Next Article in Journal
Effects of Forest Fragmentation on the Volume of Wood Resources in Managed, Pine-Dominated Forests in Poland
Previous Article in Journal
Measuring Radial Variation in Basic Density of Pendulate Oak: Comparing Increment Core Samples with the IML Power Drill
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Environmental Effects among Differently Located and Fertile Sites on Forest Basal-Area Increment in Temperate Zone

1
Department of Geology and Soil Science, Faculty of Forestry and Wood Technology, Mendel University, Zemědělská 3, CZ-613 00 Brno, Czech Republic
2
Global Change Research Institute CAS, Belidla 986/4a, CZ-603 00 Brno, Czech Republic
3
Forest Management Institute Brandýs nad Labem, Nábřežní 1326, CZ-250 01 Brandýs nad Labem-Stará Boleslav, Czech Republic
4
Institute of Physics in Opava, Silesian University in Opava, Bezručovo náměstí 1150/13, CZ-746 01 Opava, Czech Republic
5
Centre for Environment and Land Assessment—Ekotoxa, Otická 37, CZ-746 01 Opava, Czech Republic
*
Author to whom correspondence should be addressed.
Forests 2022, 13(4), 588; https://doi.org/10.3390/f13040588
Submission received: 4 February 2022 / Revised: 5 April 2022 / Accepted: 6 April 2022 / Published: 8 April 2022
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
Environmental properties differently influence the growth of forest tree species. The antagonistic effects of variable environmental properties classify the forest response according to various tree compositions among different sites. The division of the forest response was assessed in 52 stands arranged into 26 types of 13 site management populations (MPs) in 5 areas in the Czech Republic territory. The assessment was performed using time-series multiple regressions of basal-area increment from pure immature stands of Norway spruce (Picea abies), Scots pine (Pinus sylvestris), oaks (Quercus sp.), ash (Fraxinus excelsior) and willows (Salix sp.) dependent on the interpolated average temperatures, annual precipitation, atmospheric concentrations of SO2, NOx and O3 and soil properties over the period 1971–2008 at p < 0.05. Site MPs differentiated the forest response to a greater extent than tree species. The response of the forests was significantly distributed by means of the montane, upland and waterlogged sites. The multiple determination index (r2) ≥ 0.6 indicated an adaptable tree increment but an interval of r2 between 0.80–0.92 implied forest sensitivity to variability in environmental properties on non-waterlogged sites. The index r2 < 0.6 suggested a fluctuating forest increment that reflects environmental variability inconsistently. The fluctuating increment most affected the spruce and pine stands grown from upland to submontane locations. Montane spruce stands, as well as rock pines, appeared to be one of the most sensitive ones to environmental change. Floodplain forests seemed as adaptable to variable environmental properties.

1. Introduction

Environmental properties differently affect plant growth. While temperatures and water availability directly support plant growth, air pollution disrupts it. However, individual pollutants influence plant growth differently, depending on responses and physiological activity [1]. Plants partially surmount unfavourable environmental properties that do not exceed critical levels by utilising surpluses of favourable environmental properties [2]. The substitution of untoward properties is manifested by similar plant growth simultaneously in favourable and unfavourable conditions. Therefore, the estimation of the combined interaction effects between atmospheric and soil properties indicate ways of living, and community adaptation to environmental change [3].
The environmental change impacts on plant communities are most distinctively regulated by solar activity. Solar radiation and available water are essential to trigger plant photosynthesis [4]. The effects of insolation alter during solar cycles, but they also vary with the geographical location, time of day and season. Diverse insolation of variously located areas conditions albedo is related to natural organic matter production intensity [5]. However, the albedo decline is caused by both warming and ecosystem transformation. Nevertheless, plant growth is more substantially directed by insolation cycles than by temperatures. The tree species growth is slightly more temperature-dependent within colder seasons than warmer ones [6]. Warmer seasons are accompanied by the far-reaching effects of the atmospheric–oceanic circulation directed by solar activity, while within colder seasons, the distant circulations are replaced by local influences [7].
The compositions of atmosphere and soil prefigure the efficacy of the photosynthetically obtained energy used to create organic matter. The growth processes in plant communities are interconnected with the development of soil carbon content. Soil carbon development simultaneously affects soil fertility and climate variability [8]. The ecosystem carbon cycle transfers the impacts of global environmental changes to the chemical composition of organic matter due to the influence of solar activity cycles on the isotopic composition of plant bodies. The period of the reduced solar activity is manifested by an increase in 14C isotope content in plants, and vice versa [9]. Carbon interconnects the ecosystem processes by water retention and microbial interactions that reduce C/N during nitrogen access [10]. The soil C/N, bulk density (Dd) of individual soil groups, as well as a plant stand density, divide CO2 flow most. At the same time, the development of soil pH is immediately followed by alterations in the Ca/Al, Ca/Mg and Ca/Mn ratios of the wood. Decreasing soil Ca2+ availability enhances the Al3+ content of the wood, which reduces the Ca/Al ratio [11]. Similarly, decreasing soil pH reduces the Mg uptake by plants, therefore, the Ca/Mg ratio is significantly increased [12]. However, soil C/N and Dd control up to 87% of the dispersion in the CO2 flow from forests [13].
The direct effects of solar radiation on plant growth have been obscured by a cultural restructuring of ecosystems, the intensification of the greenhouse effect and the eutrophication of the environment. The restructuring of the terrestrial ecosystems has had a heavy impact on forests by fragmentation and simplification of the species composition [14]. The diminution in the Earth’s surface albedo with the retreat of montane glaciation or deforestation and greenhouse gas emissions accelerated global warming at the end of the Little Ice Age after 1850 [15]. Whereas the retreat of glaciation was naturally caused by the intensification of the solar activity and the rise in summer temperatures, other reasons for the albedo decrease ensued the increasing need of the human race to utilise landscape [16]. Current global warming is unprecedented throughout the Holocene due to human activity [17]. The present-day environmental change has accelerated forest tree species’ growth more than the previous Medieval Climatic Optimum [18]. Differences in the forest tree species’ growth between colder and warmer areas have decreased [19]. Nevertheless, environmental change does not cause a general improvement in growth conditions, but the overwhelming majority of nutrient-rich sites are threatened by erosion, drainage and acidification [20].
The threat to growth conditions affects not only the forest growth intensity but also causes an imbalance in plant nutrition. Fragmentation and transformation of the forest species’ composition have increased soil organic matter losses and CO2 emissions. The subsequent decrease in the return of nutrients from plant shedding to the soil elicits an ecosystem stability loss to mitigate the acidification process [2]. The course of forest ecosystem acidification differs after exceeding critical loads by acid deposition or owing to reduced water availability. The acid deposition load gradually disrupts the ecosystem. The disturbance initially impacts the soil substance interchange until available base cations neutralising the incoming acids are depleted. Only then, does the forest status deteriorate [21].
Restoration of reduced soil fertility as well as adaptation to global warming has become dependent on sufficient soil nitrogen supply to maintain optimally low C/N [22]. Eutrophiering growth of atmospheric CO2 content and physiologically active nitrogen compounds increases a forest increment merely in the conditions of soil nutrient excess and subdued acidification. Otherwise, eutrophication causes an imbalance between individual substance sources, reducing plant growth [23]. Nevertheless, the growth of the main management tree species at nutrient-rich sites has influenced an increase in the overall increment of European forests since the end of the Little Ice Age [3]. The diameter increment of boreal forests increases by 25% at optimal soil conditions under 2 × CO2 atmosphere [24].
The aim of the presented analysis of variable growth condition effects on the tree species diameter increment was to assess forests’ responses to environmental change among different sites. Usual observations of effects on tree species growth were focused either on one tree species in different conditions e.g., [6,11,23,25] or on several tree species under the same conditions [4,25]. The fundamental methods of the tree species growth assessment do not allow for the evaluation of responses in mixed stands to differently acting environmental properties. Variable atmospheric or soil properties influence forest tree species diversely. A diverse effect is indicated by distinctive parameters at correlation functions [8]. Similar correlation parameters suggest forest resilience under various conditions. Nevertheless, the distribution and structure of planted forests prone to acidification or living pests complicate the assessment of the impacts of alterations in growth conditions [23]. The evaluation of the development in unevenly distributed forests became dependent on the growth condition properties’ homogenisation. Homogenisation consists of determining the mean courses of variable properties in generalised growth conditions [26]. The generalisation process is based on defining a repeatable series of the landscape geographical division within the unrepeatable area. Each unique area is characterised by a uniform macroclimate; therefore, its structure develops alongside types of relief and rocks [27]. Although altitudinal transitions divide forest growth conditions most significantly, most rock types considerably predetermine soil properties [28]. Within the generalised division of growth conditions, the resilience was assessed between natural and planted stands in the different locations and fertile series. The analysis of the forest growth conditions included both a comparison of various tree species’ increments at the same site series and a comparison of the increment of one tree species among various site series [29]. The comparison among different tree species suggests adaptability to the site, while a comparison of various sites suggests a dispersion in adaptation ability.

2. Material and Methods

2.1. Geographical Pattern

