Search for: All records

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 (O₂). 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 O₂inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O₂trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O₂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×10¹⁴[mol/decade] for OMIP2, and −1.06 ± 0.68×10¹⁴[mol/decade] for CMIP6 Historical, respectively. The trend in the upper ocean global O₂inventory for the latest observations based on the World Ocean Database 2018 is −0.98×10¹⁴[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 O₂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 O₂. 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 O₂changes.

Pinatubo erupted during the first decadal survey of ocean biogeochemistry, embedding its climate fingerprint into foundational ocean biogeochemical observations and complicating the interpretation of long‐term biogeochemical change. Here, we quantify the influence of the Pinatubo climate perturbation on externally forced decadal and multi‐decadal changes in key ocean biogeochemical quantities using a large ensemble simulation of the Community Earth System Model designed to isolate the effects of Pinatubo, which cleanly captures the ocean biogeochemical response to the eruption. We find increased uptake of apparent oxygen utilization and preindustrial carbon over 1993–2003. Nearly 100% of the forced response in these quantities are attributable to Pinatubo. The eruption caused enhanced ventilation of the North Atlantic, as evidenced by deep ocean chlorofluorocarbon changes that appear 10–15 years after the eruption. Our results help contextualize observed change and contribute to improved constraints on uncertainty in the global carbon budget and ocean deoxygenation.

The oceanic response to recent tropical eruptions is examined in Large Ensemble (LE) experiments from two fully coupled global climate models, the Community Earth System Model (CESM) and the Geophysical Fluid Dynamics Laboratory Earth System Model (ESM2M), each forced by a distinct volcanic forcing dataset. Following the simulated eruptions of Agung, El Chichón, and Pinatubo, the ocean loses heat and gains oxygen and carbon, in general agreement with available observations. In both models, substantial global surface cooling is accompanied by El Niño–like equatorial Pacific surface warming a year after the volcanic forcing peaks. A mechanistic analysis of the CESM and ESM2M responses to Pinatubo identifies remote wind forcing from the western Pacific as a major driver of this El Niño–like response. Following eruption, faster cooling over the Maritime Continent than adjacent oceans suppresses convection and leads to persistent westerly wind anomalies over the western tropical Pacific. These wind anomalies excite equatorial downwelling Kelvin waves and the upwelling of warm subsurface anomalies in the eastern Pacific, promoting the development of El Niño conditions through Bjerknes feedbacks a year after eruption. This El Niño–like response drives further ocean heat loss through enhanced equatorial cloud albedo, and dominates global carbon uptake as upwelling of carbon-rich waters is suppressed in the tropical Pacific. Oxygen uptake occurs primarily at high latitudes, where surface cooling intensifies the ventilation of subtropical thermocline waters. These volcanically forced ocean responses are large enough to contribute to the observed decadal variability in oceanic heat, carbon, and oxygen.