Land-surface models (LSMs) are crucial components of the Earth system models (ESMs) that are used to make coupled climate–carbon cycle projections for the 21st century. The Joint UK Land Environment Simulator (JULES) is the land-surface model used in the climate and weather forecast models of the UK Met Office. JULES is also extensively used offline as a land-surface impacts tool, forced with climatologies into the future. In this study, JULES is automatically differentiated with respect to JULES parameters using commercial software from FastOpt, resulting in an analytical gradient, or adjoint, of the model. Using this adjoint, the adJULES parameter estimation system has been developed to search for locally optimum parameters by calibrating against observations. This paper describes adJULES in a data assimilation framework and demonstrates its ability to improve the model–data fit using eddy-covariance measurements of gross primary production (GPP) and latent heat (LE) fluxes. adJULES also has the ability to calibrate over multiple sites simultaneously. This feature is used to define new optimised parameter values for the five plant functional types (PFTs) in JULES. The optimised PFT-specific parameters improve the performance of JULES at over 85 % of the sites used in the study, at both the calibration and evaluation stages. The new improved parameters for JULES are presented along with the associated uncertainties for each parameter.

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.

The rate at which land surface soils dry following rain events is an important feature of terrestrial models. It determines, for example, the water availability for vegetation, the occurrences of droughts, and the surface heat exchanges. As such, surface soil moisture (SSM) “drydowns,” i.e., the SSM temporal dynamics following a significant rainfall event, are of particular interest when evaluating and calibrating land surface models (LSMs). By investigating drydowns, characterized by an exponential decay time scale τ, we aim to improve the representation of SSM in the ORCHIDEE global LSM. We consider τ calculated over 18 International Soil Moisture Network sites found within the footprint of FLUXNET towers, covering different vegetation types and climates. Using the ORCHIDEE LSM, we compare τ from the modeled SSM time series to values computed from in situ SSM measurements. We then assess the potential of using τ observations to constrain some water, carbon, and energy parameters of ORCHIDEE, selected using a sensitivity analysis, through a standard Bayesian optimization procedure. The impact of the SSM optimization is evaluated using FLUXNET evapotranspiration and gross primary production (GPP) data. We find that the relative drydowns of SSM can be well calibrated using observation-based τ estimates, when there is no need to match the absolute observed and modeled SSM values. When evaluated using independent data, τ-calibration parameters were able to improve drydowns for 73% of the sites. Furthermore, the fit of the model to independent fluxes was only minutely changed. We conclude by considering the potential of global satellite products to scale up the experiment to a global-scale optimization.

Peatlands at high latitudes have accumulated >400 Pg carbon (C) because saturated soil and cold temperatures suppress C decomposition. This substantial amount of C in Arctic and Boreal peatlands is potentially subject to increased decomposition if the water table (WT) decreases due to climate change, including permafrost thaw‐related drying. Here, we optimize a version of the Organizing Carbon and Hydrology In Dynamic Ecosystems model (ORCHIDEE‐PCH4) using site‐specific observations to investigate changes in CO₂ and CH₄ fluxes as well as C stock responses to an experimentally manipulated decrease of WT at six northern peatlands. The unmanipulated control peatlands, with the WT 2 kg C m⁻² over 100 years when WT is lowered by 50 cm, while permafrost peatlands temporally switched from C sinks to sources. These results highlight that reductions in C storage capacity in response to drying of northern peatlands are offset in part by reduced CH₄ emissions, thus slightly reducing the positive carbon climate feedbacks of peatlands under a warmer and drier future climate scenario.

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 CO₂ 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⁻¹. 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 pmol m⁻² s⁻¹, instead of 2 pmol m⁻² s⁻¹ 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^∘ 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⁻¹ and anoxic soils are a source of +96 GgS yr⁻¹, leading to a global net soil sink of only −30 GgS yr⁻¹, 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⁻² s⁻¹) and anoxic (up to 36.8 pmol COS m⁻² s⁻¹) 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 Météorologie 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⁻¹ 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.