Quantifying the effect of persistent dryer climates on forest productivity and implications for forest planning: a case study in northern Germany

Forest management decisions are based on expectations of future developments. For sound decisions it is essential to accurately predict the expected values in future developments and to account for their inherent uncertainty, for example the impact of climate change on forests. Changing climatic conditions affect forest productivity and alter the risk profile of forests and forest enterprises. Intensifying drought stress is seen as one major risk factor threatening forest management in the north German lowlands. Drought stress reduces tree growth and vitality and might even trigger mortality. But so far, it is not possible to quantify effects of a persistent dryer climate on forest productivity at a level suitable for forest management. We apply a well-established single-tree forest growth simulator to quantify the effect of persistent dryer climates on future forest productivity. We analyse the growth of Scots pine (Pinus sylvestris L.), European beech (Fagus sylvatica L.) and oak (Quercus robur L. and Quercus petraea (Matt.) Liebl.) in two forest regions in the north German lowlands for a time interval of 60 years until 2070. The growth response under three different climate projections is compared to a baseline scenario. The results show clear differences in volume increment to persistent dryer climates between tree species. The findings exhibit regional differences and temporal trends. While mean annual increment at biological rotation age of Scots pine and oak predominantly benefits from the projected climate conditions until 2070, beech might suffer losses of up to 3 m3·ha–1·yr–1 depending on climate scenario and region. However, in the projection period 2051 to 2070 the uncertainty ranges comprise positive as well as negative climatic effects for all species. The projected changes in forest growth serve as quantitative contributions to provide decision support in the evaluation of, for example, species future site suitability and timber supply assessments. The analysis of productivity changes under persistent dryer climate complements the drought vulnerability assessment which is applied in practical forestry in northwestern Germany today. The projected species’ productivity has strong implications for forest management and the inherent uncertainty needs to be accounted for.

