scholarly journals Variability of North Atlantic CO<sub>2</sub> fluxes for the 2000–2017 period estimated from atmospheric inverse analyses

2021 ◽  
Vol 18 (15) ◽  
pp. 4549-4570
Author(s):  
Zhaohui Chen ◽  
Parvadha Suntharalingam ◽  
Andrew J. Watson ◽  
Ute Schuster ◽  
Jiang Zhu ◽  
...  

Abstract. We present new estimates of the regional North Atlantic (15–80∘ N) CO2 flux for the 2000–2017 period using atmospheric CO2 measurements from the NOAA long-term surface site network in combination with an atmospheric carbon cycle data assimilation system (GEOS-Chem–LETKF, Local Ensemble Transform Kalman Filter). We assess the sensitivity of flux estimates to alternative ocean CO2 prior flux distributions and to the specification of uncertainties associated with ocean fluxes. We present a new scheme to characterize uncertainty in ocean prior fluxes, derived from a set of eight surface pCO2-based ocean flux products, and which reflects uncertainties associated with measurement density and pCO2-interpolation methods. This scheme provides improved model performance in comparison to fixed prior uncertainty schemes, based on metrics of model–observation differences at the network of surface sites. Long-term average posterior flux estimates for the 2000–2017 period from our GEOS-Chem–LETKF analyses are −0.255 ± 0.037 PgC yr−1 for the subtropical basin (15–50∘ N) and −0.203 ± 0.037 PgC yr−1 for the subpolar region (50–80∘ N, eastern boundary at 20∘ E). Our basin-scale estimates of interannual variability (IAV) are 0.036 ± 0.006 and 0.034 ± 0.009 PgC yr−1 for subtropical and subpolar regions, respectively. We find statistically significant trends in carbon uptake for the subtropical and subpolar North Atlantic of −0.064 ± 0.007 and −0.063 ± 0.008 PgC yr−1 decade−1; these trends are of comparable magnitude to estimates from surface ocean pCO2-based flux products, but they are larger, by a factor of 3–4, than trends estimated from global ocean biogeochemistry models.

2020 ◽  
Author(s):  
Zhaohui Chen ◽  
Parvadha Suntharalingam ◽  
Andrew J. Watson ◽  
Ute Schuster ◽  
Jiang Zhu ◽  
...  

Abstract. We present new estimates of the regional North Atlantic (15° N–80° N) CO2 flux for the 2000–2017 period using atmospheric CO2 measurements from the NOAA long term surface site network in combination with an atmospheric data assimilation system (GEOSChem–LETKF). We also assess the sensitivity of flux estimates to the representation of the prior ocean flux distribution and to the associated specification of prior flux uncertainty, including a specification that is dependent on the agreement among the multiple representations of the prior ocean flux. Long term average flux estimates for the 2000–2017 period are −0.26±0.04 PgC y−1 for the subtropical basin (15° N–50° N), and −0.25±0.04 PgC y−1 for the subpolar region (50° N–80° N, west of 20° E). Our basin–scale estimates of the amplitude of interannual variability (IAV) are 0.037±0.006 PgC y−1 and 0.025±0.009 PgC y−1 for subtropical and subpolar regions respectively. We find a statistically significant trend in carbon uptake for the subtropical North Atlantic of −0.062±0.009 PgC y−1 decade−1 over this period.


2005 ◽  
Vol 2 (4) ◽  
pp. 437-472 ◽  
Author(s):  
S. Raynaud ◽  
O. Aumont ◽  
K. B. Rodgers ◽  
P. Yiou ◽  
J. C. Orr

Abstract. The magnitude of the interannual variability of North Atlantic air-sea CO2 fluxes remains uncertain. Fluxes inferred from atmospheric inversions have large variability, whereas those simulated by ocean models have small variability. Part of the difference is that unlike typical atmospheric inversions, ocean models come with spatial resolution at the sub-basin scale. Here we explore sub-basin-scale spatiotemporal variability in the North Atlantic in one ocean model in order to better understand why the the North Atlantic basin may well contribute very little to the global variability of air-sea CO2 flux. We made two simulations with a biogeochemical model coupled to a global ocean general circulation model (OGCM), which itself was forced by 55-year NCEP reanalysis fields. In the first simulation, atmospheric CO2 was maintained at the preindustrial level (278 ppmv); in the second simulation, atmospheric CO2 followed the observed increase. Simulated air-sea CO2 fluxes and associated variables were analysed with a statistical tool known as multichannel singular spectrum analysis (MSSA). We found that the subtropical gyre is not the largest contributor to the overall, basin-wide variability, in contrast to previous suggestions. The subpolar gyre and the inter-gyre region (the transition area between subpolar and subtropical gyres) also contribute with multipolar anomalies at multiple frequencies: these tend to cancel one another in terms of the basin-wide air-sea CO2 flux. We found a strong correlation between the air-sea CO2 fluxes and the North Atlantic Oscillation (NAO), but only if one takes into account time lags as does MSSA (maximum r=0.64 for lags between 1 and 3 years). The contribution of anthropogenic CO2 to total variability was negligible at interannual time scales, whereas at the decadal (13-year) time scale, it increased variability by 30%.


