Mechanistic Drivers of Reemergence of Anthropogenic Carbon in the Equatorial Pacific

AbstractRelatively rapid reemergence of anthropogenic carbon (Cant) in the Equatorial Pacific is of potential importance for its impact on the carbonate buffering capacity of surface seawater and thereby impeding the ocean's ability to further absorb Cant from the atmosphere. We explore the mechanisms sustaining Cant reemergence (upwelling) from the thermocline to surface layers by applying water mass transformation diagnostics to a global ocean/sea ice/biogeochemistry model. We find that the upwelling rate of Cant (0.4 PgC yr−1) from the thermocline to the surface layer is almost twice as large as air‐sea Cant fluxes (0.203 PgC yr−1). The upwelling of Cant from the thermocline to the surface layer can be understood as a two‐step process: The first being due to diapycnal diffusive transformation fluxes and the second due to surface buoyancy fluxes. We also find that this reemergence of Cant decreases dramatically during the 1982/1983 and 1997/1998 El Niño events

Anthropogenic nitrogen inputs and impacts on oceanic N2O fluxes in the northern Indian Ocean: The need for an integrated observation and modelling approach

Anthropogenically-derived nitrogen input to the northern Indian Ocean has increased significantly in recent decades, based on both observational and model-derived estimates This external nutrient source is supplied by atmospheric deposition and riverine fluxes, and has the potential to affect the vulnerable biogeochemical systems of the Arabian Sea and Bay of Bengal, influencing productivity and oceanic production of the greenhouse-gas nitrous-oxide (N2O). We summarize current estimates of this external nitrogen source to the northern Indian Ocean from observations and models, highlight implications for regional marine N2O emissions using model-based analyses, and make recommendations for measurement and model needs to improve current estimates and future predictions of this impact. Current observationally-derived estimates of deposition and riverine nitrogen inputs are limited by sparse measurements and uncertainties on accurate characterization of nitrogen species composition. Ocean model assessments of the impact of external nitrogen sources on regional marine N2O production in the northern Indian Ocean estimate potentially significant changes but also have large associated uncertainties. We recommend an integrated program of basin-wide measurements combined with high-resolution modeling and more detailed characterization of nitrogen-cycle process to address these uncertainties and improve current estimates and predictions

Concentrations, ratios, and sinking fluxes of major bioelements at Ocean Station Papa

© The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Roca-Marti, M., Benitez-Nelson, C. R., Umhau, B. P., Wyatt, A. M., Clevenger, S. J., Pike, S., Horner, T. J., Estapa, M. L., Resplandy, L., & Buesseler, K. O. Concentrations, ratios, and sinking fluxes of major bioelements at Ocean Station Papa. Elementa: Science of the Anthropocene, 9(1), (2021): 00166, https://doi.org/10.1525/elementa.2020.00166.Fluxes of major bioelements associated with sinking particles were quantified in late summer 2018 as part of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) field campaign near Ocean Station Papa in the subarctic northeast Pacific. The thorium-234 method was used in conjunction with size-fractionated (1–5, 5–51, and >51 μm) concentrations of particulate nitrogen (PN), total particulate phosphorus (TPP), biogenic silica (bSi), and particulate inorganic carbon (PIC) collected using large volume filtration via in situ pumps. We build upon recent work quantifying POC fluxes during EXPORTS. Similar remineralization length scales were observed for both POC and PN across all particle size classes from depths of 50–500 m. Unlike bSi and PIC, the soft tissue–associated POC, PN, and TPP fluxes strongly attenuated from 50 m to the base of the euphotic zone (approximately 120 m). Cruise-average thorium-234-derived fluxes (mmol m–2 d–1) at 120 m were 1.7 ± 0.6 for POC, 0.22 ± 0.07 for PN, 0.019 ± 0.007 for TPP, 0.69 ± 0.26 for bSi, and 0.055 ± 0.022 for PIC. These bioelement fluxes were similar to previous observations at this site, with the exception of PIC, which was 1 to 2 orders of magnitude lower. Transfer efficiencies within the upper twilight zone (flux 220 m/flux 120 m) were highest for PIC (84%) and bSi (79%), followed by POC (61%), PN (58%), and TPP (49%). These differences indicate preferential remineralization of TPP relative to POC or PN and larger losses of soft tissue relative to biominerals in sinking particles below the euphotic zone. Comprehensive characterization of the particulate bioelement fluxes obtained here will support future efforts linking phytoplankton community composition and food-web dynamics to the composition, magnitude, and attenuation of material that sinks to deeper waters.The authors would like to acknowledge support from the National Aeronautics and Space Administration as part of the EXport Processes in the Ocean from RemoTe Sensing program awards 80NSSC17K0555 and 80NSSC17K0662. They also acknowledge the funding from the Woods Hole Oceanographic Institution’s Ocean Twilight Zone study for MRM and KOB, the National Science Foundation Graduate Research Fellowship Program for AMW, and the Ocean Frontier Institute for MRM

