Carbonyl sulfide (COS) is an atmospheric trace gas of interest for C cycle research because COS uptake by continental vegetation is strongly related to terrestrial gross primary productivity (GPP), the largest and most uncertain flux in atmospheric CO2 budgets. However, to use atmospheric COS as an additional tracer of GPP, an accurate quantification of COS exchange by soils is also needed. At present, the atmospheric COS budget is unbalanced globally, with total COS flux estimates from oxic and anoxic soils that vary between -409 and -89 GgS yr(-1). This uncertainty hampers the use of atmospheric COS concentrations to constrain GPP estimates through atmospheric transport inversions. In this study we implemented a mechanistic soil COS model in the ORCHIDEE (Organising Carbon and Hydrology In Dynamic Ecosystems) land surface model to simulate COS fluxes in oxic and anoxic soils. Evaluation of the model against flux measurements at seven sites yields a mean root mean square deviation of 1.6 pmolm(-2)s(-1), instead of 2 pmol m(-2)s(-1) when using a previous empirical approach that links soil COS uptake to soil heterotrophic respiration. However, soil COS model evaluation is still limited by the scarcity of observation sites and long-term measurement periods, with all sites located in a latitudinal band between 39 and 62 degrees N and no observations during wintertime in this study. The new model predicts that, globally and over the 2009-2016 period, oxic soils act as a net uptake of -126 GgS yr(-1) and anoxic soils are a source of +96 GgS yr(-1), leading to a global net soil sink of only -30 GgS yr(-1), i.e. much smaller than previous estimates. The small magnitude of the soil fluxes suggests that the error in the COS budget is dominated by the much larger fluxes from plants, oceans, and industrial activities. The predicted spatial distribution of soil COS fluxes, with large emissions from oxic (up to 68.2 pmol COS m(-2) s(-1)) and anoxic (up to 36.8 pmol COS m(-2) S-1) soils in the tropics, especially in India and in the Sahel region, marginally improves the latitudinal gradient of atmospheric COS concentrations, after transport by the LMDZ (Laboratoire de Meteorologie Dynamique) atmospheric transport model. The impact of different soil COS flux representations on the latitudinal gradient of the atmospheric COS concentrations is strongest in the Northern Hemisphere. We also implemented spatiotemporal variations in near-ground atmospheric COS concentrations in the modelling of biospheric COS fluxes, which helped reduce the imbalance of the atmospheric COS budget by lowering soil COS uptake by 10 % and plant COS uptake by 8 % globally (with a revised mean vegetation budget of -576 GgS yr(-1) over 2009-2016). Sensitivity analyses highlighted the different parameters to which each soil COS flux model is the most responsive, selected in a parameter optimization framework. Having both vegetation and soil COS fluxes modelled within ORCHIDEE opens the way for using observed ecosystem COS fluxes and larger-scale atmospheric COS mixing ratios to improve the simulated GPP, through data assimilation techniques.< Leer menos