2021 ◽  
Author(s):  
Geneviève Elsworth ◽  
Nicole Lovenduski ◽  
Karen McKinnon

&lt;p&gt;Internal climate variability plays an important role in the abundance and distribution of phytoplankton in the global ocean. Previous studies using large ensembles of Earth system models (ESMs) have demonstrated their utility in the study of marine phytoplankton variability. These ESM large ensembles simulate the evolution of multiple alternate realities, each with a different phasing of internal climate variability. However, ESMs may not accurately represent real world variability as recorded via satellite and in situ observations of ocean chlorophyll over the past few decades. Observational records of surface ocean chlorophyll equate to a single ensemble member in the large ensemble framework, and this can cloud the interpretation of long-term trends: are they externally forced, caused by the phasing of internal variability, or both? Here, we use a novel statistical emulation technique to place the observational record of surface ocean chlorophyll into the large ensemble framework. Much like a large initial condition ensemble generated with an ESM, the resulting synthetic ensemble represents multiple possible evolutions of ocean chlorophyll concentration, each with a different phasing of internal climate variability. We further demonstrate the validity of our statistical approach by recreating a ESM ensemble of chlorophyll using only a single ESM ensemble member. We use the synthetic ensemble to explore the interpretation of long-term trends in the presence of internal variability. Our results suggest the potential to explore this approach for other ocean biogeochemical variables.&lt;/p&gt;


2021 ◽  
Vol 13 (2) ◽  
pp. 777-808
Author(s):  
Luke Gregor ◽  
Nicolas Gruber

Abstract. Ocean acidification has profoundly altered the ocean's carbonate chemistry since preindustrial times, with potentially serious consequences for marine life. Yet, no long-term, global observation-based data set exists that allows us to study changes in ocean acidification for all carbonate system parameters over the last few decades. Here, we fill this gap and present a methodologically consistent global data set of all relevant surface ocean parameters, i.e., dissolved inorganic carbon (DIC), total alkalinity (TA), partial pressure of CO2 (pCO2), pH, and the saturation state with respect to mineral CaCO3 (Ω) at a monthly resolution over the period 1985 through 2018 at a spatial resolution of 1∘×1∘. This data set, named OceanSODA-ETHZ, was created by extrapolating in time and space the surface ocean observations of pCO2 (from the Surface Ocean CO2 Atlas, SOCAT) and total alkalinity (TA; from the Global Ocean Data Analysis Project, GLODAP) using the newly developed Geospatial Random Cluster Ensemble Regression (GRaCER) method (code available at https://doi.org/10.5281/zenodo.4455354, Gregor, 2021). This method is based on a two-step (cluster-regression) approach but extends it by considering an ensemble of such cluster regressions, leading to improved robustness. Surface ocean DIC, pH, and Ω were then computed from the globally mapped pCO2 and TA using the thermodynamic equations of the carbonate system. For the open ocean, the cluster-regression method estimates pCO2 and TA with global near-zero biases and root mean squared errors of 12 µatm and 13 µmol kg−1, respectively. Taking into account also the measurement and representation errors, the total uncertainty increases to 14 µatm and 21 µmol kg−1, respectively. We assess the fidelity of the computed parameters by comparing them to direct observations from GLODAP, finding surface ocean pH and DIC global biases of near zero, as well as root mean squared errors of 0.023 and 16 µmol kg−1, respectively. These uncertainties are very comparable to those expected by propagating the total uncertainty from pCO2 and TA through the thermodynamic computations, indicating a robust and conservative assessment of the uncertainties. We illustrate the potential of this new data set by analyzing the climatological mean seasonal cycles of the different parameters of the surface ocean carbonate system, highlighting their commonalities and differences. Further, this data set provides a novel constraint on the global- and basin-scale trends in ocean acidification for all parameters. Concretely, we find for the period 1990 through 2018 global mean trends of 8.6 ± 0.1 µmol kg−1 per decade for DIC, −0.016 ± 0.000 per decade for pH, 16.5 ± 0.1 µatm per decade for pCO2, and −0.07 ± 0.00 per decade for Ω. The OceanSODA-ETHZ data can be downloaded from https://doi.org/10.25921/m5wx-ja34 (Gregor and Gruber, 2020).


