The Arctic experiences climate changes that are among the fastest in the world and affect all Earth system components. Despite expected increase in terrigenous inputs to the Arctic Ocean, their impacts on biogeochemical cycles are currently largely neglected in IPCC-like models. We used a state-of-the-art high-resolution ocean biogeochemistry model, that includes carbon and nutrient inputs from rivers and coastal erosion, to produce twenty-first-century pan-Arctic projections. Surprisingly, even with an anticipated rise in primary production across a wide range of emission scenarios, our findings indicate that climate change leads to a counterintuitive 40% reduction in the efficiency of the Arctic's biological carbon pump by 2100, to which terrigenous inputs contribute 10%. Terrigenous inputs will also drive intense coastal CO 2 outgassing, reducing the Arctic Ocean's carbon sink by at least 10% (33 TgC yr -1 ). These unexpected positive feedbacks, mostly due to accelerated remineralization rates, lower the Arctic Ocean’s capacity for sequestering carbon.

The air-sea CO2 flux FCO2 is an important component of the global carbon budget and understanding its response to climate change is crucial to adjust mitigation pathways. Multi-linear regression supports the expectation that the balance between the CO2 partial pressures of air and the sea surface (pCO2) is the most important driver of temporal FCO2 variability. Discrepancies between state-of-the-art Earth System Models (ESMs) and gridded pCO2-products suggest that systematic biases exist across an ensemble of ESMs. In the equatorial regions, upwelling variability of carbon-rich water is biased in ESMs as modeled and observed sea surface temperature are generally uncorrelated. In the high latitudes, the climate change induced trend towards lighter sea water is overestimated in ESMs, which yields - in contrast to observations - shallower mixed layers over the contemporary period and hence a suppressed carbon supply from depth. While mixed layer depth variability and trends appear biased throughout the global ocean, this is not a determining factor for pCO2 variability in subtropical gyres. The results highlight the importance of accurately modeling hydrographic properties to obtain robust estimates of FCO2 and its variability.

The Southern Ocean is a major region of ocean carbon uptake, but its future changes remain uncertain under climate change. Here we show the projected shift in the Southern Ocean CO2 sink using a suite of Earth System Models, revealing changes in the mechanism, position and seasonality of the carbon uptake. The region of dominant CO2 uptake shifts from the Subtropical to the Antarctic region under the high-emission scenario. The warming-driven sea-ice melt, increased ocean stratification, mixed layer shoaling, and a weaker vertical carbon gradient is projected to together reduce the winter de-gassing in the future, which will trigger the switch from mixing-driven outgassing to solubility-driven uptake in the Antarctic region during the winter season. The future Southern Ocean carbon sink will be poleward-shifted, operating in a hybrid mode between biologically-driven summertime and solubility-driven wintertime uptake with further amplification of biologically-driven uptake due to the increasing Revelle Factor.

As the climate warms, the grounded ice sheet and floating ice shelves surrounding Antarctica are melting and releasing additional freshwater into the Southern Ocean. Nonetheless, almost all existing coupled climate models have fixed ice sheets and lack the physics required to represent the dominant sources of Antarctic melt. These missing ice dynamics represent a key uncertainty that is typically unaccounted for in current global climate change projections. Previous modelling studies that have imposed additional Antarctic meltwater have demonstrated regional impacts on Southern Ocean stratification, circulation, and sea ice, as well as remote changes in atmospheric circulation, tropical precipitation, and global temperature. However, these previous studies have used widely varying rates of freshwater forcing, have been conducted using different climate models and configurations, and have reached differing conclusions on the magnitude of meltwater–climate feedbacks. The Southern Ocean Freshwater Input from Antarctica (SOFIA) initiative brings together a team of scientists to quantify the climate system response to Antarctic meltwater input along with key aspects of the uncertainty. In this paper, we summarize the state of knowledge on meltwater discharge from the Antarctic ice sheet and ice shelves to the Southern Ocean and explain the scientific objectives of our initiative. We propose a series of coupled and ocean–sea ice model experiments, including idealized meltwater experiments, historical experiments with observationally consistent meltwater input, and future scenarios driven by meltwater inputs derived from stand-alone ice sheet models. Through coordinating a multi-model ensemble of simulations using a common experimental design, open data archiving, and facilitating scientific collaboration, SOFIA aims to move the community toward better constraining our understanding of the climate system response to Antarctic melt.