Background A more and more aggravating climate change shakes one of the key features of forest planning, i.e. to assign species for cultivation according to site suitability (Brang et al. 2008). Formerly, species selection depended on static parameters, e.g. soil nutrients, site water budget and silviculture region as a substitute for regional climatic conditions. In the meantime, the rapidly changing climate affects the primary growing conditions well below a classical rotation period (e.g. Boisvenue and Running 2006;Pretzsch et al. 2014;Kohnle et al. 2014;Aertsen et al. 2014). Thus, forest planning has to account for dynamic site conditions and management decisions are based on expectations of future site developments. Besides projections of site conditions species selection for future forests depends on species resilience, its productivity and ecological conduciveness to name but a few. For sound decisions it is essential to accurately predict the expected values in future developments and to account for their inherent uncertainty, for example the impact of different climate change scenarios on forests. Spittlehouse (2005) categorizes forest management actions to encounter climate change into societal adaptation, adaptation of the forest and adaptation to the forest. Societal adaptation covers adjusting people's expectations on forest services and changes in forest policy on management and conservation. Adaptation of the forest includes, for example, species selection. And adaptation to the forest comprises amongst others applying the adequate silvicultural regime, harvesting and utilization of more salvage wood and modifications in wood processing. All these measures can be proactive approaches to reduce vulnerability and enhance recovery in advance or reactive responses to diminish the impact afterwards (Ogden and Innes 2007).
According to Spittlehouse and Stewart (2003) determination of vulnerabilities and critical thresholds is an integral part of forest management adaptation to climate change. One of the core issues for forest management in the north German lowlands is the vulnerability of forests to projected changes in precipitation pattern. The most likely reduced rainfall in the growing season and increasing evapotranspiration due to rising temperatures will result in persistent dryer conditions or even drought stress on many forest sites. This will most likely affect forest productivity and alter the risk profile of forests and forest enterprises (Lindner et al. 2010(Lindner et al. , 2014. As drought stress reduces tree growth and vitality and might even trigger tree mortality (Hanson and Weltzin 2000;Allen et al. 2010), adaptation measures should account for the effects of persistent dryer climate conditions on forests (Hogg and Bernier 2005). Two main future drought types can be distinguished. First, intense and more frequent dry spells in the growing season will occur in central Europe (IPCC 2013;Lindner et al. 2014). Second, the decrease of mean precipitation in the growing season along with increasing mean temperatures will result in long-term persistent dryer conditions. Many studies quantify growth reactions after short-term extreme drought events (e.g. Orwig and Abrams 1997;Klos et al. 2009;Michelot et al. 2012;Weber et al. 2013). These exemplary studies confirm the obviously negative effect of extreme drought events on tree growth while evincing differences according to for example species, sites, forest layers or distribution range. While these growth declines after singular short-term drought events are mostly reversible on many sites (Bigler et al. 2006;Zingg and Bürgi 2008;Sutmöller et al. 2017), still with increasing frequency of drought events the short-term growth reductions might aggregate to considerable losses in the long run.
Other case studies investigate the long-term effect of persistent dryer conditions on forest productivity. Kint et al. (2012) observed a drought-induced growth decline for beech in Belgium since the 1960s. Piovesan et al. (2008) related the reduction in stand basal area increment observed for beech forests in the central Appennines, Italy, since the last three to four decades to dryer conditions. Furthermore, they found evidence for a stronger effect of temperature increase to drought than reduction of precipitation. They conclude that for beech persistent drought stress will overcompensate the positive influence of an earlier beginning of the growing season. Barber et al. (2000) assume drought stress as one possible explanation for the lessening of the positive effect of temperature increase on radial growth of white spruce in interior Alaska. Based on observed radial growth data from Scots pine, European beech and Pedunculate oak in northern Germany Bauwe et al. (2015) developed regression models and projected growth until 2100 using the emission scenario A1B for future climate conditions. Their simulation results show growth reductions for beech and oak with increasingly dryer conditions while Scots pine did not display a clear trend.
These case studies give evidence for long-term growth trends which correspond to persistent drought stress. But so far, it is not possible to satisfactorily quantify effects of persistent dryer climates on a level suitable for forest management (e.g. Adams et al. 2009;Choat et al. 2012;Anderegg et al. 2013). Until now, in Germany forest planning relies on drought vulnerability thresholds which are based on expert knowledge and/or singular observations on monitoring plots or experimental sites or current presence/absence observations in relation to various environmental factors (e.g. Spellmann et al. 2011Spellmann et al. , 2015Kölling et al. 2009;Falk et al. 2013;Hanewinkel et al. 2014). To enhance the understanding of interactions between persistent dryer climates and forest productivity and to harness quantitative drought assessment information for forest planning we pursue three goals in the presented research: (1) We quantify the effect of persistent dryer climates on forest productivity in two north German forest regions until 2070. This is pivotal as we want to complement the currently qualitative drought vulnerability assessment by quantitative effects of climate change on forest growth. For this purpose we apply a single-tree growth simulator which already serves as a prediction tool for practical forest planning in northern Germany. For our analysis we present and employ an extended site-sensitive model version.
(2) We quantify the uncertainty inherent in climate change by applying three different climate projections.
(3) Finally, we point out how forest growth simulations of actual current forests in a specified region and for a defined projection period might add information to support forest planning decisions. This analysis serves as a first step towards a quantitative drought vulnerability assessment to be considered in forest management on a regional basis, for example for silvicultural management strategies, species selection decisions or long-term timber supply studies.

Methods
This study is based on projections of forest development applying the site-sensitive version of the single-tree growth simulator WaldPlaner 2.1 (Hansen and Nagel 2014). Site-sensitivity comprises, among others, effects of climatic conditions on tree growth, i.e. the impact of precipitation and temperature. In the simulation design we vary the climatic conditions according to three climate scenarios while all other factors are held constant. Thereby, differences in forest growth in the simulation results are traced back solely to variations in climate.

Study area
In this study we use forest inventory data from two regions Uelzen and Fläming in the north German lowlands (Fig. 1).
The two regions feature different current climatic and site conditions with favourable water supply in Uelzen and already today problematic precipitation and available water capacity in the lower parts of the Fläming region at least for less tolerant tree species (Table 1). Climate change is believed to increase water shortage in the future. Standard forest management surveys (cf. Böckmann et al. 1998) provide information on forest stands at 999 locations in Uelzen and 1008 locations in Fläming. Figure 2 displays the individual components of the Wald-Planer 2.1 simulation framework to project forest development (for equations and parameters refer to Hansen and Nagel 2014). Single growth functions and a previous version of the model system have been extensively validated (Schmidt and Hansen 2007;Vospernik et al. 2010Vospernik et al. , 2015Nagel 2013;Sprauer and Nagel 2015). One growth projection step commonly covers a 5-year period. The single tree data base stores individual tree parameters as described in the initialization data section below. An ingrowth routine adds trees to the data set and the density-dependent mortality model predicts single tree mortality. The competition Fig. 1 The two study regions Uelzen and Fläming in the north German lowlands near the cities of Hamburg and Berlin (land cover data in the overview map © ESA / ESA GlobCorine Project, led by Université Catholique de Louvain) index (c66) and the change in the competition index (c66c) due to tree removal and mortality are calculated afterwards. Basal area increment (bai) is a function of the tree parameters age, crown surface area (csa), the competition index and its change. The updated diameter at breast height (dbh) is then used to derive an updated tree height applying the site-sensitive longitudinal height-diameter model by Schmidt (2010) which utilizes climate and soil data. When applying climate scenarios in impact studies it is recommended to use average values of 30-year periods (IPCC 2013(IPCC , p. 1450. Therefore, for each 5-year growth prediction the values of mean annual temperature sum in the growing season and aridity index are calculated for a 30-year period centered around the respective projection period and applied to the height-diameter model. Furthermore, development of crown surface area (csa) depends on crown width and crown length which, in turn, are functions of tree diameter and tree height. Thus, as tree height development is predicted site-sensitively soil and climate parameters also affect diameter increment. Finally, thinning and harvesting rules can be applied.

Model description
The height-diameter model by Schmidt (2010) constitutes the site-sensitive core of the single-tree growth model. It builds on a height-diameter model by Lappi (1997) which is a modified version of the Korf function (Eq. 1).
with   (1926) in stand k in the climatic reference period 1961 to 1990; lon k , lat k : Gauss-Krueger coordinates of stand k's centroid; elev k : elevation [m] above sea level of stand k's centroid; c p xa : parameter constants (with x = 1,3 and 4, respectively) of Chapman-Richards function used to approximate the age trend of the original parameter A (see formula 1); p 2a : parameter of Chapman-Richards function used to approximate the age trend of the original parameter A (see formula 1); f xa : one-dimensional smoothing regression term (with x = 2, 3, 4 and 5, respectively) to describe non-linear effects of independent variables on the original parameter A (see formula 1); f 6a : two-dimensional smoothing regression spline to account for large scale spatial autocorrelation of the original parameter A (see formula 1); p xb : parameter of linear effects (with x = 0, 1 and 2, respectively) on the original parameter B (see formula 1).
In a first step the expected tree height value is estimated using, among others, the climate variables temperature sum in the growing season (tsum) and De Martonne's (1926) aridity index (ari) as described in the climate data section below. The model is parameterized based on data of the first two German National Forest Inventories, the forest enterprise inventory of Lower Saxony and experimental plots for Pedunculate and Sessile oak, European beech, Norway spruce, Douglas-fir and Scots pine in Germany (Eq. 2).
The second stage of the height-diameter model (Eq. 3) uses the predictors available water capacity (awc), soil nutrients (nut) and ground water class (gw) to further enhance site sensitivity (cf. Fleck et al. 2015, p. 47). A separate second stage is necessary as the soil parameters were not available throughout Germany during model development but only for Lower Saxony.
with E(H kti ): expected tree height value [m] for tree i in stand k at time t with H kti~N (exp(η kti ), σ 2 ); d EðH kti Þ : estimated expected tree height value [m] for tree i in stand k at time t using Eq. 2; awc k : available water capacity in stand k [mm]; Z k : indicator variable to define ground water impact in stand k with categories "ground water impact" and "no ground water impact"; gw k : vector of the categorical variable to define ground water impact with categories "no ground water impact (ground water table > 2.5 m)", "low to medium ground water impact (ground water table 0.95 to 2.5 m)" and "medium to strong ground water impact (ground water table < 0.95 m)"; nut k : vector of the categorical variable to define soil nutrient supply in stand k (1 = very poor, 2 = poor, 3 = medium, 4 = good, 5 = rich; including further differentiation using + andvariants); lon k , lat k : Gauss-Krueger coordinates of stand k's centroid; f 1a : one-dimensional smoothing regression spline to describe the non-linear effect of awc; f 2a : two-dimensional smoothing regression spline to account for large scale spatial autocorrelation; β 1 , β 2 : coefficient vectors. The estimated expected tree height values calculated with Eq. 2 are used as an offset in the second model stage (Eq. 3).
Both stages of the site-sensitive height-diameter model are implemented in the forest simulation framework WaldPlaner. Various studies already proved the climate-sensitivity and biological plausibility of the model (Schmidt 2010;Albert et al. 2015Albert et al. , 2016Fleck et al. 2015). In a comparative analysis bias and precision of the model proved satisfying (Albrecht et al. 2017).

Soil data
The height-diameter model requires among others soil data as predictors. For each of the 2007 locations in Uelzen and Fläming information on soil nutrients, available water capacity and a classification of groundwater availability is provided. Details are given in Albert et al. (2015Albert et al. ( , 2016. In this analysis, however, soil data is of lower interest as the comparison of tree growth focusses on the climate aspect. Nevertheless, due to the interaction between climate and soil parameters the effect of soil properties may differ among the simulation runs.

Climate data
The 5 th IPCC report introduces a new generation of climate scenarios to depict the future concentration of global warming gases in the atmosphere (IPCC 2013). Four scenarios, denominated representative concentration pathways (RCP), are selected in the report. The scenarios differ in the underlying pathways, i.e. under RCP 8.5 rising CO 2 -equivalents after 2100 (>1,370 ppm in 2100), two pathways stabilizing the CO 2 -equivalents on different levels in 2100 (~850 ppm under RCP 6.0 and~650 ppm under RCP 4.5), and a peak and decline pathway with a maximum of~490 ppm before 2100 under RCP 2.6.
Based on the climate scenario The regional climate model STARS (Orlowsky et al. 2008) further downscales the realisations of the global circulation models to Germany. In the period 2051 to 2070 the projected temperature based on INM-CM4 shows a 1.1°C deviation from today's (1991 to 2010) mean annual temperature for the north German lowlands. ECHAM6 exhibits a 1.6°C deviation and AC-CESS1.0 2.7°C. Subsequently, these three climate scenarios are labelled minimum, median and maximum climate run. Furthermore, a retrospective climate scenario from 1951 to 2010 is projected using STARS based on daily weather records from 1218 climate stations run by the German Weather Service. Thus, a consistent time series from 1951 to 2070 is provided for each climate station.
The site-sensitive height-diameter model applies the climate parameters temperature sum in the growing season and aridity index as predictors. The beginning of the growing season is defined using the model LNVAR by Menzel (1997, p. 52 ff ). The end is calculated according to suggestions by Walther and Linderholm (2006) and Frich et al. (2002) using either the temperature criterion or the short day criterion, whichever is met first. The temperature criterion applies, if the moving average of mean temperature in a 7-day period between July and October falls below 5°C. We follow von Wilpert (1990) and fix the short day criterion, defined as the day length needed for xylem growth, to October 5. De Martonne (1926) defines the aridity index as the ratio of annual precipitation sum [mm] and mean annual temperature [°C] +10. We calculate the aridity index based on annual values as the soil functions as a water storage. Finally, the temperature sum in the growing season and the aridity index is provided for all 2007 locations in Uelzen and Fläming in yearly resolution from 1951 to 2070 by interpolating the projected values of all climate stations within a 20 km distance from the respective location. The interpolation is done using WaSiM-ETH's (Schulla 2015; Schulla and Jasper 2007) distance-weighted regression model. Figure 3 illustrates how the temperature sum in the growing season and aridity index develop differently in Uelzen and Fläming over time. The linear trend lines clearly display an increase in temperature over time ( Fig. 3a and b). For the aridity index, with lower values indicating dryer conditions, the findings show increasing dryer conditions with time ( Fig. 3c and d). The temperature trend in Uelzen is on a slightly lower level than in Fläming. However, the aridity index states much dryer conditions in Fläming.
The baseline climate scenario is used as a reference in this study. It requires to avoid any temporal changes in temperature sum and aridity index looking at 20-year averages. Annual fluctuations, however, are necessary to display short-term climate variability. Therefore, the annual climate values of the reference period 1991 to 2010 are randomly mixed to create new sequences for each of the three projection periods 2011 to 2030, 2031 to 2050, and 2051 to 2070. This procedure returns identical means and variances for the three future periods.

Initialization data
Using the forest simulation framework WaldPlaner 2.1 (Hansen and Nagel 2014) 2007 forest stands with individual trees are generated based on the provided stand information, i.e. for every species and for each layer values on age, quadratic mean diameter, and volume per hectare as well as species proportions and stocking degree. Single tree heights are estimated site-sensitively applying the height-diameter model. These initial tree heights are further updated in each 5-year prediction period with newly predicted values from the heightdiameter model. As the survey data has different inventory dates all stands are projected to 2010 as a uniform starting year for the growth projections.
The effect of different climatic conditions on singletree volume increment, periodic volume increment at stand level and mean annual increment at biological rotation age of the three most abundant tree species in Uelzen and Fläming, i.e. Scots pine (Pinus sylvestris L.), European beech (Fagus sylvatica L.) and Pedunculate and Sessile oak (Quercus robur L. and Quercus petraea (Matt.) Liebl.), is analysed comparing the simulation results under three different climate projections to a baseline scenario with constant climate until 2070. For the analysis we distinguish three 20-year projection periods, i.e. the first period from 2011 to 2030, the second from 2031 to 2050 and the third from 2051 to 2070. In the course of the simulations more and more individual trees are removed by the thinning and harvesting algorithm or due to mortality and new trees are introduced by the ingrowth algorithm. Ingrowth, on the other hand, is not identical on a single tree basis among the four simulation runs. Thus, although starting with quite a number of individuals in 2011 in later periods the data base for tree-wise comparison becomes too small especially in younger age classes. Therefore, tree collectives are selected for comparison (cf. Additional file 1: Tables S1 and S2). Under each scenario all trees of one species, i.e. pine, beech or oak, being alive at the end of the respective 20-year projection period or being cut or dying because of density induced mortality during the period form one collective. As a result 9 collectives under the baseline scenario (3 periods × 3 species) are compared to 27 collectives under climate change (3 periods × 3 species × 3 climate scenarios). For each tree within the collectives the mean annual volume increment is predicted.
The collectives to be compared, i.e. per species and within one period, consist of roughly the same number of individuals, most of them being identical, and have a very similar age span. The tree dimensions vary as projected tree growth depends on the assumed climate conditions. The establishment of new stands and possible species conversions, influence the age structure of all collectives as well. Under each simulation run the same settings for the business-as-usual silvicultural treatment scenario and the same rule-sets for regeneration and forest conversion apply (cf. Albert et al. 2015).

Deriving the climatic effect on volume increment
Based on the tree collectives described above the climate effect on mean single tree volume increment [m 3 ·yr -1 ] over initial volume classes is calculated separately for tree species and projection period as the deviation between the respective climate scenario and the baseline scenario: with iv ijx , iv ijb : volume increment of tree i in initial volume class j under climate scenario x and under baseline b [m 3 ·yr -1 ]; x: climate scenario (minimum, median, maximum); n jx , n jb : number of trees in initial volume class j under climate scenario x and under baseline b.
Periodic volume increment is also calculated based on the tree collectives, however, only trees of the dominant stand are taken into account as it is common in forest planning to differentiate between layers. The deviation in periodic volume increment [m 3 ·ha -1 ·yr -1 ] between the respective climate scenario and the baseline scenario is derived separately for each of the three species and for each of the three projection periods as follows: with iv ijx , iv ijb : volume increment of tree i in stand j under climate scenario x and under baseline b [m 3 ·yr -1 ]; ba jx , ba jb : proportion of respective species basal area in stand j under climate scenario x and under baseline b; x: climate scenario (minimum, median, maximum); n kjx , n kjb : number of trees of species k in stand j under climate scenario x and under baseline b; m x , m b : number of stands under climate scenario x and under baseline b.
The climate effect on mean annual increment at biological rotation age (MAI max ) is derived only for stands which are or soon will be close to rotation age, i.e. European beech and Scots pine stands aged 81 and older and oak stands aged 121 and older at the end of the respective projection period. The rationale of this restriction is that we derive MAI max based on the simulated stand's top height using static site index assessment refraining from the dynamic influence of changing site conditions which is a permissible approximation in this case. In order to obtain MAI max -values for younger stands taking climate change into consideration we would need to project forest growth climate-sensitively for younger stands beyond 2070. However, for this study no data on climate projections beyond 2070 is available. Therefore, only stands close to rotation age warrant a meaningful interpretation of MAI max -values under climate change conditions.
The site index at age 100 (SI) is estimated using the yield tables by Nagel (2017) based on the simulated top height and age of beech, oak and pine occurrence in each stand i at the end of the respective projection period. The SI-values are further transformed into MAImax -values applying following functions: Sessile and pedunculate oak : Scots pine : Finally, the climate effect on MAI max [m 3 ·ha -1 ·yr -1 ] is calculated separately for tree species and projection period as follows: with MAI max ix , MAI max i b : mean annual increment at biological rotation age in stand i under climate scenario x and under baseline b [m 3 ·ha -1 ·yr -1 ]; x: climate scenario (minimum, median, maximum); m x , m b : number of stands under climate scenario x and under baseline b.

Results
It is the primary objective of this simulation study to quantify the effect of long-term persistent dryer climates on forest productivity and to derive tangible recommendations for regional forest planning. Thus, we first compare volume increment under climate change to volume increment under constant current climate. The analysis based on mean values from single trees provides deeper insights into individual growth reactions to climate change. The innate age trend on volume increment is eliminated by looking at deviations. However, the absolute value of volume increment depends on initial tree dimension (c.f. Berrill and O'Hara 2014). Consequently, the larger the initial volume the more likely a large deviation between volume increment under climate change and constant climate will occur.
Five general findings are apparent when analyzing the climatic effects on mean single-tree volume increment over initial volume classes (Figs. 4 and 5). (1) The expected trend of larger deviations in absolute terms in larger initital volume classes is obvious for all species in both regions except for Scots pine in Uelzen. Accompanying, the uncertainty range prominently increases over initial volume classes, i.e. the differences in volume increment deviation among the three climate scenarios are relatively minor for small trees but increase to a significant extent for larger trees.
(2) Climate change has a positive effect throughout initial volume classes in most cases in the first period. This effect diminishes or even reverses over time.
(3) While the uncertainty range includes only positive deviations in most cases in the first period, the uncertainty range indicates positive as well as negative climate impacts on volume increment throughout initial volume classes in most cases in the second and third period. This is especially apparent for Scots pine in the third period and for European beech as well as oak in the second and third period in both regions. (4) Scots pine is the least sensitive species regarding the climatic effect on volume increment followed by oak and beech being affected the strongest (please observe the dashed lines in Figs. 4 and 5 indicating different scales of the vertical axes). (5) Furthermore, the climatic effect is stronger in Fläming than in Uelzen.
Deviations in mean single-tree volume increment over initial volume classes help to analyze climate effects on tree growth on the scale of biological productivity. However, the climatic effect on volume increment at the forest stand level might differ. Thus, it is also necessary to report results of deviations in periodic volume increment. On the scale of area productivity using periodic Fig. 4 Deviation in mean single-tree volume increment for different classes of initial single-tree volume between the respective climate scenario (minimum (yellow), median (orange) and maximum (dark red)) and the base line scenario for Scots pine (a), European beech (b) and Pedunculate and sessile oak (c) in Uelzen (blurred gray area indicates the uncertainty range; *) no observations under maximum climate scenario; dashed lines emphasize the -5 dm 3 •yr -1 to 5 dm 3 •yr -1 range for better comparison due to unequal scaling of the vertical axes) volume increment more general tendencies of climatic effects on forest productivity can be identified (Fig. 6). As a reference we provide distributions on absolute values of periodic volume increment in the supplement (Additional file 2: Figures S1-S6).
The overall deviations in periodic volume increment displayed in Fig. 6 point out five general tendencies. (1) In accordance with the findings on mean single-tree volume increment periodic volume increment of all species benefits from climate change in the first projection period in both regions except for European beech and oak in Uelzen under the median and maximum scenario. (2) The uncertainty range widens over time and in most cases encompasses positive as well as negative deviations. Especially European beech strongly reacts on climate scenarios with positive deviations up to 2 m 3 ·ha -1 ·yr -1 as well as negative values down to -3 m 3 ·ha -1 ·yr -1 in the third period. (3) Basically, the positive effect of climate change in the first period diminishes over time and even turns negative in many cases. While the deviation in periodic volume increment of Scots pine shows only a few small negative values in the third period, periodic volume increment of oak is negatively affected by climate change in Uelzen entirely but not so much in Fläming. The climate effect on periodic volume increment of European beech shows a disparate picture over time with positive tendencies under minimum and median climate scenario in Uelzen and a clearly negative trend under median and maximum climate scenario in Fläming.  Deviation in mean single-tree volume increment for different classes of initial single-tree volume between the respective climate scenario (minimum (yellow), median (orange) and maximum (dark red)) and the base line scenario for Scots pine (a), European beech (b) and Pedunculate and sessile oak (c) in Flaeming (blurred gray area indicates the uncertainty range; *) no observations under median and maximum climate scenario; dashed lines emphasize the -5 dm 3 •yr -1 to 5 dm 3 •yr -1 range for better comparison due to unequal scaling of the vertical axes) noticeable in Fläming looking at Scots pine and European beech. In Fläming the periodic volume increment of both species shows more positive reactions in the first period and under median and maximum scenario significant negative deviations in the third period. Periodic volume increment of pedunculate and sessile oak, however, is affected more heavily by climate change in Uelzen.
While the results displayed in Figs. 4, 5 and 6 emphasize the immediate short-term climatic effect on volume increment during the 20-year projection periods, the analysis of mean annual increment at biological rotation age (MAI max ) detects the cumulative long-term consequences of climate change on forest productivity. The differences in MAI max indicate how the gradual divergences in volume growth between the baseline and the three climate scenarios accumulate to economically meaningful additional or fewer productivity, respectively (for absolute values on MAI max please refer to Additional file 2: Figures S1-S6).
In Uelzen Scots pine and oak both benefit with additional MAI max of up to 0.4 m 3 ·ha -1 ·yr -1 depending on climate scenario and projection period (Fig. 7a). The projection only indicates a marginal loss for oak under maximum climate scenario in the third period. In Fläming the deviation in MAI max shows similar patterns for pine and oak, however with slightly larger uncertainty ranges comprising larger positive values up to 0.6 m 3 ·ha -1 ·yr -1 for pine (0.25 m 3 ·ha -1 ·yr -1 for oak) and negative values down to -0.2 m 3 ·ha -1 ·yr -1 for both species in the third period (Fig. 7b). While the uncertainty range increases significantly for both species in both regions over time, there is no obvious trend in productivity with increasing dyrer conditions until 2070. Contrary for European beech, the losses in MAI max clearly increase over time. While beech productivity benefits from climate change in the first period in most cases with deviations in MAI max up to 0.7 m 3 ·ha -1 ·yr -1 in Uelzen and Fläming, there is a strong negative effect in Uelzen of down to -1 m 3 ·ha -1 ·yr -1 in the second and even down to -3 m 3 ·ha -1 ·yr -1 in the third projection period. In Fläming a similar trend for beech is observed with negative values down to -0.8 m 3 ·ha -1 ·yr -1 in the second and down to -2 m 3 ·ha -1 ·yr -1 in the third period. However, in Fläming under minimum climate scenario positive Fig. 6 Overall deviation in periodic volume increment between the respective climate scenario (minimum (yellow), median (orange) and maximum (dark red)) and the base line scenario for Scots pine, European beech and Pedunculate and sessile oak in Uelzen (a) and Flaeming (b) (blurred gray area indicates the uncertainty range) deviations in MAI max of 0.4 m 3 ·ha -1 ·yr -1 in the second period and even of 0.8 m 3 ·ha -1 ·yr -1 in the third period are projected.

Discussion
Generally, simulation studies serve to evaluate future developments in the light of differing framework conditions (Clark et al. 2001). The main objective of simulation studies is to provide information beyond measurements and observations for forest planning and to support decision making. In contrast, resource assessments, field experiments and observational studies collect data and record forest conditions in the first place (Zhao et al. 2014). The measured data and the results and conclusions derived thereof are used for policy making, forest planning and resource allocation, hypothesis testing and modelling. Therefore, simulation studies on possible future developments and measurement data from different sources with a retrospective view may complement each other especially for forest planning issues under uncertainty.
The presented simulation study relies on current forest survey data which is then projected to assess forest development under climate change until 2070. The analysis of persistent dryer climate on forest growth and forest productivity reveals insights on three major issues.
(1) As a prerequisite for any decision support the effect of climate change on forest growth has to be quantified. Recommendations for adaptation measures cannot be derived before potential deviations between the current growth level and the projected growth level under climate change are assessed. Furthermore, it is crucial for reliable forest planning to take the uncertainty in projections into account. In view of the practical application of this simulation study the deviation in volume increment between constant climate and three climate scenarios for three species, Scots pine, European beech and Pedunculate and Sessile oak, is evaluated in two forest regions. Thus, the results show variations due to the variability in the initial forest and site conditions as well as the variability in the climate projections. The immediate, short-term effect of climate change on volume increment gives evidence of the growth model's sensitivity to the opposing climate factors, i.e. the positive influence of rising temperatures on growth and the inhibitory effect of reduced precipitation. In the near future of the first projection period the positive effect of rising temperatures ensures higher periodic volume increments of up to 2 m 3 ·ha -1 ·yr -1 depending on species, region and climate scenario (cf. Fig. 6). However, in the projection period 2051 to 2070 the effect of persistent dryer conditions finally overcompensates the positive effect of higher temperatures in most cases in Uelzen and Fläming with losses up to -3 m 3 ·ha -1 ·yr -1 . This quantitative comparison of periodic volume increment under constant and climate change conditions allows to evaluate species' productivity and future sitesuitability beyond a risk assessment alone. Specifically, in Uelzen beech will most likely suffer significant decreases in volume increment. Thus, persistent dryer conditions will strongly affect managing beech as a leading species on many sites. In the Fläming region, however, where beech predominantly serves as an admixed species already today (Albert et al. 2017), projected decreases in volume increment will be tolerable. On the other hand, Scots pine in pure and mixed stands dominates in both regions with proportions of almost 70% of the entire forest area in Uelzen and even 90% in Fläming. Therefore, even small decreases in volume increment will amount to considerable losses in the regions. Similar to beech, oak only plays a minor role in both regions looking at the current proportions of 5% in Uelzen and 4% in Fläming. However, even small potential decreases in volume increment might have a considerable impact on oak management in both regions as the economic value per cubic meter oak wood is high. While the analysis of single-tree and periodic volume increment in the three projection periods reveals the immediate impact of changing climatic conditions on forest growth, the evaluation of changes in MAI max aims at the long-term, aggregating climatic effect on forest productivity. Basically, the results on periodic volume increment and MAI max are in accordance. On average the productivity of all pine stands close to rotation age will profit from climate change until 2070, with the exception of the MAI max -value under maximum climate scenario in Fläming in the third projection period. Consequently, looking at overall productivity climate change does not negatively affect the management of Scots pine stands being harvested in Uelzen and Fläming until 2070. Contrary, persistent dryer climates heavily affect beech management with MAI max losses up to -3 m 3 ·ha -1 ·yr -1 in Uelzen and -2 m 3 ·ha -1 ·yr -1 in Fläming. Especially in stands with beech as the leading species it might be necessary to adapt the management goals, e.g. to reduce the target diameter in order not to unduly extend rotation age. As a summary, the effect of climate change on the growth of Scots pine and European beech is quite similar looking at periodic volume increment and MAI max . However, small distinctions are obvious for oak. The reported losses in volume increment for oak in the second and third period in both regions are not as pronounced looking at MAI max . Especially in comparison to beech MAI max -values for oak indicate only minor losses (Fig. 7) while the negative deviations in periodic volume increment in the third period especially in Uelzen and under maximum climate scenario in Fläming seem considerable (Fig. 6). In the case of Uelzen, we explain this discrepancy with the cultivation effort of oak even on less favourable sites in the beginning of the 1970s. These oak stands seem to suffer under dryer conditions, thus, the deviations in periodic volume increment are negative. However, these stands are not included in the data used for MAI max assessment because of the age restriction. On the other hand, the older oak stands mainly stock on sites which are not so much affected by dryer conditions, thus, deviations in MAI max are positive.
(2) The second major issue of the presented analysis is to improve drought stress assessments with quantitative data on forest productivity. One common method to account for potential drought stress limitations when evaluating species for natural or artificial regeneration in northwestern Germany is a species-specific vulnerability assessment based on site water budget according to Grier and Running (1977), i.e. the sum of climatic water balance in the growing season (difference of precipitation and potential evapotranspiration for grass reference (Allen et al. 1998)) and the available water capacity (Spellmann et al. 2007. Generally, vulnerability is defined as the extent to which a system is susceptible to disturbances, or in our case how a forest ecosystem is affected by climate change (Füssel and Klein 2006). Specifically, drought vulnerability is related to exposition, i.e. duration and severity of water shortage, species' resilience against fluctuations in water availability and species' ability to adapt to climate-induced drought. Spellmann et al. (2007 assess drought stress by transferring ranges of site water budget into species-specific vulnerability classes ( Table 2). The reasoning behind this method is that different tree species need a specific amount of water supply to achieve an economically advisable cultivation and management. And forest management is only economically promising, if the tree species show sufficient growth potential and potential risks do not exceed a certain level. The values in Table 2 were derived from observations on species distribution ranges in combination with a qualitative assessment of cultivation risks. The cultivation risks comprise effects on growth, vitality and mortality. The values classify mean site conditions, i.e. not accounting for influencing factors such as topographic site exposition. The values indicate Norway spruce as most prone to drought vulnerability, followed by European beech, oak and Douglas-fir and Scots pine as most resistant. An example illustrates how results from a drought vulnerability assessment following the method proposed by Spellmann et al. (2007 can valuably be amended by quantitative information on future species productivity. In a former study, we evaluated the necessity of adapting tree species selection in Uelzen and Fläming taking drought restrictions into account (Albert et al. 2017). For Uelzen we concluded that the number of Scots pine dominated stands will increase from 60% refraining from drought limitations to 69% under medium climate scenario and even 80% under maximum climate scenario until 2070 (under minimum climate scenario the proportion is estimated to be constant). This strong preference of Scots pine in future tree species selection based simply on drought vulnerability limits is in line with the current analysis that no severe losses neither for periodic volume increment (Fig. 6) nor for MAI max (Fig. 7) are to be expected. For oak, however, a recommended slight increase from 7% under the assumption of sufficient water supply to 8% under maximum climate scenario might be challenged looking at the projected future productivity. If the productivity goal is rated high in the decision process, potential losses in periodic volume increment of up to -2.5 m 3 ·ha -1 ·yr -1 might influence the decision process on species selection towards more productive species under climate change although oak qualifies as a species to be cultivated under the drought vulnerability limits. Evaluating European beech requires a different line of reasoning. In our species selection study we expected the proportion of beech dominated stands in Uelzen to drop from 16% with no assumed drought limitations to 14% under minimum, and an insignificant 2% under median and maximum climate scenario. While these recommendations are based solely on the intrinsic and fuzzy approach transferring expert-knowledge of species' productivity and risks into drought vulnerability classes we complement and confirm these results with sound information on projected substantial losses for beech looking at periodic volume increment and MAI max under climate change in Uelzen.
In summary, as forest growth is projected as a function of several climatic and soil properties the derived results of the presented study complement the drought vulnerability assessment of Table 2 beyond evaluating the site water budget alone. The quantitative statement on potential productivity losses, including an uncertainty range which depends on the assumed climate projections and the time horizon, enriches the information basis for forest planning decisions in the light of persistent dryer climates.
(3) It is another objective to evaluate the approach's potential to derive conclusions for forest management and to formulate adaptation strategies using the two regions Uelzen and Fläming as an example. First of all, our results correspond to the perception by Kätzel and Höppner (2011) who conclude that Scots pine will persist as an important and economically meaningful species for forest management in the northeastern German lowlands and that oak and beech will show different reactions to climate change depending on other environmental factors. Our results also confirm findings by Bauwe et al. (2015) who predicted the radial growth response to climate change of Scots pine, European beech and Pedunculate oak along a precipitation gradient in northeastern Germany. As in our analysis growth rates for beech and oak were projected to decrease under the A1B climate scenario until 2100 while Scots pine showed slight growth increase as well as decrease depending on site and location (as a reference, A1B is best comparable to RCP 6.0). A broad study on sustainable land use management by Spellmann et al. (2017) specifically recommends to promote climate-adapted tree species and species mixtures for future forest generations in order to lower and distribute potential future risks in Uelzen. Furthermore, the already completed advance planting of beech should be enriched with natural regeneration of more drought tolerant species such as oak, Scots pine and birch or even by planting Douglas-fir, Grand fir or red oak in beech gaps. Our analysis clearly supports the recommendations by Spellmann et al. to promote Scots pine and critically question the future role of beech in Uelzen. However, our results on potential productivity losses in the case of oak at least relativize its importance for forest management in Uelzen. In any case, it is crucial in forest planning to implement tree species selection for future stands as a case-by-case decision. In the same study, Spellmann et al. conclude for the Fläming region that the increasingly dryer conditions will also negatively affect Scots pine. Our results on periodic volume increment seem to support this assessment. However, the silvicultural scope in tree species selection is limited in Fläming, thus, it is most advisable to promote mixed forests, i.e. mainly leading Scots pine or Douglas-fir with admixtures of beech or oak.

Assumptions and model limitations
In this study the interaction between the site-sensitive longitudinal height-diameter model and the basal area increment model as well as other functionalities, e.g. mortality and ingrowth, within the framework of the growth simulator WaldPlaner plays a major role. The interpretation of our simulation results is subject to various assumptions and model limitations.
(1) First of all, the results are only valid for conditions specified by the applied growth functions, i.e. the model functionalities and the parameterization basis, and in the range of the assumed climate projections. The site-sensitive height-diameter model is parameterized for the entire area of Germany. The other growth functions of the WaldPlaner framework are parameterized for northwestern Germany. Thus, the spatial scope in this study is well covered. Furthermore, the height-diameter model is a statistical model fit on past observations. This implies two possible pitfalls. First, under climate change the parameter values might well be in the extrapolation range with unsuspected model behavior. However, biologically feasible results were obtained in various sensitivity studies on the model (Schmidt 2010;Albert and Schmidt 2014). Second, the model philosophy assumes that the significance of the selected model parameters has the same explanatory power under future climatic conditions. However, under more extreme climate conditions other factors could strongly influence forest growth. For example, Michelot et al. (2012) compared growth reactions of European beech, sessile oak and Scots pine to climatic variations and soil water deficits in France. Their findings clearly show that growth of beech and Scots pine will be determined by mainly summer maximum temperatures and precipitation in the growing season also under future climate conditions, whereas oak growth will be strongly affected by autumn droughts due to tree carbohydrate storage depletion. In this case, with increasing autumn drought stress and a higher frequency of extreme drought events during that season in the future a new parameter, e.g. autumn water deficit, would presumably be significant in modelling growth. These possible shifts on parameter significance and interactions under climate change are beyond the capabilities of statistical model approaches. But we assume that the broad range of observed climate values and especially different combinations on temperature and aridity index in the past throughout Germany covers or is close to projected future conditions for the large part.
The simulation results strongly rely on the climate projections, which in turn depend on the applied regional climate model STARS (Orlowsky et al. 2008). Independent from the three climate scenarios resulting from three different global circulation models (GCM), STARS always projects climate patterns with dryer summers and wetter winters in Germany under any scenario. These climate patterns stimulate the presented impact analysis. These uniform climate patterns originate from a specific temperature-driven resampling algorithm (TCR) being used for climate projections by the statistical model STARS. The GCMs define the temperature forcing as a periodic trend. TCR generates a new sequence of annual climate values according to the temperature forcing in a first step. In a second step 12-day-weather sequences are newly arranged. Wechsung (2014, 2015) criticized the underlying assumption by Gerstengarbe et al. (2013) to validate STARS based on mean statistical estimates during the period 1951 to 1975 in comparison to a new period with newly arranged values using TCR. Moreover, Wechsung (2014, 2015) prove a negative bias in STARS' projections of precipitation of 5% per 1°C temperature increase. Nevertheless, accepting the model limitations and regarding the type of climate change projection Wechsung and Wechsung (2014) recommend STARS to be used in climate impact studies. Therefore, we employ STARS' climate projections in the presented impact study in lack of alternative available data, knowing about potential pitfalls.
(2) Conclusions on the results are only permitted under the assumption that the scope of future climate conditions is within the quantified uncertainty range of the climate projections. The knowledge on the likelihood of a certain climatic development is low. Therefore, probabilities cannot be assigned whether or not certain climate conditions will be realized (Foley 2010). As a consequence scenario analysis employing climate ensembles is utilized to quantify a range of possible outcomes. This range defines the uncertainty in climate modelling. However, this approach falls short to cover the uncertainty in its entirety as beyond scenario analysis which assumes certain boundary conditions there is a considerable amount of uncertainty left due to the lack of any knowledge that originates from the long time scales and complexity, termed deep uncertainty (Walker et al. 2003;Kandlikar et al. 2005). Knutti and Sedláček 2012 reported on the robustness and uncertainties in the Coupled Model Intercomparison Project Phase 5 (CMIP5) by the World Climate Research Programme. They conclude that in spite of tremendous improvements of the climate models and higher spatial resolutions the projection uncertainties did not decrease compared to successive IPCC reports. Hence, Knutti and Sedláček (2012) advise to account for the uncertainties when deriving decision support from climate impact analysis. Our reported analysis reflects climate projection uncertainties by applying three climate scenarios which is a compromise between adequate computational workload and sufficient consideration of possible projection uncertainties. Informed recommendations for forest planning thus are possible, however, any conclusion derived from impact analysis must remain highly conditional on the underlying ranges of climate projections. But at best, conclusions drawn are robust, e.g. in Uelzen under each of the three applied climate scenarios the deviation in MAI max indicates a consistent positive climate effect for Scots pine and, contrary for European beech, a consistent negative effect in the third period (Fig.  7). Therefore, any decision in favor of Scots pine management will result in higher productivity. And the other way around, any decision to the disadvantage of beech management, e.g. turning away from the currently common practice of converting Scots pine stands into stand types with leading beech or even an earlier conversion of current beech stands to more productive species, are justified by losses in MAI max .
(3) The dynamics in climate parameters are accounted for in the height development, although the height-diameter model is strictly speaking static. The climate effect on basal area growth, however, is indirect as the predictor crown surface area is estimated based on tree height. As analyzed in this study, the growth model is capable to reproduce effects of persistent dryer climate conditions on forest productivity in a long-term perspective. One shortcoming of the growth model is that disturbances such as extreme short-term drought events are neglected. The impact of severe dry spells on tree growth was confirmed in numerous studies (e.g. Orwig and Abrams 1997;Klos et al. 2009;Michelot et al. 2012;Weber et al. 2013). In summary, both the growth reactions to long-term trends and the consequences of disturbances are important factors to adequately project forest development. Lindner et al. (2014) stated in this context that not accounting for extreme events might explain the deviations between observed growth trends and growth projections in many cases. Nevertheless, decision support for forest planning and specifically in the assessment of species suitability can be granted based on projections of species productivity as shown in the impact study in combination with a sound vulnerability analysis accounting for possible disturbances, for instance an increasing frequency of extreme drought events. Ecophysiological processes such as interactions between competition, species mixture, pathogens, nitrogen deposition or elevated CO 2 levels with drought stress are also not accounted for in the applied growth model. Evidence for these interactions are provided by several studies. For example, Giuggiola et al. (2013) reported for Scots pine on xeric sites in Switzerland that a strong reduction in basal area of 40% to 60% of the regular management regime mitigated drought effects. On the issue of species mixture Arthur and Dech (2016) contributed further findings that species mixtures have a significant influence on the effects of drought on forest growth. The hypothesis that species composition can mitigate drought effects was also supported by Pretzsch et al. (2013). On the other hand, Thurm et al. (2016) showed for Douglas-fir and European beech mixed stands that not all species equally profit from mixtures looking at drought stress release. In their study Thurm et al. (2016) provided evidence that in mixture Douglasfir is positively affected while beech shows a higher sensitivity to climate impacts. Besides the direct interaction of tree or stand properties on drought stress also external factors may alter the impact of drought on forest growth. A strong interaction between drought and pathogens and, eventually, with tree vitality, growth and even mortality exists (e.g. Desprez-Loustau et al. 2006;Overbeck and Schmidt 2012;Schoneberg 2017). There is still a great uncertainty about the direction and the magnitude in the response of insects and fungi to drought (Jactel et al. 2011). But as drought causes physiological changes in tree defenses it will alter the patterns and effects of biotic disturbances (Ayres and Lombardero 2000). Consequently, the interaction of drought and pathogens will affect tree growth and tree vitality (Dobbertin 2005).
Increased nitrogen deposition has a fertilizing effect on forest growth (Hyvönen et al. 2007;Kahle et al. 2008a). Then again, the positive effect of nitrogen does not apply to all sites equally (e.g. Näsholm et al. 2000). Furthermore, some studies suggest that nitrogen fertilization might have a negative effect on root growth and might increase root mortality, thus reducing tree vitality (Clemensson-Lindell and Persson 1995aPersson , 1995bPersson and Ahlström 2002 Recently, simulation studies indicated the high sensitivity of forest growth to assumed future CO 2 increase Lindner et al. 2014;Schelhaas et al. 2015). The assessment of interactions between elevated CO 2 -concentration and drought is even more sophisticated. Hättenschwiler et al. (1997), Ellsworth (1999) and Körner (2006) reported a possibly beneficial impact of elevated CO 2 on tree growth under drought conditions. Some studies indicated a higher wateruse efficiency of trees as the main reason for a positive effect of higher CO 2 under drought stress (e.g. Battipaglia et al. 2013;Keenan et al. 2013). Contrary, Tognetti et al. (2000) found no evidence for growth response to increasing CO 2 during drought stress for five Mediterranean forest tree species. All the above mentioned interactions could alter tree growth under extreme drought events and even under a persistent dryer climate. More research is needed to include these effects in growth models to improve informed decision support for forest management. Ciais et al. (2005) clarified the extent of the implications as they assumed that extreme drought events and persistent dryer conditions might significantly alter the long-term carbon balances and reverse forests from sinks to sources in Europe.

Conclusions
This study illustrates a method to assess the long-term effects of persistent dryer climate conditions on forest growth and productivity by comparing mean single-tree volume increment, periodic volume increment and mean annual increment at biological rotation age (MAI max ) under constant climate conditions and under three climate projections until 2070. The projected changes in forest growth within a specified uncertainty range serve as quantitative contributions to provide decision support in forest planning, e.g. adapting silvicultural management strategies, evaluating species future site suitability and assessing the impact on timber supply. The quantitative analysis of possible productivity changes under persistent dryer climate complements the expertknowledge based drought vulnerability assessment which is applied in northwestern Germany today (cf. Table 2). The drought vulnerability assessment as a synoptic approach qualitatively aggregates the effects of drought on tree growth, vitality and mortality based on thresholds for site water budget. The reported projection of productivity changes quantifies effects of persistent dryer climates on forest growth separately. In addition, the effects of drought on vitality and mortality still needs to be quantified on a level suitable for forest management.
In particular, on the grounds of the specific simulation results following conclusions for forest planning in the two regions can be drawn: (1) Although the uncertainty range in the third period (2051 to 2070) encompasses positive as well as negative effects on volume increment in many cases for European beech, pedunculated and sessile oak and Scots pine in both regions decisions on species suitability have to be made. In this case, differences in site conditions need to be taken into account and furthermore the decision also is a question of risk attitude. The results presented are average values for all sites considered. Thus, favourable site conditions, e.g. sufficient available water capacity, might well compensate the negative effect of persistent dryer climates and different recommendations for these sites apply.
(2) Scots pine is rated as the least sensitive species to dryer climate conditions. MAI max of Scots pine even increases in the time period and regions under consideration. Thus, from the growth and yield point of view Scots pine will be well suited for forest management in the future. The same but not quite as pronounced applies to the management of oak in the Fläming region while being doubtful in Uelzen.
(3) European beech faces the most severe growth reductions. In stands with beech as the leading species the negative impact of climate change might call for silvicultural adaptation measures and a change in species in the next forest generation. On sites, where a vulnerability analysis does not exclude beech from cultivation, forest management with beech as an admixture species is still advisable at the cost of lower productivity. (4) The clear temporal trend in magnitude and direction of the effect of persistent dryer climates on forest growth and productivity has to be considered.

Additional files
Additional file 1: Table S1. Characteristic numbers for the collectives under baseline (base) and minimum (min), median (med) and maximum (max) climate scenario for Uelzen at the end of the respective period (for age, dbh and height 1 st quantile, median (med) and 3 rd quantile values are given). Table S2. Characteristic numbers for the collectives under baseline (base) and minimum (min), median (med) and maximum (max) climate scenario for Fläming at the end of the respective period (for age, dbh and height 1 st quantile, median (med) and 3 rd quantile values are given). (DOCX 28 kb) Additional file 2: Figure S1. Boxplots of periodic volume increment (upper row) and mean annual increment at biological rotation age (lower row) under baseline (gray), minimum (yellow), median (orange) and maximum (red) climate scenario for European beech in Uelzen (included are only stands with beech as a leading species). Figure S2. Boxplots of periodic volume increment (upper row) and mean annual increment at biological rotation age (lower row) under baseline (gray), minimum (yellow), median (orange) and maximum (red) climate scenario for Pedunculate and sessile oak in Uelzen (included are only stands with oak as a leading species). Figure S3. Boxplots of periodic volume increment (upper row) and mean annual increment at biological rotation age (lower row) under baseline (gray), minimum (yellow), median (orange) and maximum (red) climate scenario for Scots pine in Uelzen (included are only stands with Scots pine as a leading species). Figure  S4. Boxplots of periodic volume increment (upper row) and mean annual increment at biological rotation age (lower row) under baseline (gray), minimum (yellow), median (orange) and maximum (red) climate scenario for European beech in Fläming (included are only stands with beech as a leading species). Figure S5. Boxplots of periodic volume increment (upper row) and mean annual increment at biological rotation age (lower row) under baseline (gray), minimum (yellow), median (orange) and maximum (red) climate scenario for Pedunculate and sessile oak in Fläming (included are only stands with oak as a leading species). Figure  S6. Boxplots of periodic volume increment (upper row) and mean annual increment at biological rotation age (lower row) under baseline (gray), minimum (yellow), median (orange) and maximum (red) climate scenario for Scots pine in Fläming (included are only stands with Scots pine as a leading species). (ZIP 1002 kb) Abbreviations ari: Aridity index; awc: Available water capacity; bai: Basal area increment; c66: Competition index at 66% of tree height from top; c66c: Change in competition index c66; csa: Crown surface area; cwb: Climatic water balance; dbh: Diameter at breast height; GCM: Global circulation model; gs: Growing season; gw: Groundwater class; iV: Periodic volume increment; iv: Volume increment of single tree; MAI max : Mean annual increment at biological rotation age; nut: Soil nutrient class; RCP: Representative concentration pathway; SI: Site index at reference age 100; TCR: Temperature-driven resampling algorithm; tsum: Temperature sum