Inverse modelling of carbonyl sulfide: implementation, evaluation and implications for the global budget

Abstract. Carbonyl sulfide (COS) has the potential to be used as a climate diagnostic due to its close coupling to the biospheric uptake of CO2 and its role in the formation of stratospheric aerosol. The current understanding of the COS budget, however, lacks COS sources, which have previously been allocated to the tropical ocean. This paper presents a first attempt of global inverse modelling of COS within the 4-Dimensional variational data-assimilation system of the TM5 chemistry transport model (TM5-4DVAR) and a comparison of the results with independent COS observations. We focus on the global COS budget, including COS production from its precursors carbon disulfide (CS2) and dimethyl sulfide (DMS). To this end, we implemented COS uptake by soil and vegetation from an updated biosphere model (SiB4), and new inventories for anthropogenic and biomass burning emissions. The model framework is capable of closing the COS budget by optimizing for missing emissions using NOAA observations in the period 2000–2012. The addition of 432 Gg S a−1 COS is required to obtain a good fit with NOAA observations. This missing source shows little year-to-year variations, but considerable seasonal variations. We found that the missing sources are likely located in the tropical regions, and an overestimated biospheric sink in the tropics cannot be ruled out. Moreover, high latitudes in the Northern Hemisphere require extra COS uptake or reduced emissions. HIPPO aircraft observations, NOAA airborne profiles from an ongoing monitoring program, and several satellite data sources are used to evaluate the optimized model results. This evaluation indicates that COS in the free troposphere remains underestimated after optimization. Assimilation of HIPPO observations slightly improves this model bias, which implies that additional observations are urgently required to constrain sources and sinks of COS. We finally find that the biosphere flux dependency on surface COS mixing ratio may substantially lower the fluxes of the SiB4 biosphere model over strong uptake regions. In planned further studies we will implement this biosphere dependency, and additionally assimilate satellite data with the aim to better separate the role of the oceans and the biosphere in the global COS budget.