Small-scale eddies play an important role in preconditioning and restratifying the water column before and after mixing events, thereby affecting deep water formation variability. Results from a realistic eddy-resolving (~5 km local horizontal resolution) ocean model suggest that small-scale temperature fluxes due to turbulent potential to kinetic energy conversion are the main driver of mixed layer restratification during deep convection in the Labrador Sea interior and the West Greenland Current. This resupply of heat due to turbulent upward buoyancy fluxes exhibits a large interannual variability imposed by the atmospheric forcing. Eddy fluxes only become active in periods of strong buoyancy loss, while being quiescent otherwise. In a low-resolution (~20 km) control simulation the modeled turbulence is strongly reduced and the associated modeled and parameterized heat fluxes are too weak to increase stratification.

Abstract. Abrupt changes in ocean biogeochemical variables occur as a result of human-induced climate forcing as well as those which are more gradual and occur over longer timescales. These abrupt changes have not yet been identified and quantified to the same extent as the more gradual ones. We review and synthesise abrupt changes in ocean biogeochemistry under human-induced climatic forcing. We specifically address the ocean carbon and oxygen cycles because the related processes of acidification and deoxygenation provide important ecosystem hazards. Since biogeochemical cycles depend also on the physical environment, we also describe the relevant changes in warming, circulation, and sea ice. We include an overview of the reversibility or irreversibility of abrupt marine biogeochemical changes. Important implications of abrupt biogeochemical changes for ecosystems are also discussed. We conclude that there is evidence for increasing occurrence and extent of abrupt changes in ocean biogeochemistry as a consequence of rising greenhouse gas emissions.

The Southern Ocean is a major region of ocean carbon uptake, but its future changes remain uncertain under climate warming. Here we show the projected shift in the Southern Ocean CO2 sink using a suite of Earth System Models, revealing changes in the mechanism, position and seasonality of the carbon uptake. Dominant CO2 uptake shifts from the Subtropical to the Antarctic region under the high-emission scenario by the end of the 21st century. The warming-driven sea-ice melt, increased ocean stratification, mixed layer shoaling, and a weaker vertical carbon gradient will together reduce the winter outgassing in the future, which will trigger the switch from mixing-driven outgassing to solubility-driven uptake in the Antarctic region during the winter season. The future Southern Ocean carbon sink will be poleward-shifted, operating in a hybrid mode between biologically-driven summertime and solubility-driven wintertime uptake with further amplification of biological uptake by the increasing Revelle Factor.

Transient simulations of the global fully coupled climate model COSMOS under realistic varying orbital and greenhouse gas forcings are systematically compared to diatom oxygen isotope ( δ 18 O diatom ) records from Russian lakes with focus on Eurasian Holocene climate trends. The measured δ 18 O diatom decrease and other temperature proxies are interpreted as large‐scale cooling throughout the Holocene while the model simulations are biased too warm, likely through missing radiative forcings. This large‐scale warm bias also dictates the modeled δ 18 O precipitation . Hence, at locations where the signs of model and proxy temperature/precipitation trends agree, measured δ 18 O diatom and modeled δ 18 O precipitation trends show notable accordance. An increased temporal variability of modeled δ 18 O precipitation is linked to persistent atmospheric circulation patterns. Applying the transient forcings in an accelerated way (every 10th year only) yields a similar, yet weaker or delayed model response, especially in the ocean.

