Global warming alters ocean conditions, which can have dramatic consequences for marine species. Yet, the centennial-scale effects and reversibility of habitat viability for marine species, particularly those that are important to fisheries, remain uncertain. Using the Aerobic Growth Index, we quantify the impacts of warming and deoxygenation on the contemporary habitat volume of 46 exploited marine species in novel temperature stabilization and overshoot simulations until 2500. We demonstrate that only around half of the simulated loss of contemporary (1995–2014) habitat volume is realized when warming levels are first reached. Moreover, in an overshoot scenario peaking at 2 °C global warming before stabilizing at 1.5 °C, the maximum decrease in contemporary habitat volume occurs more than 150 years post-peak warming. Species' adaptation may strongly mitigate impacts depending on adaptation rate and pressure. According to our study, marine species will be affected for centuries after temperature stabilization and overshoot, with impacts surpassing those during the transient warming phase.

Antarctic toothfish are a commercially exploited upper-level predator in the Southern Ocean. As many of its prey, the ectothermic, water-breathing Antarctic toothfish is specifically adapted to the temperature and oxygen conditions present in the high-latitude Southern Ocean. Additionally, the life cycle of Antarctic toothfish depends on sea-ice dynamics and the transport of individuals by currents between regions with different prey. To assess the impact of 21st-century climate change on potential interactions of Antarctic toothfish and its prey, we here employ the extended aerobic growth index (AGI), which quantifies the effect of ocean temperature and oxygen levels on the habitat viability of individual species. We quantify changes in predator-prey interactions by a change in viable habitat overlap as obtained with the AGI. As environmental data, we use future projections for four emission scenarios from the model FESOM-REcoM, which is specifically designed for applications on and near the Antarctic continental shelf. For the two highest-emission scenarios, we find that warming and deoxygenation in response to climate change cause a subsurface decline of up to 40% in viable habitat overlap of Antarctic toothfish with important prey species, such as Antarctic silverfish and icefish. Acknowledging regional differences, our results demonstrate that warming and deoxygenation alone can significantly perturb predator-prey habitat overlap in the Southern Ocean. Our findings highlight the need for a better quantitative understanding of climate change impacts on Antarctic species to better constrain future ecosystem impacts of climate change.

Ocean temperature and dissolved oxygen shape marine habitats in an interplay with species' physiological characteristics. Therefore, the observed and projected warming and deoxygenation of the world's oceans in the 21st century may strongly affect species' habitats. Here, we implement an extended version of the Aerobic Growth Index (AGI), which quantifies whether a viable population of a species can be sustained in a particular location. We assess the impact of projected deoxygenation and warming on the contemporary habitat of 47 representative marine species covering the epipelagic, mesopelagic, and demersal realms. AGI is calculated for these species for the historical period and into the 21st century using bias-corrected environmental data from six comprehensive Earth system models. While habitat viability decreases nearly everywhere with global warming, the impact of this decrease is strongly species dependent. Most species lose less than 5 % of their contemporary habitat volume at 2 ∘C of global warming relative to preindustrial levels, although some individual species are projected to incur losses 2–3 times greater than that. We find that the in-habitat spatiotemporal variability of O2 and temperature (and hence AGI) provides a quantifiable measure of a species' vulnerability to change. In the event of potential large habitat losses (over 5 %), species vulnerability is the most important indicator. Vulnerability is more critical than changes in habitat viability, temperature, or pO2 levels. Loss of contemporary habitat is for most epipelagic species driven by the warming of ocean water and is therefore elevated with increased levels of global warming. In the mesopelagic and demersal realms, habitat loss is also affected by pO2 decrease for some species. Our analysis is constrained by the uncertainties involved in species-specific critical thresholds, which we quantify; by data limitations on 3D species distributions; and by high uncertainty in model O2 projections in equatorial regions. A focus on these topics in future research will strengthen our confidence in assessing climate-change-driven losses of contemporary habitats across the global oceans.

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.