Decadal acidification in the water masses of the Atlantic Ocean

Global ocean acidification is caused primarily by the ocean's uptake of CO2 as a consequence of increasing atmospheric CO2 levels. We present observations of the oceanic decrease in pH at the basin scale (50°S-36°N) for the Atlantic Ocean over two decades (1993-2013). Changes in pH associated with the uptake of anthropogenic CO2 (ΔpHCant) and with variations caused by biological activity and ocean circulation (ΔpHNat) are evaluated for different water masses. Output from an Institut Pierre Simon Laplace climate model is used to place the results into a longer-term perspective and to elucidate the mechanisms responsible for pH change. The largest decreases in pH (ΔpH) were observed in central, mode, and intermediate waters, with a maximum ΔpH value in South Atlantic Central Waters of -0.042 ± 0.003. The ΔpH trended toward zero in deep and bottom waters. Observations and model results show that pH changes generally are dominated by the anthropogenic component, which accounts for rates between -0.0015 and -0.0020/y in the central waters. The anthropogenic and natural components are of the same order of magnitude and reinforce one another in mode and intermediate waters over the time period. Large negative ΔpHNat values observed in mode and intermediate waters are driven primarily by changes in CO2 content and are consistent with (i) a poleward shift of the formation region during the positive phase of the Southern Annular Mode in the South Atlantic and (ii) an increase in the rate of the water mass formation in the North Atlantic

High-resolution spatial and temporal measurements of particulate organic carbon flux using thorium-234 in the northeast Pacific Ocean during the EXport processes in the ocean from RemoTe sensing field campaign

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Buesseler, K. O., Benitez-Nelson, C. R., Roca-Marti, M., Wyatt, A. M., Resplandy, L., Clevenger, S. J., Drysdale, J. A., Estapa, M. L., Pike, S., & Umhau, B. P. High-resolution spatial and temporal measurements of particulate organic carbon flux using thorium-234 in the northeast Pacific Ocean during the EXport processes in the ocean from RemoTe sensing field campaign. Elementa: Science of the Anthropocene, 8(1), (2020): 030, https://doi.org/10.1525/elementa.030.The EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) program of National Aeronautics and Space Administration focuses on linking remotely sensed properties from satellites to the mechanisms that control the transfer of carbon from surface waters to depth. Here, the naturally occurring radionuclide thorium-234 was used as a tracer of sinking particle flux. More than 950 234Th measurements were made during August–September 2018 at Ocean Station Papa in the northeast Pacific Ocean. High-resolution vertical sampling enabled observations of the spatial and temporal evolution of particle flux in Lagrangian fashion. Thorium-234 profiles were remarkably consistent, with steady-state (SS) 234Th fluxes reaching 1,450 ± 300 dpm m−2 d−1 at 100 m. Nonetheless, 234Th increased by 6%–10% in the upper 60 m during the cruise, leading to consideration of a non-steady-state (NSS) model and/or horizontal transport, with NSS having the largest impact by decreasing SS 234Th fluxes by 30%. Below 100 m, NSS and SS models overlapped. Particulate organic carbon (POC)/234Th ratios decreased with depth in small (1–5 μm) and mid-sized (5–51 μm) particles, while large particle (>51 μm) ratios remained relatively constant, likely influenced by swimmer contamination. Using an average SS and NSS 234Th flux and the POC/234Th ratio of mid-sized particles, we determined a best estimate of POC flux. Maximum POC flux was 5.5 ± 1.7 mmol C m−2 d−1 at 50 m, decreasing by 70% at the base of the primary production zone (117 m). These results support earlier studies that this site is characterized by a modest biological carbon pump, with an export efficiency of 13% ± 5% (POC flux/net primary production at 120 m) and 39% flux attenuation in the subsequent 100 m (POC flux 220 m/POC flux 120m). This work sets the foundation for understanding controls on the biological carbon pump during this EXPORTS campaign.The authors would like to acknowledge support from the National Aeronautics and Space Administration (NASA) as part of the EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) program awards 80NSSC17K0555 and 80NSSC17K0662; the Woods Hole Oceanographic Institution’s Ocean Twilight Zone study for KOB and MRM, and the National Science Foundation Graduate Research Fellowship Program (NSF-GRFP) for funding and support of AW

Consistency and Challenges in the Ocean Carbon Sink Estimate for the Global Carbon Budget