2021 ◽  
Vol 17 (5) ◽  
pp. 2223-2254
Author(s):  
Katherine A. Crichton ◽  
Andy Ridgwell ◽  
Daniel J. Lunt ◽  
Alex Farnsworth ◽  
Paul N. Pearson

Abstract. Since the middle Miocene (15 Ma, million years ago), the Earth's climate has undergone a long-term cooling trend, characterised by a reduction in ocean temperatures of up to 7–8 ∘C. The causes of this cooling are primarily thought to be due to tectonic plate movements driving changes in large-scale ocean circulation patterns, and hence heat redistribution, in conjunction with a drop in atmospheric greenhouse gas forcing (and attendant ice-sheet growth and feedback). In this study, we assess the potential to constrain the evolving patterns of global ocean circulation and cooling over the last 15 Ma by assimilating a variety of marine sediment proxy data in an Earth system model. We do this by first compiling surface and benthic ocean temperature and benthic carbon-13 (δ13C) data in a series of seven time slices spaced at approximately 2.5 Myr intervals. We then pair this with a corresponding series of tectonic and climate boundary condition reconstructions in the cGENIE (“muffin” release) Earth system model, including alternative possibilities for an open vs. closed Central American Seaway (CAS) from 10 Ma onwards. In the cGENIE model, we explore uncertainty in greenhouse gas forcing and the magnitude of North Pacific to North Atlantic salinity flux adjustment required in the model to create an Atlantic Meridional Overturning Circulation (AMOC) of a specific strength, via a series of 12 (one for each tectonic reconstruction) 2D parameter ensembles. Each ensemble member is then tested against the observed global temperature and benthic δ13C patterns. We identify that a relatively high CO2 equivalent forcing of 1120 ppm is required at 15 Ma in cGENIE to reproduce proxy temperature estimates in the model, noting that this CO2 forcing is dependent on the cGENIE model's climate sensitivity and that it incorporates the effects of all greenhouse gases. We find that reproducing the observed long-term cooling trend requires a progressively declining greenhouse gas forcing in the model. In parallel to this, the strength of the AMOC increases with time despite a reduction in the salinity of the surface North Atlantic over the cooling period, attributable to falling intensity of the hydrological cycle and to lowering polar temperatures, both caused by CO2-driven global cooling. We also find that a closed CAS from 10 Ma to present shows better agreement between benthic δ13C patterns and our particular series of model configurations and data. A final outcome of our analysis is a pronounced ca. 1.5 ‰ decline occurring in atmospheric (and ca. 1 ‰ ocean surface) δ13C that could be used to inform future δ13C-based proxy reconstructions.


2019 ◽  
Author(s):  
Fabrice Lacroix ◽  
Tatiana Ilyina ◽  
Jens Hartmann

