Land models, which have been developed by the modeling community in the past few decades to predict future states of ecosystems and climate, have to be critically evaluated for their performance skills of simulating ecosystem responses and feedback to climate change. Benchmarking is an emerging procedure to measure performance of models against a set of defined standards. This paper proposes a benchmarking framework for evaluation of land model performances and, meanwhile, highlights major challenges at this infant stage of benchmark analysis. The framework includes (1) targeted aspects of model performance to be evaluated, (2) a set of benchmarks as defined references to test model performance, (3) metrics to measure and compare performance skills among models so as to identify model strengths and deficiencies, and (4) model improvement. Land models are required to simulate exchange of water, energy, carbon and sometimes other trace gases between the atmosphere and land surface, and should be evaluated for their simulations of biophysical processes, biogeochemical cycles, and vegetation dynamics in response to climate change across broad temporal and spatial scales. Thus, one major challenge is to select and define a limited number of benchmarks to effectively evaluate land model performance. The second challenge is to develop metrics of measuring mismatches between models and benchmarks. The metrics may include (1) a priori thresholds of acceptable model performance and (2) a scoring system to combine data–model mismatches for various processes at different temporal and spatial scales. The benchmark analyses should identify clues of weak model performance to guide future development, thus enabling improved predictions of future states of ecosystems and climate. The near-future research effort should be on development of a set of widely acceptable benchmarks that can be used to objectively, effectively, and reliably evaluate fundamental properties of land models to improve their prediction performance skills.

Forest models are instrumental for understanding and projecting the impact of climate change on forests. A considerable number of forest models have been developed in the last decades. However, few systematic and comprehensive model comparisons have been performed in Europe that combine an evaluation of modelled carbon and water fluxes and forest structure. We evaluate 13 widely used, state‐of‐the‐art, stand‐scale forest models against field measurements of forest structure and eddy‐covariance data of carbon and water fluxes over multiple decades across an environmental gradient at nine typical European forest stands. We test the models' performance in three dimensions:accuracy of local predictions(agreement of modelled and observed annual data),realism of environmental responses(agreement of modelled and observed responses of daily gross primary productivity to temperature, radiation and vapour pressure deficit) andgeneral applicability(proportion of European tree species covered). We find that multiple models are available that excel according to our three dimensions of model performance. For the accuracy of local predictions, variables related to forest structure have lower random and systematic errors than annual carbon and water flux variables. Moreover, the multi‐model ensemble mean provided overall more realistic daily productivity responses to environmental drivers across all sites than any single individual model. The general applicability of the models is high, as almost all models are currently able to cover Europe's common tree species. We show that forest models complement each other in their response to environmental drivers and that there are several cases in which individual models outperform the model ensemble. Our framework provides a first step to capturing essential differences between forest models that go beyond the most commonly used accuracy of predictions. Overall, this study provides a point of reference for future model work aimed at predicting climate impacts and supporting climate mitigation and adaptation measures in forests.

Wildfires, insect outbreaks, and storms cause large pulses of tree mortality. Climate change amplifies these forest disturbances, yet their future magnitude and extent remain uncertain. Here, we simulated future forest disturbance regimes at 100-meter resolution across Europe using a deep learning–based simulation framework. Our results show that forest disturbances will continue to increase throughout the 21st century, with disturbed areas more than doubling relative to the recent past under an unabated continuation of climate change. Wildfires are the main agent driving future disturbance change. Changing disturbances result in an increase in young forests, substantially altering Europe’s forest demography. Because of their profound implications for forest carbon storage and the habitat value of forest ecosystems, disturbances should be a priority of forest policy and management.

Stand age significantly influences the functioning of forest ecosystems by shaping structural and physiological plant traits, affecting water and carbon budgets. Forest age distribution is determined by the interplay of tree mortality and regeneration, influenced by both natural and anthropogenic disturbances. Unfortunately, human-driven alteration of tree age distribution presents an underexplored avenue for enhancing forest stability and resilience. In our study, we investigated how age impacts the stability and resilience of the forest carbon budget under both current and future climate conditions. We employed a state-of-the-science biogeochemical, biophysical, validated process-based model on historically managed forest stands, projecting their future as undisturbed systems, i.e., left at their natural evolution with no management interventions (i.e., forests are left to develop undisturbed). Such a model, forced by climate data from five Earth System Models under four representative climate scenarios and one baseline scenario to disentangle the effect of climate change, spanned several age classes as representative of the current European forests' context, for each stand. Our findings indicate that Net Primary Production (NPP) peaks in the young and middle-aged classes (16- to 50-year-old), aligning with longstanding ecological theories, regardless of the climate scenario. Under climate change, the beech forest exhibited an increase in NPP and maintained stability across all age classes, while resilience remained constant with rising atmospheric CO₂ and temperatures. However, NPP declined under climate change scenarios for the Norway spruce and Scots pine sites. In these coniferous forests, stability and resilience were more influenced. These results underscore the necessity of accounting for age class diversity -lacking in most, if not all, the current Global Vegetation Models - for reliable and robust assessments of the impacts of climate change on future forests' stability and resilience capacity. We, therefore, advocate for customized management strategies that enhance the adaptability of forests to changing climatic conditions, taking into account the diverse responses of different species and age groups to climate.

The carbon cycle component of the newly developed Earth system model of intermediate complexity CLIMBER-X is presented. The model represents the cycling of carbon through the atmosphere, vegetation, soils, seawater and marine sediments. Exchanges of carbon with geological reservoirs occur through sediment burial, rock weathering and volcanic degassing. The state-of-the-art HAMOCC6 model is employed to simulate ocean biogeochemistry and marine sediment processes. The land model PALADYN simulates the processes related to vegetation and soil carbon dynamics, including permafrost and peatlands. The dust cycle in the model allows for an interactive determination of the input of the micro-nutrient iron into the ocean. A rock weathering scheme is implemented in the model, with the weathering rate depending on lithology, runoff and soil temperature. CLIMBER-X includes a simple representation of the methane cycle, with explicitly modelled natural emissions from land and the assumption of a constant residence time of CH₄ in the atmosphere. Carbon isotopes ¹³C and ¹⁴C are tracked through all model compartments and provide a useful diagnostic for model–data comparison. A comprehensive evaluation of the model performance for the present day and the historical period shows that CLIMBER-X is capable of realistically reproducing the historical evolution of atmospheric CO₂ and CH₄ but also the spatial distribution of carbon on land and the 3D structure of biogeochemical ocean tracers. The analysis of model performance is complemented by an assessment of carbon cycle feedbacks and model sensitivities compared to state-of-the-art Coupled Model Intercomparison Project Phase 6 (CMIP6) models. Enabling an interactive carbon cycle in CLIMBER-X results in a relatively minor slow-down of model computational performance by ∼ 20 % compared to a throughput of ∼ 10 000 simulation years per day on a single node with 16 CPUs on a high-performance computer in a climate-only model set-up. CLIMBER-X is therefore well suited to investigating the feedbacks between climate and the carbon cycle on temporal scales ranging from decades to >100 000 years.