Growth condition effects on the diameter increment of pure forest stands were assessed in terms of atmospheric and soil properties. The growth conditions were characterised by the selection of dominant sites defining the natural environment of the Czech Republic in the Central European landscape. The Czech Republic is situated in a hinterland of transitional temperate climate (78,866 km2; 115–1602 m a.s.l.; 2.0–9.4 °C; 250–1470 mm). The territory of the CR consists mainly of the Bohemian Massif (84.6%) which is covered by the Outer Western Carpathians in the east (15.4%). The Bohemian Massif differs from other Central European montane systems by means of a central basin permeated by rock cities and volcanic mountain ranges [30].
The natural forest cover of the CR exceeds 98%, with the exception of high-mountain, rocky or peatland treeless areas. The current forest cover is 34%. The natural representation of tree species consisted mainly of European beech (Fagus sylvatica) (40%), oaks (Quercus sp.) (18%), Silver fir (Abies alba) (16%), Norway spruce (Picea abies) (15%) and Scots pine (Pinus sylvestris) (3%) [31], while the current forest composition is mainly composed by 44% of Norway spruce, almost 10% of Scots pine, more than 10% of European beech and 8% of oaks [32].
Forest growth conditions in the Czech Republic are divided into 14 unique areas and 27 management populations (MP) of forest types [33]. Each MP is characterised by the altitudinal zone of climatic conditions and by soil series of acidity and water availability. For the environmental change impact analysis, the forest stands were selected at MPs dominating the Outer Western Carpathians (N), the Eastern Sudetes Mountains (K), the Elbe River basin (H), the North Bohemian Sandstone Highlands (I) and the volcanic mountain ranges (B). Furthermore, the Outer Western Carpathians are a flysch highland, where mountainous areas are delimited by the Moravian-Silesian Beskids, the Eastern Sudetes Mountains are mostly formed by the crystalline mountains of Hrubý Jeseník and Králický Sněžník. The North Bohemian Sandstone Highlands and the volcanic mountain ranges are not significantly different in altitude, but they contrast with fertility (Figure 1).
The selection of the forest stands was performed in the MPs, covering >5% of the total forest soil area in each characteristic area. The stands were of managementally significant Norway spruce, Scots pine and oaks, as well as the azonal ecosystems of European ash (Fraxinus excelsior) and willows (Salix sp.) surviving Quaternary climate changes [34]. Zonal MPs were represented by the site-optimal stand-forming tree species and compared with the managementally transformed composition. Optimal stand-forming tree species for the individual MPs were determined from estimates of potential natural vegetation in the database of the Regional Forest Development Plans (RFDP) [31]. In each MP, the individual stand types were sampled in pairs, allowing for the minimised calculation of statistically representative averages and deviations [35]. The selection was limited to the single-stage closed stands older than 40 years with a representation of the investigated tree species > 80% within an area > 1 ha [36].

2.2. Field Sampling

The selected 52 forest stands were surveyed at the representative plots 30 × 30 m in May–June 2009. The representative plot was located in the geometric centre of the stand. The survey consisted of soil and dendrochronological sampling. An altitude of the surveyed plots was determined as the average through an overlay between vectors of the forest stand representative part and the contour lines spaced 1 m apart in the Basic Base of Geographical Data managed by the Czech Surveyor and Cadastral Administration [37]. The bedrock was typed according to the biogeographical division in the vector model 1:50,000 managed by the Nature Conservation Agency of the Czech Republic [38], while the soil groups were evaluated according to WRB-ISSS-ISRIC [39].
The soil survey included the determination of clay content, pH/H2O, base saturation (BS), Corg, Al2O3 and CaO in the upper horizons from five regularly located trenches in each plot. Clay content was determined sedimentographically; pH, acidometrically; BS as the proportion of the base cation contents from the total amount of exchangeable cations by extraction in 0.1M BaCl2 at pH 7.8; Corg was determined by infrared spectrophotometry; while the Al2O3 and CaO contents were determined by AAS from extraction in aqua regia [40].
A dendrochronological sampling of the representative area commenced with the selection of 30 sample trees with the dimensions of the mean stem with a deviation of the breast-height diameter (d1.3) ± 2 cm. The tree-ring series were sampled by Pressler borer with a 5-mm inner diameter at the measurement height d1.3 in a direction parallel to the slope, avoiding reaction wood [41]. The age, mean diameter, height (h), volume (Vha) and canopy of the monitored stands were determined by the forest management plans.

2.3. Time-Series Analysis

The response of each forest stand was assessed by comparison between the time series of current basal-area increment (BAI) and the series of annual average temperatures (T), annual precipitation (P), atmospheric concentrations of SO2, NOx and O3 and soil properties in the period of 1971–2008 (Figure 2). The period length was adapted to the commencement of the assessment of instrumentally measured atmospheric pollution effects on immature forests during climate change [34,42,43]. BAI is a generalisation of the diameter increment which characterises stand density at scale, indicating the effects from the external environment [23]. Climatic properties were interpolated using modified IDW from a rural measurement station system in investigated forest areas [44]. Concentrations of air pollutants were interpolated using ordinary kriging from the Information System on Air Quality of the Czech Hydrometeorological Institute stations [45]. Pollution concentrations before instrumental measuring commencement were derived linearly according to estimations of pre-industrial biogeochemistry [46].
Atmospheric SO2 concentrations were the most serious environmental load in Central Europe until 1994 [47]. Impacts of NOx and O3 loads worsened significantly, in particular after 2000, due to the intensification of climate change [1,48]. The selection of soil properties from sets of physical and chemical characteristics focused on the attributes correlating with forest health statuses. While pH and BS are fundamental physico-chemical properties indicating acidification and related changes in fertility, the soil clay content forms a matrix maintaining a constant level of base saturation, however, Corg affects soil sorption variability [49]. Al2O3 and CaO contents are components of weathering that are adversely released or bound during acidification [26]. The development of soil properties was generalised by average values characterising particular MPs within selected forest areas in the RFDP database (Figure 3).
BAI was obtained from the mean inter-year difference between ratios of standardised stand canopy and basal area (Gtab) modelled at yield tables [50]. The development of the canopy was characterised by the ratios between growing Vha and modelled Vtab. Previous Vha were derived from ratios of modelled tree number development and mean stem volume through nonlinear models of dependences between tree-ring increments (d1.3) and stand height growth [51]. The d1.3 widths were measured using a VIAS TimeTable measuring system (SCIEM) with an accuracy of 10 μM. Values of d1.3 were synchronised visually, cross-dated by PAST 4 [52] and controlled statistically using COFECHA [53].

2.4. Statistical Assessment

The statistical assessment of forest response was divided into exploratory analysis and multiple regressions. The exploratory analysis consisted of linear correlations and residue analysis. While linear correlations indicated the effects of the individual growth condition properties on stand type BAI, the residue analysis was designed to reveal inaccuracies in the linear models. Linear correlations were evaluated as significant at p < 0.05, less significant between p ≥ 0.05 and p < 0.50, and insignificant. Residue analysis was performed using the residual standard deviation (RSD), elevation (ER) and skewness (AR). The residues were calculated as differences between obtained and modelled BAIs. The low linear approximation was detected in ER > RSD4 or AR > RSD3 [54]. Multiple regressions were calculated for each plot, but merely the parameters of the most relevant function from a pair of the same stand types were considered at p < 0.05. Independent environmental variables were normalised by z-transformation expressing deviations of the measured values from the mean. The multiple model was compiled as follows:
BAI = i = 1 m a i . x i + b
where xi is the z-transformed growth condition property, ai is the slope parameter of the growth condition property, b is the intercept parameter, i is the sequence of the growth condition property and m is the number of the growth condition properties. Functions among the individual stand types were subsequently contrasted using minimal Fischer–Snedecor criteria (Fmin) and determination indices (r2) [55].

3. Results

3.1. Management Populations of Forest Growth Conditions

Management populations of forest soils are unevenly represented in the individual areas of the Czech Republic. Simultaneously, no MP is dominant in all the forest areas of the Czech Republic. Only six MPs exceeded their representation of 5% of the forest soil area. The six predominant MPs include 62.9% of the forest soils in the Czech Republic. The acidic sites predominate from the uplands to the mountains, while the nutrient-rich sites are only from the highlands to the mountains. The upland acidic series includes 5.7% of the forest soils, highlands (10.7%) and montane (11.5%). The nutrient-rich series includes 17.9% in the highlands while merely 11.1% in the montane locations. The montane sites are characteristic of a larger proportion of the gleyed series (5.9%) (Table 1). The upland acidic sites are most concentrated in the area of the Moravian floodplains (M). The highland acidic sites are most concentrated in the Beroun areas (C) and the Central Bohemian areas (D) while the highland nutrient-rich sites are most concentrated separately in the area of the volcanic mountain ranges (B) and further surroundings of the Moravian floodplains in the Bohemian-Moravian intermountain (L) and Eastern Sudetes (K), to the Carpathian (N) areas. The montane acidic sites predominate in the largest mountain ranges in the territory of the Czech Republic from the Ore Mountains (A), the Šumava Range (E), and the Bohemian-Moravian Highlands (G) to the Western Sudetes (J). The montane nutrient-rich to gleyed sites do not dominate in any of the Czech forest areas, but mostly, except for the gleyed series in the Eastern Sudetes and the Carpathians, they come after the dominant soil series representation.
Selected forest areas demonstrating the diversity of the Czech Republic are specific for four dominant site MPs. Forest soils in the selected areas are most represented by upland (9.5%), highland (28.1%) montane (16.8%) and nutrient-rich sites and rock towns (8.1%). Highland nutrient-rich sites dominate in three of the five selected areas. They cover 28% of the forests in the volcanic mountain ranges but more than 44% of the forests in the Outer Western Carpathians. Volcanic mountain ranges are typical of the common dominance of the upland (16.4%) to highland (28.0%) nutrient-rich sites, while the Elbe River basin (H) is mainly formed by acidic (21.7%) to nutrient-rich (34.7%) sites in the uplands. Contrarily, the Eastern Sudetes and the Outer Western Carpathians are similarly characterised by a predominance of nutrient-rich series from the highlands (34.7% and 44.2%) to montane altitudes (25.0% and 23.5%). The included North Bohemian Sandstone Highlands (I) differ, mainly in the predominance of rock cities (40.3%) followed by highland acidic soils (14.7%).

3.2. Properties of Assessed Ecosystems

3.2.1. Growth Conditions