Abstract. Rivers are a major source of nutrients, carbon and alkalinity for the global ocean, where the delivered compounds strongly impact biogeochemical processes. In this study, we firstly estimate pre-industrial riverine fluxes of nutrients, carbon and alkalinity based on a hierarchy of weathering and land-ocean export models, while identifying regional hotspots of the land-ocean exports. Secondly, we implement the riverine loads into a global biogeochemical ocean model and describe their implications for oceanic nutrient concentrations, the net primary production (NPP) and CO2 fluxes globally, as well as in a regional shelf analysis. Thirdly, we quantify the terrestrial origins and the long-term oceanic fate of riverine carbon in the framework, while assessing the potential implementation of riverine carbon fluxes in a fully coupled land-atmosphere-ocean model. Our approach leads to annual pre-industrial riverine exports of 3.7 Tg P, 27 Tg N, 158 Tg Si and 603 Tg C, which were derived from weathering and non-weathering sources and were fractionated into organic and inorganic compounds. We thereby identify the tropical Atlantic catchments (20 % of global C), Arctic rivers (9 % of total C) and Southeast Asian rivers (15 % of total C) as dominant providers of carbon to the ocean. The riverine exports lead to a global oceanic source of CO2 to the atmosphere (231 Tg C yr−1), which is largely a result of a source from inorganic riverine carbon loads (183 Tg C yr−1), and from organic riverine carbon inputs (128 Tg C yr−1). Additionally, a sink of 80 Tg C yr−1 is caused by the enhancement of the biological carbon uptake by dissolved inorganic nutrient inputs, resulting alkalinity production and a slight model drift. While large outgassing fluxes are mostly found in proximity to major river mouths, substantial outgassing fluxes can also be observed further offshore, most prominently in the tropical Atlantic. Furthermore, we find evidence for the interhemispheric transfer of carbon in the model; we detect a stronger relative outgassing flux (49 % of global river induced outgassing) in the southern hemisphere in comparison to the hemisphere's relative riverine inputs (33 % of global river inputs), as well as an outgassing flux of 17 Tg C yr-1 in the Southern Ocean. Riverine exports lead to a strong increase in NPP in the tropical West Atlantic, Bay of Bengal and the East China Sea (166 %, 377 % and 71 % respectively). While the NPP is not strongly sensitive to riverine loads on the light limited Arctic shelves, the CO2 flux is strongly altered due to substantial dissolved carbon supplies to the region. While our study confirms that the ocean circulation is the main driver for open ocean biogeochemical distributions, it reveals the necessity to consider riverine exports for the representation of heterogeneous features of the coastal ocean, to represent riverine-induced carbon outgassing, as well as to consider the long-term volcanic CO2 flux to close the atmospheric carbon budget in a coupled land-ocean-atmosphere setting.


2021 ◽  
Author(s):  
Mario Uchimiya ◽  
William Schroer ◽  
Malin Olofsson ◽  
Arthur S. Edison ◽  
Mary Ann Moran

AbstractOrganic carbon transfer between surface ocean photosynthetic and heterotrophic microbes is a central but poorly understood process in the global carbon cycle. In a model community in which diatom extracellular release of organic molecules sustained growth of a co-cultured bacterium, we determined quantitative changes in the diatom endometabolome and the bacterial uptake transcriptome over two diel cycles. Of the nuclear magnetic resonance (NMR) peaks in the diatom endometabolites, 38% had diel patterns with noon or mid-afternoon maxima; the remaining either increased (36%) or decreased (26%) through time. Of the genes in the bacterial uptake transcriptome, 94% had a diel pattern with a noon maximum; the remaining decreased over time (6%). Eight diatom endometabolites identified with high confidence were matched to the bacterial genes mediating their utilization. Modeling of these coupled inventories with only diffusion-based phytoplankton extracellular release could not reproduce all the patterns. Addition of active release mechanisms for physiological balance and bacterial recognition significantly improved model performance. Estimates of phytoplankton extracellular release range from only a few percent to nearly half of annual net primary production. Improved understanding of the factors that influence metabolite release and consumption by surface ocean microbes will better constrain this globally significant carbon flux.


2008 ◽  
Vol 21 (24) ◽  
pp. 6599-6615 ◽  
Author(s):  
Arne Biastoch ◽  
Claus W. Böning ◽  
Julia Getzlaff ◽  
Jean-Marc Molines ◽  
Gurvan Madec

Abstract The causes and characteristics of interannual–decadal variability of the meridional overturning circulation (MOC) in the North Atlantic are investigated with a suite of basin-scale ocean models [the Family of Linked Atlantic Model Experiments (FLAME)] and global ocean–ice models (ORCA), varying in resolution from medium to eddy resolving (½°–1/12°), using various forcing configurations built on bulk formulations invoking atmospheric reanalysis products. Comparison of the model hindcasts indicates similar MOC variability characteristics on time scales up to a decade; both model architectures also simulate an upward trend in MOC strength between the early 1970s and mid-1990s. The causes of the MOC changes are examined by perturbation experiments aimed selectively at the response to individual forcing components. The solutions emphasize an inherently linear character of the midlatitude MOC variability by demonstrating that the anomalies of a (non–eddy resolving) hindcast simulation can be understood as a superposition of decadal and longer-term signals originating from thermohaline forcing variability, and a higher-frequency wind-driven variability. The thermohaline MOC signal is linked to the variability in subarctic deep-water formation, and rapidly progressing to the tropical Atlantic. However, throughout the subtropical and midlatitude North Atlantic, this signal is effectively masked by stronger MOC variability related to wind forcing and, especially north of 30°–35°N, by internally induced (eddy) fluctuations.


