Next Article in Journal
USDA Forest Service Forest Products Laboratory: Acetylation of Wood 1945–1966
Next Article in Special Issue
Seasonal Photosynthesis and Carbon Assimilation of Dynamics in a Zelkova serrata (Thunb.) Makino Plantation
Previous Article in Journal
A Novel Tree Biomass Estimation Model Applying the Pipe Model Theory and Adaptable to UAV-Derived Canopy Height Models
Previous Article in Special Issue
Prospects of Using Tree-Ring Earlywood and Latewood Width for Reconstruction of Crops Yield on Example of South Siberia
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 12
Issue 3
10.3390/f12030259
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tree Rings of European Beech (Fagus sylvatica L.) Indicate the Relationship with Solar Cycles during Climate Change in Central and Southern Europe

by
Václav Šimůnek
1,*,
Zdeněk Vacek
1,
Stanislav Vacek
1,
Francesco Ripullone
2,
Vojtěch Hájek
1 and
Giuseppe D’Andrea
1
1
Faculty of Forestry and Wood Sciences, Czech University of Life Sciences Prague, Kamýcká 129, 16500 Praha 6-Suchdol, Czech Republic
2
School of Agricultural, Forestry and Environmental Sciences, University of Basilicata, Viale dell’Ateneo Lucano 10, I-85100 Potenza, Italy
*
Author to whom correspondence should be addressed.
Forests 2021, 12(3), 259; https://doi.org/10.3390/f12030259
Submission received: 14 January 2021 / Revised: 12 February 2021 / Accepted: 19 February 2021 / Published: 24 February 2021
(This article belongs to the Special Issue Effect of Climate Change on Forest Growth and Phenology)
Browse Figures
Versions Notes

Abstract

:
The impact of solar cycles on forest stands, while important in the development of the forest environment during climate change, has not yet been sufficiently researched. This work evaluates the radial growth of European beech (Fagus sylvatica L.) in the mountain areas of southern Italy and central Europe (Czech Republic, Poland) in correlation to solar cycles (sunspot number), extreme climatic events, air temperatures and precipitation totals. This research is focused on the evaluation of the radial growth of beech (140 dendrochronological samples with 90–247 years of age) from 1900 to 2019. The time span was divided into the following three periods: (1) a period of regular harvesting (1900–1969), (2) a period of air pollution crisis (1970–1985) and (3) a period of forest protection (1986–2019). The results indicate that the solar cycle was significantly involved in radial growth on all research plots. With regard to the evaluated precipitation totals, seasonal temperatures and the sunspot number, the latter was the most significant. Temperatures had a positive effect and precipitation had a negative effect on the radial increment of beech in central Europe, while in southern Italy, the effect of temperature and precipitation on the increment is reversed. In general, the limiting factor for beech growth is the lack of precipitation during the vegetation season. The number of negative pointer years (NPY) with an extremely low increment rose in relation to the decreasing southward latitude and the increasing influence of climate change over time, while a higher number of NPY was found in nutrient-richer habitats compared to nutrient-poorer ones. Precipitation and temperature were also reflected in the cyclical radial growth of European beech. The relationship between solar cycles and the tree ring increment was reversed in southern Italy and central Europe in the second and third (1970–2019) time periods. In the first time period (1900–1969), there was a positive relationship of the increment to solar cycles on all research plots. In the tree rings of European beech from southern Italy and central Europe, a relationship to the 11-year solar cycle has been documented. This study will attempt to describe the differences in beech growth within Europe, and also to educate forest managers about the relevant influence of solar cycles. Solar activity can play an important role in the growth of European beech in central and southern Europe, especially during the recent years of global climate change.

1. Introduction

European beech (Fagus sylvatica L.) is distributed over a large number of environmental zones throughout Europe [1]. A lack of precipitation, associated with drought, limits European beech in its distribution in southern Europe [2,3,4], while low temperatures and spring frosts curb the beech in the north [5,6]. European beech thrives in the optimal conditions of central Europe, and at the same time, we are witnessing its retreat from the southern parts of Europe [7]. Lately, due to climate change, this tree species has also been expanding into more northern areas within its distribution range [8,9,10]. The occurrence of beech forests is limited by late frosts and droughts [11]. In the past few decades, European beech has been showing better results in locations towards the north, such as southern Sweden [12], as well as in mountainous locations, e.g., in the Czech Republic [13,14,15]. In recent years, however, beech forests have been showing a greater sensitivity to drought [16], yet this tree species can adapt accommodatingly to drought in central Europe under the conditions of global climate change [6,13,17].
Concerning radial growth, European beech adapts better to climatic fluctuations at lower altitudes and in nutrient-rich habitats [18,19]. By contrast, at higher altitudes, there have been greater fluctuations in beech growth since 1975 [20]. Climatic fluctuations in the mountainous areas of central Europe are most frequently caused by lower temperatures, with annual temperatures affecting the radial growth of European beech more than precipitation totals [13,21,22]. Fluctuations in the radial growth of European beech, caused by various biotic and abiotic influences, can lead to cyclic periods [21]. For example, in the 1970s and 1980s, there was a significant negative effect of air pollution on the radial growth of European beech in the mountains of the Czech Republic, which also led to a weakening or even dying of forest stands [23,24]. Notable incremental fluctuations of European beech are also induced by frost, which can affect the increment of the entire vegetation period [21,25].
European beech stands in the Czech Republic and Italy differ in their life cycles as well as their life span. In the Krkonoše Mountains (both in the Czech Republic and Poland), generations of European beech stands are replaced after 230 to 250 years on average [26]. In contrast, old stands of European beech in Italy may live to be around 300 years or more [27,28]. The age of beech stands in Italy will shorten in response to global climate change, especially increasing average temperatures and droughts [28]. However, the radial growth of beech stands in Italy is more affected by the negative effects of drought, which limits the distribution of this tree species in the Apennines [29]. Beech stands in Italy are profoundly more vulnerable and frequently stressed by drought, and this can reduce beech increment over time [27]. In the Mediterranean forests of higher altitudes, there are also spring droughts, which negatively affect the tree rings of European beech [30].
Solar cycles affect the climate of our planet, as was cited by the theory of “solar forcing”, to name one. The theory documents that solar irradiance and geomagnetic activities are important drivers of the Earth’s atmosphere [31]. Some works also describe that solar activity, in the long run, affects the temperature of the Earth’s atmosphere, which can be accompanied by the occurrence of colder winter periods [32]. Solar activity is sometimes a stronger inducer of climate change than anthropogenic influences, also affecting other factors such as precipitation or temperature [33]. Sunspot cycles are also associated with the cycle of the Amazon River flow [34] and are directly related to flood or drought seasons in the Songhua river basin in China [35]. On the other hand, the solar cycles may be associated with drought periods and extreme weather fluctuations, for example in Kuwait [36]. The effects of 11-year solar cycles on climate patterns in Europe have been observed over the last 250 years, while the impact of sunspot activity on the climate increased during the 19th century [37].
Global climate is affected by solar cycles [38,39]; however, the imprint of this cycle might not be the same everywhere. The imprint of solar activity is different in the Indian, Pacific and Atlantic Oceans [40]. There is a large number of factors that react to solar variability differently within terrestrial climate parameters [41]. One important parameter is that solar activity affects the properties and formation of clouds in our atmosphere, which is caused by the ionization of cosmic rays that mirrors 11-year cycles of solar activity [42,43,44]. Some studies even found a link between tree ring radiocarbon production and sunspots [45,46]. The production rate of radiocarbon isotope is an indicator of cosmic radiation in the upper atmosphere, which is also connected to the solar cycle [39]. Solar cycles also affect the average monthly temperatures in Europe [47], as well as, indirectly, the North Atlantic Oscillation (NAO) and the Atlantic Multidecadal Oscillation (AMO) [48]. Solar cycles throughout the NAO are also partly associated with precipitation occurrence in Europe [49]. Solar activity may also be reflected in many other climatic cycles [50]. Tree rings of beech in the Apennines are associated with the NAO [51], and even the NAO is influenced by the solar cycle [37,49]. It is also important to mention that there is a link between fluctuations in the radial growth of European beech in the mountains of the northeast Czech Republic and solar activity [52,53]. Some research in the field of dendrochronology proves the influence of solar cycles on the radial growth of trees, which is well described by research from northwestern Russia [54,55], Chile [56] and the Tibetan Plateau [57]. The impact of solar cycles can also be negatively reflected in tree rings, which has also been investigated in cork oak (Quercus suber L.) [58].
This study focuses on the radial growth of European beech at higher altitudes of the Krkonoše Mountains in Poland and the Czech Republic and in the southern Apennines in Italy. The main objective of this study was to evaluate the initial effect of the temperature, precipitation and sunspot number on the radial growth of European beech during three time periods that cover 119 years of growth history. The time periods were chosen deliberately due to various and important differences in the development of beech stands over time. The first time period (1st period) indicates the phase of man-managed stands of European beech. The second time period (2nd period) indicates the air pollution calamity on research plots in the Czech Republic [23,52], and the same time period was also marked out for Italian research plots in order to maintain comparability of time periods. The third time period (3rd period) denotes the phase without human harvesting interventions, when the beech stands developed naturally and spontaneously. The sunspot number was chosen intentionally to find a possible relationship between the radial growth of European beech and the 11-year solar cycle, which has properties and intensities other than precipitation or temperature. Another aim of this research was to assess the cycles of radial growth of European beech using the Fourier analysis and cross-correlation in relation to the factors examined. This is a very underexplored topical issue regarding the growth process of one of the most important and promising central European tree species during global climate change.

2. Methodology

2.1. Study Area

