Ocean ecosystems are increasingly stressed by human-induced changes of their physical, chemical and biological environment. Among these changes, warming, acidification, deoxygenation and changes in primary productivity by marine phytoplankton can be considered as four of the major stressors of open ocean ecosystems. Due to rising atmospheric CO₂ in the coming decades, these changes will be amplified. Here, we use the most recent simulations performed in the framework of the Coupled Model Intercomparison Project 5 to assess how these stressors may evolve over the course of the 21st century. The 10 Earth system models used here project similar trends in ocean warming, acidification, deoxygenation and reduced primary productivity for each of the IPCC's representative concentration pathways (RCPs) over the 21st century. For the "business-as-usual" scenario RCP8.5, the model-mean changes in the 2090s (compared to the 1990s) for sea surface temperature, sea surface pH, global O₂ content and integrated primary productivity amount to +2.73 (±0.72) °C, −0.33 (±0.003) pH unit, −3.45 (±0.44)% and −8.6 (±7.9)%, respectively. For the high mitigation scenario RCP2.6, corresponding changes are +0.71 (±0.45) °C, −0.07 (±0.001) pH unit, −1.81 (±0.31)% and −2.0 (±4.1)%, respectively, illustrating the effectiveness of extreme mitigation strategies. Although these stressors operate globally, they display distinct regional patterns and thus do not change coincidentally. Large decreases in O₂ and in pH are simulated in global ocean intermediate and mode waters, whereas large reductions in primary production are simulated in the tropics and in the North Atlantic. Although temperature and pH projections are robust across models, the same does not hold for projections of subsurface O₂ concentrations in the tropics and global and regional changes in net primary productivity. These high uncertainties in projections of primary productivity and subsurface oxygen prompt us to continue inter-model comparisons to understand these model differences, while calling for caution when using the CMIP5 models to force regional impact models.

Numerical models of ocean biogeochemistry are relied upon to make projections about the impact of climate change on marine resources and test hypotheses regarding the drivers of past changes in climate and ecosystems. In large areas of the ocean, iron availability regulates the functioning of marine ecosystems and hence the ocean carbon cycle. Accordingly, our ability to quantify the drivers and impacts of fluctuations in ocean ecosystems and carbon cycling in space and time relies on first achieving an appropriate representation of the modern marine iron cycle in models. When the iron distributions from 13 global ocean biogeochemistry models are compared against the latest oceanic sections from the GEOTRACES program, we find that all models struggle to reproduce many aspects of the observed spatial patterns. Models that reflect the emerging evidence for multiple iron sources or subtleties of its internal cycling perform much better in capturing observed features than their simpler contemporaries, particularly in the ocean interior. We show that the substantial uncertainty in the input fluxes of iron results in a very wide range of residence times across models, which has implications for the response of ecosystems and global carbon cycling to perturbations. Given this large uncertainty, iron fertilization experiments based on any single current generation model should be interpreted with caution. Improvements to how such models represent iron scavenging and also biological cycling are needed to raise confidence in their projections of global biogeochemical change in the ocean.

Among the tropical oceans, the western Indian Ocean hosts one of the largest concentrations of marine phytoplankton blooms in summer. Interestingly, this is also the region with the largest warming trend in sea surface temperatures in the tropics during the past century—although the contribution of such a large warming to productivity changes has remained ambiguous. Earlier studies had described the western Indian Ocean as a region with the largest increase in phytoplankton during the recent decades. On the contrary, the current study points out an alarming decrease of up to 20% in phytoplankton in this region over the past six decades. We find that these trends in chlorophyll are driven by enhanced ocean stratification due to rapid warming in the Indian Ocean, which suppresses nutrient mixing from subsurface layers. Future climate projections suggest that the Indian Ocean will continue to warm, driving this productive region into an ecological desert.

Past model studies have projected a global decrease in marine net primary production (NPP) over the 21st century, but these studies focused on the multi-model mean rather than on the large inter-model differences. Here, we analyze model-simulated changes in NPP for the 21st century under IPCC's high-emission scenario RCP8.5. We use a suite of nine coupled carbon–climate Earth system models with embedded marine ecosystem models and focus on the spread between the different models and the underlying reasons. Globally, NPP decreases in five out of the nine models over the course of the 21st century, while three show no significant trend and one even simulates an increase. The largest model spread occurs in the low latitudes (between 30° S and 30° N), with individual models simulating relative changes between −25 and +40 %. Of the seven models diagnosing a net decrease in NPP in the low latitudes, only three simulate this to be a consequence of the classical interpretation, i.e., a stronger nutrient limitation due to increased stratification leading to reduced phytoplankton growth. In the other four, warming-induced increases in phytoplankton growth outbalance the stronger nutrient limitation. However, temperature-driven increases in grazing and other loss processes cause a net decrease in phytoplankton biomass and reduce NPP despite higher growth rates. One model projects a strong increase in NPP in the low latitudes, caused by an intensification of the microbial loop, while NPP in the remaining model changes by less than 0.5 %. While models consistently project increases NPP in the Southern Ocean, the regional inter-model range is also very substantial. In most models, this increase in NPP is driven by temperature, but it is also modulated by changes in light, macronutrients and iron as well as grazing. Overall, current projections of future changes in global marine NPP are subject to large uncertainties and necessitate a dedicated and sustained effort to improve the models and the concepts and data that guide their development.

In this paper the interplay between tropical cyclones (TCs) and the Northern Hemispheric ocean heat transport (OHT) is investigated. In particular, results from a numerical simulation of the twentieth-century and twenty-first-century climates, following the Intergovernmental Panel on Climate Change (IPCC) twentieth-century run (20C3M) and A1B scenario protocols, respectively, have been analyzed. The numerical simulations have been performed using a state-of-the-art global atmosphere–ocean–sea ice coupled general circulation model (CGCM) with relatively high-resolution (T159) in the atmosphere. The CGCM skill in reproducing a realistic TC climatology has been assessed by comparing the model results from the simulation of the twentieth century with available observations. The model simulates tropical cyclone–like vortices with many features similar to the observed TCs. Specifically, the simulated TCs exhibit realistic structure, geographical distribution, and interannual variability, indicating that the model is able to capture the basic mechanisms linking the TC activity with the large-scale circulation. The cooling of the surface ocean observed in correspondence of the TCs is well simulated by the model. TC activity is shown to significantly increase the poleward OHT out of the tropics and decrease the poleward OHT from the deep tropics on short time scales. This effect, investigated by looking at the 100 most intense Northern Hemisphere TCs, is strongly correlated with the TC-induced momentum flux at the ocean surface, where the winds associated with the TCs significantly weaken (strengthen) the trade winds in the 5°–18°N (18°–30°N) latitude belt. However, the induced perturbation does not impact the yearly averaged OHT. The frequency and intensity of the TCs appear to be substantially stationary through the entire 1950–2069 simulated period, as does the effect of the TCs on the OHT.