2013 ◽  
Vol 10 (11) ◽  
pp. 7793-7815 ◽  
Author(s):  
P. Landschützer ◽  
N. Gruber ◽  
D. C. E. Bakker ◽  
U. Schuster ◽  
S. Nakaoka ◽  
...  

Abstract. The Atlantic Ocean is one of the most important sinks for atmospheric carbon dioxide (CO2), but this sink has been shown to vary substantially in time. Here we use surface ocean CO2 observations to estimate this sink and the temporal variability from 1998 through 2007 in the Atlantic Ocean. We benefit from (i) a continuous improvement of the observations, i.e. the Surface Ocean CO2 Atlas (SOCAT) v1.5 database and (ii) a newly developed technique to interpolate the observations in space and time. In particular, we use a two-step neural network approach to reconstruct basin-wide monthly maps of the sea surface partial pressure of CO2 (pCO2) at a resolution of 1° × 1°. From those, we compute the air–sea CO2 flux maps using a standard gas exchange parameterization and high-resolution wind speeds. The neural networks fit the observed pCO2 data with a root mean square error (RMSE) of about 10 μatm and with almost no bias. A check against independent time-series data and new data from SOCAT v2 reveals a larger RMSE of 22.8 μatm for the entire Atlantic Ocean, which decreases to 16.3 μatm for data south of 40° N. We estimate a decadal mean uptake flux of −0.45 ± 0.15 Pg C yr−1 for the Atlantic between 44° S and 79° N, representing the sum of a strong uptake north of 18° N (−0.39 ± 0.10 Pg C yr−1), outgassing in the tropics (18° S–18° N, 0.11 ± 0.07 Pg C yr−1), and uptake in the subtropical/temperate South Atlantic south of 18° S (−0.16 ± 0.06 Pg C yr−1), consistent with recent studies. The strongest seasonal variability of the CO2 flux occurs in the temperature-driven subtropical North Atlantic, with uptake in winter and outgassing in summer. The seasonal cycle is antiphased in the subpolar latitudes relative to the subtropics largely as a result of the biologically driven winter-to-summer drawdown of CO2. Over the 10 yr analysis period (1998 through 2007), sea surface pCO2 increased faster than that of the atmosphere in large areas poleward of 40° N, while in other regions of the North Atlantic the sea surface pCO2 increased at a slower rate, resulting in a barely changing Atlantic carbon sink north of the Equator (−0.01 ± 0.02 Pg C yr−1 decade−1). Surface ocean pCO2 increased at a slower rate relative to atmospheric CO2 over most of the Atlantic south of the Equator, leading to a substantial trend toward a stronger CO2 sink for the entire South Atlantic (−0.14 ± 0.02 Pg C yr−1 decade−1). In contrast to the 10 yr trends, the Atlantic Ocean carbon sink varies relatively little on inter-annual timescales (±0.04 Pg C yr−1; 1 σ).


2017 ◽  
Vol 30 (6) ◽  
pp. 1985-1997 ◽  
Author(s):  
Damien Desbruyères ◽  
Elaine L. McDonagh ◽  
Brian A. King ◽  
Virginie Thierry

The early twenty-first century’s warming trend of the full-depth global ocean is calculated by combining the analysis of Argo (top 2000 m) and repeat hydrography into a blended full-depth observing system. The surface-to-bottom temperature change over the last decade of sustained observation is equivalent to a heat uptake of 0.71 ± 0.09 W m−2 applied over the surface of Earth, 90% of it being found above 2000-m depth. The authors decompose the temperature trend pointwise into changes in isopycnal depth (heave) and temperature changes along an isopycnal (spiciness) to describe the mechanisms controlling the variability. The heave component dominates the global heat content increase, with the largest trends found in the Southern Hemisphere’s extratropics (0–2000 m) highlighting a volumetric increase of subtropical mode waters. Significant heave-related warming is also found in the deep North Atlantic and Southern Oceans (2000–4000 m), reflecting a potential decrease in deep water mass renewal rates. The spiciness component shows its strongest contribution at intermediate levels (700–2000 m), with striking localized warming signals in regions of intense vertical mixing (North Atlantic and Southern Oceans). Finally, the agreement between the independent Argo and repeat hydrography temperature changes at 2000 m provides an overall good confidence in the blended heat content evaluation on global and ocean scales but also highlights basin-scale discrepancies between the two independent estimates. Those mismatches are largest in those basins with the largest heave signature (Southern Ocean) and reflect both the temporal and spatial sparseness of the hydrography sampling.


Sign in / Sign up

Export Citation Format

Share Document