The survey was conducted in 26 forest stand types at 13 MPs of four altitudinal zones and 11 soil series. The soil series was composed of 11 bedrock types and 15 soil groups. The selected forest stand types were divided in parallel between the Bohemian Massif and the Outer Western Carpathians. Conditions of the Bohemian Massif were characterised by 10 stand types in the upland areas, 7 stand types in the montane areas, and 2 in the rock cities. The Western Carpathians were defined by 2 types of floodplain forests and 5 types of montane forest stands. The montane sites were similarly defined by foothill to ravine, exposed, nutrient-rich and the peak timberline series, whereas the uplands were evaluated by means of the correlation between flat, exposed, nutrient-rich and the gleyed series.
The selected forest stands were located in the range of altitudes between 192–1173 m. The lowest-lying one was the pine stand in the upland gleyed series. The highest-lying one was the montane spruce stand in the Outer Western Carpathians. The extreme values of the growth conditions intervened several times in floodplain willow beds, gleyed upland pine and oak stands and montane spruce forests. While the lowest temperatures of +3.51 °C occurred at the sites of the Sudeten montane spruce forests; the lowest SO2 and NOx concentrations and the highest O3 load appeared in the Carpathian montane spruce forests. On the contrary, the gleyed upland sites of the pine and oak stands were most burdened by SO2 and NOx atmospheric pollution, however, only by minimal concentrations of tropospheric O3. Similarly, the lowest annual precipitation of 550 mm was detected in the nutrient-rich sites of the upland oak stands while the highest precipitation of more than 1280 mm accompanied the Carpathian ravine sites (Table 2).
The marginal values of the soil properties most often accompanied the nutrient-rich sites in contrast to the rock cities. The most frequent bedrock of the selected forests was formed by sandstones. Sandstones centrically occurred in the Carpathian forests and in the North Bohemian steep reliefs. In this context, Arenic Podzols and Skeletic Cambisols were more frequent than Haplic Cambisols, concentrated more on igneous rocks. However, Luvisols, Luvic Cambisols and Entic Podzols associations, which are predominant in the inland uplands, were similarly common. The lowest soil clay content of 1.1% was discovered in the exposed Sudetes sites of the spruce forests, while the highest soil clay content of 31.5% was determined at the nutrient-rich sites of the upland areas in oak stands. The lowest soil pH/H2O 3.7 was ascertained at the Carpathian nutrient-rich sites while the highest pH/H2O 6.1 was discovered in the submontane floodplains. Nevertheless, the highest pH corresponded to the highest soil base saturation at the floodplain sites (97.7%), whereas the lowest BS (11.3%) was identified at the foot of the Outer Western Carpathians. The values of the soils’ chemical properties followed the distribution of bedrock types rather than relief, such as the pH or BS. The lowest Corg content was demonstrated in the natural pines from the rock city areas (0.75%) while the highest Corg content was determined in the Carpathian montane spruce forests (17%). Similarly, the lowest Al2O3 content was discovered in the soils from the rock cities (only 0.07%) but the highest content of total aluminium was identified in the nutrient-rich sites of the Elbe oak stands (5.6%). On the contrary, the lowest soil CaO contents of 0.01% partially corresponded to the occurrence of the lowest pH and BS in the Outer Western Carpathians, but the highest CaO content was recorded in the Elbe exposed oak groves (2.09%) (Table 3).

3.2.2. Stand Types

The characteristics of the selected stand types were most often distributed along the uplands and mountain ranges. The floodplain forests were characterised by features similar to the submontane stands, although their location was related to the upland sites. Edge stand characteristics values occurred more often in the willow woods, ravine pines and exposed Carpathian forests. The willow stands were the youngest ones at the time of the survey (45 years) while the North Bohemian ravine pines were the oldest ones (123 years). Although most of the surveyed forest stands were younger than 100 years, natural rock pines, exposed upland pines and spruce stands, as well as nutrient-rich oak groves, were older than 100 years.
The low age foreshadowed the lowest ascertained stand height of 14 ± 8 m in the willow stands but did not significantly affect the maximum height > 30 m in the exposed Carpathian spruce stands. The mean d1.3 was smallest in the planted rock pines (22 cm), while on average, the largest d1.3 > 39 cm was detected in the spruce stands at the foot of the Carpathians. The mean d1.3 > 30 cm was most frequently discovered in the mountain forests, whereas smaller mean trunk diameters predominated in the upland forests. The stand height affected marginal increment values. The lowest basal-area increment was <2.5 m2/ha.year and occurred in stands with a height of <20 m. On the contrary, the highest increments > 8 m2/ha.year occurred in stands higher than 20 m on rich riverine or foothill sites (Table 4; Figure 4).

3.3. Response of Forest Stands

The significant correlations of more than 50% of the growth condition properties qualified multiple regressions explaining > 50% of the BAI variance. However, the significance of multiple regressions was not reduced by more frequent insignificant simple correlations.
Temperatures, atmospheric NOx and O3 concentrations and the soil clay content significantly affected the most evaluated stand types. Precipitation, temperatures, atmospheric SO2, soil clay, pH, base saturation, Al2O3 and CaO impacted BAI in more than 75% of the stand types. On the contrary, the development of soil Corg significantly influenced merely less than 12% of the compared stand types. Minor significant correlations involving < 50% of the properties were identified in 73% of the stand types while the predominant minor significant correlations of most ecosystem properties affected only 27% of the stand types. BAI of the planted rock pines, gleyed pine stands, nutrient-rich oak forests and upland gleyed spruce stands, as well as the wet montane spruce stands and Carpathian ravine forests, were influenced by the significant or low correlations covering more than 90% of the growth conditions. Significant correlations of more than 50% of the evaluated growth conditions impacted BAI not only in the Carpathian ravine forests but also in the East-Sudeten montane spruce forests, East-Sudeten gleyed spruce stands and the gleyed oak stands. The most influential growth conditions were components of acid deposition and soil physico-chemical pH and BS, partly also temperatures, precipitation, soil clay, Corg, Al2O3 and CaO (Table 5).
Multiple regressions among the BAI of the individual stand types and growth conditions included the range r2 0.16–0.92. The lowest multiple r2 occurred in the upland exposed series. The highest r2 was ascertained on ravine Carpathian slopes. The smallest difference between incremental r2 within one stand type was determined at pine stands planted in rock cities. The greatest difference in r2 was detected between the spruce stands in flat uplands. The site distributed multiple r2 more distinctively than tree species. The Carpathian forests varied from the Eastern Sudetes by a greater r2 BAI function within all soil series. The regressions of the Carpathian forests were more similar due to the positive parameters of SO2 and soil pH, while East-Sudeten montane forests were dependent on the same parameters for atmospheric precipitation, SO2 and soil pH, BS as well as Corg. In general, the residual normality was maintained in the conditions of the Bohemian Massif. The most similar residues occurred in pines, flat and exposed uplands, as well as in the montane wet and timberline sites.
On average, the highest values of r2 > 0.7 of multiple regressions were detected in all the investigated montane forests and in the North Bohemian natural pines separately. The non-waterlogged submontane natural forests were characterised by the tightness of multiple regressions of r2 < 0.6, but the nutrient-rich forests reached an r2 between 0.60–0.77. The naturally upland oak stands evinced the most significant regression of BAI at the optimal nutrient-rich sites. By contrast, the spruce stands were less dependent on the condition development than the pine and oak ones in all the intra-Bohemian soil series. The higher dependence in natural upland forests was conditioned by differentiated parameters at atmospheric properties and by equally directional parameters of soil clay, pH, BS, Corg and CaO (Table 6).

4. Discussion

Multiple dependence of the current radial increment on external conditions was most significantly distributed along different sites. On the other hand, the natural forest stands are more closely reflected by the effects of changing site conditions than unnatural tree species populations. The dependence proximity indicated the magnitude of the forest tree species increment response to multiple variabilities of growth conditions.
The forest tree species increment naturally depends on environmental properties as well as on the competition of individual species. The assessment of different tree species' responses has confirmed diversification by fertility, regional climate change and atmospheric pollution [3,24,29]. Environment significantly influenced the increment of most of the compared forest stand types. However, the parameters of individual environmental properties differed in size and direction among stand types. Atmospheric precipitation, soil clay, pH, Corg and Al2O3 influenced most forest stands to the same extent. Tropospheric O3 and soil Corg and Al2O3 significantly affected forests in the individual areas contradictorily. The opposing effect is related to the geographical location, stimulating diverse impacts of regional climate change together with age and forest stand density [56].
Air pollution has reduced tree species increments, particularly in mountain altitudes. Nonetheless, the individual pollutants acted opposite among the different nutrient-rich sites from the uplands to the mountains. Atmospheric SO2 has significantly affected the natural pine forests, intra-Bohemian gleyed forests, and East-Sudeten ravine to exposed ecosystems, floodplain ash stands and Carpathian ravine forests. NOx pollution has affected several monitored stands from the natural pines, nutrient-rich and gleyed upland forests and waterlogged East-Sudeten ecosystems to the timberline. In the Outer-Carpathian areas, the willow floodplains and exposed slopes were particularly impacted. Tropospheric O3 significantly influenced (especially) nutrient-rich to gleyed upland communities in the Eastern Sudetes. The increment of the Carpathian forests was most impacted by tropospheric O3 concentrations in the nutrient-rich exposed sites. The development of SO2 concentrations reduced forest increment most significantly in the gleyed and exposed sites from the uplands to the mountains while the NOx course contributed to the forest increment increase in lower altitudes and reduced the increment of montane forests. Contrarily, O3 in the lower altitudes mostly reduced forest increment but did not affect selected montane forests markedly.
The adverse effects of atmospheric pollution are associated with plant–climate interactions. The temperature development mostly influences montane forest growth, whereas upland forests depend more on available soil water [57]. The climate change impacts are modified by the complicated arrangement of mountain relief. Relief contributes to forest vulnerability with damage exceeding 32% of the tree species population. Forests growing at altitudes ≤ 950 m a.s.l. are affected more by drought than by pollution but the forest development > 950 m a.s.l. was impacted by the long-distance transmission of pollution and their increment was increased after the decrease in acidic load due to global warming [58]. The altitude exceeds 90% of the relief risk while the slope, exposure and relief shape include only less than 7% of the predispositions. Montane forest damage is mainly concentrated on ridges, higher slope parts, round valleys and northern to north-western exposures due to the relief effects [59]. Simultaneously, the climate change impacts on forests are more distinctive in mountain locations than in the uplands. On the one hand, climate change increases tree species increment and prolongs the growth dynamics of older stands. On the other hand, it increases the frequency of sudden weather disturbances [4,25,43].
On the one hand, the close relationship between tree increment and environmental properties indicated compliance with the site, but on the other hand, it indicated sensitivity to the change in any environmental properties. The increment of the nutrient-rich upland forests was dependent on environmental variability more than on the mountain altitudes. The unnatural pine stands in the upland locations differed slightly in their dependence on the environmental properties due to natural oak groves. Nevertheless, the pine stands grown in rock cities were significantly less dependent on the environment than natural pines. Scots pine appears to be very sensitive to growth condition variability as it is most sensitive to drought from the Central-European managementally important tree species [6]. Even though Norway spruce or European beech may also be susceptible to drought, the extended growing season supported their growth in the humid montane altitudes [29]. While the growth of the natural pines followed fluctuating atmospheric precipitation and higher O3 load, the oak stands under natural conditions reacted by means of an opposite decrease in their dependence with a more variable environment on exposed sites and by an increase on nutrient-rich sites. The influence of nutrient-rich sites on tree species increment decreased with the altitude. On the contrary, the increment of upland gleyed forests was contingent on the environment less than in mountain altitudes. The waterlogged oak groves, including floodplain forests, appear to be more resistant to environmental change in upland altitudes due to a sufficient supply of soil water to the sites [60].

