Magnitude and timing of density reduction are key for the resilience to severe drought in conifer-broadleaf mixed forests in Central Europe

We applied a modified forest gap model (ForClim) to depict changes in stand water transpiration via density reduction as a forest adaptation strategy. This approach is the key to analyzing the ecological resilience to drought, stress-induced mortality, and economic efficiency of managed mixed forest stands in Central Europe. The results show that specific geographic conditions and forest composition define the optimal stand density of drought-resilient forests. Reducing stand density has been recognized as a valid strategy to increase forest resilience to drought. Moreover, to develop adaptive management strategies (AMS) under climate change, it is crucial to consider not only drought resilience but also the economic efficiency of alternative AMS proposed to alleviate drought effects. To analyze how decreased inter-tree competition among overstorey trees affects stand vulnerability to drought and its expected yield. We integrated experimental thinning data and historical responses to drought years in a climate-sensitive forest gap model, ForClim. We tested a business as usual (BAU) and three alternative AMS (“do-nothing,” low- and high-intensity overstorey removal) in mixed stands of Norway spruce (Picea abies), silver fir (Abies alba), and European beech (Fagus sylvatica) along an elevational gradient of 520–1020 m a.s.l. in Central Europe. High-intensity overstorey removal in mixed stands of all three species considerably increased forest volume growth resilience to drought and decreased stress-induced mortality by two-thirds vis à vis a “do-nothing” strategy. In sites including only conifer species, forest resilience was equally improved by high- and low-intensity overstorey removal compared to that in the BAU strategy. Regarding the timber economy, high-intensity overstorey removal resulted in a higher economic revenue of mixed stands (~ 22% higher net present value than other strategies) on the high-elevation sites (> 1000 m a.s.l.). Modifying forest density and structure by overstorey removal is principally suitable to increase forest resilience to drought and improve its economic efficiency. The magnitude of the effect however depends on the geographical setting and forest composition.


Introduction
Droughts in combination with higher temperatures are a growing concern in terrestrial ecosystems, particularly for forests where drought-induced declines in tree growth and survival are reported. Climate change scenarios for Central Europe suggest an increase in the occurrence of consecutive dry days, as a function of the increase in temperature and decreases in precipitation and water availability (Greve et al. 2018;Hoegh-Guldberg et al. 2018). The negative impacts of drought on European temperate forests have increased in the past decades; particularly, rising temperatures and high seasonal variability in precipitation increase the likelihood of large-scale forest mortality . Patterns of declining growth rates due to region-specific climate variability have been observed, and the limited species tolerance to hotter droughts impairs the capacity of individual trees and stands to recover from drought, thus increasing their mortality risk .
Drought stress caused by high evaporative demand and low soil water content affects whole-plant hydraulic conductance from the soil to the leaves (Bréda et al. 2006). Tree vulnerability to drought depends on its individual-, population-, and speciesspecific functional properties such as the stomatal regulation strategy or hydraulic architecture (Landsberg and Sands 2011;Martínez-Vilalta and Garcia-Forner 2017). Tree size plays an important role in the drought sensitivity of trees, usually with more detrimental impacts for taller trees (cf. Bennett et al. 2015). Stand characteristics (i.e., density, structure, composition) also interact with species vulnerability to shape forest response to drought (Forrester and Albrecht 2014). Importantly, tree species diversity can enhance forest resistance to drought (e.g., specific mixtures that improve water availability of forest stands; cf. Grossiord et al. 2014). However, factors such as the type of mixture and tree species identity, depending on site aridity, basal area, and age, greatly modulate the beneficial effects of mixtures regarding drought resilience, which thus cannot be generalized (Pardos et al. 2021).
Likewise, stand structure can modulate forest resilience to drought-induced growth declines, on the one hand by lower physiological vulnerability of understorey trees, and especially when small trees are released from competition (Jones et al. 2019). On the other hand, higher stand density often is linked to exacerbated growth reductions during drought due to the decreased availability of soil moisture (Gleason et al. 2017). Novel management strategies to increase the adaptive capacity of forest ecosystems are thus required. This will necessitate an improved understanding of the interactions between drought and management strategies, as the latter are affecting competition among trees for critical limiting resources.
Management can enhance the drought resilience of forests, i.e., their capacity to return to the level of functions observed prior to the event (e.g., growth rate). Through experimental density reduction across a range of forest types and sites, the role of density reduction for resilience in past droughts was investigated (e.g., Elkin et al. 2015;Sohn et al. 2016a), suggesting that in regions with a high probability of drought stress, short and intensive thinning intervals are useful to increase drought resilience. However, forest density reduction may not be appropriate in all conditions because it may increase soil evaporation in residual stands, and thus in turn induce drought stress Gleason et al. 2017). Additionally, species-specific studies including a focus on the effects of tree size and the specific management approach are missing. Moreover, to develop adaptive management strategies in a climate change context, we need to consider not only drought resilience but also the economic efficiency of the management changes that are proposed to alleviate drought effects. For instance, a wide set of alternative management options and measures of forest resilience should be considered, depending on climate change impacts, benefits and costs expected from adaptation, and stakeholders' risk aversion (Yousefpour et al. 2017).
The use of climate-sensitive forest growth simulation models can support the development of new management strategies to increase forest resilience to drought (Fontes et al. 2010;Mäkelä et al. 2000). These models consider the species-and size-specific responses to changes in light and water availability in the stand, as well as local conditions that drive those responses (Mina et al. 2017a;Rasche et al. 2011). In landscapes experiencing drought, simulations of management intensity and harvesting types have been tested to evaluate the resilience of managed forests to climate change (e.g., Niinimäki et al. 2012;Seidl et al. 2017). Forest-enterprise, regional, and global assessments of adaptive management have also tested the effectiveness of forest growth simulation models to explore adaptive management options in decision-making (Thrippleton et al. 2020;Yousefpour et al. 2019;Yousefpour et al. 2018). Despite these recent applications of mechanistic simulation models for promoting forest resilience, only a few model applications consider both ecological (i.e., improving resilience to drought) and economic aspects.
Norway spruce (Picea abies (L.) Karst) and silver fir (Abies alba Mill.) are among the most widespread economically and ecologically important tree species in Europe. Monocultures of Norway spruce have shown less adaptability to changing climate conditions, especially to drought (Lévesque et al. 2013;Sidor et al. 2015), and vulnerability to subsequent disturbances (e.g., bark beetles; Temperli et al. 2013). In contrast, silver fir is characterized by a higher tolerance to drought (Zang et al. 2014). European beech (Fagus sylvatica L.) is the most abundant broadleaf forest tree in Central Europe (Dittmar et al. 2003) and is more drought-resistant than Norway spruce and silver fir (Niinemets and Valladares 2006;Zang et al. 2014). Although the three species can be found together in mixed mountain forests in Central Europe, their specific responses to mean annual temperature and the length of the growing season modify their performance along elevation gradients. At higher elevations, spruce is adapted well to the low prevailing temperatures, while at low elevations with high temperatures, fir and especially beech are adapted better (Kölling 2007). Accordingly, admixing silver fir and European beech is likely a good strategy to improve the resilience of Norway spruce stands in drought-prone environments.
Overstorey removal to improve the drought resilience in mountain mixed forests requires analyses that include species mixture as well as local to regional geographic conditions, e.g., topography, micro-climate, and soil formation (Lindner et al. 2010). Empirical and simulation studies analyzing the effects of contemporary and future climate change report different effects on tree growth and mortality, depending on the elevation and topography of mountain forests in Central Europe (Thrippleton et al. 2020). In the present study, we used historic drought resilience data from six managed sites along an elevational gradient (520-1020 m a.s.l. in southwestern Germany) to assess the effectiveness of overstorey removal for promoting drought resilience. Specifically, we simulated a "business as usual" (BAU) and a "do nothing" management strategy along with two adaptive strategies of high-and low-intensity overstorey removal. We applied the climatesensitive, process-based forest gap model ForClim that was modified to better capture the relationship between stand density (i.e., stand transpiration) and water availability on the one hand and the growth dynamics of these three species on the other hand. Moreover, we used simulated yield to calculate the economic efficiency of these strategies regarding timber production.
As the main objectives of this study, we investigated (1) how the magnitude and frequency of overstorey removals influenced the drought resilience of the stands; (2) how drought and stand composition interact with the resilience of forest growth; (3) how the geographical setting affects the optimal density of drought-resilient forests; (4) and how simulated management strategies affect the economic efficiency (net present value, NPV) of forests. The proposed approach allows for the short-term planning of thinning activities (i.e., over less than a decade) and, as a dynamic ecological-economic modeling system, provides the information necessary to adopt strategic changes in thinning over multiple decades.

