Abstract. Anthropogenic climate change is projected to lead to ocean warming, acidification, deoxygenation, reductions in near-surface nutrients, and changes to primary production, all of which are expected to affect marine ecosystems. Here we assess projections of these drivers of environmental change over the twenty-first century from Earth system models (ESMs) participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6) that were forced under the CMIP6 Shared Socioeconomic Pathways (SSPs). Projections are compared to those from the previous generation (CMIP5) forced under the Representative Concentration Pathways (RCPs). A total of 10 CMIP5 and 13 CMIP6 models are used in the two multi-model ensembles. Under the high-emission scenario SSP5-8.5, the multi-model global mean change (2080–2099 mean values relative to 1870–1899) ± the inter-model SD in sea surface temperature, surface pH, subsurface (100–600 m) oxygen concentration, euphotic (0–100 m) nitrate concentration, and depth-integrated primary production is +3.47±0.78 ∘C, -0.44±0.005, -13.27±5.28, -1.06±0.45 mmol m−3 and -2.99±9.11 %, respectively. Under the low-emission, high-mitigation scenario SSP1-2.6, the corresponding global changes are +1.42±0.32 ∘C, -0.16±0.002, -6.36±2.92, -0.52±0.23 mmol m−3, and -0.56±4.12 %. Projected exposure of the marine ecosystem to these drivers of ocean change depends largely on the extent of future emissions, consistent with previous studies. The ESMs in CMIP6 generally project greater warming, acidification, deoxygenation, and nitrate reductions but lesser primary production declines than those from CMIP5 under comparable radiative forcing. The increased projected ocean warming results from a general increase in the climate sensitivity of CMIP6 models relative to those of CMIP5. This enhanced warming increases upper-ocean stratification in CMIP6 projections, which contributes to greater reductions in upper-ocean nitrate and subsurface oxygen ventilation. The greater surface acidification in CMIP6 is primarily a consequence of the SSPs having higher associated atmospheric CO2 concentrations than their RCP analogues for the same radiative forcing. We find no consistent reduction in inter-model uncertainties, and even an increase in net primary production inter-model uncertainties in CMIP6, as compared to CMIP5.

Purpose of Review The changes or updates in ocean biogeochemistry component have been mapped between CMIP5 and CMIP6 model versions, and an assessment made of how far these have led to improvements in the simulated mean state of marine biogeochemical models within the current generation of Earth system models (ESMs). Recent Findings The representation of marine biogeochemistry has progressed within the current generation of Earth system models. However, it remains difficult to identify which model updates are responsible for a given improvement. In addition, the full potential of marine biogeochemistry in terms of Earth system interactions and climate feedback remains poorly examined in the current generation of Earth system models. Summary Increasing availability of ocean biogeochemical data, as well as an improved understanding of the underlying processes, allows advances in the marine biogeochemical components of the current generation of ESMs. The present study scrutinizes the extent to which marine biogeochemistry components of ESMs have progressed between the 5th and the 6th phases of the Coupled Model Intercomparison Project (CMIP).