Although both physical and biological marine changes are required to explain the 100 ppm lower atmospheric pCO2 of the Last Glacial Maximum (LGM, ∼21 ka) as compared to preindustrial (PI) times, their exact contributions are debated. Proxies of past marine carbon cycling (such as δ13C) document these changes and thus provide constraints for quantifying the drivers of long-term carbon cycle variability. This modeling study discusses the physical and biological changes in the ocean needed to simulate an LGM ocean in satisfactory agreement with proxy data, here focusing especially on δ13C. We prepared a PI and LGM equilibrium simulation using the ocean model NorESM-OC with full biogeochemistry (including the carbon isotopes δ13C and radiocarbon) and dynamic sea ice. The modeled LGM–PI differences are evaluated against a wide range of physical and biogeochemical proxy data and show agreement for key aspects of the physical ocean state within the data uncertainties. However, the lack of a simulated increase of regenerated nutrients for the LGM indicates that additional biogeochemical changes are required to simulate an LGM ocean in agreement with proxy data. In order to examine these changes, we explore the potential effects of different global mean biological pump efficiencies on the simulated marine biogeochemical tracer distributions. Through estimating which biological pump efficiency reduces LGM model–proxy biases the most, we estimate that the global mean biological pump efficiency increased from 38 % (PI) to up to 75 % (LGM). The drivers of such an increase in the biological pump efficiency may be both biological and related to circulation changes that are incompletely captured by our model – such as stronger isolation of Southern Source Water. Finally, even after considering a 75 % biological pump efficiency in the LGM ocean, a remaining model–proxy error in δ13C exists that is 0.07 ‰ larger than the 0.19 ‰ data uncertainty. This error indicates that additional changes in ocean dynamics are needed to simulate an LGM ocean in agreement with proxy data.

The dominant pacing of glacial‐interglacial cycles in deep‐ocean δ18O records changed substantially during the Mid‐Pleistocene Transition. The precessional cycle (∼23 ky) is absent during the Early Pleistocene, which we show can be explained by cancellation of the hemispherically antiphased precessional cycle in the Early Pleistocene interior ocean. Such cancellation develops due to mixing of North Atlantic and Southern Ocean δ18O signals at depth, and shows characteristic spatial patterns. We explore the cancellation potential for different North Atlantic and Southern Ocean deep‐water source δ18O values using a tracer transport ocean model. Cancellation of precession occurs for all signal strengths and is widespread for a signal strength typical for the Early Pleistocene. Early Pleistocene precessional power is therefore likely incompletely archived in deep‐sea δ18O records, concealing the true periodicity of the glacial cycles in the two hemispheres.

The ocean carbon cycle is a key player in the climate system through its role in regulating the atmospheric carbon dioxide concentration and other processes that alter the Earth's radiative balance. In the second version of the Norwegian Earth System Model (NorESM2), the oceanic carbon cycle component has gone through numerous updates that include, amongst others, improved process representations, increased interactions with the atmosphere, and additional new tracers. Oceanic dimethyl sulfide (DMS) is now prognostically simulated and its fluxes are directly coupled with the atmospheric component, leading to a direct feedback to the climate. Atmospheric nitrogen deposition and additional riverine inputs of other biogeochemical tracers have recently been included in the model. The implementation of new tracers such as “preformed” and “natural” tracers enables a separation of physical from biogeochemical drivers as well as of internal from external forcings and hence a better diagnostic of the simulated biogeochemical variability. Carbon isotope tracers have been implemented and will be relevant for studying long-term past climate changes. Here, we describe these new model implementations and present an evaluation of the model's performance in simulating the observed climatological states of water-column biogeochemistry and in simulating transient evolution over the historical period. Compared to its predecessor NorESM1, the new model's performance has improved considerably in many aspects. In the interior, the observed spatial patterns of nutrients, oxygen, and carbon chemistry are better reproduced, reducing the overall model biases. A new set of ecosystem parameters and improved mixed layer dynamics improve the representation of upper-ocean processes (biological production and air–sea CO2 fluxes) at seasonal timescale. Transient warming and air–sea CO2 fluxes over the historical period are also in good agreement with observation-based estimates. NorESM2 participates in the Coupled Model Intercomparison Project phase 6 (CMIP6) through DECK (Diagnostic, Evaluation and Characterization of Klima) and several endorsed MIP simulations.