Study sites
The data used for the simulation analysis and the design of the strategies came from experimental plots located at six sites in the Black Forest area, Southwest Germany ( Figure 1). The management experiment included an irregular group shelterwood regeneration method ("Femelschlag"), comparing strategies with different regeneration phase lengths, and thus frequency of overstorey tree removal (20, 35, and 50 years) and a control plot (removal of 50% of the 10-year volumetric increment; Weise 1996). Plots of approximately 0.25 ha were established between 1980 and 1981 in 90-to 126-year-old stands dominated by Norway spruce and silver fir, in which no harvesting had taken place for the past 10 years (Puettmann et al. 2009). In every site, we chose an experimental control plot as the basis of our simulations ("business as usual" strategy; BAU). European beech covered 2 to 20% basal area of the sites 220, 222, 223, and 225 ( Figure 1) in the first observation period (1980)(1981). Each of the six sites differs in terms of geographical setting (e.g., parental material, topography, elevation), climate, and stand characteristics, leading to a wide range of light and water availability ( Table 1). The expected mean annual increment (MAI) for spruce and fir monocultures varies between 7 and 6 m 3 ha -1 year -1 for low productivity sites and between 12.5 and 11 m 3 ha -1 year -1 for the medium productivity sites (Forrester et al. 2013). The climate is temperate, with mean annual temperature ranging from 8.3 to 5.5°C and annual precipitation ranging from 940 to 2020 mm along the elevational gradient between 520 and 1020 m a.s.l., respectively. Additional information about the sites and their management is provided in Weise (1996), Puettmann et al. (2009), andForrester et al. (2013).