The studied area is located in the national parks of the Krkonoše Mountains, both in the Czech Republic and in Poland, i.e., in the Krkonoše/Karkonosze Transboundary Biosphere Reserve. Italian research plots are located in the southern Apennines within the Appennino Lucano National Park, near a place called Sellata. A total of 4 permanent research plots were evaluated, of which the first is in the Polish portion of the Krkonoše Mountains in a place called Chojnik, not far from the town of Sobieszów. The second research plot is in Eastern Krkonoše, in Rýchory, near the town of Žacléř. The third and fourth research plots are in the southern Apennines near Sellata. The plots in Krkonoše were established and selected in 1980 for long-term research. All research plots were selected in accordance with the major representation of beech, altitude and homogeneous stand structure of European beech. Samples were taken from forest stands with a 100% share of European beech. All research plots are located in non-intervention areas, where no harvesting operations had been carried out since 1985.
The basic site and stand characteristics are given in Table 1. The worldwide Köppen classification was used for climate categorization in Table 1 [59]. Precipitation and temperature conditions differ for the Czech Republic/Poland and for Italy, as described in Figure 1. The Bedřichov meteorological station for Krkonoše (Chojnik 1 and Rýchory 2) was selected for its data availability, as it has the longest time series for precipitation and temperatures available to the Krkonoše montane region at the altitude of 780 m a.s.l. The mean annual temperature for Chojnik 1 and Rýchory 2 is 3.7 °C, and the annual total precipitation is 1300 mm. The mean number of vegetation days for Czech and Polish plots is around 120 days and the number of days with snow cover reaches 117 [21]. The mean annual temperature for the Italian plots is 13.7 °C and the mean annual precipitation is 1032 mm according to the Abriola meteorological station (1050 m a.s.l.). The mean number of vegetation days for Italian plots is 135 days (same period for snow cover) [60]. In Figure 1, a descriptive map of research plots is also displayed.

2.2. Data Collection

For dendrochronological analysis of samples from the Czech Republic, Poland and Italy, increment cores were taken from European beech using a Pressler auger, perpendicular to the trunk axis at a height of 1.3 m above the ground. Structurally homogeneous beech stands with stocking of 0.8–1 were selected for sampling. Samples were taken from randomly selected (RNG function) healthy co-dominant and dominant trees, whose average tree diameter at breast height had to be dbh > 30 cm. The height of all sampled trees was measured with a Laser Vertex hypsometer (Haglöf, Långsele, Västernorrland, Sweden). The Vertex hypsometer shows the measurement accuracy (instrumental resolution) to 0.1 m according to the manufacturer. The diameter at breast height was also measured for all trees using a Mantax Blue metal caliper (Haglöf, Långsele, Västernorrland, Sweden). This caliper provides an accuracy (instrumental resolution) of 1 mm according to the manufacturer. These research plots were monitored annually, and their detailed structural analysis was performed every 5 years. A total of 140 samples were taken for dendrochronological analysis; Table 2 offers more detailed information. Increment cores were measured using a LINTAB measuring table [61] with an Olympus microscope. The measuring table provides an accuracy (instrumental resolution) of 0.01 mm, and the TSAP-Win software [61] was used to record the cores. Measurements were made from the bark towards the heartwood, perpendicular to the center of the trunk, so that each tree ring was measured perpendicular to the trunk axis. The subsequent cross-dating of the increment cores was performed with the Cdendro program so that the cross-correlation index was CC > 25 for each sample [62].
Monthly temperature and precipitation data for Krkonoše were provided by the Czech Hydrometeorological Institute, Prague [63]. The mean monthly climatic data for the research plots in the Czech Republic and Poland were supplied by the meteorological station of Bedřichov (50°47′30.7″ N 15°08′31.7″ E) at an altitude of 780 m a.s.l.; the station is 35.5/54.3 km away from the research plots in the Czech Republic and Poland, respectively. Monthly temperature and precipitation data for the research plots in Italy were provided by the Italian Civil Protection Authority, Basilicata Region [64]. Characteristics of the Italian research plots were provided by the meteorological station in Abriola (40°30′28.0″ N 15°48′46.6″ E), at an altitude of 1050 m a.s.l. The distance from the Sellata 3 and La Lama 4 research plots to the meteorological station is 3.2/5.1 km, respectively. Annual sunspot number data were taken from the Royal Observatory of Belgium, Brussels [65].

2.3. Data Analysis

Dendrochronological data were processed in R software [66] using the “dplr” package [67]. Detrending of each tree was performed by negative exponential detrending with an inserted spline of 1/3 of the age of each tree using “dplr” instructions [68]. Such detrending removes the age trend while maintaining low-frequency climate signals [55,69]. The expressed population signal (EPS) was calculated for the 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 in relation to the climatic data was a significant EPS threshold so that EPS > 0.85 [68]. We also calculated the signal-to-noise ratio (SNR) that represents the signal strength of the chronology and R-bar (inter-series correlations) [70]. The analysis of negative pointer years was carried out [71]. For each tree, the pointer year was tested as an extremely narrow tree ring that did not reach 40% of the increment average from the four preceding years [71]. The occurrence of the negative year was proved if a strong reduction in increment occurred in at least 20% of the trees in the plot.
Spectral analyses for indexed (detrended) radial increments were created with Statistica 13 software [72]. The calculation was performed with the “Single Fourier (Spectral) Analysis” function, using the output “Periodogram” plot by “Period”. Furthermore, this software was used to calculate correlation coefficients and cross-correlations for lag −15 years (in relation to the sunspot number, vegetation season temperature and annual precipitation). With Statistica 13 software (Statsoft, Tulsa), we created cross-periodograms (real), used to study a multivariate spatial process [73]. We also used the “signal” and “dplr” wavelet plot for Krkonoše and Italy for the period 1900–2019 to make the dendrochronological signal more apparent for up to 32 periods/cycles [74,75]. For the next statistical analysis, the time span was divided into the following 3 periods: (1) a period of active forest management and regular harvesting (1900–1969), (2) a period of air pollution crisis characterized by extremely high SO2 concentrations and acid rains (1970–1985), and (3) a period of forest protection without human harvesting interventions (1986–2019).

3. Results

3.1. Tree-Ring Characteristics and Extreme Climatic Events

Dendrochronological characteristics are described in Table 2 and include the mean tree ring increment at a plot (Mean), age of the youngest and the oldest sample trees (Age), standard deviation (Std.), inter-series correlations (R-bar), expressed population signal (EPS) and signal-to-noise ratio (SNR). The research plots in Italy (Sellata 3 and La Lama 4) show a larger mean increment, on average 25% larger than the research plots in central Europe (Chojnik 1 and Rýchory 2). The age of the research plots indicates that Sellata 3 and La Lama 4 were generally older than Chojnik 1 and Rýchory 2. Even the research plot of La Lama 4 reached up to 247 years. The range between maximum and minimum age was due to age variability in the forest stands; however, all research areas showed a visually homogeneous forest structure. A greater variability of radial growth was found in both locations with lower altitude compared to a higher one, while a higher number of negative pointer years (NPY) characterized by an extremely low increment was found on richer soils compared to poorer ones.
The number of NPY ranged from 0–4 for central Europe and from 6–8 in southern Italy. A climatically significant year common for Rýchory 2 and Sellata 3 was 2016, when the share of precipitation, lower by 26% (compared to the average), was documented in the vegetation season. Another common NPY was recorded in 2013 for the Italian research plots Sellata 3 and La Lama 4, where a higher annual precipitation of 1204 mm (compared to 1024 mm per year on average) and lower seasonal temperatures of 28.3 °C (29.6 °C) were recorded in both plots. One more common NPY for the research areas in Italy was 1970, with the seasonal precipitation totals lower by 15% (240 mm, 280 mm). The year 2017 was also significant for the Italian research plots Sellata 3 and La Lama 4, with the total precipitation during the vegetation season lower by 41% (204 mm, compared to an average of 289 mm), the lowest in 20 years. Similarly, the NPY of 2011 in Krkonoše was negatively affected by an unbalanced frequency of precipitation during the year (monthly variability was higher by 46%). Overall, the NPY show that the Apennines plots (Sellata 3, La Lama 4) have a higher frequency of NPY than central Europe, while the Chojnik 1 research plot, for example, has not recorded a single NPY. Additionally, 15 of the 18 NPY were found in the second half of the solar cycle.

3.2. European Beech Tree-Ring Growth, Sunspot Number and Climate Variations

Figure 2 shows a clear difference between radial growth (a) in the Krkonoše Mountains in central Europe (Poland, Czech Republic) and (b) in the southern Apennines in Italy. The radial growth of beech differentiates itself in each time period. A positive relationship of radial growth in all research plots with the sunspot number is characteristic for the first time period. This fact is also confirmed in Table 3. A change in the trend of radial growth of European beech and the nature of its connection with sunspots occurred in the second time period on the Italian research plots. Another interesting feature of the second time period is a unique low or negative correlation with sunspots on almost all research plots, which is also confirmed in Figure 2 and Table 3. In the second period, there was also a significant decrease in radial growth on all research plots. In the last, or third, time period, the concurrence of sunspots and radial growth of beech in Italy was the opposite of that on research plots in the Czech Republic and Poland. A typical interconnection with the solar cycle (sunspot number) and the radial growth of European beech is shown in Figure 2. These are, for example, parallel cycles of radial growth with the solar cycle in the first time period, but the solar cycle is also reflected in the third time period (positively and negatively). In Figure 2, a significant reduction in radial growth can also be seen, often during the solar minimum, which can be reported for almost all time periods, except for the Italian research plots, where in the second and third periods, this process was reversed.
Table 3 describes the relationship between the radial growth of European beech in the research plots and the sunspot number, annual temperature, vegetation season temperature, annual precipitation and vegetation season precipitation in different time periods. Radial growth of beech reacts to sunspots, temperatures or precipitation differently in every time period. The most significant values were correlated with the sunspot number, then with annual precipitation, and subsequently with seasonal temperatures. Mean annual temperatures and seasonal precipitation totals did not show high correlation values to radial growth, but they also did not show significant results. According to Table 3, solar cycles correlated most significantly with the whole period, and subsequently, with the third time period. The tree ring width index (RWI) showed the most significant positive correlation with solar cycles in Poland and the Czech Republic (RWI Chojnik 1 and Rýchory 2) both for the whole period and the third period. For the Italian plots Sellata 3 and La Lama 4, the correlation results were different, with both plots having negative correlations with the sunspot number. The Sellata 3 plot showed a significant negative correlation over the whole period (r = −0.25) and in the third period (r = −0.37). La Lama 4 had a significant correlation in the second time period (r = −0.51), so all plots significantly correlated with the sunspot number, depending on the time period. All plots were positively correlated with solar cycles in the first time period, but only the RWI Chojnik 1 plot was significantly positive (r = 0.31) in the first period.
Seasonal temperatures had only one significant correlation with radial growth for the La Lama 4 plot in the third time period (r = −0.39). In the case of seasonal temperatures, positive correlations were found for the Czech and Polish research plots, while for the Italian ones, the values were mostly negative.
The annual precipitation totals had two significant values: for the Rýchory 2 plot for the whole time period (r = −0.18) and for the La Lama 4 plot in the first time period. Annual precipitation totals correlated differently to radial increment than seasonal temperatures or sunspot number. Thus, annual precipitation totals correlated negatively with radial growth in the Czech Republic and Poland in practically all the time periods. The Italian research plots Sellata 3 and La Lama 4 correlated positively with the annual total precipitation in all the time periods.

