the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
![](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-avatar-thumb150.png)
Climatic controls on metabolic constraints in the ocean
Precious Mongwe
Matthew Long
Takamitsu Ito
Curtis Deutsch
Yeray Santana-Falcón
Observations and models indicate that climate warming is associated with the loss of dissolved oxygen from the ocean. Dissolved oxygen is a fundamental requirement for heterotrophic marine organisms (except marine mammals) and, since the basal metabolism of ectotherms increases with temperature, warming increases organisms' oxygen demands. Therefore, warming and deoxygenation pose a compound threat to marine ecosystems. In this study, we leverage an ecophysiological framework and a compilation of empirical trait data quantifying the temperature sensitivity and oxygen requirements of metabolic rates for a range of marine species (“ecotypes”). Using the Community Earth System Model Large Ensemble, we investigate how natural climate variability and anthropogenic forcing impact the ability of marine environments to support aerobic metabolisms on interannual to multi-decadal timescales. Warming and deoxygenation projected over the next several decades will yield a reduction in the volume of viable ocean habitats. We find that fluctuations in temperature and oxygen associated with natural variability are distinct from those associated with anthropogenic forcing in the upper ocean. Further, the joint temperature–oxygen anthropogenic signal emerges sooner than temperature and oxygen independently from natural variability. Our results demonstrate that anthropogenic perturbations underway in the ocean will strongly exceed those associated with the natural system; in many regions, organisms will be pushed closer to or beyond their physiological limits, leaving the ecosystem more vulnerable to extreme temperature–oxygen events.
- Article
(13327 KB) - Full-text XML
-
Supplement
(477 KB) - BibTeX
- EndNote
Dissolved oxygen (O2) is a fundamental metabolic requirement for heterotrophic marine organisms, excluding marine mammals (Pörtner, 2002; Keeling et al., 2010; Tiano et al., 2014). The decline in ocean O2 due to warming is a tendency that has been long predicted by models (Keeling et al., 2010; Long et al., 2016; Oschlies et al., 2018) and that was recently found to be evident at the global scale in compilations of in situ observations (Schmidtko et al., 2017; Ito et al., 2017). Deoxygenation is driven by the direct effect of reduced oxygen solubility with warming, compounded by buoyancy-induced stratification in the upper ocean, which weakens the ventilation-mediated supply of fresh oxygen to the ocean interior. While the full ecological implications of ocean deoxygenation remain uncertain, it is clear that the physiological impacts of oxygen loss on marine organisms can be considered explicitly in the context of warming: basal metabolic rates for ectothermic organisms depend on ambient temperature and increase with warming (Gillooly et al., 2001); thus, higher temperatures impose additional demands for oxygen to sustain aerobic respiration (Deutsch et al., 2015). Consequently, as the ocean warms, even present-day oxygen distributions may be insufficient to meet the oxygen demands of organisms living near key physiological thresholds (Deutsch et al., 2022).
Model projections clearly demonstrate that warming and deoxygenation are consequences of human-driven climate change, yet natural climate variability also produces important fluctuations in these quantities. Indeed, evidence suggests that natural variability contributes to hypoxic events, such as those observed in the California Current, where fish and benthic-organism mortality has been associated with low-O2 waters impinging on the continental shelf (Pozo Buil and Di Lorenzo, 2017; Howard et al., 2020). A clear understanding of how natural climate variability drives fluctuations in metabolic state and the associated implications for organisms is a critical context in which to view long-term climate warming. Given that the natural system is highly dynamic, climate change signals are often masked by decadal-scale variability (Ito and Deutsch, 2010). While numerous authors have considered the detection and attribution of climate change for physical and biogeochemical variables (Rodgers et al., 2015; Long et al., 2016; Schlunegger et al., 2019), comparatively little attention has been devoted to explicitly characterizing the relative influence of natural and anthropogenic drivers of changes in the ocean's capacity to support aerobic life. In this study, we approach this challenge by leveraging the concept of the metabolic index (Φ) introduced by Deutsch et al. (2015). Φ is based on the notion that aerobic organisms can persist only where the ambient oxygen partial pressure (pO2) is sufficient to sustain respiration. Φ incorporates an explicit representation of the dependence of metabolic oxygen demand on temperature, thus providing a framework to consider how joint oxygen and temperature variability constrain viable habitats in the ocean.
Many ocean organisms may already be under threat from deoxygenation (Hoegh-Guldberg and Bruno, 2010; Breitburg et al., 2018); however, ongoing climate-driven loss of oxygen raises important questions about the future of marine ecosystems: how will anthropogenic changes in dissolved oxygen and temperature affect the capacity of ocean habitats to support aerobic metabolism? What are the spatial and temporal distributions of changes in the ocean's metabolic state associated with climate variability? At what point can anthropogenic change in the ocean's metabolic state be distinguished from natural variability? This study addresses these questions using a combination of metabolic theory; a dataset quantifying key physiological parameters for a collection of marine species adapted to specific environments (“ecotypes”); and the oxygen and temperature distributions as simulated in the Community Earth System Model, version 1, Large Ensemble (CESM1-LE), which includes 34 members simulating ocean biogeochemistry under climate variability and change for the period 1920–2100, forced using historical data and the Representative Concentration Pathway Scenario 8.5 (RCP85) (Kay et al., 2015; Long et al., 2016).
This paper is organized as follows. Section 2 presents a brief overview of the relevant metabolic theory and the associated empirical datasets and describes our approach to the analysis. In Sect. 3, we present results quantifying the joint temperature–oxygen variability simulated in the CESM1-LE, evaluating the spatiotemporal structure of variability in marine ecotype habitats, including long-term trends based on the RCP8.5 scenario and time of emergence (ToE). The main outcomes of the results are synthesized in Sect. 4 and summarized in Sect. 5.
2.1 Metabolic index
Empirical studies measuring thermal tolerance and oxygen requirements in the laboratory for an array of marine organisms have enabled an assessment of lethal thresholds (Vaquer-Sunyer and Duarte, 2008; Rosewarne et al., 2016). These data, coupled with recent advances in the theoretical frameworks, enable both explanatory and predictive power in the context of a dynamic environment (Deutsch et al., 2015; Penn et al., 2018; Howard et al., 2020). The fundamental insights here are that basal metabolic rates for ectothermic marine organisms depend on ambient temperature and generally increase with warming (Gillooly et al., 2001). Increasing basal metabolic rates impose additional demands for oxygen. Organisms use oxygen dissolved in seawater, and acquisition tends to be limited by diffusive processes; thus, oxygen supply is related to the ambient pO2. The ratio of oxygen supply to temperature-dependent demand provides a critical indicator of the capacity for an organism to meet its metabolic requirements. Deutsch et al. (2015) formalized these concepts into a quantity termed the metabolic index (Φ), which is defined as the ratio of oxygen supply to an organism's resting metabolic demand. Oxygen supply is parameterized according to a biomass-dependent scaling of pO2, capturing variations in the efficiency with which organisms acquire and utilize O2. This can be expressed as , where represents the gas transfer between an organism and its environment, and Bδ is the scaling of supply with biomass B (Piiper et al., 1971). Gas supply is represented as an Arrhenius function:
The resting metabolic demand is also expressed using the Arrhenius equation as follows:
where αD is a species-specific basal metabolic rate, Ed (eV) is the temperature dependence of oxygen supply, T is temperature, Tref is the reference temperature (15 °C), and kB is the Boltzmann constant (Gillooly et al., 2001). Gas transfer is kinematically slow at low temperatures, and, hence, organism viability can be limited by the energy required to acquire oxygen at low temperatures; thus, Eo varies with temperature. Here, we account for this by adding the temperature dependence () to Eo in the equations above (), using the mean value of eV, consistently with Deutsch et al. (2020). The metabolic index can thus be written as the ratio of :
where (atm−1) is the hypoxic tolerance, and (Es is the temperature dependence of oxygen supply) (Deutsch et al., 2015; Penn et al., 2018). The exponent, , is the allometric scaling of the supply-to-demand ratio with biomass, and it is typically near zero. Therefore, in the analysis that follows, we presume unit biomass and thus neglect the potential impacts of variations in biomass.
If Φ falls below a critical threshold value of 1, conditions are physiologically unsustainable: an organism cannot meet its basic resting metabolic oxygen requirements. Conversely, values of Φ above 1 enable organismal metabolic rates to increase by a factor of Φ above resting levels, permitting critical activities such as feeding, defense, growth and reproduction. Thus, for a given environment and species, Φ provides an estimate of the ratio of the maximum sustainable metabolic rate to the minimum rate necessary for basal metabolism. Deutsch et al. (2015) inferred the ratio of active to resting energetic demand by examining the biogeographic distribution of several species, finding that range boundaries coincide with values of Φ=1.5–7. This threshold, termed the critical rate (Φcrit), represents the minimum metabolic index required for an organism to sustain an active metabolic state, which is a more meaningful ecological threshold than requirements for resting metabolism. Therefore, in this study, we define a quantity Φ′, derived by dividing Φ by Φcrit so that, when Φ falls below 1, the organism can no longer sustain its active metabolic demand and will need to make physiological trade-offs. Accounting for these active metabolic requirements, we use an adjusted definition of the hypoxic tolerance trait, , where Ac is termed the “ecological hypoxia tolerance”, consistently with Howard et al. (2020). Where (i.e., Φ>Φcrit), an organism can sustain an active metabolic rate; where (i.e., Φ<Φcrit), O2 is insufficient, and an active metabolic state is not viable. Henceforth, our analysis focuses on Φ′; in the subsequent text and figures, metabatic index refers to Φ′ () throughout.
2.2 Physiological dataset
We make use of a dataset describing physiological parameters for a collection of 61 marine ecotypes spanning a range of ecological hypoxic tolerances (Ac) and temperature sensitivities (Eo) (Penn et al., 2018; Deutsch et al., 2020, Fig. 1a). The 61 species span benthic and pelagic habitats across four phyla in all ocean basins (Arthropoda, Chordata, Mollusca and Cnidaria). The dataset includes 28 malacostracans, 21 fishes; three bivalves and cephalopods; two copepods; and one of each for gastropods, ascidians, scleractinian corals and sharks with body mass spans of 8 orders of magnitude (Penn et al., 2018). We illustrate how the physiological traits Eo and Ac constrain habitat viability in the context of distributions of pO2 and temperature in the marine environment in Fig. 1b, which shows the minimum pO2 (i.e., pO2 at Φcrit) required to sustain an active metabolic state as a function of temperature for five combinations of Eo and Ac. The five combinations are derived from sampling the probability distributions of Eo and Ac (Fig. 1a) at the 10th, 50th and 90th percentile values (illustrated by colored stars in Fig. 1a and corresponding curves in Fig. 1b). We assume that the trait distributions are independent, which is a reasonably modest simplification; Eo is represented by a normal distribution, and Ac is represented by a lognormal distribution function (Fig. S1 in the Supplement). The pO2 at Φcrit curves shown in Fig. 1b delineate regions of pO2–temperature space that are habitable (above the curve) and uninhabitable (below the curve). The reversing curvature of pO2 at Φcrit in Fig. 1b at low temperatures captures the decrease in the organisms' oxygen acquisition efficiencies in cooler conditions, yielding cold intolerance. At very low temperatures, gas transfer is limited by the decrease in molecular gas diffusion; as a consequence, oxygen transfer into the organisms requires energy, yielding cold intolerance. This is well illustrated by the blue line in Fig. 1b.
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f01](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f01-thumb.png)
Figure 1Physiological traits determining hypoxic tolerance. (a) Scatter plot of 61 marine ecotypes for which empirically derived estimates of activation energy (Eo) and the ecological hypoxic tolerance (Ac) have been determined (Penn et al., 2018). The color shows the density of occurrence for the 61 marine ecotypes in the Ac–Eo trait space. (b) The minimum pO2 required to sustain an active metabolic state (i.e., pO2 at Φcrit; Deutsch et al., 2020) for five combinations of Ac and Eo, corresponding to the stars in panel (a); these are combinations of the 10th, 50th and 90th percentile values for each parameter. Below the pO2 lines shown, the organism would experience an oxygen deficit relative to its active metabolism requirements, effectively signifying the species-specific hypoxic conditions, based on physiological traits, for this range of temperatures.
To illustrate how the trait combinations of Eo and Ac exert control on the geographic distribution of organisms in the marine environment (Deutsch et al., 2020), we use observations of pO2 and T along a zonal transect of the Pacific Ocean and plot Φ′ for nine combinations of Eo and Ac percentile values (Fig. 2). The color bars in Fig. 2a–i show the metabolic index for an active state (Φ′); regions with values above 1 are habitable (colored), while regions with values below 1 are uninhabitable (white) on the basis of metabolic constraints (other ecological considerations are not considered). The subplots in the upper portion of the figure are arranged according to the same trait axes shown in Fig. 1a; Eo increases horizontally from left to right, and Ac increases from the bottom to the top. For the trait combination in the bottom left (low Eo, low Ac; Fig. 2g), metabolism is relatively insensitive to temperature, and tolerance for low pO2 is poor. Thus, ecotypes with low Eo and low Ac are restricted to high-latitude surface waters, where temperatures are cool and where pO2 is abundant (Fig. 2g). As Eo increases from left to right, metabolic rates become more sensitive to temperature. Then, habitat is gained at depth, where temperatures are cooler and where higher temperature sensitivity confers an advantage (Fig. 2g–i). From the bottom to the top, the increase in the tolerance of low-pO2 conditions increases habitability in regions of low pO2, enabling organisms to expand beyond high-latitude surface waters (Fig. 2g–a). The biogeographic range for organisms with high Ac is modulated by Eo; as temperature sensitivity increases, ecotype viability at high latitudes increases, but tropical surface waters become less viable (Fig. 2a–c). Henceforth, our analysis will utilize the metabolic index of the median ecotype (Eo=0.34, Ac=7.4; Fig. 2e) for illustrative purposes; i.e., all metabolic index figures refer to this median ecotype unless otherwise stated.
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f02](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f02-thumb.png)
Figure 2Annual mean metabolic index (Φ′) for nine combinations of the ecological traits Eo (metabolic temperature sensitivity) and Ac (ecological hypoxic tolerance) along a transect in the Pacific Ocean based on a climatology from the World Ocean Atlas dataset (Garcia et al., 2014). The percentile values of each trait are as follows: 10th (Eo=0.04, Ac=4.8), 50th (Eo=0.34, Ac=7.4) and 90th (Eo=0.79, Ac=17.0). The lower panels show pO2 and temperature from the WOA dataset. Note that the color bar range differs by panel and according to values where is omitted; thus, the color shows only areas where an active metabolic state can be sustained.
2.3 Earth system model simulations
This study is based on the CESM1-LE, described in detail by Kay et al. (2015). The CESM1-LE included 34 ensemble members integrated from 1920–2100 under historical and RCP8.5 forcing. The ensemble was generated by adding round-off level (10−14 K) perturbations to the air temperature field at initialization in 1920; this small difference yields rapidly diverging model solutions due to the chaotic dynamics intrinsic to the climate system, thus developing an ensemble spread that is representative of internal variability (Kay et al., 2015). Briefly, the CESM1-LE uses the Community Earth System Model, version 1 (Hurrell et al., 2013), with a horizontal resolution of (nominally) 1° in all components. The ocean component is the Parallel Ocean Program, version 2 (Smith et al., 2010), with sea ice simulated by the Los Alamos Sea Ice Model, version 4 (Hunke and Lipscomb, 2010). Ocean biogeochemistry was represented by the Biogeochemical Elemental Cycling (BEC) model (Moore et al., 2013; Lindsay et al., 2014).
Our analysis focuses on three depths: 50 m, representing near-surface dynamics; the epipelagic zone at 200 m; and the mesopelagic zone at 500 m. pO2 was calculated using the Garcia and Gordon (1992) solubility formulation. For convenience, we use the period 1920–1965 to define a minimally perturbed natural state as this period is prior to the development of substantial anthropogenic trends in ocean oxygen and temperature (Long et al., 2016). We also examine distributions over the last 3 decades of the 21st century (2070–2099) to evaluate the projected climate change signal under RCP8.5. We use the mean across all 34 ensemble members to quantify the deterministic, “forced” response of the climate system to anthropogenic influence (Deser et al., 2012). The ensemble spread is thus indicative of the amplitude of variations attributable to natural variability.
We compared the CESM1-LE (1920–1965) with the World Ocean Atlas, version 2013 (WOA2013), dataset (Garcia et al., 2014), an observationally based, gridded climatology (Fig. 3a–i). CESM1-LE generally provides a reasonable representation of pO2 and temperature distributions at the selected depths (Fig. 3); however, there are important biases to acknowledge in the context of interpreting the results. Temperature magnitudes are generally well simulated in the CESM1-LE, showing a root mean square error (RMSE) of <1.3 °C and a pattern correlation coefficient of (PCC) >0.98 in all three selected depths (50, 200 and 500 m) (Table 1). Temperature magnitudes are slightly underestimated at 50 and 200 m (mean bias of <0.3 °C) and are overestimated by 0.41 °C at 500 m. Note that, since our comparison uses CESM1-LE data from the period 1920–1965, some discrepancy in temperature might be expected from the signal of climate warming present in the WOA observations. pO2 is also reasonably well captured by the CESM1-LE (PCC<0.95), but magnitudes are slightly underestimated at depth, showing a mean bias of −1.63 and −2.1 kPa at 200 and 500 m with respect to WOA13 (Table 1). Regions of low-pO2 waters are too extensive in CESM1-LE (Fig. 3n and o), and there is a slight degradation of skill with depth for pO2 fields (Table 1). The underestimation of pO2 leads to a slight underestimation of Φ′ with respect to WOA13 and an overestimation of habitat loss in the future climate (Fig. 3p–r); however, Φ′ computed from the model fields demonstrates that the dominant spatial patterns are well captured by the CESM1-LE despite magnitudes that are slightly too low. This CESM pO2 bias is common among coarse-resolution ocean models, and it is attributed to a sluggish circulation and, hence, weak ventilation (Long et al., 2016). These differences ultimately matter most near the hypoxic zones and at the boundaries of habitable zones like the oxygen minimum zones (OMZs).
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f03](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f03-thumb.png)
Figure 3Mean-state comparison with observations. The climatological mean of (a–c, j–l) temperature (°C), (d–f, m–o) pO2 (kPa) and (g–i, p–r) metabolic index for active metabolism (Φ′) for the median ecotype (Eo=0.34, Ac=7.4); three depths are shown: (a, d, g, j, m, p) 50 m, (b, e, h, k, n, q) 200 m and (c, f, i, l, o, r) 500 m. Panels (a–i) show the WOA13 dataset, and panels (j–r) show CESM1-LE.
3.1 Joint temperature–pO2 natural variability and forced trends
The spatial distribution of the number of viable ecotypes is shown in Fig. 4 for the “unperturbed” climate (1920–1965). Our intention here is not to quantify the actual biogeographic range of organisms in the environment but rather to illustrate the ocean's ability to support respiration by marine ectotherms given the metabolic capacities afforded within the trait space of extant organisms. High-latitude environments do not impose strong aerobic constraints (cold-intolerance notwithstanding); thus, over much of the Southern Ocean, North Atlantic and Arctic Ocean, almost all 61 ecotypes can sustain respiration. The tropical oceans impose the strongest aerobic constraints, restricting the viability of ecotypes that do not have a high hypoxia tolerance (Ao). For example, less than 25 ecotypes are viable over much of the tropical surface ocean (Fig. 4a); low concentrations of oxygen at depth impose even stronger constraints, and no ecotypes are viable in the core of OMZs (Fig. 4b and c). The spatial patterns of the number of viable ecotypes are tightly controlled by temperature at the surface since pO2 is mostly near saturated levels; at depth, however, pO2 is the dominant driver of geographic patterns in ecotype viability (Figs. 2–4). Temperature generally decreases with depth, reducing the metabolic oxygen demand. However, since pO2 also decreases with depth and displays greater lateral heterogeneity, pO2 emerges as the dominant constraint of spatial structure in ecotype viability at depth.
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f04](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f04-thumb.png)
Figure 4Metabolic constraints on trait–space viability. (a–c) The number of ecotypes from the physiological trait database that are viable (total=61) in the CESM1-LE over the period 1920–1965. (d–f) The standard deviation (expressed as a percent of the mean) in the number of viable ecotypes, reflecting fluctuations driven by natural variability.
The standard deviation of annual anomalies using all CESM1-LE ensemble members provides insight into the amplitude of natural variability (Fig. 5, 1 standard deviation). Temperature and pO2 show similar patterns of natural variability in the upper ocean, both showing particularly large variance in the western tropical Pacific and Indian Ocean (Fig. 5a and d). Spatial variation in the magnitude of temperature variability generally decreases with depth, but pO2 displays relatively larger variability at depth with respect to the surface in some regions (Fig. 5a–f). The joint pO2–temperature variability manifests in variations of Φ′ (Fig. 5g–i). Natural variability in Φ′ computed for the median ecotype shows spatial patterns similar to those of temperature in the upper-surface ocean (50 m), but it is more similar to pO2 at depth. Thus, variations in Φ′ tend to be temperature-dominated near the surface but are more strongly controlled by pO2 variability at depth. Φ′ also shows the most extensive natural variability at 200 m, consistently with the variability of pO2. The number of viable species shows more dramatic fluctuations than variations in the median-ecotype Φ′; variations in the number of viable ecotypes exceed 30 % on annual timescales in the tropical upper ocean and near OMZ boundaries in the water column (Fig. 4c and d). This reflects the fact that interannual variability can preclude habitability for some regions of the Ac–Eo trait space, but these variations do not necessarily impact viability for the median ecotype (Fig. 1). In the tropical surface ocean, high temperatures (>25 °C) and saturated surfaces (pO2>20 kPa) require a high hypoxia tolerance (Ac) but permit a range of Eo values (Figs. 1b, 2a and b). Ecotypes with larger temperature sensitivity (high Eo) are particularly responsive to variations in temperature.
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f05](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f05-thumb.png)
Figure 5The amplitude of natural variability in the ocean's metabolic state. The panels show the standard deviation of the annual-mean anomalies of all ensemble members over the period 1920–1965 for (a–c) temperature (°C), (d–f) pO2 (kPa) and (g–i) the metabolic index (unitless) of the median ecotype (Eo=0.34, Ac=7.4).
CESM1-LE simulates nearly homogeneous warming between 1920–1965 and 2070–2099 in the surface ocean (50 m) under RCP8.5, with an exception being the so-called North Atlantic warming hole (Fig. 6a). Both modeling and observational studies have linked the North Atlantic warming hole to the slowing of the Atlantic overturning circulation with climate change (Keil et al., 2020). The magnitude of ocean warming generally diminishes with depth except in the North Atlantic, where, despite reductions, the overturning circulation effectively propagates anthropogenic heat anomalies into the ocean interior. pO2 shows heterogeneous changes between 1920–1965 and 2070–2099 (Fig. 6d–f). In the upper ocean, pO2 changes are generally small (<1 kPa) because the near-surface is kept close to saturation via photosynthetic oxygen production and air–sea equilibration. At depth, however, pO2 shows long-term changes linked to the accumulated effects of respiration and changes in circulation (Ito et al., 2017). At 200 m for example, the Pacific Ocean displays a basin-wide mean reduction in pO2 of 2 kPa (∼30 %), while the Atlantic and Indian basins gain about >2 kPa (∼10 %–35 %) by the end of the century. The largest long-term pO2 loss (>3 kPa) occurs in the North Pacific, while the largest pO2 gain (∼2 kPa) occurs in the North Atlantic gyre and western Indian Ocean (Fig. 6e and f).
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f06](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f06-thumb.png)
Figure 6Net long-term change (2070–2099 minus 1920–1965) in the CESM1-LE ensemble mean temperature (a–c) and (d–f) pO2 at 50, 200 and 500 m.
Figure 7 shows the relationship between interannual variations in pO2 versus temperature (pO2-T) in the unperturbed climate (1920–1965; top row) and for the forced trend associated with 21st-century climate change (2070–2099 minus 1920–1965; bottom row). The nature of the pO2–T relationship is an important indicator of the impacts of variability on the metabolic state. Furthermore, the extent to which the forced trend is characterized by a pO2–T relationship that is distinct from that associated with natural variability provides insight into the potential for advanced or delayed detection of signals in Φ relative to pO2 or temperature alone. Given that metabolic rates for most organisms increase with temperature (positive Eo), a positive correlation between variations in temperature and pO2 is generally indicative of compensating changes, wherein increased oxygen demand is at least partially offset by increased supply. Anticorrelation between temperature and pO2, by contrast, will generally be associated with compounding impacts on the metabolic index as a negative correlation indicates that reductions in pO2 (i.e., oxygen supply) accompany warming (i.e., increased demand). The sign of the pO2–T relationship in the natural climate varies regionally and with depth (Fig. 7, top row). The surface ocean is generally characterized by a weak, positive pO2–T relationship, which could manifest from, among other mechanisms, temperature-induced increases in photosynthetic oxygen production (Fig. 7a). The natural pO2–T relationship in the epipelagic region (200 m) is characterized by strong positive correlations in the tropics and negative correlations at high latitudes (Fig. 7b). A positive correlation between pO2 and temperature at this depth could be induced by variability associated with adiabatic vertical displacement of isopycnals or “heave”, which has the effect of translating background gradients in properties vertically in the water column (Long et al., 2016). Upward movement of a deep isopycnal surface would yield a negative temperature anomaly and a negative pO2 anomaly (positive correlation) as the deeper, colder waters have greater oxygen utilization signatures associated with longer ventilation ages. Negative correlations between pO2 and temperature could manifest from ventilation processes, where enhanced subduction of surface water yields anomalously cold water masses that are enriched in oxygen. The sign of these epipelagic pO2–T correlations shows some similarity to that associated with the externally forced climate (Fig. 6e), but the latter is characterized by a greater prevalence of anticorrelation, most notably in the North Pacific Ocean. At 500 m depth, the relationship between temperature and pO2 in the natural climate is almost a mirror image of the epipelagic region (Fig. 7c); the tropics generally display negative correlations, while polar regions show positive correlations (Fig. 7e). The pO2–T relationship in the forced trend at 500 m is dominated by broad regions of deeply negative correlations, with the most pronounced effect again being in the North Pacific. The negative relationship is consistent with a ventilation signal as buoyancy-induced stratification from warming curtails the introduction of new oxygen into the ocean interior. The predominantly negative pO2–T relationship associated with the forced trend is indicative of the compounding effects of climate change on metabolic state, increasing metabolic demand while simultaneously reducing oxygen supply.
3.2 Long-term habitat changes
Figure 8 shows the climate-driven changes in Φ′ for the median ecotype, as well as the impacts of climate change on the number of viable ecotypes. Notably, while pO2 in the near-surface ocean is relatively insensitive to climate change (Fig. 6d), there are reductions in Φ′ in the tropics (Fig. 9d) owing to the direct impacts of warming. These changes are associated with deep reductions in the number of viable ecotypes in the tropics (Fig. 8a). There are modest increases in Φ′ and ecotype viability at high latitudes; metabolic state in these regions is affected by cold intolerance, and, thus, warming broadens the viable region of trait space. Additionally, sea ice melt supports an increase in pO2 as gas exchange becomes more effective at restoring equilibrium oxygen concentrations. The number of viable ecotypes shows more intense patterns than those in the median-ecotype Φ′ in the upper ocean (Fig. 8). This is partly because ecotypes predicted to lose viability in the tropical regions (∼50 %) are at the extremes of the Ac-Eo distribution (Fig. 1) and are not captured by the median-ecotype Φ′. Nevertheless, outside the tropical regions, the median ecotype gives a good indication of the anthropogenic impact on marine ectotherms. The projected habitat loss in the epipelagic–pelagic North Pacific (>50 %) and habitat gain in the epipelagic–pelagic Southern Indian Ocean (∼40 %) and pelagic western tropical regions (∼40 %) are consistent with a decrease in the median-ecotype Φ′. Note that the most pronounced effects on habitat are associated with regions where climate change drives a strongly negative pO2–temperature relationship (Fig. 7).
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f08](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f08-thumb.png)
Figure 8Net change in the number of habitable ecotypes in percentage (a–c). Net metabolic index change (2070–2099 vs. 1920–1965) for the median ecotype [Eo=0.34, Ac=7.4] (d–f) at 50 m (a, d), 200 m (b, e) and 500 m (c, f).
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f09](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f09-thumb.png)
Figure 9Time of emergence (ToE) of the climate forcing signal for (a) temperature, (b) pO2, (c) the metabolic index of the median ecotype [Eo=0.34, Ac=7.4] and (d) the number of viable ecotypes for a single model grid in the North Pacific at 200 m. ToE (green arrows) is defined as the time when the forced trend signal (ensemble member time series) is above 2 standard deviations (dotted black line) compared to all ensemble members for the period 1920–1965.
3.3 Time of emergence
In this section, we examine the time of emergence (ToE, Hawkins and Sutton, 2012), the point when forced changes in pO2, temperature and Φ′ can be distinguished from the background natural variability. We define ToE as the time when the magnitude of change in the ensemble mean of a particular variable exceeds 2 standard deviations of the natural climate (1920–1965). This is illustrated in Fig. 9 for a single grid point in the North Pacific at 200 m. At this location, the forced trend in temperature shows a monotonic increase, while pO2 shows a monotonic decrease; as a result, Φ′ for the median ecotype and the number of viable ecotypes decrease over time. The anti-correlation between pO2 and temperature exacerbates trends in Φ′, and, hence, the forced trend of the median-ecotype Φ′ emerges from natural noise earlier than either pO2 or temperature do alone (Fig. 10a–c). Note that, although the ToE of ecotype viability change is directly derived from changes in Φ′, it is binary-counted; changes in ecotype viability are counted in whole numbers, and this creates a step-function temporal–spatial variation (Fig. 9d). Consequently, this step-function-like feature of ecotype viability creates discontinuities, even in spatial patterns of ToE (Fig. 10j–l), as also shown in the natural variance in Fig. 4d–f.
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f10](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f10-thumb.png)
Figure 10Time of emergence (ToE) of the climate forcing signal for temperature, pO2, ϕ and the number of viable ecotypes. ToE is defined as the time when the forced trend signal (ensemble member time series) is above 2 standard deviations compared to all ensemble members for the period 1920–1965.
The ToE values of pO2 and temperature are inverted with depth; temperature emerges earliest in the upper ocean, while pO2 emerges earlier at depth, while in the upper ocean, it emerges later or shows no emergence (Fig. 10a–f). This feature is consistent with larger long-term upper-ocean temperature changes and greater pO2 changes at depth. Near-surface ocean temperature is shown to have mostly already emerged by 2020 and is predicted to have almost completely emerged by the late 2060s under RCP85 (Fig. 10a–c). The early emergence of temperature from natural noise also persists for regions of relatively low natural variance at depth, e.g., the Southern Ocean and Atlantic basin gyres. However, regions of the largest natural variability (see Fig. 5) like the subtropical–subpolar Pacific do not emerge until close to the end of the century. For pO2, anthropogenic changes in the upper ocean generally do not emerge from natural noise before the end of the century, except for the Arctic Ocean and eastern Antarctic. In the Arctic Ocean and eastern Antarctic, pO2 gain is related to sea melt emergence by the mid-2050s (Fig. 10a). The median-ecotype Φ′ ToE shows spatial patterns that are coherent with temperature ToE in the upper ocean, with the exception of polar regions. In contrast, these are consistent with pO2 ToE patterns at depth; this is consistent with net long-term Φ′ changes in Fig. 9d. The emergence of the anthropogenic signal in ecotype viability closely resembles the median-ecotype Φ′ spatial patterns but with non-harmonious spatial patterns due to the step-function-like counting feature of viability changes. It shows that the predicted ∼50 % ecotype viability loss in the tropics (Fig. 6a) may already be distinguishable from natural variability by the mid-2030s. In the North Pacific, the predicted >50 % ecotype viability loss in the epipelagic–pelagic regions is predicted to start emerging in the 2040s at 500 m and in the 2080s at 200 m (Fig. 10k and l).
In summary, we showed that, because of the surface ocean’s large warming signal and the lowest pO2 loss outside of the polar regions in the RCP85 climate scenario, the surface ocean is characterized by habitat loss in the tropics and a slight habitat gain in the polar regions (Fig. 11). Sea ice melt supports oxygen gain through the enhancement of temperature-driven solubility in the surface polar regions. At depth, warming is less prevalent by the end of the 21st century; however, oxygen loss related to the weakening ventilation of the ocean interior as the ocean becomes more stratified has a stronger impact on metabolic reliance, leading to habitat loss in the tropics and subtropics. On the other hand, cooler temperatures and efficient ventilation in polar regions create an oxygen-rich environment. Thus, in contrast to tropical and subtropical regions, warming leads to a slight habitat gain (Fig. 11) as organisms escape the cold intolerance imposed by molecular gas diffusion at low temperatures.
![https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f11](https://bg.copernicus.org/articles/21/3477/2024/bg-21-3477-2024-f11-thumb.png)
Figure 11Summary figure: the distribution of ecotype viability within representative ocean temperature and pO2 boundaries for the 66 species analyzed in this study. The markers represent the subsampled regions, with polar regions denoted by circles, tropical regions denoted by squares and subtropical regions denoted by stars. The colors represent the depth levels; 50 m (gray), 200 m (purple) and 500 m (red). Each arrow shows the estimated joint temperature–pO2 climate change vector based on the net changes in temperature and pO2 (as depicted in Fig. 6).
The human-induced rapid warming of the planet has been shown to drive ocean deoxygenation (Ito et al., 2017; Schmidtko et al., 2017; Long et al., 2016). Higher metabolic oxygen demand at higher temperatures (Gillooly et al., 2001; Deutsch et al., 2015, 2022) raises concerns about the ability of marine ectotherms to support aerobic respiration in the future. This study set out to characterize the anticipated climate change signal in the ocean's metabolic state in the context of natural variability using metabolic theory as a basis to examine the capacity of the environment to support ectothermic marine heterotrophs.
The spatial variation in pO2 and temperature in the unperturbed natural climate state set biogeographic boundaries based on ectotherms' physiological performances. The resilience of these ectotherms' biogeographic structures to natural variability and long-term climate warming is perturbed by the joint pO2–temperature changes, effectively measured by the metabolic index (Φ′). An increase in the capacity of the organisms to support aerobic respiration increases Φ′ – for example, by ocean cooling or increases in oxygen supply. On the contrary, ocean warming and decreases in oxygen supply decrease Φ′. There are exceptions in extremely low-temperature environments (Fig. 11), where aerobic respiration is also limited by kinematic gas transfer into the organism in addition to environmental oxygen supply. Relative changes in pO2 and temperature in the natural variability and forced trend therefore regulate ectotherms' resilience to environmental changes. Under the RCP85 climate scenario, the ocean generally warms homogeneously, but concurrent pO2 changes are heterogeneous and vary with depth. Thus, the characteristics of these pO2–temperature forced trend changes determine when the climate change impact on marine ectotherms can be distinguishable from natural variability.
In the surface ocean, pO2 is generally abundant and relatively uniform; thus, spatial temperature variations have a dominant constraint on the spatial variations of the organismic metabolic state. The warmest parts of the surface ocean, the tropical oceans, can only support about 10–20 (∼30 %) of the 61 ecotypes, while cooler regions in the extra tropics have nearly 100 % viability. Moreover, since warming anomalies propagate from the surface, the surface tropical oceans also show the largest natural variance in temperature and ecotype viability. This is because extremely warm temperatures in the surface tropics (>25 °C) are mainly suited for organisms with a high temperature sensitivity (Eo), which are relatively fewer and are mostly close to their physiological limits (Storch et al., 2014). Large natural variability in the warmest parts of the tropical surface ocean precludes the forced trend signal from emerging from the natural variability in the ecotype viability by end of the century, although the ocean warms the most at the surface. Nevertheless, the large warming trends in the surface ocean generally emerge relatively early (the 2020s) from natural variability in both temperature and ecotype viability in most regions. Minimal changes in surface pO2 in the forced trend affirm that surface ocean marine ectotherms are mainly perturbed by temperature in the context of anthropogenic changes. In polar regions, warming has a counterintuitive effect on marine ectotherms with respect to most parts of the surface ocean. There, warming helps organisms escape extreme cold intolerances by enhancing membrane kinematic gas transfer, which enhances Φ′ and thus ecotype richness in the future (Fig. 11)
In the epipelagic and mesopelagic regions (200 and 500 m), the forced temperature trend and natural variability are broadly smaller than the surface ocean, while pO2 changes show the opposite. Thus, at depth, pO2 plays a more intricate role in perturbating marine ectotherm habitats in the context of anthropogenic warming with respect to the surface ocean, where temperature plays a dominant role. Contrasting the regression between pO2 and temperature in the natural climate, forced trends provide an instructive framework to analyze ectotherms' long-term changes. Regions showing different correlations between temperature and pO2 in the forced trends in comparison to the natural climate suggest a loss of metabolic resilience and a loss of habitat, and these regions tend to have a relatively early ToE values. For instance, in the epipelagic and mesopelagic North Pacific, temperature–pO2 regressions switched from a positive correlation in the unperturbed climate to a strong negative correlation in the forced trend (Fig. 7). The North Pacific pelagic–epipelagic region is projected to lose nearly half of the present climate ecotype viability by the end of the 21st century, with the projected habitat loss starting to emerge by the late 2030s under the RCP85 climate scenario. On the other hand, in the Arctic Ocean and some parts of the Southern Ocean, same-sign pO2–temperature correlations in the forced trends result in the preservation of the marine habitat and even slight enhancements.
The joint temperature–oxygen metabolic framework in this study provides additional insights into the impact of climate change on marine ecosystems in comparison to the independent oxygen or temperature analysis. Here, we showed that, while warming is the leading-order driving mechanism of climate change, the direct effect of warming on marine ecosystems is mostly in the upper ocean. Climate-change-related oxygen loss is a major driver of marine ecosystem stress in addition to warming at depth. Incorporating organismal physiological sensitivity to oxygen–temperature changes in the metabolic framework provides insights into how climate impacts the biogeographic structure of marine habitats. We find that forced perturbations to pO2 and temperature will strongly exceed those associated with the natural system in many parts of the upper ocean, mostly pushing organisms in these environments closer to or beyond their physiological limits. Climate warming is expected to drive significant marine habitat loss in the surface tropical oceans and epipelagic–pelagic North Pacific basin while gaining marginal habitat viability in the surface Arctic Ocean and some parts of the Southern Ocean.
Data from the CESM1-LE are available through the National Center for Atmospheric Research (NCAR) Climate Data Gateway and instructions for data access are available from https://doi.org/10.5065/d6j101d1 (Kay and Deser, 2016), last accessed 25 July 2024.
The supplement related to this article is available online at: https://doi.org/10.5194/bg-21-3477-2024-supplement.
PM and ML designed the study approach. PM developed the analysis with feedback from all the co-authors. PM prepared the paper with contributions from all the co-authors.
The contact author has declared that none of the authors has any competing interests.
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.
This article is part of the special issue “Low-oxygen environments and deoxygenation in open and coastal marine waters”. It is not associated with a conference.
We would like to acknowledge the data access and computing support provided by the NCAR Cheyenne HPC.
Precious Mongwe, Matthew Long, Curtis Deutsch and Takamitsu Ito were funded by the National Science Foundation (NSF) under grant agreement no. 1737158. Precious Mongwe and Yeray Santana-Falcón were also funded by the European Union's Horizon 2020 research and innovation programme under grant agreement no. 820989 (COMFORT).
This paper was edited by Mike Roman and reviewed by three anonymous referees.
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, 46, https://doi.org/10.1126/science.aam7240, 2018.
Deser, C., Phillips, A., Bourdette, V., and Teng, H.: Uncertainty in climate change projections: The role of internal variability, Clim. Dynam., 38, 527–546, https://doi.org/10.1007/s00382-010-0977-x, 2012.
Deutsch, C., Ferrel, A., Seibel, B., Pörtner, H. O., and Huey, R. B.: Climate change tightens a metabolic constraint on marine habitats, Science, 348, 1132–1135, https://doi.org/10.1126/science.aaa1605, 2015.
Deutsch, C., Penn, J. L., and Seibel, B.: Metabolic trait diversity shapes marine biogeography, Nature, 585, 557–562, https://doi.org/10.1038/s41586-020-2721-y, 2020.
Deutsch, C., Penn, J. L., Verberk, W. C. E. P., Inomura, K., Endress, M.-G., and Payne, J. L.: Impact of warming on aquatic body sizes explained by metabolic scaling from microbes to macrofauna, P. Natl. Acad. Sci. USA, 119, e2201345119, https://doi.org/10.1073/pnas.2201345119, 2022.
Garcia, H. E. and Gordon, L. I.: Oxygen solubility in seawater: Better fitting equations, Limnol. Oceanogr., 37, 1307–1312, https://doi.org/10.4319/lo.1992.37.6.1307, 1992.
Garcia, H. E., Boyer, T. P., Locarnini, R. A., Antonov, J. I., Mishonov, A. V., Baranova, O. K., Melissa, M. Z., Reagan, J. R., and Johnson, D. R.: World Ocean Atlas 2013, Vol. 3: Dissolved Oxygen, Apparent Oxygen Utilization, and Oxygen Saturation, NOAA Atlas NESDIS 75, https://doi.org/10.7289/V5XG9P2W, 2014.
Gillooly, J., Brown, J., West, G., Savage, V., Charnov, E., Gillooly, J. F., Brown, J. H., West, G. B., Savage, V. M., and Charnov, E. L.: Effects of size and temperature on metabolic rate, Science, 293, 2248–2251, https://doi.org/10.1126/science.1061967, 2001.
Hawkins, E. and Sutton, R.: Time of emergence of climate signals, Geophys. Res. Lett., 39, L01702, https://doi.org/10.1029/2011GL050087, 2012.
Hoegh-Guldberg, O. and Bruno, J. F.: The Impact of Climate Change on the World's Marine Ecosystems, Science, 328, 1523–1528, https://doi.org/10.1126/science.1189930, 2010.
Howard, E. M., Penn, J. L., Frenzel, H., Seibel, B. A., Bianchi, D., Renault, L., Kessouri, F., Sutula, M. A., Mcwilliams, J. C., and Deutsch, C.: Climate-driven aerobic habitat loss in the California Current System, Sci. Adv., 6, eaay3188, https://doi.org/10.1126/sciadv.aay3188, 2020.
Hunke, E. C. and Lipscomb, W. H.: CICE: the Los Alamos Sea Ice Model Documentation and Software User's Manual Version 4.1, edited by: Hunke, E. C., Lipscomb, W. H., Turner, A. K., Jeffery, N., and Elliott, S., Los Alamos National Laboratory, Los Alamos, LA-CC-06-012, 2010.
Hurrell, J. W., Holland, M. M., Gent, P. R., Ghan, S., Kay, J. E., Kushner, P. J., Lamarque, J. F., Large, W. G., Lawrence, D., Lindsay, K., Lipscomb, W. H., Long, M. C., Mahowald, N., Marsh, D. R., Neale, R. B., Rasch, P., Vavrus, S., Vertenstein, M., Bader, D., Collins, W. D., Hack, J. J., Kiehl, J., and Marshall, S.: The community earth system model: A framework for collaborative research, B. Am. Meteorol. Soc., 94, 1339–1360, https://doi.org/10.1175/BAMS-D-12-00121.1, 2013.
Ito, T. and Deutsch, C.: A conceptual model for the temporal spectrum of oceanic oxygen variability, Geophys. Res. Lett., 37, L03601, https://doi.org/10.1029/2009GL041595, 2010.
Ito, T., Minobe, S., Long, M. C., and Deutsch, C.: Upper ocean O2 trends: 1958–2015, Geophys. Res. Lett., 44, 4214–4223, https://doi.org/10.1002/2017GL073613, 2017.
Kay, J. E. and Deser, C.: Large Ensemble Community Project (LENS), NCAR [data set], https://doi.org/10.5065/d6j101d1, 2016.
Kay, J. E., Deser, C., Phillips, A., Mai, A., Hannay, C., Strand, G., Arblaster, J. M., Bates, S. C., Danabasoglu, G., Edwards, J., Holland, M., Kushner, P., Lamarque, J. F., Lawrence, D., Lindsay, K., Middleton, A., Munoz, E., Neale, R., Oleson, K., Polvani, L., and Vertenstein, M.: The community earth system model (CESM) large ensemble project: A community resource for studying climate change in the presence of internal climate variability, B. Am. Meteorol. Soc., 96, 1333–1349, https://doi.org/10.1175/BAMS-D-13-00255.1, 2015.
Keeling, R. F., Körtzinger, A., and Gruber, N.: Ocean deoxygenation in a warming world, Annu. Rev. Mar. Sci., 2, 199–229, https://doi.org/10.1146/annurev.marine.010908.163855, 2010.
Keil, P., Mauritsen, T., Jungclaus, J., Hedemann, C., Olonscheck, D., and Ghosh, R.: Multiple drivers of the North Atlantic warming hole, Nat. Clim. Change, 10, 667–671, https://doi.org/10.1038/s41558-020-0819-8, 2020.
Lindsay, K., Bonan, G. B., Doney, S. C., Hoffman, F. M., Lawrence, D. M., Long, M. C., Mahowald, N. M., Moore, J. K., Randerson, J. T., and Thornton, P. E.: Preindustrial-control and twentieth-century carbon cycle experiments with the Earth system model CESM1(BGC), J. Climate, 27, 8981–9005, https://doi.org/10.1175/JCLI-D-12-00565.1, 2014.
Long, M. C., Deutsch, C., and Ito, T.: Finding forced trends in oceanic oxygen, Global Biogeochem. Cy., 30, 381–397, https://doi.org/10.1002/2015GB005310, 2016.
Moore, J. K., Lindsay, K., Doney, S. C., Long, M. C., and Misumi, K.: Marine ecosystem dynamics and biogeochemical cycling in the community earth system model [CESM1(BGC)]: Comparison of the 1990s with the 2090s under the RCP4.5 and RCP8.5 scenarios, J. Climate, 26, 9291–9312, https://doi.org/10.1175/JCLI-D-12-00566.1, 2013.
Oschlies, A., Brandt, P., Stramma, L., and Schmidtko, S.: Drivers and mechanisms of ocean deoxygenation, Nat. Geosci., 11, 467–473, https://doi.org/10.1038/s41561-018-0152-2, 2018.
Penn, J. L., Deutsch, C., Payne, J. L., and Sperling, E. A.: Temperature-dependent hypoxia explains biogeography and severity of end-Permian marine mass extinction, Science, 362, 6419, https://doi.org/10.1126/science.aat1327, 2018.
Piiper, J., Dejours', P., Haab, P., and Rahn, H.: Concepts and Basic Quantities in Gas Exchange Physiology, Resp. Physiol., 13, 292–304, 1971.
Pörtner, H. O.: Climate variations and the physiological basis of temperature dependent biogeography: systemic to molecular hierarchy of thermal tolerance in animals, Comp. Biochem. Phys. A, 132, 739–761, 2002.
Pozo Buil, M. and Di Lorenzo, E.: Decadal dynamics and predictability of oxygen and subsurface tracers in the California Current System, Geophys. Res. Lett., 44, 4204–4213, https://doi.org/10.1002/2017GL072931, 2017.
Rodgers, K. B., Lin, J., and Frölicher, T. L.: Emergence of multiple ocean ecosystem drivers in a large ensemble suite with an Earth system model, Biogeosciences, 12, 3301–3320, https://doi.org/10.5194/bg-12-3301-2015, 2015.
Rosewarne, P. J., Wilson, J. M., and Svendsen, J. C.: Measuring maximum and standard metabolic rates using intermittent-flow respirometry: A student laboratory investigation of aerobic metabolic scope and environmental hypoxia in aquatic breathers, J. Fish Biol., 88, 265–283, https://doi.org/10.1111/jfb.12795, 2016.
Schlunegger, S., Rodgers, K. B., Sarmiento, J. L., Frölicher, T. L., Dunne, J. P., Ishii, M., and Slater, R.: Emergence of anthropogenic signals in the ocean carbon cycle, Nat. Clim. Change, 9, 719–725, https://doi.org/10.1038/s41558-019-0553-2, 2019.
Schmidtko, S., Stramma, L., and Visbeck, M.: Decline in global oceanic oxygen content during the past five decades, Nature, 542, 335–339, https://doi.org/10.1038/nature21399, 2017.
Smith, R., Jones, P., Briegleb, B., Bryan, F., Danabasoglu, G., Dennis, J., Dukowicz, J., Eden, C., Fox-Kemper, B., Gent, P., Hecht, M., Jayne, S., Jochum, M., Large, W., Lindsay, K., Maltrud, M., Norton, N., Peacock, S., Vertenstein, M., and Yeager, S.: The Parallel Ocean Program (POP) Reference Manual Ocean Component of the Community Climate System Model (CCSM) and Community Earth System Model (CESM) 1, http://n2t.net/ark:/85065/d70g3j4h (last access: 26 July 2024), 2010.
Storch, D., Menzel, L., Frickenhaus, S., and Pörtner, H. O.: Climate sensitivity across marine domains of life: Limits to evolutionary adaptation shape species interactions, Global Change Biol., 20, 3059–3067, https://doi.org/10.1111/gcb.12645, 2014.
Tiano, L., Garcia-Robledo, E., Dalsgaard, T., Devol, A. H., Ward, B. B., Ulloa, O., Canfield, D. E., and Revsbech, N. P.: Oxygen distribution and aerobic respiration in the north and south eastern tropical Pacific oxygen minimum zones, Deep-Sea Res. Pt. I, 94, 173–183, https://doi.org/10.1016/j.dsr.2014.10.001, 2014.
Vaquer-Sunyer, R. and Duarte, C. M.: Thresholds of hypoxia for marine biodiversity, P. Natl. Acad. Sci. USA, 105, 15452–15457, 2008.