skip to main content


Title: 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
NSF-PAR ID:
10311393
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Biogeosciences
Volume:
18
Issue:
9
ISSN:
1726-4189
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO2 fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5° grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (Vcmax) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) (r2 =  0.76; Nash–Sutcliffe modeling efficiency, MEF  =  0.76) and ecosystem respiration (ER, r2 =  0.78, MEF  =  0.75), with lesser accuracy for latent heat fluxes (LE, r2 =  0.42, MEF  =  0.14) and and net ecosystem CO2 exchange (NEE, r2 =  0.38, MEF  =  0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r2 values (0.57–0.86). For spatial across-site gradients of annual mean GPP, ER, NEE, and LE, r2 values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted (r2 < 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized Vcmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average Vcmax value. 
    more » « less
  4. Abstract. In the global methane budget, the largest natural sourceis attributed to wetlands, which encompass all ecosystems composed ofwaterlogged or inundated ground, capable of methane production. Among them,northern peatlands that store large amounts of soil organic carbon have beenfunctioning, since the end of the last glaciation period, as long-termsources of methane (CH4) and are one of the most significant methanesources among wetlands. To reduce uncertainty of quantifying methane flux in theglobal methane budget, it is of significance to understand the underlyingprocesses for methane production and fluxes in northern peatlands. A methanemodel that features methane production and transport by plants, ebullitionprocess and diffusion in soil, oxidation to CO2, and CH4 fluxes tothe atmosphere has been embedded in the ORCHIDEE-PEAT land surface modelthat includes an explicit representation of northern peatlands.ORCHIDEE-PCH4 was calibrated and evaluated on 14 peatland sites distributedon both the Eurasian and American continents in the northern boreal andtemperate regions. Data assimilation approaches were employed to optimizedparameters at each site and at all sites simultaneously. Results show thatmethanogenesis is sensitive to temperature and substrate availability overthe top 75 cm of soil depth. Methane emissions estimated using single siteoptimization (SSO) of model parameters are underestimated by 9 g CH4 m−2 yr−1 on average (i.e., 50 % higher than the site average ofyearly methane emissions). While using the multi-site optimization (MSO),methane emissions are overestimated by 5 g CH4 m−2 yr−1 onaverage across all investigated sites (i.e., 37 % lower than the siteaverage of yearly methane emissions). 
    more » « less
  5. 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