Effect of Climate Change on the Growth of Endangered Scree Forests in Krkonoše National Park (Czech Republic)

: Scree forests with large numbers of protected plants and wildlife are seriously threatened by climate change due to more frequent drought episodes, which cause challenges for very stony, shallow soils. The effect of environmental factors on the radial growth of ﬁve tree species—European beech ( Fagus sylvatica L.), Norway spruce ( Picea abies (L.) Karst.), sycamore maple ( Acer pseudoplatanus L.), European ash ( Fraxinus excelsior L.), and mountain elm ( Ulmus glabra Huds.)—was studied in the mixed stands (105–157 years) in the western Krkonoše Mountains (Czech Republic) concerning climate change. These are communities of maple to ﬁr beechwoods (association Aceri-Fagetum sylvaticae and Luzulo-Abietetum albae ) on ranker soils at the altitude 590–700 m a.s.l. Production, structure, and biodiversity were evaluated in seven permanent research plots and the relationships of the radial growth (150 cores) to climatic parameters (precipitation, temperature, and extreme conditions) and air pollution (SO 2 , NO X , ozone exposure). The stand volume reached 557–814 m 3 ha − 1 with high production potential of spruce and ash. The radial growth of beech and spruce growing in relatively favorable habitat conditions (deeper soil proﬁle and less skeletal soils) has increased by 16.6%–46.1% in the last 20 years. By contrast, for sycamore and ash growing in more extreme soil conditions, the radial growth decreased by 12.5%–14.6%. However, growth variability increased (12.7%–29.5%) for all tree species, as did the occurrence of negative pointer years (extremely low radial growth) in the last two decades. The most sensitive tree species to climate and air pollution were spruce and beech compared to the resilience of sycamore and ash. Spectral analysis recorded the largest cyclical ﬂuctuations (especially the 12-year solar cycle) in spruce, while ash did not show any signiﬁcant cycle processes. The limiting factors of growth were droughts with high temperatures in the vegetation period for spruce and late frosts for beech. According to the degree of extreme habitat conditions, individual tree species thus respond appropriately to advancing climate change, especially to an increase in the mean temperature (by 2.1 ◦ C), unevenness in precipitation, and occurrence of extreme climate events in the last 60 years.


Introduction
Forest management practices have a considerable impact on forest ecosystems, especially in the context of global climate change for adaptive and mitigation measures for increasing stability of forests [1,2]. Lack of precipitation and, as a result, more frequent and severe droughts will increasingly be a major risk to the growth, ecological stability, and vitality of forest ecosystems during global climate change [3,4]. Higher air temperatures, in

Study Site
The area of interest is shallow grounded rubble in the western Krkonoše Mountains (west of Harrachov) on the slopes near the Jizera Valley at an altitude of 590-700 m ( Figure 1). This area is part of the II and III nature protection zones of the Krkonoše National Park. These are richly structured mixed stands that belong to special-purpose forests and are left to spontaneous development, or are managed by the selection and shelterwood systems. Stands are most typically upland communities of maple and fir beechwoods, wherein their natural composition, European beech had 40%, sycamore maple 20%, silver fir 20%, Norway spruce 10%, and the rest consisted of European ash and mountain elm. However, silver fir almost became extinct there due to the synergism of the air pollution calamity (high concentrations of SO 2 ) in the 1980s and 1990s and the attack of the silver fir woolly aphid (Dreyfusia nordmannianae) in the tree layer [58]. These are communities belonging mainly to the association Aceri-Fagetum sylvaticae, and also to the association Luzulo-Abietetum albae.
The mean annual temperature of the locality is 4.9 • C (10.7 • C in the vegetation period) and the annual total precipitation varies around 1200 mm (600 mm in the vegetation period), with maxima in July ( Figure 1). The length of the vegetation period is 122 days and Lang's rain factor is 245. The study area was qualified as a humid continental climate zone, characterized by warm to hot, humid summers and cold to severely cold winters (Dfb) as defined by the Köppen climate classification [59]. In terms of climate change, the average annual temperature has increased by 2.1 • C over the last 58 years (1961-2019) and the total annual precipitation has decreased by 15 mm, while fluctuations of monthly climatic data increased in both cases, especially in the case of precipitation. When comparing climatic extremes, 7 of the 10 warmest months were found in the last 15 years. Annual averages of SO 2 concentrations during the air pollution disaster (in the 1970s and 1980s) ranged 120-410 µg m −3 and led to a significant decrease after 2000 to the current values of ca 4 µg m −3 [60]. The NO X concentration currently averages around 8 µg m −3 , and ozone exposure (AOT40F) shows mean values of 27,000 ppb h −1 [61]. The mean annual temperature of the locality is 4.9 °C (10.7 °C in the vegetation period) and the annual total precipitation varies around 1200 mm (600 mm in the vegetation period), with maxima in July ( Figure 1). The length of the vegetation period is 122 days and Lang's rain factor is 245. The study area was qualified as a humid continental climate zone, characterized by warm to hot, humid summers and cold to severely cold winters (Dfb) as defined by the Köppen climate classification [59]. In terms of climate change, the average annual temperature has increased by 2.1 °C over the last 58 years (1961-2019) and the total annual precipitation has decreased by 15 mm, while fluctuations of monthly climatic data increased in both cases, especially in the case of precipitation. When comparing climatic extremes, 7 of the 10 warmest months were found in the last 15 years. Annual averages of SO2 concentrations during the air pollution disaster (in the 1970s and 1980s) ranged 120-410 µg m -3 and led to a significant decrease after 2000 to the current values of ca 4 µg m -3 [60]. The NOX concentration currently averages around 8 µg m -3 , and ozone exposure (AOT40F) shows mean values of 27,000 ppb h -1 [61].
The geological basis consists mainly of quartzites and phyllites. The predominant soil types are modal and cambic Rankers and ranker Cambisols. Depth and skeletality in scree forests are quite variable in a mosaic arrangement with numerous rocky outcrops. On very shallow soils (depth up to 30-40 cm) and most skeletal soils (skeleton content over 70%), there are mostly sycamore maple, ash, and elm. On relatively deeper (depth of soil profile over 50 cm) and less skeletal soils (skeleton content up to 55%), European beech and Norway spruce dominate.
A detailed overview of the basic site and stand data of permanent research plots (PRPs) is given in Table 1. The forest stands usually consist of two to three tree layers, with the age of the upper tree layer from 105 to 130 years. Given their management, the studied stands were divided into two groups. PRPs 1, 2, 5, and 7 have been left to develop spontaneously for more than 30 years, and PRPs 3, 4, and 6 are managed by near-natural methods. It is mainly the local introduction of silver fir, the proportion of which was reduced to a minimum during the air pollution disaster in the 1980s. The geological basis consists mainly of quartzites and phyllites. The predominant soil types are modal and cambic Rankers and ranker Cambisols. Depth and skeletality in scree forests are quite variable in a mosaic arrangement with numerous rocky outcrops. On very shallow soils (depth up to 30-40 cm) and most skeletal soils (skeleton content over 70%), there are mostly sycamore maple, ash, and elm. On relatively deeper (depth of soil profile over 50 cm) and less skeletal soils (skeleton content up to 55%), European beech and Norway spruce dominate.
A detailed overview of the basic site and stand data of permanent research plots (PRPs) is given in Table 1. The forest stands usually consist of two to three tree layers, with the age of the upper tree layer from 105 to 130 years. Given their management, the studied stands were divided into two groups. PRPs 1, 2, 5, and 7 have been left to develop spontaneously for more than 30 years, and PRPs 3, 4, and 6 are managed by near-natural methods. It is mainly the local introduction of silver fir, the proportion of which was reduced to a minimum during the air pollution disaster in the 1980s.