The climate-sensitive forest gap model ForClim
ForClim (v4.0.1.2) is a climate-sensitive forest gap model developed to simulate the dynamics of temperate forests over a wide range of environmental conditions, especially under varying management regimes and climate change scenarios (Bugmann 1996;Rasche et al. 2012). It has been successfully applied to predict species composition and productivity of temperate and mountain forests in Western Europe (Mina et al. 2017a;Rasche et al. 2012) and North America Gutiérrez et al. 2016). The influence of climate and stand characteristics on demographic processes, i.e., tree growth, regeneration, and mortality, is based on specific ecological assumptions, simulated at the stand level.
The growth dynamics of every tree size cohort on small patches of land (here, 800 m 2 ) are based on species-specific maximum growth rates, which are reduced depending on the extent to which environmental factors are at suboptimal levels (Bugmann 2001;Moore 1989). Soil moisture, soil nitrogen, and degree days (sum of days with average temperature > 5.5°C) are site-specific growth reduction factors, whereas available light and crown length are growth reduction factors related to the development of the cohorts on each patch. Mortality of adult trees is determined based on a combination of a constant "background" mortality, which depends on the species-specific maximum age, and a stress-induced mortality that occurs if the annual stem diameter increment is lower than 0.3 mm or below 10% of the maximum specific diameter increment for three or more consecutive years (Bugmann 1996). Mortality is simulated individually for each tree within Mean annual increment (m 3 ha -1 yr -1 ; A. alba, P. abies) 12.5, 11 7, 6 12, 10 7, 10 11, 9 11, 10 Pre-treatment volume (m 3 ha -1 ) 5 Fig. 1 Location of the six study sites in the Black Forest (sites 221-223, 225), the Swabian-Franconian Forest (site 224), and the South-Western Swabian Alps (site 220). For every site, species composition and elevation class (low, < 700 m a.s.l.; medium, 720-830 m a.s.l.; high, > 1000 m a.s.l.) are indicated a given cohort based on a stochastic approach using Monte Carlo techniques (more details in Bircher et al. 2015;. A flexible management submodel allows for the application of a broad range of common silvicultural treatments such as clear cutting, shelterwood, thinning, or planting (Rasche et al. 2011), and has been proven suitable to represent empirical and analytical harvesting algorithms like single-stem removal or relative diameter classes (Mina et al. 2017b). ForClim is a horizontally non-explicit forest model without interactions between individual simulated patches; hence, tree removals are executed randomly within the stand, but they depend on forest properties. In each patch, the harvesting algorithm selects trees to be removed based on their diameter at breast height (dbh) according to defined classes and rules until a certain volume is reached, using a stochastic Weibull function (Mina et al. 2017a;Rasche et al. 2011).

Improvements in ForClim
Despite several changes to the structure and parameterization of ForClim to better simulate tree growth and light environment (e.g., Mina et al. 2017b), the model tends to underestimate the growth release effect after a disturbance such as a management intervention (Didion et al. 2009). More importantly, the drought intensity experienced by each cohort in ForClim v4.0.1.2 is insensitive to stand density and does not change among stands subject to different management regimes. Elkin et al. (2015) improved the representation of water demand in the stand through an empirical relationship between evapotranspiration and leaf area index in the LandClim model. Yet, competition for water is still not simulated explicitly in ForClim (Huber et al. 2018). This had to be improved to accurately predict the effects of overstorey removal on forest resilience. Additionally, during extreme droughts, the growth reduction function in ForClim can reach zero, which rarely occurs in the field for these species and at these sites, while it strongly affects the calculations of resistance before and recovery after drought. This is due to the lack of a representation of carbon reserves in the model, which are used by trees to grow despite harsh conditions (Guillemot et al. 2017). Thus, beyond the age trend, the tree-level growth time series simulated by ForClim do not feature temporal autocorrelation (Bircher et al. 2015). Additionally, recent studies show an increase in mortality of trees with larger diameters, not only old trees, in European forests (Etzold et al. 2019;Senf et al. 2018). These issues were addressed by modifications to the model as described in Appendix 1.

Identification of severe drought years
We used the mean annual drought index (uDrAn) computed in ForClim to select and characterize drought years (details in  and Elkin et al. (2015)). The calculation of the drought index depends on the monthly ratio between evapotranspiration and the evaporative demand from the soil (sensu Eq. 5 in Appendix). Five drought years were identified, i.e., 1983, 1991, 2003, with high values (uDrAn > 0.2) in the severe drought years 2003 and 2016, while the other drought years were mild (uDrAn < 0.2).

Measurement of resistance, recovery, and resilience
To analyze the effects of overstorey removal strategies on forest growth resilience of the remaining individuals, we used an assessment of engineering resilience that is normalized by system state (Ingrisch and Bahn 2018). We framed the calculation of resilience and its components (Eqs. 1-3) by the indices developed by Lloret et al. (2011). Forest growth responses to drought were measured by mean annual volume increment (MAI; m 3 ha -1 year -1 ) of the trees that were not harvested during the whole experiment, including new and advanced regeneration. MAI was calculated from the simulation results of dbh and height growth. The system state baseline for comparison was the pre-or post-disturbance growth rate (i.e., MAI), averaged over a 4-year period (cf. Anderegg et al. (2015) or DeSoto et al. (2020)): Recovery ¼ PostDr=Dr ð2Þ where PreDr is the average MAI of the 4-year period preceding the drought, Dr is the MAI of the drought year, and PostDr is the average MAI of the 4-year period following the drought. To analyze the temporal resilience of the remaining trees, we calculated the three indices for every year but focused our analysis on the identified drought years. For the drought year 1983, we used a baseline period of 2 years because the simulations start only in 1980-1981. In the case of the drought year 2016, we simulated the stand development until 2020 to compute the recovery and consequent resilience of the baseline of 4 years.
The use of growth based indices to quantify forest resilience has been criticized in the past because of the (a) inconsistency regarding the selection and characterization of drought events and (b) its calculation procedure (Schwarz et al. 2020). We consider that using the annual drought index in ForClim (cf. Subsection 2.4) is a consistent empirical process to select and characterize drought years. As for the second limitation, we agreed with the problems discussed; therefore, to corroborate our analysis of engineering resilience, we also analyzed the stress-induced mortality rate as an indicator of ecological resilience (Nikinmaa et al. 2020;Schwarz et al. 2020). To this end, we evaluated the mortality rate in a moving window of 4 years as the ratio between accumulated mortality (measured as m 3 ) and total stand productivity at the end of the 4-year period (measured as m 3 ).

Adaptive forest management strategies considering an ecological-economic approach
The ecological-economic effects of density reduction were simulated using a variable target diameter harvest of the initial overstorey trees (i.e., trees with dbh ≥ 30 cm). For this, we considered the greater sensitivity of taller trees to drought (Olson et al. 2018;Stovall et al. 2019) as well as the discounted net economic revenue for spruce-fir trees of the diameter classes > 40 cm (based on the regional averaged prices and costs for the period 2000-2016; see Appendix Table 5). We simulated the past BAU strategy for 36 years  along with three alternative management strategies, including two initial overstorey removal intensities and a "do nothing" strategy. Based on the frequency and intensity of the overstorey removal of the original experiment (Puettmann et al. 2009;Weise 1996), we simulated high-and lowintensity overstorey removal strategies. Under high intensity, the initial overstorey was removed over 20 years  in four harvesting intervals of 5 years ( Table 2). The harvesting started with the more drought-vulnerable spruce trees with dbh ≥ 60 cm (Sidor et al. 2015) and followed with silver fir and, when necessary, European beech (both ≥ 60 cm). Once the initial trees in this class had been harvested, we continued in the same species order with trees between 40 and 60 cm, and then between 30 and 40 cm until the initial overstorey was completely harvested. Afterwards, we applied control interventions (i.e., removal of 50% of the 10-year volumetric increment) every 10 years. Under the low-intensity strategy, the entire overstorey would be removed over 50 years . We applied control interventions every 10 years during 1980-2000 to foster the growth of overstorey trees. The first overstorey removal was applied in 2005 with an intensity of 20% of the current overstorey volume. We continued the lowintensity removal with 10% increases in the intensity rate every 5 years based on the overstorey volume measured in 2005. We used the same dbh class and species order as in the highintensity strategy. Finally, under the "do nothing" strategy, we only applied a control intervention in the initial year (1980)(1981) and continued the simulation without any management intervention. We used identical harvesting strategies for all sites.

Model evaluation and application
We evaluated the latest version of ForClim (v.4.0.1.2, as described in Appendix 1) against forest inventory data from the original experimental plots. For each site, forest dynamics were simulated under BAU, "do nothing," and the two adaptive management strategies. Climatic conditions were based on regionalized station-based data of precipitation and air temperature with a monthly resolution (German Weather Service; Dietrich et al. 2019). The stands were initialized using dbh data from the forest inventory conducted between 1980 and 1981. Each tree was randomly allocated to the ForClim patches, whose number was obtained by dividing the original site area by the patch size of 800 m 2 . To reduce stochastic noise in the simulations (Didion et al. 2009), individual patches were replicated 200 times. Tree regeneration was allowed to occur naturally in every simulation constrained by site conditions and browsing pressure (Appendix Table 4). From the outcomes of the model at the patch level, we upscaled the values to 1 hectare (i.e., basal area, standing volume, and stress-induced mortality). To evaluate the goodness-of-fit between BAU inventory observations and Table 2 Frequency and intensity of the harvest interventions (BAU in reality and three alternative strategies) in the simulations Year Adaptive management strategy BAU (business as usual) "Do nothing" Low-intensity initial overstorey removal (Low) High-intensity initial overstorey removal (High) *Removal of 50% of the 10-year volume increment for trees with a dbh ≥ 10 cm simulations, we used the relative root mean square error (RMSE) and the percent bias (pbias). We simulated the four forest management strategies in every site and analyzed the results of the resistance, recovery, resilience, and stress-induced mortality rates for each strategy. We looked for their differences during the observation period and identified drought years, considering all six sites individually and their forest composition (fir-spruce and fir-sprucebeech) and elevation (low, medium, and high). We evaluated pairwise differences among management strategies (t-test for paired samples). To analyze the economic implications of the four management strategies, we calculated the expected economic benefits in terms of net present value (NPV), estimating the cash flow of harvest interventions and the corresponding changes in stumpage value. Based on the aims of our study and expected low aversion to risk of the forest managers of our research area, we used a discount rate of 1% (Yousefpour et al. 2017). Averaged prices and costs by diameter classes for Baden-Württemberg for the period 2000-2016 (see Appendix Table 5) are used in the calculation of NPV: where R t is the net cash inflow-outflow during a single period t, i is the discount rate, and n is the number of periods. All calculations and analyses were performed in the R environment (R 3.5.3; Development Core Team R 2019).

Model evaluation
ForClim simulations of basal area (BA) and standing volume matched the observed data over 36 years (eight inventories from 1980 to 2016) quite well, with a percentage bias < 2% and < 5%, respectively. As Figure 2 shows, in fir-spruce mixed sites, simulations of BA marginally underestimated the observed values, ranging from −1.9 to −0.8%. There were negligible underestimation and overestimation of BA in sites where European beech is present, ranging from −0.6 to +1.5%. Standing volume showed a slightly higher bias (−4.7 to +0.6%) compared to BA (Appendix Fig. 10). Results of the relative RMSE of BA and standing volume showed similar patterns of differences among sites. Overall, ForClim simulated stand development under current management and climate conditions very well (additional information is provided in Appendix 2).

Model application
The intra-site comparison among overstorey removal strategies indicated that volume MAI (i.e., forest growth) diverged between the strategies along the simulation period ( Figure 3). During the first 20 years, the high-intensity overstorey removal showed the lowest volume MAI at all sites, while the rest of the strategies (i.e., low-intensity, BAU, and "do nothing") maintained forests with nearly complete stock and similar high MAI values. However, in 2009, 1 year before a control intervention and almost a decade after the complete initial overstorey removal under the high-intensity treatment, MAI increased in all sites by~0.7 m 3 ha -1 year -1 . During the first two decades of the experiment, the low-intensity removal strategy showed a positive but decreasing trend in MAI, bỹ 0.68 m 3 ha -1 year -1 and~0.42 m 3 ha -1 year -1 , respectively. After 2005, MAI decreased by~-4.7 m 3 ha -1 year -1 over 10 years, as more overstorey trees were removed. At the highelevation sites, the MAI of the high-intensity strategy exceeded that of the low-intensity strategy by the end of the simulation, with increases of +1.4 m 3 ha -1 year -1 at site 222 and + 2.4 m 3 ha -1 year -1 at site 223. Overall, MAI was positively affected by standing volume and tree size but negatively affected by the intensity and timing of overstorey removal.
Forest growth responses to overstorey removal during drought years differed by species composition and the geographical setting. High-elevation stands (sites 222 and 223) with moderate slope (> 12°) and higher soil water holding capacity ("bucket size," BS = 20 cm) showed a higher volume MAI under the high-intensity overstorey removal. Hence, no drought-induced negative effect on growth rates was detectable at these sites under the high-intensity strategy (~+ 12% higher MAI during the severe drought of 2003; Figure 3). Sites below 1000 m a.s.l., however, showed an increased sensitivity to drought events under all strategies (e.g., decrease of MAI in 2003). The volume MAI showed an entirely different magnitude of effect in fir-spruce-beech sites under the highintensity strategy, depending on its elevation, slope, and BS. In site 220, located at 830 m a.s.l. with a gentle slope (3°) and BS of 20 cm, MAI increased by + 4.4% during the drought of 2003. In site 225, located at 700 m a.s.l. with a moderate slope (14°) and BS of 14 cm, volume MAI decreased by −7.7%. For the same overstorey removal strategy in fir-spruce sites, MAI decreased by~−38%. Site 224, located at 520 m a.s.l. with a moderate slope and BS of 14 cm, showed the greatest decrease in volume MAI by −66.5%. Site 221, at 720 m a.s.l. with a steep-south slope of 26°and BS of 12 cm, showed a decrease of −9.3% during the drought year 2003.
Under the low-intensity strategy, across the entire elevational gradient, MAI decreased by~−14% in sites with three species and by~−70% in sites with a fir-spruce composition. Lower elevation (< 1000 m a.s.l.) and smaller bucket size (< 14 cm) affected negatively the volume MAI. For the BAU and "do nothing" strategies, the 2003 severe drought reduced volume MAI by~−15% and~−70% in sites with three and two species, respectively. Overall, across the complete elevation gradient, high-intensity overstorey removal decreased the drought effect on volume MAI, with a higher magnitude in stands including European beech and higher water availability.

Forest resilience, resistance, and recovery
3.3.1 Resilience, resistance, recovery, and management strategies We analyzed forest resilience, resistance, and recovery from drought of the trees remaining after management interventions. The stands under the "do nothing" and BAU strategies behaved similarly, but we found differences between these and stands under both the low-and highintensity strategies (Figure 4). Analyzing the five drought years, high-intensity overstorey removal had the highest average growth resilience (~1.05) and resistance (~0.93) among all strategies. It also was the lowest-performing strategy in terms of growth recovery (~1.29), but it was still larger than 1.0, indicating the capacity of stand growth to recover following drought events. The magnitude of the effects in terms of resilience and resistance to drought among the other management strategies decreased with increasing presence of overstorey trees and total stand volume. Consequently, stand performance in terms of resilience and resistance declined progressively under the low-intensity overstorey  removal, BAU, and "do nothing" strategies, respectively. The inverse order of strategies described the progressive decline in stand performance regarding recovery. Overall, during drought years, forest resilience and resistance were low under strategies maintaining a higher standing volume and more overstorey trees, whereas these strategies showed higher recovery. When overstorey removal occurred just before or after a drought year (± 4 years), forest resilience increased compared to the control interventions of BAU. Specifically, during the droughts of 1983, 1991, and 2003, high-intensity overstorey removal induced the highest resistance to drought (~+1%,+ 7%, and~+25%, respectively). Due to these increases of resistance, forest resilience was also highest during the droughts of 1991 (~+9.5%) and 2003 (~+7%). Meanwhile, during the 2006 drought, low-intensity overstorey removal showed the highest resilience (~+6%), resistance (~+5%), and recovery (~+1%) values compared to the other strategies. Finally, in the drought of 2016, only the growth rate under the low-intensity strategy was resilient to drought (~1.01) due to increased resistance (~+8%) and recovery (~+1%) relative to the other strategies. The remaining strategies with control interventions showed resilience lower than 1, indicating that the previous growth rate was not recovered in the 4 years following the drought. Thus, compared to control interventions, the application of overstorey removal interventions positively affects the resilience (mostly by increasing forest resistance to droughts; see details in Appendix Figs. 11 and 12).  Table 2). Gray vertical lines indicate drought years (1983, 1991, 2003, 2006, and 2016) 3.3.2 Resilience, resistance, recovery, and forest composition In the analysis of forest resilience under management strategies by drought year, we observed a different pattern for the two site groups with different forest compositions (i.e., firspruce vs. fir-spruce-beech). During the first decade of the experiment, both types of stand composition showed an increasing trend in forest resilience under all strategies (~+17% in fir-spruce-beech and~+10% in fir-spruce), while in both cases, the magnitude of the effect was higher under the high-intensity overstorey removal ( Figure 5). After the drought of 1991, the decreasing trend in resilience under all strategies (− 9%) was reversed only slightly in fir-spruce-beech stands, while the stands with fir-spruce composition showed a continuously decreasing trend (~−36%,~−42%,~−32%, −22% under the "do nothing," BAU, low-intensity, and highintensity overstorey removal strategies, respectively). nothing," BAU, low-intensity, and high-intensity overstorey removal strategies, respectively). During the same period, the magnitude of the effect was smaller in fir-spruce-beech stands under all strategies (~+6%). Finally, between the last two droughts (2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016), forest resilience decreased in both stand composition types, but to a lower extent in fir-sprucebeech stands (~−13%), while in fir-spruce stands, the effect was higher, especially under the "do nothing" and BAU strategies (~−33%). On average across the whole period, the resilience of fir-spruce-beech stands was higher than that of firspruce stands.
The effect of management strategies on forest resistance and recovery of the two site groups also depended on drought severity. In particular, in fir-spruce stands under all strategies, the effects of the severe droughts in 2003 and 2016 decreased on average its resistance by~−67.5% (Appendix Fig. 13). This decrease was compensated by an increase in recovery in the next 4 years (~+364%; Appendix Fig. 14). During the mild droughts of 1983, 1991, and 2006, its resistance and recovery followed the same pattern but to a lower extent (~−9.6% and 14.1%, respectively). For fir-spruce-beech stands, during severe droughts, the decrease in resistance was smaller than in the other site group (~−14.1%), while its recovery showed an inverse response, being higher than the decrease in resistance (~+65.1%). During the mild droughts, the increase in forest resistance was counteracted by a decrease in recovery (~+3.1% and~−3%, respectively). Overall, severe droughts had a larger effect on resistance and recovery than mild droughts. In addition, the presence of European beech in the mixture favored a high and more stable resistance and recovery.

