Abstract. The uptake of carbonyl sulfide (COS) by terrestrial plants is linked tophotosynthetic uptake of CO2 as these gases partly share the sameuptake pathway. Applying COS as a photosynthesis tracer in models requires anaccurate representation of biosphere COS fluxes, but these models have notbeen extensively evaluated against field observations of COS fluxes. In thispaper, the COS flux as simulated by the Simple Biosphere Model, version 4(SiB4), is updated with the latest mechanistic insights and evaluated with siteobservations from different biomes: one evergreen needleleaf forest, twodeciduous broadleaf forests, three grasslands, and two crop fields spread overEurope and North America. We improved SiB4 in several ways to improve itsrepresentation of COS. To account for the effect of atmospheric COS molefractions on COS biosphere uptake, we replaced the fixed atmospheric COS molefraction boundary condition originally used in SiB4 with spatially andtemporally varying COS mole fraction fields. Seasonal amplitudes of COS molefractions are ∼50–200 ppt at the investigated sites with aminimum mole fraction in the late growing season. Incorporating seasonalvariability into the model reduces COS uptake rates in the late growingseason, allowing better agreement with observations. We also replaced theempirical soil COS uptake model in SiB4 with a mechanistic model thatrepresents both uptake and production of COS in soils, which improves thematch with observations over agricultural fields and fertilized grasslandsoils. The improved version of SiB4 was capable of simulating the diurnal andseasonal variation in COS fluxes in the boreal, temperate, and Mediterraneanregion. Nonetheless, the daytime vegetation COS flux is underestimated onaverage by 8±27 %, albeit with large variability across sites. On aglobal scale, our model modifications decreased the modeled COS terrestrialbiosphere sink from 922 Gg S yr−1 in the original SiB4 to753 Gg S yr−1 in the updated version. The largest decrease influxes was driven by lower atmospheric COS mole fractions over regions withhigh productivity, which highlights the importance of accounting forvariations in atmospheric COS mole fractions. The change to a different soilmodel, on the other hand, had a relatively small effect on the globalbiosphere COS sink. The secondary role of the modeled soil component in theglobal COS budget supports the use of COS as a global photosynthesis tracer. Amore accurate representation of COS uptake in SiB4 should allow for improvedapplication of atmospheric COS as a tracer of local- to global-scaleterrestrial photosynthesis.
more »
« less
Carbonyl sulfide: comparing a mechanistic representation of the vegetation uptake in a land surface model and the leaf relative uptake approach
Abstract. Land surface modellers need measurable proxies toconstrain the quantity of carbon dioxide (CO2) assimilated bycontinental plants through photosynthesis, known as gross primary production(GPP). Carbonyl sulfide (COS), which is taken up by leaves through theirstomates and then hydrolysed by photosynthetic enzymes, is a candidate GPPproxy. A former study with the ORCHIDEE land surface model used a fixedratio of COS uptake to CO2 uptake normalised to respective ambientconcentrations for each vegetation type (leaf relative uptake, LRU) tocompute vegetation COS fluxes from GPP. The LRU approach is known to havelimited accuracy since the LRU ratio changes with variables such asphotosynthetically active radiation (PAR): while CO2 uptake slows underlow light, COS uptake is not light limited. However, the LRU approach hasbeen popular for COS–GPP proxy studies because of its ease of applicationand apparent low contribution to uncertainty for regional-scaleapplications. In this study we refined the COS–GPP relationship andimplemented in ORCHIDEE a mechanistic model that describes COS uptake bycontinental vegetation. We compared the simulated COS fluxes againstmeasured hourly COS fluxes at two sites and studied the model behaviour andlinks with environmental drivers. We performed simulations at a global scale,and we estimated the global COS uptake by vegetation to be −756 Gg S yr−1,in the middle range of former studies (−490 to −1335 Gg S yr−1). Basedon monthly mean fluxes simulated by the mechanistic approach in ORCHIDEE, wederived new LRU values for the different vegetation types, ranging between0.92 and 1.72, close to recently published averages for observed values of1.21 for C4 and 1.68 for C3 plants. We transported the COS using the monthlyvegetation COS fluxes derived from both the mechanistic and the LRUapproaches, and we evaluated the simulated COS concentrations at NOAA sites.Although the mechanistic approach was more appropriate when comparing tohigh-temporal-resolution COS flux measurements, both approaches gave similarresults when transporting with monthly COS fluxes and evaluating COSconcentrations at stations. In our study, uncertainties between these twoapproaches are of secondary importance compared to the uncertainties in theCOS global budget, which are currently a limiting factor to the potential ofCOS concentrations to constrain GPP simulated by land surface models on theglobal scale.
more »
« less
- Award ID(s):
- 1848618
- PAR ID:
- 10311393
- Date Published:
- Journal Name:
- Biogeosciences
- Volume:
- 18
- Issue:
- 9
- ISSN:
- 1726-4189
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. Ground-level ozone (O3) is a major air pollutant that adversely affects human health and ecosystem productivity. Removal of troposphericO3 by plant stomatal uptake can in turn cause damage to plant tissues with ramifications for ecosystem and crop health. In manyatmospheric and land surface models, the functionality of stomata opening is represented by a bulk stomatal conductance, which is oftensemi-empirically parameterized and highly fitted to historical observations. A lack of mechanistic linkage to ecophysiological processes such asphotosynthesis may render models inadequate to represent plant-mediated responses of atmospheric chemistry to long-term changes in CO2,climate, and short-lived air pollutant concentrations. A new ecophysiology module was thus developed to mechanistically simulate land−atmosphereexchange of important gas species in GEOS-Chem, a chemical transport model widely used in atmospheric chemistry studies. The implementation not onlyallows for dry deposition to be coupled with plant ecophysiology but also enables plant and crop productivity and functions to respond dynamically toatmospheric chemical changes. We conduct simulations to evaluate the effects of the ecophysiology module on simulated dry deposition velocity andconcentration of surface O3 against an observation-derived dataset known as SynFlux. Our estimated stomatal conductance and dry depositionvelocity of O3 are close to SynFlux with root-mean-squared errors (RMSEs) below 0.3 cm s−1 across different plant functionaltypes (PFTs), despite an overall positive bias in surface O3 concentration (by up to 16 ppbv). Representing ecophysiology wasfound to reduce the simulated biases in deposition fluxes from the prior model but worsen the positive biases in simulated O3concentrations. The increase in positive concentration biases is mostly attributable to the ecophysiology-based stomatal conductance being generallysmaller (and closer to SynFlux values) than that estimated by the prior semi-empirical formulation, calling for further improvements in non-stomataldepositional and non-depositional processes relevant for O3 simulations. The estimated global O3 deposition flux is864 Tg O3 yr−1 with GEOS-Chem, and the new module decreases this estimate by 92 Tg O3 yr−1. Estimated global grossprimary production (GPP) without O3 damage is 119 Pg C yr−1. O3-induced reduction in GPP is 4.2 Pg C yr−1(3.5 %). An elevated CO2 scenario (580 ppm) yields higher global GPP (+16.8 %) and lower global O3depositional sink (−3.3 %). Global isoprene emission simulated with a photosynthesis-based scheme is 317.9 Tg C yr−1, which is31.2 Tg C yr−1 (−8.9 %) less than that calculated using the MEGAN(Model of Emissions of Gases and Aerosols from Nature) emission algorithm. This new model development dynamicallyrepresents the two-way interactions between vegetation and air pollutants and thus provides a unique capability in evaluating vegetation-mediatedprocesses and feedbacks that can shape atmospheric chemistry and air quality, as well as pollutant impacts on vegetation health, especially for anytimescales shorter than the multidecadal timescale.more » « less
-
null (Ed.)Robust estimates for the rates and trends in terrestrial gross primary production (GPP; plant CO 2 uptake) are needed. Carbonyl sulfide (COS) is the major long-lived sulfur-bearing gas in the atmosphere and a promising proxy for GPP. Large uncertainties in estimating the relative magnitude of the COS sources and sinks limit this approach. Sulfur isotope measurements ( 34 S/ 32 S; δ 34 S) have been suggested as a useful tool to constrain COS sources. Yet such measurements are currently scarce for the atmosphere and absent for the marine source and the plant sink, which are two main fluxes. Here we present sulfur isotopes measurements of marine and atmospheric COS, and of plant-uptake fractionation experiments. These measurements resulted in a complete data-based tropospheric COS isotopic mass balance, which allows improved partition of the sources. We found an isotopic (δ 34 S ± SE) value of 13.9 ± 0.1‰ for the troposphere, with an isotopic seasonal cycle driven by plant uptake. This seasonality agrees with a fractionation of −1.9 ± 0.3‰ which we measured in plant-chamber experiments. Air samples with strong anthropogenic influence indicated an anthropogenic COS isotopic value of 8 ± 1‰. Samples of seawater-equilibrated-air indicate that the marine COS source has an isotopic value of 14.7 ± 1‰. Using our data-based mass balance, we constrained the relative contribution of the two main tropospheric COS sources resulting in 40 ± 17% for the anthropogenic source and 60 ± 20% for the oceanic source. This constraint is important for a better understanding of the global COS budget and its improved use for GPP determination.more » « less
-
Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) 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 synthesise datasets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy and cement production data. Emissions from land-use change (ELUC) are estimated by bookkeeping models based on land-use data. The global atmospheric CO2 growth rate (GATM) is computed from changes in concentration measured at surface stations. The global net uptake of CO2 by the ocean (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based fCO2-products. The global net uptake of CO2 by the land (SLAND) is estimated with dynamic global vegetation models. Additional lines of evidence are provided by atmospheric inversions, atmospheric oxygen measurements, ocean interior observation-based estimates, and Earth System Models. This year, we introduced corrections on the ELUC, SOCEAN and SLAND estimates. The sum of all sources and sinks results in the carbon budget imbalance (BIM), a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as ± 1σ. For the year 2024, EFOS increased by 1.1 % relative to 2023, with fossil emissions at 10.3 ± 0.5 GtC yr−1 (including the cement carbonation sink, 0.2 GtC yr−1), ELUC was 1.3 ± 0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.6 ± 0.9 GtC yr−1 (42.4 ± 3.2 GtCO2 yr−1). Also, for 2024, GATM was 7.9 ± 0.2 GtC yr−1 (3.73 ± 0.1 ppm yr−1), 2.2 GtC above the 2023 growth rate. SOCEAN was 3.4 ± 0.4 GtC yr−1 and SLAND was 1.9 ± 1.1 GtC yr−1, leaving a large negative BIM (−1.7 GtC yr−1), suggesting that the total sink or GATM is strongly overestimated in 2024. The global atmospheric CO2 concentration averaged over 2024 reached 422.8 ± 0.1 ppm. Preliminary data for 2025 suggest an increase in EFOS relative to 2024 of +1.0 % (0.2 % to 1.7 %) globally, and atmospheric CO2 concentration increasing by 2.1 ppm reaching 425.6 ppm, 53 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2024, with a near-zero overall budget imbalance, although discrepancies of up to around 1 GtC yr−1 persist for the representation of annual to decadal variability in CO2 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 CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the mean ocean sink. This living data update documents changes in methods and datasets applied to this most-recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2025 (Friedlingstein et al., 2025c).more » « less
-
Abstract. Gross primary productivity (GPP) is the largest carbon flux in the Earth system, playing a crucial role in removing atmospheric carbon dioxide and providing carbohydrates needed for ecosystem metabolism. Despite the importance of GPP, however, existing estimates present significant uncertainties and discrepancies. A key issue is the underrepresentation of the CO2 fertilization effect, a major factor contributing to the increased terrestrial carbon sink over recent decades. This omission could potentially bias our understanding of ecosystem responses to climate change. Here, we introduce CEDAR-GPP, the first global machine-learning-upscaled GPP product that incorporates the direct CO2 fertilization effect on photosynthesis. Our product is comprised of monthly GPP estimates and their uncertainty at 0.05° resolution from 1982 to 2020, generated using a comprehensive set of eddy covariance measurements, multi-source satellite observations, climate variables, and machine learning models. Importantly, we used both theoretical and data-driven approaches to incorporate the direct CO2 effects. Our machine learning models effectively predict monthly GPP (R2 ∼ 0.72), the mean seasonal cycles (R2 ∼ 0.77), and spatial variabilities (R2 ∼ 0.63) based on cross-validation at flux sites. After incorporating the direct CO2 effects, the predicted long-term GPP trend across global flux towers substantially increases from 3.1 to 4.5–5.4 gC m−2 yr−1, which aligns more closely with the 7.7 gC m−2 yr−1 trend detected from eddy covariance data. While the global patterns of annual mean GPP, seasonality, and interannual variability generally align with existing satellite-based products, CEDAR-GPP demonstrates higher long-term trends globally after incorporating CO2 fertilization and reflected a strong temperature control on direct CO2 effects. The estimated global GPP trend is 0.57–0.76 PgC yr−1 from 2001 to 2018 and 0.32–0.34 PgC yr−1 from 1982 to 2018. Estimating and validating GPP trends in data-scarce regions, such as the tropics, remains challenging, underscoring the importance of ongoing ground-based monitoring and advancements in modeling techniques. CEDAR-GPP offers a comprehensive representation of GPP temporal and spatial dynamics, providing valuable insights into ecosystem–climate interactions. The CEDAR-GPP product is available at https://doi.org/10.5281/zenodo.8212706 (Kang et al., 2024).more » « less
An official website of the United States government