The climate sensitivity is known to be mainly determined by the atmosphere model but here we discover that the ocean model can change a given transient climate response (TCR) by as much as 20% while the equilibrium climate sensitivity (ECS) change is limited to 10%. In our study, two different coupled CMIP6 models (MPI-ESM and AWI-CM) in two different resolutions each are compared. The coupled models share the same atmosphere-land component ECHAM6.3, which has been developed at the Max-Planck-Institute for Meteorology (MPI-M). However, as part of MPI-ESM and AWI-CM, ECHAM6.3 is coupled to two different ocean models, namely the MPIOM sea ice-ocean model developed at MPI-M and the FESOM sea ice-ocean model developed at the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI). A reason for the different TCR is different ocean heat uptake through greenhouse gas forcing in AWI simulations compared to MPI-M simulations. Specifically, AWI-CM simulations show stronger surface heating than MPI-ESM simulations while the MPI-M model accumulates more heat in the deeper ocean. The vertically integrated ocean heat content is increasing stronger in MPI-M model configurations compared to AWI model configurations in the high latitudes. Strong vertical mixing in MPI-M model configurations compared to AWI model configurations seems to be key for these differences. The strongest difference in vertical ocean mixing occurs inside the Weddell Gyre, but there are also important differences in another key region, the northern North Atlantic. Over the North Atlantic, these differences materialize in a lack of a warming hole in AWI model configurations and the presence of a warming hole in MPI-M model configurations. All these differences occur largely independent of the considered model resolutions.

The transient climate response (TCR) is 20% higher in the Alfred Wegener Institute Climate Model (AWI-CM) compared to the Max Planck Institute Earth System Model (MPI-ESM) whereas the equilibrium climate sensitivity (ECS) is by up to 10% higher in AWI-CM. These results are largely independent of the two considered model resolutions for each model. The two coupled CMIP6 models share the same atmosphere-land component ECHAM6.3 developed at the Max Planck Institute for Meteorology (MPI-M). However, ECHAM6.3 is coupled to two different ocean models, namely the MPIOM sea ice-ocean model developed at MPI-M and the FESOM sea ice-ocean model developed at the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI). A reason for the different TCR is related to ocean heat uptake in response to greenhouse gas forcing. Specifically, AWI-CM simulations show stronger surface heating than MPI-ESM simulations while the latter accumulate more heat in the deeper ocean. The vertically integrated ocean heat content is increasing slower in AWI-CM model configurations compared to MPI-ESM model configurations in the high latitudes. Weaker vertical mixing in AWI-CM model configurations compared to MPI-ESM model configurations seems to be key for these differences. The strongest difference in vertical ocean mixing occurs inside the Weddell and Ross Gyres and the northern North Atlantic. Over the North Atlantic, these differences materialize in a lack of a warming hole in AWI-CM model configurations and the presence of a warming hole in MPI-ESM model configurations. All these differences occur largely independent of the considered model resolutions.