5. Conclusions

The atmospheric and soil properties influenced diameter increment development of Central-European forest tree species differentially among the montane, upland and waterlogged sites. Temperatures, acid deposition and soil physico-chemical properties, Al2O3 and CaO significantly affected more than 75% of the stand types. The forests adaptable to variable growth conditions were characteristic by multiple regression r2 ≥ 0.60. Spruce and pine stands growing in upland altitudes seemed to be more susceptible to growth condition variability than natural forests. Resistant stands were defined by natural forests at sites sufficiently supplied with water. Susceptibility to the environmental change was higher at higher altitudes. High-mountain forests, as well as the natural pines, seemed to be one of the most sensitive species to environmental change.

Author Contributions

Conceptualization, P.S.; methodology, P.S., M.Z. (Miroslav Zeman) and M.Z. (Miloš Zapletal); validation, P.R. and M.Z. (Miroslav Zeman); formal analysis, P.S. and P.R.; investigation, P.S. and P.R.; resources, P.S.; data curation, P.R.; writing—original draft preparation, P.S.; writing—review and editing, P.S. and M.Z. (Miloš Zapletal).; visualization, P.S.; supervision, M.Z. (Miloš Zapletal). All authors have read and agreed to the published version of the manuscript.

Funding

The study has received funding from the European Union's Horizon 2020 Programme for Research & Innovation under grant agreement No 952314 ASFORCLIC.

Data Availability Statement

