For the past decade, observations of carbonyl sulfide (OCS or COS) have been investigated as a proxy for carbon uptake by plants. OCS is destroyed by enzymes that interact with CO₂ during photosynthesis, namely carbonic anhydrase (CA) and RuBisCO, where CA is the more important one. The majority of sources of OCS to the atmosphere are geographically separated from this large plant sink, whereas the sources and sinks of CO₂ are co-located in ecosystems. The drawdown of OCS can therefore be related to the uptake of CO₂ without the added complication of co-located emissions comparable in magnitude. Here we review the state of our understanding of the global OCS cycle and its applications to ecosystem carbon cycle science. OCS uptake is correlated well to plant carbon uptake, especially at the regional scale. OCS can be used in conjunction with other independent measures of ecosystem function, like solar-induced fluorescence and carbon and water isotope studies. More work needs to be done to generate global coverage for OCS observations and to link this powerful atmospheric tracer to systems where fundamental questions concerning the carbon and water cycle remain.

The global budget of tropospheric carbonyl sulfide (OCS) is believed to be at equilibrium because background air concentrations have remained roughly stable over at least the last decade. Since the uptake of OCS by leaves (associated with photosynthesis) and soils have been revised significantly upwards recently, an equilibrated budget can only be obtained with a compensatory source of OCS. It has been assumed that the missing source of OCS comes from the low-latitude ocean, following the incident solar flux. The present work uses parameterizations of major production and removal processes of organic compounds in the NEMO-PISCES (Nucleus for European Modelling of the Ocean, Pelagic Interaction Scheme for Carbon and Ecosystem Studies) ocean general circulation and biogeochemistry model to assess the marine source of OCS. In addition, the OCS photo-production rates computed with the NEMO-PISCES model~were evaluated independently using the UV absorption coefficient of chromophoric dissolved organic matter (derived from satellite ocean color data) and apparent quantum yields available in the literature. Our simulations show global direct marine emissions of OCS in the range of 573–3997 GgS yr−1, depending mostly on the quantification of the absorption rate of chromophoric dissolved organic matter. The high estimates of that range are unlikely, as they correspond to a formulation that most likely overestimate photo-production process. Low and medium (813 GgS yr−1) estimates derived from the NEMO-PISCES model are however consistent spatially and temporally~with the suggested missing source of Berry et al. (2013), allowing us thus to close the global budget of OCS given the recent estimates of leaf and soil OCS uptake.

Clear analogies between carbonyl sulfide (OCS) and carbon dioxide (CO2) diffusion pathways through leaves have been revealed by experimental studies, with plant uptake playing an important role for the atmospheric budget of both species. Here we use atmospheric OCS to evaluate the gross primary production (GPP) of three dynamic global vegetation models (Lund–Potsdam–Jena, LPJ; National Center for Atmospheric Research – Community Land Model 4, NCAR-CLM4; and Organising Carbon and Hydrology In Dynamic Ecosystems, ORCHIDEE). Vegetation uptake of OCS is modeled as a linear function of GPP and leaf relative uptake (LRU), the ratio of OCS to CO2 deposition velocities of plants. New parameterizations for the non-photosynthetic sinks (oxic soils, atmospheric oxidation) and biogenic sources (oceans and anoxic soils) of OCS are also provided. Despite new large oceanic emissions, global OCS budgets created with each vegetation model show exceeding sinks by several hundred Gg S yr−1. An inversion of the surface fluxes (optimization of a global scalar which accounts for flux uncertainties) led to balanced OCS global budgets, as atmospheric measurements suggest, mainly by drastic reduction (up to −50 %) in soil and vegetation uptakes. The amplitude of variations in atmospheric OCS mixing ratios is mainly dictated by the vegetation sink over the Northern Hemisphere. This allows for bias recognition in the GPP representations of the three selected models. The main bias patterns are (i) the terrestrial GPP of ORCHIDEE at high northern latitudes is currently overestimated, (ii) the seasonal variations of the GPP are out of phase in the NCAR-CLM4 model, showing a maximum carbon uptake too early in spring in the northernmost ecosystems, (iii) the overall amplitude of the seasonal variations of GPP in NCAR-CLM4 is too small, and (iv) for the LPJ model, the GPP is slightly out of phase for the northernmost ecosystems and the respiration fluxes might be too large in summer in the Northern Hemisphere. These results rely on the robustness of the OCS modeling framework and, in particular, the choice of the LRU values (assumed constant in time) and the parameterization of soil OCS uptake with small seasonal variations. Refined optimization with regional-scale and seasonally varying coefficients might help to test some of these hypothesis.