3.3. Cross-Correlation up to Fifteen Years Back in Relation to Tree Ring Growth

Cross-correlations in Figure 3 show the relationship of the sunspot number, annual precipitation and seasonal temperatures to the radial growth from up to 15 preceding seasons in different time periods. The cross-correlations of this research describe the relationship of radial growth to the studied factors up to 15 years into the past. Cross-correlations show that the most significant values are found in the sunspot number, while the number of significant values in precipitation and temperature are almost identical. In relation to the radial growth of European beech, the sunspot number correlates in waves. The most concurrent correlations with sunspots on all research plots are found in the first time period. By contrast, in the second and third time periods, the correlations turn completely against each other when the Chojnik 1 and Rýchory 2 plots are positively correlated to the sunspot number, while the Sellata 3 and La Lama 4 plots are negatively correlated to the solar cycle. The most common and highest significant correlation to sunspots in terms of relative lag years is in lag 0 to −1 years for Rýchory 2 and Chojnik 1 in both the whole and first time periods. Nevertheless, the La Lama 4 plot is significant to sunspots in lag −7 to −8 relative years in the third time period. In the second time period, the plots of Sellata 3 and La Lama 4 correlated significantly to the sunspots in lag 0 to −1 relative years. The cross-correlations of radial growth to sunspots show that all time periods contained significant values for almost all research plots.
Annual precipitation totals and seasonal temperatures cross-correlate less than the sunspot number, which is apparent in Figure 3. In the case of precipitation and temperature, it is also true that the research plots correlate conversely between the plots in Krkonoše and in Italy. For precipitation, the whole time period and the first time period correlate best, which are also the longest monitored time periods. Significant values for precipitation are found in the whole time period for the relative years 0 and −10. Another important period for precipitation is the first time period, where the values significantly correlate back from lag −6 to −10 relative years. The cross-correlation of radial beech growth and seasonal temperatures only shows significant values in the first time period. The results for precipitation in the first period show a significant value for the research plot Rýchory 2 (lag −2 relative years), La Lama 4 (lag −2 relative years) and Sellata 3 (lag −4 to −5 relative years and also −12 years).

3.4. Spectral Analysis

Cross-periodograms show common cycles between the radial growth data series in relation to the sunspot number, annual precipitation and seasonal temperatures. Cross-periodograms describe both a negative and positive cyclical relationship between our examined data. Figure 4 describes the whole time period, with an assessment of the influence of “SUNSPOT” (sunspot number), “Total prec.” (annual precipitation) and “Season temp.” (season temperature) on radial growth from all research plots in the period 1900–2019.
Solar cycles (sunspot number) are reflected in the observed positive radial growth of European beech by 11-year cycles on all research plots. The La Lama 4 research plot, where the 11-year cycle correlated both positively and negatively with the indexed radial growth and the sunspot number, shows 20- and 40-year cycles. On the Chojnik 1 research plot, there was a negative cross-period of radial growth and sunspots in 20- and 40-year cycles. The Rýchory 2 research plot formed negative cross-periods of radial growth and sunspots in 6- and 30-year cycles.
Important cycles of the cross-periodogram for all research plots in relation to precipitation and increment show that on the Italian plots (Sellata 3 and La Lama 4), there were negative values for 30- and 40-year cycles. On the research plots in Krkonoše (Chojnik 1, Rýchory 2), there were various periodic responses of radial growth of European beech to the total precipitation, but it is possible to say that the 30- to 40-year periods were not as significant as in Italy. The cross-periodogram shows that precipitation cycles were negatively correlated with radial beech growth on the Chojnik 1 and Rýchory 2 research plots. On the Sellata 3 and La Lama 4 plots, positive 60-year precipitation cycles in radial growth have occurred.
The seasonal temperature cross-periodogram shows minor differences. The Chojnik 1 and Rýchory 2 plots had the most legible 25-year positive cycles of seasonal temperatures with radial growth in the periodogram. The cross-periodogram also shows that the research plots in Italy (Sellata 3 and La Lama 4) have had 40-year positive cycles in relation to seasonal temperatures and radial growth. Furthermore, the La Lama 4 plot has had negative 60-year cycles of increment as related to temperature, on the cross-periodogram.
The Fourier analysis in Figure 5 shows the difference between the 1st and third time periods. As in the graph in Figure 5, the Chojnik 1 research plot showed 4-, 7-, 11- and 18-year cycles in the first time period. In the third time period, there were 7- to 12-year cycles on the Chojnik 1 research plot. In the first period, the Rýchory 2 research plots showed 4-, 18- and 35-year cycles in radial growth. In the third period, there were 3-, 33- and 7- to 11-year cycles in radial growth on the research plot of Rýchory 2. The Sellata 3 research plot showed 8-, 12- and 24-year cycles in the first period. In the third time period on the Sellata 3 research plot, there were 4- and 12-year cycles in radial increment. The La Lama 4 research plot in the first time period showed 5-, 10-, and 18- to 35-year cycles in relation to the radial growth of beech. In the third period on the La Lama 4 research plot, there were 4-, 5- and 7-year cycles in the radial growth of European beech.
It is evident from Figure 5 that 7- to 12-year cycles most commonly occur in radial growth, but these results are time-period-dependent. Furthermore, each research plot apparently has 9- to 12-year cycles in at least one time period.
The wavelet plot in Figure 6 describes the occurrence of cycles (periods in the graph) in time between the Apennines and the Krkonoše Mountains. It documents that 24- to 32-year cycles during the period 1900–1990 were the most important in Krkonoše. There were also less significant cycles in the period (dark gray color) from 1910 to 1960 in Krkonoše, with approximately 6- and 20-year cycles reported. However, the darkest gray color on the Krkonoše research plots shows that from 1985 to 2018, cycles of around 4 and 12 years were also important. Research plots in Italy revealed significant 16- to 22-year cycles for the period 1900–1940, and also 4-year cycles for the period 2010–2015. Periods from 1920 to 1940 with cycles of about 12 years were less significant on the Italian plots. In addition, there were less significant 3- to 8- year cycles in the period of 1960–2015. The wavelet plot in Figure 6 shows 9- to 13- year cycles, which are the most common for Krkonoše and Italy in the 1990–2019 period, but the results from Figure 6 are below the statistical significance limit.

4. Discussion

4.1. Different Growth Conditions of European Beech in Central Europe and the Southern Apennines

Beech stands in central Europe have fewer vegetation days than in the Mediterranean, which is reflected in the size of the mean radial growth. Increment on the research plots in Krkonoše is smaller by an order of magnitude (approximately by 25%) than on the research plots in the southern Apennines (Table 2). Beech stands in Krkonoše belong to a European beech provenance more sensitive to drought, which may be related to its earlier fall of the assimilation apparatus. In contrast, beech stands in the Mediterranean are more drought-resistant and defoliate later [28,76]. In the Krkonoše Mountains, temperatures were proven to have a greater effect on the radial growth of European beech than precipitation totals [21,77]. On the Chojnik 1 research plot, our results in Table 3 confirm that temperatures had a higher correlation with radial growth than with the total precipitation, while on the Rýchory 2 research plot, this fact is validated in the third time period. The total precipitation had a negative effect on the radial increment in central Europe. Similar negative responses to precipitation were also confirmed by reference sources [22,78]. For beech in Krkonoše, it is also true that with a higher altitude, the temperatures play a greater role in radial increment [79], which is verified by our results. It is evident from Table 3 that in the 1st and third time periods, the temperatures correlated more at Rýchory (760 m a.s.l.), which is located higher than the research plot of Chojnik 1 (510 m a.s.l.). Altitude also influenced growth consistency at both locations, while higher incremental variability was found in lower-situated areas compared to those with higher altitude. In general, beech stands have faced increasing drought stress in recent years [11], but in the montane areas of Krkonoše, spring frosts have more significant negative effects [21].
Beech forests in the Mediterranean are also affected by drought and spring frosts, as evidenced, for example, in the Pyrenees [80] and the Apennines [81]. Drought is a limiting factor in the European beech growth in the Apennines [29], as is confirmed by the negatively significant years (NPY) in Italy, where the impact of drought on radial growth was recorded on our plots Sellata 3 and La Lama 4 several times, e.g., in 1970 and 2017. In general, more NPY were detected in central and southern Europe in nutrient-richer habitats compared to the poor ones. On productive sites, climate change can hamper tree growth and forest productivity [82]. Also, other studies have shown that trees were more sensitive to climate on the more productive sites compared to poor ones due to differing proportions of latewood to earlywood [83,84]. Assessing tree responses to climate change without simultaneously considering soil properties and climate may be misleading, since soil nutrients can influence the growth response of trees to drought [85]. Our results also showed that 15 of the 18 NPY were observed in the second half of the second part of the studied period. An increasing number of NPY during the time may be caused by an increasing number of extreme climatic events in ongoing climate change [80,86,87].
Beech stands in the southern Apennines (Sellata 3 and La Lama 4) were negatively correlated with temperatures and positively correlated with precipitation (Table 3 and Figure 3), which is exactly the opposite from the research plots in Krkonoše. Scientific literature also confirms that these correlations of precipitation and temperature to radial beech growth in the Apennines may be due to winter accumulation of precipitation in the soil [51]. The negative correlation of the beech radial growth to the temperatures on the Sellata 3 and La Lama 4 plots could be explained by the effect of higher temperatures during dry periods [88].

