Quantification of ocean heat uptake from changes in atmospheric O2 and CO2 composition

The ocean is the main source of thermal inertia in the climate system. Ocean heat uptake during recent decades has been quantified using ocean temperature measurements. However, these estimates all use the same imperfect ocean dataset and share additional uncertainty due to sparse coverage, especially before 2007. Here, we provide an independent estimate by using measurements of atmospheric oxygen (O2) and carbon dioxide (CO2) – levels of which increase as the ocean warms and releases gases – as a whole ocean thermometer. We show that the ocean gained 1.29 ± 0.79 × 1022 Joules of heat per year between 1991 and 2016, equivalent to a planetary energy imbalance of 0.80 ± 0.49 W watts per square metre of Earth’s surface. We also find that the ocean-warming effect that led to the outgassing of O2 and CO2 can be isolated from the direct effects of anthropogenic emissions and CO2 sinks. Our result – which relies on high-precision O2 atmospheric measurements dating back to 1991 – leverages an integrative Earth system approach and provides much needed independent confirmation of heat uptake estimated from ocean data

Simulations of ocean deoxygenation in the historical era: insights from forced and coupled models

Ocean deoxygenation due to anthropogenic warming represents a major threat to marine ecosystems and fisheries. Challenges remain in simulating the modern observed changes in the dissolved oxygen (O2). Here, we present an analysis of upper ocean (0-700m) deoxygenation in recent decades from a suite of the Coupled Model Intercomparison Project phase 6 (CMIP6) ocean biogeochemical simulations. The physics and biogeochemical simulations include both ocean-only (the Ocean Model Intercomparison Project Phase 1 and 2, OMIP1 and OMIP2) and coupled Earth system (CMIP6 Historical) configurations. We examine simulated changes in the O2 inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O2 trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O2 inventory for each of the MIP simulations over the past 5 decades is 0.03 ± 0.39×1014 [mol/decade] for OMIP1, −0.37 ± 0.15×1014 [mol/decade] for OMIP2, and −1.06 ± 0.68×1014 [mol/decade] for CMIP6 Historical, respectively. The trend in the upper ocean global O2 inventory for the latest observations based on the World Ocean Database 2018 is −0.98×1014 [mol/decade], in line with the CMIP6 Historical multi-model mean, though this recent observations-based trend estimate is weaker than previously reported trends. A comparison across ocean-only simulations from OMIP1 and OMIP2 suggests that differences in atmospheric forcing such as surface wind explain the simulated divergence across configurations in O2 inventory changes. Additionally, a comparison of coupled model simulations from the CMIP6 Historical configuration indicates that differences in background mean states due to differences in spin-up duration and equilibrium states result in substantial differences in the climate change response of O2. Finally, we discuss gaps and uncertainties in both ocean biogeochemical simulations and observations and explore possible future coordinated ocean biogeochemistry simulations to fill in gaps and unravel the mechanisms controlling the O2 changes

Impacts of ENSO on air-sea oxygen exchange: observations and mechanisms

Models and observations of Atmospheric Potential Oxygen (APO ≃ O2 + 1.1*CO2) are used to investigate the influence of El Niño Southern Oscillation (ENSO) on air-sea O2 exchange. An atmospheric transport inversion of APO data from the Scripps flask network shows significant interannual variability in tropical APO fluxes that is positively correlated with the Niño3.4 index, indicating anomalous ocean outgassing of APO during El Niño. Hindcast simulations of the Community Earth System Model (CESM) and the Institut Pierre-Simon Laplace (IPSL) model show similar APO sensitivity to ENSO, differing from the Geophysical Fluid Dynamic Laboratory (GFDL) model, which shows an opposite APO response. In all models, O2 accounts for most APO flux variations. Detailed analysis in CESM shows the O2 response is driven primarily by ENSO-modulation of the source and rate of equatorial upwelling, which moderate the intensity of O2 uptake due to vertical transport of low-O2 waters. These upwelling changes dominate over counteracting effects of biological productivity and thermally-driven O2 exchange. During El Niño, shallower and weaker upwelling leads to anomalous O2 outgassing, whereas deeper and intensified upwelling during La Niña drives enhanced O2 uptake. This response is strongly localized along the central and eastern equatorial Pacific, leading to an equatorial zonal dipole in atmospheric anomalies of APO. This dipole is further intensified by ENSO-related changes in winds, reconciling apparently conflicting APO observations in the tropical Pacific. These findings suggest a substantial and complex response of the oceanic O2 cycle to climate variability that is significantly (>50%) underestimated in magnitude by ocean models

Retraction Note: Quantification of ocean heat uptake from changes in atmospheric O<SUB>2</SUB> and CO<SUB>2</SUB> composition

International audienceThis Article has been retracted; see accompanying Retraction Note