Articles | Volume 17, issue 13
https://doi.org/10.5194/bg-17-3439-2020
https://doi.org/10.5194/bg-17-3439-2020
Research article
 | Highlight paper
 | 
06 Jul 2020
Research article | Highlight paper |  | 06 Jul 2020

Twenty-first century ocean warming, acidification, deoxygenation, and upper-ocean nutrient and primary production decline from CMIP6 model projections

Lester Kwiatkowski, Olivier Torres, Laurent Bopp, Olivier Aumont, Matthew Chamberlain, James R. Christian, John P. Dunne, Marion Gehlen, Tatiana Ilyina, Jasmin G. John, Andrew Lenton, Hongmei Li, Nicole S. Lovenduski, James C. Orr, Julien Palmieri, Yeray Santana-Falcón, Jörg Schwinger, Roland Séférian, Charles A. Stock, Alessandro Tagliabue, Yohei Takano, Jerry Tjiputra, Katsuya Toyama, Hiroyuki Tsujino, Michio Watanabe, Akitomo Yamamoto, Andrew Yool, and Tilo Ziehn
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.78C, -0.44±0.005, -13.27±5.28, -1.06±0.45mmol 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.32C, -0.16±0.002, -6.36±2.92, -0.52±0.23mmol 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.

Dates
1 Introduction

1.1 Ocean warming, acidification, deoxygenation, nutrient stress, and reduced primary production

Since the pre-industrial period the global oceans have experienced fundamental changes in physical and biogeochemical conditions as a result of anthropogenic climate change. Although these changes reflect the climate services that the oceans provide through heat and carbon storage, they also have major implications for the health of marine ecosystems. Ocean ecosystems are affected by the direct and indirect consequences of climate change. Atmospheric warming and rising CO2 concentrations drives ocean warming and acidification, while these direct factors cause changes that modulate other important components of the ocean system, such as oxygenation, nutrient levels, and net primary production.

Temperature is a principal determinant of biological metabolism in the ocean (e.g. Eppley, 1972) and plays a major role in shaping the global distribution of marine species (e.g. Thomas et al., 2012; Sunagawa et al., 2015). The radiative forcing associated with greenhouse gas emissions results in an accumulation of heat in the Earth system, most of which is taken up by the oceans (Frölicher et al., 2014). Globally averaged sea surface temperature (SST) has increased by +0.7 C over the last 100 years (Bindoff et al., 2007), with observations indicating that the heat content trend in the upper 2000 m of the ocean has increased from 0.55 to 0.68 Jm-2s-1 since 1991 (Cheng et al., 2019).

Earth system models project twenty-first century increases in SST under all of the Representative Concentration Pathways (RCPs; Bopp et al., 2013). While certain marine organisms may have the potential to acclimate to rising ocean temperatures, poleward range shifts of many species have already been observed (Gregory et al., 2009; Sorte et al., 2010), with associated declines in tropical diversity projected (Thomas et al., 2012). Concurrently, the frequency, intensity, and duration of ocean heat waves has increased in the observational record and is projected to substantially increase in the future (Frölicher et al., 2018). This has already had serious impacts on marine foundation taxa such as corals, seagrasses, and kelps (Garrabou et al., 2009; Hobday et al., 2016; Smale et al., 2019).

A consequence of ocean warming is an increase in vertical density gradients and enhanced stratification. This results in a reduction in the supply of nutrients to the euphotic zone, with enhanced nutrient limitation generally leading to observed declines in net primary production (Behrenfeld et al., 2001; Behrenfeld et al., 2006). Earth system model projections consistently show enhanced stratification and associated reductions in euphotic zone nutrient concentrations under scenarios of climate change (Bopp et al., 2001; Sarmiento et al., 2004; Cabré et al., 2014; Fu et al., 2016). This generally results in projected global reductions in net primary production that are driven by enhanced phytoplankton nutrient limitation in the low-latitude oceans (Steinacher et al., 2010; Bopp et al., 2013; Krumhardt et al. 2017; Kwiatkowski et al., 2017; Moore et al., 2018). The projected magnitude of net primary production declines is highly uncertain across model ensembles (Bopp et al., 2013; Krumhardt et al., 2017), in part due to concurrent changes in phytoplankton light and temperature limitation, as well as altered top-down grazing, all of which can compensate for nutrient-driven production declines (Taucher and Oschlies, 2011; Laufkötter et al., 2015). However, declines in phytoplankton primary production are consistently amplified in higher trophic levels such as zooplankton (Chust et al., 2014; Stock et al., 2014; Kwiatkowski et al., 2018) and fish (Lotze et al., 2019).

Dissolved oxygen in the ocean exerts a strong control on marine ecosystems. At low O2 levels, marine animals are unable to sustain aerobic metabolism, which can lead to mortality (Vaquer-Sunyer and Duarte, 2008). Oxygen levels also affect many oceanic biogeochemical cycles through an impact on redox reactions and microbial metabolism (e.g. on the nitrogen cycle, Gruber, 2004). Global warming is driving a global decline of dissolved oxygen in the ocean, referred to as ocean deoxygenation, because of a warming-induced reduction in O2 solubility and increased stratification and reduced ventilation (Keeling et al., 2010; Oschlies et al., 2018). A recent assessment, based on three different analyses (Helm et al., 2011; Schmidtko et al., 2017; Ito et al., 2017), concluded that the oxygen content over the first 1000 m of the ocean has decreased by 0.5 to 3.3 % (or 0.7 to 3.5 % between 100 and 600 m) over 1970–2010 (Bindoff, et al., in press). In coastal systems, this warming effect is exacerbated by the effects of increased loading of nutrients and organic matter, which also lead to oxygen decline and an increase in coastal ocean dead zones (Breitburg et al., 2018). Earth System model projections consistently show continuing declines in oxygen over the twenty-first century as a function of the employed scenario (Bopp et al., 2013; Cocco et al., 2013), with large uncertainties in the tropics and for the evolution of oxygen minimum zones (Cabré et al., 2015).

The uptake of carbon by the oceans affects marine chemistry via ocean acidification (Gattuso and Buddemeier, 2000; Orr et al., 2005; Doney et al., 2009), a process that increases seawater concentrations of CO2, H+, and HCO3- and reduces pH and CO32- ion concentrations. The oceans have absorbed approximately 30 % of anthropogenic carbon emissions since the pre-industrial era (Sabine et al., 2004; Khatiwala et al., 2009, 2013; Gruber et al., 2019), resulting in global surface pH declines of approximately 0.1 units (Bindoff et al., 2007). Declines in global open ocean surface pH are 0.018 units per decade over 1991–2011 (Lauvset et al. 2015), with individual time series stations exhibiting declines of 0.017 to 0.027 units per decade (Bindoff, et al., in press). Earth system models have projected twenty-first century global surface ocean pH declines of up to 0.33 units under previous high-emissions scenarios (Bopp et al., 2013), with associated changes in the seasonal cycles of seawater carbonate chemistry (McNeil and Sasse, 2016; Kwiatkowski and Orr, 2018; Landschützer et al., 2018).

The impact of ocean acidification on marine species is extensive and diverse. Calcifying species, such as echinoderms, bryozoans, and cnidarians, exhibit depressed calcification, growth, and survival under acidification (Kroeker et al., 2010; Albright et al., 2016; Kwiatkowski et al., 2016), altering the competitive balance in ecosystems (Kroeker et al., 2013). In teleost fish and marine invertebrates, ion exchange is reduced under acidification, depressing protein synthesis and metabolic rates (Langenbuch et al., 2006; Pörtner, 2008). Physiological and behavioural functioning is also sensitive to acidification, with olfactory discrimination (Munday et al., 2009) and predator–prey responses (Watson et al., 2014, 2017) shown to be impaired under more acidified conditions.

Marine organisms typically experience changes in multiple physical and geochemical conditions simultaneously, with impacts determined by the interactions between potential stressors. For example, the combined effect of warming and deoxygenation is projected to force poleward and vertical contractions of metabolically viable habitat for marine ectotherms (Deutsch et al., 2015). At the physiological level, experimental studies indicate that synergistic effects between potential marine stressors are common (Gunderson et al., 2016). Compound warming and acidification has been shown to exacerbate negative impacts on photosynthesis, calcification, reproduction, and survival of marine organisms (Harvey et al., 2013), while compound exposure to acidification and low oxygen can also have synergistic effects (McBryan et al., 2013) and may reduce the thermal tolerance of certain species (Pörtner, 2010).