Mortality
The stress-induced mortality rate decreased depending on the intensity of forest management, with reduced mortality under higher overstorey removal strategies, especially compared to the "do nothing" strategy. The average reduction of the mortality rate during the experiment  was highest under high-and low-intensity overstorey removal (~−57.3% and~−61.7%). Under BAU, the reduction was considerably lower (~−21.1%), while under the "do nothing" strategy, the mortality rate increased slightly (~+5.8%; Figure 6). Under the overstorey removal strategies, high-intensity overstorey removal every 5 years resulted in a decrease of the mortality rate of~−66.7% between 1980 and 2000. For the same stands, the mortality rate increased (~+28%) when a control intervention was executed only every 10 years in the period 2001-2016. The stands with low-intensity overstorey removal showed a decreased mortality rate (~−19%) for the period 1980-2000 when control interventions were applied. In 2001-2016, when overstorey removal started under the lowintensity management, the mortality rate decreased by− 52.7%. Overall, the magnitude of the effect in the reduction of the mortality rate depended on the type and timing of the management intervention. The mean dbh of dead trees decreased over the simulation period under the "do nothing" strategy, varying by elevation. Higher stress-induced mortality in trees with a dbh < 25 cm was found at low and medium elevations, with a mean dbh of living trees of 39 cm and 37 cm, respectively. At high elevations, the mean dbh of living trees was 35 cm, and trees with a dbh < 20 cm were most affected by the increased stressinduced mortality rate (Appendix Fig. 15). Altogether, these responses indicate a high sensitivity of the stress-induced mortality rate to the timing and magnitude of overstorey removal. Thus, following a "do nothing" strategy in the face of increasing drought impacts on forests is likely to result in higher mortality, especially among understorey trees.

