DOI: 10.6084/m9.figshare.30696074.v1
title: <b>Scots pine management under climate change: a dual-model simulation experiment</b>

description: <b>Simulations</b><b>Forest management projections of </b><b><i>P</i></b><b><i>inus sylvestris</i></b><b> virtual stands across xeric, mesic, and humid climatic conditions in Catalonia, NE Spain, were conducted under three climate change scenarios (current climate RCP 0.0, RCP 4.5, and RCP 8.5) and three management alternatives (low-intensity/high-frequency thinning LIT, high-intensity/low-frequency thinning HIT, and no management NM). Simulations were run using two forest dynamics models: the process-based GOTILWA+ (Gracia et al., 1999) and the gap model SORTIE-ND (Canham et al., 2005).</b>The simulation outputs were compiled into three files:<b>A raw RDS file</b> containing the complete set of simulation outputs, including all <b>10 repetitions for each of the 54 simulation cases</b>, capturing stochasticity in stand dynamics.<b>An aggregated dataset</b> containing <b>averaged values across the 10 repetitions, CSV</b><b>Timber production metrics</b>, calculated for each <b>harvesting interval, CSV</b>The three datasets include identifiers for <b>model</b>, <b>climatic condition</b>, <b>climate scenario</b>, <b>management alternative</b>, <b>simulation year (or harvest interval)</b>, and stand-level attributes such as basal area, DBH, timber production metrics.<b>Fields in </b><code><strong>sim_means</strong></code>The dataset <code><strong>sim_means</strong></code> contains the averaged values across the 10 repetitions for each simulation case. Each row corresponds to a unique combination of simulator, site, management alternative, climate scenario, and simulation year. The fields included are:<b>simulator</b> – Forest dynamics model used in the simulation (GOTILWA+ or SORTIE-ND).<b>site</b> – Climatic condition of the virtual stand (xeric, mesic, humid).<b>mng</b> – Management alternative applied (PS00 = no management, PS08 = low-intensity/high-frequency thinning, PS09 = high-intensity/low-frequency thinning).<b>rcp</b> – Climate scenario (current climate = RCP 0.0, RCP 4.5, RCP 8.5).<b>sim_year</b> – Simulation year.<b>BA</b> – Stand basal area (m² ha⁻¹).<b>STV</b> – Standing timber volume (m³ ha⁻¹).<b>DBH</b> – Mean diameter at breast height (cm).<b>YLD</b> – Yield or timber removals derived from management interventions (m³ ha⁻¹).<b>F</b><b>R</b>– Fire risk index.<b>Background on the Simulation Design</b><b>Stand data</b><i>Pinus sylvestris</i> (Scots pine) is a key species in NE Spain due to its growing stock, carbon sequestration capacity, and harvested volume (Gracia et al., 2004). We simulated its response to current management recommendations (adapted from Pique et al., 2017) under three climate scenarios using GOTILWA+ and SORTIE-ND. Parameter estimates were embedded in GOTILWA+, while SORTIE-ND parameters followed Ameztegui et al. (2015).Stand locations from the Spanish National Forest Inventory were selected to represent <b>humid, mesic, and xeric</b> conditions. To ensure consistent inputs across sites and models, we constructed <b>virtual stands</b> with identical structure, populated with <b>1450 </b><b><i>P. sylvestris</i></b><b> trees</b> with DBH between 10–20 cm.<b>Climate change scenarios</b>Baseline (1981–2010) climate data were generated by interpolating daily records from SMC and AEMET using <i>meteoland</i> v1.0.3 R package (De Cáceres et al., 2018). Future conditions were represented by <b>RCP 4.5</b> and <b>RCP 8.5</b> projections. GOTILWA+ used daily climate series, whereas SORTIE-ND used monthly temperature and precipitation inputs.<b>Management alternatives</b>Management objectives followed the ORGEST guidelines (Piqué et al., 2017), prioritizing <b>timber production</b> and <b>fire risk reduction</b>. Strategies combined the shelterwood method with thinning from below (González-Olabarria et al., 2017). The three management alternatives simulated were:<b>PS08 or LIT </b>corresponding to the <i>Low-Intensity / High-Frequency Thinning</i> regime<b>PS09 or HIT</b> corresponding to the <i>High-Intensity / Low-Frequency Thinning</i> regime<b>PS00 or NM</b> corresponding to the <i>No Management</i> alternativeThinning schedules were <b>adaptive</b>, triggered when basal area (BA) and diameter at breast height (DBH) thresholds were reached, while shelterwood stages were implemented at fixed 10-year intervals.