4.2. Solar Cycles, Climate Change and the Possible Link to European Beech

Solar activity affects the climate on our planet [50], which is also reflected in the NAO [37,49]. Many other studies have also shown that the effect of solar activity has an impact on precipitation and temperature [33,34]. Temperatures on the European continent correlate with the NAO mainly during the winter and spring periods [47]. Precipitation totals are also linked to solar activity, but it is important to mention that precipitation is also affected by solar activity during July in Italy and during May, June and July in the Czech Republic [49]. As in Figure 1, there is also a noticeable difference between monthly precipitation totals and temperatures during the year, when precipitation totals differ significantly—in Krkonoše, the main precipitation occurs over the vegetation season, while in Italy, it occurs during the vegetation off-season.
Wind currents high in the atmosphere, or “jet streaming”, are associated with solar cycles, where these wind currents are blocked during the solar minimum [89,90], leading to a colder winter season [91]. In contrast, the solar maximum leads to an acceleration of wind currents [92], which, through “Ferrel cells” [93], increases the pressure of winds affecting the European continent. This may be confirmed by research in Spain, where the solar cycle has been found to be associated with precipitation and wind anomalies [94]. The solar cycle is associated with the occurrence of large forest fires and dry weather seasons during the solar minimum, which has been proven in Turkey [95]. Our research plots in the mountainous areas of the Mediterranean, Sellata 3 and La Lama 4 (Figure 2), have recorded higher increments in recent years (third period) during solar minimums. This fact can be supported by better conditions for growth at higher altitudes, which is confirmed by correlations presented in Table 3, showing the relationship between temperatures and radial growth as being predominantly negative.
Solar cycles are linked to climate change and temperature differences on the planet’s surface. There is a study that describes a possible association of low sunspot activity with climatically cold periods on the Earth [32]. Solar cycles do not affect temperatures and precipitation directly, and it has been shown that solar activity is not imprinted on these factors entirely [86,96]. The total contribution of solar activity to variations in the Earth’s global temperature and climate is insignificant but not negligible [97]. Other studies describe the link between cosmic rays and solar cycles [98]. It is the cosmic rays that are associated with cloud formation, and this also affects the amount of light falling on the planet’s surface, where low solar activity (low sunspots) leads to a higher amount of cosmic rays, creating a thicker cloud cover and cooling the planet [38,39]. Our results may be related to the decreasing activity of sunspots [99] due to the fact that there has been a lower radial growth in Krkonoše since around 1990 (Figure 2a), and a decrease in the radial growth of beech in the southern Apennines since 2010 (Figure 2b). In the central Apennines, there has also been a long-term reduction in the radial growth of European beech due to drought since 1970 [27], while concurrently, sunspot number has been lower since 1980 (Figure 2). Our results show the response of radial growth to solar cycles in Figure 2 and Figure 3 and Table 2, when the research areas in Krkonoše correlated positively and the areas in southern Italy correlated negatively in the third time period. The opposite reactions of the radial growth to the solar cycle were observed in our results after 1960. Simultaneously, global temperature was affected due to increasing CO2 since 1960, which disturbed the natural process of solar cycles [100]. This fact of opposite correlations could be attributed to solar activity or climate change, where both factors are accompanied by changes in the NAO [47]. This is because changes in air temperatures on the surface of the Earth can also be related to the changes in low and high pressures above Europe [91,101].

4.3. Feedback and the Relationship of Radial Growth to the Sunspot Number Across the Time Frame

Solar cycles play a role in the radial growth of European beech in Krkonoše, which has been proven in the vicinity of research plots in Rýchory [52]. Our results confirm that not only Rýchory 2, but also the northern part of Krkonoše (Chojnik 1), are associated by positive correlations (Table 3) with the solar cycle during the first and third time periods (Figure 2, Figure 3, Figure 4, Figure 5). A positive relationship of Scots pine growth has also been found in northwestern Russia [55]. It was also confirmed in Krkonoše in the 1980s, when an air pollution calamity (high SO2 concentrations) occurred, reducing the radial growth of European beech [21,24,102]. Our results from the research plots in Chojnik 1 and Rýchory 2 in the second time period also describe this event. In the second period, there was also a decrease in radial growth on the Italian research plots of Sellata 3 and La Lama 4, while Figure 3 shows a significant result of the negative cross-correlation with solar activity for both Italian plots. In the third time period, the negative relationship between the plots in the Apennines (Sellata 3 and La Lama 4) also continued, while the first and second periods could be compared to the research of solar cycles and their relation to cork oaks in the Mediterranean in Portugal, where a negative relation between tree bark increment and the solar cycle was confirmed [58]. Other studies have also confirmed a negative correlation of solar cycles with the radial growth of Pinus pinaster in northern Portugal. The same study found a minor negative correlation of radial growth and solar cycles in southern Slovakia [103].
Cross-correlations in Figure 3 show the feedback of the radial growth of European beech and solar cycles on the Italian research plots in the third time period, where Sellata 3 showed significant correlations in the 7 to 8 preceding years. This could indicate a possible shift in the effect of solar activity in the Apennines. We can even find regression correlations of increment and solar cycles in the first and second time periods for the Italian plots, which are 1- to 2-year shifts in the significant correlation coefficients. The effect of cross-correlations with solar cycles could be attributed to the effect of the NAO, while for European beech, the influence of the NAO on the radial growth in central Italy was confirmed [51]. Changes in the radial growth of European beech in Italy in the second and third time periods may have been due to the influences of the NAO over Europe, which may have caused different responses of cross-correlations between Krkonoše and the Apennines [104,105].

4.4. Recorded Cycles in Beech Tree Rings

In terms of radial growth, our results in Figure 5 and Figure 6 show 9- to 12-year cycles most frequently, recorded for each research plot in at least one of the time periods. These periods are associated with solar activity, which has 8- to 12-year cycles. These cycles have been recorded in the radial growth of various tree species, e.g., in Europe [103], Russia [55], South America [106] and Asia [57]. Furthermore, 4- to 7-year cycles have been reported for radial growth, which can be classified as the “Schwabe cycle”, the second harmonic cycle of the sun associated with the aforementioned cosmic rays [107]. The high-frequency oscillations for about three years might be associated with the frequent changing in global temperatures [108,109]. In addition, Figure 4 shows a minor share of 4- to 7-year cycles in a cross-periodogram with temperatures and precipitation, where, e.g., 8-year cycles have been recorded in Germany as being most synchronous with the vegetation phenology of beech [110]. Short (3- to 4-year) cycles can also be found in the cross-periodogram of radial growth and precipitation, and again, they may be linked to the NAO [49]. Longer (17- to 35-year) cycles could be assigned to the 22-year Hale cycle, which repeats through meteorological indicators such as temperature and ozone concentration [40] or precipitation [49]. This 22-year cycle has even been found in tree rings in Tibet [57] and northwestern Russia [55]. Our results also partly include the occurrence of 33- to 36-year cycles, which might indicate the dynamics of droughts [108]. Multi-year cycles in the cross-periodogram are among the influences of the Gleissberg cycle, which takes 80–90 years [40]; this cycle is linked to global temperature changes [97].

5. Conclusions

Solar cycles were recorded in the tree rings of European beech both in Krkonoše and in the southern Apennines. A higher number of negative significant years, characterized by extremely low radial increments, were recorded in southern compared to central Europe, while a lower number of NPY was found at nutrient-poorer sites. The frequency of NPY also increased over time with increasing frequency of extreme climatic events in response to climate change, especially droughts in vegetation periods. Radial growth on individual research plots correlated differently within each time period. The most significant values were recorded for the sunspot number (compared to temperature and precipitation), which correlated profoundly with the radial growth of beech. Each research plot correlated at least once with the sunspot number. Research plots in Krkonoše responded positively to the sunspot number, seasonal temperatures and annual precipitation totals, while the research plots in the southern Apennines correlated in exactly the opposite way than those research areas in central Europe. In the first time period (1900–1969), there were positive correlations with the sunspot number on all research plots. In the second time period (1970–1985), there was a decrease in radial growth both in the Apennines and Krkonoše. In the last, third period, time period (1986–2019), there were the most significant correlations between the radial growth of beech and the sunspot number, while precipitation and temperature correlated less significantly. In total, 4- to 7-year, 9- to 12-year and 17- to 35-year cycles were found in the radial growth of beech. Cross-periodograms of beech radial growth confirmed mainly 11-year cycles of solar activity and 30-year cycles for annual precipitation totals. This research demonstrates that solar cycles manifest differently in central and southern Europe. Research of the intensity and repetition of solar cycles can be helpful in uncovering new and unexplored processes that affect radial growth across Europe. Our findings could help in understanding the adaptation of forest management to climate change. Moreover, the solar cycles might be observed in radial growth even without radiocarbon analysis. This study should improve the forestry and climatological research understanding of natural cycles and their effect on the radial growth of European beech during global climate change.

Author Contributions