Economic assessment
Net present value (NPV) was affected by the overstorey removal strategies. Compared to BAU, NPV increased with increasing overstorey removal intensity and frequency. In high-elevation fir-spruce-beech stands (sites 222 and 223), the high-intensity removal performed best, with an averagẽ +22% increase in NPV from 1980 to 2016 (Figure 7). Medium-and low-elevation stands composed of three species showed similar results. At the medium-elevation stand (site 220), the high-intensity strategy resulted in a +49% higher NPV than the BAU strategy. At the low-elevation stand (site 225), the high-intensity strategy showed a +10% higher NPV. Fir-spruce stands at medium and low elevation resulted in a higher NPV under the high-intensity strategy than BAU, but the difference was smaller than that of fir-spruce-beech sites. At sites 221 and 224, the average expected value of the stand was~+3.7% higher under the high-intensity strategy. Overall, in stands with three species, NPV under the BAU scenario was lower than NPV under either overstorey removal strategy but still higher than that under the "do nothing" management. In fir-spruce sites, the NPV of BAU was slightly lower than that of the overstorey removal strategies but higher than under the "do nothing" management.

Discussion
Our results show a strong positive effect of overstorey removal strategies on forest resilience to drought. In contrast, a "do nothing" strategy is likely to increase tree mortality, particularly for understorey trees. Our simulations also indicate differences related to the geographical setting and stand composition, with high resilience in forest stands at low-and medium-elevation sites. However, the economic output of high-intensity strategies on high-elevation sites was higher than that under other strategies.