Carbonyl sulfide (COS), a trace gas showing striking similarity to CO₂ in terms of biochemical diffusion pathway into leaves, has been recognized as a promising indicator of the plant gross primary production (GPP), the amount of carbon dioxide that is absorbed through photosynthesis by terrestrial ecosystems. However, large uncertainties about the other components of its atmospheric budget prevent us from directly relating the atmospheric COS measurements to GPP. The largest uncertainty comes from the closure of its atmospheric budget, with a source component missing. Here, we explore the benefit of assimilating both COS and CO₂ measurements into the LMDz atmospheric transport model to obtain consistent information on GPP, plant respiration and COS budget. To this end, we develop an analytical inverse system that optimizes biospheric fluxes for the 15 plant functional types (PFTs) defined in the ORCHIDEE global land surface model. Plant uptake of COS is parameterized as a linear function of GPP and of the leaf relative uptake (LRU), which is the ratio of COS to CO₂ deposition velocities in plants. A possible scenario for the period 2008–2019 leads to a global biospheric sink of 800 GgS yr⁻¹, with higher absorption in the high latitudes and higher oceanic emissions between 400 and 600 GgS yr⁻¹ most of which is located in the tropics. As for the CO₂ budget, the inverse system increases GPP in the high latitudes by a few GtC yr⁻¹ without modifying the respiration compared to the ORCHIDEE fluxes used as a prior. In contrast, in the tropics the system tends to weaken both respiration and GPP. The optimized components of the COS and CO₂ budgets have been evaluated against independent measurements over North America, the Pacific Ocean, at three sites in Japan and at one site in France. Overall, the posterior COS concentrations are in better agreement with the COS retrievals at 250 hPa from the MIPAS satellite and with airborne measurements made over North America and the Pacific Ocean. The system seems to have rightly corrected the underestimated GPP over the high latitudes. However, the change in seasonality of GPP in the tropics disagrees with solar-induced fluorescence (SIF) data. The decline in biospheric sink in the Amazon driven by the inversion also disagrees with MIPAS COS retrievals at 250 hPa, highlighting the lack of observational constraints in this region. Moreover, the comparison with the surface measurements in Japan and France suggests misplaced sources in the prior anthropogenic inventory, emphasizing the need for an improved inventory to better partition oceanic and continental sources in Asia and Europe.

Land surface modellers need measurable proxies to constrain the quantity of carbon dioxide (CO₂) assimilated by continental plants through photosynthesis, known as gross primary production (GPP). Carbonyl sulfide (COS), which is taken up by leaves through their stomates and then hydrolysed by photosynthetic enzymes, is a candidate GPP proxy. A former study with the ORCHIDEE land surface model used a fixed ratio of COS uptake to CO₂ uptake normalised to respective ambient concentrations for each vegetation type (leaf relative uptake, LRU) to compute vegetation COS fluxes from GPP. The LRU approach is known to have limited accuracy since the LRU ratio changes with variables such as photosynthetically active radiation (PAR): while CO₂ uptake slows under low light, COS uptake is not light limited. However, the LRU approach has been popular for COS–GPP proxy studies because of its ease of application and apparent low contribution to uncertainty for regional-scale applications. In this study we refined the COS–GPP relationship and implemented in ORCHIDEE a mechanistic model that describes COS uptake by continental vegetation. We compared the simulated COS fluxes against measured hourly COS fluxes at two sites and studied the model behaviour and links 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⁻¹, in the middle range of former studies (−490 to −1335 Gg S yr⁻¹). Based on monthly mean fluxes simulated by the mechanistic approach in ORCHIDEE, we derived new LRU values for the different vegetation types, ranging between 0.92 and 1.72, close to recently published averages for observed values of 1.21 for C₄ and 1.68 for C₃ plants. We transported the COS using the monthly vegetation COS fluxes derived from both the mechanistic and the LRU approaches, and we evaluated the simulated COS concentrations at NOAA sites. Although the mechanistic approach was more appropriate when comparing to high-temporal-resolution COS flux measurements, both approaches gave similar results when transporting with monthly COS fluxes and evaluating COS concentrations at stations. In our study, uncertainties between these two approaches are of secondary importance compared to the uncertainties in the COS global budget, which are currently a limiting factor to the potential of COS concentrations to constrain GPP simulated by land surface models on the global scale.