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.


Title: Equatorial Pacific pCO 2 Interannual Variability in CMIP6 Models
Abstract The El Niño‐Southern Oscillation (ENSO) in the equatorial Pacific is the dominant mode of global air‐sea carbon dioxide (CO2) flux interannual variability (IAV). Air‐sea CO2fluxes are driven by the difference between atmospheric and surface ocean pCO2, with variability of the latter driving flux variability. Previous studies found that models in Coupled Model Intercomparison Project Phase 5 (CMIP5) failed to reproduce the observed ENSO‐related pattern of CO2fluxes and had weak pCO2IAV, which were explained by both weak upwelling IAV and weak mean vertical dissolved inorganic carbon (DIC) gradients. We assess whether the latest generation of CMIP6 models can reproduce equatorial Pacific pCO2IAV by validating models against observations‐based data products. We decompose pCO2IAV into thermally and non‐thermally driven anomalies to examine the balance between these competing anomalies, which explain the total pCO2IAV. The majority of CMIP6 models underestimate pCO2IAV, while they overestimate sea surface temperature IAV. Insufficient compensation of non‐thermal pCO2to thermal pCO2IAV in models results in weak total pCO2IAV. We compare the relative strengths of the vertical transport of temperature and DIC and evaluate their contributions to thermal and non‐thermal pCO2anomalies. Model‐to‐observations‐based product comparisons reveal that modeled mean vertical DIC gradients are biased weak relative to their mean vertical temperature gradients, but upwelling acting on these gradients is insufficient to explain the relative magnitudes of thermal and non‐thermal pCO2anomalies.  more » « less
Award ID(s):
2019625 1948624 2219829
PAR ID:
10386413
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Biogeosciences
Volume:
127
Issue:
12
ISSN:
2169-8953
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract This study examines historical simulations of ENSO in the E3SM-1-0, CESM2, and GFDL-CM4 climate models, provided by three leading U.S. modeling centers as part of the Coupled Model Intercomparison Project phase 6 (CMIP6). These new models have made substantial progress in simulating ENSO’s key features, including: amplitude; timescale; spatial patterns; phase-locking; spring persistence barrier; and recharge oscillator dynamics. However, some important features of ENSO are still a challenge to simulate. In the central and eastern equatorial Pacific, the models’ weaker-than-observed subsurface zonal current anomalies and zonal temperature gradient anomalies serve to weaken the nonlinear zonal advection of subsurface temperatures, leading to insufficient warm/cold asymmetry of ENSO’s sea surface temperature anomalies (SSTA). In the western equatorial Pacific, the models’ excessive simulated zonal SST gradients amplify their zonal temperature advection, causing their SSTA to extend farther west than observed. The models underestimate both ENSO’s positive dynamic feedbacks (due to insufficient zonal wind stress responses to SSTA) and its thermodynamic damping (due to insufficient convective cloud shading of eastern Pacific SSTA during warm events); compensation between these biases leads to realistic linear growth rates for ENSO, but for somewhat unrealistic reasons. The models also exhibit stronger-than-observed feedbacks onto eastern equatorial Pacific SSTAs from thermocline depth anomalies, which accelerates the transitions between events and shortens the simulated ENSO period relative to observations. Implications for diagnosing and simulating ENSO in climate models are discussed. 
    more » « less
  2. Abstract El Niño–Southern Oscillation (ENSO) is an important but not the only source of interannual variability over the Indo–western Pacific. Non-ENSO forced variability in the region has received recent attention because of the implications for rainy-season prediction. Using a 35-member CESM1 Large Ensemble (CESM-LE) and 30 CMIP6 models, this study shows that the ensemble means project intensified interannual variability for precipitation, low-level winds, and sea level pressure under global warming, associated with the enhanced large-scale anomalous anticyclone (AAC) over the tropical northwestern (NW) Pacific after the ENSO signal is removed. A decomposition based on the column water vapor budget reveals that enhanced precipitation variability is due to the increased background specific humidity. The resultant anomalous diabatic heating intensifies the AAC, which further strengthens the precipitation anomalies. Over the tropical NW Pacific, the wind-induced evaporative cooling on the southeastern flank of the AAC is countered by the increased shortwave radiation due to the strengthened precipitation reduction. Tropospheric temperature anomalies in the ensemble means show no significant change, suggesting no apparent change of the interbasin positive feedback between the AAC and northern Indian Ocean SST. Intermodel analysis based on CMIP6 reveals that models with a larger increase in ENSO-unrelated precipitation variability over the NW Pacific are associated with stronger background warming in the eastern equatorial Pacific, due to the modulated Walker and Hadley circulations. 
    more » « less
  3. Abstract The deep ocean releases large amounts of old, pre‐industrial carbon dioxide (CO2) to the atmosphere through upwelling in the Southern Ocean, which counters the marine carbon uptake occurring elsewhere. This Southern Ocean CO2release is relevant to the global climate because its changes could alter atmospheric CO2levels on long time scales, and also affects the present‐day potential of the Southern Ocean to take up anthropogenic CO2. Here, year‐round profiling float measurements show that this CO2release arises from a zonal band of upwelling waters between the Subantarctic Front and wintertime sea‐ice edge. This band of high CO2subsurface water coincides with the outcropping of the 27.8 kg m−3isoneutral density surface that characterizes Indo‐Pacific Deep Water (IPDW). It has a potential partial pressure of CO2exceeding current atmospheric CO2levels (∆PCO2) by 175 ± 32 μatm. Ship‐based measurements reveal that IPDW exhibits a distinct ∆PCO2maximum in the ocean, which is set by remineralization of organic carbon and originates from the northern Pacific and Indian Ocean basins. Below this IPDW layer, the carbon content increases downwards, whereas ∆PCO2decreases. Most of this vertical ∆PCO2decline results from decreasing temperatures and increasing alkalinity due to an increased fraction of calcium carbonate dissolution. These two factors limit the CO2outgassing from the high‐carbon content deep waters on more southerly surface outcrops. Our results imply that the response of Southern Ocean CO2fluxes to possible future changes in upwelling are sensitive to the subsurface carbon chemistry set by the vertical remineralization and dissolution profiles. 
    more » « less
  4. Abstract The global air‐sea CO2flux (F) impacts and is impacted by a plethora of climate‐related processes operating at multiple time scales. In bulk mass transfer formulations, F is driven by physico‐ and bio‐chemical factors such as the air‐sea partial pressure difference (∆pCO2), gas transfer velocity, sea surface temperature, and salinity–all varying at multiple time scales. To de‐convolve the impact of these factors on variability in F at different time scales, time‐resolved estimates of F were computed using a global data set assembled between 1988 and 2015. The F anomalies were defined as temporal deviations from the 28‐year time‐averaged value. Spectral analysis revealed four dominant timescales of variability in F–subseasonal, seasonal, interannual, and decadal with relative amplitude differences varying across regions. A second‐order Taylor series expansion was then conducted along these four timescales to separate drivers across differing regions. The analysis showed that on subseasonal timescales, wind speed variability explains some 66% of the global F anomaly and is the dominant driver. On seasonal, interannual, and decadal timescales, the ∆pCO2effect controlled by the ∆pCO2anomaly, explained much of the F anomaly. On decadal timescales, the F anomaly was almost entirely governed by the ∆pCO2effect with large contributions from high latitudes. The main drivers across timescales also dominate the regional F anomaly, particularly in the mid‐high latitude regions. Finally, the driver of the ∆pCO2effect was closely connected with the relative strength of atmospheric pCO2and the nonthermal component of oceanic pCO2anomaly associated with dissolved inorganic carbon and alkalinity. 
    more » « less
  5. Abstract The Chukchi Sea is an increasing CO2sink driven by rapid climate changes. Understanding the seasonal variation of air‐sea CO2exchange and the underlying mechanisms of biogeochemical dynamics is important for predicting impacts of climate change on and feedbacks by the ocean. Here, we present a unique data set of underway sea surface partial pressure of CO2(pCO2) and discrete samples of biogeochemical properties collected in five consecutive cruises in 2014 and examine the seasonal variations in air‐sea CO2flux and net community production (NCP). We found that thermal and non‐thermal effects have different impacts on sea surfacepCO2and thus the air‐sea CO2flux in different water masses. The Bering summer water combined with meltwater has a significantly greater atmospheric CO2uptake potential than that of the Alaskan Coastal Water in the southern Chukchi Sea in summer, due to stronger biological CO2removal and a weaker thermal effect. By analyzing the seasonal drawdown of dissolved inorganic carbon (DIC) and nutrients, we found that DIC‐based NCP was higher than nitrate‐based NCP by 66%–84% and attributable to partially decoupled C and N uptake because of a variable phytoplankton stoichiometry. A box model with a non‐Redfield C:N uptake ratio can adequately reproduce observedpCO2and DIC, which reveals that, during the intensive growing season (late spring to early summer), 30%–46% CO2uptake in the Chukchi Sea was supported by a flexible stoichiometry of phytoplankton. These findings have important ramification for forecasting the responses of CO2uptake of the Chukchi ecosystem to climate change. 
    more » « less