Forest structure and forest resilience
We found that a high frequency and high intensity in overstorey removals promoted forest resilience following drought events, confirming the results of Puettmann et al. (2009) who evaluated empirical indices in the same study sites and explained this increase in resilience as a consequence of release effects following density reduction. In contrast to the simulations under the "do nothing" and BAU management strategies, the gaps produced by the overstorey removal strategies positively influenced the growth of new regeneration or otherwise suppressed young trees ("advanced regeneration"; Jones et al. 2009).
The magnitude of the growth response to thinning strategies is likely to vary with tree size and light-use efficiency. For example, light absorption and light-use efficiency are strongly influenced by the structure and composition of individual trees' neighborhoods along this elevation gradient (Forrester and Albrecht 2014). Thus, by starting the harvesting with taller trees, adaptive management strategies allow for the growth of shade-tolerant silver fir trees that tolerate dense neighborhoods dominated by Norway spruce (Forrester et al. 2013). Furthermore, the harvest of taller trees increased the light share for young trees, which reduced its stress-induced mortality in combination with greater soil resource availability (Berntson and Wayne 2000). Moreover, the growth and survival of previously suppressed trees benefited from overstorey removal in short intervals (i.e., 5-10 years; cf. Elkin et al. 2015;Sohn et al. 2016b). Year Stress−induced mortality rate Do nothing BAU Low High Fig. 6 Stress-induced mortality rate among management strategies. Gray vertical lines indicate drought years (1983, 1991, 2003, 2006, and 2016) Forest resilience to drought was also enhanced by improved water availability, driven by the decreased inter-tree competition among overstorey removal strategies in coniferbroadleaf mixed forests. Greater decreases in resistance during severe drought events were associated with greater reductions in resilience. Although growth recovered in the next 4 years, only in stands with the three species were both resistance and recovery maintained during the next drought events. The stability in resistance and recovery (i.e., values close to 1.0) of stands under high-intensity overstorey removal conferred a higher resilience to both site groups. This is consistent with recent reviews linking mixing effects with growth in highly thinned stands. For example, Forrester and Bauhus (2016) suggested that mixing effects due to species-specific radial development and shade tolerance drive nutrient and water availability and determine the growth response in fir-sprucebeech stands. Additionally, other large-scale studies reported that drought resistance increases with the ratio of soil moisture to potential evapotranspiration (Gazol et al. 2017;Stuart-Haëntjens et al. 2018). In that sense, we observed an increase in drought resistance in stands with low total standing volume and fewer overstorey trees across the entire elevation gradient, indicating that higher water availability resulting from intensive overstorey removal may increase drought resistance. Drought recovery was higher in the short term (e.g., over 4 years) in stands maintaining high total standing volume and more overstorey trees. However, this enhanced growth capacity vanished fast. Year NPV (thousand € ha −1 ) Do nothing BAU Low High Fig. 7 Net present value (NPV) of the stands using 1% as the discount rate. Gray vertical lines indicate drought years (1983, 1991, 2003, 2006, and 2016)