Data Collection
FieldMap technology (IFER-Monitoring and Mapping Solutions Ltd.) was used to determine the structure of the tree layer of forest stands on seven PRPs with a size of 50 × 50 m (0.25 ha) in 2020. The position of all individuals of the tree layer with a diameter at breast height over 4 cm was measured. The diameter limit was determined according to Slanař et al. [18], Bulušek et al. [36], Král et al. [60], and other researches, where individuals under 4 cm are included to natural regeneration. Tree heights, heights of the live crown base, and crown projection area were also measured in the tree layer, in at least 4 directions perpendicular to each other. Diameters of the tree layer were measured by a Mantax Blue metal caliper (Haglöf, Sweden) with an accuracy to 1 mm, and the heights were measured using a Vertex laser hypsometer (Haglöf, Sweden) with an accuracy to 0.1 m.
For dendrochronological analysis of the samples, increment cores were taken using a Pressler auger perpendicular to the axis of the trunk at a height of 1.3 m above the ground. Thirty core samples were randomly (RNG Excel function) taken from healthy predominant and dominant trees (according to Kraft/1884/classification) of the five following tree species: European beech, Norway spruce, sycamore maple, European ash, and mountain elm. A total of 150 cores were taken for dendrochronology (135 cores were used for analysis). Increment cores were measured with an Olympus microscope with a LINTAB measuring table (Rinntech, Heidelberg, Germany). The measurement was performed with an accuracy of 0.01 mm, and TSAP-Win software (Rinntech, Heidelberg, Germany) was used to record the increment cores. Measurements were made from the bark towards the core perpendicular to the center of the trunk so that each annual ring was measured perpendicular to the trunk axis.
Data sets from meteorological stations and air pollution monitoring stations were used to derive stress factors related to climate and air pollution. Available data from the Desná-Souš station (772 m above sea level; GPS 50 • 47 23 N, 15 • 19 9 E), which is approximately 6 km away from the monitored PRPs, was used to analyze the air pollution situation according to SO 2 (1970SO 2 ( -2012 and NO x (1994-2019) concentrations, and AOT40F (1993-2019). Mean and maximum values of concentrations in µg m −3 and ppb h −1 and their 95 quantiles were used for evaluation. Climate behavior in terms of temperature and precipitation ratios was evaluated based on the data from the meteorological station Harrachov from the years 1961-2019 (690 m above sea level; GPS 50 • 46 29 N, 15 • 25 39 E, located about 3 km from the monitored PRPs). The development of temperature and precipitation conditions was studied on the basis of data on mean annual temperature, temperature in the vegetation period, temperature outside the vegetation period, temperature in individual months, temperature variability, minima and maxima, and an annual sum of precipitation, precipitation in the vegetation period, precipitation outside the vegetation period, precipitation in individual months, precipitation variability, minima, and maxima. The study area have been long-term monitored by various specialists in the fields, and damage by other climatic factors (wind, late frost) and pests (bark beetles, Cryptococcus fagisuga, Ophiostoma ulmi, Rhytisma acerinum) has been recorded.

Data Analysis
The basic structure, biodiversity, and production characteristics of the tree layer were evaluated by the SIBYLA 5.1. software [62]. The stand volume was calculated according to Petráš and Pajtík [63]. The relative stand density index (SDI) [64], the crown closure (CC) [65], and the crown projection area (CPA) were observed for each plot. The maximum SDI value was derived from the model of the yield tables [66]. Slenderness coefficient (height to diameter ratio-HDR) was evaluated as an indicator of static stand stability [67,68] according to the following range: excellent (HDR < 82), good (83)(84)(85)(86)(87)(88)(89)(90)(91)(92), satisfactory (93)(94)(95)(96)(97)(98)(99)(100)(101), and dissatisfactory (>102) [69]. Species diversity was evaluated by species richness [70], species heterogeneity [71], and species evenness [72]. In terms of structure, horizontal structure [73], Arten-profile index [74], vertical diversity, crown differentiation [75], and diameter and height differentiation [76] were calculated. Moreover, total stand diversity [75] was determined based on tree species diversity, diversity of vertical structure, diversity of tree spatial distribution, and diversity of crown differentiation ( Table 2). For calculation of horizontal structure, the software PointPro 2.1 (CULS Prague, Zahradník) was used. The Monte Carlo technique was used to simulate random data sets. For each investigated plot, 4000 simulations of the Poisson forest with the same stand density as the investigated stand were performed. When the value of the aggregation index exceeds the upper or lower limit of the random distribution interval (at significance 0.05), the horizontal structure of trees is regular or aggregated [77]. Cross-dating the increment cores was performed in the CDendro program so that the cross-correlation index was CC > 25 for each sample [78]. Dendrochronological data were processed using R software [79] using the "dplr" package [80][81][82]. The detrending of each tree was performed negatively by exponential detrending with an interleaved spline of 1/3 of the age of each tree. This detrending removes the age trend while maintaining low-frequency climate signals [83]. The EPS (expressed population signal) was calculated for detrended data. The EPS represents the reliability of a chronology as a fraction of the joint variance of the theoretical infinite tree population. The limit for using the data for comparison against climatic data was a significant EPS limit so that EPS > 0.85 [80]. The SNR (signal to noise ratio), which represents the signal strength of chronology, and R-bar (inter-series correlations) [84] were also calculated. The first-order autocorrelation (Ar1) was also carried out. The EPS, SNR, R-bar, Ar1 were calculated by the instructions to "dplr" [80] based on common dendrochronological theories [84,85].
The pointer years analysis of the relative growth change was carried out [86] to determine changes in tree-rings. The pointer year is created by normalizing the tree-ring series in a moving window with each year placed as a central point. The pointer years reflect the number of standard deviations from the local mean of the average ring width of the previous four years. The pointer years identify event years where pointer year > 0.75 standard deviation of the previous four years. The threshold on the percentage of trees in a negative or positive event year was used for better clarity [87]. The pointer years or mean annual growth deviation (in percentage) are distinguished for the most common event-year class.
Spectral analyses for the detrended radial increment of each studied tree species were performed with Statistica 13 software (StatSoft, Tulsa, OK, USA). The calculation was performed with the "Single Fourier (Spectral) Analysis" function, using the output "Periodogram" plot by Period".
The average ring series from individual tree species were correlated with climatic data (precipitation, temperature) and air pollution data (AOT40F, SO 2 , and NO x concentrations) in Statistica 13 software (StatSoft, Tulsa, OK, USA). DendroClim software [88] was used to model the radial increment depending on the monthly climatic characteristics.
Unconstrained principal component analysis (PCA) in Canoco 5 [89] was used to analyze the relationships between climate factors, air pollution load, and the radial growth of tree species. Data were log-transformed, centered, and standardized before the analysis. The results of PCA were exported into the form of an ordination diagram. The situation map was made in ArcGIS 10 software (Esri, West Redlands, CA, USA).

Stand Structure and Production
The number of live trees in the tree layer in 2020 ranged between 260 and 480 trees ha −1 (on average 328 trees ha −1 ); for stands left to spontaneous development, the average was 296 trees ha −1 , and for managed stands, 372 trees ha −1 ( Table 3). The stand density index ranged from 0.60 to 0.81 (on average 0.68). The canopy closure ranged from 76.6 to 97.4%. Considering the rather extreme habitat, the stand production capacity was at a relatively high level. The average basal area reached 46.2 m 2 ha −1 for stands left to spontaneous development and 45.8 m 2 ha −1 for managed stands. The highest stand volume was 814 m 3 ha −1 on PRP 7, and the smallest was 597 m 3 ha −1 on PRP 3, where overhead release felling was performed. In stands left to spontaneous development, the average stand volume was 730 m 3 ha −1 , and in managed stands, it was lower by 17.9% (599 m 3 ha −1 ). The mean annual increment varies in the range of 4.28-6.06 m 3 ha −1 yr −1 ; for unmanaged stands, it was 5.78 m 3 ha −1 yr −1 , and in managed stands, 4.81 m 3 ha −1 yr −1 . The slenderness ratio-within the limits for the given habitat and stand conditions-showed excellent static stability of these stands (HDR 62.6-74.2). Table 3. Structural characteristics of tree layer on permanent research plots 1-7 in the western Krkonoše Mountains in 2020.

PRP
Age Notes: Age-mean stand age, DBH-mean quadratic diameter at breast height, h-mean tree height, f-form factor, v-mean tree volume, N-number of trees per hectare, BA-basal area, V-stand volume, HDR-height to diameter ratio (slenderness quotient), MAI-mean annual increment, CC-canopy closure, CPA-crown projection area, SDI-stand density index.
The diameter distribution and relationship between DBH and height of the tree layer are clearly shown in Figure 2. Three to four stories were measured in stands left to spontaneous development and two to three levels in managed stands. The tree layers of all PRPs were dominated by beech followed by spruce and maple. On PRP 1, 2, and 7, elm was admixed, while it was ash on PRP 4 and 6. Fir was interspersed only on PRP 2 and 3. All tree species except the fir are represented in all stories. Only the fir is represented in the lowest story. The diameter distribution on the left side is closer to the shape of the Gaussian curve with the highest frequency of trees in diameter class 48-56 cm. The highest tree was observed in the case of spruce (42.6 m) followed by beech (41.0 m), while it was lowest in the case of sycamore (33.0 m) and ash (33.7 m).
PRPs were dominated by beech followed by spruce and maple. On PRP 1, 2, and 7, elm was admixed, while it was ash on PRP 4 and 6. Fir was interspersed only on PRP 2 and 3. All tree species except the fir are represented in all stories. Only the fir is represented in the lowest story. The diameter distribution on the left side is closer to the shape of the Gaussian curve with the highest frequency of trees in diameter class 48-56 cm. The highest tree was observed in the case of spruce (42.6 m) followed by beech (41.0 m), while it was lowest in the case of sycamore (33.0 m) and ash (33.7 m).

Tree Layer Biodiversity
In terms of species diversity, the species richness of the tree layer was low to medium, as was species heterogeneity ( Table 4). The species evenness ranged from medium to very high. The horizontal structure of the tree layer was significantly (α = 0.05) aggregated on PRP 1 and 7, on other PRPs, the spatial pattern was random with a prevailing tendency to aggregated distribution of trees. According to the standardized A index, the vertical structure of the tree layer was moderately diversified on most PRPs, only on PRP 6 was it strongly diversified. However, the diversity of the vertical structure according to B index was high on all PRPs except PRP 2, where it was medium. The diameter differentiation was medium on all PRPs, and the height differentiation was low to medium. Crown differentiation was medium on all PRPs. In terms of the total diversity of the tree layer, uneven stand structure was observed on all PRPs, while diversity was higher by 1.3% (B = 6.555) in unmanaged compared to managed stands (B = 6.469).

Tree Layer Biodiversity
In terms of species diversity, the species richness of the tree layer was low to medium, as was species heterogeneity ( Table 4). The species evenness ranged from medium to very high. The horizontal structure of the tree layer was significantly (α = 0.05) aggregated on PRP 1 and 7, on other PRPs, the spatial pattern was random with a prevailing tendency to aggregated distribution of trees. According to the standardized A index, the vertical structure of the tree layer was moderately diversified on most PRPs, only on PRP 6 was it strongly diversified. However, the diversity of the vertical structure according to B index was high on all PRPs except PRP 2, where it was medium. The diameter differentiation was medium on all PRPs, and the height differentiation was low to medium. Crown differentiation was medium on all PRPs. In terms of the total diversity of the tree layer, uneven stand structure was observed on all PRPs, while diversity was higher by 1.3% (B = 6.555) in unmanaged compared to managed stands (B = 6.469).

Dynamics of Radial Growth
The dendrochronological analyses show that Norway spruce (RW = 2.41 mm) reached the largest average diameter increment (Table 5). Tree species such as mountain elm, European ash, and sycamore maple had similar RW ranging from 1.96 to 2.09 mm. The lowest diameter increment was recorded for European beech (RW = 1.43 mm). Comparing the period 2010-2019 to 1990-2009, radial growth increased in spruce by 46.1% and in beech by 20.6%, while it decreased in ash by 14.9% and sycamore by 12.5% (no changes in the case of elm-0.1%). On the other hand, the annual variability in growth increased in all tree species, from 12.7% (in elm) to 29.5% (in ash) in the last 20 years. The oldest tree recorded in the studied area was a 152-year-old beech. Other tree species showed a maximum age from 120 to 145 years. Spruce reached higher autocorrelation values (Ar1 = 0.76) than European beech (Ar1 = 0.68). However, other admixed and less represented tree species  The standardized annual ring curves showed different courses of growth in tree species from 1960 to 2019 ( Figure 3). The largest fluctuations in the radial growth were in spruce and elm, followed by beech compared to low variability in growth as in the case of ash and sycamore. In spruce, the negative decrease of the radial growth in 1972 was caused by a lack of precipitation (62% of precipitation of the long-term average 1961-2019) until 1974. In the period from 1970 to 1985, the growth of spruce was negatively affected by air pollution load. The effect of air pollution in synergism with winter drying in spruce was most prominent in the years 1980-1984. In early spring 1995, due to a thick layer of wet snow, some of the spruce trees were damaged by top breaks, and the following year, 1996, was very cold (75% of the average), which influenced the decrease in the radial growth. Long-term droughts in the vegetation period were recorded in 2006 (47% of the average in June and July). The significantly declining radial growth in 2018-2019 was due to a substantially reduced sum of precipitation (74% of the average) and high average temperature (126% of the average) during both years, and especially in the vegetation period.
In beech, the years with a high decrease in the radial growth were 1979, 1981, 2011, and 2016. On a smaller scale, European beech also reacted to air pollution in the RWI between 1979 and 1986. In the spring of 1979 and 1981, the assimilation apparatus was damaged by late frosts. In 1985, there was a strong outbreak of beech scale (Cryptococcus fagisuga). In 1995, beech crowns were severely damaged by rime, and in 1996, late frosts recurred. Very strong damage to the beech assimilation apparatus occurred at the beginning of the vegetation period in 2011, and milder in 2016. The sharply declining radial growth in 2016-2019 was the same in the case of spruce. For elm, in the period 1972-1974, there was a weaker incidence of Ophiostoma ulmi, and in 1985-1995, there was a strong incidence of Ophiostoma ulmi in synergism with increased air pollution. In 2004-2007 and 2009-2011, Ophiostoma ulmi was recorded there again. In 2014, precipitation was significantly below normal (74% of the long-term average). In ash, a decrease in the radial growth in 1996 was caused by a very cold year in combination with damage to the assimilation apparatus by a late frost, as in the case of beech. In sycamore, the reduced radial growth in 1993 and 1996 was caused by infestation with Rhytisma acerinum, and in 2011, there was damage to the assimilation apparatus by late frosts.
wet snow, some of the spruce trees were damaged by top breaks, and the following year, 1996, was very cold (75% of the average), which influenced the decrease in the radial growth. Long-term droughts in the vegetation period were recorded in 2006 (47% of the average in June and July). The significantly declining radial growth in 2018-2019 was due to a substantially reduced sum of precipitation (74% of the average) and high average temperature (126% of the average) during both years, and especially in the vegetation period. In beech, the years with a high decrease in the radial growth were 1979, 1981, 2011, and 2016. On a smaller scale, European beech also reacted to air pollution in the RWI between 1979 and 1986. In the spring of 1979 and 1981, the assimilation apparatus was damaged by late frosts. In 1985, there was a strong outbreak of beech scale (Cryptococcus Relative changes in the growth of annual rings and negative pointer years are recorded in Figure 4. The most stable fluctuations in the mean growth deviation are recorded for spruce, where positive and negative growth periods alternate after four to seven years. Spruce recorded two significantly negative years in 2018 and 2019. It is 2018 and 2019 that are abnormally high in terms of temperature averages. The average temperature over the last ten years is 6.5 • C, except in the years 2018 and 2019, it shows 7.3 and 7.2 • C, which is almost 0.9 • C above normal. At the same time, in 2018 and 2019, the total precipitation decreased. The average precipitation for the last 10 years was 1190 mm, but the years 2018 and 2019 show only 916 mm and 1008 mm. In contrast to the growth reactions of spruce, beech showed a more positive average growth deviation, and only one significant negative year was recorded in 2011, caused by the aforementioned significant spring frosts, which damaged beech sprouts. Pointer years of spruce and beech are easily visible from the decline in growth in Figure 3. Other tree species (elm, ash, sycamore) do not indicate any significant pointer years in growth.
Spruce recorded two significantly negative years in 2018 and 2019. It is 2018 and 2019 that are abnormally high in terms of temperature averages. The average temperature over the last ten years is 6.5 °C, except in the years 2018 and 2019, it shows 7.3 and 7.2 °C, which is almost 0.9 °C above normal. At the same time, in 2018 and 2019, the total precipitation decreased. The average precipitation for the last 10 years was 1190 mm, but the years 2018 and 2019 show only 916 mm and 1008 mm. In contrast to the growth reactions of spruce, beech showed a more positive average growth deviation, and only one significant negative year was recorded in 2011, caused by the aforementioned significant spring frosts, which damaged beech sprouts. Pointer years of spruce and beech are easily visible from the decline in growth in Figure 3. Other tree species (elm, ash, sycamore) do not indicate any significant pointer years in growth.

Spectral Analysis in Radial Growth
Recurring processes in the radial growth of tree species are clearly characterized by a single spectral analysis ( Figure 5), which describes the repeated period and the relative

Spectral Analysis in Radial Growth
Recurring processes in the radial growth of tree species are clearly characterized by a single spectral analysis ( Figure 5), which describes the repeated period and the relative magnitude of individual cycles over the periodogram values. Norway spruce, the most represented tree species, goes through the most prominent cycles reaching 1.2 periodogram values. Mountain elm also had large cycles in radial growth (up to 0.6 periodogram values). European beech indicated smaller cycles that reach up to 0.3 periodogram values. The smallest cyclic processes in growth (periodogram values up to 0.1) were recorded for European ash and sycamore maple. Spruce recorded 12, 20, and 60-year cycles in the radial growth, beech, 7 and 20-year cycles, elm, 15-year cycles, and sycamore maple, 20-year cycles. Ash did not show any significant cyclic processes in spectral analysis.
Spectral analysis recorded the largest cyclical fluctuations in spruce, while it was also found that with a lower proportion of tree species, cyclical fluctuations in the radial increment decreased. Across tree species, 7, 12, 20, and 60-year cycles in the radial increment were recorded in ring widths. gram values. Mountain elm also had large cycles in radial growth (up to 0.6 periodogram values). European beech indicated smaller cycles that reach up to 0.3 periodogram values. The smallest cyclic processes in growth (periodogram values up to 0.1) were recorded for European ash and sycamore maple. Spruce recorded 12, 20, and 60-year cycles in the radial growth, beech, 7 and 20-year cycles, elm, 15-year cycles, and sycamore maple, 20-year cycles. Ash did not show any significant cyclic processes in spectral analysis. Spectral analysis recorded the largest cyclical fluctuations in spruce, while it was also found that with a lower proportion of tree species, cyclical fluctuations in the radial increment decreased. Across tree species, 7, 12, 20, and 60-year cycles in the radial increment were recorded in ring widths.

Effect of Climate on Radial Growth
The highest effect of climate factor (monthly temperature and precipitation) was observed in European beech (10 significant months) followed by Norway spruce (8 months), while the most climate-resistant tree species were European ash (3 months) and sycamore maple (4 months; Figure 6). In beech, the temperature had only a positive effect and precipitation only a negative effect on the radial growth. It was not so obvious for other tree species. In spruce, the monthly climate factors in the previous year had a higher effect on the radial growth compared to the current year, while there was no significant (p < 0.05) effect of the previous year in the case of ash.
while the most climate-resistant tree species were European ash (3 months) and sycamore maple (4 months; Figure 6). In beech, the temperature had only a positive effect and precipitation only a negative effect on the radial growth. It was not so obvious for other tree species. In spruce, the monthly climate factors in the previous year had a higher effect on the radial growth compared to the current year, while there was no significant (p < 0.05) effect of the previous year in the case of ash. Generally, temperature had a higher effect on the radial growth compared to precipitation in relation to individual months ( Figure 6). The most significant months affecting the radial growth were April and July in the current years. July was also the most important month in the previous year. The highest negative correlation (r = −0.36) with the temperature was observed in the radial growth of ash in December of the previous year, or in July of the current year (r = −0.35) in the case of precipitation.

Interaction Between Climate, Air pollution, and Radial Growth of Tree Species
The results of PCA are presented in an ordination diagram showing relationships among radial growth, climatic factors (temperature and precipitation), air pollution, and tree species (Figure 7). The x-axis illustrates the concentrations of SO2 and radial growth of spruce, and the y-axis represents the temperature variability (standard deviation of monthly data). The mean SO2 concentration and mean temperature outside the vegetation Generally, temperature had a higher effect on the radial growth compared to precipitation in relation to individual months ( Figure 6). The most significant months affecting the radial growth were April and July in the current years. July was also the most important month in the previous year. The highest negative correlation (r = −0.36) with the temperature was observed in the radial growth of ash in December of the previous year, or in July of the current year (r = −0.35) in the case of precipitation.

Interaction between Climate, Air pollution, and Radial Growth of Tree Species
The results of PCA are presented in an ordination diagram showing relationships among radial growth, climatic factors (temperature and precipitation), air pollution, and tree species (Figure 7). The x-axis illustrates the concentrations of SO 2 and radial growth of spruce, and the y-axis represents the temperature variability (standard deviation of monthly data). The mean SO 2 concentration and mean temperature outside the vegetation period were the highest explanatory variables concerning the radial growth of tree species, while it was the lowest for the sum of precipitation outside the vegetation period and an annual sum of precipitation. In terms of climatic factors, temperature had a higher effect on the radial growth of the studied tree species compared to precipitation. The growth of spruce and elm was the most sensitive tree species to climatic factors, while a low effect was found for maple and ash. period were the highest explanatory variables concerning the radial growth of tree species, while it was the lowest for the sum of precipitation outside the vegetation period and an annual sum of precipitation. In terms of climatic factors, temperature had a higher effect on the radial growth of the studied tree species compared to precipitation. The growth of spruce and elm was the most sensitive tree species to climatic factors, while a low effect was found for maple and ash. According to individual tree species, the growth of beech is strongly affected by temperatures: significant (p < 0.01) positive effect of the mean annual temperature in the previous year and the mean temperature outside vegetation period of the current year (Table  6). Similarly, a significant (p < 0.05) positive effect of the mean temperature in June and July was observed in ash, while a significant (p < 0.05) negative effect of the mean temperature outside the vegetation period was observed only in elm. In terms of climatic variability, high differences in monthly values in temperature had a significant positive effect on growth of (p < 0.01) ash and (p < 0.05) elm, and negative in (p < 0.05) ash in the case of precipitation variability. No significant effect of cumulative climatic factors was observed in maple and spruce; however, spruce was a sensitive tree species to air pollution load. This tree species was significantly (p < 0.05-0.01) negatively correlated with mean and max SO2 concentrations, max NOX concentrations, and also with both ozone exposure indicators. However, to a lesser extent, beech showed a significant (p < 0.05) negative effect of the mean and max AOT40F and max NOX concentrations. Other tree species were not influenced by air pollution load. Table 6. Correlation matrix describing interactions between the radial growth of Norway spruce, European beech, mountain elm, European ash, and sycamore maple, as well as precipitation and According to individual tree species, the growth of beech is strongly affected by temperatures: significant (p < 0.01) positive effect of the mean annual temperature in the previous year and the mean temperature outside vegetation period of the current year (Table 6). Similarly, a significant (p < 0.05) positive effect of the mean temperature in June and July was observed in ash, while a significant (p < 0.05) negative effect of the mean temperature outside the vegetation period was observed only in elm. In terms of climatic variability, high differences in monthly values in temperature had a significant positive effect on growth of (p < 0.01) ash and (p < 0.05) elm, and negative in (p < 0.05) ash in the case of precipitation variability. No significant effect of cumulative climatic factors was observed in maple and spruce; however, spruce was a sensitive tree species to air pollution load. This tree species was significantly (p < 0.05-0.01) negatively correlated with mean and max SO 2 concentrations, max NO X concentrations, and also with both ozone exposure indicators. However, to a lesser extent, beech showed a significant (p < 0.05) negative effect of the mean and max AOT40F and max NO X concentrations. Other tree species were not influenced by air pollution load. Table 6. Correlation matrix describing interactions between the radial growth of Norway spruce, European beech, mountain elm, European ash, and sycamore maple, as well as precipitation and temperature (1961-2019), and concentrations of SO 2 , NO X (1994-2019), and AOT40F (1993-2019) in the western Krkonoše Mountains; significant correlations are designated by * (p < 0.05) and ** (p < 0.01).

Stand Structure and Biodiversity
In the studied mixed forests, aged 105-151, the stand volume ranged from 597 to 814 m 3 ha −1 , whereas the mean stand volume was 730 m 3 ha −1 in unmanaged stands and 599 m 3 ha −1 in managed stands (−17.9%). A lower stand volume was observed in mixed scree forests in the eastern (591-699 m 3 ha −1 ; aged 104-158 years) and middle parts (145-195 m 3 ha −1 ; aged 83-107 years) of the Czech Republic [33,90]. Concerning the connection of the stand volume and climate change, many researchers supported an assumption that mixed forest stands provide higher timber production and more resistance to increasing drought periods [61,91], whereas lack of precipitation was also a problem in recent years in our study area. Conversely, comparing the total diversity of unmanaged (B = 6.555) and managed (B = 6.469) forests, we have found no difference (±1.3%). Hájek et al. [34] documented a significantly higher total diversity in a comparable site in the Broumovsko Protected Area (B = 7.088-11.277).
In terms of species diversity, the species richness and heterogeneity of the tree layer was low to medium; however, the species evenness ranged from medium to very high in study forest stands. The number of tree species ranged from three to five species per PRPs.
Hájek et al. [34] found a significantly higher species diversity in mixed scree forests and counted from 5 to 13 tree species in lower altitudes. In terms of horizontal structure, the distribution of the tree layer was significantly aggregated on two PRPs, whereas on other PRPs, the spatial pattern was random with a prevailing tendency to aggregation. In natural forests, trees are arranged randomly or in aggregates [92,93], while in managed forests, trees are more regularly distributed [94,95]. Clustering of trees in natural forests can be caused by an uneven distribution of resources (especially soil depth, skeletality, light availability), interspecific competition, seed dispersal, or management and disturbance histories [96,97]. Even after many years since the end of stand management, the structure of forest stands may retain some consequential characteristics from previous management methods [98,99], as changes in stands caused by management lead to changes in interactions among trees and within the forest structure [100,101]. Studies of stand structure in scree and ravine forests have documented both random and aggregated tree arrangements [33,56]. Clustered tree arrangements have been observed, especially in extreme, strongly skeletal, and shallow soils [57], as in our case. A key factor in these forest stands is slope erosion caused by unstable debris beds [41]. For these reasons, most scree forests in the Krkonoše National Park have been left to spontaneous development or managed by fine near-nature practices.
In our case, in the managed stands, the number of trees was higher by 25.7% (372 trees ha −1 ) compared to stands both left to spontaneous development (296 trees ha −1 ) and managed stands. For example, Heiri et al. [98] showed that the number of large trees increased within a few decades since the end of management, but the total number of trees had decreased. Burrascano et al. [102] found that the higher density of large live trees and lower total tree density is conclusive in natural forests compared to mature managed forests as in our case. The density of trees in scree forests is generally lower than in the surrounding stands of a given altitudinal vegetation zone [54,103]. The lack of large trees in managed forests is often compensated by the higher overall density of trees in the stand [104,105], but this is not always the case [106,107].

Effect of Non-Climatic Factors to Radial Growth
On a European scale, the increased growth of forest ecosystems has been evident since the 1970s, often attributed to rising temperatures combined with increased nitrogen deposition and increasing CO 2 in the atmosphere [22,108]. On the other hand, other researchers have shown a decrease in growth [109,110]. From our results, the increase in growth was evident in Norway spruce (46.1%) and European beech (20.6%) in the last 20 years, but was not demonstrated in mountain elm (0.1%), sycamore maple (−12.5%), and European ash (−14.9%), which grow prevailingly in the most extreme conditions of shallow grounded rubble (with the largest skeletality and the shallowest soil profile). In the Czech Republic, however, during this period, the increment of stands was also significantly limited by several negative factors [39,58] including long-term droughts [111,112]. From our results, the influence of negative factors on the radial growth of the studied tree species was clearly visible.
In the Czech Republic, beech is negatively influenced by beech scale, which was reported, e.g., from eastern Krkonoše Mts. by Šimůnek et al. [113] or from Orlické hory Mts. [97]. Spruce is mostly infected by numerous insect pests (especially spruce bark beetle, larch bud moth), hoofed game, and fungal pathogens [114,115]. It was similar in our study, especially in the case of significant negative effects of SO 2 , NO x , and AOT40F concentrations on the growth of spruce, which was the most sensitive tree species of all of the studied samples. In sycamore, we proved the occurrence of only Rhytisma acerinum. The same was true in the scree forests in the Broumov region [33]. In our study, sycamore showed the highest stability in radial growth and resistance to adverse factors. In ash, the situation is similar in the eastern Krkonoše Mts., where fungal pathogens were confirmed [116] as in our study area. Moreover, Ophiostoma ulmi was documented in elm where the most significant damage occurred in the years 1985 to 1997, which is also evidenced by the significantly reduced radial growth. A similar development was reported from various parts of Europe [117,118].
Mina et al. [119] then showed that the growth of individual tree species in mixtures is greater than in monocultures. In general, it can be said that mixed stands are also considered to be more stable because the risk of their disturbance is distributed over a wider range of species-specific properties [120,121]. Similar results regarding the positive effect of mixing on growth were reported by Toïgo et al. [26]. On the other hand, Nothdurft and Engel [122] focused on individual trees in their research and do not share these conclusions. Vannoppen et al. [108] then state that mixed beech and oak stands have a positive effect on beech growth and, conversely, there is a negative effect on oak growth, which justifies the differences in the ecological demands of the tree species. Forrester et al. [123] claim better growth in mixed stands compared to monocultures by more efficient use of underground resources rather than above-ground ones. In mixed stands, better decomposition of the litter as well as the cycle of mineral nutrients may also be a positive factor [124,125]. The positive effects of mixing are enhanced by the ability of individual tree species to create different root systems [126], as in our case-the root structures vary from superficial through substitute tap to tap roots. The variability of the radial growth of the studied tree species corresponds in general to the type of root system and the character of soil conditions on shallow grounded rubble. This thesis is supported by the results from monospecific stands, where the variability of increment fluctuates considerably depending on different habitat and stand ratios [58]. The positive effect of mixing on the growth depends heavily on climatic and habitat conditions [127].

Effect of Climate to Radial Growth
Climatic factors showed considerable variability during the studied period. In the beginning, the increment was influenced more by cold or frost damage and vice versa at the end of the observed period by higher temperatures and lack of precipitation, which is basically in line with climate dynamics. From the southern regions of Europe, Begović et al. [46] offer a parallel in the positive relationship to summer precipitation and the negative effect of summer temperature. Similar results were obtained by other works [128,129] that deal with climate change. Zimmermann et al. [130], based on their analysis of the sensitivity of growth to climate, showed that both high summer temperatures and low summer precipitation are among the main factors limiting growth. Our results show a large variability affected by winter temperatures, although other studies have shown a positive relationship between winter temperature and the increment [131,132]. Carrer et al. [133] showed a decreasing effect of maximum winter temperatures on the radial growth from alpine to subtropical conditions in Italy.
The effects of climate change depend on the tree species [134]. Pretzsch et al. [135] showed that spruce responds to drought by a more pronounced decrease in increment, and significantly negative effects of high temperatures on spruce increment have also been demonstrated [136]. In addition, beech can tolerate sporadic dry summers relatively well, but once they repeat immediately, one after the other, its growth decreases significantly as well [130]. From our results, this is especially evident in recent years (2017-2019), which is in line with the results of Ols et al. [137]. In 1961-2010, the precipitation balance in the individual decades was relatively balanced in our case. Contrary to this, in the years 2010-2019, the average annual total precipitation decreased by 12%, but increased during the vegetation period. At the same time, each decennium warmed by an average of 0.3-0.4 • C (for 1961-2019, it was a total of 1.8 • C).
In our study, the monthly climate factors in the previous year had a higher effect on the radial growth in spruce and the current year in ash. Based upon long-term research, Vitali et al. [138] confirmed a positive correlation of increment with summer precipitation at lower altitudes. Conversely, their study does not present significant negative correlations with summer temperatures at higher altitudes. At higher altitudes, on the other hand, higher temperatures can have a positive effect on growth due to the earlier onset of cambial activity and xylem cell differentiation [139,140]. It is generally accepted that the onset of tree growth after winter rest is markedly sensitive to temperature [141,142]. Moreover, severe frosts can cause a significant reduction in growth [131].
Individual species are affected differently by lack of precipitation, high temperatures, or a combination of both, keeping in mind that they also have different growth reactions [143]. In addition, growth reactions may be specific in site and species [144]. Based on this, it is possible to explain the relatively significant differentiation of our results both by the microhabitat conditions, the species composition, and the different types of root systems. We observed the highest effect of climate factors in beech and spruce, while relatively climate-resistant tree species were ash and sycamore. Tree species may also have different growth responses in terms of their water management strategy, where isohydric species such as spruce close their stomata immediately during drought to maintain a consistent minimum water leaf potential-to prevent high water losses and disruption of water balance in the plant [91]. However, this reaction can lead to reduced photosynthetic activity and an overall reduction in growth [145,146]. On the other hand, anisohydric species such as beech close their stomata slowly and thereby may be exposed to hydrologi-cal deficiency [147], although photosynthetic activity may continue [148]. In addition, at higher temperatures, there are greater respiratory carbon losses, which lead to reduced carbohydrate availability and may be the cause of reduced growth [149]. In general, trees can handle high temperatures for some time if they have enough water resources to reduce stress conditions [150].
Drought adaptation strategies of trees also differ by tree species [151]. It can lead to better maintenance of osmoregulatory and hydrological functions [152]. Also, mixing deciduous and non-deciduous tree species can lead to better growth, thanks to the diversification of environmental resources, which, eventually, can lead to reduced stress [153]. In our case, for example, spruce can benefit from a greater transmission of light and water through the crowns of other tree species in the spring and winter months, when deciduous trees do not have leaves. Rötzer et al. [154] describe this profitable alliance between spruce and beech. Last but not least, the different type of root system also contributes to increased resource efficiency in stands [155]. Higher spring temperatures are also advantageous for fir and beech due to the reduction in the occurrence of late frosts to which they are sensitive [22,131] and the acceleration of cambial activity [156,157]. Mild winter conditions can also contribute to improved growth [158]. Conversely, the correlation of growth and precipitation at the middle and higher altitudes is usually not as significant as the correlation of growth and temperature, especially in the vegetation period [133,159]. Generally, the most significant months (relating to climatic factors) affecting the cambium formation and growth of annual rings are June and July [160,161]. Similarly, we concluded that the months of April and July were the most significant in our study.
In the radial growth, 12-year cycles were observed in Norway spruce. These 12-year periods in radial growth may be caused by solar activity that repeats at 8 to 12-year intervals [162]. These cyclical processes have been recorded in the eastern Krkonoše Mountains on European beech stands [7], and have also been confirmed for beech increment in southern Italy [163]. Even Scots pine (Pinus sylvestris L.) in Russia has recorded these 11-year cycles in the annual growth [164,165].
Moreover, four to seven-year cycles were recorded for spruce, beech, and mountain elm. We classify these cycles as the "Schwabe cycle", which is the second harmonic cycle of the sun [166]. These four to seven-year periods are associated with shorter-term temperature and precipitation cycles [163], where, for example, eight-year cycles have been recorded in Germany as most synchronous with vegetation phenology of beech [167].
Longer 20-year cycles were classified as the 22-year Hale cycle, which repeats through meteorological indicators such as temperature and ozone concentration [168] or precipitation [169]. This 22-year cycle has also been found in annual rings in eastern Russia [170] or Slovakia [171]. Longer 60-year cycles in the spectral analysis of spruce and beech are among the influences of the Gleissberg cycle, which is repeated from 70 to 100 years [168], and this cycle is associated with global temperature changes [172].

Conclusions
The study of radial growth in mixed stands with dominant European beech, Norway spruce, and admixed sycamore maple, European ash, and mountain elm on very stony shallow soils in the Krkonoše Mts. shows that individual tree species respond appropriately to ongoing climate change (increase of temperature and fluctuation in precipitation). While beech and spruce predominantly growing in more favorable habitat conditions (deeper soil profile and less skeletal soils) have increased in the last two decades, the radial growth of sycamore and ash growing in the most extreme soil conditions (very shallow and strongly skeletal soils) decreased. This indicates that a more severe impaction of climate change on forest stands in extreme habitats can be expected. Moreover, the occurrence of negative pointer years with extremely low radial growth and growth variability has increased in the last 20 years. On the other hand, the most vulnerable tree species to growth cyclicity, climate factors, and air pollution load were spruce and beech compared to the resiliency in ash and especially in sycamore. These durable tree species should be supported by silviculture practice on these endangered, specific habitats. Their significant stabilizing potential in the most extreme soil conditions of shallow grounded rubble during global climate change is very important. In these habitats and stand conditions, it is essential to fulfill the ecological functions of the forest, especially the soil, climate, and water protection functions, which are becoming more important due to climate change. These mixed scree forests should be systematically maintained by adaptive forest management. These measures should include the supporting of vital deep-rooted deciduous trees (beech, maple, and ash) during sanitary silviculture interventions and the introduction of other native tree species (silver fir, rowan, limes) to increase stand diversity, resilience and sustainability.

Data Availability Statement:
Monthly climate data and air pollution data for Krkonoše Mountains are available from the Czech Hydrometeorological Institute, Prague (www.chmi.cz). The tree ring data and stand data presented in this study are available on request from the corresponding author.