Restrictions apply to the availability of data used. Data was obtained from third party and are available from the authors with the permission of data owner.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Augustin, S.; Bolte, A.; Holzhausen, M.; Wolff, B. Exceedance of critical loads of nitrogen and sulphur and its relation to forest conditions. Eur. J. For. Res. 2005, 124, 289–300. [Google Scholar] [CrossRef]
  2. Yu, L.; Zanchi, G.; Akselsson, C.; Wallander, H.; Belyazid, S. Modeling the forest phosphorus nutrition in a southwestern Swedish forest site. Ecol. Modell. 2018, 369, 88–100. [Google Scholar] [CrossRef] [Green Version]
  3. Pretzsch, H.; Biber, P.; Schütze, G.; Uhl, E.; Rötzer, T. Forest stand growth dynamics in Central Europe have accelerated since 1870. Nat. Commun. 2014, 5, 4967. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Rossi, S.; Deslauriers, A.; Anfodillo, T.; Morin, H.; Saracino, A.; Motta, R.; Borghetti, M. Conifers in cold environment synchronize maximum growth rate of tree-ring formation with day lenght. New Phytol. 2006, 170, 301–310. [Google Scholar] [CrossRef] [Green Version]
  5. Tomášková, I.; Pokorný, R.; Marek, M.V. Influence of stand density, thinning and elevated CO2 on stem wood density on spruce. J. For. Sci. 2007, 53, 400–405. [Google Scholar] [CrossRef] [Green Version]
  6. Esper, J.; Frank, D.C.; Timonen, M.; Zorita, E.; Wilson, R.J.S.; Luterbacher, J.; Holzkämper, S.; Fischer, N.; Wagner, S.; Nievergelt, D.; et al. Orbital forcing of tree-ring data. Nat. Clim. Change 2012, 2, 862–866. [Google Scholar] [CrossRef]
  7. Samec, P. Analysis of paleoclimatic time series for evaluation of recent climate change on European North Atlantic Region and Bohemian Massif selected localities. Rep. For. Res. 2006, 51, 272–278. [Google Scholar]
  8. Ganopolski, A.; Winhelmann, R.; Schellnhuber, H.J. Critical insolation–CO2 relation for diagnosing past and future glacial inception. Nature 2016, 529, 200–203. [Google Scholar] [CrossRef]
  9. Miyake, F.; Panyushkina, I.P.; Jull, A.J.T.; Adolphi, F.; Brehm, N.; Helama, S.; Kanzawa, K.; Moriya, T.; Muscheler, R.; Nicolussi, K.; et al. A Single-Year Cosmic Ray Event at 5410 BCE Registered in 14C of Tree Rings. Geophys. Res. Lett. 2021, 48, e2021GL093419. [Google Scholar] [CrossRef]
  10. Liu, S.; Bond-Lamberty, B.; Hicke, J.A.; Vargas, R.; Zhao, S.; Chen, J.; Edburg, S.L.; Hu, Y.; Liu, J.; McGuire, A.D.; et al. Simulating the impacts of disturbances on forest carbon cycling in North America: Processes, data, models, and challenges. J. Geophys. Res. 2011, 116, G00K08. [Google Scholar] [CrossRef] [Green Version]
  11. Lee, K.-S.; Hung, D.V.; Kwak, J.-H.; Lim, S.-S.; Lee, K.-H.; Choi, W.-J. Tree Ring Ca/Al as an indicator of Historical Soil Acidification of Pinus densiflora Forest in Southern Korea. Korean J. Environ. Agric. 2011, 30, 229–233. [Google Scholar] [CrossRef]
  12. Chen, L.; Wu, F.-H.; Liu, T.-W.; Chen, J.; Li, Z.-J.; Pei, Z.-M.; Zheng, H.-L. Soil acidity reconstruction nased on tree ring information of a dominant species Abies fabri in the subalpine forest ecosystems in southwest China. Environ. Pollut. 2010, 158, 3219–3224. [Google Scholar] [CrossRef]
  13. Ngao, J.; Epron, D.; Delpierre, N.; Bréda, N.; Granier, A.; Longdoz, B. Spatial variability of soil CO2 efflux linked to soil parameters and ecosystem characteristics in a temperate beech forest. Agric. For. Meteorol. 2012, 154–155, 136–146. [Google Scholar] [CrossRef]
  14. Hansen, M.C.; Potapov, P.V.; Moore, R.; Hancher, M.; Turubanova, S.A.; Tyukavina, A.; Thau, D.; Stehman, S.V.; Goetz, S.J.; Loveland, T.R.; et al. High-resolution Global Maps of 21st-Centrury Forest Cover Change. Science 2013, 342, 850–853. [Google Scholar] [CrossRef] [Green Version]
  15. Vincent, C.; Le Meur, E.; Six, D. Solving the paradox of the end of the Little Ice Age in the Alps. Geophys. Res. Lett. 2005, 32, L09706. [Google Scholar] [CrossRef]
  16. Jones, P.D.; Mann, M.E. Climate Over Past Millennia. Rev. Geophys. 2004, 42, 143. [Google Scholar] [CrossRef] [Green Version]
  17. Mann, M.E.; Jones, P.D. Global surface temperatures over the past two millennia. Geophys. Res. Lett. 2003, 30, 17814. [Google Scholar] [CrossRef] [Green Version]
  18. Briffa, K.R.; Osborn, T.J.; Schweingruber, F.H.; Harris, I.C.; Jones, P.D.; Shiyatov, S.G.; Vaganov, E.A. Low-frequency Temperature Variations from a Nothern Tree Ring Density Network. J. Geophys. Res. 2001, 106, 2629–2941. [Google Scholar] [CrossRef]
  19. Osborn, T.J.; Briffa, K.R. The spatial extent of 20th century warmth in the context of the past 1200 years. Science 2006, 311, 841–844. [Google Scholar] [CrossRef]
  20. Gül, E.; Erşahin, S. Evaluating the desertification vulnerability of a semiarid landscape under different land uses with the environmental sensitivity index. Land Degrad. Dev. 2019, 30, 811–823. [Google Scholar] [CrossRef]
  21. Augustin, S. Linking critical thresholds for acidity to forest condition by using element contents in tree rings: Consequences for the development of an integrated ecosystem assessment. iForest 2009, 2, 51–53. [Google Scholar] [CrossRef]
  22. Srivastava, P.; Singh, R.; Tripathi, S.; Singh, P.; Singh, S.; Singh, H.; Raghubanshi, A.S.; Mishra, P.K. Soil Carbon Dynamics Under Changing Climate—A Research Transition from Absolute to Relative Roles of Inorganic Nitrogen Pools and Associated Microbial Processes: A Review. Pedoshepere 2017, 27, 792–806. [Google Scholar] [CrossRef]
  23. Ďurský, J.; Škvarenina, J.; Mindáš, J.; Miková, A. Regional analysis of climate change impact on Norway spruce (Picea abies L. Karst) growth in Slovak mountain forests. J. For. Sci. 2006, 51, 306–315. [Google Scholar] [CrossRef] [Green Version]
  24. Sigurdsson, B.D.; Medhurst, J.L.; Wallin, G.; Eggertsson, O.; Linder, S. Growth of mature boreal Norway spruce was not affected by elevated [CO2] and/or air temperature unless nutrient availability was improved. Tree Physiol. 2013, 33, 1192–1205. [Google Scholar] [CrossRef] [Green Version]
  25. Rossi, S.; Deslauriers, A.; Anfodillo, T.; Carrer, M. Age-dependent xylogenesis in timberline conifers. New Phytol. 2008, 177, 199–208. [Google Scholar] [CrossRef]
  26. Ziche, D.; Seidling, W. Homogenisation of climate time series from ICP Forests Level II monitoring sites in Germany based on interpolated climate data. Ann. For. Sci. 2010, 67, 804. [Google Scholar] [CrossRef] [Green Version]
  27. Zeng, C.; Yang, L.; Zhu, A.-X.; Rossiter, D.G.; Liu, J.; Liu, J.; Qin, C.; Wang, D. Mapping soil organic matter concentration at different scales using a mixed geographically weighted regression method. Geoderma 2016, 281, 69–82. [Google Scholar] [CrossRef]
  28. Minasny, B.; McBratney, A.B.; Salvador-Blanes, S. Quantitative models for pedogenesis—A review. Geoderma 2008, 144, 140–157. [Google Scholar] [CrossRef]
  29. Debel, A.; Meier, W.J.-H.; Bräuning, A. Climate Signals for Growth Variations of F. sylvatica, P. abies and P. sylvestris in Southeast Germany over the Past 50 Years. Forests 2021, 12, 1433. [Google Scholar] [CrossRef]
  30. Neuhäuslová, Z.; Blažková, D.; Grulich, V.; Husová, M.; Chytrý, M.; Jeník, J.; Jirásek, J.; Kolbek, J.; Kropáč, Z.; Ložek, V.; et al. Map of Potential Natural Vegetation of the Czech Republic /Text part/; Academia: Prague, Czech Republic, 1998; pp. 1–341. [Google Scholar]
  31. Šindelář, J. A preliminary report on autochthonous species composition of forests in the territory of the Czech Republic. Lesnictví For. 1995, 41, 293–299. [Google Scholar]
  32. Zouhar, V. Dřevinná skladba. In Oblastní Plány Rozvoje Lesů 2; ÚHÚL Brandýs nad Labem, Česká Lesnická Společnost, Ministerstvo Zemědělství České Republiky: Prague, Czech Republic, 2018; pp. 20–22. [Google Scholar]
  33. Samec, P.; Caha, J.; Tuček, P.; Zapletal, M.; Cudlín, P.; Kučera, M. Discrimination between acute and chronic decline of Central European forests using map algebra of the growth condition and forest biomass fuzzy sets: A case study. Sci. Total Environ. 2017, 599–600, 899–909. [Google Scholar] [CrossRef]
  34. Mikšovský, J.; Brázdil, R.; Štěpánek, P.; Zahradníček, P.; Pišoft, P. Long-term variability of temperature and precipitation in the Czech Lands: An attribution analysis. Clim. Chang. 2014, 125, 253–264. [Google Scholar] [CrossRef]
  35. Horn, P.S.; Pesce, A.J.; Copeland, B.E. A robust approach to reference interval estimation and evaluation. Clin. Chem. 1998, 44, 622–631. [Google Scholar] [CrossRef] [Green Version]
  36. Peter, M.; Ayer, F.; Cudlín, P.; Egli, S. Belowground ectomycorrhizal communities in three Norway spruce stands with different degrees of decline in the Czech Republic. Mycorrhiza 2008, 18, 157–169. [Google Scholar] [CrossRef] [Green Version]
  37. Neteler, M.; Mitášova, H. Open Source GIS: A GRASS GIS Approach; Kluwer Academic Publishers: Oxford, UK; Springer: Amsterdam, The Netherlands, 2002. [Google Scholar]
  38. Culek, M.; Grulich, V. Biogeographical division. 1:500,000. In Landscape Atlas of the Czech Republic; Hrnčianová, T., Mackovčin, P., Zvara, I., Eds.; Ministry of Environment, The Silva Tarouca Research Institute for Landscape and Ornamental Gardering: Prague, Czech Republic, 2009; pp. 134–135. [Google Scholar]
  39. Schad, P.; van Huysteen, C.; Michéli, E.; Vargas, R. (Eds.) World Reference Base for Soil Resources 2014. International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; World Soil Resources Reports No. 106; FAO: Rome, Italy, 2014. [Google Scholar]
  40. Cools, N.; De Vos, B. Availability and evaluation of European forest soil monitoring data in the study on the effects of air pollution on forests. iForest 2011, 4, 205–211. [Google Scholar] [CrossRef]
  41. Gut, U.; Árvai, M.; Bijak, S.; Camarero, J.; Cedro, A.; Cruz-García, R.; Garamszegi, B.; Hacket-Pain, A.; Hevia, A.; Huang, W.; et al. No systematic effects of sampling direction on climate-growth relationships in a large-scale, multi-species tree-ring data set. Dendrochronologia 2019, 57, 125624. [Google Scholar] [CrossRef]
  42. Mann, M.E.; Bradley, R.S.; Hughes, M.K. Global-Scale Temperature Patterns and Climate Forcing Over the Past Six Centuries. Nature 1998, 392, 779–787. [Google Scholar] [CrossRef]
  43. Feliksik, E.; Wilczyński, S.; Durło, G. Variability of air temperatures of the temperate climatic beltof the Beskid Sądecki Mountains in the period 1971–2000 as exemplified by the climatological station in Kopciowa. Electron. J. Pol. Agric. Univ. 2001, 4, 2. [Google Scholar]
  44. Hadas, P. Analysis of the Climate as a Stress Factor to Forest Ecosystems. Ekologia Bratislava 2000, 19, 162–176. [Google Scholar]
  45. Zapletal, M. Atmospheric deposition of nitrogen and sulphur in relation to critical loads of nitrogen and acidity in the Czech Republic. J. For. Sci. 2006, 52, 92–100. [Google Scholar] [CrossRef] [Green Version]
  46. Oulehle, F.; Hleb, R.; Houška, J.; Šamonil, P.; Hofmeister, J.; Hruška, J. Anthropogenic acidification effects in primeval forests in the Transcarpathian Mts., western Ukraine. Sci. Total Environ. 2010, 408, 856–864. [Google Scholar] [CrossRef]
  47. Keller, K.; Černý, M.; Bytnerowitz, A.; Smith, L.; Šesták, M.; Michalec, M.; Pernegr, V.; Kučera, J. Air Pollution Status of a Representative Site in the Czech Republic Brdy Mountains. Environ. Pollut. 1997, 98, 291–297. [Google Scholar]
  48. Percy, K.E.; Ferretti, M. Air pollution and forest health: Toward new monitoring concepts. Environ. Pollut. 2004, 130, 113–126. [Google Scholar] [CrossRef]
  49. Schröder, W.; Nickel, S.; Jenssen, M.; Riediger, J. Methodology to assess and map the potential development of forest ecosystems exposed to climate change and atmospheric nitrogen deposition: A pilot study in Germany. Sci. Total Environ. 2015, 521–522, 108–122. [Google Scholar] [CrossRef]
  50. Monserud, R.A.; Sterba, H. A basal area increment model for individual tree growing in even- and uneven-aged forest stands in Austria. For. Ecol. Manag. 1996, 80, 57–80. [Google Scholar] [CrossRef]
  51. Korf, V. Analysis and pattern of forest stand increment. Lesnictví 1972, 18, 585–604. [Google Scholar]
  52. Knibbe, B. PAST4—Personal Analysis System for Tree Ring Research Version 4. Instruction Manual; SCIEM/Bernhard Knibbe: Vienna, Austria, 2004; pp. 1–140. [Google Scholar]
  53. Holmes, R.L. Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bull. 1983, 43, 69–78. [Google Scholar]
  54. Meloun, M.; Militký, J.; Hill, M.; Brereton, R.G. Crucial problems in regression modelling and their solutions. Analyst 2002, 127, 433–450. [Google Scholar] [CrossRef]
  55. Chicco, D.; Warrens, M.J.; Jurman, G. The coefficient of determination R-squared is more informative than SMAPE, MAE, MAPE, MSE and RMSE in regression analysis evaluation. Comput. Sci. 2021, 7, e623. [Google Scholar] [CrossRef]
  56. 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 2013, 27, 97–110. [Google Scholar] [CrossRef]
  57. Carrer, M.; Urbinati, C. Long-term change in the sensitivity of tree-ring growth to climate forcing in Larix decidua. New Phytol. 2006, 170, 861–872. [Google Scholar] [CrossRef] [PubMed]
  58. Wilczyński, S.; Feliksik, E. The dendrochronological monitoring of the Western Beskid Mountains (Southern Poland) on the basis of radial increments of Norway spruce (Picea abies /L./ Karst.). Electron. J. Pol. Agric. Univ. 2004, 7, 7. [Google Scholar]
  59. Wang, Z.; Lyu, L.; Liu, W.; Liang, H.; Huang, J.; Zhang, Q.-B. Topographic patterns of forest decline as detected from tree rings and NDVI. Catena 2021, 198, 105011. [Google Scholar] [CrossRef]
  60. Pechanec, V.; Machar, I.; Kilianová, H.; Vyvlečka, P.; Seják, J.; Pokorný, J.; Štěrbová, L.; Prokopová, M.; Cudlín, P. Ranking the Key Forest Habitats in Ecosystem Function Provision: Case Study from Morava River Basin. Forests 2021, 12, 138. [Google Scholar] [CrossRef]