Modelling adaptive forest management
Process-based models are crucial to analyze the effects of forest structural changes on forest resilience to drought. Simulated growth dynamics should be sensitive to stand structural changes to test forest resilience. In ForClim, the structure representation satisfactorily captures the growth sensitivity of young trees to increases in light availability resulting from overstorey removal. Additionally, the simulations consider a vegetation-dependent water demand factor that relates the leaf area index to stand evapotranspiration (see Appendix 1;Elkin et al. 2015, Thrippleton et al. 2018. Thereby, overstorey removal reduced stand transpiration by diminishing foliage area. This, in combination with increased light availability for understorey trees, was particularly beneficial for shadetolerant species like silver fir (Grassi and Umberto Bagnaressi 2001). Furthermore, the integration of adaptive management strategies in a climate-sensitive forest gap model allowed us to test the effects of decreased resource competition and larger growth rate increment, as demonstrated by the differences among overstorey removal strategies. Although we have improved the representation of short-term drought effects on forest resilience and mortality (in ForClim), for longer drought periods, other issues (e.g., representation of carbon reserve pools; Guillemot et al. 2017) ought to be considered explicitly and by consideration of mechanisms, rather than heuristically as in our approach (cf. Appendix 1).

The effect of forest composition and the geographical setting
As mentioned above, mixtures of conifer and broadleaf tree species can improve forest resilience to drought. Reviewing thinning effects on temperate forests, Sohn et al. (2016b) found that individual-tree radial growth was on average 125% higher in heavily thinned compared to that in unthinned stands. They also found that soil water availability was 54% higher in heavily thinned compared to that in unthinned stands during a drought year. The positive effect of increased drought resilience in mixtures compared to monocultures has also been measured recently in closed even-aged, well to fully stock stands (Pardos et al. 2021). They found a strong positive effect of conifer-broadleaf mixed stands on drought resilience and resistance, and suggested that the considerably different root morphology and architecture among species leads to trait complementarity, thereby improving the exploitation of belowground resources. Conifer-conifer mixtures did not differ from monospecific conifer stands in terms of resistance and resilience, which could be linked to their greater similarities in traits and strategies to cope with drought (i.e., isohydric behavior). Similarly, in our simulations, the fir-spruce-beech combination was found to be more resistant and resilient to severe droughts than the fir-spruce stands.
The long-term resilience of mixed stands through overstorey removal was greatly affected by the elevational gradient and soil water holding capacity. Negative impacts from climate change are likely to occur in forests at low elevation (≤ 1000 m a.s.l). In contrast, high-elevation forests are more likely to benefit from warmer conditions in the absence of disturbances (Thrippleton et al. 2020). We observed that drought strongly affects volume growth across management strategies on lowand medium-elevation sites, especially those with low soil water holding capacity. Meanwhile, at high-elevation sites with high soil water holding capacity, only a slight effect was measured. These findings are consistent with the understanding that growth responses vary greatly according to topography, which plays a key role in altering stand water availability and determining local climate conditions, also according to aspect and slope (Zou et al. 2007). Individual trees feature specific mechanisms to maintain the integrity of their hydraulic system and avoid cavitation. In the short term, trees utilize dynamic and reversible regulation processes via these mechanisms, such as stomatal closure to reduce transpiration (Bréda et al. 2006). Our study demonstrates the short-term positive effects of adaptive forest management on the growth resilience of mixed forests, but the magnitude of the effect depends greatly on forest composition and the geographical setting. Reduced stand transpiration by high-intensity overstorey removal may help to regulate the negative impacts of climate change on low-elevation forests (≤ 1000 m a.s.l), increasing resistance, recovery, and resilience to drought conditions over 4-year periods (according to our state system baseline; Figure 3). Overall, an analysis of more experimental sites is necessary to support this conclusion across Central European temperate forests.

Forest resilience economy
The net present value was found to increase concurrently with the frequency and intensity of overstorey removal, offering a useful short-term management strategy to increase the economic efficiency of forest management. However, over longer periods and under changing growing conditions and different site qualities, density reduction responses may be overestimated (Pretzsch 2020). Elkin et al. (2015) showed that the effects of management are quickly disappearing in the growth response and would be rather costly to maintain via frequent interventions. The range of growth responses under different overstorey removal strategies in conifer-broadleaf stands ( Figure 5) allows for the future exploration of different levels of frequency and intensity of overstorey reductions. This can provide insights to optimize drought-resilient forest economic efficiency. In comparison, overstorey removals in fir-spruce stands did not considerably affect their economic efficiency compared to BAU management, as lower MAI counterbalanced the effects of the early intensive interventions on NPV. The results of our economic efficiency analysis are in line with those of other authors such as Knoke et al. (2021). They used a framework of economic resilience and found that a forest management system with higher structural mixture (e.g., continuous cover forestry) performs economically and ecologically better than a monospecific system (e.g., clear fell). The economic resilience focuses on the time needed for a forest to gain the previous value after a disturbance. Another concept analyzed by Pukkala (2021) is using resilience and economic benefit in multi-criteria evaluation of the social benefit of forest management to rank alternative silvicultural systems. Notwithstanding, analyzing ecological and economic resilience is a useful framework for decisionmakers to evaluate the robustness of management regimes to deal with droughts.
It should be noted that we assumed fixed timber prices and a constant discount rate for the entire period to calculate the NPV, ignoring that volatile prices and changes of the discount rate are an important source of risk for management investments (Knoke et al. 2001;Yousefpour et al. 2012). However, recent studies like Messerer et al. (2020) that used a simple modeling approach found that including empirical growth responses after harvest interventions greatly reduces the economic risk of optimized, uneven-aged forest management. Therefore, we recommend that future studies consider price and discount rate changes due to drought events or other disturbances in addition to process-based models.
This study contributes to the clarification that overstorey removal as an adaptation strategy may improve the ecological and economic efficiency in Central European conifer-broadleaf forests. Brèteau-Amores et al. (2019) illustrated the importance of considering ecological and economic benefits to evaluate drought-adaptation strategies. Moreover, Radke et al. (2020) recently demonstrated how combining economic metrics with ecological criteria (i.e., volume growth) can be used to delineate robust adaptive management practices. However, such joint analyses are largely lacking in the literature. Notwithstanding, our research further highlights the importance of increasing forest resilience for long-term adaptive management planning. Additional studies that simultaneously analyze resilience indices, mortality rates, and economic efficiency are crucial for the multi-dimensional evaluation of drought impacts and to identify the best drought management strategies.

Conclusion
In conclusion, the tested adaptive forest management strategies to remove overstorey trees are a promising strategy to increase ecological and economic resilience to severe drought events. Our study demonstrates the short-term positive effects of adaptive forest management on the growth resilience of mixed forests, but the magnitude of the effect depends greatly on forest composition and the geographical setting. However, the enhanced growth capacity vanished fast; therefore, the evaluation of adaptive forest management in the long term to compare the benefits from different strategies is required. Nevertheless, our study provides the highly relevant suggestion that under exacerbated drought conditions, maintaining unmanaged forests (i.e., the "do nothing" strategy) may not be an efficient strategy, especially because new and "advanced" regeneration will be negatively affected by drought conditions. This study shows an innovative combination of diverse historical data and forest models to address the potential consequences of drought in forest management under a framework of engineering, ecological, and economic resilience.

Appendix 1
Summary for Appendix 1

Section in manuscript Subsection
Item Description 2. Methods 2.3. Improvements in ForClim Figure 8 Effect of stand density on the evapotranspirative demand from vegetation *** Growth function Table 3 Minimum and maximum annual tree-ring width (TRW) ratios *** Alternative background mortality 2.4. Identification of severe drought years *** Description of the process Figure 9 Mean annual drought index Additional information on model parameters Table 4 Site-specific parameters used in ForClim Additional data for the calculation of economic indicators Effect of stand density on the evapotranspirative demand from vegetation The evapotranspirative demand (D) and water supply from the soil (S) are the base of the drought index calculation in ForClim . Elkin et al. (2015) functionally coupled the vegetation leaf area index (LAI) to the water supply side to expand the framework and explicitly consider the effect of vegetation in calculating the drought index. That approach considered empirical data from Moreno and Cubera (2008), which estimated the lowest water use at an LAI = 0 (unvegetated condition) and increased by 25% at an LAI = 8 (dense stand). However, the new approach did not consider the evapotranspirative demand from different vegetation density levels like our research goals needed. To solve this issue, we use a function linking stand demand (D stand ) to vegetation LAI via a vegetation-dependent demand factor (F vegetation ) (see Eq. 5). The new equation was initially developed by Thrippleton et al. (2018) to account for the contribution of different stand vertical layers to the stand water demand in the spatially explicit model LandClim. The new demand factor is based on the original formulation in  and in the empirical linear relationship between stand water use and LAI used in Elkin et al. (2015) (see Figure 8).
where PET is potential evapotranspiration (based on the approach of Thornthwaite and Mather, 1957) and P i is the amount of water intercepted by the canopy (see  and Thrippleton et al. (2018) for further details). Here, we present the linear equation used to calculate F vegetation , dependent on LAI (Eq. 6). In ForClim, LAI is sensitive to management interventions. Therefore, the vegetation-dependent demand factor for potential evapotranspiration in the stand is sensitive to stand density variations as needed for our research goals.
Growth function To avoid zero diameter increments in the model and constrain the estimations of resistance and recovery, we derived a minimum (and maximum) change of the diameter increment from one year to the next for each species using field data from research plots and implemented those thresholds in ForClim. Our calculation was based on the pancontinental tree-ring width (TRW) database by Cailleret et al. (2017) and on TRW measurements carried out at the study sites and analyzed following standard dendrochronological procedures (cf. Speer 2010). We then calculated the speciesspecific annual ratios of any two subsequent diameter increments, considering a minimum and maximum envelope of 99.9% of the resulting data cloud to avoid extreme values ( Table 3). The modified version of ForClim (v4.0.1.2) uses these measured TRW ratios to constrain the variability of simulated diameter increments from 1 year to the next when the simulated increment ratio is smaller or higher than the measured TRW ratio.
Alternative background mortality We used a size-dependent background mortality function rather than an age-dependent function to capture the vulnerability of large trees to climate warming (Stegen et al. 2011;Xu et al. 2017). The alternative theoretical function had already been tested in ForClim by , with beneficial results.
Identification of drought years using ForClim v4.0.1.2 We used the mean annual drought index (uDrAn) in ForClim to identify drought years. To calculate the uDrAn, we used the  Elkin et al. 2015) monthly climate data from Dietrich et al. (2019). The index's positive values indicate an increase in yearly potential evapotranspiration and/or decreased precipitation (Figure 9). Strong drought effects were simulated in sites with lower soil water potential (i.e., BS < 14 cm, steep to moderate slope; sites 221, 224, 225). Sites with higher soil water potential (i.e., BS = 20 cm, moderate to flat slope; sites 220, 222, 223) showed smaller drought effects. For the identification of drought years, we identified years with uDrAn > 0.2 as severe drought events (i.e., years 2003 and 2016); then, we selected years with values   , 1983, 1991, 2003, 2006, and 2016) 0.1 < uDrAn < 0.2 as mild droughts (i.e., years 1983 and 1991) and finally for values 0 < uDrAn < 0.1 consistent at all sites as mild droughts (i.e., year 2006).

