Search

Current climate warming is accelerating mass loss from glaciers and ice sheets. In Greenland, the rates of mass changes are now dominated by changes in surface mass balance (SMB) due to increased surface melting. To improve the future sea-level rise projections, it is therefore critical to have an accurate estimate of the SMB, which depends on the representation of the processes occurring within the snowpack. The Explicit Snow (ES) scheme implemented in the land surface model Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) has not yet been adapted to ice-covered areas. Here, we present the preliminary developments we made to apply the ES model to glaciers and ice sheets. Our analysis mainly concerns the model's ability to represent ablation-related processes. At the regional scale, our results are compared to the MAR regional atmospheric model outputs and to MODIS albedo retrievals. Using different albedo parameterizations, we performed offline ES simulations forced by the MAR model over the 2000–2019 period. Our results reveal a strong sensitivity of the modelled SMB components to the albedo parameterization. Results inferred with albedo parameters obtained using a manual tuning approach present very good agreement with the MAR outputs. Conversely, with the albedo parameterization used in the standard ORCHIDEE version, runoff and sublimation were underestimated. We also tested parameters found in a previous data assimilation experiment, calibrating the ablation processes using MODIS snow albedo. While these parameters greatly improve the modelled albedo over the entire ice sheet, they degrade the other model outputs compared to those obtained with the manually tuned approach. This is likely due to the model overfitting to the calibration albedo dataset without any constraint applied to the other processes controlling the state of the snowpack. This underlines the need to perform a “multi-objective” optimization using auxiliary observations related to internal snowpack processes. Although there is still room for further improvements, the developments reported in the present study constitute an important advance in assessing the Greenland SMB with possible extension to mountain glaciers or the Antarctic ice sheet.

With the growing complexity of land surface models used to represent the terrestrial part of wider Earth system models, the need for sophisticated and robust parameter optimisation techniques is paramount. Quantifying parameter uncertainty is essential for both model development and more accurate projections. In this study, we assess the power of history matching by comparing results to the variational data assimilation approach commonly used in land surface models for parameter estimation. Although both approaches have different setups and goals, we can extract posterior parameter distributions from both methods and test the model–data fit of ensembles sampled from these distributions. Using a twin experiment, we test whether we can recover known parameter values. Through variational data assimilation, we closely match the observations. However, the known parameter values are not always contained in the posterior parameter distribution, highlighting the equifinality of the parameter space. In contrast, while more conservative, history matching still gives a reasonably good fit and provides more information about the model structure by allowing for non-Gaussian parameter distributions. Furthermore, the true parameters are contained in the posterior distributions. We then consider history matching's ability to ingest different metrics targeting different physical parts of the model, thus helping to reduce the parameter space further and improve the model–data fit. We find the best results when history matching is used with multiple metrics; not only is the model–data fit improved, but we also gain a deeper understanding of the model and how the different parameters constrain different parts of the seasonal cycle. We conclude by discussing the potential of history matching in future studies.

Predicting the responses of terrestrial ecosystem carbon to future global change strongly relies on our ability to model accurately the underlying processes at a global scale. However, terrestrial biosphere models representing the carbon and nitrogen cycles and their interactions remain subject to large uncertainties, partly because of unknown or poorly constrained parameters. Parameter estimation is a powerful tool that can be used to optimise these parameters by confronting the model with observations. In this paper, we identify sensitive model parameters from a recent version of the ORgainzing Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE) land surface model that includes the nitrogen cycle. These sensitive parameters include ones involved in parameterisations controlling the impact of the nitrogen cycle on the carbon cycle and, in particular, the limitation of photosynthesis due to leaf nitrogen availability. We optimise these ORCHIDEE parameters against carbon flux data collected on sites from the FLUXNET network. However, optimising against present-day observations does not automatically give us confidence in future projections of the model, given that environmental conditions are likely to shift compared to the present day. Manipulation experiments give us a unique look into how the ecosystem may respond to future environmental changes. One such type of manipulation experiment, the Free Air CO2 Enrichment (FACE) experiment, provides a unique opportunity to assess vegetation response to increasing CO2 by providing data under ambient and elevated CO2 conditions. Therefore, to better capture the ecosystem response to increased CO2, we add the data from two FACE sites to our optimisations, in addition to the FLUXNET data. We use data from both CO2 conditions of FACE, which allows us to gain extra confidence in the model simulations using this set of parameters. We find that we are able to improve the magnitude of modelled productivity. Although we are unable to correct the interannual variability fully, we start to simulate possible progressive nitrogen limitation at one of the sites. Using an idealised simulation experiment based on increasing atmospheric CO2 by 1 % yr−1 over 100 years, we find that optimising against only FLUXNET data tends to imply a large fertilisation effect, whereas optimising against FLUXNET and FACE data (with information about nutrient limitation and acclimation of plants) decreases it significantly.