Stratospheric ozone and water vapour are key components of the Earth system, and past and future changes to both have important impacts on global and regional climate. Here, we evaluate long-term changes in these species from the pre-industrial period (1850) to the end of the 21st century in Coupled Model Intercomparison Project phase 6 (CMIP6) models under a range of future emissions scenarios. There is good agreement between the CMIP multi-model mean and observations for total column ozone (TCO), although there is substantial variation between the individual CMIP6 models. For the CMIP6 multi-model mean, global mean TCO has increased from ∼ 300 DU in 1850 to ∼ 305 DU in 1960, before rapidly declining in the 1970s and 1980s following the use and emission of halogenated ozone-depleting substances (ODSs). TCO is projected to return to 1960s values by the middle of the 21st century under the SSP2-4.5, SSP3-7.0, SSP4-3.4, SSP4-6.0, and SSP5-8.5 scenarios, and under the SSP3-7.0 and SSP5-8.5 scenarios TCO values are projected to be ∼ 10 DU higher than the 1960s values by 2100. However, under the SSP1-1.9 and SSP1-1.6 scenarios, TCO is not projected to return to the 1960s values despite reductions in halogenated ODSs due to decreases in tropospheric ozone mixing ratios. This global pattern is similar to regional patterns, except in the tropics where TCO under most scenarios is not projected to return to 1960s values, either through reductions in tropospheric ozone under SSP1-1.9 and SSP1-2.6, or through reductions in lower stratospheric ozone resulting from an acceleration of the Brewer–Dobson circulation under other Shared Socioeconomic Pathways (SSPs). In contrast to TCO, there is poorer agreement between the CMIP6 multi-model mean and observed lower stratospheric water vapour mixing ratios, with the CMIP6 multi-model mean underestimating observed water vapour mixing ratios by ∼ 0.5 ppmv at 70 hPa. CMIP6 multi-model mean stratospheric water vapour mixing ratios in the tropical lower stratosphere have increased by ∼ 0.5 ppmv from the pre-industrial to the present-day period and are projected to increase further by the end of the 21st century. The largest increases (∼ 2 ppmv) are simulated under the future scenarios with the highest assumed forcing pathway (e.g. SSP5-8.5). Tropical lower stratospheric water vapour, and to a lesser extent TCO, shows large variations following explosive volcanic eruptions.

We simulate the two Coupled Model Intercomparison Project scenarios RCP4.5 and RCP8.5, to assess the effects of melt‐induced fresh water on the Atlantic meridional overturning circulation (AMOC). We use a newly developed climate model with high resolution at the coasts, resolving the complex ocean dynamics. Our results show an AMOC recovery in simulations run with and without an included ice sheet model. We find that the ice sheet adds a strong decadal variability on the freshwater release, resulting in intervals in which it reduces the surface runoff by high accumulation rates. This compensating effect is missing in climate models without dynamic ice sheets. Therefore, we argue to assess those freshwater hosing experiments critically, which aim to parameterize Greenland's freshwater release. We assume the increasing net evaporation over the Atlantic and the resulting increase in ocean salinity, to be the main driver of the AMOC recovery.

The global-mean surface air temperature change due to a doubled carbon dioxide concentration in the atmosphere (equilibrium climate sensitivity, ECS) is an important measure to quantify the impact of predicted anthropogenic climate change. The latest climate modeling intercomparison project (CMIP6) exhibits a higher ECS compared to the previous climate model generation (1.8 to 5.6 K for CMIP6 versus 1.5 to 4.5 K for CMIP5). The increase in ECS is likely due to decreases in extratropical low cloud coverage and albedo, caused by improvements in the numerical aerosol schemes. Our state-of-the-art Earth system model AWI-ESM, developed at the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, yields an ECS of 3.59-3.62 K, which is close to the CMIP5 mean. Using a set of varying model configurations, we identify dynamic vegetation and model resolution as the primary driving factors which influence the modeled global response to an increased greenhouse gas forcing.

The influence of a high horizontal resolution (5–15 km) on the general circulation and hydrography in the North Atlantic is investigated using the Finite Element Sea Ice–Ocean Model (FESOM). We find a stronger shift of the upper-ocean circulation and water mass properties during the model spinup in the high-resolution model version compared to the low-resolution (~1°) control run. In quasi equilibrium, the high-resolution model is able to reduce typical low-resolution model biases. Especially, it exhibits a weaker salinification of the North Atlantic subpolar gyre and a reduced mixed layer depth in the Labrador Sea. However, during the spinup adjustment, we see that initially improved high-resolution features partially reduce over time: the strength of the Atlantic overturning and the path of the North Atlantic Current are not maintained, and hence hydrographic biases known from low-resolution ocean models return in the high-resolution quasi-equilibrium state. We identify long baroclinic Rossby waves as a potential cause for the strong upper-ocean adjustment of the high-resolution model and conclude that a high horizontal resolution improves the state of the modeled ocean but the model integration length should be chosen carefully.