Based on the 2019 assessment of the Global Carbon Project, the ocean took up on average, 2.5 ± 0.6 PgC yr−1 or 23 ± 5% of the total anthropogenic CO2 emissions over the decade 2009–2018. This sink estimate is based on simulation results from global ocean biogeochemical models (GOBMs) and is compared to data-products based on observations of surface ocean pCO2 (partial pressure of CO2) accounting for the outgassing of river-derived CO2. Here we evaluate the GOBM simulations by comparing the simulated surface ocean pCO2 to observations. Based on this comparison, the simulations are well-suited for quantifying the global ocean carbon sink on the time-scale of the annual mean and its multi-decadal trend (RMSE <20 μatm), as well as on the time-scale of multi-year variability (RMSE <10 μatm), despite the large model-data mismatch on the seasonal time-scale (RMSE of 20–80 μatm). Biases in GOBMs have a small effect on the global mean ocean sink (0.05 PgC yr−1), but need to be addressed to improve the regional budgets and model-data comparison. Accounting for non-mapped areas in the data-products reduces their spread as measured by the standard deviation by a third. There is growing evidence and consistency among methods with regard to the patterns of the multi-year variability of the ocean carbon sink, with a global stagnation in the 1990s and an extra-tropical strengthening in the 2000s. GOBMs and data-products point consistently to a shift from a tropical CO2 source to a CO2 sink in recent years. On average, the GOBMs reveal less variations in the sink than the data-based products. Despite the reasonable simulation of surface ocean pCO2 by the GOBMs, there are discrepancies between the resulting sink estimate from GOBMs and data-products. These discrepancies are within the uncertainty of the river flux adjustment, increase over time, and largely stem from the Southern Ocean. Progress in our understanding of the global ocean carbon sink necessitates significant advancement in modeling and observing the Southern Ocean carbon sink including (i) a game-changing increase in high-quality pCO2 observations, and (ii) a critical re-evaluation of the regional river flux adjustment

Inconsistent strategies to spin up models in CMIP5: implications for ocean biogeochemical model performance assessment

International audienceDuring the fifth phase of the Coupled Model Inter-comparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hind-casts were compared with available modern biogeochemi-cal observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to-model uncertainty

Global carbon budget 2019

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) are based on energy statistics and cement production data, while emissions from land use change (ELUC), mainly deforestation, are based on land use and land use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2009–2018), EFF was 9.5±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.9±0.02 GtC yr−1 (2.3±0.01 ppm yr−1), SOCEAN 2.5±0.6 GtC yr−1, and SLAND 3.2±0.6 GtC yr−1, with a budget imbalance BIM of 0.4 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2018 alone, the growth in EFF was about 2.1 % and fossil emissions increased to 10.0±0.5 GtC yr−1, reaching 10 GtC yr−1 for the first time in history, ELUC was 1.5±0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5±0.9 GtC yr−1 (42.5±3.3 GtCO2). Also for 2018, GATM was 5.1±0.2 GtC yr−1 (2.4±0.1 ppm yr−1), SOCEAN was 2.6±0.6 GtC yr−1, and SLAND was 3.5±0.7 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in EFF of +0.6 % (range of −0.2 % to 1.5 %) based on national emissions projections for China, the USA, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. Overall, the mean and trend in the five components of the global carbon budget are consistently estimated over the period 1959–2018, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations shows (1) no consensus in the mean and trend in land use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018a, b, 2016, 2015a, b, 2014, 2013). The data generated by this work are available at https://doi.org/10.18160/gcp-2019 (Friedlingstein et al., 2019)

Global Carbon Budget 2022

Accurate assessment of anthropogenic carbon dioxide (CO$_2$) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO$_2$ emissions (E$_{FOS}$) are based on energy statistics and cement production data, while emissions from land-use change (E$_{LUC}$), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO$_2$ concentration is measured directly, and its growth rate (G$_{ATM}$) is computed from the annual changes in concentration. The ocean CO$_2$ sink (S$_{OCEAN}$) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO$_2$ sink (S$_{LAND}$) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (B$_{IM}$), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2021, E$_{FOS}$ increased by 5.1 % relative to 2020, with fossil emissions at 10.1 ± 0.5 GtC yr$^{−1}$ (9.9 ± 0.5 GtC yr$^{−1}$ when the cement carbonation sink is included), and E$_{LUC}$ was 1.1 ± 0.7 GtC yr$^{−1}$, for a total anthropogenic CO$_2$ emission (including the cement carbonation sink) of 10.9 ± 0.8 GtC yr$^{−1}$ (40.0 ± 2.9 GtCO$_2$). Also, for 2021, G$_{ATM}$ was 5.2 ± 0.2 GtC yr$^{−1}$ (2.5 ± 0.1 ppm yr$^{−1}$), S$_{OCEAN}$ was 2.9  ± 0.4 GtC yr$^{−1}$, and S$_{LAND}$ was 3.5 ± 0.9 GtC yr$^{−1}$, with a B$_{IM}$ of −0.6 GtC yr$^{−1}$ (i.e. the total estimated sources were too low or sinks were too high). The global atmospheric CO$_2$ concentration averaged over 2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest an increase in E$_{FOS}$ relative to 2021 of +1.0 % (0.1 % to 1.9 %) globally and atmospheric CO$_2$ concentration reaching 417.2 ppm, more than 50 % above pre-industrial levels (around 278 ppm). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2021, but discrepancies of up to 1 GtC yr$^{−1}$ persist for the representation of annual to semi-decadal variability in CO$_2$ fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO$_2$ flux in the northern extratropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set. The data presented in this work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b)