The role of the land carbon cycle in climate change remains highly uncertain. A key source of the projection spread is related to the assumed response of photosynthesis to warming, especially in the tropics. The optimum temperature for photosynthesis determines whether warming positively or negatively impacts photosynthesis, thereby amplifying or suppressing CO2 fertilisation of photosynthesis under CO2-induced global warming. Land carbon cycle models have been extensively calibrated against local eddy flux measurements, but this has not previously been clearly translated into a reduced uncertainty in terms of how the tropical land carbon sink will respond to warming. Using a previous parameter perturbation ensemble carried out with version 3 of the Hadley Centre coupled climate–carbon cycle model (HadCM3C), we identify an emergent relationship between the optimal temperature for photosynthesis, which is especially relevant in tropical forests, and the projected amount of atmospheric CO2 at the end of the century. We combine this with a constraint on the optimum temperature for photosynthesis, derived from eddy covariance measurements using the adjoint of the Joint UK Land Environment Simulator (JULES) land surface model. Taken together, the emergent relationship from the coupled model and the constraint on the optimum temperature for photosynthesis define an emergent constraint on future atmospheric CO2 in the HadCM3C coupled climate–carbon cycle under a common emissions scenario (A1B). The emergent constraint sharpens the probability density of simulated CO2 change (2100–1900) and moves its peak to a lower value of 497 ± 91 compared to 607 ± 128 ppmv (parts per million by volume) when using the equal-weight prior. Although this result is likely to be model and scenario dependent, it demonstrates the potential of combining the large-scale emergent constraint approach with a parameter estimation using detailed local measurements.

Greenland ice sheet mass loss continues to accelerate as global temperatures increase. The surface albedo of the ice sheet determines the amount of absorbed solar energy, which is a key factor in driving surface snow and ice melting. Satellite-retrieved snow albedo allows us to compare and optimise modelled albedo over the entirety of the ice sheet. We optimise the parameters of the albedo scheme in the ORCHIDEE (Organizing Carbon and Hydrology in Dynamic Ecosystems) land surface model for 3 random years taken over the 2000–2017 period and validate over the remaining years. In particular, we want to improve the albedo at the edges of the ice sheet, since they correspond to ablation areas and show the greatest variations in runoff and surface mass balance. By giving a larger weight to points at the ice sheet's edge, we improve the model–data fit by reducing the root-mean-square deviation by over 25 % for the whole ice sheet for the summer months. This improvement is consistent for all years, even those not used in the calibration step. We also show the optimisation successfully improves the model–data fit at 87.5 % of in situ sites from the PROMICE (Programme for Monitoring of the Greenland Ice Sheet) network. We conclude by showing which additional model outputs are impacted by changes to the albedo parameters, encouraging future work using multiple data streams when optimising these parameters.

Land–atmosphere (L–A) interactions encompass the co-evolution of the land surface and overlying planetary boundary layer, primarily during daylight hours. However, many studies have been conducted using monthly or entire-day mean time series due to the lack of subdaily data. It is unclear whether the inclusion of nighttime data alters the assessment of L–A coupling or obscures L–A interactive processes. To address this question, we generate monthly (M), entire-day mean (E), and daytime-only mean (D) data based on the ERA5 (5th European Centre for Medium-Range Weather Forecasts reanalysis) product and evaluate the strength of L–A coupling through two-legged metrics, which partition the impact of the land states on surface fluxes (the land leg) from the impact of surface fluxes on the atmospheric states (the atmospheric leg). Here we show that the spatial patterns of strong L–A coupling regions among the M-, D-, and E-based diagnoses can differ by more than 80 %. The signal loss from E- to M-based diagnoses is determined by the memory of local L–A states. The differences between E- and D-based diagnoses can be driven by physical mechanisms or averaging algorithms. To improve understanding of L–A interactions, we call attention to the urgent need for more high-frequency data from both simulations and observations for relevant diagnoses. Regarding model outputs, two approaches are proposed to resolve the storage dilemma for high-frequency data: (1) integration of L–A metrics within Earth system models, and (2) producing alternative daily datasets based on different averaging algorithms.

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 pmol m−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∘ 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 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−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.

Land surface modellers need measurable proxies to constrain the quantity of carbon dioxide (CO2) 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 CO2 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 CO2 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−1, in the middle range of former studies (−490 to −1335 Gg S yr−1). 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 C4 and 1.68 for C3 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.

Earth System Models (ESMs) continue to diagnose a wide range of carbon budgets for each level of global warming. Here, we present emergent constraints on the carbon budget as a function of global warming, which combine the available ESM historical simulations and future projections for a range of scenarios, with observational estimates of global warming and anthropogenic CO 2 emissions to the present day. We estimate mean and likely ranges for cumulative carbon budgets for the Paris targets of 1.5 °C and 2 °C of global warming of 812 [691, 933] PgC and 1048 [881, 1216] PgC, which are more than 10% larger than the ensemble mean values from the CMIP6 models. The linearity between cumulative emissions and global warming is found to be maintained at least until 4 °C, and is consistent with an effective Transient Climate Response to Emissions (eTCRE) of 2.1 [1.8, 2.6] °C/1000PgC, from a global warming of 1.2 °C onwards., (© 2024. The Author(s).)