Here we assess future projections of climate-related drivers of marine impacts within the Coupled Model Intercomparison Project Phase 6 (CMIP6; Eyring et al., 2016; O'Neill et al., 2016) simulations, evaluating how these differ from previous CMIP5 (Taylor et al., 2011) simulations. We focus on projected changes in ocean temperature, pH, dissolved O2 and NO3- concentration, and net primary production across 13 CMIP6 and 10 CMIP5 Earth system models. Expanding on previous CMIP intercomparison studies (e.g. Steinacher et al., 2010; Bopp et al., 2013; Cocco et al., 2013), which typically concentrate on upper-ocean mean state changes in biogeochemistry, we also assess benthic impacts and changes in seasonal cycles.

1.2 Ocean biogeochemical model development since CMIP5

A comprehensive assessment of changes between CMIP5 and CMIP6 in the ocean biogeochemical components of Earth system models (ESMs) and their associated skill is provided in Séférian et al. (2020). Since CMIP5, CMIP6 has seen a general increase in the horizontal grid resolution of physical ocean models and a limited increase in vertical resolution. The latter may be particularly important for ecosystem projections as it directly affects simulated stratification, a key factor influencing changes in ocean impact drivers (Capotondi et al., 2012; Bopp et al., 2013; Laufkötter et al., 2015; Kwiatkowski et al., 2017) and their impact on higher trophic levels (Stock et al., 2014; Chust et al., 2014; Kwiatkowski et al., 2018; Lotze et al., 2019). Updates in the representation of ocean biogeochemical processes between CMIP5 and CMIP6 have generally included increases in model complexity (Séférian et al., 2020). Specifically, CMIP6 models provide more widespread inclusion of dissolved oxygen, micronutrients, such as iron, variable stoichiometric ratios, and improved representation of lower trophic levels including heterotrophic bacteria and the cycling and sinking of organic matter (Séférian et al., 2020).

Relative to CMIP5, the CMIP6 Earth system models display an improved ability to reproduce the modern mean-state distribution of a number of key biogeochemical tracers (Séférian et al., 2020). Although global-scale totals of ocean carbon flux and net primary production estimates have been improved in CMIP6 with respect to CMIP5, the simulated geographical distribution of present-day mean state air–sea CO2 fluxes and surface chlorophyll concentrations show only moderate improvements between CMIP5 and CMIP6. There are also moderate improvements in the representation of subsurface dissolved oxygen concentrations in most ocean basins. Model skill in the representation of surface macronutrient concentrations in CMIP6 has improved for dissolved silicate but declined slightly for nitrate.

2 Methodology

The analysis of projected multiple ocean impact drivers presented here focuses on three key depth levels: the upper ocean, the thermocline, and the benthic zone. The surface zone is where most biological activity is concentrated in the oceans and where impacts from climate change are typically the greatest. Specifically, we assess projections of surface ocean temperature, surface ocean pH, subsurface dissolved O2 concentration (averaged between 100 and 600 m), upper-ocean NO3- concentration (averaged between 0 and 100 m) and net primary production (depth integrated over the full water column). The choice of vertical range for O2 reflects the potential importance of the expansion of oxygen minimum zones, which are more prominent at such depths. The choice of vertical range for NO3- reflects its importance as a critical macronutrient supporting primary production in the euphotic zone. Both vertical ranges are chosen to be compatible with the recent assessment of marine drivers in the IPCC Special Report on the Ocean and Cryosphere (Bindoff, et al., in press). Additionally, for the CMIP6 models we assess benthic ecosystem drivers, focussing on projections of bottom temperature, pH, and O2 concentration. The benthic level is defined as the bottom-ocean model layer at each grid point. As such, its exact depth depends on vertical discretisation and bathymetry, which differs across the CMIP6 ensemble. All benthic model outputs were corrected for potential drift at the grid cell level (e.g. Gehlen et al., 2014; Séférian et al., 2016) using coincident pre-industrial control simulations.

2.1 Processing and analysis of model outputs

All ESMs assessed in the CMIP5 and CMIP6 ensembles (Tables 1 and 2) include physical ocean models and coupled ocean biogeochemistry schemes that account for some or all of the potential ocean impact drivers: temperature, pH, O2, NO3-, and net primary production. A total of 10 CMIP5 and 13 CMIP6 models are assessed with the model ensemble size differing among scenarios depending on contributions from each model group. The CMIP5 ensemble is the same as that used in the comprehensive assessment of projected ocean drivers provided by Bopp et al. (2013). Only one ensemble member per model is used for a given scenario. That is, in CMIP terminology we typically use ensemble member “r1i1p1” from each CMIP5 model and “r1i1p1fx” from each CMIP6 model (where “fx” is the recommended set of external forcings employed by the various modelling groups). Consequently, we do not assess the role of internal variability in the emergence of climate-related changes in marine ecosystems drivers (e.g. Frölicher et al., 2016; Lovenduski et al., 2016; Krumhardt et al. 2017; Freeman et al., 2018). Two of the CMIP6 models included in our analysis (GFDL-CM4 and ACCESS-ESM1.5) do not include NO3- as a prognostic tracer. Hence, their NO3- concentrations were calculated from modelled total dissolved inorganic phosphorus assuming a constant Redfield ratio of 16:1.

Table 1The CMIP6 Earth system models used in this study; their individual components used to represent ocean, sea ice, and marine biogeochemistry; and the ocean impact drivers and simulations that were assessed.

Download Print Version | Download XLSX

Table 2The CMIP5 Earth system models used in this study; their individual components used to represent ocean, sea ice, and marine biogeochemistry; and the simulations that were assessed. All models provided temperature, pH, oxygen, nitrate, and NPP outputs, as in Bopp et al. (2013).

Download Print Version | Download XLSX

To facilitate intercomparison, model output on each native grid was regridded to the same regular 1×1 horizontal grid using distance weighted average remapping (climate data operators; remapdis). Model outputs were kept on their native vertical grids, with vertical discretisation ranging from 40 (MPI-ESM1.2) to 75 (IPSL-CM6A-LR, CNRM-ESM2-1 and UKESM1-0-LL) levels, except for models using hybrid or isopycnic vertical coordinates for which model outputs were vertically regridded to 35 (GFDL-ESM4, GFDL-CM4) and 70 (NorESM2-LM) levels. Following generally adopted practice (e.g. Bopp et al., 2013), all models were given equal weighting in the respective CMIP6 and CMIP5 ensemble mean. However, within the CMIP6 ensemble two modelling groups contributed two ESMs and within the CMIP5 ensemble three modelling groups contributed two ESMs, which is likely to influence the extent of model independence (Masson and Knutti, 2011; Knutti et al., 2015; Sanderson et al., 2015; Lovenduski and Bonan, 2017).

The CMIP5 historical simulations had variable start dates between 1850 and 1861, all of which finished in 2005; the subsequent RCP simulations started in 2006 and were run until at least 2099. In CMIP6, there is greater temporal consistency. All CMIP6 historical simulations were made over 1850–2014, while the subsequent SSP scenarios started in 2015 and ran until at least 2100. To facilitate comparison between CMIP5 and CMIP6, the historical and future projections of ocean impact drivers in both phases of CMIP are presented as anomalies relative to 1870–1899 mean values of their respective historical simulations. When solely evaluating twenty-first century projections in the SSPs, however, the last 20 years of the CMIP6 historical simulations (1995–2014) are used as a baseline period. Throughout the analysis, the uncertainty associated with global mean projections is assessed using the inter-model SD (given ± uncertainties). At regional scales, projection robustness is evaluated using previously adopted approaches (e.g. Bopp et al., 2013), including whether the magnitude of the multi-model mean anomaly exceeds the inter-model SD or if there is at least 80 % model sign agreement. The interquartile range of regional projections is given in the annexes.

2.2 From Representative Concentration Pathways to Shared Socioeconomic Pathways

Aside from changes in ESMs, a fundamental difference between CMIP5 and CMIP6 is that they differ in the future scenarios used for anthropogenic emissions and land use change. Those scenarios are derived from integrated assessment models and based on plausible future pathways of societal development. In CMIP6, the Shared Socioeconomic Pathways (SSPs) provided via the Scenario Model Intercomparison Project (ScenarioMIP) are used instead of the RCPs that were used in CMIP5 (O'Neill et al., 2016). The SSPs provide revised emission and land use scenarios relative to the RCPs (Riahi et al., 2017).

In this study, we confine our assessment of ocean impact drivers to concentration-driven simulations, focussing on Tier 1 SSPs of ScenarioMIP (O'Neill et al., 2016), that is SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5, which result in end-of-century approximate total radiative forcing levels of 2.6, 4.5, 7.0, and 8.5 W m−2, respectively. The SSPs have generally higher associated concentrations of atmospheric CO2 and lower associated atmospheric concentrations of CH4 and N2O relative to their RCP counterparts (O'Neill et al., 2016; Meinshausen et al., 2019). This is particularly the case for SSP5-8.5, which in comparison to RCP8.5 assumes that coal constitutes a greater proportion of the primary energy mix in the second half of the twenty-first century (Kriegler et al., 2017). Given that differences among projections of surface ocean acidification are dominated by scenario uncertainty, with relatively little inter-model uncertainty and internal variability (e.g. Bopp et al., 2013; Frölicher et al., 2016), such changes in atmospheric concentrations of CO2 are expected to have a large impact on projections of ocean pH and related carbonate system variables.

Alongside the assessment of the SSP concentration-driven model outputs, outputs from models forced under the four RCPs (RCP2.6, RCP4.5, RCP6.0 and RCP8.5) are also assessed in parallel. This allows some comparison with past CMIP5 assessments (e.g. Bopp et al., 2013). However, RCP6.0 has no direct SSP analogue, while SSP3-7.0 has no direct RCP analogue.

3 Results and discussion

3.1 Comparison with historical global trends

Observed historical trends in global mean surface ocean temperature, pH, and subsurface oxygen were compared with the multi-model mean of the CMIP6 ensemble over the corresponding years of historical simulations (Table 3). Global observations of historical trends in euphotic-zone nitrate concentrations and integrated primary production were deemed insufficiently robust, given the associated interannual–decadal variability, to be assessed in the models (Elsworth et al. 2020). The observed 1901–2012 SST warming of +0.06C per decade is well reproduced in the CMIP6 ensemble, although SST warming between 1979 and 2012 is warm-biased in the multi-model mean. Global surface ocean acidification of −0.018 units per decade between 1991 and 2011 is also well reproduced by the CMIP6 models, particularly considering that the observed trend is a reconstruction based on discrete surface ocean fCO2 measurements and alkalinity estimates that are derived from temperature and salinity (Lauvset et al., 2015). Finally, with respect to subsurface deoxygenation, the observed dissolved oxygen trend of −0.30 to −1.52mmol m−3 per decade from 1970 to 2010 (90 % confidence range; Bindoff, et al., in press) encompasses the CMIP6 multi-model mean response over the corresponding years. Given the performance of the CMIP6 models at reproducing ocean biogeochemical mean conditions (Séférian et al., 2020) and trends, they are deemed appropriate to project future trends in biogeochemistry under the SSPs.

Table 3Historical trends in ocean impact drivers in observations and CMIP6 models. Observed multi-decadal historical trends in global mean sea surface temperature, surface pH, and subsurface O2 (averaged between 100 and 600 m) and the corresponding trend in the CMIP6 ensemble. Uncertainty estimates are the inter-model SD.

Download Print Version | Download XLSX

3.2 Global upper-ocean projections

Under all SSPs, global multi-model mean sea surface temperature is projected to increase, while surface pH, subsurface dissolved oxygen concentration, euphotic-zone nitrate concentration and net primary production are projected to decline during the twenty-first century (Fig. 1). The projected change in the five ocean impact drivers increases with associated radiative forcing across the four SSPs. Under the high-mitigation SSP1-2.6 scenario, the end-of-century model mean changes (2080–2099 mean values relative to 1870–1899) in sea surface temperature, surface pH, subsurface oxygen concentration, euphotic nitrate concentration, and net primary production are +1.42±0.32C, -0.16±0.002, -6.36±2.92, -0.52±0.23mmol m−3, and -0.56±4.12 %, respectively. Under the high-emissions scenario SSP5-8.5, the corresponding changes are 3.47± 0.78 C, -0.44±0.005, -13.27±5.28, -1.06±0.45mmol m−3, and -2.99±9.11 % (Table 4), respectively. Across these two scenarios, the separation between CMIP6 projections of sea surface temperature and pH and to a lesser extent oxygen and nitrate, further demonstrate the effectiveness of intense mitigation strategies in limiting twenty-first century marine ecosystem exposure to potential stress. This is in agreement with assessments of previous multi-model projections (e.g. CMIP5; Bopp et al., 2013).

Table 4Global mean changes in multiple ocean impact drivers across the CMIP6 and CMIP5 ensembles. Global mean anomalies of sea surface temperature, surface ocean pH, subsurface dissolved O2 concentration (averaged between 100 and 600 m), upper-ocean NO3- (averaged between 0–100 m), and depth-integrated net primary production for the CMIP6 SSPs and CMIP5 RCPs. Anomalies are given as 2080–2099 mean values relative to the 1870–1899 mean. Uncertainty estimates are the inter-model standard deviation.

Download Print Version | Download XLSX

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f01

Figure 1Global mean projections of upper-ocean impact drivers. Global mean projections of (a) sea surface temperature (C), (b) surface ocean pH, (c) subsurface dissolved O2 concentration (averaged between 100–600 m;mmolm-3), (d) euphotic-zone NO3- (averaged between 0–100 m;mmolm-3), and (e) depth-integrated net primary production ( %). Values are anomalies relative to the 1870–1899 reference period. CMIP6 mean anomalies for the historical and SSP simulations are shown as solid lines with shading representing the inter-model SD. CMIP5 projections only show the multi-model mean. The model ensemble size for each scenario is given in parentheses.

Download

Following previous assessments (Bopp et al., 2013), inter-model uncertainty is estimated as the inter-model SD. Although some of this model spread is due to internal variability, this contribution is limited for global averages and its relative contribution to inter-model uncertainty is expected to decline throughout the twenty-first century (Frölicher et al., 2016). Relative to scenario uncertainty, which is estimated as the maximum difference between mean SSP projections, inter-model uncertainty is extremely low for surface pH projections, which show distinct separation between the SSPs prior to 2050. The low inter-model uncertainty associated with projections of surface ocean pH is well characterised and associated with the identical CO2 forcing used by all ESMs in concentration-driven SSP and RCP projections (Lovenduski et al., 2016), a weak climate–pH feedback (Orr et al., 2005; McNeil and Matear, 2007), limited interannual variability, and consistently adopted standards for ESM ocean carbonate chemistry equations (Orr et al., 2017). Surface ocean pCO2 and corresponding carbonate chemistry generally follow changes in atmospheric CO2 with a global mean equilibration time of approximately 8 months (Gattuso and Hansson, 2011). The differences between projected surface pH across the SSPs therefore reflect the divergence of prescribed atmospheric CO2 concentrations, i.e. the different scenarios.

In comparison to pH, projections of SST exhibit greater inter-model uncertainty (Fig. 1). This uncertainty is likely to result from differences in climate sensitivity between models. Historically, such differences have been attributed to diversity in cloud feedbacks and to a lesser extent water vapour and lapse-rate feedbacks (Andrews et al., 2012; Vial et al., 2013). For projections of subsurface oxygen and euphotic-zone nitrate concentrations, inter-model uncertainty is greater still and can exceed scenario uncertainty. This greater inter-model uncertainty is a result of oxygen and nitrate concentrations being strongly influenced by both physical changes (e.g. changes in solubility, circulation, and mixing) and changes in biological sources and sinks (Stramma et al., 2012; Fu et al., 2016; Bopp et al., 2017; Oschlies et al., 2018).

The inter-model uncertainty associated with CMIP6 net primary production projections is consistently larger than the scenario uncertainty. Indeed, for each SSP, individual models project both increases and decreases in global primary production, with the inter-model SD encompassing positive and negative anomalies (Fig. 1). This is a consequence of net primary production changes reflecting a diverse and delicately balanced suite of bottom-up and top-down ecological processes, which are variously parameterised across models. Bottom-up changes in phytoplankton growth rates are typically driven by a combination of enhanced nutrient limitation and reduced temperature and light limitation (Doney, 2006; Steinacher et al., 2010; Taucher and Oschlies, 2011), while top-down changes in zooplankton grazing rates can simultaneously influence the stock of phytoplankton biomass (Laufkötter et al., 2015). The accurate simulation of many of the biogeochemical tracers upon which net primary production (NPP) depends (e.g. the distribution of iron; Tagliabue et al., 2016) represents a significant and ongoing challenge to ESMs (Séférian et al., 2020).

3.3 Regional patterns of upper-ocean change

Global-scale projections of end-of-century upper-ocean impact drivers (2080–2099 anomalies relative to 1995–2014 mean values) exhibit spatial variability that is both ocean impact driver and SSP dependent (Fig. 2). CMIP6 projections of SST show near global increases under both SSP1-2.6 and SSP5-8.5, with generally high regional robustness across the model ensemble. The greatest warming is evident in the Northern Hemisphere, particularly the Arctic Ocean and high-latitude North Pacific, where multi-model mean warming can exceed 2 C in SSP1-2.6 and 5 C in SSP5-8.5. This Arctic amplification is well established in both observations (Bekryaev et al., 2010) and models, and thought to be primarily driven by temperature and surface albedo feedbacks (Screen and Simmonds, 2010; Pithan and Mauritsen, 2014). The notable exception to warming is in the subpolar North Atlantic where there is minor cooling in SSP1-2.6 and limited warming in SSP5-8.5. Although it is associated with high inter-model uncertainty (Figs. 2b,c and A1), this “warming hole” is well documented in both observations and models and typically related to a slow down in the Atlantic meridional overturning circulation (Drijfhout et al., 2012; Menary and Wood, 2018). Spatial patterns of SST anomalies are broadly consistent with those of the CMIP5 ensemble (Bopp et al., 2013).

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f02

Figure 2Projections of multiple upper-ocean impact drivers under SSP1-2.6 and SSP5-8.5. CMIP6 mean historical climatologies and twenty-first century anomalies in (a–c) sea surface temperature (C), (d–f) surface ocean pH, (g–i) subsurface dissolved O2 concentration (averaged between 100 and 600 m;mmolm-3), (j–l) euphotic-zone NO3- (averaged between 0 and 100 m;mmolm-3), and (m–o) depth-integrated net primary production (gCm-2y-1). Anomalies are 2080–2099 mean values relative to the 1995–2014 baseline period. Stippling designates areas of projection robustness. For temperature and pH this is defined as the magnitude of the mean anomaly exceeding the inter-model SD. For O2, NO3, and NPP this is defined as at least 80 % model sign agreement.

Anomalies in surface ocean pH are ubiquitously negative under both SSP1-2.6 and SSP5-8.5, with very low associated inter-model uncertainty (Figs. 2e and f and A1). Consistent with past model projections (McNeil and Matear, 2007; Steinacher et al., 2009; Bopp et al., 2013), the greatest declines in pH are projected in the higher latitudes and the Arctic Ocean in particular, where model mean declines can exceed 0.45 in SSP5-8.5 (2080–2099 anomalies relative to 1995–2014; Fig. 2e, f). The enhanced Arctic Ocean acidification reflects both the large surface ocean warming in the Arctic, which acts to decrease pH, as well as the related loss of permanent and semi-permanent sea ice (McNeil and Matear, 2007; Steinacher et al., 2009). Sea ice melt increases anthropogenic carbon uptake and decreases pH by both providing a greater surface area for air–sea gas exchange and simultaneously enhancing air–sea CO2 fluxes by dilution of dissolved inorganic carbon concentrations with freshwater (Yamamoto-Kawai et al., 2009; Yamamoto et al., 2012).

Although global mean subsurface (100–600 m) O2 concentration is projected to decline under all SSPs, there is a high degree of variability in projections at regional scales (Fig. 2h,i). The largest declines in subsurface O2 generally occur at higher latitudes and in particular in the North Pacific, where declines in the multi-model mean can exceed 45 mmol m−3 in SSP5-8.5. In equatorial regions of the Atlantic and Indian Ocean and upwelling regions of the Pacific, increases in subsurface O2 concentration are projected under both SSP1-2.6 and SSP5-8.5. However, these increases have high inter-model uncertainty and are at odds with historical observations of OMZ expansion (Stramma et al., 2008; Helm et al., 2011; Schmidtko et al., 2017; Ito et al., 2017), with an assessed range of 3.0 to 8.3 % between 1970 and 2010 (Bindoff et al., in press), though it has been suggested that such observations are the result of climate variability (Deutsch et al., 2011; Bindoff et al., in press). These subsurface O2 increases are, however, consistent with previous projections, including those from CMIP5, which have highlighted that coarse-resolution models struggle to reproduce subsurface ventilation pathways in these regions (Stramma et al., 2012; Andrews et al., 2013; Busecke et al., 2019).

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f03

Figure 3Change in subsurface oxygen saturation and apparent oxygen utilisation under SSP5-8.5. CMIP6 multi-model mean changes in (a) subsurface dissolved O2 concentration (averaged between 100 and 600 m;mmolm-3), (b) subsurface O2 saturation (O2 sat), and (c) subsurface apparent oxygen utilisation (AOU) in 2080–2099 of SSP5-8.5 relative to 1995–2014. Stippling designates robustness, as defined by at least 80 % model sign agreement (only a subset of CMIP6 models were considered).

For a subset of the CMIP6 models (CanESM5, CNRM-ESM2-1, GFDL-CM4, IPSL-CM6A-LR, MIROC-ES2L, MPI-ESM1.2-HR and UKESM1-0-LL), projected changes in subsurface O2 concentration under SSP5-8.5 were decomposed into changes in O2 saturation (O2 sat) and apparent oxygen utilisation (AOU), where ΔO2=ΔO2 sat-ΔAOU (Fig. 3). O2 sat was computed from model temperature and salinity outputs and represents the effect of oxygen solubility changes on dissolved O2 concentration, while AOU was calculated as ΔAOU=ΔO2 sat-ΔO2, and is affected by both changes in biological consumption of O2 and in ventilation and stratification. The heightened reductions in subsurface O2 in the North Pacific and North Atlantic are shown to be the result of consistent reductions in O2 sat and increases in AOU, which act to reinforce O2 concentration declines. In contrast, the projected increases in O2 in the tropical Indian and Atlantic Oceans are shown to be the result of reductions in AOU that lower oxygen demand more than the concurrent reductions in O2 sat. These tropical reductions in AOU are generally robust across the model ensemble despite this not being the case for the coincident increases in O2. The spatial patterns of CMIP6 projected changes in subsurface O2 sat and AOU under SSP5-8.5 are similar to that of the CMIP5 models under RCP8.5 (Bopp et al., 2017). The general reduction in O2 sat has been shown to be predominantly due to warming-driven reductions in solubility, while the heightened AOU declines in the North Pacific and North Atlantic have been primarily attributed to reductions in ventilation and an increase in the age of these waters (Bopp et al., 2017; Tjiputra et al., 2018).

CMIP6 multi-model mean projections of NO3- concentrations in the euphotic zone (0–100 m) show variable regional declines under SSP1-2.6 and SSP5-8.5 (Fig. 2k,l). These projected declines are robust and largest in the Arctic Ocean, equatorial eastern Pacific, North Atlantic, and North Pacific, where they can exceed 3 mmol m−3 in SSP5-8.5. NO3- concentrations show limited anomalies in the subtropical gyres where concentrations are already very low. The CMIP6 spatial pattern of euphotic-zone NO3- anomalies is in broad agreement with CMIP5 projections (Fu et al., 2016).

Projections of primary production anomalies are highly diverse across regions (Fig. 2n,o). The global CMIP6 multi-model mean decline in primary production is shown to be primarily driven by declines in the North Atlantic and the western equatorial Pacific, which can exceed 40 gCm-2y-1 under SSP5-8.5. In the high latitudes, primary production generally increases, with anomalies approaching 20 gCm-2y-1 in parts of the Arctic and Southern Oceans under SSP5-8.5. Such changes have historically been associated with enhanced stratification as the upper ocean warms (Doney, 2006). In tropical and mid-latitude waters, where phytoplankton are nutrient limited, this tends to reduce the vertical nutrient supply and exacerbate nutrient stress. In contrast, in high-latitude waters, where phytoplankton are typically light limited, enhanced stratification can reduce mixing below the euphotic depth and therefore result in reduced light stress. However, the simplicity of this paradigm has been challenged by observations, which show regionally variable coupling between changes in stratification and primary production on interannual timescales (Lozier et al., 2011; Dave and Lozier, 2013), with recent studies demonstrating the additional importance of changes to the horizontal advection of nutrients (Whitt, 2019) and zooplankton grazing (Laufkötter et al., 2015) in shaping regional primary production responses.

Although the general pattern of NPP changes is similar in CMIP6 compared to CMIP5, regional declines are reduced in magnitude, less spatially extensive, and are typically less robust. This is particularly notable in the Indian Ocean and subtropical North Pacific, which were regions of consistent NPP decline in CMIP5 projections (Bopp et al., 2013) but exhibit limited robust trends in CMIP6. The projected increase in NPP in the high latitudes is, however, broadly consistent with previous model intercomparisons (Steinacher et al., 2010; Bopp et al., 2013), which have attributed these increases to the retreat of sea ice, reduced deep mixing and corresponding reductions in light limitation.

3.4 Stratification and mixed-layer depth

The difference between densities at the surface, 200, and 1000 m are used as stratification indices in the CMIP6 projections and given alongside changes in the maximum mixed-layer depth (Fig. 4). Global multi-model mean stratification is consistently projected to increase over the twenty-first century, with greater increases under SSPs that have higher radiative forcing. Under SSP1-2.6, the global change in stratification index is +0.13±0.05kg m−3 between the surface and 200 m and +0.26±0.08kg m−3 between the surface and 1000 m, while under SSP5-8.5 this increase to +0.58±0.11 and +0.90±0.20kg m−3, respectively. Over the same period, the global mean maximum annual mixed-layer depth shoals by 7.0±3.3 m in SSP1-2.6 and 19.5±2.6 m in SSP5-8.5. With the exception of the Arctic Ocean, multi-model projections of increasing stratification are typically consistent across most regions, with robustness increasing under SSP5-8.5 and when the upper 1000 m of the water column is considered. Projected shoaling of the maximum mixed-layer depth is also generally robust across the multi-model ensemble, albeit with slightly less model consistency, as would be expected given the coincident climatic changes in wind stress and surface heat fluxes (Fig. 4h,i).

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f04

Figure 4Change in stratification and mixed-layer depth in SSP1-2.6 and SSP5-8.5. CMIP6 mean historical climatologies and anomalies in (a–c) stratification index between 200 m and the surface (kg m−3), (d–f) stratification index between 1000 m and the surface (kg m−3), and (g–i) maximum annual mixed-layer depth (m). Anomalies are 2080–2099 mean values relative to 1995–2014. The stratification index is defined as the difference in density between given depths and the surface. Stippling designates robustness, as defined by the magnitude of the mean anomaly exceeding the inter-model SD.

The regions of projected primary production decline in the North Atlantic and western equatorial Pacific are typically associated with heightened increases in stratification, notably in the upper 200 m of the water column. In these regions, particularly the North Atlantic, the maximum mixed-layer depth also shoals. As such, there is strong evidence that reduced vertical mixing and entrainment of nutrients into the upper ocean is, at least partially, responsible for these regional declines in primary production. However, similar increases in stratification and reductions in mixed-layer depth occur in regions such as the North Pacific and Indian Ocean, where declines in primary production are largely absent. Therefore, further assessment of simultaneous changes in processes such as nutrient advection (e.g. Whitt, 2019), nitrogen fixation (Riche and Christian, 2018), the microbial loop (e.g. Schmittner et al., 2008; Taucher and Oschlies, 2011), and top-down grazing pressure (e.g. Laufkötter et al., 2015) are required to fully understand the regional primary production response in CMIP6.

3.5 Compound stressors

The projected occurrence of multiple potential ecosystem stressors in the upper ocean was determined across the SSPs using prescribed thresholds of surface warming (>+2C), surface acidification (<-0.2 units), subsurface deoxygenation (<-30mmol m−3) and euphotic-zone NO3- decline (<-1mmol m−3), with anomalies calculated as 2080–2099 mean values relative to 1995–2014 (Fig. 5). It should be noted that our choice of stressor thresholds, based on the magnitude of biogeochemical anomalies, are somewhat arbitrary. Indeed, it could be argued that absolute biogeochemical thresholds, for example as defined for hypoxia or oligotrophy, may better reflect potential ecosystem stress. Moreover, the thresholds take no account of regional differences in natural temperature, pH, O2, and NO3- variability, which may mediate ecosystem responses to changes in mean conditions (e.g. Kroeker et al., 2020). That being said, a single threshold that encompasses the variety of ecosystem responses to a particular stressor likely does not exist.

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f05

Figure 5Compound upper-ocean impact drivers. Regions where projected CMIP6 sea surface warming exceeds 2 C, euphotic-zone (0–100 m) NO3- decline exceeds 1 mmol m−3, surface ocean pH decline exceeds 0.2, and subsurface (100–600 m) dissolved O2 concentration decline exceeds 30 mmol m−3 in (a) SSP1-2.6, (b) SSP2-4.5, (c) SSP3-7.0, and (d) SSP5-8.5. The exceedance of driver thresholds is determined from 2080 to 2099 anomalies relative to 1995–2014 values.

The concurrent exceedance of multiple thresholds increases with associated radiative forcing across the SSPs, indicative of greater potential compound ecosystem stressors. The tropical and subtropical oceans are generally characterised by projected compound warming and acidification under SSP3-7.0 and SSP5-8.5, with additional nutrient thresholds exceeded in regions of equatorial upwelling. The North Pacific is characterised by high sensitivity to potential compound stressors, with all thresholds of warming, acidification, deoxygenation, and nutrient decline exceeded under SSP5-8.5. In contrast, the projected occurrence of compound stressors is limited in the Southern Ocean, where only the acidification threshold is consistently exceeded. The North Atlantic is characterised by sensitivity to combined acidification and nutrient stress, while the Arctic Ocean is sensitive to compound warming, acidification, and nutrient stress.

3.6 CMIP6 vs. CMIP5 projections

While the temporal behaviour of changes in ocean impact drivers is similar across the CMIP5 and CMIP6 model ensembles (Fig. 1), the CMIP6 Earth system models generally project greater global surface ocean warming, acidification, subsurface deoxygenation, and euphotic zone NO3 reductions than the CMIP5 projections performed with comparable radiative forcing (Fig. 6, Table 4). The CMIP6 models, however, project reduced global primary production declines relative to comparable CMIP5 simulations. There is no consistent reduction in inter-model uncertainty in CMIP6. In fact, with respect to projections of primary production, inter-model uncertainty is substantially increased in CMIP6.

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f06

Figure 6Comparison between CMIP6 and CMIP5 end-of-century changes in upper-ocean impact drivers. (a) Atmospheric CO2 concentration and radiative forcing derived from the MAGICC6 model in year 2100 for the CMIP6 SSPs and CMIP5 RCPs. (b–e) Global mean anomalies of (b) surface ocean pH, (c) subsurface dissolved O2 concentration (averaged between 100 and 600 m;mmolm-3), (d) NO3- concentration (averaged between 0 and 100 m;mmolm-3), and (e) integrated net primary production ( %) against anomalies of sea surface temperature (C). Anomalies are 2080–2099 mean values relative to the 1870–1899 baseline period. Error bars represent the inter-model SD.

Download

The projected end-of-century SST increase (2080–2099 minus 1870–1899) in SSP1-2.6, SSP2-4.5, and SSP5-8.5 is higher than in RCP2.6, RCP4.5 and RCP8.5, respectively. This enhanced CMIP6 warming is attributable to generally greater climate sensitivity in the CMIP6 model ensemble relative to the CMIP5 ensemble (Forster et al., 2019). Indeed, the same version of a reduced complexity climate model (MAGICC7.0) run with CMIP5 and CMIP6 forcings, projects marginally greater warming of near-surface air temperatures in the RCPs than comparative SSPs (Meinshausen et al., 2019), underlining that the greater SST increases in CMIP6 are likely driven by changes to models and not forcing datasets.

The enhanced acidification in CMIP6 relative to CMIP5 is consistent across models (Fig. 6, Table 4) and attributable to higher prescribed atmospheric CO2 levels in the forcing of the SSP scenarios relative to the RCP scenarios with equivalent radiative forcing (Meinshausen et al., 2019). Year 2100 atmospheric CO2 levels are 1135.2, 602.8, and 445.6 ppm in SSP5-8.5, SSP2-4.5, and SSP1-2.6, respectively. The corresponding levels in RCP8.5, RCP4.5, and RCP2.6 are 936, 538, and 421 ppm, respectively (Meinshausen et al., 2011). Therefore, although the SSP and RCP simulation pairs have analogous end-of-century radiative forcing, the higher CO2 levels in the SSPs result in greater acidification for the CMIP6 projections.

The greater euphotic-zone NO3- concentration declines in SSPs compared to their RCP analogues are likely a consequence of the enhanced surface warming in CMIP6 models. This warming results in a greater increase in upper-ocean stratification than that projected in CMIP5 models (Cabré et al., 2014; Fu et al., 2016). At the global scale, models have been shown to consistently exhibit strong negative correlations between relative stratification anomalies and relative NO3- anomalies on interannual timescales (Fu et al., 2016). The greater increases in stratification in CMIP6 therefore result in greater reductions in mixing and entrainment of nutrient-rich deep waters into the euphotic zone in comparison with CMIP5.

The enhanced subsurface deoxygenation in the SSPs relative to comparable RCPs is likely the consequence of both physical and biogeochemical processes (e.g. Bopp et al., 2017; Oschlies et al., 2018). The greater warming in CMIP6 projections results in a greater reduction in O2 solubility, while also affecting the ventilation and transport of O2 within the ocean interior. In addition, concurrent changes in biological production, export and respiration can either mitigate or exacerbate physically driven subsurface deoxygenation (Oschlies et al., 2018).

The reduced declines in global net primary production projected for the twenty-first century in the SSPs, relative to comparable RCPs, combined with large increases in the associated inter-model uncertainty, is striking (Fig. 6e), particularly given that declines in euphotic zone NO3- concentrations are typically greater in the SSPs. This suggests that the temporal evolution of phytoplankton resource limitation and grazing pressure under climate change may have significantly altered between CMIP5 and CMIP6. In previous CMIP biogeochemistry intercomparisons, all models projected global primary production declines, albeit with large inter-model uncertainty (Steinacher et al., 2010; Bopp et al., 2013). However, in CMIP6, four of the models (CESM2, CESM2-WACCM, CNRM-ESM2-1, and IPSL-CM6A-LR) consistently project global increases in primary production across the SSPs and are primarily responsible for both the reduced multi-model mean declines and the large increase in inter-model SD. Global increases in net primary production have been previously documented in other model studies (Schmittner et al., 2008; Taucher and Oschlies, 2011; Laufkötter et al., 2015) and attributed to temperature-dependent intensification of the microbial loop increasing regenerated production. Further analysis of the CMIP6 models that project primary production increases is clearly required to determine whether this is also the case or additional processes (e.g. the temporal evolution of nitrogen fixation or iron limitation) explain why they differ from previous generations of the same Earth system model family.

3.7 Global benthic ocean projections

On average, bottom waters are consistently projected to warm, acidify, and deoxygenate across the twenty-first century (Fig. 7). Under SSP1-2.6, the end-of-century model mean changes (2080–2099 relative to 1870–1899) in bottom-water temperature, pH, and dissolved O2 are +0.12±0.03C, -0.018±0.001, and -5.14±2.04mmol m−3, respectively. Under SSP5-8.5 the corresponding changes are +0.22±0.04C, -0.030±0.002, and -6.04±2.19mmol m−3 (Table 5). Thus, even for bottom waters, CMIP6 projections highlight that intense mitigation strategies can limit ecosystem exposure to potential warming and acidification stress during the twenty-first century (e.g. Tittensor et al., 2010; Levin and Le Bris, 2015).

Table 5Global mean projected changes in benthic ocean impact drivers in CMIP6. Global mean anomalies of bottom-water temperature (C), pH, and dissolved O2 concentration for the CMIP6 SSPs. Anomalies are 2080–2099 mean values relative to the 1870–1899 baseline period. Uncertainty estimates are the inter-model SD.

Download Print Version | Download XLSX

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f07

Figure 7Global mean projections of benthic ocean impact drivers. CMIP6 global mean projections of benthic (a) temperature (C), (b) pH, and (c) dissolved O2 concentration (mmol m−3). Values are anomalies relative to the 1870–1899 reference period. Mean anomalies for the historical and SSP simulations are shown as solid lines, with shading representing the inter-model SD. The model ensemble size for each scenario is given in parentheses.

Download

The magnitude of projected changes in bottom waters is less than in surface and upper-ocean waters, while bottom-water uncertainties for a given scenario are larger (Fig. 7). This contrast is particularly evident for pH projections with the SSPs, whose ranges of uncertainty fully separate before 2050 in the surface ocean (Fig. 1) but still have some overlap in 2080 for bottom waters. This relative increase in inter-model uncertainty results from surface ocean chemistry being in equilibrium with the same atmospheric CO2 concentrations for all models. Conversely, benthic pH changes are strongly influenced by ocean circulation, which transports anthropogenic carbon from the upper ocean to the seafloor and is variably impacted by climate change across models (e.g. Gregory et al., 2005; Cheng et al., 2013). The increased uncertainty in pH projections with depth has been previously noted for CMIP5 projections in the North Atlantic (Gehlen et al., 2014) and Arctic Ocean (Steiner et al., 2014). For projected global deoxygenation in bottom waters, inter-model uncertainty is substantially larger than scenario uncertainty in CMIP6. As with projections of subsurface dissolved O2, this larger model uncertainty results from the isolation of bottom waters from the atmosphere. Thus, bottom waters at a given temperature and salinity may deviate substantially from the value that would be determined by their solubility and air–sea equilibrium due to effects from other physical and biogeochemical processes (e.g. Oschlies et al., 2018).

3.8 Regional patterns of benthic ocean change

In bottom waters, the end-of-century spatial distributions of changes in temperature, pH, and dissolved O2 are similar between SSPs (Fig. 8) and in broad agreement with CMIP5 projections (Sweetman et al., 2017). The intensity of warming, acidification, and deoxygenation is generally greater in SSP5-8.5 than SSP1-2.6 in benthic waters above 2000 m, while at greater depths the impact is similar for both SSPs.

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f08

Figure 8Projections of multiple benthic ocean impact drivers under SSP1-2.6 and SSP5-8.5. CMIP6 mean historical climatologies and anomalies in benthic (a–d) temperature (C), (e–h) pH, and (j–l) dissolved O2 concentration (mmol m−3). Anomalies are 2080–2099 mean values relative to 1995–2014. Stippling designates areas of projection robustness. For temperature and pH this is defined as the magnitude of the mean anomaly exceeding the inter-model SD. For O2 this is defined as at least 80 % model sign agreement. Vertical profiles show global mean benthic anomalies in 200 m depth intervals, with error bars denoting the inter-model SD.

The largest projected benthic warming in SSP1-2.6 and SSP5-8.5 occurs in continental shelf waters, the Arctic seas, and the Southern Ocean, where temperature increases can exceed 0.5 C by the end-of-century (2080–2099 average relative to the 1995–2014 baseline). In contrast, for most of the abyssal benthic ocean projected increases in temperature are less than 0.2 C. The characteristic North Atlantic “warming hole” present in projections of the surface ocean (Fig. 2) is also evident in benthic layers above 1000 m, such as the mid-Atlantic ridge (Fig. 8b,c). This represents the only major region of multi-model mean benthic cooling across SSP1-2.6 and SSP5-8.5 with high associated uncertainty. As in the surface ocean, this cooling is likely associated with a slow down of the Atlantic meridional overturning circulation (Drijfhout et al., 2012; Menary and Wood, 2018).

Projected end-of-century acidification is highly limited in most bottom waters. However, in the North Atlantic, Arctic seas, and certain continental shelf waters, pH changes can exceed −0.1 in SSP1-2.6 and −0.2 in SSP5-8.5. For shelf waters, the greater bottom-water pH declines can be the result of coupling between surface waters, which experience large changes in carbonate chemistry, and bottom waters (e.g. through mixing and entrainment), as well as benthic remineralisation of organic matter (Bates et al., 2009). In contrast, enhanced bottom-water acidification in the North Atlantic is associated with deep-water formation and high uptake of anthropogenic carbon (Sabine et al., 2004), which rapidly propagates anomalies in surface ocean chemistry to depth. Bottom-water acidification has been previously projected in the North Atlantic by an ensemble of CMIP5 models under RCP8.5 (Gehlen et al., 2014; Sweetman et al., 2017).

In contrast to temperature and pH, projections of benthic dissolved O2 concentration show changes that are not predominantly confined to shelf waters and specific regions. Most of the global benthic ocean is projected to experience deoxygenation under both SSP1-2.6 and SSP5-8.5, even at depths below 2000 m (Fig. 8). Bottom-water deoxygenation is higher in the Southern Ocean, equatorial Pacific, and North Atlantic, where declines in the multi-model mean can exceed 20 mmol m−3.

It should be noted that Earth system models are not explicitly designed to explore the benthic biogeochemical response to climate change and that certain caveats should be considered. Model spin-up simulations, although longer in CMIP6 than CMIP5 (Séférian et al., 2020), are typically insufficient in length to equilibrate biogeochemical conditions in the deep ocean (Séférian et al., 2016), and therefore contemporaneous pre-industrial control simulations are required to correct biogeochemical drift. Vertical thickness of bottom-ocean layers is also highly variable across the CMIP6 ensemble, although for a given model resolution it is typically highest near the surface and decreases dramatically with depth. As such, the extent to which continental shelves are resolved greatly differs and uncertainties associated with resolution are pronounced in the abyssal ocean. Moreover, the representation of biogeochemical processes associated with ocean sediments and benthic ecosystems is typically absent or highly limited in ESMs.

3.9 Depth of maximum acidification

The depth of maximum end-of-century pH and [H+] change is often below the surface, and it varies regionally in CMIP6 projections (Fig. 9). Although the maximum pH change is usually found in surface waters in the high latitudes and upwelling regions, it is typically located between 200–400 m in subtropical mode and intermediate waters. Because of its log scale, if the change in pH were identical in surface and subsurface waters it would imply a larger absolute change in [H+] in the subsurface, where the mean [H+] is higher. Indeed, a change in pH represents a relative change in [H+], not an absolute change in that quantity. That relationship, combined with higher [H+] at depth, means that the maximum change in [H+] is usually deeper than it is for pH. Furthermore, the spatial distribution of the maximum change in pH and [H+] also differs.

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f09

Figure 9Magnitude and depth of maximum pH and [H+] change under SSP5-8.5. The CMIP6 ensemble mean maximum change in (a) pH and (b) [H+] in 2080–2099 of SSP5-8.5 relative to 1995–2014. The mean depth at which the maximum (c) pH and (d) [H+] change is projected. Stippling designates robustness, as defined by the mean anomaly exceeding the inter-model SD.

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f10

Figure 10Change in the seasonal amplitude of surface ocean [H+] and pH. The CMIP6 historical climatologies and SSP5-8.5 anomalies in peak-to-peak seasonal amplitude of surface ocean (a, b) [H+] and (c, d) pH. Anomalies are calculated from the mean seasonal amplitude in 2080–2099 relative to that in 1995–2014. Stippling designates robustness, as defined by at least 80 % model sign agreement.

Enhanced acidification in subsurface mode and intermediate waters has been observed at time series stations (Dore et al., 2009; Byrne et al., 2010; Bates et al., 2012) and in CMIP5 model projections (Resplandy et al., 2013; Bopp et al., 2013; Watanabe and Kawamiya, 2017). Although observational studies have suggested that this enhancement results from changes in circulation and biological activity (Dore et al., 2009; Byrne et al., 2010), model results indicate that it can be explained by the geochemical effect of rising atmospheric CO2 and the particular carbonate chemistry of these waters (Orr, 2011; Resplandy et al., 2013). Specifically, the enhanced acidification sensitivity in mode and intermediate waters has been attributed to their lower temperatures and their higher ratio of dissolved inorganic carbon to total alkalinity relative to that found in surface waters of the same regions (Orr, 2011; Resplandy et al., 2013).

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f11

Figure 11The seasonal amplitude of sea surface temperature. The CMIP6 multi-model mean peak-to-peak seasonal amplitude of sea surface temperature (C) in (a) 1995–2014 of the historical simulations, (b) 2080–2099 of the SSP5-8.5 simulations, and (c) the change in seasonal amplitude between the two periods. Stippling designates robustness, as defined by at least 80 % model sign agreement.

3.10 Surface ocean seasonality

Changes in the seasonal amplitude of surface ocean temperature, pH, and hydrogen ion concentration ([H+]) were determined after detrending by subtracting a cubic spline fit from the monthly time series in each grid cell and then calculating the annual peak-to-peak amplitude for each year of the detrended dataset, following the approach of Kwiatkowski and Orr (2018). Under SSP5-8.5, the seasonal amplitude of global surface ocean [H+] is projected to increase by +73± 12 % across the CMIP6 ensemble (2080–2099 average relative to 1995–2014; Fig. 10). Concurrently, the seasonal amplitude of global surface ocean pH is projected to decrease by -10±5 %. Increases in the seasonal amplitude of [H+] are robust across all regions but are generally highest in the high latitudes. In contrast, declines in the seasonal amplitude of pH are typically only robust in the low latitudes and mid-latitudes, with inconsistent projections of pH seasonal amplitude projected in the Arctic and Southern Oceans.

The simultaneous amplification of [H+] and attenuation of pH seasonal cycles is consistent with previous assessments of CMIP5 projections, with Kwiatkowski and Orr (2018) showing [H+] seasonal amplification of +81± 16 % and pH seasonal attenuation of -16±7 % under RCP8.5 (2090–2099 anomalies relative to 1990–1999). Although counterintuitive, this results from the log scale of pH, which means that the seasonal amplitude of pH depends not only on the seasonal amplitude of [H+] but also on the inverse of the annual mean [H+]. As the projected increase in annual mean [H+] is usually greater than the corresponding increase in the seasonal amplitude of [H+], the seasonal amplitude of pH declines as a result. Increases in the seasonal cycle of [H+] have been shown to be primarily driven by the geochemical effect of increasing atmospheric CO2. This affects both the seasonal amplitude of the controlling variables dissolved inorganic carbon and alkalinity, as well the sensitivity of [H+] to seasonal changes in temperature, dissolved inorganic carbon, and alkalinity (Kwiatkowski and Orr, 2018). Given the near-linear relationship between [H+] and pCO2 on annual timescales (Orr, 2011), projected increases in the seasonal amplitude of [H+] are in agreement with historical observations (Landschützer et al., 2018) and twenty-first century projections from CMIP5 models (McNeil and Sasse, 2016; Gallego et al., 2018) of increasing pCO2 seasonal amplitude.

The multi-model mean seasonal amplitude of global surface ocean temperature is projected to increase by +0.59±0.21C across SSP5-8.5 (Fig. 11). Over most of the ocean, the seasonal amplitude of sea surface temperature is projected to show limited but robust increases (<+0.5C). However, in the North Atlantic, North Pacific, and Southern Ocean, increases in the seasonal amplitude of SST can exceed +2C, and in the Arctic Ocean, the seasonal amplitude typically increases by >+5C.

The CMIP6 projections of the changing seasonal amplitude of SST under SSP5-8.5 are consistent with previous projections from the CMIP5 models (Carton et al., 2015; Alexander et al., 2018). The limited increases in SST seasonal amplitude for most of the global ocean have been attributed to greater relative shoaling of the mixed-layer depth in summer than in winter (Alexander et al., 2018). However, in the Arctic Ocean the large increase in SST seasonal amplitude is primarily due to the loss of sea ice. The seasonal melting and refreezing of sea ice accounts for approximately half of the present-day seasonal Arctic Ocean net surface heat flux, buffering seasonal variability in Arctic Ocean heat content and SSTs (Serreze et al., 2007; Fig. 11). The loss of this seasonal melting–freezing cycle under high-emissions scenarios such as RCP8.5 has been shown to account for a doubling of seasonal Arctic Ocean heat content variability. Ice loss further amplifies the seasonal cycle of SSTs by increasing the seasonal cycle of net surface heat fluxes. The net downward radiative flux increases in summer as albedo declines, while the net upward radiative flux increases in winter due to greater evaporative and sensible heat loss (Carton et al., 2015).

4 Conclusions

The latest CMIP6 Earth system models consistently project global surface ocean warming and acidification, subsurface deoxygenation, and euphotic-zone nitrate reductions in the twenty-first century. Multi-model mean projections of global net primary production show declines in the twenty-first century, although this is with large inter-model uncertainty. The projected change in these ocean impact drivers is shown to increase with radiative forcing across the SSPs, highlighting the benefit of emissions reductions to upper-ocean ecosystems. The magnitude of projected warming, acidification, and deoxygenation is lower in the benthic ocean, with greater inter-model uncertainty relative to scenario uncertainty. However, the extent of warming and acidification is still limited under lower-emissions scenarios, demonstrating the potential benefits of mitigation to benthic ecosystems.

In addition to changing mean-state conditions, the CMIP6 models also project changes to the seasonal cycles of temperature and carbonate chemistry under the SSPs. The seasonal amplitude of surface ocean acidity ([H+]) nearly doubles over the twenty-first century under SSP5-8.5, with a concurrent reduction in the seasonal amplitude of pH. Over the same period, the seasonal amplitude of temperature is projected to increase, particularly in the Arctic Ocean.

The CMIP6 projections of warming, acidification, deoxygenation, and nutrient reduction are greater than those of previous CMIP5 models under comparable radiative forcing. The enhanced acidification is a consequence of higher atmospheric CO2 concentrations in the SSPs than their RCP analogues. The enhanced warming, however, reflects the greater climate sensitivity of the CMIP6 models. This increased warming results in greater increases in upper-ocean stratification, which contributes to greater reductions in euphotic nitrate and subsurface oxygen concentration. The CMIP6 multi-model mean projections of primary production declines are less than those of previous CMIP5 models under comparable radiative forcing; however, there is a large increase in inter-model uncertainty that requires further assessment.

Projected changes to the mean state and seasonality of biogeochemical ocean conditions are likely to present major challenges to diverse marine ecosystems from the surface ocean to abyssal depths. Potential organism stress is likely to be exacerbated by simultaneous exposure to multiple biogeochemical changes, emphasising the need for extensive emissions reductions.

Appendix A
https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f12

Figure A1The interquartile range associated with CMIP6 projections of upper-ocean impact drivers. The CMIP6 interquartile range (Q3Q1) of projected anomalies in (a, b) sea surface temperature (C), (c, d) surface ocean pH, (e, f) subsurface dissolved O2 concentration (averaged between 100 and 600 m;mmolm-3), (g, h) euphotic-zone NO3- (averaged between 0 and 100 m;mmolm-3), and (i, j) depth-integrated net primary production (gCm-2y-1). The corresponding mean projections are given in Fig. 2.

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f13

Figure A2The interquartile range associated with CMIP6 projections of stratification and mixed-layer depth. The CMIP6 interquartile range (Q3Q1) of projected anomalies in (a, b) the stratification index between 200 m and the surface (kg m−3), (c, d) the stratification index between 1000 m and the surface (kg m−3), and (e, f) the maximum annual mixed-layer depth (m). The corresponding mean projections are given in Fig. 4.

https://www.biogeosciences.net/17/3439/2020/bg-17-3439-2020-f14

Figure A3The interquartile range associated with CMIP6 projections of benthic impact drivers. The CMIP6 interquartile range (Q3Q1) of projected anomalies in benthic (a, b) temperature (C), (c, d) pH, and (e, f) dissolved O2 concentration (mmol m−3). The corresponding mean projections are given in Fig. 8.

Data availability

The Earth system model output used in this study is available via the Earth System Grid Federation (https://esgf-node.ipsl.upmc.fr/projects/esgf-ipsl/).

Author contributions

LK and LB conceived and designed this study. LK, LB, and OT processed model outputs and performed the analysis. All authors contributed to the ocean biogeochemistry development of the CMIP6 ESMs and/or the manuscript text.

Competing interests

The authors declare that they have no conflict of interest.

Disclaimer

This article reflects only the authors' views; the funding agencies and their executive agencies are not responsible for any use that may be made of the information that the article contains.

Acknowledgements

We thank Hyung-Gyu Lim for providing an early review of our analysis. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP. For CMIP ,the US Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provided coordinating support and led the development of software infrastructure in partnership with the Global Organisation for Earth System Science Portals. The authors also thank the IPSL modelling group for the software infrastructure, which facilitated CMIP analysis. Lester Kwiatkowski and Laurent Bopp thank the 4C project under grant agreement no. 821003, the MTES/FRB Acidoscope project, and the ENS-Chanel research chair. Roland Séférian thanks the TRIATLAS project under grant agreement no. 817578. Laurent Bopp, Jörg Schwinger, Jerry Tjiputra, Yeray Santana-Falcón, and Roland Séférian thank the COMFORT project under grant agreement no. 820989. Michio Watanabe and Akitomo Yamamoto thank the Integrated Research Programme for Advancing Climate Models (TOUGOU) under grant agreement no. JPMXD0717935715 from the Ministry of Education, Culture, Sports, Science, and Technology (MEXT), Japan. Nicole S. Lovenduski is grateful for support from the US National Science Foundation (OCE-1752724, OCE-1558225). Jörg Schwinger and Jerry Tjiputra acknowledge support from the Research Council of Norway through projects INES (no. 270061), KeyClim (no. 295046), and COLUMBIA (no. 275268), as well as support from the Bjerknes Centre for Climate Research through the project LOES. Alessandro Tagliabue received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 724289).

Financial support

This research has been supported by the Horizon 2020 (CRESCENDO – 641816, COMFORT – 820989, and TRIATLAS – 817578), the National Science Foundation (grant nos. OCE-1752724 and OCE-1558225), and the Agence Nationale de la Recherche (grant no. ANR-18-ERC2-0001-01).

Review statement

This paper was edited by Fortunat Joos and reviewed by two anonymous referees.

References

Albright, R., Caldeira, L., Hosfelt, J., Kwiatkowski, L., Maclaren, J. K., Mason, B. M., Nebuchina, Y., Ninokawa, A., Pongratz, J., Ricke, K. L., Rivlin, T., Schneider, K., Sesboüé, M., Shamberger, K., Silverman, J., Wolfe, K., Zhu, K., and Caldeira, K.: Reversal of ocean acidification enhances net coral reef calcification, Nature, 531, 362–365, https://doi.org/10.1038/nature17155, 2016. 

Alexander, M. A., Scott, J. D., Friedland, K. D., Mills, K. E., Nye, J. A., Pershing, A. J., and Thomas, A. C.: Projected sea surface temperatures over the 21st century: Changes in the mean, variability and extremes for large marine ecosystem regions of Northern Oceans, Elem. Sci. Anth., 6, 9, https://doi.org/10.1525/elementa.191, 2018. 

Andrews, O. D., Bindoff, N. L., Halloran, P. R., Ilyina, T., and Le Quéré, C.: Detecting an external influence on recent changes in oceanic oxygen using an optimal fingerprinting method, Biogeosciences, 10, 1799–1813, https://doi.org/10.5194/bg-10-1799-2013, 2013. 

Andrews, T., Gregory, J. M., Webb, M. J., and Taylor, K. E.: Forcing, feedbacks and climate sensitivity in CMIP5 coupled atmosphere-ocean climate models, Geophys. Res. Lett., 39, L09712, https://doi.org/10.1029/2012GL051607, 2012. 

Bates, N. R., Mathis, J. T., and Cooper, L. W.: Ocean acidification and biologically induced seasonality of carbonate mineral saturation states in the western Arctic Ocean, J. Geophys. Res.-Oceans, 114, C11009, https://doi.org/10.1029/2008JC004862, 2009. 

Bates, N. R., Best, M. H. P., Neely, K., Garley, R., Dickson, A. G., and Johnson, R. J.: Detecting anthropogenic carbon dioxide uptake and ocean acidification in the North Atlantic Ocean, Biogeosciences, 9, 2509–2522, https://doi.org/10.5194/bg-9-2509-2012, 2012. 

Behrenfeld, M. J., Randerson, J. T., McClain, C. R., Feldman, G. C., Los, S. O., Tucker, C. J., Falkowski, P. G., Field, C. B., Frouin, R., Esaias, W. E., Kolber, D. D., and Pollack, N. H.: Biospheric Primary Production During an ENSO Transition, Science, 291, 2594–2597, https://doi.org/10.1126/science.1055071, 2001. 

Behrenfeld, M. J., O'Malley, R. T., Siegel, D. A., McClain, C. R., Sarmiento, J. L., Feldman, G. C., Milligan, A. J., Falkowski, P. G., Letelier, R. M., and Boss, E. S.: Climate-driven trends in contemporary ocean productivity, Nature, 444, 752–755, https://doi.org/10.1038/nature05317, 2006. 

Bekryaev, R. V., Polyakov, I. V., and Alexeev, V. A.: Role of Polar Amplification in Long-Term Surface Air Temperature Variations and Modern Arctic Warming, J. Climate, 23, 3888–3906, https://doi.org/10.1175/2010JCLI3297.1, 2010. 

Bentsen, M., Bethke, I., Debernard, J. B., Iversen, T., Kirkevåg, A., Seland, Ø., Drange, H., Roelandt, C., Seierstad, I. A., Hoose, C., and Kristjánsson, J. E.: The Norwegian Earth System Model, NorESM1-M – Part 1: Description and basic evaluation of the physical climate, Geosci. Model Dev., 6, 687–720, https://doi.org/10.5194/gmd-6-687-2013, 2013. 

Bindoff, N. L., Cheung, W. W. L., Kairo, J. G., Arístegui, J., Guinder, V. A., Hallberg, R., Hilmi, N., Jiao, N., Karim, M. S., Levin, L., O'Donoghue, S., Purca Cuicapusa, S. R., Rinkevich, B., Suga, T., Tagliabue, A., and Williamson, P.: Changing Ocean, Marine Ecosystems, and Dependent Communities. In: IPCC Special Report on the Ocean and Cryosphere in a Changing Climate, edited by: Pörtner, H.-O., Roberts, D. C., Masson-Delmotte, V., Zhai, P., Tignor, M., Poloczanska, E., Mintenbeck, K., Alegría, A., Nicolai, M., Okem, A., Petzold, J., Rama, B., Weyer, N. M., in press., 2020. 

Bindoff, N. L., Willebrand, J., Artale, V., Cazenave, A., Gregory, J. M., Gulev, S., Hanawa, K., Le Quere, C., Levitus, S., Nojiri, Y., Shum, C. K., Talley, L. D., Unnikrishnan, A. S., Josey, S. A., Tamisiea, M., Tsimplis, M., and Woodworth, P.: Observations: oceanic climate change and sea level, in Climate change 2007: the physical science basis. Contribution of Working Group I, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., Miller, H. L., Cambridge University Press, Cambridge, pp. 385–428, 2007. 

Bopp, L., Monfray, P., Aumont, O., Dufresne, J.-L., Treut, H. L., Madec, G., Terray, L., and Orr, J. C.: Potential impact of climate change on marine export production, Global Biogeochem. Cy., 15, 81–99, https://doi.org/10.1029/1999GB001256, 2001. 

Bopp, L., Resplandy, L., Orr, J. C., Doney, S. C., Dunne, J. P., Gehlen, M., Halloran, P., Heinze, C., Ilyina, T., Séférian, R., Tjiputra, J., and Vichi, M.: Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models, Biogeosciences, 10, 6225–6245, https://doi.org/10.5194/bg-10-6225-2013, 2013. 

Bopp, L., Resplandy, L., Untersee, A., Le Mezo, P., and Kageyama, M.: Ocean (de)oxygenation from the Last Glacial Maximum to the twenty-first century: insights from Earth System models, Philos. Trans. A Math. Phys. Eng. Sci., 375, 20160323, https://doi.org/10.1098/rsta.2016.0323, 2017. 

Boucher, O., Denvil, S., Caubel, Arnaud; F., and Marie A.: IPSL IPSL-CM6A-LR model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.5195, 2018. 

Boucher, O., Denvil, S., Caubel, Arnaud; F., and Marie A.: IPSL IPSL-CM6A-LR model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.5264, 2019. 

Boucher, O., Servonnat, J., Albright, A. L., Aumont, O., Balkanski, Y., Bastrikov, V., Bekki, S., Bonnet, R., Bony, S., Bopp, L., Braconnot, P., Brockmann, P., Cadule, P., Caubel, A., Cheruy, F., Codron, F., Cozic, A., Cugnet, D., D’Andrea, F., Davini, P., Lavergne, C. de, Denvil, S., Deshayes, J., Devilliers, M., Ducharne, A., Dufresne, J.-L., Dupont, E., Éthé, C., Fairhead, L., Falletti, L., Flavoni, S., Foujols, M.-A., Gardoll, S., Gastineau, G., Ghattas, J., Grandpeix, J.-Y., Guenet, B., Guez, L., Guilyardi, É., Guimberteau, M., Hauglustaine, D., Hourdin, F., Idelkadi, A., Joussaume, S., Kageyama, M., Khodri, M., Krinner, G., Lebas, N., Levavasseur, G., Lévy, C., Li, L., Lott, F., Lurton, T., Luyssaert, S., Madec, G., Madeleine, J.-B., Maignan, F., Marchand, M., Marti, O., Mellul, L., Meurdesoif, Y., Mignot, J., Musat, I., Ottlé, C., Peylin, P., Planton, Y., Polcher, J., Rio, C., Rochetin, N., Rousset, C., Sepulchre, P., Sima, A., Swingedouw, D., Thiéblemont, R., Traore, A. K., Vancoppenolle, M., Vial, J., Vialard, J., Viovy, N. and Vuichard, N.: Presentation and evaluation of the IPSL‐CM6A‐LR climate model, J. Adv. Model. Earth Sy., 12, e2019MS002010, https://doi.org/10.1029/2019MS002010, 2020. 

Breitburg, D., Levin, L. A., Oschlies, A., Grégoire, M., Chavez, F. P., Conley, D. J., Garçon, V., Gilbert, D., Gutiérrez, D., Isensee, K., Jacinto, G. S., Limburg, K. E., Montes, I., Naqvi, S. W. A., Pitcher, G. C., Rabalais, N. N., Roman, M. R., Rose, K. A., Seibel, B. A., Telszewski, M., Yasuhara, M., and Zhang, J.: Declining oxygen in the global ocean and coastal waters, Science, 359, https://doi.org/10.1126/science.aam7240, 2018. 

Busecke, J. J. M., Resplandy, L., and Dunne, J. P.: The Equatorial Undercurrent and the Oxygen Minimum Zone in the Pacific, Geophys. Res. Lett., 46, 6716–6725, https://doi.org/10.1029/2019GL082692, 2019. 

Byrne, R. H., Mecking, S., Feely, R. A., and Liu, X.: Direct observations of basin-wide acidification of the North Pacific Ocean, Geophys. Res. Lett., 37, L02601, https://doi.org/10.1029/2009GL040999, 2010. 

Cabré, A., Marinov, I., and Leung, S.: Consistent global responses of marine ecosystems to future climate change across the IPCC AR5 earth system models, Clim. Dynam., 45, 1253–1280, https://doi.org/10.1007/s00382-014-2374-3, 2014. 

Cabré, A., Marinov, I., Bernardello, R., and Bianchi, D.: Oxygen minimum zones in the tropical Pacific across CMIP5 models: mean state differences and climate change trends, Biogeosciences, 12, 5429–5454, https://doi.org/10.5194/bg-12-5429-2015, 2015. 

Cagnazzo, C., Manzini, E., Fogli, P. G., Vichi, M., and Davini, P.: Role of stratospheric dynamics in the ozone–carbon connection in the Southern Hemisphere, Clim. Dynam., 41, 3039–3054, https://doi.org/10.1007/s00382-013-1745-5, 2013. 

Capotondi, A., Alexander, M. A., Bond, N. A., Curchitser, E. N., and Scott, J. D.: Enhanced upper ocean stratification with climate change in the CMIP3 models, J. Geophys. Res.-Oceans, 117, C04031, https://doi.org/10.1029/2011JC007409, 2012. 

Carton, J. A., Ding, Y., and Arrigo, K. R.: The seasonal cycle of the Arctic Ocean under climate change, Geophys. Res. Lett., 42, 7681–7686, https://doi.org/10.1002/2015GL064514, 2015. 

Cheng, L., Abraham, J., Hausfather, Z., and Trenberth, K. E.: How fast are the oceans warming?, Science, 363, 128–129, https://doi.org/10.1126/science.aav7619, 2019. 

Cheng, W., Chiang, J. C. H., and Zhang, D.: Atlantic Meridional Overturning Circulation (AMOC) in CMIP5 Models: RCP and Historical Simulations, J. Climate, 26, 7187–7197, https://doi.org/10.1175/JCLI-D-12-00496.1, 2013. 

Chust, G., Allen, J. I., Bopp, L., Schrum, C., Holt, J., Tsiaras, K., Zavatarelli, M., Chifflet, M., Cannaby, H., Dadou, I., Daewel, U., Wakelin, S. L., Machu, E., Pushpadas, D., Butenschon, M., Artioli, Y., Petihakis, G., Smith, C., Garçon, V., Goubanova, K., Le Vu, B., Fach, B. A., Salihoglu, B., Clementi, E., and Irigoien, X.: Biomass changes and trophic amplification of plankton in a warmer ocean, Glob. Change Biol., 20, 2124–2139, https://doi.org/10.1111/gcb.12562, 2014. 

Cocco, V., Joos, F., Steinacher, M., Frölicher, T. L., Bopp, L., Dunne, J., Gehlen, M., Heinze, C., Orr, J., Oschlies, A., Schneider, B., Segschneider, J., and Tjiputra, J.: Oxygen and indicators of stress for marine life in multi-model global warming projections, Biogeosciences, 10, 1849–1868, https://doi.org/10.5194/bg-10-1849-2013, 2013. 

Collins, W. J., Bellouin, N., Doutriaux-Boucher, M., Gedney, N., Halloran, P., Hinton, T., Hughes, J., Jones, C. D., Joshi, M., Liddicoat, S., Martin, G., O'Connor, F., Rae, J., Senior, C., Sitch, S., Totterdell, I., Wiltshire, A., and Woodward, S.: Development and evaluation of an Earth-System model – HadGEM2, Geosci. Model Dev., 4, 1051–1075, https://doi.org/10.5194/gmd-4-1051-2011, 2011. 

Danabasoglu, G.: NCAR CESM2 model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.7522, 2019a. 

Danabasoglu, G.: NCAR CESM2 model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.7753, 2019b. 

Danabasoglu, G.: NCAR CESM2-WACCM model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.10028, 2019c. 

Danabasoglu, G.: NCAR CESM2-WACCM model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.10101, 2019d. 

Dave, A. C. and Lozier, M. S.: Examining the global record of interannual variability in stratification and marine productivity in the low-latitude and mid-latitude ocean, J. Geophys. Res.-Oceans, 118, 3114–3127, https://doi.org/10.1002/jgrc.20224, 2013. 

Deutsch, C., Brix, H., Ito, T., Frenzel, H., and Thompson, L.: Climate-Forced Variability of Ocean Hypoxia, Science, 333, 336–339, https://doi.org/10.1126/science.1202422, 2011. 

Deutsch, C., Ferrel, A., Seibel, B., Pörtner, H.-O., and Huey, R. B.: Climate change tightens a metabolic constraint on marine habitats, Science, 348, 1132–1135, https://doi.org/10.1126/science.aaa1605, 2015. 

Doney, S. C.: Plankton in a warmer world, Nature, 444, 695–696, https://doi.org/10.1038/444695a, 2006. 

Doney, S. C., Fabry, V. J., Feely, R. A., and Kleypas, J. A.: Ocean Acidification: The Other CO2 Problem, Annu. Rev. Marine Sci., 1, 169–192, https://doi.org/10.1146/annurev.marine.010908.163834, 2009. 

Dore, J. E., Lukas, R., Sadler, D. W., Church, M. J., and Karl, D. M.: Physical and biogeochemical modulation of ocean acidification in the central North Pacific, Proc. Natl. Acad. Sci. USA, 106, 12235–12240, https://doi.org/10.1073/pnas.0906044106, 2009. 

Drijfhout, S., van Oldenborgh, G. J., and Cimatoribus, A.: Is a Decline of AMOC Causing the Warming Hole above the North Atlantic in Observed and Modeled Warming Patterns?, J. Climate, 25, 8373–8379, https://doi.org/10.1175/JCLI-D-12-00490.1, 2012. 

Dufresne, J.-L., Foujols, M.-A., Denvil, S., Caubel, A., Marti, O., Aumont, O., Balkanski, Y., Bekki, S., Bellenger, H., Benshila, R., Bony, S., Bopp, L., Braconnot, P., Brockmann, P., Cadule, P., Cheruy, F., Codron, F., Cozic, A., Cugnet, D., Noblet, N., Duvel, J.-P., Ethé, C., Fairhead, L., Fichefet, T., Flavoni, S., Friedlingstein, P., Grandpeix, J.-Y., Guez, L., Guilyardi, E., Hauglustaine, D., Hourdin, F., Idelkadi, A., Ghattas, J., Joussaume, S., Kageyama, M., Krinner, G., Labetoulle, S., Lahellec, A., Lefebvre, M.-P., Lefevre, F., Levy, C., Li, Z. X., Lloyd, J., Lott, F., Madec, G., Mancip, M., Marchand, M., Masson, S., Meurdesoif, Y., Mignot, J., Musat, I., Parouty, S., Polcher, J., Rio, C., Schulz, M., Swingedouw, D., Szopa, S., Talandier, C., Terray, P., Viovy, N., and Vuichard, N.: Climate change projections using the IPSL-CM5 Earth System Model: from CMIP3 to CMIP5, Clim. Dynam., 40, 2123–2165, https://doi.org/10.1007/s00382-012-1636-1, 2013. 

Dunne, J. P., John, J. G., Adcroft, A. J., Griffies, S. M., Hallberg, R. W., Shevliakova, E., Stouffer, R. J., Cooke, W., Dunne, K. A., Harrison, M. J., Krasting, J. P., Malyshev, S. L., Milly, P. C. D., Phillipps, P. J., Sentman, L. T., Samuels, B. L., Spelman, M. J., Winton, M., Wittenberg, A. T., and Zadeh, N.: GFDL's ESM2 Global Coupled Climate–Carbon Earth System Models. Part I: Physical Formulation and Baseline Simulation Characteristics, J. Climate, 25, 6646–6665, https://doi.org/10.1175/JCLI-D-11-00560.1, 2012. 

Dunne, J. P., Horowitz, L. W. , Adcroft, A. J., Ginoux, P. , Held, I. M., John, J. G., Krasting, J. P., Malyshev, S., Naik, V., Paulot, F., Shevliakova, E., Stock, C. A., Zadeh, N., Balaji, V., Blanton, C., Dunne, K. A., Dupuis, C., Durachta, J., Dussin, R., Gauthier, P. P. G., Griffies, S. M., Guo, H., Hallberg, R. W., Harrison, M., He, J., Hurlin, W., McHugh, C., Menzel, R., Milly, P. C. D., Nikonov, S., Paynter, D. J., Ploshay, J., Radhakrishnan, A., Rand, K., Reichl, B. G., Robinson, T., Schwarzkopf, D. M., Sentman, L. A., Underwood, S., Vahlenkamp, H., Winton, M., Wittenberg, A. T., Wyman, B., Zeng, Y., and Zhao, M.: The GFDL Earth System Model version 4.1 (GFDL-ESM4.1): Model description and simulation characteristics, J. Adv. Model. Earth Syst., 2019MS002008, 2020a. 

Dunne, J. P., Bociu, I., Bronselaer, B., Guo, H., John, J. G., Krasting, J. P., Stock, C. A., Winton, M., and Zadeh, N.: Simple Global Ocean Biogeochemistry with Light, Iron, Nutrients and Gas version 2 (BLINGv2): Model description and simulation characteristics in GFDL's CM4.0, J. Adv. Model. Earth Syst., 2019MS002008, 2020b. 

Elsworth, G. W., Lovenduski, N. S., McKinnon, K. A., Krumhardt, K. M., and Brady, R. X.: Finding the fingerprint of anthropogenic climate change in marine phytoplankton abundance, Curr. Clim. Change Rep., 6, https://doi.org/10.1007/s40641-020-00156-w, 2020. 

Eppley, R. W.: Temperature and phytoplankton growth in the sea, Fish Bull., 70, 1063–85, 1972. 

Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958, https://doi.org/10.5194/gmd-9-1937-2016, 2016. 

Freeman, N. M., Lovenduski, N. S., Munro, D. R., Krumhardt, K. M., Lindsay, K., Long, M. C., and Maclennan, M.: The variable and changing Southern Ocean Silicate Front: Insights from the CESM Large Ensemble, Global Biogeochem. Cy., 32, 752–768, https://doi.org/10.1029/2017GB005816, 2018. 

Forster, P. M., Maycock, A. C., McKenna, C. M., and Smith, C. J.: Latest climate models confirm need for urgent mitigation, Nat. Clim. Change, 1–4, https://doi.org/10.1038/s41558-019-0660-0, 2019. 

Frölicher, T. L., Sarmiento, J. L., Paynter, D. J., Dunne, J. P., Krasting, J. P., and Winton, M.: Dominance of the Southern Ocean in Anthropogenic Carbon and Heat Uptake in CMIP5 Models, J. Climate, 28, 862–886, https://doi.org/10.1175/JCLI-D-14-00117.1, 2014. 

Frölicher, T. L., Rodgers, K. B., Stock, C. A., and Cheung, W. W. L.: Sources of uncertainties in 21st century projections of potential ocean ecosystem stressors, Global Biogeochem. Cy., 30, 2015GB005338, https://doi.org/10.1002/2015GB005338, 2016. 

Frölicher, T. L., Fischer, E. M., and Gruber, N.: Marine heatwaves under global warming, Nature, 560, 360–364, https://doi.org/10.1038/s41586-018-0383-9, 2018. 

Fu, W., Randerson, J. T., and Moore, J. K.: Climate change impacts on net primary production (NPP) and export production (EP) regulated by increasing stratification and phytoplankton community structure in the CMIP5 models, Biogeosciences, 13, 5151–5170, https://doi.org/10.5194/bg-13-5151-2016, 2016. 

Hartmann, D. L., Tank, A. M. K., Rusticucci, M., Alexander, L. V., Brönnimann, S., Charabi, Y. A. R., Dentener, F. J., Dlugokencky, E. J., Easterling, D. R., Kaplan, A., and Soden, B. J.: Observations: atmosphere and surface, in: Climate Change 2013 – the physical science basis: Working group I contribution to the fifth assessment report of the intergovernmental panel on climate change, Cambridge University Press, Cambridge, pp. 159–254, 2013. 

Gallego, M. A., Timmermann, A., Friedrich, T., and Zeebe, R. E.: Drivers of future seasonal cycle changes in oceanic pCO2, Biogeosciences, 15, 5315–5327, https://doi.org/10.5194/bg-15-5315-2018, 2018. 

Garrabou, J., Coma, R., Bensoussan, N., Bally, M., Chevaldonné, P., Cigliano, M., Diaz, D., Harmelin, J. G., Gambi, M. C., Kersting, D. K., Ledoux, J. B., Lejeusne, C., Linares, C., Marschal, C., Pérez, T., Ribes, M., Romano, J. C., Serrano, E., Teixido, N., Torrents, O., Zabala, M., Zuberer, F., and Cerrano, C.: Mass mortality in Northwestern Mediterranean rocky benthic communities: effects of the 2003 heat wave, Glob. Change Biol., 15, 1090–1103, https://doi.org/10.1111/j.1365-2486.2008.01823.x, 2009. 

Gattuso, J.-P. and Buddemeier, R. W.: Ocean biogeochemistry: Calcification and CO2, Nature, 407, 311–314, 2000. 

Gattuso, J.-P. and Hansson, L.: Ocean Acidification, OUP, Oxford, 2011. 

Gehlen, M., Séférian, R., Jones, D. O. B., Roy, T., Roth, R., Barry, J., Bopp, L., Doney, S. C., Dunne, J. P., Heinze, C., Joos, F., Orr, J. C., Resplandy, L., Segschneider, J., and Tjiputra, J.: Projected pH reductions by 2100 might put deep North Atlantic biodiversity at risk, Biogeosciences, 11, 6955–6967, https://doi.org/10.5194/bg-11-6955-2014, 2014. 

Gent, P. R., Danabasoglu, G., Donner, L. J., Holland, M. M., Hunke, E. C., Jayne, S. R., Lawrence, D. M., Neale, R. B., Rasch, P. J., Vertenstein, M., Worley, P. H., Yang, Z.-L., and Zhang, M.: The Community Climate System Model Version 4, J. Climate, 24, 4973–4991, https://doi.org/10.1175/2011JCLI4083.1, 2011. 

Giorgetta, M. A., Jungclaus, J., Reick, C. H., Legutke, S., Bader, J., Böttinger, M., Brovkin, V., Crueger, T., Esch, M., Fieg, K., Glushak, K., Gayler, V., Haak, H., Hollweg, H.-D., Ilyina, T., Kinne, S., Kornblueh, L., Matei, D., Mauritsen, T., Mikolajewicz, U., Mueller, W., Notz, D., Pithan, F., Raddatz, T., Rast, S., Redler, R., Roeckner, E., Schmidt, H., Schnur, R., Segschneider, J., Six, K. D., Stockhause, M., Timmreck, C., Wegner, J., Widmann, H., Wieners, K.-H., Claussen, M., Marotzke, J., and Stevens, B.: Climate and carbon cycle changes from 1850 to 2100 in MPI-ESM simulations for the Coupled Model Intercomparison Project phase 5, J. Adv. Model. Earth Syst., 5, 572–597, https://doi.org/10.1002/jame.20038, 2013. 

Good, P., Sellar, A., Tang, Y., Rumbold, S., Ellis, R., Kelley, D., and Kuhlbrodt, T.: MOHC UKESM1.0-LL model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.6405, 2019. 

Gregory, B., Christophe, L., and Martin, E.: Rapid biogeographical plankton shifts in the North Atlantic Ocean, Glob. Change Biol., 15, 1790–1803, https://doi.org/10.1111/j.1365-2486.2009.01848.x, 2009. 

Gregory, J. M., Dixon, K. W., Stouffer, R. J., Weaver, A. J., Driesschaert, E., Eby, M., Fichefet, T., Hasumi, H., Hu, A., Jungclaus, J. H., Kamenkovich, I. V., Levermann, A., Montoya, M., Murakami, S., Nawrath, S., Oka, A., Sokolov, A. P., and Thorpe, R. B.: A model intercomparison of changes in the Atlantic thermohaline circulation in response to increasing atmospheric CO2 concentration, Geophys. Res. Lett., 32, L12703, https://doi.org/10.1029/2005GL023209, 2005. 

Gruber, N.: The Dynamics of the Marine Nitrogen Cycle and its Influence on Atmospheric CO2 Variations, in The Ocean Carbon Cycle and Climate, edited by: Follows, M., Oguz, T., Springer Netherlands, Dordrecht, pp. 97–148, 2004. 

Gruber, N., Clement, D., Carter, B. R., Feely, R. A., van Heuven, S., Hoppema, M., Ishii, M., Key, R. M., Kozyr, A., Lauvset, S. K., Lo Monaco, C., Mathis, J. T., Murata, A., Olsen, A., Perez, F. F., Sabine, C. L., Tanhua, T., and Wanninkhof, R.: The oceanic sink for anthropogenic CO2 from 1994 to 2007, Science, 363, 1193–1199, https://doi.org/10.1126/science.aau5153, 2019. 

Gunderson, A. R., Armstrong, E. J., and Stillman, J. H.: Multiple Stressors in a Changing World: The Need for an Improved Perspective on Physiological Responses to the Dynamic Marine Environment, Annu. Rev. Marine Sci., 8, 357–378, https://doi.org/10.1146/annurev-marine-122414-033953, 2016. 

Guo, H., John, J. G., Blanton, C., McHugh, C., Nikonov, S., Radhakrishnan, A., Rand, K., Zadeh, N. T., Balaji, V., Durachta, J., Dupuis, C., Menzel, R., Robinson, T., Underwood, S., Vahlenkamp, H., Bushuk, M., Dunne, K. A., Dussin, R., Gauthier, P. P., Ginoux, P., Griffies, S. M., Hallberg, R., Harrison, M., Hurlin, W., Lin, P., Malyshev, S., Naik, V., Paulot, F., Paynter, D. J., Ploshay, J., Reichl, B., Schwarzkopf, D. M., Seman, C. J., Shao, A., Silvers, L., Wyman, B., Yan, X., Zeng, Y., Adcroft, A., Dunne, J. P., Held, I. M., Krasting, J. P., Horowitz, L. W., Milly, P. C. D., Shevliakova, E., Winton, M., Zhao, M., and Zhang, R.: NOAA-GFDL GFDL-CM4 model output prepared for CMIP6 CMIP. Version 20180701. Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.1402, 2018a. 

Guo, H., John, J. G., Blanton, C., McHugh, C., Nikonov, S., Radhakrishnan, A., Rand, K., Zadeh, N. T., Balaji, V., Durachta, J., Dupuis, C., Menzel, R., Robinson, T., Underwood, S., Vahlenkamp, H., Dunne, K. A., Gauthier, P. P., Ginoux, P., Griffies, S. M., Hallberg, R., Harrison, M., Hurlin, W., Lin, P., Malyshev, S., Naik, V., Paulot, F., Paynter, D. J., Ploshay, J., Schwarzkopf, D. M., Seman, C. J., Shao, A., Silvers, L., Wyman, B., Yan, X., Zeng, Y., Adcroft, A., Dunne, J. P., Held, I. M., Krasting, J. P., Horowitz, L. W., Milly, C., Shevliakova, E., Winton, M., Zhao, M., and Zhang, R.: NOAA-GFDL GFDL-CM4 model output prepared for CMIP6 ScenarioMIP. Version 20180701. Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.9242, 2018b. 

Hajima, T., Abe, M., Arakawa, O., Suzuki, T., Komuro, Y., Ogura, T., Ogochi, K., Watanabe, M., Yamamoto, A., Tatebe, H., Noguchi, M. A., Ohgaito, R., Ito, A., Yamazaki, D., Ito, A., Takata, K., Watanabe, S., Kawamiya, M., and Tachiiri, K.: MIROC MIROC-ES2L model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.5602, 2019. 

Hajima, T., Watanabe, M., Yamamoto, A., Tatebe, H., Noguchi, M. A., Abe, M., Ohgaito, R., Ito, A., Yamazaki, D., Okajima, H., Ito, A., Takata, K., Ogochi, K., Watanabe, S., and Kawamiya, M.: Development of the MIROC-ES2L Earth system model and the evaluation of biogeochemical processes and feedbacks, Geosci. Model Dev., 13, 2197–2244, https://doi.org/10.5194/gmd-13-2197-2020, 2020. 

Harvey, B. P., Gwynn-Jones, D., and Moore, P. J.: Meta-analysis reveals complex marine biological responses to the interactive effects of ocean acidification and warming, Ecol. Evol., 3, 1016–1030, https://doi.org/10.1002/ece3.516, 2013. 

Held, I. M., Guo, H., Adcroft, A., Dunne, J. P., Horowitz, L. W., Krasting, J., Shevliakova, E., Winton, M., Zhao, M., Bushuk, M., Wittenberg, A. T., Wyman, B., Xiang, B., Zhang, R., Anderson, W., Balaji, V., Donner, L., Dunne, K., Durachta, J., Gauthier, P. P. G., Ginoux, P., Golaz, J.-C., Griffies, S. M., Hallberg, R., Harris, L., Harrison, M., Hurlin, W., John, J., Lin, P., Lin, S.-J., Malyshev, S., Menzel, R., Milly, P. C. D., Ming, Y., Naik, V., Paynter, D., Paulot, F., Rammaswamy, V., Reichl, B., Robinson, T., Rosati, A., Seman, C., Silvers, L. G., Underwood, S., and Zadeh, N.: Structure and Performance of GFDL's CM4.0 Climate Model, J. Adv. Model. Earth Syst., 11, 3691–3727, https://doi.org/10.1029/2019MS001829, 2019. 

Helm, K. P., Bindoff, N. L., and Church, J. A.: Observed decreases in oxygen content of the global ocean, Geophys. Res. Lett., 38, L23602, https://doi.org/10.1029/2011GL049513, 2011. 

Hobday, A. J., Alexander, L. V., Perkins, S. E., Smale, D. A., Straub, S. C., Oliver, E. C. J., Benthuysen, J. A., Burrows, M. T., Donat, M. G., Feng, M., Holbrook, N. J., Moore, P. J., Scannell, H. A., Sen Gupta, A., and Wernberg, T.: A hierarchical approach to defining marine heatwaves, Progress Oceanogr., 141, 227–238, https://doi.org/10.1016/j.pocean.2015.12.014, 2016. 

Ito, T., Minobe, S., Long, M. C., and Deutsch, C.: Upper ocean O2 trends: 1958–2015, Geophys. Res. Lett., 44, 4214–4223, https://doi.org/10.1002/2017GL073613, 2017. 

John, J. G., Blanton, C., McHugh, C., Radhakrishnan, A., Rand, K., Vahlenkamp, H., Wilson, C., Zadeh, N. T., Gauthier, P. P., Dunne, J.P., Dussin, R., Horowitz, L. W., Lin, P., Malyshev, S., Naik, V., , Ploshay, J., Silvers, L., Stock, C., Winton, M., and Zeng, Y.: NOAA-GFDL GFDL-ESM4 model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.1414, 2018. 

Jungclaus, J., Bittner, M., Wieners, K.-H., Wachsmann, F., Schupfner, M., Legutke, S., Giorgetta, M., Reick, C., Gayler, V., Haak, H., de Vrese, P., Raddatz, T., Esch, M., Mauritsen, T., von Storch, J.-S., Behrens, J., Brovkin, V., Claussen, M., Crueger, T., Fast, I., Fiedler, S., Hagemann, S., Hohenegger, C., Jahns, T., Kloster, S., Kinne, S., Lasslop, G., Kornblueh, L., Marotzke, J., Matei, D., Meraner, K., Mikolajewicz, U., Modali, K., Müller, W., Nabel, J., Notz, D., Peters, K., Pincus, R., Pohlmann, H., Pongratz, J., Rast, S., Schmidt, H., Schnur, R., Schulzweida, U., Six, K., Stevens, B., Voigt, A., and Roeckner, E.: MPI-M MPI-ESM1.2-HR model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.6594, 2019. 

Keeling, R. F., Körtzinger, A., and Gruber, N.: Ocean Deoxygenation in a Warming World, Annu. Rev. Marine Sci., 2, 199–229, https://doi.org/10.1146/annurev.marine.010908.163855, 2010. 

Khatiwala, S., Primeau, F., and Hall, T.: Reconstruction of the history of anthropogenic CO2 concentrations in the ocean, Nature, 462, 346–349, https://doi.org/10.1038/nature08526, 2009. 

Khatiwala, S., Tanhua, T., Mikaloff Fletcher, S., Gerber, M., Doney, S. C., Graven, H. D., Gruber, N., McKinley, G. A., Murata, A., Ríos, A. F., and Sabine, C. L.: Global ocean storage of anthropogenic carbon, Biogeosciences, 10, 2169–2191, https://doi.org/10.5194/bg-10-2169-2013, 2013. 

Knutti, R., Masson, D., and Gettelman, A.: Climate model genealogy: Generation CMIP5 and how we got there, Geophys. Res. Lett., 40, 1194–1199, 2015. 

Krasting, J. P., John, J. G., Blanton, C., McHugh, C., Nikonov, S., Radhakrishnan, A., Rand, K., Zadeh, N. T., Balaji, V., Durachta, J., Dupuis, C., Menzel, R., Robinson, T., Underwood, S., Vahlenkamp, H., Dunne, K. A., Gauthier, P. P., Ginoux, P., Griffies, S. M., Hallberg, R., Harrison, M., Hurlin, W., Malyshev, S., Naik, V., Paulot, F., Paynter, D. J., Ploshay, J., Schwarzkopf, D. M., Seman, C. J., Silvers, L., Wyman, B., Zeng, Y., Adcroft, A., Dunne, J. P., Dussin, R., Guo, H., He, J., Held, I. M., Horowitz, L. W., Lin, P., Milly, P. C. D., Shevliakova, E., Stock, C., Winton, M., Xie, Y., and Zhao, M.: NOAA-GFDL GFDL-ESM4 model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.1407, 2018. 

Kriegler, E., Bauer, N., Popp, A., Humpenöder, F., Leimbach, M., Strefler, J., Baumstark, L., Bodirsky, B. L., Hilaire, J., Klein, D., Mouratiadou, I., Weindl, I., Bertram, C., Dietrich, J.-P., Luderer, G., Pehl, M., Pietzcker, R., Piontek, F., Lotze-Campen, H., Biewald, A., Bonsch, M., Giannousakis, A., Kreidenweis, U., Müller, C., Rolinski, S., Schultes, A., Schwanitz, J., Stevanovic, M., Calvin, K., Emmerling, J., Fujimori, S., and Edenhofer, O.: Fossil-fueled development (SSP5): An energy and resource intensive scenario for the 21st century, Global Environ. Change, 42, 297–315, https://doi.org/10.1016/j.gloenvcha.2016.05.015, 2017. 

Kroeker, K. J., Kordas, R. L., Crim, R. N., and Singh, G. G.: Meta-analysis reveals negative yet variable effects of ocean acidification on marine organisms, Ecol. Lett., 13, 1419–1434, https://doi.org/10.1111/j.1461-0248.2010.01518.x, 2010. 

Kroeker, K. J., Micheli, F., and Gambi, M. C.: Ocean acidification causes ecosystem shifts via altered competitive interactions, Nat. Clim. Change, 3, 156–159, https://doi.org/10.1038/nclimate1680, 2013. 

Kroeker, K. J., Bell, L. E., Donham, E. M., Hoshijima, U., Lummis, S., Toy, J. A., and Willis-Norton, E.: Ecological change in dynamic environments: Accounting for temporal environmental variability in studies of ocean change biology, Glob. Change Biol., 26, 54–67, https://doi.org/10.1111/gcb.14868, 2020. 

Krumhardt, K. M., Lovenduski, N. S., Long, M. C., and Lindsay, K.: Avoidable impacts of ocean warming on marine primary production: Insights from the CESM ensembles, Global Biogeochem. Cy., 31, 114–133, https://doi.org/10.1002/2016GB005528, 2017. 

Kwiatkowski, L. and Orr, J. C.: Diverging seasonal extremes for ocean acidification during the twenty-first century, Nat. Clim. Change, 8, 141–145, https://doi.org/10.1038/s41558-017-0054-0, 2018. 

Kwiatkowski, L., Gaylord, B., Hill, T., Hosfelt, J., Kroeker, K. J., Nebuchina, Y., Ninokawa, A., Russell, A. D., Rivest, E. B., Sesboüé, M., and Caldeira, K.: Nighttime dissolution in a temperate coastal ocean ecosystem increases under acidification, Sci. Rep., 6, 22984, https://doi.org/10.1038/srep22984, 2016. 

Kwiatkowski, L., Bopp, L., Aumont, O., Ciais, P., Cox, P. M., Laufkötter, C., Li, Y., and Séférian, R.: Emergent constraints on projections of declining primary production in the tropical oceans, Nat. Clim. Change, 7, 355–358, https://doi.org/10.1038/nclimate3265, 2017. 

Kwiatkowski, L., Aumont, O., and Bopp, L.: Consistent trophic amplification of marine biomass declines under climate change, Glob. Change Biol., 25, 218–229, https://doi.org/10.1111/gcb.14468, 2018. 

Landschützer, P., Gruber, N., Bakker, D. C. E., Stemmler, I., and Six, K. D.: Strengthening seasonal marine CO2 variations due to increasing atmospheric CO2, Nat. Clim. Change, 8, 146–150, https://doi.org/10.1038/s41558-017-0057-x, 2018. 

Langenbuch, M., Bock, C., Leibfritz, D., and Pörtner, H. O.: Effects of environmental hypercapnia on animal physiology: A 13C NMR study of protein synthesis rates in the marine invertebrate Sipunculus nudus, Comp. Biochem. Physiol. A, 144, 479–484, https://doi.org/10.1016/j.cbpa.2006.04.017, 2006. 

Laufkötter, C., Vogt, M., Gruber, N., Aita-Noguchi, M., Aumont, O., Bopp, L., Buitenhuis, E., Doney, S. C., Dunne, J., Hashioka, T., Hauck, J., Hirata, T., John, J., Le Quéré, C., Lima, I. D., Nakano, H., Séférian, R., Totterdell, I., Vichi, M., and Völker, C.: Drivers and uncertainties of future global marine primary production in marine ecosystem models, Biogeosciences, 12, 6955–6984, https://doi.org/10.5194/bg-12-6955-2015, 2015. 

Lauvset, S. K., Gruber, N., Landschützer, P., Olsen, A., and Tjiputra, J.: Trends and drivers in global surface ocean pH over the past 3 decades, Biogeosciences, 12, 1285–1298, https://doi.org/10.5194/bg-12-1285-2015, 2015. 

Levin, L. A. and Bris, N. L.: The deep ocean under climate change, Science, 350, 766–768, https://doi.org/10.1126/science.aad0126, 2015. 

Lotze, H. K., Tittensor, D. P., Bryndum-Buchholz, A., Eddy, T. D., Cheung, W. W. L., Galbraith, E. D., Barange, M., Barrier, N., Bianchi, D., Blanchard, J. L., Bopp, L., Büchner, M., Bulman, C. M., Carozza, D. A., Christensen, V., Coll, M., Dunne, J. P., Fulton, E. A., Jennings, S., Jones, M. C., Mackinson, S., Maury, O., Niiranen, S., Oliveros-Ramos, R., Roy, T., Fernandes, J. A., Schewe, J., Shin, Y.-J., Silva, T. A. M., Steenbeek, J., Stock, C. A., Verley, P., Volkholz, J., Walker, N. D., and Worm, B.: Global ensemble projections reveal trophic amplification of ocean biomass declines with climate change, Proc. Natl. Acad. Sci. USA, 116, 12907–12912, https://doi.org/10.1073/pnas.1900194116, 2019. 

Lovenduski, N. S., McKinley, G. A., Fay, A. R., Lindsay, K., and Long, M. S.: Partitioning uncertainty in ocean carbon uptake projections: Internal variability, emission scenario, and model structure, Global Biogeochem. Cy., 30, 1276–1287, https://doi.org/10.1002/2016GB005426, 2016. 

Lovenduski, N. S. and Bonan, G. B.: Reducing uncertainty in projections of terrestrial carbon uptake, Env. Res. Lett, 12, 044020, https://doi.org/10.1088/1748-9326/aa66b8, 2017. 

Lozier, M. S., Dave, A. C., Palter, J. B., Gerber, L. M., and Barber, R. T.: On the relationship between stratification and primary productivity in the North Atlantic, Geophys. Res. Lett., 38, L18609, https://doi.org/10.1029/2011GL049414, 2011. 

Masson, D. and Knutti, R.: Climate model genealogy, Geophys. Res. Lett., 38, L08703, https://doi.org/10.1029/2011GL046864, 2011. 

Mauritsen, T., Bader, J., Becker, T., Behrens, J., Bittner, M., Brokopf, R., Brovkin, V., Claussen, M., Crueger, T., Esch, M., and Fast, I.: Developments in the MPI-M Earth System Model version 1.2 (MPI-ESM1.2) and its response to increasing CO2, J. Adv. Model. Earth Syst., 11, 998–1038, https://doi.org/10.1029/2018MS001400, 2019. 

McBryan, T. L., Anttila, K., Healy, T. M., and Schulte, P. M.: Responses to Temperature and Hypoxia as Interacting Stressors in Fish: Implications for Adaptation to Environmental Change, Integr. Comp. Biol., 53, 648–659, https://doi.org/10.1093/icb/ict066, 2013. 

McNeil, B. I. and Matear, R. J.: Climate change feedbacks on future oceanic acidification, Tellus B, 59, 191–198, https://doi.org/10.1111/j.1600-0889.2006.00241, 2007. 

McNeil, B. I. and Sasse, T. P.: Future ocean hypercapnia driven by anthropogenic amplification of the natural CO2 cycle, Nature, 529, 383–386, https://doi.org/10.1038/nature16156, 2016. 

Meinshausen, M., Smith, S. J., Calvin, K., Daniel, J. S., Kainuma, M. L. T., Lamarque, J.-F., Matsumoto, K., Montzka, S. A., Raper, S. C. B., Riahi, K., Thomson, A., Velders, G. J. M., and van Vuuren, D. P. P.: The RCP greenhouse gas concentrations and their extensions from 1765 to 2300, Clim. Change, 109, 213, https://doi.org/10.1007/s10584-011-0156-z, 2011. 

Meinshausen, M., Nicholls, Z., Lewis, J., Gidden, M. J., Vogel, E., Freund, M., Beyerle, U., Gessner, C., Nauels, A., Bauer, N., Canadell, J. G., Daniel, J. S., John, A., Krummel, P., Luderer, G., Meinshausen, N., Montzka, S. A., Rayner, P., Reimann, S., Smith, S. J., van den Berg, M., Velders, G. J. M., Vollmer, M., and Wang, H. J.: The SSP greenhouse gas concentrations and their extensions to 2500, Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2019-222, in review, 2019. 

Menary, M. B. and Wood, R. A.: An anatomy of the projected North Atlantic warming hole in CMIP5 models, Clim. Dynam., 50, 3063–3080, https://doi.org/10.1007/s00382-017-3793-8, 2018. 

Moore, J. K., Fu, W., Primeau, F., Britten, G. L., Lindsay, K., Long, M., Doney, S. C., Mahowald, N., Hoffman, F., and Randerson, J. T.: Sustained climate warming drives declining marine biological productivity, Science, 359, 1139–1143, https://doi.org/10.1126/science.aao6379, 2018. 

Munday, P. L., Dixson, D. L., Donelson, J. M., Jones, G. P., Pratchett, M. S., Devitsina, G. V., and Døving, K. B.: Ocean acidification impairs olfactory discrimination and homing ability of a marine fish, Proc. Natl. Acad. Sci. USA, 106, 1848–1852, https://doi.org/10.1073/pnas.0809996106, 2009. 

Müller, W. A., Jungclaus, J. H., Mauritsen, T., Baehr, J., Bittner, M., Budich, R., Bunzel, F., Esch, M., Ghosh, R., Haak, H., and Ilyina, T.: A higher-resolution version of the Max Planck Institute Earth System Model (MPI-ESM1.2-HR), J. Adv. Model. Earth Syst., 10, 1383–1413, https://doi.org/10.1029/2017MS001217, 2018. 

O'Neill, B. C., Tebaldi, C., van Vuuren, D. P., Eyring, V., Friedlingstein, P., Hurtt, G., Knutti, R., Kriegler, E., Lamarque, J.-F., Lowe, J., Meehl, G. A., Moss, R., Riahi, K., and Sanderson, B. M.: The Scenario Model Intercomparison Project (ScenarioMIP) for CMIP6, Geosci. Model Dev., 9, 3461–3482, https://doi.org/10.5194/gmd-9-3461-2016, 2016. 

Orr, J. C.: Recent and future changes in ocean carbonate chemistry, Ocean Acidification, 1, 41–66, 2011. 

Orr, J. C., Fabry, V. J., Aumont, O., Bopp, L., Doney, S. C., Feely, R. A., Gnanadesikan, A., Gruber, N., Ishida, A., Joos, F., Key, R. M., Lindsay, K., Maier-Reimer, E., Matear, R., Monfray, P., Mouchet, A., Najjar, R. G., Plattner, G.-K., Rodgers, K. B., Sabine, C. L., Sarmiento, J. L., Schlitzer, R., Slater, R. D., Totterdell, I. J., Weirig, M.-F., Yamanaka, Y., and Yool, A.: Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms, Nature, 437, 681–686, https://doi.org/10.1038/nature04095, 2005. 

Orr, J. C., Najjar, R. G., Aumont, O., Bopp, L., Bullister, J. L., Danabasoglu, G., Doney, S. C., Dunne, J. P., Dutay, J.-C., Graven, H., Griffies, S. M., John, J. G., Joos, F., Levin, I., Lindsay, K., Matear, R. J., McKinley, G. A., Mouchet, A., Oschlies, A., Romanou, A., Schlitzer, R., Tagliabue, A., Tanhua, T., and Yool, A.: Biogeochemical protocols and diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP), Geosci. Model Dev., 10, 2169–2199, https://doi.org/10.5194/gmd-10-2169-2017, 2017. 

Oschlies, A., Brandt, P., Stramma, L., and Schmidtko, S.: Drivers and mechanisms of ocean deoxygenation, Nat. Geosci., 11, 467–473, https://doi.org/10.1038/s41561-018-0152-2, 2018. 

Pithan, F. and Mauritsen, T.: Arctic amplification dominated by temperature feedbacks in contemporary climate models, Nat. Geosci., 7, 181–184, https://doi.org/10.1038/ngeo2071, 2014. 

Pörtner, H.-O.: Ecosystem effects of ocean acidification in times of ocean warming: a physiologists view, Mar. Ecol. Prog. Ser., 373, 203–217, https://doi.org/10.3354/meps07768, 2008. 

Pörtner, H.-O.: Oxygen- and capacity-limitation of thermal tolerance: a matrix for integrating climate-related stressor effects in marine ecosystems, J. Exp. Biol., 213, 881–893, https://doi.org/10.1242/jeb.037523, 2010. 

Resplandy, L., Bopp, L., Orr, J. C., and Dunne, J. P.: Role of mode and intermediate waters in future ocean acidification: Analysis of CMIP5 models, Geophys. Res. Lett., 40, 3091–3095, https://doi.org/10.1002/grl.50414, 2013. 

Riahi, K., van Vuuren, D. P., Kriegler, E., Edmonds, J., O'Neill, B. C., Fujimori, S., Bauer, N., Calvin, K., Dellink, R., Fricko, O., Lutz, W., Popp, A., Cuaresma, J. C., Kc, S., Leimbach, M., Jiang, L., Kram, T., Rao, S., Emmerling, J., Ebi, K., Hasegawa, T., Havlik, P., Humpenöder, F., Da Silva, L. A., Smith, S., Stehfest, E., Bosetti, V., Eom, J., Gernaat, D., Masui, T., Rogelj, J., Strefler, J., Drouet, L., Krey, V., Luderer, G., Harmsen, M., Takahashi, K., Baumstark, L., Doelman, J. C., Kainuma, M., Klimont, Z., Marangoni, G., Lotze-Campen, H., Obersteiner, M., Tabeau, A., and Tavoni, M.: The Shared Socioeconomic Pathways and their energy, land use, and greenhouse gas emissions implications: An overview, Global Environ. Change, 42, 153–168, https://doi.org/10.1016/j.gloenvcha.2016.05.009, 2017. 

Riche, O. G. J. and Christian, J. R.: Ocean dinitrogen fixation and its potential effects on ocean primary production in Earth system model simulations of anthropogenic warming, Elem. Sci. Anth., 6, 16, https://doi.org/10.1525/elementa.277, 2018. 

Sabine, C. L., Feely, R. A., Gruber, N., Key, R. M., Lee, K., Bullister, J. L., Wanninkhof, R., Wong, C. S., Wallace, D. W. R., Tilbrook, B., Millero, F. J., Peng, T.-H., Kozyr, A., Ono, T., and Rios, A. F.: The Oceanic Sink for Anthropogenic CO2, Science, 305, 367–371, https://doi.org/10.1126/science.1097403, 2004. 

Sanderson, B. M., Knutti, R., and Caldwell, P.: A Representative Democracy to Reduce Interdependency in a Multimodel Ensemble, J. Climate, 28, 5171–5194, https://doi.org/10.1175/JCLI-D-14-00362.1, 2015. 

Sarmiento, J. L., Slater, R., Barber, R., Bopp, L., Doney, S. C., Hirst, A. C., Kleypas, J., Matear, R., Mikolajewicz, U., Monfray, P., Soldatov, V., Spall, S. A., and Stouffer, R.: Response of ocean ecosystems to climate warming, Global Biogeochem. Cy., 18, GB3003, https://doi.org/10.1029/2003GB002134, 2004. 

Schmidtko, S., Stramma, L., and Visbeck, M.: Decline in global oceanic oxygen content during the past five decades, Nature, 542, 335–339, https://doi.org/10.1038/nature21399, 2017. 

Schmittner, A., Oschlies, A., Matthews, H. D., and Galbraith, E. D.: Future changes in climate, ocean circulation, ecosystems, and biogeochemical cycling simulated for a business-as-usual CO2 emission scenario until year 4000 AD, Global Biogeochem. Cy., 22, GB1013, https://doi.org/10.1029/2007GB002953, 2008. 

Schupfner, M., Wieners, K.-H., Wachsmann, F., Steger, C., Bittner, M., Jungclaus, J., Früh, B., Pankatz, K., Giorgetta, M., Reick, C., Legutke, S., Esch, M., Gayler, V., Haak, H., de Vrese, P., Raddatz, T., Mauritsen, T., von Storch, J.-S., Behrens, J., Brovkin, V., Claussen, M., Crueger, T., Fast, I., Fiedler, S., Hagemann, S., Hohenegger, C., Jahns, T., Kloster, S., Kinne, S., Lasslop, G., Kornblueh, L., Marotzke, J., Matei, D., Meraner, K., Mikolajewicz, U., Modali, K., Müller, W., Nabel, J., Notz, D., Peters, K., Pincus, R., Pohlmann, H., Pongratz, J., Rast, S., Schmidt, H., Schnur, R., Schulzweida, U., Six, K., Stevens, B., Voigt, A., and Roeckner, E.: DKRZ MPI-ESM1.2-HR model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.4398, 2019. 

Screen, J. A. and Simmonds, I.: The central role of diminishing sea ice in recent Arctic temperature amplification, Nature, 464, 1334–1337, https://doi.org/10.1038/nature09051, 2010. 

Séférian, R., Gehlen, M., Bopp, L., Resplandy, L., Orr, J. C., Marti, O., Dunne, J. P., Christian, J. R., Doney, S. C., Ilyina, T., Lindsay, K., Halloran, P. R., Heinze, C., Segschneider, J., Tjiputra, J., Aumont, O., and Romanou, A.: Inconsistent strategies to spin up models in CMIP5: implications for ocean biogeochemical model performance assessment, Geosci. Model Dev., 9, 1827–1851, https://doi.org/10.5194/gmd-9-1827-2016, 2016. 

Séférian, R.: CNRM-CERFACS CNRM-ESM2-1 model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.4068, 2018. 

Séférian, R.: CNRM-CERFACS CNRM-ESM2-1 model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.4191, 2019. 

Séférian, R., Nabat, P., Michou, M., Saint-Martin, D., Voldoire, A., Colin, J., Decharme, B., Delire, C., Berthet, S., Chevallier, M., Sénési, S., Franchistéguy, L., Vial, J., Mallet, M., Joetzjer, E., Geoffroy, O., Guérémy, J.-F., Moine, M.-P., M'Sadek, R., Ribes, A., Rocher, M., Roehrig, R., Salas-y-Mélia, D., Sanchez, E., Terray, L., Valcke, S., Waldman, R., Aumont, O., Bopp, L., Deshayes, J., Éthé, C., and Madec, G.: Evaluation of CNRM Earth System model, CNRM-ESM2-1: role of Earth system processes in present-day and future climate, J. Adv. Model. Earth Syst., 11, 4182–4227, https://doi.org/10.1029/2019MS001791, 2019. 

Séférian, R., Berthet, S., Yool, A., Palmiéri, J., Bopp, L., Tagliabue, A., Kwiatkowski, L., Aumont, O., Christian, J., Dunne, J., Gehlen, M., Ilyina, T., John, J. G., Li, H., Long, M., Luo, J. Y., Nakano, H., Romanou, A., Schwinger, J., Stock, C., Santana-Falcón, Y., Takano, Y., Tjiputra, J., Tsujino, H., Watanabe, M., Wu, T., Wu, F., Yamamoto, A. : Tracking improvement in simulated marine biogeochemistry between CMIP5 and CMIP6, Current Climate Change Reports, doi:10.1007/s40641-020-00160-0, 2020. 

Seland, Ø., Bentsen, M., Oliviè, D. J. L., Toniazzo, T., Gjermundsen, A., Graff, L. S., Debernard, J. B., Gupta, A. K., He, Y., Kirkevåg, A., Schwinger, J., Tjiputra, J., Aas, Kjetil S., Bethke, I., Fan, Y., Griesfeller, J., Grini, A., Guo, C., Ilicak, M., Karset, I. H. H., Landgren, O. A., Liakka, J., Moseid, K. O., Nummelin, A., Spensberger, C., Tang, H., Zhang, Z., Heinze, C., Iversen, T., and Schulz, M.: NCC NorESM2-LM model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.8036, 2019a. 

Seland, Ø., Bentsen, M., Oliviè, D. J. L., Toniazzo, T., Gjermundsen, A., Graff, L. S., Debernard, J. B., Gupta, A. K., He, Y., Kirkevåg, A., Schwinger, J., Tjiputra, J., Aas, Kjetil S., Bethke, I., Fan, Y., Griesfeller, J., Grini, A., Guo, C., Ilicak, M., Karset, I. H. H., Landgren, O. A., Liakka, J., Moseid, K. O., Nummelin, A., Spensberger, C., Tang, H., Zhang, Z., Heinze, C., Iversen, T., and Schulz, M.: NCC NorESM2-LM model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.8319, 2019b. 

Sellar, A. A., Jones, C. G., Mulcahy, J., Tang, Y., Yool, A., Wiltshire, A., O'connor, F. M., Stringer, M., Hill, R., Palmieri, J., Woodward, S., Mora, L., Kuhlbrodt, T., Rumbold, S., Kelley, D. I., Ellis, R., Johnson, C. E., Walton, J., Abraham, N. L., Andrews, M. B., Andrews, T., Archibald, A. T., Berthou, S., Burke, E., Blockley, E., Carslaw, K., Dalvi, M., Edwards, J., Folberth, G. A., Gedney, N., Griffiths, P. T., Harper, A. B., Hendry, M. A., Hewitt, A. J., Johnson, B., Jones, A., Jones, C. D., Keeble, J., Liddicoat, S., Morgenstern, O., Parker, R. J., Predoi, V., Robertson, E., Siahaan, A., Smith, R. S., Swaminathan, R., Woodhouse, M. T., Zeng, G., and Zerroukat, M.: UKESM1: Description and evaluation of the UK Earth System Model, J. Adv. Model. Earth Syst., 11, 4513–4558, https://doi.org/10.1029/2019MS001739, 2019. 

Serreze, M. C., Barrett, A. P., Slater, A. G., Steele, M., Zhang, J., and Trenberth, K. E.: The large-scale energy budget of the Arctic, J. Geophys. Res.-Atmos., 112, D11122, https://doi.org/10.1029/2006JD008230, 2007. 

Smale, D. A., Wernberg, T., Oliver, E. C. J., Thomsen, M., Harvey, B. P., Straub, S. C., Burrows, M. T., Alexander, L. V., Benthuysen, J. A., Donat, M. G., Feng, M., Hobday, A. J., Holbrook, N. J., Perkins-Kirkpatrick, S. E., Scannell, H. A., Gupta, A. S., Payne, B. L., and Moore, P. J.: Marine heatwaves threaten global biodiversity and the provision of ecosystem services, Nat. Clim. Change, 1, 306–312, https://doi.org/10.1038/s41558-019-0412-1, 2019. 

Sorte, C. J. B., Williams, S. L., and Carlton, J. T.: Marine range shifts and species introductions: comparative spread rates and community impacts, Global Ecol. Biogeogr., 19, 303–316, https://doi.org/10.1111/j.1466-8238.2009.00519.x, 2010. 

Steinacher, M., Joos, F., Frölicher, T. L., Plattner, G.-K., and Doney, S. C.: Imminent ocean acidification in the Arctic projected with the NCAR global coupled carbon cycle-climate model, Biogeosciences, 6, 515–533, https://doi.org/10.5194/bg-6-515-2009, 2009. 

Steinacher, M., Joos, F., Frölicher, T. L., Bopp, L., Cadule, P., Cocco, V., Doney, S. C., Gehlen, M., Lindsay, K., Moore, J. K., Schneider, B., and Segschneider, J.: Projected 21st century decrease in marine productivity: a multi-model analysis, Biogeosciences, 7, 979–1005, https://doi.org/10.5194/bg-7-979-2010, 2010. 

Steiner, N. S., Christian, J. R., Six, K. D., Yamamoto, A., and Yamamoto-Kawai, M.: Future ocean acidification in the Canada Basin and surrounding Arctic Ocean from CMIP5 earth system models, J. Geophys. Res.-Oceans, 119, 332–347, https://doi.org/10.1002/2013JC009069, 2014. 

Stock, C. A., Dunne, J. P., and John, J. G.: Drivers of trophic amplification of ocean productivity trends in a changing climate, Biogeosciences, 11, 7125–7135, https://doi.org/10.5194/bg-11-7125-2014, 2014. 

Stock, C. A., Dunne, J. P., Fan, S., Ginoux, P., John, J. G., Krasting, J. P., Laufkötter, C., Paulot, F., and Zadeh, N.,: Ocean Biogeochemistry in GFDL's Earth System Model 4.1 and its Response to Increasing atmospheric CO2, J. Adv. Model. Earth Syst., 2019MS002043, 2020. 

Stramma, L., Johnson, G. C., Sprintall, J., and Mohrholz, V.: Expanding Oxygen-Minimum Zones in the Tropical Oceans, Science, 320, 655–658, https://doi.org/10.1126/science.1153847, 2008. 

Stramma, L., Oschlies, A., and Schmidtko, S.: Mismatch between observed and modeled trends in dissolved upper-ocean oxygen over the last 50 yr, Biogeosciences, 9, 4045–4057, https://doi.org/10.5194/bg-9-4045-2012, 2012. 

Sunagawa, S., Coelho, L. P., Chaffron, S., Kultima, J. R., Labadie, K., Salazar, G., Djahanschiri, B., Zeller, G., Mende, D. R., Alberti, A., Cornejo-Castillo, F. M., Costea, P. I., Cruaud, C., d'Ovidio, F., Engelen, S., Ferrera, I., Gasol, J. M., Guidi, L., Hildebrand, F., Kokoszka, F., Lepoivre, C., Lima-Mendez, G., Poulain, J., Poulos, B. T., Royo-Llonch, M., Sarmento, H., Vieira-Silva, S., Dimier, C., Picheral, M., Searson, S., Kandels-Lewis, S., Coordinators, T. O., Bowler, C., Vargas, C. de, Gorsky, G., Grimsley, N., Hingamp, P., Iudicone, D., Jaillon, O., Not, F., Ogata, H., Pesant, S., Speich, S., Stemmann, L., Sullivan, M. B., Weissenbach, J., Wincker, P., Karsenti, E., Raes, J., Acinas, S. G., and Bork, P.: Structure and function of the global ocean microbiome, Science, 348, 6237, https://doi.org/10.1126/science.1261359, 2015. 

Swart, N. C., Cole, J. N. S., Kharin, V. V., Lazare, M., Scinocca, J. F., Gillett, N. P., Anstey, J., Arora, V., Christian, J. R., Hanna, S., Jiao, Y., Lee, W. G., Majaess, F., Saenko, O. A., Seiler, C., Seinen, C., Shao, A., Sigmond, M., Solheim, L., von Salzen, K., Yang, D., and Winter, B.: The Canadian Earth System Model version 5 (CanESM5.0.3), Geosci. Model Dev., 12, 4823–4873, https://doi.org/10.5194/gmd-12-4823-2019, 2019a. 

Swart, N. C., Cole, J. N. S., Kharin, V. V., Lazare, M., Scinocca, J. F., Gillett, N. P., Anstey, J., Arora, V., Christian, J. R., Jiao, Y., Lee, W. G., Majaess, F., Saenko, O. A., Seiler, C., Seinen, C., Shao, A., Solheim, L., von Salzen, K., Yang, D., Winter, B., and Sigmond, M.: CCCma CanESM5 model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.3673, 2019b. 

Swart, N. C., Cole, J. N. S., Kharin, V. V., Lazare, M., Scinocca, J. F., Gillett, N. P., Anstey, J., Arora, V., Christian, J. R., Jiao, Y., Lee, W. G., Majaess, F., Saenko, O. A., Seiler, C., Seinen, C., Shao, A., Solheim, L., von Salzen, K., Yang, D., Winter, B., and Sigmond, M.: CCCma CanESM5 model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.3683, 2019c. 

Sweetman, A. K., Thurber, A. R., Smith, C. R., Levin, L. A., Mora, C., Wei, C.-L., Gooday, A. J., Jones, D. O. B., Rex, M., Yasuhara, M., Ingels, J., Ruhl, H. A., Frieder, C. A., Danovaro, R., Würzberg, L., Baco, A., Grupe, B. M., Pasulka, A., Meyer, K. S., Dunlop, K. M., Henry, L.-A., and Roberts, J. M.: Major impacts of climate change on deep-sea benthic ecosystems, Elementa, 5, 4, https://doi.org/10.1525/elementa.203, 2017. 

Tachiiri, K., Abe, M., Hajima, T., Arakawa, O., Suzuki, T., Komuro, Y., Ogochi, K., Watanabe, M., Yamamoto, A., Tatebe, H., Noguchi, M. A., Ohgaito, R., Ito, A., Yamazaki, D., Ito, A., Takata, K., Watanabe, S., and Kawamiya, M.: MIROC MIROC-ES2L model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.5745, 2019. 

Tagliabue, A., Aumont, O., DeAth, R., Dunne, J. P., Dutkiewicz, S., Galbraith, E., Misumi, K., Moore, J. K., Ridgwell, A., Sherman, E., Stock, C., Vichi, M., Völker, C. and Yool, A.: How well do global ocean biogeochemistry models simulate dissolved iron distributions?, Global Biogeochem. Cy., 30, 149–174, https://doi.org/10.1002/2015GB005289, 2016. 

Tang, Y., Rumbold, S., Ellis, R., Kelley, D., Mulcahy, J., Sellar, A., Walton, J. and Jones, C.: MOHC UKESM1.0-LL model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.6113, 2019. 

Taucher, J. and Oschlies, A.: Can we predict the direction of marine primary production change under global warming?, Geophys. Res. Lett., 38, L02603, https://doi.org/10.1029/2010GL045934, 2011. 

Taylor, K. E., Stouffer, R. J., and Meehl, G. A.: An Overview of CMIP5 and the Experiment Design, B. Am. Meteorol. Soc., 93, 485–498, https://doi.org/10.1175/BAMS-D-11-00094.1, 2011. 

Thomas, M. K., Kremer, C. T., Klausmeier, C. A., and Litchman, E.: A Global Pattern of Thermal Adaptation in Marine Phytoplankton, Science, 338, 1085–1088, https://doi.org/10.1126/science.1224836, 2012. 

Tittensor, D. P., Baco, A. R., Hall-Spencer, J. M., Orr, J. C., and Rogers, A. D.: Seamounts as refugia from ocean acidification for cold-water stony corals, Mar. Ecol., 31, 212–225, https://doi.org/10.1111/j.1439-0485.2010.00393.x, 2010. 

Tjiputra, J. F., Goris, N., Lauvset, S. K., Heinze, C., Olsen, A.,Schwinger, J., and Steinfeldt, R.: Mechanisms and Early Detections of Multidecadal Oxygen Changes in the Interior Subpolar North Atlantic, Geophys. Res. Lett., 45, 4218–4229, https://doi.org/10.1029/2018GL077096, 2018. 

Tjiputra, J. F., Schwinger, J., Bentsen, M., Morée, A. L., Gao, S., Bethke, I., Heinze, C., Goris, N., Gupta, A., He, Y.-C., Olivié, D., Seland, Ø., and Schulz, M.: Ocean biogeochemistry in the Norwegian Earth System Model version 2 (NorESM2), Geosci. Model Dev., 13, 2393–2431, https://doi.org/10.5194/gmd-13-2393-2020, 2020. 

Vaquer-Sunyer, R. and Duarte, C. M.: Thresholds of hypoxia for marine biodiversity, P. Natl. Acad. Sci. USA, 105, 15452–15457, https://doi.org/10.1073/pnas.0803833105, 2008. 

Vial, J., Dufresne, J.-L., and Bony, S.: On the interpretation of inter-model spread in CMIP5 climate sensitivity estimates, Clim. Dynam., 41, 3339–3362, https://doi.org/10.1007/s00382-013-1725-9, 2013. 

Vichi, M., Manzini, E., Fogli, P. G., Alessandri, A., Patara, L., Scoccimarro, E., Masina, S., and Navarra, A.: Global and regional ocean carbon uptake and climate change: sensitivity to a substantial mitigation scenario, Clim. Dynam., 37, 1929–1947, https://doi.org/10.1007/s00382-011-1079-0, 2011. 

Watanabe, M. and Kawamiya, M.: Remote effects of mixed layer development on ocean acidification in the subsurface layers of the North Pacific, J. Oceanogr., 73, 771–784, https://doi.org/10.1007/s10872-017-0431-3, 2017. 

Watson, S.-A., Lefevre, S., McCormick, M. I., Domenici, P., Nilsson, G. E. and Munday, P. L.: Marine mollusc predator-escape behaviour altered by near-future carbon dioxide levels, Proc. Roy. Soc. Lond. B-Biol. Sci., 281, 20132377, https://doi.org/10.1098/rspb.2013.2377, 2014. 

Watson, S.-A., Fields, J. B., and Munday, P. L.: Ocean acidification alters predator behaviour and reduces predation rate, Biol. Lett., 13, 20160797, https://doi.org/10.1098/rsbl.2016.0797, 2017. 

Whitt, D. B.: On the Role of the Gulf Stream in the Changing Atlantic Nutrient Circulation During the 21st Century, in Kuroshio Current: Physical, Biogeochemical, and Ecosystem Dynamics, 51–82, American Geophysical Union (AGU), 2019. 

Yamamoto, A., Kawamiya, M., Ishida, A., Yamanaka, Y., and Watanabe, S.: Impact of rapid sea-ice reduction in the Arctic Ocean on the rate of ocean acidification, Biogeosciences, 9, 2365–2375, https://doi.org/10.5194/bg-9-2365-2012, 2012. 

Yamamoto-Kawai, M., McLaughlin, F. A., Carmack, E. C., Nishino, S., and Shimada, K.: Aragonite Undersaturation in the Arctic Ocean: Effects of Ocean Acidification and Sea Ice Melt, Science, 326, 1098–1100, https://doi.org/10.1126/science.1174190, 2009. 

Yukimoto, S., Kawai, H., Koshiro, T., Oshima, N., Yoshida, K., Urakawa, S., Tsujino, H., Deushi, M., Tanaka, T., Hosaka, M., Yabu, S., Yoshimura, H., Shindo, E., Mizuta, R., Obata, A., Adachi, Y., and Ishii, M.: The Meteorological Research Institute Earth System Model version 2.0, MRI-ESM2.0: Description and basic evaluation of the physical component, J. Meteorol. Soc. Jpn., 97, 931–965, https://doi.org/10.2151/jmsj.2019-051, 2019a. 

Yukimoto, S., Koshiro, T., Kawai, H., Oshima, N., Yoshida, K., Urakawa, S., Tsujino, H., Deushi, M., Tanaka, T., Hosaka, M., Yoshimura, H., Shindo, E., Mizuta, R., Ishii, M., Obata, A., and Adachi, Y.: MRI MRI-ESM2.0 model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.6842, 2019b. 

Yukimoto, S., Koshiro, T., Kawai, H., Oshima, N., Yoshida, K., Urakawa, S., Tsujino, H., Deushi, M., Tanaka, T., Hosaka, M., Yoshimura, H., Shindo, E., Mizuta, R., Ishii, M., Obata, A., and Adachi, Y.: MRI MRI-ESM2.0 model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.6929, 2019c. 

Ziehn, T., Chamberlain, M., Law, R., Lenton, A., Bodman, R., Dix, M., Stevens, L., Wang, Y.-P., and Srbinovsky, J.: The Australian Earth System Model:ACCESS-ESM1.5, Journal of Southern Hemisphere Earth Systems Science (JSHESS), 2020. 

Ziehn, T., Chamberlain, M., Lenton, A., Law, R., Bodman, R., Dix, M., Wang, Y., Dobrohotoff, P., Srbinovsky, J., Stevens, L., Vohralik, P., Mackallah, C., Sullivan, A., O'Farrell, S., and Druken, K.: CSIRO ACCESS-ESM1.5 model output prepared for CMIP6 CMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.4231, 2019a.  

Ziehn, T., Chamberlain, M., Lenton, A., Law, R., Bodman, R., Dix, M., Wang, Y., Dobrohotoff, P., Srbinovsky, J., Stevens, L., Vohralik, P., Mackallah, C., Sullivan, A., O'Farrell, S., and Druken, K.: CSIRO ACCESS-ESM1.5 model output prepared for CMIP6 ScenarioMIP, Earth System Grid Federation, available at: https://doi.org/10.22033/ESGF/CMIP6.4320, 2019b. 

Download
Short summary
We assess 21st century projections of marine biogeochemistry in the CMIP6 Earth system models. These models represent the most up-to-date understanding of climate change. The models generally project greater surface ocean warming, acidification, subsurface deoxygenation, and euphotic nitrate reductions but lesser primary production declines than the previous generation of models. This has major implications for the impact of anthropogenic climate change on marine ecosystems.
Altmetrics
Final-revised paper
Preprint