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 more » 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. « less
Authors:
; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
1848618
Publication Date:
NSF-PAR ID:
10311393
Journal Name:
Biogeosciences
Volume:
18
Issue:
9
ISSN:
1726-4189
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 productionmore »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.« less
  2. 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 meanmore »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.« less
  3. 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 thanmore »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).« less
  4. Elger, Kirsten ; Carlson, David ; Klump, Jens ; Peng, Ge (Ed.)
    Air-sea flux of carbon dioxide (CO2) is a critical component of the global carbon cycle and the climate system with the ocean removing about a quarter of the CO2 emitted into the atmosphere by human activities over the last decade. A common approach to estimate this net flux of CO2 across the air-sea interface is the use of surface ocean CO2 observations and the computation of the flux through a bulk parameterization approach. Yet, the details for how this is done in order to arrive at a global ocean CO2 uptake estimate varies greatly, unnecessarily enhancing the uncertainties. Here we reduce some of these uncertainties by harmonizing an ensemble of products that interpolate surface ocean CO2 bservations to near global coverage. We propose a common methodology to fill in missing areas in the products and to calculate fluxes and present a new estimate of the net flux. The ensemble data product, SeaFlux (Gregor & Fay (2021), doi.org/10.5281/zenodo.4133802, https://github.com/luke-gregor/SeaFlux), accounts for the diversity of the underlying mapping methodologies. Utilizing six 30 global observation-based mapping products (CMEMS-FFNN, CSIR-ML6, JENA-MLS, JMA-MLR, MPI-SOMFFN, NIESFNN), the SeaFlux ensemble approach adjusts for methodological inconsistencies in flux calculations that can result in an average error of 15%more »in global mean flux estimates. We address differences in spatial coverage of the surface ocean CO2 between the mapping products which ultimately yields an increase in CO2 uptake of up to 19% for some products. Fluxes are calculated using three wind products (CCMPv2, ERA5, and JRA55). Application of an appropriately scaled gas exchange 35 coefficient has a greater impact on the resulting flux than solely the choice of wind product. With these adjustments, we derive an improved ensemble of surface ocean pCO2 and air-sea carbon flux estimates. The SeaFlux ensemble suggests a global mean uptake of CO2 from the atmosphere of 1.92 +/- 0.35 PgC yr-1. This work aims to support the community effort to perform model-data intercomparisons which will help to identify missing fluxes as we strive to close the global carbon budget.« less
  5. Abstract. Plant activity in semi-arid ecosystems is largely controlled by pulses of precipitation, making them particularly vulnerable to increased aridity expected with climate change. Simple bucket-model hydrology schemes in land surface models (LSMs) have had limited ability in accurately capturing semi-arid water stores and fluxes. Recent, more complex, LSM hydrology models have not been widely evaluated against semi-arid ecosystem in situ data. We hypothesize that the failure of older LSM versions to represent evapotranspiration, ET, in arid lands is because simple bucket models do not capture realistic fluctuations in upper layer soil moisture. We therefore predict that including a discretized soil hydrology scheme based on a mechanistic description of moisture diffusion will result in an improvement in model ET when compared to data because the temporal variability of upper layer soil moisture content better corresponds to that of precipitation inputs. To test this prediction, we compared ORCHIDEE LSM simulations from (1) a simple conceptual 2-layer bucket scheme with fixed hydrological parameters; and (2) a 11-layer discretized mechanistic scheme of moisture diffusion in unsaturated soil based on Richards equations against daily and monthly soil moisture and ET observations, together with data-derived transpiration / evaporation, T / ET, ratios, from six semi-arid grass, shrub and forestmore »sites in the southwestern USA. The 11-layer scheme also has modified calculations of surface runoff, bare soil evaporation, and water limitation to be compatible with the more complex hydrology configuration. To diagnose remaining discrepancies in the 11-layer model, we tested two further configurations: (i) the addition of a term that captures bare soil evaporation resistance to dry soil; and (ii) reduced bare soil fraction. We found that the more mechanistic 11-layer model results better representation of the daily and monthly ET observations. We show that is likely because of improved simulation of soil moisture in the upper layers of soil (top 5 cm). Some discrepancies between observed and modelled soil moisture and ET may allow us to prioritize future model development. Adding a soil resistance term generally decreased simulated E and increased soil moisture content, thus increasing T and T / ET ratios and reducing the negative T / ET model-data bias. By reducing the bare soil fraction in the model, we illustrated that modelled leaf T is too low at sparsely vegetated sites. We conclude that a discretized soil hydrology scheme and associated developments improves estimates of ET by allowing the model to more closely match the pulse precipitation dynamics of these semi-arid ecosystems; however, the partitioning of T from bare soil evaporation is not solved by this modification alone.« less