Abstract:Terrestrial ecosystem models can provide major insights into the responses of Earth's ecosystems to environmental changes and rising levels of atmospheric CO2. To achieve this goal, biosphere models need mechanistic formulations of the processes that drive the ecosystem functioning from diurnal to decadal timescales. However, the subsequent complexity of model equations is associated with unknown or poorly calibrated parameters that limit the accuracy of long-term simulations of carbon or water fluxes and their interannual variations. In this study, we develop a data assimilation framework to constrain the parameters of a mechanistic land surface model (ORCHIDEE) with eddy-covariance observations of CO2 and latent heat fluxes made during the years 2001-2004 at the temperate beech forest site of Hesse, in eastern France. As a first technical issue, we show that for a complex process-based model such as ORCHIDEE with many (28) parameters to be retrieved, a Monte Carlo approach (genetic algorithm, GA) provides more reliable optimal parameter values than a gradient-based minimization algorithm (variational scheme). The GA allows the global minimum to be found more efficiently, whilst the variational scheme often provides values relative to local minima. The ORCHIDEE model is then optimized for each year, and for the whole 2001-2004 period. We first find that a reduced (< 10) set of parameters can be tightly constrained by the eddy-covariance observations, with a typical error reduction of 90 %. We then show that including contrasted weather regimes (dry in 2003 and wet in 2002) is necessary to optimize a few specific parameters (like the temperature dependence of the photosynthetic activity). Furthermore, we find that parameters inverted from 4 years of flux measurements are successful at enhancing the model fit to the data on several timescales (from monthly to interannual), resulting in a typical modeling efficiency of 92% over the 2001-2004 period (Nash-Sutcliffe coefficient). This suggests that ORCHIDEE is able robustly to predict, after optimization, the fluxes of CO2 and the latent heat of a specific temperate beech forest (Hesse site). Finally, it is shown that using only 1 year of data does not produce robust enough optimized parameter sets in order to simulate properly the year-to-year flux variability. This emphasizes the need to assimilate data over several years, including contrasted weather regimes, to improve the simulated flux interannual variability.

rsofun v4.4.1: A model-data integration framework for simulating ecosystem processes

pnetr: An R package for the PnET family of forest ecosystem models

A differentiable, physics-informed ecosystem modeling and learning framework for large-scale inverse problems: demonstration with photosynthesis simulations

P-model V1.0: an Optimality-Based Light Use Efficiency Model for Simulating Ecosystem Gross Primary Production

Ecosystem model optimization using in situ flux observations: benefit of Monte Carlo versus variational schemes and analyses of the year-to-year model performances

A model-data intercomparison of CO2 exchange across North America: Results from the North American Carbon Program site synthesis

The value of soil respiration measurements for interpreting and modeling terrestrial carbon cycling

Forest Carbon Flux Simulation Using Multi-Source Data and Incorporation of Remotely Sensed Model with Process-Based Model

Optimization of ecosystem model parameters through assimilating eddy covariance flux data with an ensemble Kalman filter

R3DFEM: an R package for running the 3D-CMCC-FEM model

Ecosystem C and N cycle interactions – diverse model representations and divergent model predictions versus collective empirical constraints

Scientific modelling can be accessible, interoperable and user friendly: A case study for pasture and livestock modelling in Spain

ss3sim: An R package for fisheries stock assessment simulation with Stock Synthesis

Identification of Vegetation and Soil Carbon Pools out of Equilibrium in a Process Model Via Eddy Covariance and Biometric Constraints

Optimizing a Process‐based Ecosystem Model with Eddy‐covariance Flux Measurements: A Pine Forest in Southern France

Parameter optimization, sensitivity, and uncertainty analysis of an ecosystem model at a forest flux tower site in the United States

Assimilating Sentinel-2 Data in a Modified Vegetation Photosynthesis and Respiration Model (VPRM) to Improve the Simulation of Croplands CO2 Fluxes in Europe

CN-model: A dynamic model for the coupled carbon and nitrogen cycles in terrestrial ecosystems

The XSO framework (v0.1) and Phydra library (v0.1) for a flexible, reproducible, and integrated plankton community modeling environment in Python

EEMtoolbox: A user-friendly R package for flexible ensemble ecosystem modeling

Assimilating Sentinel-2 data in a modified vegetation photosynthesis and respiration model (VPRM) to improve the simulation of croplands CO 2 fluxes in Europe