V.Š. designed the research with its methodology, collected samples, measured, analyzed the data, wrote the first draft of the manuscript and made review & editing. Z.V. collaborated on manuscript, wrote the first draft and made review & editing of the manuscript. S.V. made substantial contribution with design manuscript and wrote the first draft. F.R. sourced data collection in terrain. V.H. measured data. G.D. collected samples in terrain. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Czech University of Life Sciences Prague, Faculty of Forestry and Wood Sciences (No. IGA A19/22).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Monthly temperature and precipitation data for Krkonoše Mountains are available from the Czech Hydrometeorological Institute, Prague (www.chmi.cz (accessed on 10 January 2021)). Monthly climatic data for Italy are available from the Italian Civil Protection Authority, Basilicata Region (http://www.centrofunzionalebasilicata.it/ (accessed on 10 January 2021). Annual sunspot number data are available from the Royal Observatory of Belgium, Brussels (www.astro.oma.be (accessed on 10 January 2021)). The tree ring data presented in this study are available on request from the corresponding author.

Acknowledgments

Acknowledgement goes to the Czech Hydrometeorological Institute, to Italian Civil Protection Authority and to Royal Observatory of Belgium, for providing the data sets. We would also like to thank both Richard Lee Manore, a native speaker, and Jitka Šišáková, an expert in the field, for checking the English. Finally, we thank two anonymous reviewers for their constructive comments and suggestions that helped improve the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. EUFORGEN. Fagus sylvatica, European Beech. Available online: http://www.euforgen.org/species/fagus-sylvatica/ (accessed on 10 January 2021).
  2. Dorado-Liñán, I.; Piovesan, G.; Martínez-Sancho, E.; Gea-Izquierdo, G.; Zang, C.; Cañellas, I.; Castagneri, D.; Di Filippo, A.; Gutiérrez, E.; Ewald, J.; et al. Geographical adaptation prevails over species-specific determinism in trees’ vulnerability to climate change at Mediterranean rear-edge forests. Glob. Chang. Biol. 2019, 25, 1296–1314. [Google Scholar] [CrossRef]
  3. Geßler, A.; Keitel, C.; Matyssek, R.; Seiler, W.; Rennenberg, H. Potential risks for European beech (Fagus sylvatica L.) in a changing climate Potential risks for European beech (Fagus sylvatica L.). Trees 2007, 21, 1–11. [Google Scholar] [CrossRef]
  4. Tegel, W.; Seim, A.; Hakelberg, D.; Hoffmann, S.; Panev, M.; Westphal, T.; Büntgen, U. A recent growth increase of European beech (Fagus sylvatica L.) at its Mediterranean distribution limit contradicts drought stress. Eur. J. For. Res. 2014, 133, 61–71. [Google Scholar] [CrossRef]
  5. Eilmann, B.; Sterck, F.; Wegner, L.; De Vries, S.M.G.; Von Arx, G.; Mohren, G.M.J.; Den Ouden, J. Wood structural differences between northern and southern beech provenances growing at a moderate site. Tree Physiol. 2014, 34, 882–893. [Google Scholar] [CrossRef] [PubMed]
  6. Lenz, A.; Hoch, G.; Vitasse, Y. Fast acclimation of freezing resistance suggests no influence of winter minimum temperature on the range limit of European beech. Tree Physiol. 2016, 36, 490–501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Magri, D. Patterns of post-glacial spread and the extent of glacial refugia of European beech (Fagus sylvatica). J. Biogeogr. 2008, 35, 450–463. [Google Scholar] [CrossRef]
  8. Kramer, K.; Degen, B.; Buschbom, J.; Hickler, T.; Thuiller, W.; Sykes, M.T.; de Winter, W. Modelling exploration of the future of European beech (Fagus sylvatica L.) under climate change-Range, abundance, genetic diversity and adaptive response. For. Ecol. Manag. 2010, 259, 2213–2222. [Google Scholar] [CrossRef]
  9. Penuelas, J.; Boada, M. A global change-induced biome shift in the Montseny mountains (NE Spain). Glob. Chang. Biol. 2003, 9, 131–140. [Google Scholar] [CrossRef] [Green Version]
  10. Kolář, T.; Čermák, P.; Trnka, M.; Žid, T.; Rybníček, M. Temporal changes in the climate sensitivity of Norway spruce and European beech along an elevation gradient in Central Europe. Agric. For. Meteorol. 2017, 239, 24–33. [Google Scholar] [CrossRef]
  11. Muffler, L.; Weigel, R.; Hacket-Pain, A.J.; Klisz, M.; van der Maaten, E.; Wilmking, M.; Kreyling, J.; van der Maaten-Theunissen, M. Lowest drought sensitivity and decreasing growth synchrony towards the dry distribution margin of European beech. J. Biogeogr. 2020, 47, 1910–1921. [Google Scholar] [CrossRef]
  12. Bolte, A.; Hilbrig, L.; Grundmann, B.; Kampf, F.; Brunet, J.; Roloff, A. Climate change impacts on stand structure and competitive interactions in a southern Swedish spruce-beech forest. Eur. J. For. Res. 2010, 129, 261–276. [Google Scholar] [CrossRef] [Green Version]
  13. Vacek, Z.; Vacek, S.; Slanař, J.; Bílek, L.; Bulušek, D.; Štefančík, I.; Králíček, I.; Vančura, K. Adaption of Norway spruce and European beech forests under climate change: From resistance to close-to-nature silviculture. Cent. Eur. For. J. 2019, 65, 129–144. [Google Scholar] [CrossRef] [Green Version]
  14. Tumajer, J.; Altman, J.; Štěpánek, P.; Treml, V.; Doležal, J.; Cienciala, E. Increasing moisture limitation of Norway spruce in Central Europe revealed by forward modelling of tree growth in tree-ring network. Agric. For. Meteorol. 2017, 247, 56–64. [Google Scholar] [CrossRef]
  15. Vacek, S.; Moucha, P.; Bílek, L.; Mikeska, M.; Remeš, J.; Simon, J.; Hynek, V.; Šrůtka, P.; Schwarz, O.; Mánek, J.; et al. Péče o lesní ekosystémy v chráněných územích ČR. [Management of Forest Ecosystems in the Protected Areas of the Czech Republic]; Ministry of the Environment of the Czech Republic: Prague, Czech Republic, 2012. [Google Scholar]
  16. Conte, E.; Lombardi, F.; Battipaglia, G.; Palombo, C.; Altieri, S.; La Porta, N.; Marchetti, M.; Tognetti, R. Growth dynamics, climate sensitivity and water use efficiency in pure vs. mixed pine and beech stands in Trentino (Italy). For. Ecol. Manag. 2018, 409, 707–718. [Google Scholar] [CrossRef]
  17. Etzold, S.; Ziemińska, K.; Rohner, B.; Bottero, A.; Bose, A.K.; Ruehr, N.K.; Zingg, A.; Rigling, A. One Century of Forest Monitoring Data in Switzerland Reveals Species- and Site-Specific Trends of Climate-Induced Tree Mortality. Front. Plant Sci. 2019, 10, 307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Aertsen, W.; Janssen, E.; Kint, V.; Bontemps, J.D.; Van Orshoven, J.; Muys, B. Long-term growth changes of common beech (Fagus sylvatica L.) are less pronounced on highly productive sites. For. Ecol. Manag. 2014, 312, 252–259. [Google Scholar] [CrossRef]
  19. Vacek, S.; Prokůpková, A.; Vacek, Z.; Bulušek, D.; Šimůnek, V.; Králíček, I.; Prausová, R.; Hájek, V. Growth response of mixed beech forests to climate change, various management and game pressure in Central Europe. J. For. Sci. 2019, 65, 331–345. [Google Scholar] [CrossRef]
  20. Dittmar, C.; Zech, W.; Elling, W. Growth variations of Common beech (Fagus sylvatica L.) under different climatic and environmental conditions in Europe—A dendroecological study. For. Ecol. Manag. 2003, 173, 63–78. [Google Scholar] [CrossRef]
  21. Šimůnek, V.; Vacek, Z.; Vacek, S.; Králíček, I.; Vančura, K. Growth variability of European beech (Fagus sylvatica L.) natural forests: Dendroclimatic study from Krkonoše National Park. Cent. Eur. For. J. 2019, 65, 3–11. [Google Scholar] [CrossRef] [Green Version]
  22. Králíček, I.; Vacek, Z.; Vacek, S.; Remeš, J.; Bulušek, D.; Král, J.; Štefančík, I.; Putalová, T. Dynamics and structure of mountain autochthonous spruce-beech forests: Impact of hilltop phenomenon, air pollutants and climate. Dendrobiology 2017, 77, 119–137. [Google Scholar] [CrossRef]
  23. Vacek, S.; Hůnová, I.; Vacek, Z.; Hejcmanová, P.; Podrázský, V.; Král, J.; Putalová, T.; Moser, W.K. Effects of air pollution and climatic factors on Norway spruce forests in the Orlické hory Mts. (Czech Republic), 1979–2014. Eur. J. For. Res. 2015, 134, 1127–1142. [Google Scholar] [CrossRef]
  24. Špulák, O.; Souček, J. The Sibyla model and development of beech forests affected by air pollution. Cent. Eur. J. Biol. 2010, 5, 317–383. [Google Scholar] [CrossRef]
  25. Kooijman, A.M.; Emmer, I.M.; Fanta, J.; Sevink, J. Natural regeneration potential of the degraded Krkonoše forests. Land Degrad. Dev. 2000, 11, 459–473. [Google Scholar] [CrossRef]
  26. Vacek, S.; Vacek, Z.; Bílek, L.; Nosková, I.; Schwarz, O. Structure and development of forest stands on permanent research plots in the Krkonoše Mts. J. For. Sci. 2010, 56, 518–530. [Google Scholar] [CrossRef] [Green Version]
  27. Piovesan, G.; Biondi, F.; Di Filippo, A.; Alessandrini, A.; Maugeri, M. Drought-driven growth reduction in old beech (Fagus sylvatica L.) forests of the central Apennines, Italy. Glob. Chang. Biol. 2008, 14, 1265–1281. [Google Scholar] [CrossRef]
  28. Filippo, A.D.I.; Biondi, F.; Maugeri, M. Bioclimate and growth history affect beech lifespan in the Italian Alps and Apennines. Glob. Chang. Biol. 2012, 960–972. [Google Scholar] [CrossRef]
  29. Tognetti, R.; Lasserre, B.; Di Febbraro, M.; Marchetti, M. Modeling regional drought-stress indices for beech forests in Mediterranean mountains based on tree-ring data. Agric. For. Meteorol. 2019, 265, 110–120. [Google Scholar] [CrossRef]
  30. Di Filippo, A.; Biondi, F.; Čufar, K.; De Luis, M.; Grabner, M.; Maugeri, M.; Presutti Saba, E.; Schirone, B.; Piovesan, G. Bioclimatology of beech (Fagus sylvatica L.) in the Eastern Alps: Spatial and altitudinal climatic signals identified through a tree-ring network. J. Biogeogr. 2007, 34, 1873–1892. [Google Scholar] [CrossRef]
  31. Tartaglione, N.; Toniazzo, T.; Orsolini, Y.; Otterå, O.H. Impact of solar irradiance and geomagnetic activity on polar NOx, ozone and temperature in WACCM simulations. J. Atmos. Sol. Terr. Phys. 2020, 209, 105398. [Google Scholar] [CrossRef]
  32. Lockwood, M.; Owens, M.; Hawkins, E.; Jones, G.S.; Usoskin, I. Frost fairs, sunspots and the Little Ice Age. Astron. Geophys. 2017, 58, 2.17–2.23. [Google Scholar] [CrossRef] [Green Version]
  33. Mauas, P.J.D.; Buccino, A.P.; Flamenco, E. Solar activity forcing of terrestrial hydrological phenomena. Proc. Int. Astron. Union 2016, 12, 180–191. [Google Scholar] [CrossRef] [Green Version]
  34. Antico, A.; Torres, M.E. Evidence of a decadal solar signal in the Amazon River: 1903 to 2013. Geophys. Res. Lett. 2016, 42, 10–782. [Google Scholar] [CrossRef] [Green Version]
  35. Li, H.Y.; Xue, L.J.; Wang, X.J. Relationship between solar activity and flood/drought disasters of the Second Songhua river basin. J. Water Clim. Chang. 2015, 6, 578–585. [Google Scholar] [CrossRef] [Green Version]
  36. Almedeij, J. Long-term periodic drought modeling. Stoch. Environ. Res. Risk Assess. 2016, 30, 901–910. [Google Scholar] [CrossRef]
  37. Brugnara, Y.; Brönnimann, S.; Luterbacher, J.; Rozanov, E. Influence of the sunspot cycle on the Northern Hemisphere wintertime circulation from long upper-air data sets. Atmos. Chem. Phys. 2013, 13, 6275–6288. [Google Scholar] [CrossRef] [Green Version]
  38. Ormes, J.F. Cosmic rays and climate. Adv. Sp. Res. 2018, 62, 2880–2891. [Google Scholar] [CrossRef]
  39. Easterbrook, D.J. Cause of global climate changes: Correlation of global temperature, sunspots, solar irradiance, cosmic rays, and radiocarbon and berylium production rates. In Evidence-Based Climate Science Data Opposing CO2 Emissions as the Primary Source of Global Warming, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2016; pp. 245–262. [Google Scholar] [CrossRef]
  40. Tsiropoula, G. Signatures of solar activity variability in meteorological parameters. J. Atmos. Sol. Terr. Phys. 2003, 65, 469–482. [Google Scholar] [CrossRef]
  41. Singh, A.K.; Bhargawa, A. Delineation of possible influence of solar variability and galactic cosmic rays on terrestrial climate parameters. Adv. Sp. Res. 2020, 65, 1831–1842. [Google Scholar] [CrossRef]
  42. Haywood, J.; Boucher, O. Estimates of the direct and indirect radiative forcing due to tropospheric aerosols: A review. Rev. Geophys. 2000, 38, 513–543. [Google Scholar] [CrossRef]
  43. Jayaraman, A.; Lubin, D.; Ramachandran, S.; Ramanathan, V.; Woodbridge, E. Direct observations of aerosol radiative forcing over the tropical Indian Ocean during the January-February 1996 pre-INDOEX cruise Z 10-. J. Geophys. Res. 1998, 103, 827–836. [Google Scholar] [CrossRef]
  44. Maghrabi, A.; Kudela, K. Relationship between time series cosmic ray data and aerosol optical properties: 1999–2015. J. Atmos. Sol. Terr. Phys. 2019, 190, 36–44. [Google Scholar] [CrossRef]
  45. Cecchini, S.; Galli, M.; Nanni, T.; Ruggiero, L. Solar variability and ring widths in fossil trees. Il Nuovo Cim. C 1996, 19, 527–536. [Google Scholar] [CrossRef]
  46. Muraki, Y.; Nakamura, T.; Kitagawa, H.; Masuda, K.; Nagaya, K.; Miyahara, H. Transition of solar cycle length in association with the occurrence of grand solar minima indicated by radiocarbon content in tree-rings. Quat. Geochronol. 2008, 3, 208–212. [Google Scholar] [CrossRef]
  47. Lüdecke, H.J.; Cina, R.; Dammschneider, H.J.; Lüning, S. Decadal and multidecadal natural variability in European temperature. J. Atmos. Sol. Terr. Phys. 2020, 205, 105294. [Google Scholar] [CrossRef]
  48. Bice, D.; Montanari, A.; Vučetić, V.; Vučetić, M. The influence of regional and global climatic oscillations on Croatian climate. Int. J. Climatol. 2012, 32, 1537–1557. [Google Scholar] [CrossRef]
  49. Laurenz, L.; Lüdecke, H.J.; Lüning, S. Influence of solar activity changes on European rainfall. J. Atmos. Sol. Terr. Phys. 2019, 185, 29–42. [Google Scholar] [CrossRef]
  50. Le Mouël, J.L.; Lopes, F.; Courtillot, V. A Solar Signature in Many Climate Indices. J. Geophys. Res. Atmos. 2019, 124, 2600–2619. [Google Scholar] [CrossRef]
  51. Piovesan, G.; Schirone, B. Winter North Atlantic oscillation effects on the tree rings of the Italian beech (Fagus sylvatica L.). Int. J. Biometeorol. 2000, 44, 121–127. [Google Scholar] [CrossRef] [PubMed]
  52. Šimůnek, V.; Sharma, R.P.; Vacek, Z.; Vacek, S.; Hůnová, I. Sunspot area as unexplored trend inside radial growth of European beech in Krkonoše Mountains: A forest science from different perspective. Eur. J. For. Res. 2020. [Google Scholar] [CrossRef]
  53. Šimůnek, V.; Vacek, Z.; Sharma, R.; Vacek, S. Sunspot Cycles as the Underlaying Factor of European Beech Natural Developmen. In Proceedings of Central European Silviculture; Houšková, K., Černý, J., Eds.; Mendel University: Brno, Czech Republic, 2019; pp. 267–277. [Google Scholar]
  54. Kasatkina, E.A.; Shumilov, O.I.; Timonen, M. Solar activity imprints in tree ring-data from northwestern Russia. J. Atmos. Sol. Terr. Phys. 2019, 193, 105075. [Google Scholar] [CrossRef]
  55. Shumilov, O.I.; Kasatkina, E.A.; Mielikainen, K.; Timonen, M.; Kanatjev, A.G. Palaeovolcanos, Solar activity and pine tree-rings from the Kola Peninsula (northwestern Russia) over the last 560 years Palaeovolcanos. Int. J. Environ. Res. 2011, 5, 855–864. [Google Scholar]
  56. Rigozo, N.R.; Nordemann, D.J.R.; Echer, E.; Zanandrea, A.; Gonzalez, W.D. Solar variability effects studied by tree-ring data wavelet analysis. Adv. Sp. Res. 2002, 29, 1985–1988. [Google Scholar] [CrossRef]
  57. Wang, X.; Zhang, Q. Bin Evidence of solar signals in tree rings of Smith fir from Sygera Mountain in southeast Tibet. J. Atmos. Sol. Terr. Phys. 2011, 73, 1959–1966. [Google Scholar] [CrossRef]
  58. Surový, P.; Ribeiro, N.A.; De Evora, U.; Pereira, J.S.; Superior, I.; Lisbon, D.A. Influence of solar activity cycles on cork growth—A hypothesis. In Proceedings of the 19th National Solar Physics Meeting Papradno; Dorotovič, I., Ed.; SÚH: Hurbanovo, Slovakia, 2008; pp. 67–72. [Google Scholar]
  59. Tolazs, R. Atlas Podnebí Česka; Český Hydrometeorologický Ústav: Praha, Czech Republic, 2007. [Google Scholar]
  60. Piovesan, G.; Bernabei, M.; Di Filippo, A.; Romagnoli, M. A long-term tree ring beech chronology from a high-elevation old-growth forest of Central Italy. Dendrochronologia 2003. [CrossRef]
  61. Rinntech. TSAP-Win: Time Series Analysis and Presentation for Dendrochronology and Related Applications; Rinntech: Heidelberg, Germany, 2003; Available online: http://www.rimatech.com (accessed on 20 June 2019).
  62. Larsson, L.-A. Cybis Elektronik & Data AB, Saltsjöbaden, Sweden. 2010. Available online: http://www.cybis.se/ (accessed on 20 June 2019).
  63. ČHMÚ Czech Hydrometeorological Institute. Available online: http://portal.chmi.cz/historicka-data/pocasi/uzemni-srazky (accessed on 25 June 2020).
  64. Protezione Civile—Regione Basilicata Abriola a Sellata Pierfaone. Available online: http://www.centrofunzionalebasilicata.it/it/ (accessed on 10 October 2020).
  65. WDC-SILSO Sunspot Data from the World Data Center SILSO, Royal Observatory of Belgium, Brussels. Available online: http://www.sidc.be/silso/datafiles (accessed on 20 October 2020).
  66. Team R Core. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  67. Zang, C.; Buras, A.; Cecile, J.; Mudelsee, M.; Schulz, M.; Pucha-cofrep, D. Package ‘dplR’ R, Dendrochronology Program Library in R Version 2018. Available online: https://r-forge.r-project.org/projects/dplr/ (accessed on 25 June 2020).
  68. Bunn, A.; Mikko, K. Chronology Building in dplR; R Foundation for Statistical Computing: Vienna, Austria, 2018; pp. 1–13. [Google Scholar]
  69. Cook, E.R.; Shiyatov, S.G.; Mazepa, V.S.; Ecology, A.; Branch, U. Methods of Dendrochronology Applications; Cook, E.R., Kairiukstis, L.A., Eds.; Tree-Ring Laboratory, Lamont-Ooherty Geological Observatory, Columbia University: New York, NY, USA, 1990; ISBN 9789048140602. [Google Scholar]
  70. Fritts, H.C. Tree Rings and Climate; Academic Press Inc.: Tucson, AZ, USA, 1976. [Google Scholar]
  71. Schweingrub, F.H.; Eckstein, D.; Serre-Bachet, F.; Braker, O.U. Identification, Presentation and Interpretation of Event Years and Pointer Years in Dendrochronology. Dendrochronologia 1990, 8, 9–38. [Google Scholar]
  72. StatSoft. Statistica Electronic Manual; StatSoft: Tulsa, OK, USA, 2013. [Google Scholar]
  73. Lim, C.Y.; Stein, M. Properties of spatial cross-periodograms using fixed-domain asymptotics. J. Multivar. Anal. 2008, 99, 1962–1984. [Google Scholar] [CrossRef] [Green Version]
  74. Bunn, A.; Mikko, K.; Biondi, F.; Campelo, F.; Merian, P.; Qeadan, F.; Zang, C.; Pucha-Cofrep, D.; Wernicke, J. Dendrochronology Program Library in R. R package version 1.6.8. Dendrochronologia 2018, 26, 115–124. [Google Scholar] [CrossRef]
  75. Ligges, U.; Short, T.; Kienzle, P.; Schnackenberg, S.; Billinghurst, D.; Borchers, H.-W.; Carezia, A.; Dupuis, P.; Eaton, J.W.; Farhi, E.; et al. Package ‘Signal’ 2015. Available online: http://r-forge.r-project.org/projects/signal/ (accessed on 25 June 2020).
  76. Robson, T.M.; Rasztovits, E.; Aphalo, P.J.; Alia, R.; Aranda, I. Flushing phenology and fitness of European beech (Fagus sylvatica L.) provenances from a trial in La Rioja, Spain, segregate according to their climate of origin. Agric. For. Meteorol. 2013, 180, 76–85. [Google Scholar] [CrossRef]
  77. Vacek, S.; Hejcmanová, P.; Hejcman, M.; Vacek, Z. Growth, healthy status and seed production of differently aged allochtonous and autochtonous Pinus mugo stands in the Giant Mts. over 30 years. Eur. J. For. Res. 2013, 132, 801–813. [Google Scholar] [CrossRef]
  78. Putalová, T.; Vacek, Z.; Vacek, S.; Štefančík, I.; Bulušek, D.; Král, J. Tree-ring widths as an indicator of air pollution stress and climate conditions in different Norway spruce forest stands in the Krkonoše Mts. Cent. Eur. For. J. 2019, 65, 21–33. [Google Scholar] [CrossRef]
  79. Dulamsuren, C.; Hauck, M.; Kopp, G.; Ruff, M.; Leuschner, C. European beech responds to climate change with growth decline at lower, and growth increase at higher elevations in the center of its distribution range (SW Germany). Trees Struct. Funct. 2017, 31, 673–686. [Google Scholar] [CrossRef]
  80. Gazol, A.; Camarero, J.J.; Colangelo, M.; de Luis, M.; Martínez del Castillo, E.; Serra-Maluquer, X. Summer drought and spring frost, but not their interaction, constrain European beech and Silver fir growth in their southern distribution limits. Agric. For. Meteorol. 2019, 278, 107695. [Google Scholar] [CrossRef]
  81. Nolè, A.; Rita, A.; Ferrara, A.M.S.; Borghetti, M. Effects of a large-scale late spring frost on a beech (Fagus sylvatica L.) dominated Mediterranean mountain forest derived from the spatio-temporal variations of NDVI. Ann. For. Sci. 2018, 75, 1–11. [Google Scholar] [CrossRef] [Green Version]
  82. Weemstra, M.; Eilmann, B.; Sass-Klaassen, U.G.W.; Sterck, F.J. Summer droughts limit tree growth across 10 temperate species on a productive forest site. For. Ecol. Manag. 2013, 306, 142–149. [Google Scholar] [CrossRef]
  83. Leonelli, G.; Denneler, B.; Bergeron, Y. Climate sensitivity of trembling aspen radial growth along a productivity gradient in northeastern British Columbia, Canada. Can. J. For. Res. 2008, 38, 1211–1222. [Google Scholar] [CrossRef]
  84. Urban, J.; Holušová, K.; Menšík, L.; Čermák, J.; Kantor, P. Tree allometry of Douglas fir and Norway spruce on a nutrient-poor and a nutrient-rich site. Trees Struct. Funct. 2013, 27, 97–110. [Google Scholar] [CrossRef]
  85. Lévesque, M.; Walthert, L.; Weber, P. Soil nutrients influence growth response of temperate tree species to drought. J. Ecol. 2016, 104, 377–387. [Google Scholar] [CrossRef]
  86. Šimůnek, V.; Vacek, Z.; Vacek, S. Solar cycles in salvage logging: National data from the Czech Republic confirm significant correlation. Forests 2020, 11, 973. [Google Scholar] [CrossRef]
  87. Seidl, R.; Thom, D.; Kautz, M.; Martin-Benito, D.; Peltoniemi, M.; Vacchiano, G.; Wild, J.; Ascoli, D.; Petr, M.; Honkaniemi, J.; et al. Forest disturbances under climate change. Nat. Clim. Chang. 2017, 7, 395–402. [Google Scholar] [CrossRef] [Green Version]
  88. Ferrara, C.; Marchi, M.; Fabbio, G.; Fares, S.; Bertini, G.; Piovosi, M.; Salvati, L. Exploring nonlinear intra-annual growth dynamics in Fagus sylvatica L. trees at the Italian ICP-forests level II network. Forests 2019, 10, 584. [Google Scholar] [CrossRef] [Green Version]
  89. Gray, L.J.; Woollings, T.J.; Andrews, M.; Knight, J. Eleven-year solar cycle signal in the NAO and Atlantic/European blocking. Q. J. R. Meteorol. Soc. 2016, 142, 1890–1903. [Google Scholar] [CrossRef] [Green Version]
  90. Adolphi, F.; Muscheler, R.; Svensson, A.; Aldahan, A.; Matthes, K. Persistent link between solar activity and Greenland climate during the Last Glacial Maximum. Nat. Geosci. 2014, 7, 662–666. [Google Scholar] [CrossRef]
  91. Ma, H.; Chen, H.; Gray, L.; Zhou, L.; Li, X.; Wang, R.; Zhu, S. Changing response of the North Atlantic/European winter climate to the 11 year solar cycle. Environ. Res. Lett. 2018, 13. [Google Scholar] [CrossRef]
  92. Kodera, K.; Kuroda, Y. Dynamical response to the solar cycle. J. Geophys. Res. Atmos. 2002, 107, ACL 5-1–ACL 5-12. [Google Scholar] [CrossRef] [Green Version]
  93. Brönnimann, S.; Ewen, T.; Griesser, T.; Jenne, R. Multidecadal signal of solar variability in the upper troposphere during the 20th century. Space Sci. Rev. 2006, 125, 305–317. [Google Scholar] [CrossRef]
  94. Lopez-Bustins, J.A.; Esteban, P.; Labitzke, K.; Langematz, U. The role of the stratosphere in Iberian Peninsula rainfall: A preliminary approach in February. J. Atmos. Sol. Terr. Phys. 2007, 69, 1471–1484. [Google Scholar] [CrossRef] [Green Version]
  95. Uğur, B.; Feriha, Y. Forecasting risky years for forest fires depending on sunspot cycle. J. For. Res. 2017, 4, 133–142. [Google Scholar]
  96. Kumar, V.; Dhaka, S.K.; Panwar, V.; Singh, N.; Rao, A.S.; Malik, S.; Yoden, S. Detection of solar cycle signal in the tropospheric temperature using COSMIC data. Curr. Sci. 2018, 115, 2232–2239. [Google Scholar] [CrossRef]
  97. Nagovitsyn, Y.A. Specific features in the effect of solar activity on the Earth’s climate changes. Geomagn. Aeron. 2014, 54, 1010–1013. [Google Scholar] [CrossRef]
  98. Hathaway, D.H. The solar cycle. Living Rev. Sol. Phys. 2015, 12, 83. [Google Scholar] [CrossRef] [PubMed]
  99. Singh, A.K.; Bhargawa, A. Prediction of declining solar activity trends during solar cycles 25 and 26 and indication of other solar minimum. Astrophys. Space Sci. 2019, 364, 12. [Google Scholar] [CrossRef]
  100. Kristoufek, L. Has global warming modified the relationship between sunspot numbers and global temperatures? Phys. A Stat. Mech. Appl. 2017, 468, 351–358. [Google Scholar] [CrossRef] [Green Version]
  101. Ding, Z.; Wu, R. Quantifying the internal variability in multi-decadal trends of spring surface air temperature over mid-to-high latitudes of Eurasia. Clim. Dyn. 2020, 55, 2013–2030. [Google Scholar] [CrossRef]
  102. Vacek, S.; Bílek, L.; Schwarz, O.; Hejcmanová, P.; Mikeska, M. Effect of Air Pollution on the Health Status of Spruce Stands Effect of Air Pollution on the Health Status of Spruce Stands. Mt. Res. Dev. 2013, 33, 40–50. [Google Scholar] [CrossRef]
  103. Dorotovič, I.; Louzada, J.L.; Rodrigues, J.C.; Karlovský, V. Impact of Solar Activity on the Growth of Pine Trees: Case Study. Eur. J. For. Res. 2014, 133, 639–648. [Google Scholar] [CrossRef]
  104. Kotsias, G.; Lolis, C.J.; Hatzianastassiou, N.; Levizzani, V.; Bartzokas, A. On the connection between large-scale atmospheric circulation and winter GPCP precipitation over the Mediterranean region for the period 1980–2017. Atmos. Res. 2020, 233, 104714. [Google Scholar] [CrossRef]
  105. Wibig, J.; Piotrowski, P. Impact of the air temperature and atmospheric circulation on extreme precipitation in Poland. Int. J. Climatol. 2018, 38, 4533–4549. [Google Scholar] [CrossRef]
  106. Jean, D.; Nordemann, R.; Echer, E.; Evangelista, H.; Pereira, M.; Echer, D.S. Solar activity imprints in tree ring width from Chile (1610–1991). J. Atmos. Sol. Terr. Phys. 2007, 69, 1049–1056. [Google Scholar]
  107. McCracken, K.G.; Beer, J.; McDonald, F.B. A five-year variability in the modulation of the galactic cosmic radiation over epochs of low solar activity. Geophys. Res. Lett. 2002, 29, 3–6. [Google Scholar] [CrossRef] [Green Version]
  108. Matveev, S.M.; Chendev, Y.G.; Lupo, A.R.; Hubbart, J.A.; Timashchuk, D.A. Climatic Changes in the East-European Forest-Steppe and Effects on Scots Pine Productivity. Pure Appl. Geophys. 2017, 174, 427–443. [Google Scholar] [CrossRef]
  109. Tsonis, A.A.; Swanson, K.; Kravtsov, S. A new dynamical mechanism for major climate shifts. Geophys. Res. Lett. 2007, 34, 1–5. [Google Scholar] [CrossRef]
  110. Carl, G.; Doktor, D.; Koslowsky, D. Phase difference analysis of temperature and vegetation phenology for beech forest: A wavelet approach. Stoch. Environ. Res. Risk Assess. 2013, 27, 1221–1230. [Google Scholar] [CrossRef]
Forests 12 00259 g001 550
Figure 1. Location of permanent research plots (black dots) with monthly climatic characteristics of the mean precipitation and mean temperature in the meteorological stations (black stars) in Bedřichov, the Czech Republic, and Abriola, Italy, over the whole time period of 1900–2019.
Figure 1. Location of permanent research plots (black dots) with monthly climatic characteristics of the mean precipitation and mean temperature in the meteorological stations (black stars) in Bedřichov, the Czech Republic, and Abriola, Italy, over the whole time period of 1900–2019.
Forests 12 00259 g001
Forests 12 00259 g002 550
Figure 2. Standardized ring width chronologies of European beech on research plots in the Czech Republic/Poland (a) and in Italy (b). RWI: ring width index; spl8: eight-year spline. The periods described in the time window are as follows: 1st period: 1900–1969 (regular harvesting period); 2nd period: 1970–1985 (air pollution crisis period); 3rd period: 1986–2019 (forest protection period).
Figure 2. Standardized ring width chronologies of European beech on research plots in the Czech Republic/Poland (a) and in Italy (b). RWI: ring width index; spl8: eight-year spline. The periods described in the time window are as follows: 1st period: 1900–1969 (regular harvesting period); 2nd period: 1970–1985 (air pollution crisis period); 3rd period: 1986–2019 (forest protection period).
Forests 12 00259 g002
Forests 12 00259 g003 550
Figure 3. Cross-correlation of ring width index with sunspot number, precipitation and temperature in the three time periods and the whole time period. Values are significant at p < 0.05 and marked with a circle. Lag: year (Lag −1 is the correlation of the 1 preceding year).
Figure 3. Cross-correlation of ring width index with sunspot number, precipitation and temperature in the three time periods and the whole time period. Values are significant at p < 0.05 and marked with a circle. Lag: year (Lag −1 is the correlation of the 1 preceding year).
Forests 12 00259 g003
Forests 12 00259 g004 550
Figure 4. Cross-periodogram of ring width index (RWI) with sunspot number (SUNSPOT), annual precipitation (Total prec.) and seasonal temperature (Season temp.) in the whole time period, 1900–2019.
Figure 4. Cross-periodogram of ring width index (RWI) with sunspot number (SUNSPOT), annual precipitation (Total prec.) and seasonal temperature (Season temp.) in the whole time period, 1900–2019.
Forests 12 00259 g004
Forests 12 00259 g005 550
Figure 5. Single spectral analysis of the indexed ring-width chronology for European beech in the first and third periods.
Figure 5. Single spectral analysis of the indexed ring-width chronology for European beech in the first and third periods.
Forests 12 00259 g005
Forests 12 00259 g006 550
Figure 6. Wavelet plot calculated from mean RWI (ring width index) of European beech stands in the Krkonoše mountains in Poland and the Czech Republic (left) and plots in Italy in the Apennines mountains (right). Black color represents the 90% significance level against the dark gray noise level.
Figure 6. Wavelet plot calculated from mean RWI (ring width index) of European beech stands in the Krkonoše mountains in Poland and the Czech Republic (left) and plots in Italy in the Apennines mountains (right). Black color represents the 90% significance level against the dark gray noise level.
Forests 12 00259 g006
Table
Table 1. Overview of basic site and stand characteristics of research plots in 2019.
Table 1. Overview of basic site and stand characteristics of research plots in 2019.
Name of PlotGPSAltitude
(m)		Exposure *Slope (°)Height
(m)		Diameter (cm)Volume
(m3 ha−1)		Soil TypeKöppen Classification
Chojnik 150°50′12.1″ N 15°38′27.8″ E510NW162339380Modal
Cambisol		Dfb
Rýchory 250°39′57.7″ N 15°53′05.2″ E760NE272944540Eutrophic CambisolsDfb
Sellata 340°32′21.5″ N 15°47′39.9″ E1275E263360720Epileptic Phaeozems Csa
La Lama 440°28′22.5″ N 15°45′35.2″ E1340SE72952560Haplic PhaeozemsCsa
Notes: * NW—northwest; NE—northeast; E—east; SE—southeast; Dfb—warm summer humid continental climate with large seasonal temperature differences, with warm to hot summers and cold winters; Csa—hot summer Mediterranean climate with dry summers and mild, wet winters.
Table
Table 2. Characteristics of tree-ring chronologies for European beech in research plots for 1900–2019.
Table 2. Characteristics of tree-ring chronologies for European beech in research plots for 1900–2019.
Plot NameNo. TreesAge (Min/Max)Mean RWStd.R-barESPSNRNPY
(Samples)(mm)(mm)
Chojnik 12990/1231.990.780.410.9417−
Rýchory 233103/1821.050.530.260.92111913, 1953, 2011, 2016
Sellata 34094/1522.330.690.250.95191962, 1970, 1988, 1994, 2012, 2013, 2016, 2017
La Lama 438125/2471.730.580.340.92121933, 1970, 1971, 1981, 2013, 2017
Notes: No. trees—number of trees; Age—age of youngest and oldest sample tree; Mean RW—mean ring width in mm, Std.—standard deviation in mm; R-bar—inter-series correlation; EPS—expressed population signal; SNR—signal-to-noise ratio; NPY—negative pointer years with significantly extreme low radial growth.
Table
Table 3. Correlation coefficients for the tree-ring width index (RWI) of research plots in different time periods to the sunspot number, annual temperature, vegetation season temperature, annual precipitation and vegetation season precipitation. Significant correlation values are in bold; the correlations are significant at p < 0.05.
Table 3. Correlation coefficients for the tree-ring width index (RWI) of research plots in different time periods to the sunspot number, annual temperature, vegetation season temperature, annual precipitation and vegetation season precipitation. Significant correlation values are in bold; the correlations are significant at p < 0.05.
Years1900–20191986–20191970–19851900–1969
Time PeriodWhole Period3rd Period2nd Period1st Period
Sunspot number
RWI Chojnik 10.260.370.070.31
RWI Rýchory 20.230.54−0.120.19
RWI Sellata 3−0.25−0.37−0.430.11
RWI La Lama 4−0.07−0.02−0.510.11
Annual temperature
RWI Chojnik 10.140.180.370.10
RWI Rýchory 20.050.280.370.07
RWI Sellata 3−0.040.02−0.07−0.02
RWI La Lama 40.08−0.230.12−0.08
Vegetation season temperature
RWI Chojnik 10.090.090.390.11
RWI Rýchory 20.080.280.320.15
RWI Sellata 3−0.060.04−0.21−0.02
RWI La Lama 4−0.12−0.390.12−0.14
Annual precipitation
RWI Chojnik 1−0.100.00−0.28−0.06
RWI Rýchory 2−0.18−0.10−0.48−0.19
RWI Sellata 30.100.050.320.06
RWI La Lama 40.150.020.110.24
Vegetation season precipitation
RWI Chojnik 1−0.080.04−0.21−0.07
RWI Rýchory 2−0.15−0.16−0.45−0.17
RWI Sellata 30.11−0.010.450.13
RWI La Lama 40.120.060.200.14
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Šimůnek, V.; Vacek, Z.; Vacek, S.; Ripullone, F.; Hájek, V.; D’Andrea, G. Tree Rings of European Beech (Fagus sylvatica L.) Indicate the Relationship with Solar Cycles during Climate Change in Central and Southern Europe. Forests 2021, 12, 259. https://doi.org/10.3390/f12030259

AMA Style

Šimůnek V, Vacek Z, Vacek S, Ripullone F, Hájek V, D’Andrea G. Tree Rings of European Beech (Fagus sylvatica L.) Indicate the Relationship with Solar Cycles during Climate Change in Central and Southern Europe. Forests. 2021; 12(3):259. https://doi.org/10.3390/f12030259

Chicago/Turabian Style

Šimůnek, Václav, Zdeněk Vacek, Stanislav Vacek, Francesco Ripullone, Vojtěch Hájek, and Giuseppe D’Andrea. 2021. "Tree Rings of European Beech (Fagus sylvatica L.) Indicate the Relationship with Solar Cycles during Climate Change in Central and Southern Europe" Forests 12, no. 3: 259. https://doi.org/10.3390/f12030259

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop