Accurately predicting terrestrial ecosystem responses to climate change is crucial for addressing global challenges. This relies on mechanistic modelling of ecosystem processes through Land Surface Models (LSMs). Despite their importance, LSMs face significant uncertainties due to poorly constrained parameters, especially in carbon cycle predictions. This paper reviews the progress made in using data assimilation (DA) for LSM parameter optimisation, focusing on carbon-water-vegetation interactions, as well as discussing the technical challenges faced by the community. These challenges include identifying sensitive model parameters and their prior distributions, characterising errors due to observation biases and model-data inconsistencies, developing observation operators to interface between the model and the observations, tackling spatial and temporal heterogeneity as well as dealing with large and multiple datasets, and including the spin-up and historical period in the assimilation window. We then outline how machine learning (ML) can help address these issues, proposing different avenues for future work that integrate ML and DA to reduce uncertainties in LSMs. We conclude by highlighting future priorities, including the need for international collaborations, to fully leverage the wealth of available Earth observation datasets, harness machine learning advances, and enhance the predictive capabilities of LSMs.

Comprehensive land models are subject to significant parametric uncertainty, which can be hard to quantify due to the large number of parameters and high model computational costs. We constructed a large parameter perturbation ensemble (PPE) for the Community Land Model vesion 5.1 with biogeochemistry configuration (CLM5.1-BGC). We performed more than 2000 simulations perturbing 211 parameters across six forcing scenarios. This provides an expansive dataset, which can be used to identify the most influential parameters on a wide range of output variables globally, by biome, or by plant functional type. We found that parameter effects can exceed scenario effects and that a small number of parameters explains a large fraction of variance across our ensemble. The most important parameters can differ regionally and also based on the forcing scenario. The software infrastructure developed for this experiment has greatly reduced the human and computer time needed for CLM PPEs, which can facilitate routine investigation of parameter sensitivity and uncertainty, as well as automated calibration.

The Arctic hydrological system is an interconnected system that is experiencing rapid change. It is comprised of permafrost, snow, glacier, frozen soils, and inland river systems. Permafrost degradation, trends towards earlier snow melt, a lengthening snow-free season, soil ice melt, and warming frozen soils all challenge hydrologic simulation under climate change in the Arctic. In this study, we provide an improved representation of the hydrologic cycle across a regional Arctic domain using a generalizable optimization methodology and workflow for the community. We applied the Community Terrestrial Systems Model (CTSM) across the US state of Alaska and the Yukon River Basin at 4-km spatial resolution. We highlight several potentially useful high-resolution CTSM configuration changes. Additionally, we performed a multi-objective optimization using snow and river flow metrics within an adaptive surrogate-based model optimization scheme. Four representative river basins across our study domain were selected for optimization based on observed streamflow and snow water equivalent observations at ten SNOTEL sites. Fourteen sensitive parameters were identified for optimization with half of them not directly related to hydrology or snow processes. Across fifteen out-of-sample river basins, thirteen had improved flow simulations after optimization and the median Kling-Gupta Efficiency of daily flow increased from 0.40 to 0.63. In addition, we adapted the Shapley Decomposition to disentangle each parameter’s contribution to streamflow performance changes, with the seven non-hydrological parameters providing a non-negligible contribution to performance gains. The snow simulation had limited improvement, likely because snow simulation is influenced more by meteorological forcing than model parameter choices.