skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: 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 O2inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O2trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O2inventory 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 O2inventory 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 O2inventory 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 O2changes.  more » « less
Award ID(s):
1948599
PAR ID:
10494115
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; « less
Publisher / Repository:
Frontiers Media SA
Date Published:
Journal Name:
Frontiers in Marine Science
Volume:
10
ISSN:
2296-7745
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, Journal of Geophysical Research: Machine Learning and Computation)
    Abstract The ocean oxygen (O2) inventory has declined in recent decades but the estimates of O2trend are uncertain due to its sparse and irregular sampling. A refined estimate of deoxygenation rate is developed using machine learning techniques and biogeochemical Argo array. The source data includes historical shipboard (bottle and CTD‐O2) profiles from 1965 to 2020 and biogeochemical Argo profiles after 2005. Neural network and random forest algorithms were trained using approximately 80% of this data and the remaining 20% for validation. The training data is further divided into 5‐fold decadal groups to perform cross validation and hyperparameter tuning. Through different combinations of algorithm types and predictor variable sets, an ensemble of gridded monthly O2data sets was generated with similar skills (root‐mean‐square error ∼13–18 μmol/kg and R2 ∼ 0.9). The largest errors are found in the oxycline and frontal regions with strong lateral and vertical gradients. The mapping was repeated with shipboard data only and with both shipboard and Argo data. The effect of including Argo data on the estimated global deoxygenation trends has a major impact with an 56% increase while reducing the uncertainty by 40% as measured by the ensemble spread. This study demonstrates the importance of new biogeochemical Argo arrays in relatively data‐poor regions such as the Southern Ocean. 
    more » « less
  2. ; ; ; ; ; ; ; ; ; ; et al (, Biogeosciences)
    Abstract. The global ocean’s oxygen content has declined significantly over the past several decades and is expected to continue decreasing under global warming with far reaching impacts on marine ecosystems and biogeochemical cycling. Determining the oxygen trend, its spatial pattern and uncertainties from observations is fundamental to our understanding of20 the changing ocean environment. This study uses a suite of CMIP6 Earth System Models to evaluate the biases and uncertainties in oxygen distribution and trends due to sampling sparseness. Model outputs are sub-sampled according to the spatial and temporal distribution of the historical shipboard measurements, and an optimal interpolation method is applied to fill data gaps. Sub-sampled results are compared to full model output, revealing the biases in global and basin-wise oxygen content trends. The optimal interpolation underestimates the modeled global deoxygenation trends, capturing approximately25 two-thirds of the full model trends. North Atlantic and Subpolar North Pacific are relatively well sampled, and the optimal interpolation is capable of reconstructing more than 80% of the oxygen trend. In contrast, pronounced biases are found in the equatorial oceans and the Southern Ocean, where the sampling density is relatively low. Optimal interpolation of the historical dataset estimated the global oxygen loss of 1.5% over the past 50 years. However, the ratio of global oxygen trend between the subsampled and full model output, increases the estimated loss rate to 1.7 to 3.1% over the past 50 years, which partially30 overlaps with previous studies. The approach taken in this study can provide a framework for the intercomparison of different statistical gap-fill methods to estimate oxygen content trends and its uncertainties due to sampling sparseness. 
    more » « less
  3. ; ; ; ; ; ; ; ; ; ; et al (, Biogeosciences)
    The global ocean's oxygen content has declined significantly over the past several decades and is expected to continue decreasing under global warming, with far-reaching impacts on marine ecosystems and biogeochemical cycling. Determining the oxygen trend, its spatial pattern, and uncertainties from observations is fundamental to our understanding of the changing ocean environment. This study uses a suite of CMIP6 Earth system models to evaluate the biases and uncertainties in oxygen distribution and trends due to sampling sparseness. Model outputs are sub-sampled according to the spatial and temporal distribution of the historical shipboard measurements, and the data gaps are filled by a simple optimal interpolation method using Gaussian covariance with a constant e-folding length scale. Sub-sampled results are compared to full model output, revealing the biases in global and basin-wise oxygen content trends. The simple optimal interpolation underestimates the modeled global deoxygenation trends, capturing approximately two-thirds of the full model trends. The North Atlantic and subpolar North Pacific are relatively well sampled, and the simple optimal interpolation is capable of reconstructing more than 80% of the oxygen trend in the non-eddying CMIP models. In contrast, pronounced biases are found in the equatorial oceans and the Southern Ocean, where the sampling density is relatively low. The application of the simple optimal interpolation method to the historical dataset estimated the global oxygen loss to be 1.5% over the past 50 years. However, the ratio of the global oxygen trend between the sub-sampled and full model output has increased the estimated loss rate in the range of 1.7% to 3.1% over the past 50 years, which partially overlaps with previous studies. The approach taken in this study can provide a framework for the intercomparison of different statistical gap-filling methods to estimate oxygen content trends and their uncertainties due to sampling sparseness. 
    more » « less
  4. ; ; ; ; (, Journal of Geophysical Research: Oceans)
    Abstract The southeastern Atlantic Ocean is a crucial yet understudied region for the ocean absorption of anthropogenic carbon (Canth). Data from the A12 (2020) and A13.5 (2010) cruises offer an opportunity to examine changes in dissolved inorganic carbon (DIC), its stable isotope (δ13C), and Canthover the past decade within a limited region (1∼3°E, 32∼42°S). For the decade of 2010–2020, Canthinvasion was observed from the sea surface down to 1,200 m based on both DIC and δ13C data. The mean Canthincrease rate (1.08 ± 0.26 mol m−2 yr−1) during this period accelerated from 0.87 ± 0.05 mol m−2 yr−1during the previous period (1983/84–2010). The δ13C‐based Canthincrease closely matches the DIC‐based estimation below 500 m but is 26% higher in the upper ocean. This discrepancy is likely due to δ13C's longer air‐sea exchange timescale, seasonal variability in the upper ocean, and the chosen ratio of anthropogenically induced changes in δ13C and DIC. Finally, column inventory changes based on the two methods also exhibit very similar mean Canthuptake rates. The paired DIC concentration and stable isotope dataset may enhance our ability to constrain Canthaccumulation and its controlling mechanisms in the ocean. 
    more » « less
  5. ; ; ; ; ; ; ; ; (, Geophysical Research Letters)
    Abstract This study assesses the effective climate sensitivity (EffCS) and transient climate response (TCR) derived from global energy budget constraints within historical simulations of eight CMIP6 global climate models (GCMs). These calculations are enabled by use of the Radiative Forcing Model Intercomparison Project (RFMIP) simulations, which permit accurate quantification of the radiative forcing. Long‐term historical energy budget constraints generally underestimate EffCS from CO2quadrupling and TCR from CO2ramping, owing to changes in radiative feedbacks and changes in ocean heat uptake efficiency. Atmospheric GCMs forced by observed warming patterns produce lower values of EffCS that are more in line with those inferred from observed historical energy budget changes. The differences in the EffCS estimates from historical energy budget constraints of models and observations are traced to discrepancies between modeled and observed historical surface warming patterns. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email