Additional information on model parameters in
The values in Table 4 indicate site-specific environmental and biotic conditions that act as species-specific limitations for sapling (dbh = 1.27 cm) establishment (e.g., browsing pressure, soil moisture derived from the bucket size), and tree growth reduction factors (e.g., nitrogen availability, soil moisture derived from the bucket size). The values are based on data provided by Forstliche Versuchs-und Forschungsanstalt Baden-Württemberg (FVA).   Table 6 Percentage of bias Table 7 Root mean squared error 3.2. Model application Figure 10 Standing volume under BAU 3.3.1. Resilience, resistance, recovery, and management strategies Figure 11 Differences in forest resistance  Table 8 Averaged NPV, discount rate 1% Table 9 Averaged LEV, discount rate 1% Table 10 Averaged NPV, discount rate 2% Table 11 Averaged LEV, discount rate 2%

Additional data for the calculation of economic indicators
PBIAS ¼ ∑ sim−obs ð Þ =n ∑obs=n * 100 (8) where sim is the result of the model for the inventory year and obs is the measurement in the inventory. We considered the total number of measurements in the inventory (the inventory includes measurements every 5 years approximately).  Fig. 14 Forest recovery of remaining trees after management strategies separated by species composition. Recovery = 1: the average growth rate before the drought was maintained after the drought event. Recovery > 1: the growth rate was increased after the drought event. Recovery < 1: the growth rate was decreased after the drought event. Gray vertical lines indicate drought years (1983, 1991, 2003, 2006, and 2016). The values were averaged over stands composed of fir-spruce-beech and fir-spruce trees  Fig. 13 Forest resistance to drought of the remaining trees after management strategies separated by species composition. Resistance = 1: the average growth rate before the drought was maintained after the drought event. Resistance > 1: the growth rate was increased after the drought event. Resistance < 1: the growth rate was decreased after the drought event. Gray vertical lines indicate drought years (1983, 1991, 2003, 2006, and 2016). The values were averaged over stands composed of firspruce-beech and fir-spruce trees   (1983, 1991, 2003, 2006, and 2016) Declarations Consent for publication All the authors gave their informed consent to this publication and its content.

Conflict of interest
The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's