Figure 1. Distribution of selected forest stand types in relation to climatological stations including attribution of time series at investigated altitudinal zones. Forest area ranges: A—Ore Mts.; B—Volcanic Mts.; C—Beroun; D—Central Bohemia; E—Bohemian Forest; F—South-Bohemian basins; G—Bohemian-Moravian Highland; H—Elbe basin; I—North-Bohemian sandstone; J—Western Sudetes; K—Eastern Sudetes; L—Moravian Foothills; M—Moravian Floodplains; N—Outer Western Carpathians.
Figure 1. Distribution of selected forest stand types in relation to climatological stations including attribution of time series at investigated altitudinal zones. Forest area ranges: A—Ore Mts.; B—Volcanic Mts.; C—Beroun; D—Central Bohemia; E—Bohemian Forest; F—South-Bohemian basins; G—Bohemian-Moravian Highland; H—Elbe basin; I—North-Bohemian sandstone; J—Western Sudetes; K—Eastern Sudetes; L—Moravian Foothills; M—Moravian Floodplains; N—Outer Western Carpathians.
Forests 13 00588 g001
Figure 2. Development of year average temperatures (T), annual precipitation (P) and total atmospheric pollution concentrations at altitudinal zones covering sampled forest stand types.
Figure 2. Development of year average temperatures (T), annual precipitation (P) and total atmospheric pollution concentrations at altitudinal zones covering sampled forest stand types.
Forests 13 00588 g002
Figure 3. Development of soil series characteristics clay content, pH/H2O, base saturation (BS), organic carbon (Corg) and total A2O3 and CaO contents.
Figure 3. Development of soil series characteristics clay content, pH/H2O, base saturation (BS), organic carbon (Corg) and total A2O3 and CaO contents.
Forests 13 00588 g003
Figure 4. Current basal-area increment of forest stands selected on various soil series.
Figure 4. Current basal-area increment of forest stands selected on various soil series.
Forests 13 00588 g004
Table 1. Proportion of forest type management populations composed of relief types and soil series in forest area ranges of the Czech Republic (%). For an explanation of forest area ranges locations see Figure 1.
Table 1. Proportion of forest type management populations composed of relief types and soil series in forest area ranges of the Czech Republic (%). For an explanation of forest area ranges locations see Figure 1.
ReliefSeriesABCDEFGHIJKLMNCR
FloodplainRiverine0.20.10.10.70.11.10.04.20.40.11.00.544.12.01.5
Swamp1.11.61.10.90.71.01.00.90.50.61.30.50.90.00.8
RockyPine0.70.47.80.40.117.70.215.340.31.2-0.40.4-4.0
UplandExposed0.17.85.14.20.10.30.02.42.40.10.04.40.70.31.8
Acidic1.13.414.310.90.43.20.121.75.20.70.411.530.30.55.7
Nutrient-rich0.016.44.54.10.01.20.034.72.41.80.36.121.09.54.8
Gleyed0.60.02.90.70.110.40.217.70.30.10.10.20.80.71.6
KarsticExposed0.23.30.10.70.10.00.2---0.51.2-0.10.4
Nutrient-rich---0.00.10.00.0---0.00.6--0.1
Wet0.0-0.00.00.06.60.1-0.5--0.0--0.2
Foothill/highlandRavine3.04.62.51.93.03.40.50.66.82.31.22.11.10.62.1
Exposed1.010.44.45.31.30.81.60.02.93.16.18.00.04.43.7
Acidic7.15.722.531.15.89.17.21.414.714.72.712.7-1.310.7
Nutrient-rich0.228.013.424.23.95.46.50.54.18.634.740.50.744.217.9
Gleyed0.33.511.214.21.422.73.40.12.41.32.32.20.12.04.4
MontaneExposed9.93.10.90.07.0-3.5-1.08.66.50.7-9.74.4
Acidic28.10.22.20.029.00.029.7-9.521.95.62.8-0.411.5
Nutrient-rich10.69.10.50.212.80.122.7-3.715.225.03.2-23.511.1
Gleyed7.62.05.40.114.50.120.2-1.63.94.12.0-0.45.9
Wet0.50.30.90.42.117.01.90.31.41.00.60.3-0.11.3
High-mountainExposed0.7---0.9----2.11.3--0.10.5
Acidic15.2---5.6----5.31.7--0.02.2
Nutrient-rich0.4---1.1----0.11.6--0.10.3
Gleyed1.9---3.9-0.1--0.70.3--0.00.7
Wet9.1-0.2-3.7-0.9--1.70.1--0.01.3
Timberline0.4---2.1----3.52.2--0.00.8
Supalpine0.0--------1.50.3---0.1
Table 2. Air condition characteristics of sampled forest stand types (mean ± standard deviation).
Table 2. Air condition characteristics of sampled forest stand types (mean ± standard deviation).
GeotectonicsReliefSeriesTree-SpeciesOriginAltitudeTPSO2NOxO3
Bohemian MassifRockyPinePinus sylvestrisNatural280–3307.52 ± 0.50687 ± 506.67 ± 1.4321.03 ± 1.2054.42 ± 2.07
Pinus sylvestrisPlanted265–2967.59 ± 0.21648 ± 36.72 ± 1.4621.43 ± 1.2253.65 ± 2.02
UplandFlatPicea abiesPlanted255–2938.29 ± 0.13690 ± 106.75 ± 1.4821.74 ± 1.2553.03 ± 1.99
ExposedQuercus petraeaNatural516–5887.04 ± 0.16582 ± 146.19 ± 1.1817.04 ± 1.0062.26 ± 2.74
Pinus sylvestrisPlanted245–2678.33 ± 0.05559 ± 446.77 ± 1.4821.83 ± 1.2552.86 ± 1.98
Picea abiesPlanted457–5037.17 ± 0.28643 ± 356.33 ± 1.2418.20 ± 1.0459.97 ± 2.52
Nutrient-richQuercus petraeaNatural386–3987.69 ± 0.04550 ± 46.50 ± 1.3319.63 ± 1.1157.18 ± 2.27
Pinus sylvestrisPlanted197–2628.56 ± 0.02578 ± 56.82 ± 1.5222.26 ± 1.2952.02 ± 1.94
Picea abiesPlanted236–2538.47 ± 0.01572 ± 96.79 ± 1.5022.02 ± 1.2752.49 ± 1.96
GleydedQuercus roburNatural195–2008.58 ± 0.01559 ± 26.88 ± 1.5622.78 ± 1.3351.00 ± 1.89
Pinus sylvestrisPlanted192–2038.59 ± 0.06556 ± 76.88 ± 1.5622.78 ± 1.3351.00 ± 1.89
Picea abiesPlanted277–3198.33 ± 0.29570 ± 66.68 ± 1.4421.15 ± 1.2054.20 ± 2.06
MontaneFoothillPicea abiesPlanted500–7046.49 ± 0.61940 ± 656.15 ± 1.1716.70 ± 1.0062.91 ± 2.82
RavinePinus sylvestrisNatural321–3716.88 ± 0.63624 ± 446.36 ± 1.2718.46 ± 1.0759.47 ± 2.50
ExposedPicea abiesPlanted805–8405.05 ± 0.31976 ± 635.68 ± 1.0212.84 ± 1.0170.49 ± 3.65
Nutrient-richPicea abiesPlanted741–8205.17 ± 0.43963 ± 565.72 ± 1.0313.12 ± 1.0069.93 ± 3.58
GleyedPicea abiesPlanted800–8404.99 ± 0.16873 ± 55.66 ± 1.0112.70 ± 1.0170.76 ± 3.68
WetPicea abiesPlanted849–8604.82 ± 0.04871 ± 25.60 ± 1.0012.14 ± 1.0371.86 ± 3.81
TimberlinePicea abiesNatural1046–11233.51 ± 0.341005 ± 215.15 ± 1.028.42 ± 1.2679.16 ± 4.70
CarpathiansFloodplainSwampSalix fragilisNatural235–3357.90 ± 0.40677 ± 136.71 ± 1.4521.36 ± 1.2253.79 ± 2.04
RiverineFraxinus excelsiorNatural223–3137.96 ± 0.37796 ± 1456.68 ± 1.4421.14 ± 1.2154.22 ± 2.07
MontaneFoothillPicea abiesPlanted581–6266.61 ± 0.331107 ± 196.17 ± 1.1716.88 ± 0.9962.57 ± 2.77
RavinePicea abiesPlanted1050–11503.71 ± 0.441281 ± 885.12 ± 1.038.17 ± 1.2879.65 ± 4.76
ExposedPicea abiesPlanted559–8115.83 ± 0.621160 ± 775.89 ± 1.0714.54 ± 0.9967.16 ± 3.27
Nutrient-richPicea abiesPlanted697–10625.31 ± 1.021191 ± 825.65 ± 1.0412.61 ± 1.0570.93 ± 3.71
TimberlinePicea abiesNatural1136–11733.59 ± 0.191209 ± 115.01 ± 1.057.28 ± 1.3481.39 ± 4.98
T—temperature (°C); P—precipitation (mm); SO2—sulphur dioxide (g/m3); NOx—nitrogen oxides (g/m3); O3—tropospheric ozone (g/m3).
Table 3. Soil characteristics of sampled forest stand types in relation to bedrock type (mean ± standard deviation).
Table 3. Soil characteristics of sampled forest stand types in relation to bedrock type (mean ± standard deviation).
GeotectonicsReliefSeriesTree-SpeciesOriginBedrockSoil GroupclaypHBSCorgAl2O3CaO
Bohemian MassifRockyPinePinus sylvestrisNaturalsandstoneArenic Podzol1.17 ± 0.664.19 ± 0.2525.46 ± 14.400.75 ± 0.580.07 ± 0.030.02 ± 0.01
Pinus sylvestrisPlantedsandstoneArenic Podzol2.05 ± 1.024.07 ± 0.2619.14 ± 7.700.83 ± 0.690.13 ± 0.070.02 ± 0.03
UplandFlatPicea abiesPlantedloessHaplic Luvisol2.73 ± 2.003.94 ± 0.2412.65 ± 6.864.08 ± 1.381.75 ± 0.150.11 ± 0.02
ExposedQuercus petraeaNaturalphonoliteHaplic Cambisol17.97 ± 11.285.27 ± 0.6570.47 ± 32.885.61 ± 1.814.10 ± 0.332.09 ± 0.25
Pinus sylvestrisPlantedsandstoneLuvic Cambisol15.44 ± 10.545.85 ± 1.7560.98 ± 42.311.22 ± 0.321.27 ± 0.391.21 ± 1.92
Picea abiesPlantedgraniteHaplic Cambisol12.79 ± 26.224.22 ± 0.2218.29 ± 5.562.46 ± 0.632.99 ± 0.360.18 ± 0.11
Nutrient-richQuercus petraeaNaturalphonoliteHaplic Luvisol31.47 ± 20.885.48 ± 0.4592.56 ± 7.096.46 ± 2.955.60 ± 1.610.80 ± 0.22
Pinus sylvestrisPlantedphonoliteEpidystric Stagnosol2.13 ± 0.354.34 ± 0.6032.78 ± 19.622.65 ± 0.831.03 ± 0.130.06 ± 0.04
Picea abiesPlantedmarlstoneLuvic Cambisol19.1 ± 15.553.87 ± 0.1120.63 ± 3.962.95 ± 2.431.75 ± 1.330.04 ± 0.02
GleydedQuercus roburNaturalgravel-sandEndogleyic Stagnosol9.85 ± 10.474.35 ± 0.3322.2 ± 21.231.75 ± 0.440.72 ± 0.240.03 ± 0.03
Pinus sylvestrisPlantedgravel-sandEntic Podzol4.86 ± 4.754.46 ± 0.9434.07 ± 32.611.21 ± 1.060.61 ± 0.200.05 ± 0.05
Picea abiesPlantedgravel-sandStagni-luvic Cambisol19.62 ± 23.463.85 ± 0.1011.37 ± 3.414.43 ± 2.352.02 ± 1.920.05 ± 0.05
MontaneFoothillPicea abiesPlantedamphiboliteHaplic Cambisol2.41 ± 3.224.40 ± 0.6741.09 ± 32.079.81 ± 7.174.16 ± 0.460.37 ± 0.14
RavinePinus sylvestrisNaturalsandstoneArenic Podzol6.38 ± 8.093.72 ± 0.0916.62 ± 6.771.67 ± 0.700.33 ± 0.110.02 ± 0.01
ExposedPicea abiesPlantedorthogneissSkeletic Cambisol1.10 ± 0.543.91 ± 0.1416.44 ± 2.968.96 ± 4.063.35 ± 0.540.13 ± 0.06
Nutrient-richPicea abiesPlantedphylliteEntic Podzol0.98 ± 0.404.04 ± 0.3117.65 ± 10.17.15 ± 2.183.28 ± 1.240.15 ± 0.04
GleyedPicea abiesPlantedparagneissStagnic Cambisol1.68 ± 0.733.97 ± 0.3212.05 ± 1.335.85 ± 1.813.62 ± 1.190.16 ± 0.06
WetPicea abiesPlantedparagneissHaplic Gleysol2.53 ± 0.634.02 ± 0.2015.5 ± 6.416.08 ± 1.944.15 ± 1.310.14 ± 0.07
TimberlinePicea abiesNaturalorthogneissHaplic Podzol2.60 ± 0.903.86 ± 0.1816.56 ± 5.89.19 ± 8.111.25 ± 0.490.07 ± 0.02
CarpathiansFloodplainSwampSalix fragilisNaturalfluvial loamStagnic Fluvisol16.50 ± 2.156.10 ± 0.0397.72 ± 2.474.60 ± 1.462.29 ± 0.330.25 ± 0.05
RiverineFraxinus excelsiorNaturalfluvial loamHaplic Fluvisol3.18 ± 2.345.34 ± 1.0666.77 ± 36.384.07 ± 1.353.18 ± 1.070.39 ± 0.32
MontaneFoothillPicea abiesPlantedsandstoneHaplic Cambisol1.29 ± 0.323.78 ± 0.0811.34 ± 4.5210.64 ± 5.612.68 ± 0.290.08 ± 0.03
RavinePicea abiesPlantedsandstoneUmbric Cambisol15.66 ± 8.113.95 ± 0.135.40 ± 1.144.40 ± 1.012.55 ± 0.420.01 ± 0.00
ExposedPicea abiesPlantedsandstoneSkeletic Cambisol1.24 ± 0.723.81 ± 0.1715.63 ± 7.0813.75 ± 3.462.65 ± 0.380.10 ± 0.05
Nutrient-richPicea abiesPlantedsandstoneSkeletic Cambisol1.59 ± 0.563.67 ± 0.0711.63 ± 3.6016.23 ± 5.822.45 ± 0.400.06 ± 0.01
TimberlinePicea abiesNaturalsandstoneHaplic Podzol1.22 ± 0.743.72 ± 0.1017.85 ± 6.4617.02 ± 14.251.38 ± 0.460.07 ± 0.05
Clay—soil clay (<2 μM) content (%); BS—base saturation (%); Corg—organic carbon (%); Al2O3—aluminium oxide (%); CaO—calcium oxide (%).
Table 4. Grow characteristics of sampled forest stand types (mean ± standard deviation).
Table 4. Grow characteristics of sampled forest stand types (mean ± standard deviation).
GeotectonicsReliefSeriesTree-SpeciesOriginAgeVhahd1.3BAI
Bohemian MassifRockyPinePinus sylvestrisNatural115 ± 2296 ± 423 ± 124.00 ± 2.074.48 ± 1.62
Pinus sylvestrisPlanted89 ± 3225 ± 419 ± 422.00 ± 2.023.79 ± 1.53
UplandFlatPicea abiesPlanted66 ± 13278 ± 5421 ± 1023.00 ± 1.997.45 ± 4.4
ExposedQuercus petraeaNatural77 ± 1151 ± 119 ± 122.50 ± 2.743.17 ± 1.30
Pinus sylvestrisPlanted103 ± 1256 ± 221 ± 425.00 ± 1.983.72 ± 1.37
Picea abiesPlanted109 ± 1342 ± 2623 ± 124.50 ± 2.525.28 ± 1.62
Nutrient-richQuercus petraeaNatural100 ± 2203 ± 2519 ± 327.00 ± 2.272.36 ± 0.92
Pinus sylvestrisPlanted86 ± 3328 ± 424 ± 327.00 ± 1.944.31 ± 1.28
Picea abiesPlanted85 ± 4446 ± 625 ± 527.50 ± 1.965.94 ± 1.81
GleydedQuercus roburNatural97 ± 0241 ± 523 ± 228.50 ± 1.894.08 ± 1.20
Pinus sylvestrisPlanted86 ± 1334 ± 1123 ± 027.00 ± 1.894.97 ± 1.28
Picea abiesPlanted84 ± 1376 ± 1023 ± 325.00 ± 2.065.48 ± 2.21
MontaneFoothillPicea abiesPlanted75 ± 1365 ± 3024 ± 529.50 ± 2.828.18 ± 6.2
RavinePinus sylvestrisNatural123 ± 0256 ± 622 ± 428.50 ± 2.503.39 ± 2.38
ExposedPicea abiesPlanted96 ± 3505 ± 1628 ± 636.50 ± 3.655.89 ± 1.74
Nutrient-richPicea abiesPlanted76 ± 8458 ± 3726 ± 1430.00 ± 3.586.55 ± 2.04
GleyedPicea abiesPlanted78 ± 6391 ± 3527 ± 933.00 ± 3.688.67 ± 7.56
WetPicea abiesPlanted63 ± 8252 ± 2220 ± 1123.50 ± 3.816.52 ± 2.16
TimberlinePicea abiesNatural90 ± 6346 ± 1521 ± 824.50 ± 4.704.03 ± 1.24
CarpathiansFloodplainSwampSalix fragilisNatural45 ± 2220 ± 1014 ± 836.40 ± 2.042.02 ± 1.76
RiverineFraxinus excelsiorNatural74 ± 3402 ± 1926 ± 631.50 ± 2.078.47 ± 2.87
MontaneFoothillPicea abiesPlanted79 ± 5227 ± 1228 ± 239.06 ± 2.777.13 ± 5.81
RavinePicea abiesPlanted72 ± 4216 ± 1318 ± 126.50 ± 4.763.45 ± 1.66
ExposedPicea abiesPlanted89 ± 4555 ± 2230 ± 833.50 ± 3.277.14 ± 1.86
Nutrient-richPicea abiesPlanted79 ± 4462 ± 1526 ± 132.50 ± 3.717.32 ± 2.63
TimberlinePicea abiesNatural86 ± 6240 ± 2020 ± 827.50 ± 4.983.25 ± 1.33
Vha—volume stock (m3/ha); h—height (m); d1.3—diameter in breast height (cm); BAI—basal-area increment (m2/ha).
Table 5. Linear correlations and analysis of residues between tree current basal-area increment and properties of growth conditions at sampled forest stand type. Bold statistically significant correlation at p < 0.05, normal less significant between p ≥ 0.05 and p < 0.50 and light insignificant at p ≥ 0.50. RSD—residual standard deviation; ER—residual elevation; AR—residual skewness (for detailed information about properties of growth conditions see Table 2 and Table 3).
Table 5. Linear correlations and analysis of residues between tree current basal-area increment and properties of growth conditions at sampled forest stand type. Bold statistically significant correlation at p < 0.05, normal less significant between p ≥ 0.05 and p < 0.50 and light insignificant at p ≥ 0.50. RSD—residual standard deviation; ER—residual elevation; AR—residual skewness (for detailed information about properties of growth conditions see Table 2 and Table 3).
GeotectonicsReliefSeriesTree-SpeciesOriginTPSO2NOxO3claypHBSCorgAl2O3CaORSDERAR
Bohemian MassifRockyPinePinus sylvestrisNatural−0.090.020.390.08−0.030.290.16−0.08−0.290.18−0.080.24−0.14−0.72
Pinus sylvestrisPlanted0.12−0.120.290.20−0.280.10−0.08−0.10−0.130.13−0.020.24−0.020.39
UplandFlatPicea abiesPlanted−0.010.230.07−0.210.150.240.240.22−0.14−0.010.121.37−1.440.05
ExposedQuercus petraeaNatural0.10−0.17−0.030.01−0.110.00−0.32−0.100.060.04−0.280.301.811.29
Pinus sylvestrisPlanted−0.280.070.090.23−0.070.030.070.05−0.040.16−0.030.21−0.47−0.26
Picea abiesPlanted−0.260.170.13−0.410.19−0.050.02−0.11−0.06−0.01−0.088.51−0.50−0.28
Nutrient−richQuercus petraeaNatural0.110.09−0.18−0.310.32−0.16−0.36−0.100.06−0.09−0.180.221.220.88
Pinus sylvestrisPlanted−0.36−0.110.010.51−0.340.350.17−0.04−0.150.180.120.24−0.13−0.18
Picea abiesPlanted−0.05−0.05−0.120.27−0.150.270.17−0.01−0.080.000.110.43−0.310.16
GleydedQuercus roburNatural−0.09−0.290.380.34−0.290.04−0.12−0.270.06−0.13−0.240.21−0.760.22
Pinus sylvestrisPlanted0.00−0.01−0.230.19−0.070.390.23−0.17−0.090.300.020.19−0.260.32
Picea abiesPlanted−0.12−0.120.100.41−0.350.130.12−0.120.090.00−0.130.450.700.99
MontaneFoothillPicea abiesPlanted0.51−0.300.100.10−0.67−0.090.030.000.070.00−0.010.98−0.640.45
RavinePinus sylvestrisNatural0.070.100.32−0.03−0.060.100.000.02−0.030.130.100.170.430.03
ExposedPicea abiesPlanted−0.03−0.01−0.300.040.00−0.060.050.100.13−0.140.030.281.360.18
Nutrient−richPicea abiesPlanted0.29−0.260.16−0.070.550.08−0.18−0.10−0.010.160.140.894.381.28
GleyedPicea abiesPlanted−0.22−0.120.270.40−0.45−0.250.210.33−0.070.120.000.921.60−0.81
WetPicea abiesPlanted−0.06−0.140.180.31−0.29−0.170.380.48−0.130.24−0.110.54−0.81−0.22
TimberlinePicea abiesNatural0.02−0.430.060.28−0.380.400.330.160.33−0.19−0.010.42−0.79−0.44
CarpathiansFloodplainSwampSalix fragilisNatural0.140.20−0.06−0.230.03−0.24−0.030.060.03−0.24−0.041.033.31−0.02
RiverineFraxinus excelsiorNatural−0.32−0.150.240.020.190.310.16−0.19−0.230.25−0.020.712.181.64
MontaneFoothillPicea abiesPlanted−0.29−0.280.200.03−0.01−0.210.05−0.10−0.180.310.010.524.062.00
RavinePicea abiesPlanted0.27−0.13−0.36−0.09−0.140.190.370.38−0.020.420.350.251.371.33
ExposedPicea abiesPlanted0.42−0.40−0.01−0.46−0.390.070.300.14−0.060.220.250.682.81−0.03
Nutrient−richPicea abiesPlanted0.14−0.490.130.00−0.430.030.430.14−0.110.320.330.330.160.45
TimberlinePicea abiesNatural0.03−0.08−0.170.080.050.100.260.32−0.010.350.360.382.351.46
Table 6. Multiple regressions between current basal-area increment and properties of growth conditions at forest stand types (for detail information about properties of growth conditions see Table 2 and Table 3).
Table 6. Multiple regressions between current basal-area increment and properties of growth conditions at forest stand types (for detail information about properties of growth conditions see Table 2 and Table 3).
GeotectonicsReliefSeriesTree-SpeciesOriginFminr2TPSO2NOxO3claypHBSCorgAl2O3CaOb
Bohemian MassifRockyPinePinus sylvestrisNatural5.070.61–0.850.341−0.1790.552−0.587−0.5460.248−0.0490.469−0.0410.300−0.5214.406
Pinus sylvestrisPlanted1.580.40–0.450.1400.1740.544−0.081−0.059−0.021−0.025−0.080−0.0500.1550.2723.017
UplandFlatPicea abiesPlanted0.570.20–0.640.3601.9745.069−2.4221.813−0.537−0.552−0.153−1.830−1.4980.38510.233
ExposedQuercus petraeaNatural0.830.26–0.55−0.1480.0280.118−0.1090.227−0.018−0.2420.131−0.0300.554−0.5493.097
Pinus sylvestrisPlanted1.860.44–0.690.060−0.050−0.3951.4680.104−0.164−0.0970.248−0.1500.479−0.6424.400
Picea abiesPlanted0.440.16–0.43−0.058−0.209−0.105−0.249−0.055−0.141−0.6720.862−0.064−0.045−0.1604.625
Nutrient−richQuercus petraeaNatural1.590.40–0.70−0.0790.198−0.7040.1990.134−0.347−0.5200.514−0.293−0.158−0.2932.814
Pinus sylvestrisPlanted4.130.64–0.71−0.133−0.022−0.5261.1160.0900.006−0.077−0.0400.0470.383−0.4574.865
Picea abiesPlanted2.060.47–0.580.590−0.0010.507−0.710−1.229−0.717−0.2780.718−0.1490.595−0.7236.887
GleydedQuercus roburNatural2.910.55–0.640.136−0.1610.2380.229−0.148−0.272−0.2350.072−0.0300.177−0.1024.125
Pinus sylvestrisPlanted3.520.47–0.550.365−0.026−0.3230.178−0.111−0.473−0.4050.1170.4000.2360.3645.810
Picea abiesPlanted2.960.40–0.48−0.304−0.441−0.061−0.944−0.304−0.770−1.4580.4310.2251.4120.7064.887
MontaneFoothillPicea abiesPlanted1.170.33–0.50−0.2060.2430.154−0.834−0.544−1.233−1.3500.0370.0450.921−0.0889.057
RavinePinus sylvestrisNatural3.040.43–0.52−0.0250.1520.4090.351−0.0630.911−0.181−0.5450.2990.1790.4422.858
ExposedPicea abiesPlanted4.230.56–0.720.382−0.223−0.090−0.1310.457−0.447−0.5560.025−0.682−0.7170.6106.117
Nutrient-richPicea abiesPlanted0.980.23–0.42−0.301−0.444−0.292−0.081−0.653−0.603−0.6280.508−0.275−0.4410.7916.777
GleyedPicea abiesPlanted3.990.55–0.852.567−0.051−8.8860.066−1.3050.472−25.07123.376−14.9384.9434.65510.945
WetPicea abiesPlanted6.450.70–0.761.278−3.053−14.22419.04621.6294.129−18.89912.463−18.37719.398−3.73819.110
TimberlinePicea abiesNatural11.600.78–0.840.019−0.141−0.583−0.6310.795−0.381−0.3920.405−0.597−0.2750.0463.875
CarpathiansFloodplainSwampSalix fragilisNatural2.690.53–0.600.245−0.1851.0480.3360.0742.5980.6830.0001.107−1.858−0.3631.758
RiverineFraxinus excelsiorNatural1.630.41–0.77−0.5041.055−0.8600.713−0.024−2.3971.598−0.888−1.6020.808−0.5188.951
MontaneFoothillPicea abiesPlanted4.110.59–0.662.635−0.0030.2261.2510.0250.562−18.6710.0330.5850.0730.64041.357
RavinePicea abiesPlanted14.890.82–0.92−0.3031.2540.2430.844−0.1480.0810.305−0.0710.7980.652−0.7043.698
ExposedPicea abiesPlanted8.710.73–0.84−0.0150.116−0.862−0.2920.122−0.2900.1320.021−0.2300.303−0.4027.097
Nutrient−richPicea abiesPlanted5.860.64–0.70−0.6440.7930.976−1.0670.3260.591−0.361−0.7930.2450.850−0.0278.381
TimberlinePicea abiesNatural6.120.72–0.91−0.032−0.197−0.4620.113−0.2000.2090.6040.689−0.872−0.7480.8644.013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Samec, P.; Rychtecká, P.; Zeman, M.; Zapletal, M. Environmental Effects among Differently Located and Fertile Sites on Forest Basal-Area Increment in Temperate Zone. Forests 2022, 13, 588. https://doi.org/10.3390/f13040588

AMA Style

Samec P, Rychtecká P, Zeman M, Zapletal M. Environmental Effects among Differently Located and Fertile Sites on Forest Basal-Area Increment in Temperate Zone. Forests. 2022; 13(4):588. https://doi.org/10.3390/f13040588

Chicago/Turabian Style

Samec, Pavel, Petra Rychtecká, Miroslav Zeman, and Miloš Zapletal. 2022. "Environmental Effects among Differently Located and Fertile Sites on Forest Basal-Area Increment in Temperate Zone" Forests 13, no. 4: 588. https://doi.org/10.3390/f13040588

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