Sea surface salinity restoringWe compare the ocean circulation after 1000 years of simulation in a LGM spinup with and without extra adjustments to the sea surface salinity (SSS) reference field used for salinity restoring.The adjusted spinup is the one used in our study and described in the main text.If unadjusted, the LGM SSS reference field is derived from the PI SSS reference field through the addition of a PMIP3-based sea surface salinity anomaly (Morée and Schwinger, 2019a; Morée and Schwinger, 2019b) (Fig. S1a).For the unadjusted setup NorESM-OC simulated an AMOC of 21 Sv (17 Sv for the adjusted spinup), and a Drake Passage through flow of 112 Sv (134 Sv for the adjusted spinup).The SSW profoundly decrease their volume in the Atlantic (Fig. S2).As especially the retreat of SSW is in disagreement with proxy-based reconstructions (see main text), we applied an adjustment to the SSS relaxation of -0.5 psu in the North Atlantic in the region where the anomaly was largest, as well as in the Southern Ocean (Fig. S1a).The simulated SSS for the adjusted run is shown in Fig. S1b. Figure S11 Simulated LGM O2 concentrations for a 75 % ̅̅̅̅̅̅̅̅ , with a zero contour line.Atlantic Pacific

Deep-sea δ18O records show a pronounced difference in Milankovitch periodicity between the Early and Late Pleistocene. δ18O is interpreted as a proxy for ice sheet volume and temperature, which led to the conclusion that glacial-interglacial cycles considerably changed their rhythm during the mid-Pleistocene. This transition is referred to as the mid-Pleistocene Transition (MPT). Specifically, the precessional component of the Milankovitch cycles is absent in Early Pleistocene δ18O records, despite its continuous presence in solar insolation forcing to the ice sheets. Climate feedbacks involving (sea) ice, geological processes and carbon and nutrient cycling have been proposed as causes of this marked change. We however show that the absence of an Early Pleistocene precession signal in deep-sea δ18O records could be the result of destructive interference of the precessional cycle in the interior ocean. Such cancellation is caused by the anti-phasing of the precessional cycle between the North Atlantic and Southern Ocean deep-water sources (see Figure). We explore the potential for cancellation in the transient setup of the Total Matrix Intercomparison model for a wide range of source signal strengths. Our results show that cancellation can cause the absence of the precessional signal due to cancellation in the interior South-Atlantic, Indian and Pacific basins. Cancellation is especially widespread for a relative end-member contribution typical for the Early Pleistocene. We therefore conclude that the precessional component is likely incompletely archived in Early Pleistocene δ18O records, and appears as an actual change in Milankovitch periodicity across the MPT. Proxies not susceptible to cancellation of precession (such as those currently retrieved across the MPT from Antarctica) would be able to verify to what extent deep-sea δ18O correctly represents Pleistocene climate.

The authors discussed about the glacial changes in d13C distribution in the ocean by comparing LGM ocean (NorESM-OC model) simulations with proxy data. The model significantly underestimates the glacial d13C changes compared with the proxy data; for example, negative signal of d13C in the deep Atlantic Ocean inferred from the proxy data is not reproduced in the model. At the same time, the model shows the decrease of the ocean biological pump efficiency in the LGM (33%) compared with the PI (38%), opposite to the fact that this is believed to be increased from the proxy data. The authors discussed the response of d13C by artificially increasing the ocean biological pump efficiency. The authors concluded that an approximate doubling of the global

