The CEROBI project developed a robust framework for simulating forest drought responses in the ORCHIDEE Land Surface Model through substantial advances in initialization and validation. These methodological improvements revealed that significant issues affecting forest carbon sink predictions must be addressed before isotopic constraints can be meaningfully applied to the model.
1. Model Initialization and Demographic State Variables
Hydraulic architecture calculations are fundamentally dependent on tree height and stand density, as these structural variables determine hydraulic resistance within the soil-plant-atmosphere continuum and govern water transport dynamics. If tree heights and stand structure are unrealistic, the simulated hydraulic architecture becomes flawed, leading to inaccurate representations of plant water stress. Because isotopic fractionation is closely linked to stomatal conductance and transpiration—both controlled by water stress—realistic hydraulic architecture is essential for any subsequent isotopic analysis to be physiologically plausible. To address this, an initialization procedure called "blended spin-up" was developed to match simulated tree diameter distributions with European forest inventory data from 2002-2010. The method was refined to blend tree diameters separately for each Plant Functional Type (PFT) to ensure realistic species-specific size structures. However, this refinement initially violated carbon mass conservation because ORCHIDEE simulates wood product pools at the grid cell level rather than the PFT level. Code modifications were required to restore mass conservation while maintaining demographic realism, ensuring diameter adjustments remained synchronized with carbon and nitrogen pools in the vegetation, litter, soil, and wood products.
2. Impact on Carbon Flux Trajectories
To illustrate the impact of this model spin-up, two parallel 100-year simulations of European forests were conducted with identical climate forcing, comparing standard initialization against the refined PFT-specific diameter blended spin-up. The results revealed a substantial impact on long-term carbon cycle projections, where the non-initialized simulation showed European forests as a net carbon source over the century, while the simulation starting from the blended spin-up showed them acting as a net carbon sink. This reversal in the sign of the carbon flux has profound implications for climate mitigation assessments. This difference occurs because initial tree size distributions determine current biomass stocks, mortality vulnerability, growth potential, and demographic turnover rates. Consequently, incorrect initialization creates unrealistic demographic structures that propagate through decades of simulation, producing fundamentally wrong trajectories and highlighting that demographic initialization is a first-order control on model predictions rather than a technical detail.
3. Model Evaluation and Observational Uncertainty
Rigorous assessment of model performance was conducted to validate the initial model states and subsequent forest dynamics over Europe. Since a reasonable relationship is expected between forest age and structure in managed forests, the model's capability to reproduce observed age structure was evaluated against two forest age datasets. Spatial patterns of age were well captured, correctly representing the geographic distribution of young regenerating forests versus older stands across Northern Europe, though specific age values were sometimes over- or underestimated. Discrepancies between the two observational datasets were comparable in magnitude to the differences between the model and observations, emphasizing that model validation is constrained by substantial uncertainty in the benchmark data themselves. Furthermore, while the blended spin-up utilized ground-based diameter products, it was necessary to demonstrate that simulated diameters aligned with acceptable tree heights and stand densities. Aboveground biomass was used as an independent constraint since it is the product of diameter squared, height, stand density, and wood density. The initialization procedure resulted in a 50% reduction in errors for aboveground biomass across Europe compared to non-initialized simulations, though residual biases remained in regions like Scandinavia. These discrepancies stem from uncertainties in the forest inventory, inconsistencies between in-situ and satellite-derived estimates, and potential structural limitations in the model's allometry and biomass allocation schemes. To independently validate growth dynamics, forest state variables were initialized for 1990 and run forward to 2010, showing good agreement between simulated and observed diameter distributions.
4. Assessment of Forest Response to Drought
Forest resilience was evaluated by testing the ability of the ORCHIDEE model to simulate successional biomass recovery following drought, which was represented here as a clear-cut. The model demonstrated high accuracy in temperate forests, where stomatal closure and cavitation thresholds were identified as the primary physiological parameters governing the trade-off between carbon starvation and hydraulic failure. Ultimately, CEROBI established a foundation for future research by providing a validated model infrastructure with realistic initial conditions and identifying the key mechanisms governing how forests respond to drought.