Abstract. This article describes the new Earth system model (ESM), the Model for Interdisciplinary Research on Climate, Earth System version 2 for Long-term simulations (MIROC-ES2L), using a state-of-the-art climate model as the physical core. This model embeds a terrestrial biogeochemical component with explicit carbon–nitrogen interaction to account for soil nutrient control on plant growth and the land carbon sink. The model's ocean biogeochemical component is largely updated to simulate the biogeochemical cycles of carbon, nitrogen, phosphorus, iron, and oxygen such that oceanic primary productivity can be controlled by multiple nutrient limitations. The ocean nitrogen cycle is coupled with the land component via river discharge processes, and external inputs of iron from pyrogenic and lithogenic sources are considered. Comparison of a historical simulation with observation studies showed that the model could reproduce the transient global climate change and carbon cycle as well as the observed large-scale spatial patterns of the land carbon cycle and upper-ocean biogeochemistry. The model demonstrated historical human perturbation of the nitrogen cycle through land use and agriculture and simulated the resultant impact on the terrestrial carbon cycle. Sensitivity analyses under preindustrial conditions revealed that the simulated ocean biogeochemistry could be altered regionally (and substantially) by nutrient input from the atmosphere and rivers. Based on an idealized experiment in which CO2 was prescribed to increase at a rate of 1 % yr−1, the transient climate response (TCR) is estimated to be 1.5 K, i.e., approximately 70 % of that from our previous ESM used in the Coupled Model Intercomparison Project Phase 5 (CMIP5). The cumulative airborne fraction (AF) is also reduced by 15 % because of the intensified land carbon sink, which results in an airborne fraction close to the multimodel mean of the CMIP5 ESMs. The transient climate response to cumulative carbon emissions (TCRE) is 1.3 K EgC−1, i.e., slightly smaller than the average of the CMIP5 ESMs, which suggests that “optimistic” future climate projections will be made by the model. This model and the simulation results contribute to CMIP6. The MIROC-ES2L could further improve our understanding of climate–biogeochemical interaction mechanisms, projections of future environmental changes, and exploration of our future options regarding sustainable development by evolving the processes of climate, biogeochemistry, and human activities in a holistic and interactive manner.

We present a model intercomparison project, LongRunMIP, the first collection of millennial-length (1,000+ years) simulations of complex coupled climate models with a representation of ocean, atmosphere, sea ice, and land surface, and their interactions. Standard model simulations are generally only a few hundred years long. However, modeling the long-term equilibration in response to radiative forcing perturbation is important for understanding many climate phenomena, such as the evolution of ocean circulation, time- and temperature-dependent feedbacks, and the differentiation of forced signal and internal variability. The aim of LongRunMIP is to facilitate research into these questions by serving as an archive for simulations that capture as much of this equilibration as possible. The only requirement to participate in LongRunMIP is to contribute a simulation with elevated, constant CO2 forcing that lasts at least 1,000 years. LongRunMIP is an MIP of opportunity in that the simulations were mostly performed prior to the conception of the archive without an agreed-upon set of experiments. For most models, the archive contains a preindustrial control simulation and simulations with an idealized (typically abrupt) CO2 forcing. We collect 2D surface and top-of-atmosphere fields and 3D ocean temperature and salinity fields. Here, we document the collection of simulations and discuss initial results, including the evolution of surface and deep ocean temperature and cloud radiative effects. As of October 2019, the collection includes 50 simulations of 15 models by 10 modeling centers. The data of LongRunMIP are publicly available. We encourage submissions of more simulations in the future.

Global warming is expected to decrease ocean oxygen concentrations by less solubility of surface ocean and change in ocean circulation. The associated expansion of the oxygen minimum zone would have adverse impacts on marine organisms and ocean biogeochemical cycles. Oxygen reduction is expected to persist for a thousand years or more, even after atmospheric carbon dioxide stops rising. However, long-term changes in ocean oxygen and circulation are still unclear. Here we simulate multimillennium changes in ocean circulation and oxygen under doubling and quadrupling of atmospheric carbon dioxide, using a fully coupled atmosphere-ocean general circulation model and an offline biogeochemical model. In the first 500 years, global oxygen concentration decreases, consistent with previous studies. Thereafter, however, the oxygen concentration in the deep ocean globally recovers and overshoots at the end of the simulations, despite surface oxygen decrease and weaker Atlantic meridional overturning circulation. This is because, after the initial cessation, the recovery and overshooting of deep ocean convection in the Weddell Sea enhance ventilation and supply oxygen-rich surface waters to deep ocean. Another contributor to deep ocean oxygenation is seawater warming, which reduces the export production and shifts the organic matter remineralization to the upper water column. Our results indicate that the change in ocean circulation in the Southern Ocean potentially drives millennial-scale oxygenation in deep ocean, which is opposite to the centennial-scale global oxygen reduction and general expectation.