Model simulations of the Last Glacial Maximum (LGM; ∼ 21 000 years before present) can aid the interpretation of proxy records, can help to gain an improved mechanistic understanding of the LGM climate system, and are valuable for the evaluation of model performance in a different climate state. Ocean-ice only model configurations forced by prescribed atmospheric data (referred to as “forced ocean models”) drastically reduce the computational cost of palaeoclimate modelling compared to fully coupled model frameworks. While feedbacks between the atmosphere and ocean and sea-ice compartments of the Earth system are not present in such model configurations, many scientific questions can be addressed with models of this type. Our dataset supports simulations of the LGM in a forced ocean model set-up while still taking advantage of the complexity of fully coupled model set-ups. The data presented here are derived from fully coupled palaeoclimate simulations of the Palaeoclimate Modelling Intercomparison Project phase 3 (PMIP3). The data are publicly accessible at the National Infrastructure for Research Data (NIRD) Research Data Archive at https://doi.org/10.11582/2020.00052 (Morée and Schwinger, 2020). They consist of 2-D anomaly forcing fields suitable for use in ocean models that employ a bulk forcing approach and are optimized for use with CORE forcing fields. The data include specific humidity, downwelling long-wave and short-wave radiation, precipitation, wind (v and u components), temperature, and sea surface salinity (SSS). All fields are provided as climatological mean anomalies between LGM and pre-industrial (PI) simulations. These anomaly data can therefore be added to any pre-industrial ocean forcing dataset in order to obtain forcing fields representative of LGM conditions as simulated by PMIP3 models. Furthermore, the dataset can be easily updated to reflect results from upcoming and future palaeo-model intercomparison activities.

δ13C, the standardised 13C ∕ 12C ratio expressed in per mille, is a widely used ocean tracer to study changes in ocean circulation, water mass ventilation, atmospheric pCO2, and the biological carbon pump on timescales ranging from decades to tens of millions of years. δ13C data derived from ocean sediment core analysis provide information on δ13C of dissolved inorganic carbon and the vertical δ13C gradient (i.e. Δδ13C) in past oceans. In order to correctly interpret δ13C and Δδ13C variations, a good understanding is needed of the influence from ocean circulation, air–sea gas exchange and biological productivity on these variations. The Southern Ocean is a key region for these processes, and we show here that Δδ13C in all ocean basins is sensitive to changes in the biogeochemical state of the Southern Ocean. We conduct a set of idealised sensitivity experiments with the ocean biogeochemistry general circulation model HAMOCC2s to explore the effect of biogeochemical state changes of the Southern and Global Ocean on atmospheric δ13C, pCO2, and marine δ13C and Δδ13C. The experiments cover changes in air–sea gas exchange rates, particulate organic carbon sinking rates, sea ice cover, and nutrient uptake efficiency in an unchanged ocean circulation field. Our experiments show that global mean Δδ13C varies by up to about ±0.35 ‰ around the pre-industrial model reference (1.2 ‰) in response to biogeochemical change. The amplitude of this sensitivity can be larger at smaller scales, as seen from a maximum sensitivity of about −0.6 ‰ on ocean basin scale. The ocean's oldest water (North Pacific) responds most to biological changes, the young deep water (North Atlantic) responds strongly to air–sea gas exchange changes, and the vertically well-mixed water (SO) has a low or even reversed Δδ13C sensitivity compared to the other basins. This local Δδ13C sensitivity depends on the local thermodynamic disequilibrium and the Δδ13C sensitivity to local POC export production changes. The direction of both glacial (intensification of Δδ13C) and interglacial (weakening of Δδ13C) Δδ13C change matches the direction of the sensitivity of biogeochemical processes associated with these periods. This supports the idea that biogeochemistry likely explains part of the reconstructed variations in Δδ13C, in addition to changes in ocean circulation.