# A harmonized database of European forest simulations under climate change

Marc Grünig, Werner Rammer, Katharina Albrich, Frédéric André, Andrey L.D. Augustynczik, Friedrich Bohn, Meike Bouwman, Harald Bugmann, Alessio Collalti, Irina Cristal, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Marc Grünig, Werner Rammer, Katharina Albrich, Frédéric André, Andrey L.D. Augustynczik, et al. A harmonized database of European forest simulations under climate change. Data in Brief, 2024, 54, pp.110384. 10.1016/j.dib.2024.110384. hal-04771698

HAL Id: hal-04771698

https://hal.science/hal-04771698v1

Submitted on 7 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Data Article

A harmonized database of European forest
simulations under climate change

Check for
updates

Marc Grünig<sup>a,*</sup>, Werner Rammer<sup>a</sup>, Katharina Albrich<sup>b</sup>, Frédéric André<sup>c</sup>, Andrey L.D. Augustynczik<sup>d</sup>, Friedrich Bohn<sup>e</sup>, Meike Bouwman<sup>f</sup>, Harald Bugmann<sup>g</sup>, Alessio Collalti<sup>h,i</sup>, Irina Cristal<sup>j</sup>, Daniela Dalmonech<sup>h,i</sup>, Miquel De Caceres<sup>k</sup>, Francois De Coligny<sup>l</sup>, Laura Dobor<sup>m</sup>, Christina Dollinger<sup>a</sup>, David I. Forrester<sup>n</sup>, Jordi Garcia-Gonzalo<sup>j</sup>, José Ramón González<sup>j</sup>, Ulrike Hiltner<sup>g</sup>, Tomáš Hlásny<sup>m</sup>, Juha Honkaniemi<sup>b</sup>, Nica Huber<sup>g,w</sup>, Mathieu Jonard<sup>c</sup>, Anna Maria Jönsson<sup>o</sup>, Fredrik Lagergren<sup>o</sup>, Mats Nieberg<sup>p,q,r</sup>, Marco Mina<sup>s</sup>, Frits Mohren<sup>f</sup>, Christine Moos<sup>t</sup>, Xaxier Morin<sup>u</sup>, Bart Muys<sup>v</sup>, Mikko Peltoniemi<sup>b</sup>, Christopher PO Reyer<sup>p</sup>, Ilié Storms<sup>v</sup>, Dominik Thom<sup>a</sup>, Maude Toïgo<sup>u,x</sup>, Rupert Seidl<sup>a</sup>

<figure><img src="image_3.png" /></figure>

<footnote><sup>a</sup> TUM School of Life Sciences, Ecosystem Dynamics and Forest Management, Technical University of Munich, Hans-Carl-von-Carlowitz-Platz 2, 85354 Freising, Germany</footnote>

<footnote><sup>b</sup> Natural Resources Institute Finland, Forest Health and Biodiversity Group, Latokartanonkaari 9, 00790 Helsinki, Finland</footnote>

<footnote><sup>c</sup> Earth and Life Institute, Université catholique de Louvain, Croix du S, 1348 Ottignies-Louvain-la-Neuve, Belgium</footnote>

<footnote><sup>d</sup> International Institute for Applied Systems Analysis, Integrated Biosphere Futures Research Group, Schlossplatz 1, A-2361 Laxenburg, Austria</footnote>

<footnote><sup>e</sup> Helmholtz Centre for Environmental Research UFZ, Permoserstraße 15, 04318 Leipzig, Germany</footnote>

<footnote><sup>f</sup> Wageningen University & Research, Forest Ecology and Forest Management Group, Droevendaalsesteeg 3a, 6708 PB Wageningen, the Netherlands</footnote>

<footnote><sup>g</sup> ETH Zürich, Forest Ecology, Institute of Terrestrial Ecosystems, Universitätstrasse 16, 8006 Zürich, Switzerland</footnote>

<footnote><sup>h</sup> National Research Council of Italy (CNR-ISAFOM), Institute for Agriculture and Forestry Systems in the Mediterranean, Forest Modelling Lab., Via Madonna Alta 128, 06128 Perugia, Italy</footnote>

<footnote><sup>i</sup> National Biodiversity Future Center (NBFC), Piazza Marina, 61 90133 Palermo, Italy</footnote>

<footnote><sup>j</sup> Forest Science and Technology Center of Catalonia (CTFC), Crta. de St. Llorenç de Morunys, 25280 Solsona, Spain</footnote>

<footnote><sup>k</sup> CREAF, E08193 Bellaterra (Cerdanyola del Vallès), Catalonia, Spain</footnote>

<footnote><sup>l</sup> AMAP, INRAE-CIRAD-CNRS-IRD-Univers Montpellier, 34398 Montpellier cedex 5, France</footnote>

<footnote><sup>m</sup> Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, 165 21 Prague 6, Kamýcká 129, Czech Republic</footnote>

<footnote>* Corresponding author.</footnote>

<footnote>E-mail address: marc.gruenig@tum.de (M. Grünig).</footnote>

<footnote>Social media: 👄 @GruenigMarc (M. Grünig), 👄 @RammerWerner (W. Rammer), 👄 @RupertSeidl (R. Seidl)</footnote>
## ARTICLE INFO

**Article history:**

Received 1 March 2024

Revised 27 March 2024

Accepted 29 March 2024

Available online 3 April 2024

**Dataset link:** Data for: A harmonized database of European forest simulations under climate change (Original data)

**Keywords:**

Process-based models

Vegetation dynamics

Europe's forests

Forest development

Forest structure

Forest composition

Forest functioning

## ABSTRACT

Process-based forest models combine biological, physical, and chemical process understanding to simulate forest dynamics as an emergent property of the system. As such, they are valuable tools to investigate the effects of climate change on forest ecosystems. Specifically, they allow testing of hypotheses regarding long-term ecosystem dynamics and provide means to assess the impacts of climate scenarios on future forest development. As a consequence, numerous local-scale simulation studies have been conducted over the past decades to assess the impacts of climate change on forests. These studies apply the best available models tailored to local conditions, parameterized and evaluated by local experts. However, this treasure trove of knowledge on climate change responses remains underexplored to date, as a consistent and harmonized dataset of local model simulations is missing.

