Biogeochemical versus ecological consequences of modeled ocean
physics
Abstract. Regional and idealized modeling studies have shown that increasing the physical resolution of biogeochemical models to include mesoscale and submesoscale dynamics can result in both increases and decreases in phytoplankton biomass and primary production, as well as changes in phytoplankton community structure. Here we present a systematic study of the differences generated by coupling the same ecological-biogeochemical model to a 1°, coarse-resolution, and 1/6°, eddy-permitting, global ocean circulation model. Surprisingly, we find that the modeled phytoplankton community is largely unchanged, with the same phenotypes dominating in both cases. Conversely, there are large regional variations in integrated primary production, phytoplankton and zooplankton biomass. In the subtropics, mixed layer depths are, on average, deeper in the eddy-permitting model, resulting in higher nutrient supply driving increases in primary production and phytoplankton biomass. In the higher latitudes, deeper spring mixed layer depths in the eddy-permitting model result in increased light limitation during the spring bloom. Counter-intuitively, this does not drive a decrease in phytoplankton biomass, but is reflected in decreased primary production and zooplankton biomass. We explain these similarities and differences in the model using the framework of resource competition theory, and find that they are the consequence of changes in the regional and seasonal nutrient supply and light environment, mediated by differences in the modeled mixed layer depths. Although previous work has suggested that complex models may respond chaotically and unpredictably to changes in forcing, we find that our model responds in a predictable way to different ocean circulation forcing, despite its complexity.