Jukuri, open repository of the Natural Resources Institute Finland (Luke) All material supplied via Jukuri is protected by copyright and other intellectual property rights. Duplication or sale, in electronic or print form, of any part of the repository collections is prohibited. Making electronic or print copies of the material is permitted only for your own personal use or for educational purposes. For other purposes, this article may be used in accordance with the publisher’s terms. There may be differences between this version and the publisher’s version. You are advised to cite the publisher’s version. This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): Mladen Ognjenović, Ivan Seletković, Nenad Potočić, Mia Marušić, Melita P. Tadić, Mathieu Jonard, Pasi Rautio, Volkmar Timmermann, Lucija Lovreškov and Damir Ugarković. Title: Defoliation Change of European Beech (Fagus sylvatica L.) Depends on Previous Year Drought Year: 2022 Version: Published version Copyright: The Author(s) 2022 Rights: CC BY 4.0 Rights url: http://creativecommons.org/licenses/by/4.0/ Please cite the original version: Ognjenović, M.; Seletković, I.; Potočić, N.; Marušić, M.; Tadić, M.P.; Jonard, M.; Rautio, P.; Timmermann, V.; Lovreškov, L.; Ugarković, D. Defoliation Change of European Beech (Fagus sylvatica L.) Depends on Previous Year Drought. Plants 2022, 11, 730. https://doi.org/10.3390/plants11060730. ���������� ������� Citation: Ognjenović, M.; Seletković, I.; Potočić, N.; Marušić, M.; Tadić, M.P.; Jonard, M.; Rautio, P.; Timmermann, V.; Lovreškov, L.; Ugarković, D. Defoliation Change of European Beech (Fagus sylvatica L.) Depends on Previous Year Drought. Plants 2022, 11, 730. https://doi.org/ 10.3390/plants11060730 Academic Editors: Frank M. Thomas, Nenad Potočić and Oleg Chertov Received: 25 January 2022 Accepted: 7 March 2022 Published: 9 March 2022 Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affil- iations. Copyright: © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). plants Article Defoliation Change of European Beech (Fagus sylvatica L.) Depends on Previous Year Drought Mladen Ognjenović 1 , Ivan Seletković 1, Nenad Potočić 1,* , Mia Marušić 1, Melita Perčec Tadić 2, Mathieu Jonard 3, Pasi Rautio 4 , Volkmar Timmermann 5, Lucija Lovreškov 1 and Damir Ugarković 6 1 Division for Forest Ecology, Croatian Forest Research Institute, 10450 Jastrebarsko, Croatia; mladeno@sumins.hr (M.O.); ivans@sumins.hr (I.S.); miam@sumins.hr (M.M.); lucijal@sumins.hr (L.L.) 2 Meteorological and Hydrological Service, 10000 Zagreb, Croatia; melita.percec.tadic@cirus.dhz.hr 3 Earth and Life Institute, Université Catholique de Louvain, 1348 Louvain-La-Neuve, Belgium; mathieu.jonard@uclouvain.be 4 Natural Resources Institute Finland, 00790 Helsinki, Finland; pasi.rautio@luke.fi 5 Norwegian Institute of Bioeconomy Research, 1433 Ås, Norway; volkmar.timmermann@nibio.no 6 Faculty of Forestry and Wood Technology, University of Zagreb, 10000 Zagreb, Croatia; dugarkovic@sumfak.hr * Correspondence: nenadp@sumins.hr Abstract: European beech (Fagus sylvatica L.) forests provide multiple essential ecosystem goods and services. The projected climatic conditions for the current century will significantly affect the vitality of European beech. The expected impact of climate change on forest ecosystems will be potentially stronger in southeast Europe than on the rest of the continent. Therefore, our aim was to use the long-term monitoring data of crown vitality indicators in Croatia to identify long-term trends, and to investigate the influence of current and previous year climate conditions and available site factors using defoliation (DEF) and defoliation change (∆DEF) as response variables. The results reveal an increasing trend of DEF during the study period from 1996 to 2017. In contrast, no significant trend in annual ∆DEF was observed. The applied linear mixed effects models indicate a very strong influence of previous year drought on ∆DEF, while climate conditions have a weak or insignificant effect on DEF. The results suggest that site factors explain 25 to 30% DEF variance, while similar values of conditional and marginal R2 show a uniform influence of drought on ∆DEF. These results suggest that DEF represents the accumulated impact of location-specific stressful environmental conditions on tree vitality, while ∆DEF reflects intense stress and represents the current or recent status of tree vitality that could be more appropriate for analysing the effect of climate conditions on forest trees. Keywords: defoliation; monitoring; tree vitality; drought; climate change 1. Introduction Climate conditions influence the structure and function of forest ecosystems, and play an essential role in forest health [1,2]. Global warming has indisputably caused climate change, which is a significant threat to forest ecosystems [3]. The effects of climate change are generally expected to reduce tree growth and survival, predispose forests to disturbances, and ultimately change forest structure and composition at the landscape scale [4]. Therefore, there is an increasing concern in Europe over the sustainability of forest ecosystems under climate change [5]. Although vitality is a theoretical concept, it can be defined as the ability of a tree to assimilate, to survive stress, to react to changing conditions, and to reproduce [6]. As vitality cannot be measured directly, various indicators can be used to describe it [7]. Crown defoliation is a commonly used tree vitality indicator [8–10], which can be obtained cost- effectively and relatively quickly in field surveys [11]. Defoliation is defined as leaf loss in the assessable crown, as compared to a reference tree, and is observed regardless of the Plants 2022, 11, 730. https://doi.org/10.3390/plants11060730 https://www.mdpi.com/journal/plants https://doi.org/10.3390/plants11060730 https://doi.org/10.3390/plants11060730 https://creativecommons.org/ https://creativecommons.org/licenses/by/4.0/ https://creativecommons.org/licenses/by/4.0/ https://www.mdpi.com/journal/plants https://www.mdpi.com https://orcid.org/0000-0002-1335-6423 https://orcid.org/0000-0002-5581-2976 https://orcid.org/0000-0003-0559-7531 https://orcid.org/0000-0002-5192-9974 https://doi.org/10.3390/plants11060730 https://www.mdpi.com/journal/plants https://www.mdpi.com/article/10.3390/plants11060730?type=check_update&version=1 Plants 2022, 11, 730 2 of 13 cause of foliage loss [12]. Landmann [13] states that defoliation is an indicator of acute stress and subsequent recovery of forest ecosystems. However, defoliation has been criticized due to the subjectivity of the assessment, as well as it being a non-specific indicator affected by several biotic and abiotic factors [14–16]. To ensure data quality, training courses and repeated control assessments are regularly carried out on a national [17–19] and international level [20]. European beech (Fagus sylvatica L.) is a dominant broadleaved tree species in European forests that forms forest communities over a broad range of habitat conditions [21]. These forests provide multiple ecosystem goods and services [22]. Despite being adapted to a wide range of environmental conditions, the projected effects of climate change, particularly drought, will significantly affect the vitality of European beech [23,24]. Southeast Europe represents one of the most vulnerable regions with expected intensi- fication of severity and duration of droughts and heat waves. As the impacts of climate change on forests in southeast Europe will be potentially stronger and faster than on the rest of the continent [25,26], this region is an ideal model for studying the impacts of changing climatic conditions. A trend of decreasing precipitation and increasing temperatures has al- ready been observed in Croatia [27,28]. In the decade 2001–2010 alone, four drought events occurred, while only 13 took place between 1961 and 2010 [29]. In the future, the climate in Croatia is expected to be hotter and drier, with considerable impacts to be expected for the forest ecosystems. Consequently, continuous long-term forest monitoring is crucial in order to measure and assess these impacts and their consequences on ecosystem functioning. In Europe, the International Co-operative Program on Assessment and Monitoring of Air Pollution Effects on Forests (ICP Forests) is the most comprehensive European program for the large-scale assessment of forest ecosystem health [30]. The defoliation data obtained from the ICP Forests monitoring network have led to the publication of numerous studies of climate influence in several European countries, such as Switzerland [31], Germany [32], France [33], and Spain [34]. Depending on the investigated region, different climate parameters were found to have a negative impact on crown defoliation. Studies of this kind in Croatia have previously only been regional, which has limited the applicability of results [35,36]. A recent pilot study found a pronounced lag effect of both temperature and precipitation on beech defoliation [37]. Based on these studies, we hypothesize that previous and current year droughts as a consequence of high temperatures and low precipitation contribute the most to beech defoliation across Croatia. Given the ecological and economic importance of beech, it is necessary to understand the impact of climate change on beech defoliation. Therefore, by using the ICP Forests monitoring network in Croatia, our aim was to (i) identify European beech long-term defoliation trends, and to (ii) investigate the influence of current and previous year climate conditions and various site factors on beech defoliation. 2. Results 2.1. Temporal Trends in Tree Vitality There was a significant trend of increasing mean defoliation (DEF) by 0.39% annually over the study period (Figure 1). Annual overall mean defoliation values peaked in 2001 and 2014. The latter year is marked by the highest observed overall mean defoliation of 18.4%. In contrast, no significant trend in annual defoliation change (∆DEF) was observed. Time series of annual defoliation change values exhibit stationarity (Dickey–Fuller = −5.81, p < 0.01), and generally stay close to the neutral trend line. Plants 2022, 11, 730 3 of 13Plants 2022, 11, x FOR PEER REVIEW 3 of 13 Figure 1. Overall trend of crown vitality parameters from 1996–2017. (a) Defoliation trend of European beech (Tau = 0.78, Sen’s slope = 0.39, p < 0.001, orange line) and annual overall mean defoliation (DEFi, black line). (b) Defoliation change trend of European beech (Tau = −0.05, Sen’s slope = −0.01, p = 0.73, orange line) and annual overall change in mean defoliation (ΔDEFi, black line). Points represent annual plot mean values. Plot-scale trend analysis of DEF showed that 53% of plots have a significant and increasing trend of defoliation at an annual rate ranging from 0.25% to 1.29%. On the remaining 47% of plots, we did not observe any significant trend. The results of the q- statistics and Morans’ I spatial autocorrelation coefficient did not indicate spatial stratification of plot-scale DEF trends (data not shown). We did not detect any significant ΔDEF trend at the plot-scale. Therefore, the subsequent spatial stratification tests were not conducted. 2.2. Influence of Environmental Conditions on Tree Vitality During the selection procedure for DEF and ΔDEF models, additional variables were considered: stand age, site factors (altitude and orientation), and soil properties (soil pH, total nitrogen content in the soil, content of available phosphorus and available potassium in the soil) were tested. However, only altitude showed a significant impact. Although the influence of stand age was not significant, it was retained in the selection process due to improved model performance and the reported influence of stand age on defoliation in other studies [38,39]. The two linear mixed effects models (LMM) used to assess the impact of climatic variables on DEF revealed different influences of current and previous year SPEI (Standardized Precipitation Evapotranspiration Index). Current year SPEI is positively correlated with defoliation, while previous year SPEI has a negative effect. Although significant, this divergent effect of drought should be regarded with caution, since the estimated effects are weak (Table 1). Both the current and previous year temperature and precipitation did not have significant effects on defoliation. Plots located on higher altitudes had significantly higher mean defoliation. The annual increase in defoliation estimated in the LMM models is very similar to the positive defoliation trend assessed by the Mann–Kendall test (Figure 1). The marginal coefficient of determination (R2) was lower than the conditional R2 in both models, which suggests that the model random effects i.e., plot location, accounts for a high proportion of the explained DEF variance (Table 2). Table 1. Estimated model parameters, standard errors, t-values, and p-values for mean defoliation LMM models with current year (DEF-I) and previous year climate variables (DEF-II). T—mean annual temperature, P—annual sum of precipitation, SPEI3—mean annual Standardized Figure 1. Overall trend of crown vitality parameters from 1996–2017. (a) Defoliation trend of European beech (Tau = 0.78, Sen’s slope = 0.39, p < 0.001, orange line) and annual overall mean defoliation (DEFi, black line). (b) Defoliation change trend of European beech (Tau = −0.05, Sen’s slope = −0.01, p = 0.73, orange line) and annual overall change in mean defoliation (∆DEFi, black line). Points represent annual plot mean values. Plot-scale trend analysis of DEF showed that 53% of plots have a significant and increasing trend of defoliation at an annual rate ranging from 0.25% to 1.29%. On the re- maining 47% of plots, we did not observe any significant trend. The results of the q-statistics and Morans’ I spatial autocorrelation coefficient did not indicate spatial stratification of plot-scale DEF trends (data not shown). We did not detect any significant ∆DEF trend at the plot-scale. Therefore, the subsequent spatial stratification tests were not conducted. 2.2. Influence of Environmental Conditions on Tree Vitality During the selection procedure for DEF and ∆DEF models, additional variables were considered: stand age, site factors (altitude and orientation), and soil properties (soil pH, total nitrogen content in the soil, content of available phosphorus and available potassium in the soil) were tested. However, only altitude showed a significant impact. Although the influence of stand age was not significant, it was retained in the selection process due to improved model performance and the reported influence of stand age on defoliation in other studies [38,39]. The two linear mixed effects models (LMM) used to assess the impact of climatic variables on DEF revealed different influences of current and previous year SPEI (Standard- ized Precipitation Evapotranspiration Index). Current year SPEI is positively correlated with defoliation, while previous year SPEI has a negative effect. Although significant, this divergent effect of drought should be regarded with caution, since the estimated effects are weak (Table 1). Both the current and previous year temperature and precipitation did not have significant effects on defoliation. Plots located on higher altitudes had significantly higher mean defoliation. The annual increase in defoliation estimated in the LMM models is very similar to the positive defoliation trend assessed by the Mann–Kendall test (Figure 1). The marginal coefficient of determination (R2) was lower than the conditional R2 in both models, which suggests that the model random effects i.e., plot location, accounts for a high proportion of the explained DEF variance (Table 2). Plants 2022, 11, 730 4 of 13 Table 1. Estimated model parameters, standard errors, t-values, and p-values for mean defoliation LMM models with current year (DEF-I) and previous year climate variables (DEF-II). T—mean annual temperature, P—annual sum of precipitation, SPEI3—mean annual Standardized Precipitation Evapotranspiration Index calculated on a three-month time scale, lag—denotes previous year values. Estimate Std. Error t Value p Value DEF-I Intercept −7.924 1.591 −4.980 <0.001 Year 0.004 8.03 × 10−4 4.934 <0.001 Stand age −4.70 × 10−4 3.16 × 10−4 −1.486 0.138 Altitude 1.58 × 10−4 3.86 × 10−5 4.095 <0.01 T 0.007 0.004 1.954 0.051 P −1.07 × 10−5 1.19 × 10−5 −0.902 0.367 SPEI3 0.016 0.006 2.834 <0.01 DEF-II Intercept −7.542 1.571 −4.799 <0.001 Year 0.004 7.92 × 10−4 4.767 <0.001 Stand age −3.68 × 10−4 3.09 × 10−4 −1.191 0.234 Altitude 1.18 × 10−4 3.82 × 10−5 3.099 <0.01 T_lag 0.005 0.003 1.553 0.121 P_lag 9.67 × 10−6 1.19 × 10−5 0.813 0.417 SPEI3_lag −0.014 0.006 −2.501 <0.05 Table 2. Defoliation (DEF) and defoliation change (∆DEF) model performance indices: Akaike Information Criterion (AIC); Conditional coefficient of determination (Conditional R2); Marginal coefficient of determination (Marginal R2); Intraclass Correlation Coefficient (ICC); Root Mean Square Error (RMSE). AIC Conditional R2 Marginal R2 ICC RMSE DEF-I −1906 0.89 0.60 0.71 0.061 DEF-II −1915 0.88 0.62 0.68 0.062 ∆DEF-I 1824 0.09 0.09 1.20 × 10−8 0.996 ∆DEF-II 1806 0.32 0.32 7.07 × 10−9 0.984 Defoliation change LMM indicates a very strong and significant negative influence of previous year SPEI (Table 3). Positive defoliation changes were observed in years preceded by low SPEI, while negative defoliation changes were associated with high SPEI the previous year. This inverse relationship where the increase of defoliation is preceded by drought was most notable in 2001 (Figure 2). Other climate variables, including current year SPEI, did not have a significant effect and neither did stand age or elevation. Equal values of conditional and marginal R2 suggest a uniform effect of previous year drought on European beech ∆DEF (Table 2). As expected, the models did not reveal a significant trend in the change of defoliation. Plants 2022, 11, 730 5 of 13 Table 3. Estimated regression parameters, standard errors, t-values, and p-values for mean defoliation change LMM models with current year (∆DEF-I) and previous year climate variables (∆DEF-II). T—mean annual temperature, P—annual sum of precipitation, SPEI3—mean annual Standardized Precipitation Evapotranspiration Index calculated on a three-month time scale, lag—denotes previous year values. Estimate Std. Error t Value p Value ∆DEF-I Intercept 25.114 15.641 1.606 0.109 Year −0.013 0.008 −1.595 0.111 Stand age −1.91 × 10−4 0.002 −0.104 0.917 Altitude −1.50 × 10−4 3.87 × 10−4 −0.388 0.701 T 0.012 0.057 0.207 0.836 P 1.47 × 10−4 1.35 × 10−4 1.094 0.274 SPEI3 −0.006 0.091 −0.066 0.947 ∆DEF-II Intercept 15.966 15.833 1.008 0.314 Year −0.008 0.008 −0.964 0.336 Stand age −2.50 × 10−5 0.002 −0.014 0.989 Altitude −0.001 3.78 × 10−4 −1.388 0.176 T_lag −0.047 0.055 −0.857 0.392 P_lag 1.98 × 10−4 1.33 × 10−4 1.487 0.138 SPEI3_lag −0.388 0.089 −4.338 <0.01 Plants 2022, 11, x FOR PEER REVIEW 5 of 13 Figure 2. Inverse relationship between previous year SPEI (Standardized Precipitation Evapotranspiration Index) and defoliation change (ΔDEFi) in Croatia from 1996–2017. (a) Previous year mean SPEI (orange line) and standard error (orange area), calculated on a time scale of three months. (b) Annual overall defoliation change (ΔDEFi) mean (green line) and standard error (green area). 3. Discussion 3.1. Temporal Trends in Crown Vitality Defoliation is widely accepted as a proxy indicator of tree vitality and forest health, able to provide useful information on its status and trends. Long-term defoliation data series are an important asset to explore the changes in forest ecosystem health across Europe over the past 30 years [8]. This study revealed a statistically significant trend of increasing defoliation of European beech in Croatia over time. The distinctive impulses of increasing mean defoliation in 2001 and from 2012 to 2014 (Figure 1) are preceded by the dry 2000 and the extreme drought period recorded in 2011/2012 (Figure S1) [40]. The severe drought recorded in 2003 did not result in an increase in mean defoliation in the following years, which is consistent with our previous findings [37]. Sudden increases in mean defoliation preceded by years or periods of pronounced precipitation deficit were also recorded in studies conducted in France [33] and the Iberian Peninsula [41]. A significant, but weaker trend of increasing beech defoliation was also observed at a European level [42]. A study comparing the general defoliation trends between geographical regions found that southern Europe, including Croatia, has a more pronounced trend of increasing defoliation compared to central and northern Europe [41]. P 1.47 × 10−4 1.35 × 10−4 1.094 0.274 SPEI3 −0.006 0.091 −0.066 0.947 ΔDEF-II Intercept 15.966 15.833 1.008 0.314 Year −0.008 0.008 −0.964 0.336 Stand age −2.50 × 10−5 0.002 −0.014 0.989 Altitude −0.001 3.78 × 10−4 −1.388 0.176 T_lag −0.047 0.055 −0.857 0.392 P_lag 1.98 × 10−4 1.33 × 10−4 1.487 0.138 SPEI3_lag −0.388 0.089 −4.338 <0.01 Figure 2. Inverse relationship between previous year SPEI (Standardized Precipitation Evapotran- spiration Index) and defoliation change (∆DEFi) in Croatia from 1996–2017. (a) Previous year mean SPEI (orange line) and standard error (orange area), calculated on a time scale of three months. (b) Annual overall defoliation change (∆DEFi) mean (green line) and standard error (green area). 3. Discussion 3.1. Temporal Trends in Crown Vitality Defoliation is widely accepted as a proxy indicator of tree vitality and forest health, able to provide useful information on its status and trends. Long-term defoliation data series are an important asset to explore the changes in forest ecosystem health across Europe over the past 30 years [8]. This study revealed a statistically significant trend of increasing defoliation of European beech in Croatia over time. The distinctive impulses of increasing mean defoliation in 2001 and from 2012 to 2014 (Figure 1) are preceded by the dry 2000 and the extreme drought period recorded in 2011/2012 (Figure S1) [40]. The severe drought recorded in 2003 did not result in an increase in mean defoliation in the Plants 2022, 11, 730 6 of 13 following years, which is consistent with our previous findings [37]. Sudden increases in mean defoliation preceded by years or periods of pronounced precipitation deficit were also recorded in studies conducted in France [33] and the Iberian Peninsula [41]. A significant, but weaker trend of increasing beech defoliation was also observed at a European level [42]. A study comparing the general defoliation trends between geographical regions found that southern Europe, including Croatia, has a more pronounced trend of increasing defoliation compared to central and northern Europe [41]. Plot-scale analysis of defoliation revealed a statistically significant increase on 53% of plots in the 1996 to 2017 period. A study applying a similar plot-scale approach in France found that as many as 70% of beech plots showed an increasing defoliation trend from 1996 to 2009 [33]. Our results also suggest that there is no spatial grouping of the defoliation trend, which is in line with other studies [9,31]. Obviously, the trend of defoliation is not influenced by geographical position, but rather by specific environmental conditions present on a particular plot. Defoliation change (∆DEF), defined as the difference between the defoliation as- sessed in the current and previous year has so far been used in only a few studies [31,43]. Unfortunately, detailed results from these studies were not provided, and therefore a straightforward comparison with our results was not possible. While differences in as- sessment of absolute defoliation values can be expected due to national adjustments of the methods, the differences in assessing the relative change of defoliation from year to year should be negligible, and could potentially reduce the influence of possible subjec- tivity of the assessment [43]. Additionally, the absence of serial correlation of defoliation change, i.e., its stationarity, enables easier development of impact models compared to using defoliation data. 3.2. Influence of Environmental Conditions on Tree Vitality A key task at the European level is to study the impact of climate change on crown defoliation and, consequently, on forest health [44], taking into consideration a wide range of natural and anthropogenic environmental factors [45]. The established difference in marginal and conditional R2 in crown defoliation models suggests that site factors explain 25 to 30% variance of European beech defoliation over time. A high influence of specific site attributes on defoliation was also found in other studies, e.g., lower defoliation values were observed at higher nitrogen supply levels and higher pH levels [46,47]. However, soil properties did not show a significant impact during the model selection process in our study, which is in line with several studies that did not confirm the importance of soil properties for European beech defoliation [10,43]. It is possible that the applied approach to soil sampling and analysis does not provide a sufficient level of detail to detect the significance of soil variables, given that soil properties have nevertheless been shown as significant factors in some studies of European beech [48] and Norway spruce defoliation [49]. On the other hand, effects from environmental factors that change on a long-term scale, like soil properties, will not likely be detected through annual variation of defoliation [46]. Stand age was identified as a significant predictor in different approaches to modelling defoliation [39,40,45,50]. However, the established relationship between defoliation and stand age may, in many cases, represent an interaction between various stress factors and age [4]. The absence of a significant impact of stand age on European beech defoliation in Croatia (Tables 1 and 3) can be explained by a relatively small number of plots in our sample, where stands are older than 80 years (Figure S2), while the abovementioned studies were not limited by irregular age distribution. Numerous studies have observed a significant effect of drought on increasing defo- liation [32,38,39,42,45], and increased leaf loss following spring and summer heat waves was recorded both in European [50,51] and North American forests [52]. In contrast to the clear influence of spring and summer temperatures on European beech defoliation found in Spain [34], the results of this study indicate a very weak influence of all examined tempera- Plants 2022, 11, 730 7 of 13 ture variables, which is consistent with results of a beech study conducted in Germany [46]. Furthermore, we did not find a significant influence of precipitation on defoliation. This is contrary to the results of a French study, where precipitation and precipitation deficit correlated with defoliation [33]. While we could not detect any effect of air temperature and precipitation when considered independently, these climate variables showed to have a clear impact on defoliation change when combined in the SPEI drought index. The basic mechanism for regulating water loss in dry conditions is stomatal closure in plants [53]. Under conditions of increased water deficit, plants also respond by increasing water use efficiency [54], reduced growth [55], and conservative mechanisms such as limit- ing their photosynthetic activity [56]. Due to drought, plants can adapt their morphological structure by increasing the carbon allocation to the root system [57], reducing their leaf size [58], decreasing leaf area index [59], and ultimately shedding leaves [10]. During long lasting drought events, stomatal closure can significantly reduce carbon fixation by trees as well as their carbon reserves, which weakens trees and makes them more vulnerable to biotic and abiotic stresses. In extreme cases, this can lead to mortality by carbon starva- tion [60]. Severe drought during the year of bud formation, in our study indicated by low previous year SPEI and its impact on ∆DEF, decreases the number of new leaves formed in the bud thus influencing the number of leaves, leaf surface area, and twig extension in the following year [51]. In Fagus species, all leaves are pre-formed in winter buds [61,62] during late summer and early autumn [59]. Hydraulic failure may also occur during severe droughts leading to twig and leaf abscission, which can be seen as a drastic adaptation strategy to reduce evapotranspiration [63]. This effect is visible from the values of defolia- tion rising in the period from 2011 until 2014 (Figure 1), which coincides with low SPEI for the years 2011 and 2012 (Figure 2). Equal values of conditional and marginal R2 and the low value of ICC (intraclass correlation coefficient) suggests a uniform influence of previous year drought on the European beech ∆DEF throughout Croatia. On the other hand, the increase in the marginal R2 of defoliation models after the inclusion of climatic variables was slight, and the influence of all observed climatic variables on DEF was weak. This indicates that defoliation is influenced by site-specific environmental or stand factors that have not been identified in this study. Mean beech defoliation shows fluctuations that coincide with the occurrence of common to abundant fructification [48,64]. However, we were not able to include fruiting as a factor due to the lack of data. The lack of a clear ∆DEF trend, as well as the pronounced impact of drought in the previous year, may indicate that this response variable reflects intense stress, while the positive DEF trend represents the accumulated impact of location- specific stressful environmental conditions on tree vitality. Since defoliation change shows the current or recent status of tree vitality, while defoliation is an integrated indicator resulting from cumulated biotic and abiotic pressures on tree vitality over many years, defoliation change could be a more appropriate indicator for analysing the effect of recent climate conditions on tree vitality. Increasing temperatures may lead to drought thus affecting forest vitality in the region. Forest monitoring activities in southeast Europe should be intensified to determine the unknown site-specific environmental and/or stand factors that may explain a part of the variance in the present data. This could help develop adequate and locally applicable mitigation strategies to secure the future of beech forests in the region. 4. Materials and Methods 4.1. Study Area and Plot Selection The ICP Forests Level I monitoring plots in Croatia are established on intersections of a 16 × 16 km grid that contain forest cover. These plots do not have a fixed area; rather, 24 trees are chosen for defoliation assessments using a cross-cluster system with six trees in each cluster [65]. Only plots with a minimum of five European beech trees were selected to ensure that European beech was significantly represented in the mixture of tree species. To ensure defoliation data consistency over the investigated period from 1996 to 2017, Plants 2022, 11, 730 8 of 13 defoliation assessment on selected plots had to have been carried out for at least 80% of the investigated period. This resulted in the selection of 28 research plots (Figure 3). In addition to defoliation, the ICP Forests database contains information on several site factors (Table S1). Plants 2022, 11, x FOR PEER REVIEW 8 of 13 Figure 3. Location of research plots in the context of the distribution of European beech in Europe [66]. Selected Level I ICP Forests monitoring plots (orange dots); remaining Level I plots in Croatian (black dots); distribution of European beech forests in Croatia [67] (green polygon). 4.2. Defoliation Assessment and Crown Vitality Indicators Defoliation of European beech trees on the selected plots was assessed annually between mid-July and mid-August from 1996 to 2017, in 5% classes from 0 to 100%, according to the ICP Forests Manual [12]. Assessments of tree crowns was performed in comparison with the absolute reference tree. For this study, two crown vitality indicators were calculated and used as response variables: (i) the mean current year crown defoliation DEFi on plot i, and (ii) the change in the mean current year crown defoliation on plot icompared to the previous year assessment ΔDEFi. Mean values of crown vitality parameters at the plot level were used, since the values of all predictor variables could only be obtained at the plot level. Additionally, comparison of defoliation variability within and between plots showed that it was lower within plots than between plots. 4.3. Soil Sampling and Analysis Soil sampling was performed during the summer of 2019 on five points located within each of the research plots. One point was located within each of the four groups of trees that are assessed for defoliation, and an additional fifth point was located in the centre of each research plot. Soil samples were taken with a pedological drill from a depth of 0–10 cm, 10–20 cm, 20–40 cm, and 40–80 cm. Collected samples were pooled according to the sampling depth. Soil chemical parameters were analysed according to standard protocols and methods (Table S2). 4.4. Climate Data Climate monitoring stations are generally situated at considerable distances from the research plots. Therefore, the data they provide are not always representative of the research locations. To overcome this, we used gridded data produced by regression kriging (RK), which is a hybrid method of interpolation carried out in four steps [68]. The method was validated with leave-one-out cross-validation, while the root mean square Figure 3. Location of research plots in the context of the distribution of European beech in Europe [66]. Selected Level I ICP Forests monitoring plots (orange dots); remaining Level I plots in Croatian (black dots); distribution of European beech forests in Croatia [67] (green polygon). 4.2. Defoliation Assessment and Crown Vitality Indicators Defoliation of European beech trees on the selected plots was assessed annually between mid-July and mid-August from 1996 to 2017, in 5% classes from 0 to 100%, according to the ICP Forests Manual [12]. Assessments of tree crowns was performed in comparison with the absolute reference tree. For this study, two crown vitality indicators were calculated and used as response variables: (i) the mean current year crown defoliation DEFi on plot i, and (ii) the change in the mean current year crown defoliation on plot i compared to the previous year assessment ∆DEFi. Mean values of crown vitality parameters at the plot level were used, since the values of all predictor variables could only be obtained at the plot level. Additionally, comparison of defoliation variability within and between plots showed that it was lower within plots than between plots. 4.3. Soil Sampling and Analysis Soil sampling was performed during the summer of 2019 on five points located within each of the research plots. One point was located within each of the four groups of trees that are assessed for defoliation, and an additional fifth point was located in the centre of each research plot. Soil samples were taken with a pedological drill from a depth of 0–10 cm, 10–20 cm, 20–40 cm, and 40–80 cm. Collected samples were pooled according to the sampling depth. Soil chemical parameters were analysed according to standard protocols and methods (Table S2). Plants 2022, 11, 730 9 of 13 4.4. Climate Data Climate monitoring stations are generally situated at considerable distances from the research plots. Therefore, the data they provide are not always representative of the research locations. To overcome this, we used gridded data produced by regression kriging (RK), which is a hybrid method of interpolation carried out in four steps [68]. The method was validated with leave-one-out cross-validation, while the root mean square error (RMSE) was calculated between observed and interpolated values. Mean RMSEs are for mean monthly temperature from 0.5 ◦C to 0.9 ◦C, for minimum temperature from 1.1 ◦C to 1.5 ◦C, for maximum temperature from 0.7 ◦C to 1.1 ◦C, and for precipitation from 18 to 30 mm, averaged by months. Mean monthly temperature (T), minimum (Tmin) and maximum (Tmax) monthly temperature, and monthly sum of precipitation (P) from the gridded dataset on 1 km spatial resolution for Croatia [69] were used to calculate yearly values, as well as the Palmer Drought Severity Index (scPDSI) [70] and Standardized Precipitation Evapotranspiration Index (SPEI) [71]. Lower values of scPDSI and SPEI indicate a stronger drought intensity while higher values indicate a higher degree of humidity. SPEI was calculated on a time scale of 3, 6, and 12 months. 4.5. Data Analysis The trend of defoliation and defoliation change was estimated according to Sen’s slope [72], while the significance of a trend was tested by the Mann–Kendall test [73,74] with a significance level of p ≤ 0.05. Both methods are suitable for data with asymmetric distribution and in this case are significantly more accurate compared to the simple linear regression model [75]. Spatial stratification of defoliation plot-wise trends were examined by calculating the degree of spatial stratified heterogeneity using the q-statistics method [76] and the spatial autocorrelation coefficient, Moran’s I [77]. To model the impact of site factors and climate conditions on crown vitality indicators, we used linear mixed effects models (LMM) [78]. Prior to adding climate predictors, a default model was fitted: DEFit = Yeart + StandAgeit + Elevationi SamplePloti ~ N(0, σ2) (1) where DEFit is the mean crown defoliation of European beech trees for sample plot i = 1, . . . , n and for year t = 1, . . . , 22, averaged over all trees at sample plot i. SamplePloti is the random intercept, which is assumed to be normally distributed with mean 0 and variance σ2. To account for different number of trees on each plot, weights 1/αit were introduced, where αit is the number of trees assessed at plot i and year t. Since defoliation represents an estimated percentage and due to the LMM requirements, mean defoliation values were divided by 100 before model fitting. First order autocorrelative term was introduced to account for temporal autocorrelation in the model. Seidling [46] states that the serial correlation of European beech that appears over a five-year period is not as distinct as in other species studied, which is contrary to our data. Ignoring serial correlation in model fitting leads to overestimated random effects and to the inflation of the empirical Type I error rates [79], therefore it is crucial to account for this during model fitting of defoliation data. Due to high kurtosis, The Lambert W × F function [80] was applied to transform ∆DEF data to a normal distribution. Afterward, the same approach to the base model build up was applied as with DEF data, except that the first order autocorrelative term was left out since the data did not display serial correlation. The final model selection process was based on diagnostic diagrams and a procedure defined by Johnson and Omland [81]. Spearman’s correlation coefficient (rho) was calcu- lated between the crown vitality indicators and each quantitative environmental variable in order to obtain an overview of possible impacts. Of the potential independent variables, those that explain most of the variation of crown vitality indicators were selected with Plants 2022, 11, 730 10 of 13 a recursive feature elimination approach (RFE) implemented within the random forest algorithm [82]. Selected variables which were linearly correlated with other variables and had a variance inflation factors VIF > 5, a commonly used threshold in detecting multicollinearity [83], were identified. The identified collinear variables with the lower value according to the Akaike Information Criteria (AIC) [84] were retained for further model development. This subset of uncorrelated environmental variables was used as predictor variables for developing the final models. All analyses were conducted in an R programming environment [85]. Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/plants11060730/s1, References [86–91] are cited in the Supplementary Materials. Figure S1. Lin- ear regression fits (black lines) and kernel smoothing functions (colour lines) for one-year smoothed mean monthly temperature, monthly precipitation sum and mean monthly SPEI data (grey lines) on the research plots during the vegetation period., Figure S2. Distribution of 28 research plots by stand age, altitude, and orientation, Table S1. Additional site factor variables and data sources, Table S2. Descriptive statistics of soil chemical properties from 28 research plots and the applied methods of analysis Author Contributions: Conceptualization, M.O., I.S. and N.P.; Formal analysis, M.O.; Investigation, M.O.; Methodology, M.O.; Project administration, M.O., I.S. and N.P.; Resources, N.P. and M.P.T.; Supervision, N.P. and D.U.; Visualization, M.O.; Writing—original draft, M.O., I.S. and N.P.; Writing— review and editing, M.O., I.S., N.P., M.M., M.P.T., M.J., P.R., V.T., L.L. and D.U. All authors have read and agreed to the published version of the manuscript. Funding: This work has been fully supported by Croatian Science Foundation under the project VitaClim (IP-2018-01-5222). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: The data presented in this study are available on request from the corresponding author. The data are not publicly available due to legal reasons. Acknowledgments: We highly appreciate the work of researchers and technicians of the Croatian Forest Research Institute that participated in field research and laboratory analysis for this study. Conflicts of Interest: The authors declare no conflict of interest. References 1. Centritto, M.; Tognetti, R.; Leitgeb, E.; Střelcová, K.; Cohen, S. Above Ground Processes: Anticipating Climate Change Influences. In Forest Management and the Water Cycle: An Ecosystem-Based Approach; Bredemeier, M., Cohen, S., Godbold, D.L., Lode, E., Pichler, V., Schleppi, P., Eds.; Springer: Dordrecht, The Netherlands, 2011; pp. 31–64. 2. Trumbore, S.; Brando, P.; Hartmann, H. Forest health and global change. Science 2015, 349, 814–818. [CrossRef] [PubMed] 3. Stocker, T.F.; Qin, D.; Plattner, G.-K.; Alexander, L.V.; Allen, S.K.; Bindoff, N.L.; Bréon, F.-M.; Church, J.A.; Cubasch, U.; Emori, S.; et al. Technical Summary. In Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Changep; Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013. 4. De Vries, W.; Dobbertin, M.H.; Solberg, S.; van Dobben, H.F.; Schaub, M. Impacts of acid deposition, ozone exposure and weather conditions on forest ecosystems in Europe: An overview. Plant Soil 2014, 380, 1–45. [CrossRef] 5. Ferretti, M.; Waldner, P.; Verstraeten, A.; Schmitz, A.; Michel, A.; Žlindra, D.; Marchetto, A.; Hansen, K.; Pitar, D.; Gottardini, E.; et al. Criterion 2: Maintenance of Forest Ecosystem Health and Vitality. In FOREST EUROPE, 2020: State of Europe’s Forests 2020; Ministerial Conference on the Protection of Forests in Europe—Liaison Unit Bratislava: Zvolen, Slovak Republic, 2020. 6. Brang, P. Sanasilva-Bericht 1997: Zustand und Gefährdung des Schweizer Waldes-eine Zwischenbilanz nach 15 Jahren Waldschaden- forschung; Bundesamt für Umwelt Wald und Landschaft; Eidgenössische Forschungsanstalt: Birmensdorf, Switzerland, 1998. 7. Cherubini, P.; Battipaglia, G.; Innes, J.L. Tree Vitality and Forest Health: Can Tree-Ring Stable Isotopes Be Used as Indicators? Curr. For. Rep. 2021, 7, 69–80. [CrossRef] 8. Gottardini, E.; Cristofolini, F.; Cristofori, A.; Pollastrini, M.; Camin, F.; Ferretti, M. A multi-proxy approach reveals common and species-specific features associated with tree defoliation in broadleaved species. For. Ecol. Manag. 2020, 467, 118151. [CrossRef] 9. De Vries, W.; Klap, J.M.; Erisman, J.W. Effects of environmental stress on forest crown condition in Europe. Part I: Hypotheses and approach to the study. Water Air Soil Pollut. 2000, 119, 317–333. [CrossRef] https://www.mdpi.com/article/10.3390/plants11060730/s1 https://www.mdpi.com/article/10.3390/plants11060730/s1 http://doi.org/10.1126/science.aac6759 http://www.ncbi.nlm.nih.gov/pubmed/26293952 http://doi.org/10.1007/s11104-014-2056-2 http://doi.org/10.1007/s40725-021-00137-8 http://doi.org/10.1016/j.foreco.2020.118151 http://doi.org/10.1023/A:1005157509454 Plants 2022, 11, 730 11 of 13 10. De Marco, A.; Proietti, C.; Cionni, I.; Fischer, R.; Screpanti, A.; Vitale, M. Future impacts of nitrogen deposition and climate change scenarios on forest crown defoliation. Environ. Pollut. 2014, 194, 171–180. [CrossRef] 11. Dobbertin, M. Tree growth as indicator of tree vitality and of tree reaction to environmental stress: A review. Eur. J. For. Res. 2005, 124, 319–333. [CrossRef] 12. Eichhorn, J.; Roskams, P.; Potočić, N.; Timmermann, V.; Ferretti, M.; Mues, V.; Szepesi, A.; Durrant, D.; Seletković, I.; Schroeck, H.-W.; et al. Part IV: Visual Assessment of Crown Condition and Damaging Agents. Version 2020-3. In Manual on Methods and Criteria for Harmonized Sampling, Assessment, Monitoring and Analysis of the Effects of Air Pollution on Forests; Centre, U.I.F.P.C., Ed.; Thünen Institute of Forest Ecosystems: Eberswalde, Germany, 2020; p. 49 + Annex. 13. Landmann, G. Forest decline and air pollution effects in the French mountains: A synthesis. In Forest Decline and Atmospheric Deposition Effects in the French Mountains; Springer: Berlin/Heidelberg, Germany, 1995; pp. 407–452. 14. Ghosh, S.; Innes, J.L.; Hoffmann, C. Observer Variation as a Source of Error in Assessments of Crown Condition Through Time. For. Sci. 1995, 41, 235–254. [CrossRef] 15. Innes, J.L.; Landmann, G.; Mettendorf, B. Consistency of observations of forest tree defoliation in three European countries. Environ. Monit. Assess. 1993, 25, 29–40. [CrossRef] 16. Johnson, J.; Jacob, M. Monitoring the effects of air pollution on forest condition in Europe: Is crown defoliation an adequate indicator? Iforest—Biogeosci. For. 2010, 3, 86–88. [CrossRef] 17. Eickenscheidt, N.; Wellbrock, N. Consistency of defoliation data of the national training courses for the forest condition survey in Germany from 1992 to 2012. Environ. Monit. Assess. 2014, 186, 257–275. [CrossRef] [PubMed] 18. Ferretti, M.; Bussotti, F.; Cenni, E.; Cozzi, A. Implementation of Quality Assurance Procedures in the Italian Programs of Forest Condition Monitoring. Water Air Soil Pollut. 1999, 116, 371–376. [CrossRef] 19. Wulff, S. The Accuracy of Forest Damage Assessments—Experiences from Sweden. Environ. Monit. Assess. 2002, 74, 295–309. [CrossRef] 20. Ferretti, M.; König, N.; Rautio, P.; Sase, H. Quality assurance (QA) in international forest monitoring programmes: Activity, problems and perspectives from East Asia and Europe. Ann. For. Sci. 2009, 66, 403. [CrossRef] 21. Leuschner, C.; Ellenberg, H. Ecology of Central European Forests: Vegetation Ecology of Central Europe; Springer International Publishing: Cham, Switzerland, 2017; Volume 1, pp. 1–971. 22. Duncker, P.S.; Raulund-Rasmussen, K.; Gundersen, P.; Katzensteiner, K.; De Jong, J.; Ravn, H.P.; Smith, M.; Eckmullner, O.; Spiecker, H. How Forest Management affects Ecosystem Services, including Timber Production and Economic Return: Synergies and Trade-Offs. Ecol. Soc. 2012, 17, 17. [CrossRef] 23. Zimmermann, J.; Hauck, M.; Dulamsuren, C.; Leuschner, C. Climate Warming-Related Growth Decline Affects Fagus sylvatica, But Not Other Broad-Leaved Tree Species in Central European Mixed Forests. Ecosystems 2015, 18, 560–572. [CrossRef] 24. Simon, J.; Dannenmann, M.; Pena, R.; Gessler, A.; Rennenberg, H. Nitrogen nutrition of beech forests in a changing climate: Importance of plant-soil-microbe water, carbon, and nitrogen interactions. Plant Soil 2017, 418, 89–114. [CrossRef] 25. Beniston, M.; Stephenson, D.B.; Christensen, O.B.; Ferro, C.A.T.; Frei, C.; Goyette, S.; Halsnaes, K.; Holt, T.; Jylhä, K.; Koffi, B.; et al. Future extreme events in European climate: An exploration of regional climate model projections. Clim. Chang. 2007, 81, 71–95. [CrossRef] 26. Giorgi, F. Climate change hot-spots. Geophys. Res. Lett. 2006, 33, 113–128. [CrossRef] 27. Gajić-Čapka, M.; Cindrić, K.; Pasarić, Z. Trends in precipitation indices in Croatia, 1961–2010. Theor. Appl. Climatol. 2015, 121, 167–177. [CrossRef] 28. Zaninović, K.; Gajić-Čapka, M. Changes in Components of the Water Balance in the Croatian Lowlands. Theor. Appl. Climatol. 2000, 65, 111–117. [CrossRef] 29. Spinoni, J.; Antofie, T.; Barbosa, P.; Bihari, Z.; Lakatos, M.; Szalai, S.; Szentimrey, T.; Vogt, J. An overview of drought events in the Carpathian Region in 1961–2010. Adv. Sci. Res. 2013, 10, 21–32. [CrossRef] 30. Lorenz, M.; Fischer, R. Pan-European Forest Monitoring: An Overview. In Developments in Environmental Science; Volume Forest Monitoring—Methods for Terrestrial Investigations in Europe with an Overview of North America and Asia; Marco Ferretti, R.F., Ed.; Elsevier: Amsterdam, The Netherlands, 2013; p. 14. 31. Zierl, B. A simulation study to analyse the relations between crown condition and drought in Switzerland. For. Ecol. Manag. 2004, 188, 25–38. [CrossRef] 32. Seidling, W.; Ziche, D.; Beck, W. Climate responses and interrelations of stem increment and crown transparency in Norway spruce, Scots pine, and Common beech. For. Ecol. Manag. 2012, 284, 196–204. [CrossRef] 33. Ferretti, M.; Nicolas, M.; Bacaro, G.; Brunialti, G.; Calderisi, M.; Croisé, L.; Frati, L.; Lanier, M.; Maccherini, S.; Santi, E.; et al. Plot-scale modelling to detect size, extent, and correlates of changes in tree defoliation in French high forests. For. Ecol. Manag. 2014, 311, 56–69. [CrossRef] 34. De la Cruz, A.C.; Gil, P.M.; Fernández-Cancio, Á.; Minaya, M.; Navarro-Cerrillo, R.M.; Sánchez-Salguero, R.; Grau, J.M. Defoliation triggered by climate induced effects in Spanish ICP Forests monitoring plots. For. Ecol. Manag. 2014, 331, 245–255. [CrossRef] 35. Seletković, I.; Potočić, N.; Ugarković, D.; Jazbec, A.; Pernar, R.; Seletković, A.; Benko, M. Climate and relief properties influence crown condition of common beech (Fagus sylvatica L.) on the Medvednica massif. Period. Biol. 2009, 111, 435–441. 36. Potočić, N.; Seletković, I.; Ugarković, D.; Jazbec, A.; Mikac, S. The influence of climate properties on crown condition of Common beech (Fagus sylvatica L.) and Silver fir (Abies alba Mill.) on Velebit. Period. Biol. 2008, 110, 145–150. http://doi.org/10.1016/j.envpol.2014.07.027 http://doi.org/10.1007/s10342-005-0085-3 http://doi.org/10.1093/forestscience/41.2.235 http://doi.org/10.1007/BF00549790 http://doi.org/10.3832/ifor0538-003 http://doi.org/10.1007/s10661-013-3372-3 http://www.ncbi.nlm.nih.gov/pubmed/23955498 http://doi.org/10.1023/A:1005240000294 http://doi.org/10.1023/A:1014287228807 http://doi.org/10.1051/forest/2009025 http://doi.org/10.5751/ES-05066-170450 http://doi.org/10.1007/s10021-015-9849-x http://doi.org/10.1007/s11104-017-3293-y http://doi.org/10.1007/s10584-006-9226-z http://doi.org/10.1029/2006GL025734 http://doi.org/10.1007/s00704-014-1217-9 http://doi.org/10.1007/s007040050008 http://doi.org/10.5194/asr-10-21-2013 http://doi.org/10.1016/j.foreco.2003.07.019 http://doi.org/10.1016/j.foreco.2012.07.015 http://doi.org/10.1016/j.foreco.2013.05.009 http://doi.org/10.1016/j.foreco.2014.08.010 Plants 2022, 11, 730 12 of 13 37. Ognjenović, M.; Levanič, T.; Potočić, N.; Ugarković, D.; Indir, K.; Seletković, I. Interrelations of various tree vitality indicators and their reaction to climatic conditions on a european beech (Fagus sylvatica L.) plot. Šumar. List 2020, 144, 351–365. [CrossRef] 38. Sousa-Silva, R.; Verheyen, K.; Ponette, Q.; Bay, E.; Sioen, G.; Titeux, H.; Van de Peer, T.; Van Meerbeek, K.; Muys, B. Tree diversity mitigates defoliation after a drought-induced tipping point. Glob. Chang. Biol. 2018, 24, 4304–4315. [CrossRef] 39. Augustin, N.H.; Musio, M.; von Wilpert, K.; Kublin, E.; Wood, S.N.; Schumacher, M. Modeling Spatiotemporal Forest Health Monitoring Data. J. Am. Stat. Assoc. 2009, 104, 899–911. [CrossRef] 40. Cindrić, K.; Telišman Prtenjak, M.; Herceg-Bulić, I.; Mihajlović, D.; Pasarić, Z. Analysis of the extraordinary 2011/2012 drought in Croatia. Theor. Appl. Climatol. 2016, 123, 503–522. [CrossRef] 41. Carnicer, J.; Coll, M.; Ninyerola, M.; Pons, X.; Sánchez, G.; Peñuelas, J. Widespread crown condition decline, food web disruption, and amplified tree mortality with increased climate change-type drought. Proc. Natl. Acad. Sci. USA 2011, 108, 1474–1478. [CrossRef] [PubMed] 42. Timmermann, V.; Potočić, N.; Ognjenović, M.; Kirchner, T. Tree Crown Condition in 2019; Thünen Institute: Eberswalde, Germany, 2020. 43. Klap, J.M.; Oude Voshaar, J.H.; De Vries, W.; Erisman, J.W. Effects of Environmental Stress on Forest Crown Condition in Europe. Part IV: Statistical Analysis of Relationships. Water Air Soil Pollut. 2000, 119, 387–420. [CrossRef] 44. Vitale, M.; Proietti, C.; Cionni, I.; Fischer, R.; De Marco, A. Random Forests Analysis: A Useful Tool for Defining the Relative Importance of Environmental Conditions on Crown Defoliation. Water Air Soil Pollut. 2014, 225, 1992. [CrossRef] 45. Seidling, W. Signals of summer drought in crown condition data from the German Level I network. Eur. J. For. Res. 2007, 126, 529–544. [CrossRef] 46. Seidling, W. Crown condition within integrated evaluations of Level II monitoring data at the German level. Eur. J. For. Res. 2004, 123, 63–74. [CrossRef] 47. Ling, K.A.; Power, S.A.; Ashmore, M.R. A Survey of the Health of Fagus sylvatica in Southern Britain. J. Appl. Ecol. 1993, 30, 295–306. [CrossRef] 48. Eichhorn, J.; Icke, R.; Isenberg, A.; Paar, U.; Schönfelder, E. Temporal development of crown condition of beech and oak as a response variable for integrated evaluations. Eur. J. For. Res. 2005, 124, 335–347. [CrossRef] 49. Ewald, J. Ecological background of crown condition, growth and nutritional status of Picea abies (L.) Karst. in the Bavarian Alps. Eur. J. For. Res. 2005, 124, 9–18. [CrossRef] 50. Graf Pannatier, E.; Dobbertin, M.; Schmitt, M.; Thimonier, A.; Waldner, P. Effects of the drought 2003 on forests in Swiss Level II plots. In Proceedings of the Symposium: Forests in a Changing Environment—Results of 20 Years ICP Forests Monitoring, Göttingen, Germany, 25–28 October 2006; pp. 125–135. 51. Bréda, N.; Huc, R.; Granier, A.; Dreyer, E. Temperate forest trees and stands under severe drought: A review of ecophysiological responses, adaptation processes and long-term consequences. Ann. For. Sci. 2006, 63, 625–644. [CrossRef] 52. Filewod, B.; Thomas, S.C. Impacts of a spring heat wave on canopy processes in a northern hardwood forest. Glob. Chang. Biol. 2014, 20, 360–371. [CrossRef] [PubMed] 53. Zweifel, R.; Steppe, K.; Sterck, F.J. Stomatal regulation by microclimate and tree water relations: Interpreting ecophysiological field data with a hydraulic plant model. J. Exp. Bot. 2007, 58, 2113–2131. [CrossRef] [PubMed] 54. Peñuelas, J.; Hunt, J.M.; Ogaya, R.; Jump, A.S. Twentieth century changes of tree-ring δ13C at the southern range-edge of Fagus sylvatica: Increasing water-use efficiency does not avoid the growth decline induced by warming at low altitudes. Glob. Chang. Biol. 2008, 14, 1076–1088. [CrossRef] 55. Allen, C.D.; Macalady, A.K.; Chenchouni, H.; Bachelet, D.; McDowell, N.; Vennetier, M.; Kitzberger, T.; Rigling, A.; Breshears, D.D.; Hogg, E.H.; et al. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For. Ecol. Manag. 2010, 259, 660–684. [CrossRef] 56. Ogaya, R.; Peñuelas, J.; Asensio, D.; Llusià, J. Chlorophyll fluorescence responses to temperature and water availability in two co-dominant Mediterranean shrub and tree species in a long-term field experiment simulating climate change. Environ. Exp. Bot. 2011, 73, 89–93. [CrossRef] 57. Dreesen, F.E.; De Boeck, H.J.; Janssens, I.A.; Nijs, I. Summer heat and drought extremes trigger unexpected changes in productivity of a temperate annual/biannual plant community. Environ. Exp. Bot. 2012, 79, 21–30. [CrossRef] 58. Shao, H.-B.; Chu, L.-Y.; Jaleel, C.A.; Zhao, C.-X. Water-deficit stress-induced anatomical changes in higher plants. Comptes Rendus Biol. 2008, 331, 215–225. [CrossRef] 59. Perry, T.O. Dormancy of Trees in Winter. Science 1971, 171, 29–36. [CrossRef] 60. McDowell, N.G.; Allen, C.D. Darcy’s law predicts widespread forest mortality under climate warming. Nat. Clim. Chang. 2015, 5, 669–672. [CrossRef] 61. Eschrich, W.; Burchardt, R.; Essiamah, S. The induction of sun and shade leaves of the European beech (Fagus sylvatica L.): Anatomical studies. Trees 1989, 3, 1–10. [CrossRef] 62. Uemura, A.; Ishida, A.; Nakano, T.; Terashima, I.; Tanabe, H.; Matsumoto, Y. Acclimation of leaf characteristics of Fagus species to previous-year and current-year solar irradiances. Tree Physiol. 2000, 20, 945–951. [CrossRef] [PubMed] 63. McDowell, N.; Pockman, W.T.; Allen, C.D.; Breshears, D.D.; Cobb, N.; Kolb, T.; Plaut, J.; Sperry, J.; West, A.; Williams, D.G.; et al. Mechanisms of plant survival and mortality during drought: Why do some plants survive while others succumb to drought? New Phytol. 2008, 178, 719–739. [CrossRef] http://doi.org/10.31298/sl.144.7-8.2 http://doi.org/10.1111/gcb.14326 http://doi.org/10.1198/jasa.2009.ap07058 http://doi.org/10.1007/s00704-014-1368-8 http://doi.org/10.1073/pnas.1010070108 http://www.ncbi.nlm.nih.gov/pubmed/21220333 http://doi.org/10.1023/A:1005157208701 http://doi.org/10.1007/s11270-014-1992-z http://doi.org/10.1007/s10342-007-0174-6 http://doi.org/10.1007/s10342-004-0014-x http://doi.org/10.2307/2404631 http://doi.org/10.1007/s10342-005-0097-z http://doi.org/10.1007/s10342-004-0051-5 http://doi.org/10.1051/forest:2006042 http://doi.org/10.1111/gcb.12354 http://www.ncbi.nlm.nih.gov/pubmed/24038752 http://doi.org/10.1093/jxb/erm050 http://www.ncbi.nlm.nih.gov/pubmed/17490998 http://doi.org/10.1111/j.1365-2486.2008.01563.x http://doi.org/10.1016/j.foreco.2009.09.001 http://doi.org/10.1016/j.envexpbot.2011.08.004 http://doi.org/10.1016/j.envexpbot.2012.01.005 http://doi.org/10.1016/j.crvi.2008.01.002 http://doi.org/10.1126/science.171.3966.29 http://doi.org/10.1038/nclimate2641 http://doi.org/10.1007/BF00202394 http://doi.org/10.1093/treephys/20.14.945 http://www.ncbi.nlm.nih.gov/pubmed/11303569 http://doi.org/10.1111/j.1469-8137.2008.02436.x Plants 2022, 11, 730 13 of 13 64. Eickenscheidt, N.; Augustin Nicole, H.; Wellbrock, N. Spatio-temporal modelling of forest monitoring data: Modelling German tree defoliation data collected between 1989 and 2015 for trend estimation and survey grid examination using GAMMs. iForest— Biogeosci. For. 2019, 12, 338–348. [CrossRef] 65. Ferretti, M.; Fischer, R.; Mues, V.; Granke, O.; Lorenz, M.; Seidling, W.; Nicolas, M. Part II: Basic design principles for the ICP Forests Monitoring Networks. Version 2020-2. In Manual on Methods and Criteria for Harmonized Sampling, Assessment, Onitoring and Analysis of the Effects of Air Pollution on Forests; Centre, U.I.F.P.C., Ed.; Thünen Institute of Forest Ecosystems: Eberswalde, Germany, 2020; p. 33 + Annex. 66. Caudullo, G.; Welk, E.; San-Miguel-Ayanz, J. Chorological maps for the main European woody species. Data Brief 2017, 12, 662–666. [CrossRef] [PubMed] 67. Pravilnik o Vrstama Stanišnih Tipova, Karti Staništa, Ugroženim i Rijetkim Stanišnim Tipovima te o Mjerama za Očuvanje Stanišnih Tipova NN 7/2006. Available online: https://narodne-novine.nn.hr/clanci/sluzbeni/2006_01_7_156.html (accessed on 24 January 2022). 68. Perčec Tadić, M. Gridded Croatian climatology for 1961–1990. Theor. Appl. Climatol. 2010, 102, 87–103. [CrossRef] 69. Perčec Tadić, M.; Pasarić, Z.; Guijarro, J.A. Croatian High-Resolution Monthly Gridded Data Set of Homogenised Surface Air Temperature (Manuscript submitted). Theor. Appl. Climatol. 2022. 70. Palmer, W.C. Meteorological Drought; US Department of Commerce, Weather Bureau: Washington, DC, USA, 1965; Volume 45, p. 58. 71. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A Multiscalar Drought Index Sensitive to Global Warming: The Standard- ized Precipitation Evapotranspiration Index. J. Clim. 2010, 23, 1696–1718. [CrossRef] 72. Sen, P.K. Estimates of the Regression Coefficient Based on Kendall’s Tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [CrossRef] 73. Mann, H.B. Nonparametric tests against trend. Econometrica 1945, 13, 245–259. [CrossRef] 74. Kendall, M.G. Rank Correlation Methods; Griffin: Oxford, UK, 1948. 75. Wilcox, R.R. Fundamentals of Modern Statistical Methods: Substantially Improving Power and Accuracy; Springer: New York, NY, USA, 2010. 76. Wang, J.-F.; Zhang, T.-L.; Fu, B.-J. A measure of spatial stratified heterogeneity. Ecol. Indic. 2016, 67, 250–256. [CrossRef] 77. Bivand, R.S.; Pebesma, E.J.; Gomez-Rubio, V.; Pebesma, E.J. Applied Spatial Data Analysis with R; Springer: New York, NY, USA, 2013; Volume 2. 78. Laird, N.M.; Ware, J.H. Random-Effects Models for Longitudinal Data. Biometrics 1982, 38, 963–974. [CrossRef] [PubMed] 79. LeBeau, B. Impact of serial correlation misspecification with the linear mixed model. J. Mod. Appl. Stat. Methods 2016, 15, 21. [CrossRef] 80. Goerg, G.M. The Lambert Way to Gaussianize Heavy-Tailed Data with the Inverse of Tukey’s h Transformation as a Special Case. Sci. World J. 2015, 2015, 909231. [CrossRef] 81. Johnson, J.B.; Omland, K.S. Model selection in ecology and evolution. Trends Ecol. Evol. 2004, 19, 101–108. [CrossRef] 82. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [CrossRef] 83. Thompson, C.G.; Kim, R.S.; Aloe, A.M.; Becker, B.J. Extracting the Variance Inflation Factor and Other Multicollinearity Diagnostics from Typical Regression Results. Basic Appl. Soc. Psychol. 2017, 39, 81–90. [CrossRef] 84. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control. 1974, 19, 716–723. [CrossRef] 85. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2016. 86. Shumway, R.H.; Stoffer, D.S.; Stoffer, D.S. Time Series Analysis and Its Applications; Springer: Berlin/Heidelberg, Germany, 2000; Volume 3. 87. EEA. European Digital Elevation Model (EU-DEM), Version 1.1, European Environment Agency (EEA). 2016. Available online: https://land.copernicus.eu/imagery-in-situ/eu-dem/eu-dem-v1.1?tab=metadata (accessed on 24 January 2022). 88. ISO-13878; Soil Quality—Determination of Total Nitrogen Content by Dry Combustion (“Elemental Analysis”). ISO: Geneva, Switzerland, 1998. 89. Egnér, H.; Riehm, H.; Domingo, W.R. Untersuchungen über die chemische Bodenanalyse als Grundlage für die Beurteilung des Nährstoffzustandes der Böden. II. Chemische Extraktionsmethoden zur Phosphor-und Kaliumbestimmung. K. Lantbr. Ann. 1960, 26, 199–215. 90. Škorić, A. Priručnik za Pedološka Istraživanja; Fakultet Poljoprivrednih Znanosti: Zagreb, Croatia, 1985. 91. ISO-10390; Soil Quality—Determination of pH. ISO: Geneva, Switzerland, 2005. http://doi.org/10.3832/ifor2932-012 http://doi.org/10.1016/j.dib.2017.05.007 http://www.ncbi.nlm.nih.gov/pubmed/28560272 https://narodne-novine.nn.hr/clanci/sluzbeni/2006_01_7_156.html http://doi.org/10.1007/s00704-009-0237-3 http://doi.org/10.1175/2009JCLI2909.1 http://doi.org/10.1080/01621459.1968.10480934 http://doi.org/10.2307/1907187 http://doi.org/10.1016/j.ecolind.2016.02.052 http://doi.org/10.2307/2529876 http://www.ncbi.nlm.nih.gov/pubmed/7168798 http://doi.org/10.22237/jmasm/1462076400 http://doi.org/10.1155/2015/909231 http://doi.org/10.1016/j.tree.2003.10.013 http://doi.org/10.1023/A:1010933404324 http://doi.org/10.1080/01973533.2016.1277529 http://doi.org/10.1109/TAC.1974.1100705 https://land.copernicus.eu/imagery-in-situ/eu-dem/eu-dem-v1.1?tab=metadata Ognjenovic et al 2021.pdf Ognjenović et al 2022 Plants Introduction Results Temporal Trends in Tree Vitality Influence of Environmental Conditions on Tree Vitality Discussion Temporal Trends in Crown Vitality Influence of Environmental Conditions on Tree Vitality Materials and Methods Study Area and Plot Selection Defoliation Assessment and Crown Vitality Indicators Soil Sampling and Analysis Climate Data Data Analysis References