Here, our objectives were (i) to compile existing local simulations on forest development under climate change in Europe in a common database, (ii) to harmonize them to a common suite of output variables, and (iii) to provide a standardized vector of auxiliary environmental variables for each simulated location to aid subsequent investigations. Our dataset of European stand- and landscape-level forest simulations contains over 1.1 million simulation runs representing 135 million simulation years for more than 13,000 unique locations spread across Europe. The data were harmonized to consistently describe forest development in terms of stand structure (dominant height), composition (dominant species, admixed species), and functioning (leaf area index). Auxiliary variables provided include consistent daily climate information (temperature, precipitation, radiation, vapor pressure deficit) as well as information on local site conditions (soil depth, soil physical properties, soil water holding capacity, plant-available nitrogen). The present dataset facilitates analyses across models and locations, with the aim to better harness the valuable information contained in local simulations for large-scale policy support, and for fostering a deeper understanding of the effects of climate change on forest ecosystems in Europe.

© 2024 The Authors. Published by Elsevier Inc.

This is an open access article under the CC BY license

(http://creativecommons.org/licenses/by/4.0/)

**Article history:**

Received 1 March 2024

Revised 27 March 2024

Accepted 29 March 2024

Available online 3 April 2024

**Dataset link:** Data for: A harmonized database of European forest simulations under climate change (Original data)

**Keywords:**

Process-based models

Vegetation dynamics

Europe's forests

Forest development

Forest structure

Forest composition

Forest functioning

<footnote>ª CSIRO Environment, GPO Box 1700, ACT 2601, Australia</footnote>

<footnote>ᵇ Department of Physical Geography and Ecosystem Science, Lund University, Sölvegatan 12, 223 62 Lund, Sweden</footnote>

<footnote>ᵇ Potsdam Institute for Climate Impact Research (PIK), Member of the Leibniz Association, Telegrafenberg A 31, Potsdam, Germany</footnote>

<footnote>ᶜ European Forest Institute, Platz der Vereinten Nationen 7, 53113 Bonn, Germany</footnote>

<footnote>ᵈ Technische Universität Dresden, Chair of Forest Growth and Woody Biomass Production, Pienner Straße 8, 01737 Tharandt, Germany</footnote>

<footnote>ᵉ Institute for Alpine Environment, Eurac Research, Via Alessandro Volta, 13A, 39100 Bolzano, BZ, Italy</footnote>

<footnote>ᵉ² Bern University of Applied Sciences, BFH-HAFL, Länggasse 85, 3052 Zollikofen, Switzerland</footnote>

<footnote>ᵉ³ Université de Montpellier Université Paul-Valéry Montpellier - EPHE-IRD, CEFE UMR 5175, CNRS, 1919 Route de Mende, F-34293 Montpellier, France</footnote>

<footnote>ᵉ⁴ KU Leuven, Department of Earth and Environmental Sciences, Celestijnenlaan 200E, 3001 Leuven, Belgium</footnote>

<footnote>ᵉ⁵ Swiss Federal Research Institute WSL, Remote Sensing, Zürcherstrasse 111, CH-8903 Birmensdorf, Switzerland</footnote>

<footnote>ᵉ⁶ Université Bordeaux, Bordeaux Sciences Agro, INRAE, Biogeco, 69 route d'Arcachon, F-33612 Cestas, France</footnote>
Specifications Table

<table><tr><td>Subject</td><td>Environmental Sciences: Ecological modeling</td></tr><tr><td>Specific subject area</td><td>Harmonizing forest modeling simulations of climate change effects, process-based forest simulation models</td></tr><tr><td>Data format</td><td>Raw, harmonized</td></tr><tr><td>Type of data</td><td>Table, Database</td></tr><tr><td>Data collection</td><td>Data were contributed from the European forest modeling community. Each contributor uploaded simulation data files and a metadata file containing information on the design and drivers of the simulation as CSV files to an R Shiny application. Upon submission, the data were stored in a designated folder with a unique identifier assigned to each contributor. Contributions were required to follow specific criteria, including process-based simulations and annual output information on key vegetation development indicators (i.e. proportion of tree species, canopy height, and leaf area index). Other requirements were that simulation outputs were on the level of tree species, as well as in the absence of disturbances and management (or business-as-usual management).</td></tr><tr><td>Data source location</td><td>Technical University of Munich, TUM School of Life Sciences, Ecosystem Dynamics and Forest Management Group</td></tr><tr><td>Data accessibility</td><td>Repository name: Data for: A harmonized database of European forest simulations under climate change<br>Data identification number: 10.5281/zenodo.10730807<br>Direct URL to data: https://zenodo.org/records/10730807</td></tr></table>

# 1. Value of the Data

* Forest simulation model outputs from 17 different models were collected and harmonized. The dataset contains 1.1 million individual simulation runs over 135 million simulation years across 13,599 unique locations in Europe, covering large proportions of the climate and soil conditions in Europe's forests.

* The database contains standardized output variables across all models. Specifically, harmonized simulation outputs are available for canopy height (structure), leaf area index (LAI; functioning) and tree species proportions (composition) at annual time step. To provide harmonized layers of context information for simulation results, we collated daily climate data for historic climate (1981–2005) and a set of climate change scenarios (2006–2100) and consistent soil properties for all simulations.

* This is the first harmonized dataset of local forest model simulations at continental scale. The data collected here will support synthetic analyses of climate change impact on Europe's forests, and will facilitate comparative analyses across locations and models. Further, our dataset also helps to identify regions that remain underrepresented in model-based climate impact assessments and should thus be the focus of future studies.

# 2. Background

The objective was to collate projections on forest development under climate change derived from simulation models. Specifically, we compiled existing simulation data from previously conducted analyses using published models at the stand- to landscape-scale. Contributions to the dataset were made by several experts of the European forest modeling community, and all contributors are co-authors of this paper. Model outputs were compiled for three common state variables describing complementary aspects of forest ecosystems. By harmonizing the output variables across different models and adding standardized climate and soil data we created a novel, bottom-up dataset for broad-scale, multi-model assessments of climate change impacts of Europe's forests.
<caption>**Table 1**</caption>

Structure of the harmonized simulation database including standardized auxiliary data for all simulations. Discrete vegetation states were created by combining the three state variables and binning the respective continuous variables – see methods for details. Note that each row in Table 1 describes a column in the simulation database.

<table><thead><tr><th>Section</th><th>Column name</th><th>Description</th></tr></thead><tbody><tr><td rowspan="3">General</td><td>SourceID</td><td>ID to identify the source of the data (i.e. contributor).</td></tr><tr><td>SimulationID</td><td>Unique numeric identifier of the simulation created by the contributor.</td></tr><tr><td>Year</td><td>Year of the simulation starting with 1 or the calendar year (e.g., 2000).</td></tr><tr><td rowspan="2">Vegetation</td><td>Discrete vegetation state</td><td>Discrete state derived by combining the three state variables species composition, LAI class and dominant canopy height class (e.g. PIAB_3_20_22; see methods for details).</td></tr><tr><td>Dominant height</td><td>In meters, harmonized (calculated from min, mean and max heights, details on the calculations are shown in the methods).</td></tr><tr><td rowspan="6">Soil</td><td>LAI</td><td>Leaf Area Index (one-sided or projected) in m²/m².</td></tr><tr><td>WHC</td><td>Water holding capacity of the site (mm).</td></tr><tr><td>TextureSand</td><td>% sand content of the soil.</td></tr><tr><td>TextureSilt</td><td>% silt content of the soil.</td></tr><tr><td>TextureClay</td><td>% clay content of the soil.</td></tr><tr><td>SoilDepth</td><td>Depth of the plant-accessible soil (mm) without rocks (&gt; 2 mm diameter).</td></tr><tr><td rowspan="7">Climate</td><td>AvailableNitrogen</td><td>Plant available nitrogen (kg/ha/year).</td></tr><tr><td>Scenario</td><td>Combination of GCM and RCP from which daily data was obtained.</td></tr><tr><td>Temperature</td><td>Columns "tas_1" to "tas_365" with daily mean temperature [°C].</td></tr><tr><td>Precipitation</td><td>Columns "prec_1" to "prec_365" with daily precipitation [mm].</td></tr><tr><td>Radiation</td><td>Columns "rad_1" to "rad_365" with daily radiation [W/m²].</td></tr><tr><td>Vapor pressure deficit (VPD)</td><td>Columns "vpd_1" to "vpd_365" with daily vapor pressure deficit [kPa].</td></tr></tbody></table>

## 3. Data Description

The data are collected and stored in SQLite format. SQLite is a widely used open-source database format and can be accessed from all major data analysis platforms. One SQLite database contains the raw simulation outputs and a metadata table of all simulations including information about locations and harmonized soil conditions for those locations. This data follows the structure described in detail in the supplementary information (Tables S1 and S2). Simulation outputs with harmonized climate data are stored in one SQLite database per climate scenario. Tables in those databases follow the structure shown in Table 1. Further, a metadata table of all simulations, including information about locations and soil conditions for those locations is provided (Table S2).

The database contains 1,117,453 simulation runs that together contain 135,375,583 simulation years. Simulations cover 13,599 unique locations across Europe and represent 92 tree species. Simulation data were provided from 19 research groups, using 17 different forest models. All simulations were created with locally tested and evaluated models that are well-documented and published in the peer-reviewed literature (Table 2). Note that as models are further developed over time, model versions used for the simulations may vary from the cited references in some cases. While all models contributed to the coverage of climate and soil conditions across Europe, their individual contributions varied in terms of geographic range and number of simulations provided. Likewise, for some models more data were available than for others. While the majority of simulations were run with iLand, MEDFATE simulations covered the largest climate and soil gradient. The simulations in the database consist of 90.5% climate change runs and
<caption>**Table 2**</caption>

Models, number of simulations per model and their coverage of current climate and soil space of Europe. Percentage of climate and soil space refers to the area covered by the climate and soil space in which the model simulations are located. For this, we stratified the climate and soil space and checked which of the classes are covered by the simulations of each model. We then calculated the area covered by the classes that are represented by each model. Further details are described in the Experimental Design, Materials and Methods section. The Model type column distinguishes between stand-level (S) models and landscape-level (L) models.

<table><thead><tr><th rowspan="2">MODEL</th><th colspan="4">MODEL OVERVIEW</th><th rowspan="2">Model type</th><th rowspan="2">Model reference</th></tr><tr><th>Simulations</th><th>% of all</th><th>% clim</th><th>% soil</th></tr></thead><tbody><tr><td>ILAND</td><td>821,979</td><td>73.2</td><td>21.4</td><td>43.1</td><td>L</td><td>Seidl et al. 2012 [2]</td></tr><tr><td>4C</td><td>250,979</td><td>22.4</td><td>18.3</td><td>54.3</td><td>S</td><td>Lasch-Born et al., 2020 [3]</td></tr><tr><td>MEDFATE</td><td>22,464</td><td>2.0</td><td>40.8</td><td>64.5</td><td>S</td><td>De Cáceres et al., 2021 [4]</td></tr><tr><td>FORCLIM</td><td>9210</td><td>0.8</td><td>10.3</td><td>46.7</td><td>S</td><td>Bugmann, 1996 [5]; Huber et al., 2021 [6]</td></tr><tr><td>LPJ-GUESS2.1</td><td>6156</td><td>0.6</td><td>15.5</td><td>46.6</td><td>L</td><td>Smith et al., 2008 [7]</td></tr><tr><td>FORCEEPS</td><td>5040</td><td>0.5</td><td>4.8</td><td>6.7</td><td>S</td><td>Morin et al., 2021 [8]</td></tr><tr><td>TREEMIG</td><td>2820</td><td>0.3</td><td>0.6</td><td>0.3</td><td>L</td><td>Lischke et al., 2006 [9]</td></tr><tr><td>3PGN-BW</td><td>1428</td><td>0.1</td><td>40.0</td><td>59.8</td><td>S</td><td>Augustynczik and Yousefpour, 2021 [10]</td></tr><tr><td>FORMIND</td><td>1008</td><td>0.1</td><td>24.2</td><td>6.0</td><td>S</td><td>Fischer et al., 2016 [11]</td></tr><tr><td>3D-CMCC-FEM</td><td>367</td><td>&lt; 0.1</td><td>13.1</td><td>7.3</td><td>S</td><td>Collalti et al., 2018 [12]; Dalmonech et al., 2022 [13]</td></tr><tr><td>GOTILWA+</td><td>336</td><td>&lt; 0.1</td><td>10.1</td><td>7.2</td><td>S</td><td>Nadal-Sala et al., 2013 [14]</td></tr><tr><td>SORTIE-ND</td><td>278</td><td>&lt; 0.1</td><td>21.3</td><td>22.0</td><td>S</td><td>Canham et al., 2005 [15]</td></tr><tr><td>PREBAS</td><td>252</td><td>&lt; 0.1</td><td>6.6</td><td>6.2</td><td>S</td><td>Minunno et al., 2019 [16]</td></tr><tr><td>3PGMIX</td><td>173</td><td>&lt; 0.1</td><td>10.1</td><td>34.0</td><td>S</td><td>Forrester and Tang, 2016 [17]</td></tr><tr><td>LANDSCAPE</td><td>38</td><td>&lt; 0.1</td><td>11.3</td><td>7.2</td><td>L</td><td>Haas et al., 2013 [18]</td></tr><tr><td>DNDC</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HETEROFOR</td><td>12</td><td>&lt; 0.1</td><td>6.7</td><td>18.6</td><td>S</td><td>Jonard et al., 2020 [19]</td></tr></tbody></table>

9.5% of the simulations were run under baseline or observed climate conditions. The database contains simulations from 12 stand-level models and four landscape-level models (Table 2, see Bugmann & Seidl, 2022 [1] for a review on modeling approaches). For balance between stand- and landscape-level simulations, a subset of 1-ha stands from the full landscape was sampled and used as individual simulation runs.

In geographic space, the simulation runs cluster in Central Europe, Spain (Catalonia), Finland and Sweden, i.e. areas that were analyzed particularly intensively in previous modeling studies (Fig. 1). To evaluate the proportion of the geographic area of Europe that was covered by the climate and soil conditions represented in the dataset, the climate and soil space (using the variables shown in Fig. 2) of the entire continent were stratified and strata in which simulations were located were obtained. As some of the strata cover more geographic area (i.e. climate and soil conditions that occur more often), a simulation in that strata could cover a larger percentage of climate and soil space than other simulations. The climate space covered by all simulation data spans 79% of Europe's geographic area and the covered soil space spanned 75.4% of the geographic area (Fig. 2). Areas not covered are mostly unforested regions in highly continental parts of eastern Europe (i.e. parts of Ukraine, western Romania), very warm regions in the Mediterranean (i.e. southern Spain, parts of Greece and Italy), and wet and very oceanic regions including large parts of the British Isles and southern Norway. For soil conditions, we find that mainly soils with very coarse or very fine texture, low nitrogen availability and low water holding capacity are not covered by simulations.

The database contains simulations for 92 species. While most species are represented in less than 1000 simulations (median 728), *Fagus sylvatica*, *Picea abies*, *Larix decidua*, and *Pinus sylvestris* are the most prevalent, each occurring in over 400,000 simulations (Fig. 3). Furthermore, the dataset contains simulations without forest management (39.2%) and simulations implementing common practices (60.8%).
<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Locations of the simulations in geographical space. Simulations from different models are displayed with different colors and symbols. Note that the projection of the map is Lambert Azimuthal Equal-Area (LAEA), the true north can be identified by following the longitudinal lines plotted in grey at 0° and 20° East.</figcaption></figure>

# 4. Experimental Design, Materials and Methods

## 4.1. Data collection

A central database of forest simulations was created from existing model outputs from previous simulation studies, collating information from a variety of stand- to landscape-scale forest simulation models across Europe. The focus lay on models that were locally evaluated, thus encapsulating the best available bottom-up understanding of quantitative forest ecosystem dynamics in Europe. Furthermore, the dataset was restricted to results from process-based models, as these are expected to be more robust under changing environmental conditions compared to purely empirical models [20]. Simulation outputs were provided by modeling groups in two files, a simulation data file and a metadata file containing information on the design and drivers of the simulation (Table S1 & S2). Forest simulations under climate change conditions were of particular interest, but simulation runs under baseline climate conditions were also included. The minimum requirements for contribution to the dataset were that the model conducting the simulation was process-based, and that simulation outputs provide annual information on basic indicators of vegetation development (i.e. proportion of tree species, canopy height, and leaf area index). While some models provide a broader range of output indicators, these three variables were chosen as least common denominator for describing forest composition, structure, and functioning. Furthermore, the simulations had to provide information at the level of individ-
a

<figure><img src="image_5.png" /></figure>

c

<figure><img src="image_6.png" /></figure>

e

<figure><img src="image_7.png" /></figure>

b

<figure><img src="image_8.png" /></figure>

d

<figure><img src="image_9.png" /></figure>

f

<figure><img src="image_10.png" /><figcaption>**Fig. 2.** Climate and soil space covered by simulations. Climate space was stratified into unique combinations of mean annual temperature (MAT; ten classes; panel a), annual precipitation (ANP; 7 classes; panel c) and temperature seasonality (Seas; seven classes; panel e). Soil space was stratified to unique combinations of plant-available nitrogen (N; five classes, but 0–20 kg/ha/year is not covered; panel b), soil texture (classes from coarse (=1) to very fine (=5) according to the European Soil Database classification scheme for soil texture; panel d) and water holding capacity (WHC; five classes; panel f). Areas not covered by the data in this database are hatched. For more details and data sources see section 3.2.1 Climate and 3.2.2 Soil.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>Fig. 3. Number of simulations in which the 20 most prevalent tree species in the dataset occur.</figcaption></figure>

ual tree species, models simulating plant functional types were not included. With regard to the initial conditions of the simulation, both generic initial conditions (such as simulations starting from bare ground) and runs initialized with the current state of the vegetation were included. Further, simulation runs were conducted either in the absence of management interventions or assumed business-as-usual management for the area under study. Natural disturbances were not considered in the simulations. Simulation data were compiled and analyzed at stand level, hence information from landscape-scale simulation models were considered as unique data vectors at the level of simulated stands (i.e., areas of homogeneous climate and soil conditions) within a landscape. The temporal extent of simulations was variable, and both historical and future time series were included.

## 4.2. Data harmonization

Simulation data went through initial checks to ensure that the metadata for each simulation were complete and IDs in the metadata table (Table S2) were matching the simulation data (Table S1). To harmonize simulation outputs across the different models and simulation studies we used discrete vegetation states as described by Rammer & Seidl (2019) [21]. Discrete
vegetation states condense the complexity of forest vegetation by describing the structure, composition, and function of vegetation as distinct states. To that end, continuous variables were discretized, and a state was derived for every simulation year in the output of each model. With regard to stand composition, a species was categorized as dominant if it held more than 66% of a stand's basal area. Admixed species were explicitly considered if their share was ≥20%. Based on these categories, a unique string was derived to describe species composition by concatenating the first two characters of the genus name and the first two characters of the species name (e.g. *Pinus sylvestris* to PISY), and combining all dominant and admixed species occurring in a stand. Letters for dominant species were capitalized, while admixed species were in lowercase letters (e.g. PISYfasy means *Pinus sylvestris* is the dominant species and has a share of more than 66%, with *Fagus sylvatica* as an admixed species with more than 20% of the basal area). To categorize ecosystem functioning, LAI (as a key indicator for the exchange of carbon, water, and energy of the system with the atmosphere) was grouped into three classes: class 1 for sparse (LAI < 2), class 2 for moderate (2 ≤ LAI ≤ 4), and class 3 for dense forest canopies (LAI > 4). For structure, the dominant height of the canopy was employed, as the vertical utilization of space is a key element of forest structure. Canopy height can further indicate the development stage of forest stands. The height information from different models (mean, minimum, maximum) was harmonized to dominant stand height. For stands dominated by *Pinus sylvestris*, *Abies alba*, *Larix decidua*, *Picea abies*, *Fagus sylvatica* and *Quercus robur*, allometric factors calibrated by Kahn (1994) [22] were applied to describe the relationship of stand dominant height to maximum tree height. For other species, a statistical model based on yield table data was derived to estimate dominant height from mean and maximum tree height [22]. Specifically, for 14 species, a linear mixed model with dominant height as dependent variable, mean and maximum height as predictor variable, and random slopes for the different species was calibrated (conditional R² = 0.999; marginal R² = 0.996). To make the analysis robust to individual outliers, dominant height was limited to lie between 0.8 times maximum height (lower limit) and maximum height (upper limit). Some models only provided information for mean height (but not maximum height), or for minimum and maximum height only (but not mean height). In the latter case, mean height was calculated as the arithmetic mean of minimum and maximum height. For simulations for which only mean height was available, maximum height was estimated using a linear mixed model. This model was calibrated on a random subsample of the dataset of all simulation data across all models that contained both maximum and mean height. In the linear mixed model, maximum height was used as the dependent variable, mean height as predictor variable and random effects for species were included (conditional R² = 0.54; marginal R² = 0.25). Maximum height was then predicted from mean heights with this model and subsequently calculated dominant height from mean and maximum height as described before. Subsequently, continuous information on dominant height was grouped into 2-m bins. Finally, the individual states for forest composition, structure, and function were combined into a unique string describing the state of the vegetation (e.g., PIABfasy_3_20_22 representing a stand dominated by *Picea abies* (PIAB) with admixed *Fagus sylvatica* (fasy) that has a dense forest canopy (LAI class 3) and a dominant height of between 20 and 22 m). This harmonization and discretization of the underlying simulation data resulted in a total of 18,598 distinct vegetation states being recorded in our database.

## 4.3. Auxiliary data

In addition to the harmonized simulation data and the respective metadata, a common data vector of auxiliary data was compiled. This vector contained standardized climate and soil data for each simulation, in order to facilitate the analysis and interpretation of the dataset.

### 4.3.1. Climate

To derive common and coherent climate data for all simulation runs included in the database, the simulation-specific climate time series (which was restricted to annual values for tempera-
ture and precipitation) provided by the modeling groups was matched with climate data from a common climate dataset. This approach allowed to address problems of different data resolution and representation in the data used for individual simulations (e.g., some models operate with monthly climate data while others use daily climate data, some models use maximum and minimum temperature as drivers while others use mean temperature, some models include climatic drivers beyond temperature and precipitation such as vapor pressure deficit (VPD) while others do not, some simulation runs were driven by detailed downscaled climate data while others used coarser-resolution climate information as input). As common climate database EURO-CORDEX climate data was harnessed for historical conditions as well as three RCP scenarios (RCP2.6, RCP4.5 and RCP8.5), each simulated with three global circulation models (MPI-MMPI-ESM-LR, ICHEC-EC-EARTH, NCC-NorESM1-M, all downscaled with the SMHI-RCA4 RCM), resulting in 12 climate scenarios. The climate data was obtained in a 0.11° x 0.11° spatial resolution and daily temporal resolution from the Copernicus Climate Data Store. (https://cds.climate.copernicus.eu/cdsapp#!/dataset/projections-cordex-domains-single-levels?tab=overview). GCMs were selected to cover a broad gradient of temperature and precipitation conditions [23] representative for the current and future climatic conditions in Europe.

For the harmonized dataset, the goal was to obtain daily climate data from EURO-CORDEX for all simulations. To achieve this, the best matching EURO-CORDEX trajectory was assigned to each simulation and for each simulation year the daily data from the best matching year from that scenario was extracted. In more detail, information on the type of climate trajectory used for the simulation run (baseline run without climate change, SRES, or RCP scenario family) was obtained from the metadata of each simulation. Baseline conditions refer to historical time series, while SRES and RCP scenarios refer to simulations under future climate scenarios. The best matching EURO-CORDEX scenario was assigned to each simulation based on a comparison between slopes (i.e. temporal change) of temperature and precipitation used in the simulations with those calculated from EURO-CORDEX scenarios for the simulation location. The difference of the slope for temperature and precipitation was obtained and the overall difference was calculated as $ΔT + 0.1* ΔP$ to choose the trajectory with the smallest overall difference. Whenever simulation metadata contained specific information on the GCM and/or RCP used for the climate forcing used in the simulation, this information was harnessed to limit selection options from the EURO-CORDEX database. For instance, if the climate trajectory of a simulation was run with an RCP2.6 scenario, the temperature and precipitation slopes were only compared with the three RCP2.6 scenarios contained in our EURO-CORDEX selection. The best matching trajectory was then adjusted to the mean temperature and precipitation level of the simulation. Specifically, the difference between mean temperatures of the simulation and climate scenario trajectory and the multiplicative difference for precipitation was added. This was necessary to make sure that the adjusted time series represented the climate used in the simulations, as they may differ due to the relatively coarse resolution of the gridded climate data. Next, a time series from the adjusted scenario data was constructed, matching the individual years to each year in the simulation. For this, an index combining the additive difference between mean annual temperatures and the multiplicative difference of annual precipitation (again with $ΔT + 0.1* ΔP$) was calculated for each pair of simulation data year and climate scenario year. To prevent multiple occurrences of the same year, one of the three best matching years was randomly sampled with replacement. Finally, the thus constructed daily climate time series data for all variables (temperature, precipitation, radiation and VPD) was stored along the simulation it represents.

### 4.3.2. Soil

Mirroring the approach taken with climate data, the simulation metadata provided by the modelling groups was facilitated to derive a consistent and quantitative soil data set across all simulations. While some simulation metadata contained exact numbers on all relevant soil variables considered (i.e. soil depth, soil texture (sand, silt, clay percentage), water holding capacity (WHC) and available nitrogen), others provided more descriptive values for water and nutrient conditions, such as soil water or fertility ratings (see Table S3). To complete missing data and convert descriptive ratings to quantitative soil characteristics pan-European datasets on soil in-
formation were leveraged. For soil depth, soil texture and WHC gridded data (1 km resolution) from the European Soil Data Center (ESDAC [24]) was obtained. Soil fertility was approximated by means of plant-available nitrogen, considering that nitrogen is the most important macro-nutrient for forests in Europe at the continental scale. Since data on plant-available nitrogen (i.e., the annual flux of mineralized nitrogen) is not available for large scales, we used a sta-tistical approach to estimate mineralization rates based on fertility and climate data and com-bined them with data on nitrogen stocks. Specifically, the SQ1 nutrient availability map from the harmonized European soil database [25] was reclassified from the four classes of soil fertil-ity considered therein to values of 100, 65, 40 and 20 kg/ha/year of plant-available nitrogen. In a next step, from the Soilgrids dataset [26] the nitrogen pool (in kg/m²) across Europe based on relative nitrogen content (g/kg) and bulk density (kg/m³) for soil layers up to 30 cm soil depth were obtained. The reclassified soil fertility map was divided by the nitrogen pool layer to ap-proximate a coarse pseudo-mineralization rate. This pseudo-mineralization rate was then used as a dependent variable to calibrate a generalized linear model (GLM) with mean annual tem-perature, annual precipitation, temperature seasonality (obtained from CHELSA [27]; we used CHELSA data instead of EURO-CORDEX to create a layer of 1 km spatial resolution, matching the other soil layers) and soil pH (from Soilgrids) as predictor variables (D² = 0.64). The calibrated model was used to project pseudo-mineralization rates for continental Europe at 1 × 1 km res-olution providing more consistent estimates based on local soil and climate conditions. Finally, the nitrogen pool was multiplied with the modeled pseudo-mineralization rate to get an ap-proximation of plant-available nitrogen.

## 4.4. Coverage of simulated data

The coverage of the current climate and soil space of Europe's forests by simulation data compiled here cover was investigated using the harmonized data (see Fig. 2). The historical cli-mate was categorized into ten stratified bins of mean annual temperature (< -2 °C, -2 - 2 °C, 2-4 °C, 4 - 6 °C, 6 - 8 °C, 8 - 10 °C, 10 - 12 °C, 12 - 14 °C, 14 - 18 °C, > 18 °C), seven stratified bins of annual precipitation sum (< 400 mm, 400 - 550 mm, 550 - 750 mm, 750 - 1000 mm, 1000 - 1400 mm, 1400 - 2000 mm, > 2000 mm) and seven stratified bins of temper-ature seasonality (calculated as standard deviation of monthly mean temperature* 100, binned to ( 200, 200 - 400, 400 - 500, 500 - 600, 600 - 700, 700 - 900, > 900). Simulations were assigned to their bins based on the climate grid cell (0.11°) of the simulation's location, and the area that is covered by the occupied bins was calculated. This approach likely underestimates the true climatic coverage of the simulations, as individual simulations (both stand-level and landscape-level) often encompass climatic gradients within a single cell. Soil conditions were categorized into five unique combinations of stratified values along the dimensions of soil tex-ture (six classes from fine to coarse calculated with sand, silt and clay content), water holding capacity (<50 mm, 50 - 100 mm, 100 - 150 mm, 150 - 180 mm, >180 mm) and soil fertility (plant-available nitrogen of < 20, 20 - 40, 40 - 65, 65 - 100, >100 kg/ha/year). Again, each simulated location was located in the three-dimensional space of soil conditions to assess how well the compiled simulation data represented the soil conditions of Europe's forests.

## Limitations

The majority of simulations were sourced from iLand model outputs (73.2%, as shown in Table 2), totaling over 800,000 entries stored in large tables, which were unwieldy to manage. To mitigate this, we divided the simulation data into smaller chunks with different unique identi-fiers. Furthermore, the imbalance in the number of simulations from the various models creates a bias that needs consideration during database utilization for analytical purposes.
## Ethics Statement

We confirm to have read the ethical requirements for publication in Data in Brief and confirm that the current work does not involve human subjects, animal experiments, or any data collected from social media platforms.

## Declaration of Generative AI and AI-Assisted Technologies in the Writing Process

During the preparation of this work the authors used ChatGPT in order to improve readability. After using this tool, the authors reviewed and edited the content as needed and take full responsibility for the content of the publication.

## Data Availability

Data for: A harmonized database of European forest simulations under climate change
(Original data) (Zenodo)

## CRediT Author Statement

**Marc Grünig:** Conceptualization, Methodology, Software, Formal analysis, Data curation, Writing – original draft, Writing – review & editing, Visualization; **Werner Rammer:** Conceptualization, Methodology, Software, Formal analysis, Data curation, Writing – original draft, Writing – review & editing, Funding acquisition; **Katharina Albrich:** Data curation, Writing – review & editing; **Frédéric André:** Data curation, Writing – review & editing; **Andrey L.D. Augustynczik:** Data curation, Writing – review & editing; **Friedrich Bohn:** Data curation, Writing – review & editing; **Meike Bouwman:** Data curation, Writing – review & editing; **Harald Bugmann:** Data curation, Writing – review & editing; **Alessio Collalti:** Data curation, Writing – review & editing; **Irina Cristal:** Data curation, Writing – review & editing; **Daniela Dalmonech:** Data curation, Writing – review & editing; **Miquel De Caceres:** Data curation, Writing – review & editing; **Francois De Coligny:** Data curation, Writing – review & editing; **Laura Dobor:** Data curation, Writing – review & editing; **Christina Dollinger:** Data curation, Writing – review & editing; **David I. Forrester:** Data curation, Writing – review & editing; **Jordi Garcia-Gonzalo:** Data curation, Writing – review & editing; **José Ramón González:** Data curation, Writing – review & editing; **Ulrike Hiltner:** Data curation, Writing – review & editing; **Tomáš Hlásny:** Data curation, Writing – review & editing; **Juha Honkaniemi:** Data curation, Writing – review & editing; **Nica Huber:** Data curation, Writing – review & editing; **Mathieu Jonard:** Data curation, Writing – review & editing; **Anna Maria Jönsson:** Data curation, Writing – review & editing; **Fredrik Lagergren:** Data curation, Writing – review & editing; **Mats Nieberg:** Data curation, Writing – review & editing; **Marco Mina:** Data curation, Writing – review & editing; **Frits Mohren:** Data curation, Writing – review & editing; **Christine Moos:** Data curation, Writing – review & editing; **Xavier Morin:** Data curation, Writing – review & editing; **Bart Muys:** Data curation, Writing – review & editing; **Mikko Peltoniemi:** Data curation, Writing – review & editing; **Christopher PO Reyer:** Data curation, Writing – review & editing; **Ilié Storms:** Data curation, Writing – review & editing; **Dominik Thom:** Data curation, Writing – review & editing; **Maude Toigo:** Data curation, Writing – review & editing; **Rupert Seidl:** Conceptualization, Methodology, Writing – original draft, Writing – review & editing, Funding acquisition.

## Acknowledgments

We acknowledge funding received from the European Union's Horizon 2020 research and innovation program under grant agreement no. 101000574 (RESONATE: Resilient forest value
chains – enhancing resilience through natural and socio-economic responses). We further acknowledge support from the European Research Council under the European Union's Horizon 2020 research and innovation program (Grant Agreement 101001905).

## Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

## Supplementary Materials

Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.dib.2024.110384.

## References

[1] H. Bugmann, R. Seidl, The evolution, complexity and diversity of models of long-term forest dynamics, J. Ecol. 110 (10) (2022) 2288-2307.

[2] R. Seidl, W. Rammer, R.M. Scheller, T.A. Spies, An individual-based process model to simulate landscape-scale forest ecosystem dynamics, Ecol. Modell. 231 (2012) 87-100.

[3] P. Lasch-Born, F. Suckow, C.P. Reyer, M. Gutsch, C. Kollas, F.W. Badeck, H.K. Bugmann, R. Grote, C. Fürstenau, M. Lindner, J. Schaber, Description and evaluation of the process-based forest model 4C v2. 2 at four European forest sites, Geosci. Model Dev. 13 (11) (2020) 5311-5343.

[4] M. De Cáceres, M. Mencuccini, N. Martin-StPaul, J.M. Limousin, L. Coll, R. Poyatos, A. Cabon, V. Granda, A. Forner, F. Valladares, J. Martínez-Vilalta, Unravelling the effect of species mixing on water use and drought stress in Mediterranean forests: a modelling approach, Agric. For. Meteorol. 296 (2021) 108233.

[5] H.K. Bugmann, A simplified forest model to study species composition along climate gradients, Ecology 77 (7) (1996) 2055-2074.

[6] N. Huber, H. Bugmann, M. Cailleret, N. Bircher, V. Lafond, Stand-scale climate change impacts on forests over large areas: transient responses and projection uncertainties, Ecol. Appl. 31 (4) (2021) e02313.

[7] B. Smith, I.C. Prentice, M.T. Sykes, Representation of vegetation dynamics in the modelling of terrestrial ecosystems: comparing two contrasting approaches within European climate space, Glob. Ecol. Biogeogr. (2001) 621-637.

[8] X. Morin, H. Bugmann, F. De Coligny, N. Martin-StPaul, M. Cailleret, J.M. Limousin, J.M. Ourcival, B. Prevosto, G. Simioni, M. Toigo, M. Vennetier, E. Catteau, J. Guillemot, Beyond forest succession: a gap model to study ecosystem functioning and tree community composition under climate change, Funct. Ecol. 35 (4) (2021) 955-975.

[9] H. Lischke, N.E. Zimmermann, J. Bolliger, S. Rickebusch, T.J. Löffler, TreeMig: a forest-landscape model for simulating spatio-temporal patterns from stand to landscape scale, Ecol. Modell. 199 (4) (2006) 409-420.

[10] A.L.D. Augustyniczik, R. Yousefpour, Assessing the synergistic value of ecosystem services in European beech forests, Ecosyst. Serv. 49 (2021) 101264.

[11] R. Fischer, F. Bohn, M.D. de Paula, C. Dislich, J. Groeneveld, A.G. Gutiérrez, ... A. Huth, Lessons learned from applying a forest gap model to understand ecosystem and carbon dynamics of complex tropical forests, Ecol. Modell. 326 (2016) 124-133.

[12] A. Collalti, P.E. Thornton, A. Cescatti, A. Rita, M. Borghetti, A. Nolè, G. Matteucci, The sensitivity of the forest carbon budget shifts across processes along with stand development and climate change, Ecol. Appl. 29 (2) (2019) e01837.

[13] D. Dalmonech, G. Marano, J.S. Amthor, A. Cescatti, M. Lindner, C. Trotta, A. Collalti, Feasibility of enhancing carbon sequestration and stock capacity in temperate and boreal European forests via changes to management regimes, Agric. For. Meteorol. 327 (2022) 109203.

[14] D. Nadal-Sala, S. Sabaté, C. Gracia, GOTILWA+: un modelo de procesos que evalúa efectos del cambio climático en los bosques y explora alternativas de gestión para su mitigación, Ecosistemas 22 (3) (2013) 29-36.

[15] C.D. Canham, L.E. Murphy, M.J. Papaik, SORTIE-ND: software for spatially-explicit simulation of forest dynamics, Institute of Ecosystem Studies, Millbrook, NY, 2005.

[16] F. Minunno, M. Peltoniemi, S. Härkönen, T. Kalliokoski, H. Makinen, A. Mäkelä, Bayesian calibration of a carbon balance model PREBAS using data from permanent growth experiments and national forest inventory, For. Ecol. Manag. 440 (2019) 208-257.

[17] D.I. Forrester, X. Tang, Analysing the spatial and temporal dynamics of species interactions in mixed-species forests and the effects of stand density using the 3-PG model, Ecol. Modell. 319 (2016) 233-254.

[18] E. Haas, S. Klatt, A. Fröhlich, P. Kraft, C. Werner, R. Kiese, ... K. Butterbach-Bahl, LandscapeDNDC: a process model for simulation of biosphere-atmosphere-hydrosphere exchange processes at site and regional scale, Landsc. Ecol. 28 (4) (2013) 615-636.

[19] M. Jonard, F. André, F. De Coligny, L. De Wergifosse, N. Beudez, H. Davi, ... C. Vincke, HETEROFOR 1.0: a spatially explicit model for exploring the response of structurally complex forests to uncertain future conditions-Part 1: carbon fluxes and tree dimensional growth, Geosci. Model Dev. 13 (3) (2020) 905-935.
[20] E.J. Gustafson, When relationships estimated in the past cannot be used to predict the future: using mechanistic models to predict landscape ecological dynamics in a changing world, Landsc. Ecol. 28 (2013) 1429-1437.

[21] W. Rammer, R. Seidl, A scalable model of vegetation transitions using deep neural networks, Methods Ecol. Evol. 10 (6) (2019) 879-890.

[22] M. Kahn, Modellierung der Höhenentwicklung ausgewählter Baumarten in Abhängigkeit bom Standort, Forstwissenschaftliche Fakultet der Universität der Universität München und Bayerische forstliche Versuchs-und Forschungsanstalt (1994).

[23] K.M. Parding, A. Dobler, C.F. McSweeney, O.A. Landgren, R. Benestad, H.B. Erlandsen, ... H. Loukos, GCMeval-An interactive tool for evaluation and selection of climate model ensembles, Clim. Serv. 18 (2020) 100167.

[24] R. Hiederer, Mapping Soil Properties for Europe - Spatial Representation of Soil Database Attributes, Luxembourg: Publications Office of the European Union -2013 -47pp. - EUR26082EN Scientific and Technical Research series, 2013 ISSN 1831-9424, doi:10.2788/94128.

[25] G. Fischer, F. Nachtergaele, S. Prieler, H.T. van Velthuizen, L. Verelst, D. Wiberg, Global Agro-ecological Zones Assessment for Agriculture (GAEZ 2008), 2008.

[26] T. Hengl, J. Mendes de Jesus, G.B. Heuvelink, M. Ruiperez Gonzalez, M. Kilibarda, A. Blagotić, ... B Kempen, Soil-Grids250m: global gridded soil information based on machine learning, PLoS One 12 (2) (2017) e0169748.

[27] D.N. Karger, O. Conrad, J. Böhner, T. Kawohl, H. Kreft, R.W. Soria-Auza, ... M. Kessler, Climatologies at high resolution for the earth's land surface areas, Sci. Data 4 (1) (2017) 1-20.