Recent studies investigating future warming scenarios have shown that the ocean oxygen content will continue to decrease over the coming century due to ocean warming and changes in oceanic circulation. However, significant uncertainties remain regarding the magnitude and patterns of future ocean deoxygenation. Here, we simulate ocean oxygenation with the ACCESS ESM1.5 model during two past interglacials that were warmer than the preindustrial (PI) climate, the Last Interglacial (Marine Isotope Stage (MIS) 5e, ∼ 129–115 ka) and MIS 9e (∼ 336–321 ka). While orbital parameters were similar during MIS 5e and MIS 9e, with lower precession, higher eccentricity, and higher obliquity than PI, greenhouse gas radiative forcing was highest during MIS 9e and lowest during MIS 5e. We find that the global ocean is overall less oxygenated in the MIS 5e and MIS 9e simulations compared to the PI control run and that oxygen concentrations are more sensitive to changes in the distribution of incoming solar radiation than to differences in greenhouse gas concentrations. Large regions in the Mediterranean Sea are hypoxic in the MIS 5e simulation, and to a lesser extent in the MIS 9e simulation, due to an intensification and expansion of the African monsoon, enhanced river runoff and resulting freshening of surface waters and stratification. Upwelling zones off the coast of North America and North Africa are weaker in both simulations compared to the preindustrial control run, leading to less primary productivity and export production. Antarctic Bottom Water is less oxygenated, while North Atlantic Deep Water and the North Pacific Ocean at intermediate depths are higher in oxygen content. All changes in oxygen concentrations are primarily caused by changes in ocean circulation and export production and secondarily by changes in temperature and solubility.
Australian Research CouncilDP180100048SR200100008CE170100023Australian GovernmentNational Environmental Science Program (NESP)Introduction
The global ocean oxygen inventory has decreased by over 2 % since 1960, and the volume of anoxic waters has more than quadrupled over the same time period . This decrease was caused by a decline in oxygen solubility due to higher ocean temperatures , but changes in ocean circulation and biological consumption also contributed to the observed changes . Model simulations project that this decline will continue into the future, although the spatial patterns and magnitude are model-dependent and strongly influenced by the models' ocean diffusivity parameters . Future ocean deoxygenation could eventually pose severe problems for ocean ecosystems and human societies , and there has been a growing interest in recent years to better quantify ocean deoxygenation and understand its drivers.
Ocean oxygenation has also changed in the past. For example, find sulfidic conditions in the near-surface sediments of the Peruvian upwelling area during the Last Interglacial (LIG) period (Marine Isotope Stage (MIS) 5e, ∼ 129–116 kyr ago), based on sedimentary molybdenum accumulation. Records of sedimentary δ15N and redox-sensitive metals show that interglacials in the late Quaternary led to an expansion and/or intensification of near-surface and intermediate-depth suboxic zones in the eastern Pacific margins and the Arabian Sea compared to glacials . These changes in oxygenation have been generally assigned to changes in oxygen supply to the global thermocline . For large regions of the deep ocean this trend was inverted , and deep and bottom waters were better ventilated during interglacials compared to glacials. report close to suboxic conditions during glacial intervals of the past 150 kyr compared to better ventilated conditions during MIS 5e and the Holocene (the past 11 700 years) in the deep subarctic Pacific. In the deep equatorial Pacific there is multi-proxy evidence for reduced oxygen concentrations in all deep Pacific Ocean water masses below 1 km during glacials compared to MIS 5e and the Holocene . Bottom waters in the northern Arabian Sea, on the Portuguese margin, and in the North and South Atlantic, as well as South Atlantic Circumpolar Deep Water (CDW), were also less oxygenated during glacial episodes compared to MIS 5e and the Holocene . It is noteworthy that most reconstructions of ocean oxygenation in the late Quaternary compare glacials with interglacials; there is very little evidence on how oxygen concentrations varied between interglacials. This is an important caveat, as better understanding of ocean oxygen variability in the past, and, in particular, ocean oxygenation during past warm episodes, could inform on the main processes influencing current ocean deoxygenation.
One region with notably large changes in past ocean oxygen concentrations is the semi-enclosed Mediterranean Sea, where there is evidence of intervals with severe anoxia during past interglacials . These intervals are characterized by sediment layers with elevated organic carbon concentrations (sapropels) and coincide with astronomically timed episodes of monsoon intensification, causing enhanced runoff from North Africa into the Mediterranean Sea and leading to stratification . The resulting anoxia was often more intensely developed in the eastern Mediterranean than in the western Mediterranean.
Here, we analyze simulations integrated with the Australian earth system model ACCESS-ESM1.5 for two interglacials, MIS 5e, or the Last Interglacial (LIG), and MIS 9e. The LIG is an interesting time period because it was globally the warmest interglacial of the past 800 kyr . Being the most recent interglacial, the spatial and temporal resolution of climate proxies is also much better than for earlier interglacials. The LIG climate has been recently simulated by a variety of climate models under the Paleoclimate Model Intercomparison Project 4 (PMIP4) lig127k experiment .
MIS 9e (∼ 336–321 kyr ago) stands out by its high greenhouse gas forcing, with the highest carbon dioxide and methane concentrations over the past 800 kyr . It is also one of the three interglacials (together with MIS 11 and MIS 17) with the lowest value of seawater δ18O, which might imply high sea levels and small continental ice volumes , although sea-level records from the Red Sea show that MIS 5e had higher sea levels than MIS 9e . MIS 9e was an exceptionally short interglacial with the warmest temperatures in Antarctica over the past 800 kyr . We will focus on how changes in orbital parameters and greenhouse gas concentrations influenced simulated ocean oxygenation for these two interglacials.
This paper is structured as follows. Section describes the model, the experimental setup, and the time series for each of the analyzed experiments. Section first analyses changes in large-scale circulation patterns and oxygenation and then focuses on changes in the Mediterranean Sea. In Sect. we discuss uncertainties and situate our results within a broader context. Section summarizes the main results.
Methods
We use the Australian Community Climate and Earth System Simulator Earth System Model, ACCESS-ESM1.5 , to integrate three equilibrium simulations under preindustrial, MIS 5e, and MIS 9e boundary conditions. The model includes the atmosphere UK Met Office Unified Model (UM) version 7.3 , the Community Atmosphere Biosphere Land Exchange model (CABLE) version 2.4 , the NOAA/GFDL Modular Ocean Model (MOM) version 5 , and the Los Alamos National Laboratory sea ice model (LANL CICE) version 4.1 . The components communicate with each other through the Ocean Atmosphere Sea Ice Soil Model Coupling Toolkit (OASIS-MCT; ). The land and atmosphere components have a horizontal resolution of 1.875° × 1.25°, with 38 vertical levels for the atmosphere model. The ocean and sea ice components have a resolution of 1° × 1°, with 50 vertical levels for the ocean model. The horizontal resolution of the ocean model is higher near the Equator (0.33°) and in the Southern Ocean (∼ 0.4° at 70° S). Ocean biogeochemistry is represented by the Whole Ocean Model of Biogeochemistry And Trophic-dynamics (WOMBAT; ).
WOMBAT is a nutrient–phytoplankton–zooplankton–detritus (NPZD) model . It includes one functional type of phytoplankton and zooplankton. As prognostic tracers it simulates dissolved inorganic carbon (DIC), alkalinity (ALK), phosphate (PO4), oxygen (O2), and iron. The stoichiometry is fixed at a C : N : P : O2 ratio of 106 : 16 : 1 : -172. CaCO3 export from the photic zone is set at ∼ 8 % of the organic carbon export. Detrital decomposition is a function of temperature and is allowed to occur when oxygen is zero. Even though nitrification and denitrification are not explicitly included in the model and the global nitrogen budget is kept constant , this formulation emulates the effect of denitrification. The dissolution of CaCO3 occurs at a constant rate. All organic and inorganic particles reaching the bottom are remineralized, given that ACCESS-ESM1.5 does not include burial of sediments.
The preindustrial 1850 equilibrium simulation (PI) is integrated following the CMIP6 protocol with the exception of using the CMIP5 solar constant (1365.65 W m−2). The MIS 5e simulation (LIG) is integrated following the PMIP4 protocol with the solar constant of CMIP5-PMIP3 (1365.65 W m−2). Both PI and LIG simulations have been evaluated extensively in the recent literature .
Experimental setup. LIG boundary conditions follow the PMIP4 protocol . MIS 9e boundary conditions are based on for orbital parameters and peak concentrations for methane , carbon dioxide , and nitrous oxide . Please note that the first 372 years of the LIG simulation have erroneous forcing and were integrated with PI greenhouse gas concentrations.
PILIGMIS 9eOrbital parametersEccentricity0.0167640.0393780.03229469Obliquity (degrees)23.45924.04024.2397Perihelion -180100.33275.41297.994Vertical equinox21 March at noon21 March at noon21 March at noonGreenhouse gasesCarbon dioxide (ppm)284.3275298.6Methane (ppb)808.2685797Nitrous oxide (ppb)273255287.3PaleogeographyModernModernModernIce sheetsModernModernModernVegetationPIPIPIAerosolsCMIP Deck piControlCMIP Deck piControlCMIP Deck piControlIntegration length1000 years1823 years1750 yearsRemaining drift below 3000 m per 100 yearsPotential temperature (° C per 100 years)+0.004+0.040+0.055Dissolved O2 (mmol m−3 per 100 years)+0.38-0.85-1.28
The MIS 9e simulation was set up with boundary conditions corresponding to 333 ka and the same solar constant as for PI and LIG (see Table ). The Greenland and Antarctic ice sheets and vegetation distribution are the same in all three simulations, but the leaf area index is calculated prognostically. Anomalies in incoming solar radiation are shown in Fig. .
Time series of dissolved O2 concentration in mmol m−3 averaged (a) globally, (b) in the Mediterranean Sea, (c) between 500 and 1200 m depth and 40–60° N (North Atlantic Deep Water), (d) between 250 and 1200 m depth and 40–60° N in the North Pacific, and (e) at 3000 m depth between 60 and 35° S (Antarctic Bottom Water) for PI (solid lines), LIG (dotted lines), and MIS 9e (dashed lines) simulations. Red rectangles designate the last 100 years of each simulation that are analyzed in this paper. The first 372 years of LIG were integrated with PI greenhouse gas concentrations and have been deleted from our file servers; we show therefore interpolations (long dashed lines) over this time span.
Dissolved oxygen is a tracer with very long equilibration times. Figure shows time series of dissolved oxygen for the three simulations analyzed in this study. While oxygen concentrations have reached quasi-equilibrium for North Atlantic Deep Water (NADW), in the intermediate waters of the North Pacific, and in Antarctic Bottom Water (AABW) south of 35° S in all runs, there is still a slight drift in global mean dissolved oxygen in our LIG and MIS 9e simulations and a more substantial drift in the Mediterranean Sea.
In Sect. we partition changes in dissolved oxygen into two components, changes in the saturated concentration of oxygen O2sat and changes in apparent oxygen utilization (AOU). AOU estimates the oxygen consumed during respiration and can be calculated as the difference between dissolved oxygen concentration (O2) and O2sat:
AOU=O2sat(T,S)-O2,
where T is the potential temperature and S salinity. Changes in AOU are therefore a combination of changes in circulation (with sluggish water masses tending to have higher AOU) and changes in remineralization rates, which depend on the vertical export of organic matter (export production) and temperature. Please note that the here used metric AOU assumes that dissolved oxygen in surface waters is in equilibrium with the atmosphere, which might lead to an overestimation of the true oxygen utilization (TOU) .
ResultsLarge-scale oxygenation
All results reported in this section are based on 100-year means, calculated over the last 100 years of each simulation (see rectangles in Fig. ). The global temperature patterns simulated by ACCESS ESM1.5 under LIG boundary conditions are described elsewhere . Figure shows annual-mean surface air temperature (SAT) anomalies and sea surface temperature (SST) anomalies for the simulations described in this study. Due to different orbital parameters, and in particular a positive summer insolation anomaly at high northern latitudes and a positive spring anomaly at high southern latitudes (; Fig. ), the annual-mean surface air temperature in our LIG simulation is 1.11 and 1.40 °C higher north and south of 40° N/S, respectively, compared to PI. While the insolation during MIS 9e is similar to LIG , greenhouse gas concentrations are higher (Table ), with anomalies of ∼ 24 ppm for CO2, ∼ 112 ppb for CH4, and ∼ 32 ppb for N2O, thus leading to globally warmer conditions in our MIS 9e simulation compared to the LIG simulation (Fig. ). The annual-mean surface air temperature is 0.49 and 1.20 °C higher in the LIG and MIS 9e runs, respectively, compared to PI.
Annual-mean near-surface air temperature (SAT) anomalies (a, c, e) and sea surface temperature (SST) anomalies (b, d, f) in degrees Celsius (°C). Statistically insignificant differences (p value > 0.05) are shaded in gray. (a, b) LIG–PI, (c, d) MIS 9e–PI, (e, f) MIS 9e–LIG.
The PI simulation reproduces observed patterns of dissolved oxygen (World Ocean Atlas, 1965–2020) reasonably well (; Fig. c and d) but underestimates O2 in the eastern Pacific Ocean and southeastern Atlantic Ocean. It overestimates O2 in the Arabian Sea, the Southern Ocean, and the Arctic. The global mean ocean O2 concentration equals 176.3 mmol m−3 in our PI simulation compared to 158.0 and 154.5 mmol m−3 in the LIG and MIS 9e simulations, respectively. The largest differences between the simulations are located in the deep ocean (below 2500 m depth) and follow the pathways of Antarctic Bottom Water (AABW; Figs. and ). Overall, the anomalies in LIG and MIS 9e compared to PI are quite similar, with the loss in AABW O2 being slightly less pronounced in LIG (Figs. and ). Significant differences can also be seen in the Mediterranean Sea (Fig. , Sect. ).
Vertically averaged dissolved O2 concentration anomalies in mmol m−3. Statistically insignificant differences (p value > 0.01) are shaded in gray. (a) LIG–PI, (b) MIS 9e–PI, (c) vertically averaged dissolved O2 concentration for the PI simulation, and (d) vertically averaged O2 concentration from World Ocean Atlas (WOA, 1965–2022; ).
Zonally averaged dissolved O2 concentration in mmol m−3 (left), saturated O2 concentration (middle), and -1⋅AOU (apparent oxygen utilization, right). Statistically insignificant differences (p value > 0.01) are hatched and shaded in gray. Atlantic Ocean anomalies: LIG–PI (a–c) and MIS 9e–PI (d–f). Pacific Ocean anomalies: LIG–PI (g–i) and MIS 9e–PI (j–l). To facilitate the comparison between oxygen solubility and oxygen utilization, we multiplied AOU by -1.
The next two subsections analyze changes in the main large-scale bottom and deep water masses, AABW and NADW. We also see significant changes in the oxygenation of North Pacific Intermediate Water, which are described in Sect. . Changes in low-oxygen zones can have a large impact on marine life and biogeochemical cycles. We therefore analyze changes in the simulated oxygen minimum zones in Sect. .
Antarctic Bottom Water – warmer and less ventilated
AABW is warmer and weaker in our LIG and MIS 9e simulations compared to PI (Fig. ). At 30° S, the maximum absolute value of the streamfunction reaches 3.8 Sv (sverdrups) in LIG and 2.9 Sv in MIS 9e compared to 7.9 Sv in PI (Fig. ). Bottom water convection sites in the Weddell, Lazarev, and Ross seas are smaller in spatial extent and are also less intense (Fig. ). This weakening of deep-ocean convection is mainly due to reduced sea ice formation and leads to warmer, less ventilated, and therefore less oxygenated AABW. The relative contribution of changes in apparent oxygen utilization (AOU) is higher than the contribution of changes in solubility (compare middle and right columns in Fig. ). Export production is enhanced in the Southern Ocean in LIG and MIS 9e compared to PI (; Fig. ), and remineralization rates are higher due to higher temperatures . Changes in dissolved oxygen concentrations in AABW are therefore primarily due to the weakening of Antarctic Bottom Water circulation, together with the increase in export production and higher remineralization rates, and secondarily caused by the temperature-dependent solubility effect.
Meridional overturning streamfunction in sverdrups (Sv) for the Atlantic Ocean (a, c, e) and global (b, d, f). (a, b) LIG, (c, d) MIS 9e, and (e, f) PI.
Detrital organic carbon concentration and anomalies in µmol C m−3 at 250 m depth. Statistically insignificant differences (p value > 0.01) are hatched and shaded in gray. Regions that are shallower than 250 m are shown in white. (a) LIG–PI, (b) MIS 9e–PI, and (c) PI.
North Atlantic Deep Water – colder and better ventilated
Figure a and d show that North Atlantic Deep Water (NADW) is more oxygenated in both MIS 9e and LIG compared to PI, with MIS 9e being more oxygenated than LIG at mid-depth between 40° S and 20° N. This is due to enhanced deepwater formation (; Fig. ) and stronger advection of NADW at depth leading to a decrease in AOU (Fig. c and f) and to lower temperatures of NADW in LIG and MIS 9e compared to PI (Figs. and , middle column). The maximum value of the streamfunction in the Atlantic at 45° N reaches 23.6 Sv in our LIG simulation and 24.3 Sv in the MIS 9e simulation, compared to 19.9 Sv in the PI simulation. These simulated changes in NADW transport are associated with an advective–convective feedback mechanism involving the Subpolar Gyre (SPG), first described by . The SPG is stronger in our MIS 9e and LIG simulations than in PI (Fig. , left column), with an establishment of a third region of deep convection in the Labrador Sea in addition to the deep convection sites in the Norwegian Sea and south of the Greenland–Scotland Ridge (Fig. ). During boreal winter (December to February), a stronger SPG results in stronger outcropping of isopycnals and in higher sea surface salinities (Fig. , right column) due to enhanced transport of subtropical saline waters. Both are preconditioning the surface waters for convection. The resulting horizontal surface density gradients maintain the strong gyre. The new convection site in the Labrador Sea contributes colder waters to the resulting North Atlantic Deep Water mass, leading to an overall colder NADW compared to PI (Fig. ). Export production is enhanced over most of the North Atlantic Ocean in LIG and MIS 9e compared to PI and is therefore not one of the drivers of the increase in dissolved oxygen (Fig. ).
North Pacific Intermediate Water – warmer and better ventilated
The North Pacific ocean above ∼ 2000 m is significantly more oxygenated in our LIG and MIS 9e simulations than in our PI simulation (Fig. g and j), with LIG being even more oxygenated than MIS 9e. The potential temperatures of the affected water masses are higher in both simulations compared to PI (Figs. and h and k). The oxygenation is therefore solely due to changes in AOU, which shows strong negative anomalies (Fig. i and l). We can see a slight increase in ventilation during winter months in both simulations due to higher subduction rates in the western Bering Sea – along the coast of the Kamchatka Peninsula (Fig. ). This leads to younger and better ventilated waters at 250 m depth circulating across the Pacific and then southward following the California Current (shown in Fig. ). There is also a significant decrease in export production off the west coast of North America in LIG and MIS 9e compared to PI (Fig. ). This is due to a weakening of the north-westerly winds causing a weakening of the coastal upwelling zones (Fig. , left column) contributing to the gain in oxygen at deeper layers (Fig. a and b). Furthermore, better ventilated subsurface waters can lead to a decrease in nutrient concentrations and therefore an overall reduction in nutrients transported to the surface via upwelling .
Oxygen minimum zones (OMZs) – no significant change
Figure (left column) shows dissolved O2 concentrations in hypoxic zones at 300 m depth. In our PI simulation, the model represents the large-scale patterns reasonably well in the tropical and subtropical Pacific but overestimates oxygen depletion within these zones (Fig. e and g). The simulated OMZ in the Atlantic ocean off the African coast is also overestimated in extent and intensity. Contrarily, oxygen depletion in the North Pacific and the Arabian Sea is underestimated at 300 m depth, but the North Pacific ocean becomes hypoxic at deeper levels (Fig. , right column).
O2 concentration and anomalies at 300 m depth in hypoxic zones (left) and vertical minimum of O2 concentration and anomalies in hypoxic zones (right) in (a, b) LIG–PI, (c, d) MIS 9e–PI, (e, f) PI, and (g, h) World Ocean Atlas (WOA, 1965–2022). Black contour lines in panels (a)–(d) indicate the 62 mmol m−3 isolines for PI (solid), LIG (dotted), and MIS 9e (dashed). Statistically insignificant differences (p value > 0.01) are shaded in gray. Hypoxic zones are defined as zones where O2 concentration at 300 m is below 62 mmol m−3 for panels (a), (c), (e), and (g) and where the vertical minimum of O2 concentration is below 62 mmol m−3 for panels (b), (d), (f), and (h).
The extent and intensity of OMZs is very similar between the three simulations at 300 m depth (Fig. a and c). The Eastern Boundary Upwelling System off the coast of North America is weakened in our LIG and MIS 9e simulations compared to PI (Fig. ), which leads to a decrease in nutrient availability at the surface (not shown), causing lower primary productivity and export production (Fig. ). Consequently, AOU is significantly reduced in our LIG and MIS 9e simulations off the coast of California and Baja California (Fig. b and d), leading to an increase in oxygen in the northern hemispheric tropical Pacific OMZ that is partially compensated by a decrease in solubility due to higher water temperatures (Fig. a and c). The Eastern Boundary Upwelling System off Africa is also weakened (Fig. , right column) but does not affect the OMZ, which is situated further south. The OMZs in the southern hemispheric tropical Pacific and South Atlantic are slightly intensified in our LIG and MIS 9e simulations compared to PI, which is due to warmer waters and partially compensated by a decrease in AOU in these regions (Fig. ).
The right column of Fig. shows the minimum of dissolved oxygen in the water column. The model represents the extent of hypoxic waters well in the Pacific Ocean but underestimates the extent in the Indian Ocean and, again, misses the OMZ in the Arabian Sea. The extent of the South Atlantic OMZ is overestimated. The intensity of oxygen loss is overestimated in the eastern regions of all OMZs (Fig. f and h). Differences between the simulations show that the mid- to high-latitude North Pacific is significantly more oxygenated in LIG and, to a lesser extent, in MIS 9e, as discussed in Sect. . The southern and western parts of the Pacific hypoxic zones as well as most of the hypoxic zones in the Atlantic Ocean become more oxygen depleted in LIG and MIS 9e compared to PI. The minimum oxygen concentration in the water column shifts to deeper depths in the northern hemispheric Pacific and to shallower depth in most of the southern hemispheric Pacific (Fig. ).
Total volume (percentage) of oxygen-depleted water masses.
The global ocean is overall less oxygenated in MIS 9e and LIG compared to PI, and this is also reflected in the total volume of hypoxic waters (< 62 mmol m−3), which is 5 % larger in these simulations (Table ). However, the suboxic (< 10 mmol m−3) and anoxic (< 1 mmol m−3) zones are smaller in our LIG and MIS 9e simulations, mainly due to the oxygenation of the North Pacific Ocean (Sect. ).
Oxygenation of the Mediterranean Sea – large-scale hypoxia
The Mediterranean Sea is a semi-enclosed sea which is linked to the Atlantic Ocean via the relatively narrow and shallow (< 900 m) Strait of Gibraltar. Changes in large-scale circulation therefore impact the Mediterranean Sea only moderately, although the opposite is not true, as the Mediterranean outflow impacts the circulation in the Atlantic Ocean . There have been numerous time intervals during past interglacials when some regions of the Mediterranean Sea became anoxic, prompting us to analyze changes in simulated dissolved O2 in the Mediterranean Sea in this section.
Volume (percentage) of oxygen-depleted water masses in the Mediterranean Sea.
The mean O2 concentration in the Mediterranean Sea equals 218.7 mmol m−3 in our PI simulation compared to 106.2 and 61.6 mmol m−3 in our MIS 9e and LIG simulations, respectively. In the LIG simulation 63 % of Mediterranean waters are hypoxic, with dissolved oxygen concentrations below 62 mmol m−3, and 27 % are anoxic (less than 1 mmol m−3). The MIS 9e simulation is not quite as depleted in oxygen in the Mediterranean Sea, with 34 % of the total volume being hypoxic and 5.6 % anoxic (Table ).
Zonal and meridional mean vertical profiles in the Mediterranean Sea for PI (solid lines), LIG (dotted lines), and MIS 9e (dashed lines). (a) Potential temperature in degrees Celsius (°C), (b) salinity in practical salinity units (psu), (c) O2 concentration in mmol m−3, and (d) seawater age since last surface contact in years. Note that the maximum possible water age equals the length of the simulation and might therefore underestimate the water age for very old water masses.
Vertical mean O2 concentration (a, c, e) and meridional mean O2 concentration (b, d, f) in the Mediterranean Sea in mmol m−3 for PI (a, b), LIG (c, d), and MIS 9e (e, f). Symbols represent core sites discussed in the text.
The depletion of dissolved oxygen in the MIS 9e and LIG simulations is greatest in intermediate water masses (Figs. c and , right column), with horizontally averaged O2 concentrations dropping to 50 mmol m−3 in MIS 9e and reaching anoxic values in LIG between 500 and 1000 m depth (Fig. c). This deoxygenation is caused by a strong pycnocline due to a significant freshening and warming of the surface layers (Fig. a and b) that prevents mixing between surface and deeper layers and increases the water ages of intermediate and deep water masses (Fig. d). While there is deep water formation in winter months in the Ionian Sea with seasonal mean winter ventilation depth reaching 489 m, and to a lesser extent in the southern Adriatic Sea and east of the Gulf of Lion (west of Corsica) in our PI simulation, winter-mean mixed layer depths do not exceed 97 and 106 m in our LIG and MIS 9e simulations, respectively (Fig. ).
Figure shows that the central Mediterranean Sea is severely oxygen depleted in both LIG and MIS 9e simulations, with O2 concentrations dropping to zero in the Tyrrhenian Sea below 100 m depth. This is in qualitative agreement with the appearance of organic-rich dark layers in core DED 87-08 (39°42′ N, 13°35′ E; 2965 m; ) during eastern Mediterranean Sapropel events S5 (121.5–128.3 ka; ) and S10 (∼ 332 ka; ). Core KET 80-04 (39°40′ N, 13°34′ E; 2909 m) also shows evidence of anoxic conditions during S5 in the Tyrrhenian Sea .
The eastern Mediterranean basin is also oxygen depleted, especially in our LIG simulation. The simulated oxygen depletion in the Levantine Sea is in agreement with core MD 84-641 (33°02′ N, 32°38′ E; 1375 m), which recorded a strong sapropel event (S5) at the LIG (4.52 % total organic carbon (TOC)) and a less strong sapropel event (S10) at MIS 9e (2.5 % TOC) . Our results are also in agreement with , who find water mass ages of 1030 (+820/-520 years) in the deep eastern Mediterranean Sea (compared to ∼ 100 years today) at the end of S5, and with , who find euxinic conditions during S5 in the Levantine Sea.
The Ionian Sea is better oxygenated than the Levantine Sea in our simulations, but conditions at intermediate depths are still hypoxic in the MIS 9e simulation, dropping to 52.4 mmol m−3 and reaching anoxic values in the LIG simulation. Core KC07 (20°53′ N, 36°34′ E; 3410 m) records rich organic layers during both S5 and S10 but is situated at a depth where our simulations show higher oxygen concentrations. It should be noted that oxygen concentrations in the deep waters in the Ionian Sea are still drifting in our simulations and might reach suboxic values if integrated for longer.
In the western Mediterranean, sapropels are less distinctive than in the eastern Mediterranean, but they can still be detected in terms of higher-than-background organic carbon content. ODP leg 161 drilled five sites throughout the western Mediterranean and found that sapropel TOC values decrease to the west and are lowest in the Alboran Sea . Cores 975 (38°53.795′ N, 04°30.595′ E; 2416 m), 976 (36°12.32′ N, 04°18.760′ W; 1108 m), 977 (36°01.907′ N, 01°957.319′ W; 1984 m), and 979 (35°43.427′ N, 03°12.353′ W; 1062 m) all record higher organic carbon layers between 122 and 128 ka and between 328 and 331 ka, and the S5 event is color-banded in site 975, while the S10 event is homogeneous . This east–west gradient is reflected in our simulations, where waters near the Strait of Gibraltar are better oxygenated than waters in the Alboran Sea.
Anomalies of precipitation in North Africa during boreal summer (JJA) in mm d−1 for (a) LIG–PI and (b) MIS 9e–PI; (c) seasonal freshwater discharge from all rivers into the Mediterranean Sea for PI (solid line), LIG (dotted line), and MIS 9e (dashed line) in 106 kg s−1 and (d) seasonal freshwater surface fluxes (precipitation + runoff - evaporation) into the Mediterranean Sea in 106 kg s−1 for PI (solid line), LIG (dotted line), and MIS 9e (dashed line). Statistically insignificant differences (p value > 0.01) are shaded in gray in panels (a) and (b). Monthly values are processed following paleo-seasonal adjustments with reference dates of 127 and 333 ka for LIG and MIS 9e, respectively. The x axis in panels (c) and (d) represents months according to the present-day calendar.
The North African monsoon expands significantly into the Sahara region in our LIG and MIS 9e simulations, hence causing an intensification of precipitation during the wet season in North Africa . Rainfall in the Sahara region increases locally by up to +12.0 mm d−1 in LIG and +10.9 mm d−1 in MIS 9e during boreal summer (JJA) compared to PI (Fig. a and b). This excess rainfall causes excess freshwater runoff that enters the Mediterranean basin mainly through the Nile Delta. Although North African monsoon intensity peaks between June and August, outflow into the Mediterranean Sea reaches its maximum in autumn between September and November (Fig. c and d). This time lag of about 3 months is caused by the river routing timescale embedded in the ACCESS-ESM1.5 model (Total Runoff Integrating Pathways (TRIP) river routing; ) and is similar to the observed timescale of river routing.
The annual-mean surface freshwater flux into the Mediterranean Sea (precipitation + river runoff - evaporation) is equal to -344 × 106 kg s−1 in our PI simulation compared to +248 × 106 and +130 × 106 kg s−1 in our LIG and MIS 9e simulations, respectively. These positive fluxes cause freshening of the surface layers, stratification, and eventually deoxygenation of the Mediterranean Sea in our simulations (Fig. b). Compared to MIS 9e, the North African monsoon intensification is more pronounced in the LIG simulation (Fig. a and b), leading to higher river runoff (Fig. c). Sea surface salinities decrease at a faster rate and lead to stratification and deoxygenation earlier than in the MIS 9e simulation (Fig. b). The difference in oxygenation between LIG and MIS 9e is thus mainly due to the difference in the length of time during which the Mediterranean Sea is stratified, although small changes in biological productivity and export production (not shown) also contribute. It should be noted that oxygen levels in the Mediterranean Sea are still drifting in both simulations and that the equilibrium values are lower than what is presented here.
Discussion
The oxygenation of the world's oceans is very similar in our LIG and MIS 9e simulations and quite different compared to our PI simulation. These differences are primarily due to circulation changes and changes in export production and secondarily due to the solubility effect. The large-scale ocean circulation patterns, including AABW, NADW, and the ventilation of the North Pacific Ocean, are therefore very sensitive to the latitudinal and seasonal distribution of incoming solar radiation in the ACCESS ESM1.5 and less sensitive to changes in greenhouse gas concentrations within the range of these three interglacials. The insolation anomalies are indeed quite similar for LIG and MIS 9e (Fig. ), while the greenhouse gas forcing is highest for MIS 9e and lowest for LIG, with PI being in the middle (Table ).
While quantitative records of past oxygenation are difficult to reconstruct with proxies , qualitative records exist. For the more recent Pleistocene period, most records looking at past ocean oxygenation have focused on glacial–interglacial variations; for summaries see and . used epibenthic foraminifera carbon isotopes to make inferences of intermediate water ventilation and oxygenation in the northeast Pacific and suggested that waters above 2000 m depth were less oxygenated during MIS 5e. If the proposed relationship between benthic foraminifera carbon isotopes and apparent oxygen utilization holds, then overall reduced carbon isotopes during MIS 5e would imply an overall decrease in global ocean oxygenation, agreeing with our simulations. There is also evidence for reduced AABW formation and a decrease in dissolved oxygen in the Southern Ocean based on redox elements .
Oxygen minimum zones (OMZs) have expanded and contracted repeatedly in the past in response to global changes in climate . To the authors' knowledge, there is no global modeling study analyzing the extent and magnitude of OMZs during MIS 5e and MIS 9e. Simulating OMZs with earth system models is very challenging, because low oxygen concentrations are the result of a subtle balance between two large opposing processes: the physical transport of oxygen-rich waters into the region of interest and the local flux of particulate matter and remineralization. Climate models are not very skilled in representing these dynamics, mostly because their resolution does not resolve equatorial dynamics well enough. Another problem is the simple representation of biology and geochemistry in the majority of current state-of-the-art models. The ACCESS ESM1.5 is not an exception. While the extent of the OMZs in the tropical Pacific is reasonably well represented, there are deficits in other regions. Our results show that relatively small changes in circulation, such as the subtle increase in ventilation of the upper North Pacific Ocean in our LIG and MIS 9e simulations, can have a significant impact on oxygen concentrations, emphasizing the need of a realistic representation of physical circulation when simulating oxygen. For example, cores LV-28-44-3 (52°02.514′ N, 153°05.949′ E; 684 m) in the eastern Sea of Okhotsk and MR0604-PC07A (51°16′56′′ N, 149°12′60′′ E; 1247 m) in the central Sea of Okhotsk recorded intervals of suboxic conditions during MIS 5e. While our simulation points to suboxic waters east of the Kamchatka Peninsula, the waters west of the peninsula are well oxygenated (Fig. ). It should be noted that these cores were retrieved from a semi-isolated marginal basin and are influenced by small-scale local circulation patterns that are not well resolved in the ACCESS ESM1.5.
We find that the main Eastern Boundary Upwelling Systems in the North Pacific and Atlantic Oceans are weaker in the LIG and MIS 9e simulations compared to PI. This is consistent with high sea surface temperatures and faunal composition off the Canary Islands during MIS 5e and microfossil assemblages from MIS 5e off the coast of northern California . However, Si / Ti, Cd / Al, and Ni / Al records from core MD02-2508 (23°27.91′ N, 111°35.74′ W; 606 m) suggest that biological productivity was high during MIS 5e further south, off the coast of Baja California , which is contrary to our results, although some upwelling is still present in our LIG and MIS 9e simulations. Coastal upwelling may have been similar or enhanced during MIS 5e and MIS 9e along the southeast African margin , expanded along the Brazil margin during MIS 5e , and increased at higher latitudes during MIS 5e and MIS 9e . This agrees broadly with the simulated patterns of detrital organic carbon concentration anomalies.
The Mediterranean's sedimentary record is punctuated by periodic deep-sea anoxic events that manifest themselves in layers with elevated organic carbon concentrations (sapropels) and are strongly associated with times of African monsoon intensification. We find monsoon intensification in our LIG and MIS 9e simulations, enhanced river runoff, and a freshening of the surface layers of the Mediterranean Sea. The resulting vertical stratification leads to anoxic conditions in the east Mediterranean Sea in our LIG simulation and hypoxic conditions in our MIS 9e simulation. The Tyrrhenian Sea is anoxic in both simulations, and the west Mediterranean shows an east–west gradient of oxygenation that is qualitatively in agreement with sediment data. It should be noted, however, that the ACCESS ESM1.5 is a global model with relatively coarse horizontal resolution (1° × 1° in the Mediterranean Sea). Such coarse-resolution model cannot be expected to simulate Mediterranean circulation skillfully. There are three main sites of deep water formation in the Mediterranean Sea in today's climate. Western Mediterranean deep water is formed during winter in the Gulf of Lion. Eastern Mediterranean deep water formation occurs in two sites, the Adriatic Sea, where bottom water is formed in winter and exits through the Strait of Otranto, and the Rhodes Gyre in the Levantine Basin. In our PI simulation, the main deep water formation occurs in the Ionian Sea, and to a lesser extent in the southern Adriatic Sea and east of the Gulf of Lion. Our results should therefore be taken as a qualitative indication of enhanced stratification and oxygen loss without putting too much trust into simulated regional changes. For example, find shallow water euxinia during S5 on the Adriatic Shelf due to a shutdown of North Adriatic Dense Water formation. This region remains well ventilated in our simulations. In addition, note that Mediterranean oxygen content is still drifting in our LIG and MIS 9e runs. As a result, MIS 9e could reach an oxygen loss similar to LIG on a longer timescale (Fig. ).
Our study shows that even relatively small changes in boundary conditions can lead to large changes in ocean circulation, upwelling systems, export production, and ocean oxygenation. While our results cannot directly inform on future changes in ocean oxygenation, there might be some similarities. Ocean temperatures and ocean stratification are projected to continue to increase until atmospheric CO2 concentrations finally plateau. Current ocean deoxygenation is therefore not easily reversible and will persist for centuries . There will be physiological and morphological impacts on organisms, including reduced growth for a vast range of taxonomic groups . Exposure to low-oxygen conditions has also been associated with a delay in when fish produce eggs, a reduction of the number of eggs fish produce, and blindness . The metabolic demand of oxygen increases with water temperature, and when combined with deoxygenation, this can lead to respiratory distress, followed by respiratory failure and death . Over 50 mass mortality events due to hypoxia have been recorded in the tropics to date . The consequences of deoxygenation for fisheries and the world's future food supply could thus be serious .
Here we perform a snapshot experiment of the peak LIG and MIS 9e climates with a single model. These simulations thus only represent an equilibrium response to peak interglacial conditions and do not represent potential transient climatic changes associated with the end of deglaciations or abrupt climate changes linked to continental ice-sheet melting. In addition, the dissolved O2 content simulated in these experiments is very dependent on the simulated oceanic circulation. The state of the ocean circulation during these interglacials is, however, unclear. Transient simulations through whole interglacial and multi-model intercomparisons would be preferable and give a better indication of the robustness and uncertainties of past oxygenation patterns and their variability . We hope that other modeling groups will consider analyzing ocean oxygenation during past interglacials. Oxygen concentration is a variable with extremely long equilibration times, so caution should be taken when climate models are not integrated for at least 1000 years, preferably longer, and the drift in the deep ocean is still significant.
Conclusions
We integrated three equilibrium simulations with Australia's earth system model ACCESS ESM1.5 under preindustrial (PI), Last Interglacial (LIG, 127 ka), and MIS 9e (333 ka) boundary conditions. The LIG and MIS 9e simulations show similar anomalies in large-scale oxygenation pointing to the fact that circulation patterns and oxygen concentrations are more sensitive to the distribution of incoming solar radiation than to greenhouse gas concentrations within the range of these three interglacials. Antarctic Bottom Water (AABW) is weaker and warmer, leading to deoxygenation mostly due to higher oxygen utilization and, to a lesser extent, the solubility effect. North Atlantic Deep Water (NADW) is overall cooler and better ventilated, leading to higher in situ oxygen concentrations. Water masses in the upper 2000 m of the North Pacific are both warmer and higher in oxygen content, due to enhanced subduction in the Bering Sea and a reduction in export production off the coast of North America. These anomalies are more pronounced in the MIS 9e simulation than in the LIG simulation, with the exception of the North Pacific Ocean, where the oxygenation anomaly is strongest in the LIG simulation. While the global ocean is overall less oxygenated in MIS 9e and LIG compared to PI, the volume of suboxic and anoxic waters is smaller in MIS 9e and LIG, mainly due to the oxygenation of the North Pacific Ocean. The Mediterranean Sea is the exception, where oxygen is significantly more depleted in MIS 9e and LIG compared to PI, due to an intensification and expansion of the African monsoon, enhanced runoff, and resulting freshening of surface waters and stratification. In the LIG simulation, 30 % of the total seawater volume in the Mediterranean Sea is anoxic (< 1 mm m−3) and 63 % is hypoxic (< 62 mm m−3). This deoxygenation is not quite as pronounced in the MIS 9e simulation with 6 % of anoxic waters and 34 % of hypoxic waters. Dissolved oxygen concentrations are still drifting in both simulations, even after integration times of well over 1500 years.
The simulation of ocean oxygen concentrations is challenging, as oxygen levels depend on complex physical and biogeochemical processes that must be represented correctly in a climate model to simulate oxygen realistically. State-of-the-art earth system climate models systematically underestimate the observed rates of oxygen loss and are also not very skillful in reproducing the observed patterns of deoxygenation . In this study, we confirmed that subtle changes in ocean circulation can have significant impact on oxygen concentrations. For trustworthy future projections, high-resolution ocean models are therefore needed to represent the relevant ocean circulation patterns as well as possible, but these models are computationally expensive. Oceanic oxygen takes centuries to adjust, and we currently lack computer power to run high-resolution models for the timescales needed for initialization and equilibrium responses. Another important challenge is the representation of biogeochemical processes in models, as our knowledge of coupled biogeochemical processes is still incomplete.
Insolation anomalies (W m−2) across latitudes and days of year. (a) LIG–PI, (b) MIS 9e–PI, and (c) MIS 9e–LIG. The bottom x axis represents the day of year, and the top x axis represents months according to the present-day calendar.
Zonally averaged potential temperature in degrees Celsius (°C) (g, h) and anomalies (a–f) in the Atlantic Ocean (left) and in the Pacific Ocean (right). Statistically insignificant differences (p value > 0.01) are hatched and shaded in gray. (a, b) LIG–PI, (c, d) MIS 9e–PI, (e, f) MIS 9e–LIG, and (g, h) PI.
Annual-mean mixed layer depth in meters (m). (a) PI, (b) LIG, and (c) MIS 9e.
Depth-integrated barotropic stream function in the North Atlantic Ocean averaged during boreal winter in sverdrups (Sv) (DJF; a–c) and sea surface salinity in practical salinity units (psu) (f) and anomalies (d, e) in the North Atlantic Ocean averaged over boreal winter months (DJF). Statistically insignificant differences (p value > 0.01) are shaded in gray. Mean over DJF months is processed following paleo-seasonal adjustments with reference dates of 127 and 333 ka for LIG and MIS 9e, respectively. (a) PI, (b) LIG, (c) MIS 9e, (d) LIG–PI, (e) MIS 9e–PI, and (f) PI.
Mixed layer depth and anomalies in the North Pacific ocean averaged during boreal winter (DJF) in meters (m). Statistically insignificant differences (p value > 0.05) are shaded in gray. Mean over DJF months is processed following paleo-seasonal adjustments with reference dates of 127 and 333 ka for LIG and MIS 9e, respectively. (a) LIG–PI, (b) MIS 9e–PI, (c) MIS 9e–LIG, and (d) PI.
Water age since last surface contact in years (g) and anomalies (a–f) for LIG–PI (a, c, e) and MIS 9e–PI (b, d, f). (a, b) 250 m depth, (c, d) 500 m depth, (e, f) 1000 m depth, and (g) 500 m PI reference. Statistically insignificant differences (p value > 0.01) are shaded in gray.
Annual-mean water age since last surface contact at 50 m depth overlaid with annual-mean surface wind at 10 m (e, f) and anomalies (a–d) off the west coast of North America (a, c, e) and off West Africa (b, d, f). (a, b) LIG–PI, (c, d) MIS 9e–PI, and (e, f) PI. Statistically insignificant differences of water age (p value > 0.01) are shaded in gray, and statistically insignificant wind difference vectors are not represented.
Saturated O2 concentration (left) and -1⋅AOU (apparent oxygen utilization, right) at 300 m depth. Statistically insignificant differences (p value > 0.01) are shaded in gray. Anomalies of LIG–PI (a, b), MIS 9e–PI (c, d), and MIS 9e–LIG (e, f). Full fields for PI are shown in panels (g) and (h). To facilitate the comparison between oxygen solubility and oxygen utilization, we multiplied AOU by -1.
Depth of minimum O2 concentration in hypoxic zones (< 62 mmol m−3). (a) PI, (b) LIG, and (c) MIS 9e.
Mixed layer depth in the Mediterranean Sea averaged over boreal winter (DJF) in meters (m). Mean over DJF months is processed following paleo-seasonal adjustments with reference dates of 127 and 333 ka for LIG and MIS 9e, respectively. (a) PI, (b) LIG, and (c) MIS 9e.
Data availability
The model data analyzed in this paper are published on the UNSW ResData repository at 10.26190/unsworks/30420.
Author contributions
BD performed the analyses and drafted the figures under the guidance of KM and LM. KM wrote the manuscript together with BD and LM. NY integrated the LIG and MIS 9e simulations. TZ and MC contributed to the model setup and integrated the PI simulation. BH commented on an advanced draft of the manuscript.
Competing interests
At least one of the (co-)authors is a member of the editorial board of Climate of the Past. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Acknowledgements
Bartholomé Duboc acknowledges funding from the ARC Centre of Excellence for Climate Extremes (CE170100023) for an undergraduate intern scholarship. All experiments were performed on the computational facility of the National Computational Infrastructure (NCI) owned by the Australian National University through awards under the Merit Allocation Scheme and the UNSW HPC at NCI Scheme. Katrin J. Meissner and Laurie Menviel acknowledge support from the Australian Research Council (DP180100048 and SR200100008). Nicholas K. H. Yeung acknowledges the Research Training Program provided by the Australian government, a top-up scholarship provided by the Climate Change Research Centre, and support from the ARC Centre of Excellence for Climate Extremes. Tilo Ziehn and Matthew Chamberlain received funding from the Australian Government under the National Environmental Science Program (NESP). This study was undertaken by PO2, a working group of the Past Global Changes (PAGES) project, which in turn received support from the Swiss Academy of Sciences and the Chinese Academy of Sciences.
Financial support
This research has been supported by the Australian Research Council (grant nos. CE170100023, DP180100048 and SR200100008) and the Australian Government under the National Environmental Science Program (NESP).
Review statement
This paper was edited by Zhongshi Zhang and reviewed by Christian Stepanek and two anonymous referees.
References Altieri, A., Harrison, S., Seemann, J., Collin, R., Diaz, R., and Knowlton, N.: Tropical dead zones and mass mortalities on coral reefs, P. Natl. Acad. Sci. USA, 114, 3660–3665, 2017. Andersen, M., Matthews, A., Vance, D., Bar-Matthews, M., Archer, C., and de Souza, G.: A 10-fold decline in the deep Eastern Mediterranean thermohaline overturning circulation during the last interglacial period, Earth Planet. Sc. Lett., 503, 58–67, 2018. Anderson, R. F., Sachs, J. P., Fleisher, M. Q., Allen, K. A., Yu, J., Koutavas, A., and Jaccard, S. L.: Deep-Sea Oxygen Depletion and Ocean Carbon Sequestration During the Last Ice Age, Global Biogeochem. Cy., 33, 301–317, 2019. Bahl, A., Gnanadesikan, A., and Pradal, M.-A.: Variations in Ocean Deoxygenation Across Earth System Models: Isolating the Role of Parameterized Lateral Mixing, Global Biogeochem. Cy., 33, 703–724, 2019.Barbosa Aguiar, A., Peliz, A., Neves, F., Bashmachnikov, I., and Carton, X.: Mediterranean outflow transports and entrainment estimates from observations and high-resolution modelling, Progr. Oceanogr., 131, 33–45, 10.1016/j.pocean.2014.11.008, 2015.Bartlein, P. J. and Shafer, S. L.: Paleo calendar-effect adjustments in time-slice and transient climate-model simulations (PaleoCalAdjust v1.0): impact and strategies for data analysis, Geosci. Model Dev., 12, 3889–3913, 10.5194/gmd-12-3889-2019, 2019.Bengtson, S. A., Menviel, L. C., Meissner, K. J., Missiaen, L., Peterson, C. D., Lisiecki, L. E., and Joos, F.: Lower oceanic δ13C during the last interglacial period compared to the Holocene, Clim. Past, 17, 507–528, 10.5194/cp-17-507-2021, 2021.Bereiter, B., Eggleston, S., Schmitt, J., Nehrbass-Ahles, C., Stocker, T. F., Fischer, H., Kipfstuhl, S., and Chappellaz, J.: Revision of the EPICA Dome C CO2 record from 800 to 600 kyr before present, Geophys. Res. Lett., 42, 542–549, 10.1002/2014GL061957, 2015.Berger, A.: Long-term variations of caloric insolation resulting from the earth's orbital elements, Quaternary Res., 9, 139–167, 10.1016/0033-5894(78)90064-9, 1978.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, 10.5194/bg-10-6225-2013, 2013.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, eaam7240, 10.1126/science.aam7240, 2018.Cartapanis, O., Tachikawa, K., Romero, O. E., and Bard, E.: Persistent millennial-scale link between Greenland climate and northern Pacific Oxygen Minimum Zone under interglacial conditions, Clim. Past, 10, 405–418, 10.5194/cp-10-405-2014, 2014.Chamberlain, M. A., Ziehn, T., and Law, R. M.: The Southern Ocean as the climate's freight train – driving ongoing global warming under zero-emission scenarios with ACCESS-ESM1.5, Biogeosciences, 21, 3053–3073, 10.5194/bg-21-3053-2024, 2024.Choudhury, D., Menviel, L., Meissner, K. J., Yeung, N. K. H., Chamberlain, M., and Ziehn, T.: Marine carbon cycle response to a warmer Southern Ocean: the case of the last interglacial, Clim. Past, 18, 507–523, 10.5194/cp-18-507-2022, 2022.Clarke, T. M., Wabnitz, C. C., Striegel, S., Frölicher, T. L., Reygondeau, G., and Cheung, W. W.: Aerobic growth index (AGI): An index to understand the impacts of ocean warming and deoxygenation on global marine fisheries resources, Progr. Oceanogr., 195, 102588, 10.1016/j.pocean.2021.102588, 2021.Craig, A., Valcke, S., and Coquart, L.: Development and performance of a new version of the OASIS coupler, OASIS3-MCT_3.0, Geosci. Model Dev., 10, 3297–3308, 10.5194/gmd-10-3297-2017, 2017. den Dulk, M., Reichart, G., Memon, G., Roelofs, E., Zachariasse, W., and van der Zwaan, G.: Benthic foraminiferal response to variations in surface water productivity and oxygenation in the northern Arabian Sea, Mar. Micropaleontol., 35, 43–66, 1998. Diaz, R. J. and Rosenberg, R.: Spreading Dead Zones and Consequences for Marine Ecosystems, Science, 321, 926–929, 2008.Donald, H. K., Swann, G. E. A., Rae, J. W. B., and Foster, G. L.: Carbon Cycle and Circulation Change in the North Pacific Ocean at the Initiation of Northern Hemisphere Glaciation Constrained by Boron-Based Proxies in Diatoms, Paleoceanography and Paleoclimatology, 39, e2024PA004968, 10.1029/2024PA004968, 2024.Duteil, O. and Oschlies, A.: Sensitivity of simulated extent and future evolution of marine suboxia to mixing intensity, Geophys. Res. Lett., 38, L06607, 10.1029/2011GL046877, 2011.Duteil, O., Koeve, W., Oschlies, A., Bianchi, D., Galbraith, E., Kriest, I., and Matear, R.: A novel estimate of ocean oxygen utilisation points to a reduced rate of respiration in the ocean interior, Biogeosciences, 10, 7723–7738, 10.5194/bg-10-7723-2013, 2013. Elderfield, H., Ferretti, P., Greaves, M., Crowhurst, S., McCave, I. N., Hodell, D., and Piotrowski, A. M.: Evolution of Ocean Temperature and Ice Volume Through the Mid-Pleistocene Climate Transition, Science, 337, 704–709, 2012.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, 10.5194/gmd-9-1937-2016, 2016.Fontugne, M. R. and Calvert, S. E.: Late Pleistocene Variability of the Carbon Isotopic Composition of Organic Matter in the Eastern Mediterranean: Monitor of Changes in Carbon Sources and Atmospheric CO2 Concentrations, Paleoceanography, 7, 1–20, 10.1029/91PA02674, 1992.Frölicher, T. L., Aschwanden, M. T., Gruber, N., Jaccard, S. L., Dunne, J. P., and Paynter, D.: Contrasting Upper and Deep Ocean Oxygen Response to Protracted Global Warming, Global Biogeochem. Cy., 34, e2020GB006601, 10.1029/2020GB006601, 2020.Galbraith, E. D., Kienast, M., Pedersen, T. F., and Calvert, S. E.: Glacial-interglacial modulation of the marine nitrogen cycle by high-latitude O2 supply to the global thermocline, Paleoceanography, 19, PA4007, 10.1029/2003PA001000, 2004. Ganeshram, R. S., Pedersen, T. F., Calvert, S. E., McNeill, G. W., and Fontugne, M. R.: Glacial-interglacial variability in denitrification in the World's Oceans: Causes and consequences, Paleoceanography, 15, 361–376, 2000.Garcia, H. E., Wang, Z., Bouchard, C., Cross, S. L., Paver, C. R., Reagan, J. R., Boyer, T. P., Locarnini, R. A., Mishonov, A. V., Baranova, O., Seidov, D., and Dukhovskoy D.: World Ocean Atlas 2023, Volume 3: Dissolved Oxygen, Apparent Oxygen Utilization, and Oxygen Saturation, edited by: Mishonov, A., NOAA Atlas NESDIS 91, National Centers for Environmental Information (U.S.), 10.25923/rb67-ns53, 2024.Glasscock, S. K., Hayes, C. T., Redmond, N., and Rohde, E.: Changes in Antarctic Bottom Water Formation During Interglacial Periods, Paleoceanography and Paleoclimatology, 35, e2020PA003867, 10.1029/2020PA003867, 2020.Glock, N., Erdem, Z., and Schönfeld, J.: The Peruvian oxygen minimum zone was similar in extent but weaker during the Last Glacial Maximum than Late Holocene, Communications Earth & Environment, 3, 307, 10.1038/s43247-022-00635-y, 2022.Gnanadesikan, A., Dunne, J. P., and John, J.: Understanding why the volume of suboxic waters does not increase over centuries of global warming in an Earth System Model, Biogeosciences, 9, 1159–1172, 10.5194/bg-9-1159-2012, 2012.Gottschalk, J., Skinner, L. C., Jaccard, S. L., Menviel, L., Nehrbass-Ahles, C., and Waelbroeck, C.: Southern Ocean link between changes in atmospheric CO2 levels and northern-hemisphere climate anomalies during the last two glacial periods, Quaternary Sci. Rev., 230, 106067, 10.1016/j.quascirev.2019.106067, 2020.Grant, K., Grimm, R., Mikolajewicz, U., Marino, G., Ziegler, M., and Rohling, E.: The timing of Mediterranean sapropel deposition relative to insolation, sea-level and African monsoon changes, Quaternary Sci. Rev., 140, 125–141, 10.1016/j.quascirev.2016.03.026, 2016. Griffies, S.: Elements of the Modular Ocean Model (MOM) 5: (2012 release with updates), Tech. rep., GFDL Ocean Group Tech. Rep. 7, NOAA, Geophysical Fluid Dynamics Laboratory, Princeton, NJ, 618, 2014.Hayes, C. T., Martínez-García, A., Hasenfratz, A. P., Jaccard, S. L., Hodell, D. A., Sigman, D. M., Haug, G. H., and Anderson, R. F.: A stagnation event in the deep South Atlantic during the last interglacial period, Science, 346, 1514–1517, 10.1126/science.1256620, 2014. Hoogakker, B., Elderfield, H., Schmiedl, G., McCave, I. N., and Rickaby, R. E. M.: Glacial–interglacial changes in bottom-water oxygen content on the Portuguese margin, Nat. Geosci., 8, 40–43, 2015. Hoogakker, B., Lu, Z., Umling, N., Jones, L., Zhou, X., Rickaby, R. E. M., Thunell, R., Cartapanis, O., and Galbraith, E.: Glacial expansion of oxygen-depleted seawater in the eastern tropical Pacific, Nature, 562, 410–413, 2018.Hoogakker, B. A. A., Thornalley, D. J. R., and Barker, S.: Millennial changes in North Atlantic oxygen concentrations, Biogeosciences, 13, 211–221, 10.5194/bg-13-211-2016, 2016.Hoogakker, B. A. A., Davis, C., Wang, Y., Kusch, S., Nilsson-Kerr, K., Hardisty, D. S., Jacobel, A., Reyes Macaya, D., Glock, N., Ni, S., Sepúlveda, J., Ren, A., Auderset, A., Hess, A. V., Meissner, K. J., Cardich, J., Anderson, R., Barras, C., Basak, C., Bradbury, H. J., Brinkmann, I., Castillo, A., Cook, M., Costa, K., Choquel, C., Diz, P., Donnenfield, J., Elling, F. J., Erdem, Z., Filipsson, H. L., Garrido, S., Gottschalk, J., Govindankutty Menon, A., Groeneveld, J., Hallmann, C., Hendy, I., Hennekam, R., Lu, W., Lynch-Stieglitz, J., Matos, L., Martínez-García, A., Molina, G., Muñoz, P., Moretti, S., Morford, J., Nuber, S., Radionovskaya, S., Raven, M. R., Somes, C. J., Studer, A. S., Tachikawa, K., Tapia, R., Tetard, M., Vollmer, T., Wang, X., Wu, S., Zhang, Y., Zheng, X.-Y., and Zhou, Y.: Reviews and syntheses: Review of proxies for low-oxygen paleoceanographic reconstructions, Biogeosciences, 22, 863–957, 10.5194/bg-22-863-2025, 2025. Hunke, E. C., Lipscomb, W. H., Turner, A. K., Jeffery, N., and Elliott, S.: CICE: the Los Alamos sea ice model documentation and software user's manual version 4.1 LA-CC-01-012), Tech. rep., T-3 Fluid Dynamics Group, Los Alamos National Laboratory, 675, 500, 2010.Ito, T., Minobe, S., Long, M. C., and Deutsch, C.: Upper ocean O2 trends: 1958–2015, Geophys. Res. Lett., 44, 4214–4223, 2017. Jaccard, S. L. and Galbraith, E. D.: Large climate-driven changes of oceanic oxygen concentrations during the last deglaciation, Nat. Geosci., 5, 151–156, 2012. Jaccard, S. L., Galbraith, E., Sigman, D., Haug, G., Francois, R., Pedersen, T., Dulski, P., and Thierstein, H.: Subarctic Pacific evidence for a glacial deepening of the oceanic respired carbon pool, Earth Planet. Sc. Lett., 277, 156–165, 2009.Jacobel, A., Anderson, R., Jaccard, S., McManus, J., Pavia, F., and Winckler, G.: Deep Pacific storage of respired carbon during the last ice age: Perspectives from bottom water oxygen reconstructions, Quaternary Sci. Rev., 230, 106065, 10.1016/j.quascirev.2019.106065, 2020.Jimenez-Espejo, F. J., García-Alix, A., Harada, N., Bahr, A., Sakai, S., Iijima, K., Chang, Q., Sato, K., Suzuki, K., and Ohkouchi, N.: Changes in detrital input, ventilation and productivity in the central Okhotsk Sea during the marine isotope stage 5e, penultimate interglacial period, J. Asian Earth Sci., 156, 189–200, 10.1016/j.jseaes.2018.01.032, 2018.Kageyama, M., Sime, L. C., Sicard, M., Guarino, M.-V., de Vernal, A., Stein, R., Schroeder, D., Malmierca-Vallet, I., Abe-Ouchi, A., Bitz, C., Braconnot, P., Brady, E. C., Cao, J., Chamberlain, M. A., Feltham, D., Guo, C., LeGrande, A. N., Lohmann, G., Meissner, K. J., Menviel, L., Morozova, P., Nisancioglu, K. H., Otto-Bliesner, B. L., O'ishi, R., Ramos Buarque, S., Salas y Melia, D., Sherriff-Tadano, S., Stroeve, J., Shi, X., Sun, B., Tomas, R. A., Volodin, E., Yeung, N. K. H., Zhang, Q., Zhang, Z., Zheng, W., and Ziehn, T.: A multi-model CMIP6-PMIP4 study of Arctic sea ice at 127 ka: sea ice data compilation and model differences, Clim. Past, 17, 37–62, 10.5194/cp-17-37-2021, 2021.Kallel, N., Duplessy, J.-C., Labeyrie, L., Fontugne, M., Paterne, M., and Montacer, M.: Mediterranean pluvial periods and sapropel formation over the last 200 000 years, Palaeogeogr. Palaeocl., 157, 45–58, 10.1016/S0031-0182(99)00149-2, 2000.Köng, E., Zaragosi, S., Schneider, J.-L., Garlan, T., Bachèlery, P., Sabine, M., and San Pedro, L.: Gravity-Driven Deposits in an Active Margin (Ionian Sea) Over the Last 330,000 Years, Geochem. Geophy. Geosy., 18, 4186–4210, 10.1002/2017GC006950, 2017.Konijnendijk, T., Ziegler, M., and Lourens, L.: Chronological constraints on Pleistocene sapropel depositions from high resolution geochemical records of ODP Sites 967 and 968, Newsl. Stratigr., 47, 263–282, 10.1127/0078-0421/2014/0047, 2014. Kowalczyk, E., Stevens, L., Law, R., Dix, M., Wang, Y., Harman, I., Haynes, K., Srbinovsky, J., Pak, B., and Ziehn, T.: The land surface model component of ACCESS: description and impact on the simulated surface climatology, Aust. Meteorol. Ocean., 63, 65–82, 2013.Kwiatkowski, L., Torres, O., Bopp, L., Aumont, O., Chamberlain, M., Christian, J. R., Dunne, J. P., Gehlen, M., Ilyina, T., John, J. G., Lenton, A., Li, H., Lovenduski, N. S., Orr, J. C., Palmieri, J., Santana-Falcón, Y., Schwinger, J., Séférian, R., Stock, C. A., Tagliabue, A., Takano, Y., Tjiputra, J., Toyama, K., Tsujino, H., Watanabe, M., Yamamoto, A., Yool, A., and Ziehn, T.: Twenty-first century ocean warming, acidification, deoxygenation, and upper-ocean nutrient and primary production decline from CMIP6 model projections, Biogeosciences, 17, 3439–3470, 10.5194/bg-17-3439-2020, 2020. Landry, C., Steele, S., Manning, S., and Cheek, A.: Long term hypoxia suppresses reproductive capacity in the estuarine fish, Fundulus grandis, Comp. Biochem. Physiol. A Mol. Integr. Physiol., 148, 317–23, 2007.Law, R. M., Ziehn, T., Matear, R. J., Lenton, A., Chamberlain, M. A., Stevens, L. E., Wang, Y.-P., Srbinovsky, J., Bi, D., Yan, H., and Vohralik, P. F.: The carbon cycle in the Australian Community Climate and Earth System Simulator (ACCESS-ESM1) – Part 1: Model description and pre-industrial simulation, Geosci. Model Dev., 10, 2567–2590, 10.5194/gmd-10-2567-2017, 2017.Lessa, D. V., Santos, T. P., Venancio, I. M., and Albuquerque, A. L. S.: Offshore expansion of the Brazilian coastal upwelling zones during Marine Isotope Stage 5, Global Planet. Change, 158, 13–20, 10.1016/j.gloplacha.2017.09.006, 2017.Levermann, A. and Born, A.: Bistability of the Atlantic subpolar gyre in a coarse-resolution climate model, Geophys. Res. Lett., 34, L24605, 10.1029/2007GL031732, 2007.Levin, L. A.: Manifestation, Drivers, and Emergence of Open Ocean Deoxygenation, Annu. Rev. Mar. Sci., 10, 229–260, 10.1146/annurev-marine-121916-063359, 2018. Long, M. C., Deutsch, C., and Ito, T.: Finding forced trends in oceanic oxygen, Global Biogeochem. Cy., 30, 381–397, 2016.Loulergue, L., Schilt, A., Spahni, R., Masson-Delmotte, V., Blunier, T., Lemieux, B., Barnola, J.-M., Raynaud, D., Stocker, T. F., and Chappellaz, J.: Orbital and millennial-scale features of atmospheric CH4 over the past 800 000 years, Nature, 453, 383–386, 10.1038/nature06950, 2008.Lu, Z., Hoogakker, B., Hillenbrand, C.-D., Zhou, X., Thomas, E., Gutchess, K., Lu, W., Jones, L., and Rickaby, R.: Oxygen depletion recorded in upper waters of the glacial Southern Ocean, Nat. Commun., 7, 11146, 10.1038/ncomms11146, 2016. Mackallah, C., Chamberlain, M. A., Law, R. M., Dix, M., Ziehn, T., Bi, D., Bodman, R., Brown, J. R., Dobrohotoff, P., Druken, K., Evans, B., Harman, I. N., Hayashida, H., Holmes, R., Kiss, A. E., Lenton, A., Liu, Y., Marsland, S., Meissner, K., Menviel, L., O'Farrell, S., Rashid, H. A., Ridzwan, S., Savita, A., Srbinovsky, J., Sullivan, A., Trenham, C., Vohralik, P. F., Wang, Y. P., Williams, G., Woodhouse, M. T., and Yeung, N.: ACCESS datasets for CMIP6: methodology and idealised experiments, Journal of Southern Hemisphere Earth Systems Science, 72, 93–116, 2022.Maréchal, C., Boutier, A., Mélières, M.-A., Clauzel, T., Betancort, J. F., Lomoschitz, A., Meco, J., Fourel, F., Barral, A., Amiot, R., and Lécuyer, C.: Last Interglacial sea surface warming during the sea-level highstand in the Canary Islands: Implications for the Canary Current and the upwelling off African coast, Quaternary Sci. Rev., 234, 106246, 10.1016/j.quascirev.2020.106246, 2020. Martin, G. M., Milton, S. F., Senior, C. A., Brooks, M. E., Ineson, S., Reichler, T., and Kim, J.: Analysis and Reduction of Systematic Errors through a Seamless Approach to Modeling Weather and Climate, J. Climate, 23, 5933–5957, 2010.Matul, A., Abelmann, A., Khusid, T., Chekhovskaya, M., Kaiser, A., Nürnberg, D., and Tiedemann, R.: Late Quaternary changes of the oxygen conditions in the bottom and intermediate waters on the western Kamchatka continental slope, the Sea of Okhotsk, Deep-Sea Res. Pt. II, 125–126, 184–190, 10.1016/j.dsr2.2013.03.023, 2016.McCormick, L. and Levin, L.: Physiological and ecological implications of ocean deoxygenation for vision in marine organisms, Philos. Trans. A Math. Phys. Eng. Sci., 375, 20160322, 10.1098/rsta.2016.0322, 2017.Meissner, K. J., Galbraith, E. D., and Völker, C.: Denitrification under glacial and interglacial conditions: A physical approach, Paleoceanography, 20, PA3001, 10.1029/2004PA001083, 2005.Melki, T., Kallel, N., and Fontugne, M.: The nature of transitions from dry to wet condition during sapropel events in the Eastern Mediterranean Sea, Palaeogeogr. Palaeocl., 291, 267–285, 10.1016/j.palaeo.2010.02.039, 2010.Menviel, L., Govin, A., Avenas, A., Meissner, K., Grant, K., and Tzedakis, C.: Drivers of the evolution and amplitude of African Humid Periods, Communications Earth & Environment, 2, 237, 10.1038/s43247-021-00309-1, 2021.Mitsui, T., Tzedakis, P. C., and Wolff, E. W.: Insolation evolution and ice volume legacies determine interglacial and glacial intensity, Clim. Past, 18, 1983–1996, 10.5194/cp-18-1983-2022, 2022.Moffitt, S., Moffitt, R., Sauthoff, W., Davis, C., Hewett, K., and Hill, T. M.: Paleoceanographic Insights on Recent Oxygen Minimum Zone Expansion: Lessons for Modern Oceanography, PLoS ONE, 10, e0115246, 10.1371/journal.pone.0115246, 2015.Muhs, D. R., Meco, J., and Simmons, K. R.: Uranium-series ages of corals, sea level history, and palaeozoogeography, Canary Islands, Spain: An exploratory study for two Quaternary interglacial periods, Palaeogeogr. Palaeocl., 394, 99–118, 10.1016/j.palaeo.2013.11.015, 2014.Murat, A.: Pliocene-Pleistocene occurrence of sapropels in the Western Mediterranean Sea and their relation to eastern Mediterranean sapropels, in: Proceedings of the Ocean Drilling Program, 161 Scientific Results, Texas A&M University, Ocean Drilling Program, College Station, TX, United States, 10.2973/odp.proc.sr.161.244.1999, 1999. Muratli, J. M., Chase, Z., Mix, A. C., and McManus, J.: Increased glacial-age ventilation of the Chilean margin by Antarctic Intermediate Water, Nat. Geosci., 3, 23–26, 2010.Nameroff, T. J., Calvert, S. E., and Murray, J. W.: Glacial-interglacial variability in the eastern tropical North Pacific oxygen minimum zone recorded by redox-sensitive trace metals, Paleoceanography, 19, PA1010, 10.1029/2003PA000912, 2004.Oke, P. R., Griffin, D. A., Schiller, A., Matear, R. J., Fiedler, R., Mansbridge, J., Lenton, A., Cahill, M., Chamberlain, M. A., and Ridgway, K.: Evaluation of a near-global eddy-resolving ocean model, Geosci. Model Dev., 6, 591–615, 10.5194/gmd-6-591-2013, 2013.Oki, T., Nishimura, T., and Dirmeyer, P.: Assessment of Annual Runoff from Land Surface Models Using Total Runoff Integrating Pathways (TRIP), J. Meteorol. Soc. Jpn. Ser. II, 77, 235–255, 10.2151/jmsj1965.77.1B_235, 1999.Oschlies, A.: A committed fourfold increase in ocean oxygen loss, Nat. Commun., 12, 2307, 10.1038/s41467-021-22584-4, 2021.Oschlies, A., Duteil, O., Getzlaff, J., Koeve, W., Landolfi, A., and Schmidtko, S.: Patterns of deoxygenation: sensitivity to natural and anthropogenic drivers, Philos. Trans. A Math. Phys. Eng. Sci., 375, 20160325, 10.1098/rsta.2016.0325, 2017.Oschlies, A., Brandt, P., Stramma, L., and Schmidtko, S.: Drivers and mechanisms of ocean deoxygenation, Nat. Geosci., 11, 467–473, 10.1038/s41561-018-0152-2, 2018.Otto-Bliesner, B. L., Braconnot, P., Harrison, S. P., Lunt, D. J., Abe-Ouchi, A., Albani, S., Bartlein, P. J., Capron, E., Carlson, A. E., Dutton, A., Fischer, H., Goelzer, H., Govin, A., Haywood, A., Joos, F., LeGrande, A. N., Lipscomb, W. H., Lohmann, G., Mahowald, N., Nehrbass-Ahles, C., Pausata, F. S. R., Peterschmitt, J.-Y., Phipps, S. J., Renssen, H., and Zhang, Q.: The PMIP4 contribution to CMIP6 – Part 2: Two interglacials, scientific objective and experimental design for Holocene and Last Interglacial simulations, Geosci. Model Dev., 10, 3979–4003, 10.5194/gmd-10-3979-2017, 2017.Otto-Bliesner, B. L., Brady, E. C., Zhao, A., Brierley, C. M., Axford, Y., Capron, E., Govin, A., Hoffman, J. S., Isaacs, E., Kageyama, M., Scussolini, P., Tzedakis, P. C., Williams, C. J. R., Wolff, E., Abe-Ouchi, A., Braconnot, P., Ramos Buarque, S., Cao, J., de Vernal, A., Guarino, M. V., Guo, C., LeGrande, A. N., Lohmann, G., Meissner, K. J., Menviel, L., Morozova, P. A., Nisancioglu, K. H., O'ishi, R., Salas y Mélia, D., Shi, X., Sicard, M., Sime, L., Stepanek, C., Tomas, R., Volodin, E., Yeung, N. K. H., Zhang, Q., Zhang, Z., and Zheng, W.: Large-scale features of Last Interglacial climate: results from evaluating the lig127k simulations for the Coupled Model Intercomparison Project (CMIP6)–Paleoclimate Modeling Intercomparison Project (PMIP4), Clim. Past, 17, 63–94, 10.5194/cp-17-63-2021, 2021. Past Interglacials Working Group of PAGES: Interglacials of the last 800 000 years, Rev. Geophys., 54, 162–219, 2016. Penn, J. L. and Deutsch, C.: Avoiding ocean mass extinction from climate warming, Science, 376, 524–526, 2022.Pichevin, L., Martinez, P., Bertrand, P., Schneider, R., Giraudeau, J., and Emeis, K.: Nitrogen cycling on the Namibian shelf and slope over the last two climatic cycles: Local and global forcings, Paleoceanography, 20, PA2006, 10.1029/2004PA001001, 2005.Poore, R., Dowsett, H., Barron, J., Heusser, L., Ravelo, A., and Mix, A.: Multiproxy Record of the Last Interglacial (MIS 5e) off Central and Northern California, U.S.A., from Ocean Drilling Program Sites 1018 and 1020, U.S. Geological Survey, Professional Paper 1632, 10.3133/pp1632, 2000. Rohling, E., Grant, K., Bolshaw, M., Roberts, A., Siddall, M., Hemleben, C., and Kucera, M.: Antarctic temperature and global sea level closely coupled over the past five glacial cycles, Nat. Geosci., 2, 500–504, 2009. Rohling, E., Marino, G., and Grant, K.: Mediterranean climate and oceanography, and the periodic development of anoxic events (sapropels), Earth-Sci. Rev., 143, 62–97, 2015. Rose, K., Gutiérrez, D., Breitburg, D., Conley, D., Craig, K., Froehlich, H., Jeyabaskaran, R., Kripa, V., Mbaye, B., Mohamed, K., Padua, S., and Prema, D.: Chapter 10: Impacts of ocean deoxygenation on fisheries, in: Ocean deoxygenation, edited by: Laffoley, D. and Baxter, J. M., International Union for Conservation of Nature and Natural Resources, ISBN 978-2-8317-2013-5, 2019.Rush, D., Talbot, H. M., van der Meer, M. T. J., Hopmans, E. C., Douglas, B., and Sinninghe Damsté, J. S.: Biomarker evidence for the occurrence of anaerobic ammonium oxidation in the eastern Mediterranean Sea during Quaternary and Pliocene sapropel formation, Biogeosciences, 16, 2467–2479, 10.5194/bg-16-2467-2019, 2019. Sachs, J. P. and Repeta, D. J.: Oligotrophy and Nitrogen Fixation During Eastern Mediterranean Sapropel Events, Science, 286, 2485–2488, 1999. Sampaio, E., Santos, C., Rosa, I. C., Ferreira, V., Pörtner, H. O., Duarte, C. M., Levin, L. A., and Rosa, R. P.: Impacts of hypoxic events surpass those of future ocean warming and acidification, Nature Ecology and Evolution, 5, 311–321, 2021.Santana-Falcón, Y., Yamamoto, A., Lenton, A., Jones, C. D., Burger, F. A., John, J. G., Tjiputra, J., Schwinger, J., Kawamiya, M., Frölicher, T. L., Ziehn, T., and Séférian, R.: Irreversible loss in marine ecosystem habitability after a temperature overshoot, Communications Earth & Environment, 4, 343, 10.1038/s43247-023-01002-1, 2023.Schilt, A., Baumgartner, M., Blunier, T., Schwander, J., Spahni, R., Fischer, H., and Stocker, T. F.: Glacial–interglacial and millennial-scale variations in the atmospheric nitrous oxide concentration during the last 800 000 years, Quaternary Sci. Rev., 29, 182–192, 10.1016/j.quascirev.2009.03.011, 2010. Schmidtko, S., Stramma, L., and Visbeck, M.: Decline in global oceanic oxygen content during the past five decades, Nature, 542, 335–339, 2017. Scholz, F., McManus, J., Mix, A. C., Hensen, C., and Schneider, R. R.: The impact of ocean deoxygenation on iron release from continental margin sediments, Nat. Geosci., 7, 433–437, 2014.Stott, L. D., Neumann, M., and Hammond, D.: Intermediate water ventilation on the northeastern Pacific Margin during the Late Pleistocene inferred from benthic foraminiferal δ13C, Paleoceanography, 15, 161–169, 10.1029/1999PA000375, 2000.Sweere, T., Hennekam, R., Vance, D., and Reichart, G.-J.: Molybdenum isotope constraints on the temporal development of sulfidic conditions during Mediterranean sapropel intervals, Geochemical Perspectives Letters, 17, 125–141, 10.7185/geochemlet.2108, 2021.Tesi, T., Asioli, A., Massara, E. P., Montagna, P., Pellegrini, C., Nogarotto, A., Cipriani, A., Piva, A., Muschitiello, F., Rovere, M., Viaggi, P., and Trincardi, F.: Abrupt and persistent shutdown of the thermohaline forcing during MIS5e in the Adriatic Sea: Insights from shallow-water sapropel sediments, Quaternary Science Advances, 13, 100134, 10.1016/j.qsa.2023.100134, 2024.The HadGEM2 Development Team: G. M. Martin, Bellouin, N., Collins, W. J., Culverwell, I. D., Halloran, P. R., Hardiman, S. C., Hinton, T. J., Jones, C. D., McDonald, R. E., McLaren, A. J., O'Connor, F. M., Roberts, M. J., Rodriguez, J. M., Woodward, S., Best, M. J., Brooks, M. E., Brown, A. R., Butchart, N., Dearden, C., Derbyshire, S. H., Dharssi, I., Doutriaux-Boucher, M., Edwards, J. M., Falloon, P. D., Gedney, N., Gray, L. J., Hewitt, H. T., Hobson, M., Huddleston, M. R., Hughes, J., Ineson, S., Ingram, W. J., James, P. M., Johns, T. C., Johnson, C. E., Jones, A., Jones, C. P., Joshi, M. M., Keen, A. B., Liddicoat, S., Lock, A. P., Maidens, A. V., Manners, J. C., Milton, S. F., Rae, J. G. L., Ridley, J. K., Sellar, A., Senior, C. A., Totterdell, I. J., Verhoef, A., Vidale, P. L., and Wiltshire, A.: The HadGEM2 family of Met Office Unified Model climate configurations, Geosci. Model Dev., 4, 723–757, 10.5194/gmd-4-723-2011, 2011. Ufkes, E. and Kroon, D.: Sensitivity of south-east Atlantic planktonic foraminifera to mid-Pleistocene climate change, Palaeontology, 55, 183–204, 10.1111/j.1475-4983.2011.01119.x, 2012.Vaquer-Sunyer, R. and Duarte, C. M.: Thresholds of hypoxia for marine biodiversity, P. Natl. Acad. Sci. USA, 105, 15452–15457, 10.1073/pnas.0803833105, 2008.Yao, Z., Shi, X., Yin, Q., Jaccard, S., Liu, Y., Guo, Z., Gorbarenko, S. A., Wang, K., Chen, T., Wu, Z., Nan, Q., Zou, J., Wang, H., Cui, J., Wang, A., Yang, G., Zhu, A., Bosin, A., Vasilenko, Y., and Yu, Y.: Ice sheet and precession controlled subarctic Pacific productivity and upwelling over the last 550 000 years, Nat. Commun., 15, 3489, 10.1038/s41467-024-47871-8, 2024.Yeung, N.: Simulated ocean oxygenation during the interglacials MIS 5e and MIS 9e, UNSW ResData [data set], 10.26190/unsworks/30420, 2024.Yeung, N. K.-H., Menviel, L., Meissner, K. J., Taschetto, A. S., Ziehn, T., and Chamberlain, M.: Land–sea temperature contrasts at the Last Interglacial and their impact on the hydrological cycle, Clim. Past, 17, 869–885, 10.5194/cp-17-869-2021, 2021.Yeung, N. K.-H., Menviel, L., Meissner, K. J., Choudhury, D., Ziehn, T., and Chamberlain, M. A.: Last Interglacial subsurface warming on the Antarctic shelf triggered by reduced deep-ocean convection, Communications Earth & Environment, 5, 212, 10.1038/s43247-024-01383-x, 2024. Ziehn, T., Chamberlain, M. A., Law, R. M., Lenton, A., Bodman, R. W., 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, 70, 193–214, 2020.