Next Article in Journal
Characterization of Water in Wood by Time-Domain Nuclear Magnetic Resonance Spectroscopy (TD-NMR): A Review
Next Article in Special Issue
Applications of TLS and ALS in Evaluating Forest Ecosystem Services: A Southern Carpathians Case Study
Previous Article in Journal
Urban Forest Sustainability in Residential Areas in the City of Santo Domingo
Previous Article in Special Issue
SWAT Model Adaptability to a Small Mountainous Forested Watershed in Central Romania
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Past and Future of Temperate Forests State under Climate Change Effects in the Romanian Southern Carpathians

1
National Institute for Research and Development in Forestry “Marin Drăcea”, 077191 Voluntari, Romania
2
Faculty of Silviculture and Forest Engineering, “Transilvania” University of Brașov, 500123 Brașov, Romania
*
Author to whom correspondence should be addressed.
Forests 2021, 12(7), 885; https://doi.org/10.3390/f12070885
Submission received: 21 May 2021 / Revised: 25 June 2021 / Accepted: 4 July 2021 / Published: 7 July 2021
(This article belongs to the Special Issue Climate Change and Air Pollution Effects on Forest Ecosystems)

Abstract

:
Research Highlights: Carpathian forests hold high ecological and economic value while generating conservation concerns, with some of these forests being among the few remaining temperate virgin forests in Europe. Carpathian forests partially lost their original integrity due to their management. Climate change has also gradually contributed to forest changes due to its modification of the environmental conditions. Background and Objectives: Understanding trees’ responses to past climates and forms of management is critical in foreseeing the responses of forests to future conditions. This study aims (1) to determine the sensitivity of Carpathian forests to past climates using dendrochronological records and (2) to describe the effects that climate change and management will have on the attributes of Carpathian forests, with a particular focus on the different response of pure and mixed forests. Materials and Methods: To this end, we first analysed the past climate-induced growth change in a dendrochronological reference series generated for virgin forests in the Romanian Curvature Carpathians and then used the obtained information to calibrate spatially explicit forest Landis-II models for the same region. The model was used to project forest change under four climate change scenarios, from mild to extreme. Results: The dendrochronological analysis revealed a climate-driven increase in forest growth over time. Landis-II model simulations also indicate that the amount of aboveground forest biomass will tend to increase with climate change. Conclusions: There are differences in the response of pure and mixed forests. Therefore, suitable forest management is required when forests change with the climate.

1. Introduction

Forests have a specific structure and function that determine their capacity to provide multiple ecosystem services, including C (carbon) sequestration [1,2], which has a main role in atmospheric CO2 (carbon dioxide) balance and, therefore, climate change mitigation strategies and agreements. Environmental conditions, particularly with regard to the soil and climate, play a significant role in forest structure development, forest succession, productivity, etc. Since the climate is continuously changing [3,4], forests also do [5,6]. Even when climate extremes do not have strong effects at the local level, forest changes are forced by direct [6] and indirect [7] climate change effects [8].
However, in terms of structure and composition, forests are often characterized by a high diversity and spatial heterogeneity [9,10]. Carpathian temperate forests, for instance, host both pure and mixed altitudinal forests [11,12], which have very different ecological preferences and properties [11,13,14]. Mixed stands proved to have high resistance and resilience and, often, a high level of productivity, compared to pure stands [15,16,17,18,19]. On the other hand, even-aged monospecific stands instead are less complex functionally and structurally complex [20]. Functional diversity, particularly in mixed forests, provides redundancy [21,22], resistance [15,16,22] and resilience [16,17] to different disturbances. Having a relatively high forest integrity compared to other temperate regions in the world [10,18,19], the mixed forests in the Carpathians are usually constituted by uneven-aged stands with high structural and functional diversity [9,23,24].
Characterized by their high stability and adaptability to climate change and stress factors [25] and by their high productivity [26], the mixed forest management model of the Carpathians might be a useful reference for other managed forests [23,27] in the region. Forest management, including its effect on forest structure and functioning, can be used to promote long-term forest resilience and resistance [22,28,29], contributing to C sequestration, as part of the atmospheric CO2 balance and climate change mitigation and adaptation strategies and specific programs [22,30]. However, a careful assessment of the forest’s state is required.
Even though mixed stands benefit from functional diversity [31], they are also sensitive to management and changes in environmental conditions [31]. Proven drivers of mixed forest change include temperature rises [16] and drought [17]. Functional diversity guarantees greater stability at the stand and forest levels and can increase stress resilience and productivity under drought events [16], but climate change can trigger modifications in species abundance [16,32,33] and can, therefore, alter the structure and functioning of mixed forests, as well as the services they provide [1]. There is still high uncertainty regarding the local impacts of climate change and their potential changes in temperate forests [30], either pure or mixed, as well as their capacity for stocking biomass under regular management. There is also the need to reconstruct and to reanalyse past changes [8] in the context of climate change effects in order to explore future changes.
The complexity of the ecosystem processes and the need for reliable forecasts of climate change mitigation and adaptation actions make the use of a model approach that integrates plant–soil–climate relationships necessary. Landis-II models have been proven to be robust in simulating ecophysiological processes at the cohort level, including natural disturbances and management at the landscape spatial scale [34]. Landis-II models, coupled with the PnET model [35], cover local climate effects on vegetation development [35,36], and they are suitable tools for climate change forecasting [36]. However, their implementation is complex, because many parameters are intercorrelated, and many data inputs must be measured locally [36]. In addition, their calibration is complex and potentially affected by operator decisions [36]. Error propagation can affect the entire model, and the calibrated model can rarely be exported to other contexts, either geographical or temporal. Therefore, new alternatives for models’ implementation are required, and this may be achieved by integrating ecological models with machine learning (ML) tools [37,38,39] and by utilizing the best data available, such as those from long-term ecological research databases [24], e.g., dendrochronological records.
It is well known that tree ring width is significantly influenced by climate variation [40], and therefore, tree ring cores provide a historical record of the climate variability (dendroclimatology) [41], forest growth and age structure, disturbances, site conditions (dendroecology) [42], etc. Climate reconstruction is commonly achieved through the study of the climate signals found in the tree rings [40,41,42,43] of open-grown trees [14,40,44], because they are more exposed to climate events and because isolated trees are less affected by competition and other ecological processes [14,45]. Non-isolated trees, mainly in virgin mixed forests, are the best option for tackling intra- and inter-specific interactions. Carpathian virgin forests, internationally recognized for their conservation value, for research and for the services they provide [22,25,46], together with other old-grown forests in the Curvature Carpathians, give an opportunity to deepen our understanding of the roles of mixed vs. pure forests role and their response to long-term processes [23,40,41,42,43], such as climate change [8,44,47], as well as their effects on biodiversity, productivity and C cycling [9,22,29]. Understanding species interactions and responses to past climates and forms of management is critical in foreseeing forest responses to future conditions and in creating optimal strategies for climate change mitigation and adaptation. Within this context, this study aims (1) to determine the sensitivity of Carpathian forests to past climates using dendrochronological records and (2) to describe the future effects of forest management and climate change on Carpathian Forest biomass, with a particular focus on the different response of pure and mixed forests.

2. Materials and Methods

2.1. Study Area

The research area is located in the Penteleu Mountains within the Penteleu Forest District (Figure 1). It is a very representative area of Eastern Europe temperate forests specific to the Romanian Curvature Carpathians, at altitudes between 400 and 1800 m.a.s.l. The average annual temperature ranges between 4 °C and 6 °C. The minimum and maximum absolute temperatures recorded are −33.5 °C and +38 °C. The frost period ranges from 120 days at 800 m.a.s.l. to 220 days at the highest peaks. The average annual precipitation is 830 mm, mainly concentrated in summertime. The prevailing winds come from the NW, W and NE directions, with a 6–7 m∙s−1 annual average speed in high-altitude areas and 2.5–3.5 m∙s−1 in low-altitude areas.
Most of the forest stands included in the research area are temperate montane mixed conifer–broadleaved forests. The main species in terms of timber volume are Fagus sylvatica L. (beech—34.6%), Picea abies (L.) H. Karst. (Norway spruce—38.1%) and Abies alba Mill. (silver fir—22.6%). Other deciduous and softwood species provide less than 4.7% of the overall volume [48]. There is no historical evidence that extreme events (such as windthrows, drought, insect damage, fires, etc.) had a relevant role in the development of the Penteleu forest over time.
Penteleu forest, similar to most of the Romanian Carpathian forests, suffered the effects of profound changes in its management. Being initially private property owned by local landowners and local communities, logging at the beginning of the 20th century was carried out at the pleasure of the owners and in accordance with market demand. However, forest integrity was not heavily affected, regardless of the massive harvests that occurred around a forest railway. Within this period, conifer stands were clear-cut and artificially regenerated afterwards. Mixed stands were subjected to selective cutting, extracting only conifers with a diameter greater than 30 cm, while beech was not extracted for commercial purposes, because it was considered, at that time, a species of low economic value. The Romanian forest nationalization of 1948 was an important step toward obtaining quality timber in maximum quantities, while, at the same time, satisfying the protection functions. Later, in 1975, the forest protection status generated by a change in the forest legislation (especially functional zoning of forests in particular) contributed to forest conservation. Nowadays, protected forests, in which no silvicultural prescriptions or only special conservation prescriptions are allowed, occupy 39% of the total area [48]. The remaining areas are only eligible for selection cuttings (99.6%) and clear-cut prescriptions in a very small amount (0.4%). As a result, around 60% of the forest is older than 60 years old, with almost 20% of it being more than 100 years old [48], including over 200 hectares of virgin and quasi-virgin forests [49,50]. Within the context of temperate European forests, an important fraction of the Penteleu forest district area is occupied by relatively old forest.
The historical management of the different stands in Penteleu forest was reconstructed for the modelling approach following the ecological conditions and actual stand age structure and composition. In general, given the historical management changes and the information contained in the management plans [48], there is a high level of confidence in the historical management of the vast majority of the stands. An estimative historical management of Penteleu forest can be found in Supplementary Material S1.
The climate variables taken from van Oldenborgh et al.’s [51] multi-model mean climate change CMIP5 scenarios (https://climexp.knmi.nl/start.cgi; accessed on 1 July 2020) and additional climate series [43,52] for the region 25–26 E 45–46 N were statistically downscaled to 50 m resolution following Zorita et al. [53] using WorldClim 1.4 and 2.0 [54]. The downscaled monthly climate variables are the Tmin (monthly average of minimum temperature), Tmax (monthly average of maximum temperatures), PET (average monthly evapotranspiration), Prec (average monthly precipitation) and BAL (average monthly water balance). Spectral analysis of the climate variables is provided in Supplementary Material S2. Monthly radiation was also taken from WorldClim [54]. These climate products with past observations and future climate change forecasts until 2100 were used to analyse climate signal in dendrochronological series and in Landis-II simulations.

2.2. Dendrochronological Series and Climate Signal

In order to understand the principles of forests’ structure and functioning, a deep analysis of the past and present structure and functionality of virgin forests was performed in a 1 ha control plot inside the Penteleu forest. This control plot (circular, with a 56.41 m radius) was installed and measured in 2014 in a representative temperate montane conifer–broadleaved mixed–virgin forest stand. The stand is unevenly aged, and the reference plot includes 439 trees from three species (beech—63%, Norway spruce—26% and silver fir—11%) [48]. Being unaffected by human activity, we assumed that this stand is optimal for detecting climate signals on tree growth and species interactions.
Tree spatial distribution, breast height diameter (DBH) and height (h) were measured for all living trees over 8 cm DBH. The aboveground tree volume (v) of each tree was calculated with the following regression equation: log v = b0 + b1∙log d + b2∙log2 d +b3∙log h + b4∙log2h, where b0, b1, b2, b3 and b4 are the species’ national regression coefficients [55]. The overall stand volume measured was 803 m3. Tree aboveground biomass was computed using species wood density [55].
To build the reference dendrochronological series of silver fir, beech and Norway spruce, 406 radial core samples of living trees collected from the control plot were analysed and processed. The radial growth was calculated with CooRecorder 7.4 image analysis techniques [47,56]. Measurement checking and cross-dating was performed with COFECHA [42,57], and the growth series standardization was carried out with ASTRANwin [58]. The age structure was intensively analysed in an exploratory analysis to ascertain how stand development might have affected tree growth. DBH, measured through tree ring radial growth, was analysed by linear mixed models and analysis of variance in R software [59].
The statistical analysis of the dendrochronological series was performed with dplR [60], treeclim [61] and waveslim [62] packages in R software [59], and wavelet coherence among species dendrochronological series was explored with biwavelet [63]. The overall and moving correlations between climate variables (i.e., Tmin, Tmax, PET, Prec, BAL) and tree ring growth for the main species were calculated by treeclim dcc response and correlation function analysis [61]. In order to find interactions among species in mixed stands, the moving correlation matrices were inspected. This was completed with a factorial analysis of the of the species’ moving.
In addition, since there was a net change in the slope, the 1 ha plot was divided into 2 sectors during the fieldwork in order to explore and compare the relative species abundance in both sectors as a function of the local conditions. Being representative of 2 ecoregions, the reconstruction of the tree aboveground biomass in the 2 sectors was later used in the calibration of the Landis-II models.

2.3. Climate Change Projections

Landis-II [34,64], coupled with PnET models [35,65,66], was developed to simulate forest succession and the effect of climate [66] and management [67] under different climate change scenarios [66]. We built Landis-II models for the Penteleu forest district, aiming to assess management effects and the climate change impacts on Carpathian forest ecoregions, i.e., forests with homogeneous local conditions.
The model parameterization was performed following Gustafson [53] recommendations, using data extracted from scientific and technical works [12,68,69,70,71]. Since the full dataset had missing data issues, inconsistency, etc., the values were not used directly in the Landis-II model, but as input data for Landis-II calibration using the genetic algorithms (GA) [39,72] machine learning (ML) technique, with GA R package [38,39,59]. For calibration purposes, the GA fitness function for a given set of chromosome values used the difference in the biomass estimated for the 1 ha plot sectors (using the dendrochronological series) and the biomass predicted by the Landis-II models for those sectors (introduced in the model as fully independent ecoregions). In this, the supervised species calibration was bypassed, and species parameter values that are valid for the whole study area were set. Only the last period of the dendrochronological series was taken into account in order to minimize the effects of the stand dynamics effects on the calibration process.
After the calibration, the 50 m × 50 m resolution landscape was segmented in ecoregions combining the site type, forest types and relative productivity class in the management plans [48]. Soil properties for the ecoregions’ parameterization were also extracted from the management plan [48]. Overall, 22 ecoregions, covering an altitudinal gradient and a wide range of site properties, were established. Every ecoregion was provided for monthly photosynthetically active radiation (PAR) from WorldClim 2.0 radiation [54], downscaled CMIP5 Tmax, Tmin and Prec data (https://climexp.knmi.nl/start.cgi; accessed on 1 July 2020), and atmospheric CO2 concentrations for the 1500–2100 period [3] under the 2.6, 4.5, 6.0 and 8.5 greenhouse gas (GHG) atmospheric concentration RCP scenarios [3].
The 1685 forest stands in the Penteleu forest district were aggregated into 568 management system types, from which 1369 treatment prescriptions (silvicultural regeneration interventions) were sequentially combined and implemented over time at the stand and management system levels, according to historical changes in forest management, species composition and abundance, silvicultural system and age structure. All of the spatial information was managed in R software [59] spatial libraries [73,74,75,76].
The silvicultural systems implemented in Landis-II are as follows: (i) strict protection (mostly mixed uneven-aged forests; all virgin forests, quasi-virgin forests and some pure stands are covered by this management type); (ii) old selection forests (predominantly mixed uneven-aged forests); (iii) selection forests (natural composition, often mixed stands, with an uneven-aged structure); (iv) quasi-selection systems (transformation to selection forests 60–70 years; natural composition with relative even to relative uneven-aged structures); (v) group selection system (even aged with pure and mixed composition); (vi) clear-cut system (mostly even-age pure stands and some mixed stands with non-native species).
The timeframe covered by the Landis-II models comprises four main management types and periods: (i) virgin forests and traditional selection forests (1900–1940); (ii) salvage logging (1941–1973); (iii) selection forests (1973–2020); (iv) increased protection of forests (2021–2100).
Disturbances (e.g., windthrows) were modelled through harvest extension, because (a) management plans show evidence of their low relevance and (b) the salvage logging, which is commonly practiced in disturbed stands, has similar effects to unplanned clearcuttings.
The Landis-II models, calibrated and validated for the entire research area, were computed for past conditions in the 1900–2010 period and for the 2010–2100 period under the RCP 2.6, 4.5, 6.0 and 8.5 climate change scenarios. Biomass and, implicitly, C stock, capitalizing on the energy and mass fluxes within the modelling approach, were the main target variables in this study. The ecoregions’ biomass over time and the biomass in the natural distribution areas of pure and mixed forests were fitted and analysed using linear mixed models and an analysis of variance. Three key periods (1901–1910, 2001–2010 and 2091–2100) were targeted and compared. Comparisons between those periods were performed for every ecoregion and climate change scenario using Benjamini–Hochberg FDR correction [77]. The familywise and post hoc comparisons were carried out with nlme [78] and multcomp packages [79] in R software [59].

3. Results

3.1. Forest Sensitivity to Past Climate

Older trees in the 1 ha plot mixed–virgin forest stand are more than 350 years old and the maximum tree DBH is 110 cm. Exploratory analysis of the age structure and recruitment suggest that, despite the truncated diameter-based sampling, the age structure follows a geometric distribution, with a strong predominance of young age classes. This age structure has been quite stable over time, emphasized by the small peaks of recruitment detected over the last 300 years.
The analysis of the climate signal in a mixed–virgin forest stand indicates that tree radial growth depends on climate and interspecific relationships. Tree DBH, reconstructed from radial core samples collected in the Penteleu 1 ha plot, increased approximately linearly over time (Figure 2) at the stand level. The mean growth of the central 90% (0.05–0.95 quantiles) of trees in the plot ranged from 0.9 to 4.2 cm∙year−1. At the tree level, radial growth has the general tendency to be constant over time, with a few exceptions. Radial growth changed strongly over time for a reduced number of trees: some trees with intense radial growth in the early stages of development gradually reduced their growth rate when they reached physiological maturity, while other trees with reduced radial growth in the early stages suddenly increased their growth rate in later stages. The dendrochronological dataset analysis showed a non-significant effect of the tree age on the tree ring width (p-value > 0.05).
All three biwavelets among pairs of species (Figure 3) contain regions with significant wavelet coherence among dendrochronological series of pairs of species. Anti-phase relationships dominate in the beech–silver fir biwavelet, while beech–Norway spruce and silver fir–Norway spruce biwavelets tend to be in phase when there is a lag in the phase.
Climate has a significant effect on tree ring growth. At the stand level, it was summarised the relationship between some of the main climate variables (i.e., Tmin, Tmax, PET, Prec and BAL) and the radial growth extracted from tree cores through moving correlations (details can be found in Supplementary Material S3), and then assessed the similarities between the responses of the three main species in the study area (silver fir, beech and Norway spruce) to the climate variables (Figure 4). The results showed synergic and antagonistic responses to the climate. Silver fir showed a similar trend to the Norway spruce, with the first not being as significant as the second. Silver fir tended to show negative correlations with climate variables in periods in which beech and Norway spruce growth were favoured. Beginning with the end of the 19th century and until the beginning of the 20th century, a consistently positive growth rate was recorded for beech, correlated with higher-than-average summer temperatures. Beech’s correlation with September Prec and BAL varied over time, from a negative correlation to a positive one. Beech growth also had a positive correlation with late spring Prec. Norway spruce, in general, tended to show negative correlations with climate variables in the same periods in which beech was favoured by the climate. Beech moving correlations around 1750 showed significant negative correlations with September Prec and BAL, as in 1840. Additionally, around 1840, positive growth was associated with spring Prec and BAL. Later, around 1900, positive correlations were detected with spring Tmin, Tmax and PET. Around 1965, growth was positively correlated with September Prec and BAL and negatively correlated with July and September Tmin. Silver fir growth correlated negatively with spring Tmin and PET and positively with late-summer Prec and BAL around 1970. During the 19th century, silver fir growth was negatively correlated with May Tmin, Tmax and PET, although it was positively correlated with September and April Tmin and Tmax and also with May Prec and BAL. At the beginning of the 20th century, silver fir growth was negatively correlated with late-summer Tmin, Tmax and Prec and positively correlated with April Prec and BAL. Later, occasional positive significant positive correlations were detected with Tmin, Tmax and PET, while negative correlations were detected with July and August Prec and August BAL. Finally, no correlations were found for the 21st century. Norway spruce growth showed positive significant correlations with September Tmin, Tmax and PET around 1825, though negative correlations were found with April Prec and BAL. Later, until 1990, positive correlations tended to occur with spring and summer Tmin, Tmax and PET, particularly around 1875, 1925 and 1970. Significant positive correlations with April and June Prec were detected in the last decades of the 19th century and at the beginning of the 20th century. Finally, around 1945 and 1985, Norway spruce growth was eventually negatively correlated with late-summer Prec and BAL.
Thus, the first and second components of the factorial analysis of the moving correlation matrices were built based on the antagonistic behaviour of beech and silver fir. Norway spruce displayed an intermediate behaviour, and consequently, less weight was placed on the construction of the first and second components.
Tree biomass reconstructed from tree cores showed that, at the stand level, around one third of the current biomass in the Penteleu 1 ha plot was produced more than 150 years ago and approximately half was formed in the last 50 years (Figure 5a). Regarding species biomass, before 1850, most of it belonged to beech. After 1900, only half of the biomass found was attributed to beech, with the remaining biomass being provided by conifers, particularly silver fir.
Due to the plot’s heterogeneity and microtopography, we compared the biomass accumulation on steep and flat terrain. Before 1850, the accumulated tree biomass was higher in the flattest sector of the plot compared to the steepest one. After 1950, the biomass per hectare in the steepest sector (Figure 5b) was equal to the biomass in the flat sector. As extracted from the dendrochronological series, beech and silver fir are the dominant and codominant species in both sectors, with silver fir being dominant in the steep sector and beech in the flat. Norway spruce, having the lowest contribution in terms of biomass, accumulated more than 1.5 times the biomass per hectare in the flat sector, as compared to the amount of biomass measured in the steep sector.

3.2. Aboveground Living Biomass under Climate Change Scenarios

Landis-II simulations show that the forest biomass changed from 1900 to 2010 in almost all of the ecoregions (Figure 6). Some ecoregions showed a moderate increase in biomass throughout the period. Other ecoregions hardly changed, while the remaining few displayed a reduced amount of biomass. Biomass reductions in the ecoregions were often associated with the harvesting process.
Harvesting caused very strong biomass reductions in some ecoregions, particularly through clearcutting during the 1941–1973 period. However, biomass reductions were also found in other periods, harvest types and management systems, such as selection cuttings in shelterwood systems. The rise in the level of forest protection, together with the aging of the forest, contributed to a reduction in biomass loss and to moderate biomass accumulation in recent decades. In ecoregions with null or very low harvest intensity, the amount of biomass tends to increase.
Management, which, in the past, caused strong reductions in biomass through intense harvest, will contribute to biomass accumulation in the future through protection and low harvest pressure. The Landis-II projections show that, from 2010 to 2100, under a more protective management, the amount of biomass in most of the ecoregions will increase, particularly after 2025.
Apart from harvest effects, the different climate change scenarios (RCP 2.6 to RCP 8.5) had a significant role on forest biomass accumulation after 2010 (p-value ≤ 0.05). This increase in biomass was particularly noticeable in periods where harvesting was less intense. The intensity of the growth differed strongly among the ecoregions. After 2010, the accumulation of biomass is a result of a combination of both local ecoregion conditions, management and climate change scenarios. There is no common response to climate scenarios for all ecoregions, even though intense climate change scenarios, particularly RCP 8.5, tend to increase biomass accumulation. Indeed, some ecoregions have the lowest accumulated biomass under intermediate climate change scenarios (RCP 4.5 and RCP 6.0) and the highest under RCP 8.5. The amount of accumulated biomass barely changed over time in some ecoregions.
The ecoregions, aggregated according to their domination by pure and mixed forests (Figure 7), had different biomass accumulated at the beginning of the study period and over time. In 1900, when salvage logging had not begun, the biomass per hectare in pure stands was higher than in mixed stands. Around 2010, the amount of biomass was still higher in pure forest ecoregions. However, in the last period of the simulations (up to 2100), the biomass per hectare in mixed forests was similar to or even higher than that of pure forests, particularly in RCP 8.5, where the biomass per hectare in mixed forests broadly surpasses that of pure forests.

4. Discussion

This study shows that, even though virgin forests have a proven stability, both tree radial growth and biomass are affected by climate change and species interactions, thus providing long-term faithful records that are unaffected by management. In addition, it explores the reconstruction of the forest stand aboveground biomass in the past, projection in the future under different climate change scenarios and the C sequestration, as part of the atmospheric CO2 balance.
Tree DBH reconstructed from Penteleu dendrochronological series was approximately linear, as also found in other similar studies [8,31,80]. Thus, radial growth is not particularly affected by tree age. Radial growth and, thus, the DBH–age relationship, is conditioned by the environmentally limiting conditions [80] and by biological stress. Additionally, the change in the growth rates of individual trees often depends on their interactions with neighbouring trees (e.g., facilitation and competition) and other processes (e.g., mortality) [11,81], meaning that they are likely the cause of the changes found in the measured trees in Penteleu.
We found that small differences in the ecological preferences, local conditions and climate lead to several spatiotemporal patterns in the biomass accumulation over time. In mixed–virgin forest stands, intra- and inter-specific interactions have a strong impact on tree growth, while also having further consequences at the stand level. Such interactions differ in relation to local conditions and species dominance and abundance, which can lead to modifications even at short distances. The results showed that there are differences between growth rates of some trees in early vs. late stages, which (presumably) strongly depend on the physiological activity and dynamics of the in-stand competition and mortality [11,81] of the species. In the Penteleu virgin forest 1 ha plot, the steepest areas are dominated by silver fir, which is less susceptible to warmer and drier conditions than beech [13], while flat areas are occupied by beech, which is also accompanied with a higher abundance of Norway spruce. Having a higher abundance in the flattest area, beech, compared to silver fir and Norway spruce shows better resistance to different stress factors and also has a significant growth in the understory [23,81]. Due to withstand such conditions (including drought), beech growth usually has a low variability [11] compared to coniferous species. These heterogeneous responses to the climate are produced in mixed forests even when the climate conditions are favourable for all three species.
The dendrochronological series analysis proved that silver fir and beech tend to have opposite behaviours in response to climate in mixed stands, while Norway spruce, displaying intermediate behaviour, tends to negatively correlate with the climate when beech correlations are positive. Having different ecological preferences and life traits, the beech–silver fir association in mixed stands is advantageous in terms of their ability to endure drought [13,17,31,33,82] and also because of their efficient light use [83,84]. There is also evidence that beech–Norway spruce mixed stands perform better than beech pure stands [13,33] due to the complementarity of the two species.
The overall biomass, capitalizing all ecological processes, is homogeneously distributed over the 1 ha plot, highlighting the stability and species complementarity in mixed forests, as other studies have suggested [13,18,19,33]. However, despite their high efficiency in terms of the use of resources, there are concerns on the response of mixed and pure forests to climate change (e.g., Piovesan et al. [32]). It was emphasized that pure beech forests have already been affected by climate-change-induced drought in recent decades [32], as Norway spruce shows high drought resistance when compared to beech [17,82,85]. Additionally, silver fir is expected to be less susceptible to warmer and drier conditions than beech [13].
Since tree ring measurements include the effects of past conditions, such types of long-term data can be used in a modelling approach to analyse past changes and to foresee, anticipate and manage future changes [86,87]. However, the complexity, variability and information gaps contained in biological records [87,88] often require solutions such as ML methods [37,86]. GA [38,39] provided a suitable mechanism to speed up and objectivize the modelling approach and was successfully used in this study for Landis-II model calibration in order to simulate landscape changes and for forest change forecasting under different climate change scenarios, which proved to be robust.
The Landis-II models demonstrate that management had a strong impact on forest structure and functioning over time, in agreement with official records and statistics [48]. Forest resilience to management in terms of relative biomass compared to the pre-harvest level is high at the ecoregion level. Biomass commonly recovers in less than 50 years, although impacts on other forest attributes, e.g., diversity, last much longer [48]. Forest resistance to management is very low, as denoted by both (a) the impact of punctual harvest events at the ecoregion level and (b) sustained biomass stock changes due to the change in the silvicultural system, following stakeholders and policymakers’ decisions in the past and expectations after 2025, when even more protective management will be implemented.
Biomass change in the different ecological regions is not only related to management, but also to climate change effects. The most important impacts of climate change in the Carpathians will be caused by rising temperature and CO2 levels. Precipitation will be similar to the current conditions in some areas and may even slightly increase [4]; therefore, drought events are not expected to have a significant role, other than in a few ecological regions and marginal areas of the Carpathians [4]. Temperate forests will benefit from relative stability and better conditions for growth [2,4], contrary to Mediterranean forests, where droughts and fires will be the main drivers of change [89].
Some areas in the Carpathians will endure harsher conditions [4], and our results indicate that impacts will also differ at the local scale, with some ecoregions being more affected than others.
Uncertainty of climate change impacts and characteristics on forest ecosystems is relatively high [30,90], and our results show that temperate forest resistance to climate change is relatively low; however, forest biomass is predicted to be relatively similar under different climate change scenarios. Despite the findings of Nabuurs [2], projections from our Landis-II model also show that the Penteleu forest’s C sequestration capacity is still far from saturation, and the aboveground biomass will continue to grow at least until the end of the current century, even when signs of deceleration appear at the end of the period in some ecoregions.
It is important to note that this C pool-growing trend is also related to a rise in the level of forest protection provided, forest aging and the decreasing harvest pressure throughout the period. Thus, the management of temperate forests has the potential to increase carbon sinks and mitigate climate change [1].
On the other hand, adaptation strategies can focus on sustainable management. The general interest is to promote C sequestration and management should be oriented toward forest adaptation to climate change, by reducing potential negative climate change impacts, e.g., drought events [2,89], etc.
In this regard, management can be used to promote mixed forests, which reach higher biomass accumulation under intense climate change scenarios, and thus, to maximize C sequestration. Mixed stands and their particularities contribute to forest resistance and resilience at the landscape level, in comparison to pure stands [15,16,17].

5. Conclusions

The dendrochronological series in mixed stands not only contain climate signal, but species interactions, and they are useful for past events and projections. Beech, Norway spruce and silver fir, whose radial growth was not dependent of tree age, have complementary responses in mixed stands.
PnET-Landis-II models, GA calibrated using dendrochronological records, applied to the Penteleu forest, which hosts virgin and managed stands, brings to light (i) the relevance of virgin forests in vegetation long-term monitoring and as a reference data source for all forests in the Carpathians, (ii) the similarities and differences in the structure and functioning of mixed and pure forest stands and their adaptative management and (iii) the effect of climate on tree growth and climate change impact on forests, in which mitigation and adaptation strategies, programs and actions must be implemented.
Landis-II models provide valuable insights into the medium- to long-term succession for forest management and forest planning. They provide a good approach to climate change impacts, which are becoming increasingly accentuated. Landis-II simulations revealed that aboveground forest biomass changed in almost all ecological regions in the Carpathian Curvature, often associated with harvesting. Forest projections demonstrate that climate change in the Carpathian Curvature will promote forest aboveground biomass and C sequestration. Climate change will modify species and cohort interactions, leading to changes in the forest structure and functioning. Climate change impacts will depend on climate change intensity and local conditions and management. Among forest ecosystems, mixed forests are capable of higher C sequestration and higher biomass accumulation under intense climate change.
In this study, management, which is one of the main factors modifying forest structure and succession, can be used to promote resilient mixed forests, which are expected to accumulate higher biomass quantity under intense climate change, and thus, contributes to climate change mitigation and adaptation.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/f12070885/s1, Figure S1: Estimative historical management, Figure S2: Spectral analysis of climate variables (i.e., Tmin, Tmax, PET, Prec, BAL), Figure S3: Moving correlations between beech (F. sylvatica L.), Norway spruce (P. abies (L.) H. Karst.) and silver fir (A. alba Mill.) tree ring growth and climate variables (i.e., Tmin, Tmax, PET, Prec, BAL).

Author Contributions

Conceptualization, J.G.-D. and S.C.; methodology, J.G.-D., S.C. and Ș.L.; software, J.G.-D. and S.C.; validation, J.G.-D., S.C., D.P. and O.B.; formal analysis, J.G.-D.; investigation, J.G.-D. and S.C.; resources, S.C. and J.G.-D.; data curation, J.G.-D.; writing—original draft preparation, J.G.-D. and S.C.; writing—review and editing, J.G.-D., S.C., D.P., Ș.L. and O.B.; visualization, J.G.-D.; supervision, S.C. and O.B.; project administration, S.C. and O.B.; funding acquisition, S.C. and O.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the BIOSERV Programme, Project IDs PN19070103, PN19070102 and EO-ROFORMON project, ID P_37_651/SMIS 105058.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dymond, C.C.; Beukema, S.; Nitschke, C.R.; Coates, K.D.; Scheller, R.M. Carbon Sequestration in Managed Temperate Coniferous Forests under Climate Change. Biogeosciences 2016, 13, 1933–1947. [Google Scholar] [CrossRef] [Green Version]
  2. Nabuurs, G.-J.; Lindner, M.; Verkerk, P.J.; Gunia, K.; Deda, P.; Michalak, R.; Grassi, G. First Signs of Carbon Sink Saturation in European Forest Biomass. Nat. Clim. Chang. 2013, 3, 792–796. [Google Scholar] [CrossRef]
  3. van Vuuren, D.P.; Edmonds, J.; Kainuma, M.; Riahi, K.; Thomson, A.; Hibbard, K.; Hurtt, G.C.; Kram, T.; Krey, V.; Lamarque, J.-F.; et al. The Representative Concentration Pathways: An Overview. Clim. Chang. 2011, 109, 5–31. [Google Scholar] [CrossRef]
  4. Micu, D.M.; Dumitrescu, A.; Cheval, S.; Birsan, M.-V. Climate of the Romanian Carpathians; Springer Atmospheric Sciences; Springer International Publishing: Cham, Switzerlands, 2015; ISBN 978-3-319-02885-9. [Google Scholar]
  5. Huntley, B.; Baxter, R. Vegetation Ecology and Global Change. In Vegetation Ecology; van der Maarel, E., Franklin, J., Eds.; John Wiley & Sons, Ltd.: Oxford, UK, 2013; pp. 509–530. ISBN 978-1-118-45259-2. [Google Scholar]
  6. McDowell, N.G.; Allen, C.D.; Anderson-Teixeira, K.; Aukema, B.H.; Bond-Lamberty, B.; Chini, L.; Clark, J.S.; Dietze, M.; Grossiord, C.; Hanbury-Brown, A.; et al. Pervasive Shifts in Forest Dynamics in a Changing World. Science 2020, 368, eaaz9463. [Google Scholar] [CrossRef]
  7. 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. [Google Scholar] [CrossRef] [Green Version]
  8. Pretzsch, H. The Course of Tree Growth. Theory and Reality. For. Ecol. Manag. 2020, 478, 118508. [Google Scholar] [CrossRef]
  9. Chivulescu, S.; Leca, S.; Silaghi, D.; Cristea, V. Structural Biodiversity and Dead Wood in Virgin Forests from Eastern Carpathians. J. Agric. For. 2018, 64. [Google Scholar] [CrossRef] [Green Version]
  10. Grantham, H.S.; Duncan, A.; Evans, T.D.; Jones, K.R.; Beyer, H.L.; Schuster, R.; Walston, J.; Ray, J.C.; Robinson, J.G.; Callow, M.; et al. Anthropogenic Modification of Forests Means Only 40% of Remaining Forests Have High Ecosystem Integrity. Nat. Commun. 2020, 11, 5978. [Google Scholar] [CrossRef]
  11. Vrška, T.; Adam, D.; Hort, L.; Kolář, T.; Janík, D. European Beech (Fagus Sylvatica L.) and Silver Fir (Abies Alba Mill.) Rotation in the Carpathians—A Developmental Cycle or a Linear Trend Induced by Man? For. Ecol. Manag. 2009, 258, 347–356. [Google Scholar] [CrossRef]
  12. Şofletea, N.; Curtu, L. Dendrologie; Pentru Viaţă: Braşov, Romania, 2001; Volume II, ISBN 978-973-99456-1-5. [Google Scholar]
  13. Bosela, M.; Lukac, M.; Castagneri, D.; Sedmák, R.; Biber, P.; Carrer, M.; Konôpka, B.; Nola, P.; Nagel, T.A.; Popa, I.; et al. Contrasting Effects of Environmental Change on the Radial Growth of Co-Occurring Beech and Fir Trees across Europe. Sci. Total Environ. 2018, 615, 1460–1469. [Google Scholar] [CrossRef]
  14. Primicia, I.; Camarero, J.J.; Janda, P.; Čada, V.; Morrissey, R.C.; Trotsiuk, V.; Bače, R.; Teodosiu, M.; Svoboda, M. Age, Competition, Disturbance and Elevation Effects on Tree and Stand Growth Response of Primary Picea Abies Forest to Climate. For. Ecol. Manag. 2015, 354, 77–86. [Google Scholar] [CrossRef]
  15. Griess, V.C.; Acevedo, R.; Härtl, F.; Staupendahl, K.; Knoke, T. Does Mixing Tree Species Enhance Stand Resistance against Natural Hazards? A Case Study for Spruce. For. Ecol. Manag. 2012, 267, 284–296. [Google Scholar] [CrossRef]
  16. Pretzsch, H.; Schütze, G.; Uhl, E. Resistance of European Tree Species to Drought Stress in Mixed versus Pure Forests: Evidence of Stress Release by Inter-Specific Facilitation: Drought Stress Release by Inter-Specific Facilitation. Plant Biol. 2013, 15, 483–495. [Google Scholar] [CrossRef] [PubMed]
  17. Magh, R.-K.; Bonn, B.; Grote, R.; Burzlaff, T.; Pfautsch, S.; Rennenberg, H. Drought Superimposes the Positive Effect of Silver Fir on Water Relations of European Beech in Mature Forest Stands. Forests 2019, 10, 897. [Google Scholar] [CrossRef] [Green Version]
  18. MacPherson, D.M.; Lieffers, V.J.; Blenis, P.V. Productivity of Aspen Stands with and without a Spruce Understory in Alberta’s Boreal Mixedwood Forests. For. Chron. 2001, 77, 351–356. [Google Scholar] [CrossRef] [Green Version]
  19. Man, R.; Lieffers, V.J. Are Mixtures of Aspen and White Spruce More Productive than Single Species Stands? For. Chron. 1999, 75, 505–513. [Google Scholar] [CrossRef]
  20. Laiho, O.; Lahde, E.; Pukkala, T. Uneven- vs Even-Aged Management in Finnish Boreal Forests. Forestry 2011, 84, 547–556. [Google Scholar] [CrossRef]
  21. Barsoum, N.; Coote, L.; Eycott, A.E.; Fuller, L.; Kiewitt, A.; Davies, R.G. Diversity, Functional Structure and Functional Redundancy of Woodland Plant Communities: How Do Mixed Tree Species Plantations Compare with Monocultures? For. Ecol. Manag. 2016, 382, 244–256. [Google Scholar] [CrossRef] [Green Version]
  22. Petritan, I.C.; Commarmot, B.; Hobi, M.L.; Petritan, A.M.; Bigler, C.; Abrudan, I.V.; Rigling, A. Structural Patterns of Beech and Silver Fir Suggest Stability and Resilience of the Virgin Forest Sinca in the Southern Carpathians, Romania. For. Ecol. Manag. 2015, 356, 184–195. [Google Scholar] [CrossRef]
  23. Chivulescu, S.; Ciceu, A.; Leca, S.; Apostol, B.; Popescu, O.; Badea, O. Development Phases and Structural Characteristics of the Penteleu-Viforta Virgin Forest in the Curvature Carpathians. IForest Biogeosci. For. 2020, 13, 389–395. [Google Scholar] [CrossRef]
  24. Mirtl, M. Introducing the Next Generation of Ecosystem Research in Europe: LTER-Europe’s Multi-Functional and Multi-Scale Approach. In Long-Term Ecological Research; Müller, F., Baessler, C., Schubert, H., Klotz, S., Eds.; Springer: Dordrecht, The Netherlands, 2010; pp. 75–93. ISBN 978-90-481-8781-2. [Google Scholar]
  25. Parviainen, J. Virgin and Natural Forests in the Temperate Zone of Europe. For. Snow Landsc. Res. 2005, 79, 9–18. [Google Scholar]
  26. Dănescu, A.; Albrecht, A.T.; Bauhus, J. Structural Diversity Promotes Productivity of Mixed, Uneven-Aged Forests in Southwestern Germany. Oecologia 2016, 182, 319–333. [Google Scholar] [CrossRef] [PubMed]
  27. Vacek, S.; Vacek, Z.; Podrázský, V.; Bílek, L.; Bulušek, D.; Štefančík, I.; Remeš, J.; Štícha, V.; Ambrož, R. Structural Diversity of Autochthonous Beech Forests in Broumovské Stěny National Nature Reserve, Czech Republic. Austrian J. For. Sci. 2014, 131, 191–214. [Google Scholar]
  28. Duveneck, M.J.; Scheller, R.M. Measuring and Managing Resistance and Resilience under Climate Change in Northern Great Lake Forests (USA). Landsc. Ecol. 2016, 31, 669–686. [Google Scholar] [CrossRef]
  29. Curovic, M.; Spalevic, V.; Sestras, P.; Motta, R.; Dan, C.; Garbarino, M.; Vitali, M.; Urbinati, C. Structural and Ecological Characteristics of Mixed Broadleaved Old-Growth Forest (Biogradska Gora-Montenegro). Turk. J. Agric. For. 2020, 44, 428–438. [Google Scholar] [CrossRef]
  30. Bouriaud, L.; Bouriaud, O.; Elkin, C.; Temperli, C.; Reyer, C.; Duduman, G.; Barnoaiea, I.; Nichiforel, L.; Zimmermann, N.; Bugmann, H. Age-Class Disequilibrium as an Opportunity for Adaptive Forest Management in the Carpathian Mountains, Romania. Reg. Environ. Chang. 2015, 15, 1557–1568. [Google Scholar] [CrossRef]
  31. Pretzsch, H.; Hilmers, T.; Uhl, E.; Bielak, K.; Bosela, M.; del Rio, M.; Dobor, L.; Forrester, D.I.; Nagel, T.A.; Pach, M.; et al. European Beech Stem Diameter Grows Better in Mixed than in Mono-Specific Stands at the Edge of Its Distribution in Mountain Forests. Eur. J. For. Res. 2020. [Google Scholar] [CrossRef]
  32. Piovesan, G.; Biondi, F.; Filippo, A.D.; 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]
  33. Schwarz, J.A.; Bauhus, J. Benefits of Mixtures on Growth Performance of Silver Fir (Abies Alba) and European Beech (Fagus Sylvatica) Increase with Tree Size Without Reducing Drought Tolerance. Front. For. Glob. Chang. 2019, 2, 79. [Google Scholar] [CrossRef]
  34. Scheller, R.M.; Domingo, J.B.; Sturtevant, B.R.; Williams, J.S.; Rudy, A.; Gustafson, E.J.; Mladenoff, D.J. Design, Development, and Application of LANDIS-II, a Spatial Landscape Simulation Model with Flexible Temporal and Spatial Resolution. Ecol. Model. 2007, 201, 409–419. [Google Scholar] [CrossRef]
  35. Aber, J.D.; Ollinger, S.V.; Federer, C.A.; Reich, P.B.; Goulden, M.L.; Sicklighter, D.W.; Melillo, J.M.; Lathrop, R.G. Predicting the Effects of Climate Change on Water Yield and Forest Production in the Northeastern United States. Clim. Res. 1995, 5, 207–222. [Google Scholar] [CrossRef]
  36. Gustafson, E.J.; Sturtevant, B.R.; de Bruijn, A.M.G.; Lichti, N.; Jacobs, D.F.; Kashian, D.M.; Miranda, B.R.; Townsend, P.A. Forecasting Effects of Tree Species Reintroduction Strategies on Carbon Stocks in a Future without Historical Analog. Glob. Chang. Biol. 2018, 24, 5500–5517. [Google Scholar] [CrossRef] [PubMed]
  37. Bishop, C.M. Pattern Recognition and Machine Learning; Information Science and Statistics; Springer: New York, NY, USA, 2006; ISBN 978-0-387-31073-2. [Google Scholar]
  38. Scrucca, L. On Some Extensions to GA Package: Hybrid Optimisation, Parallelisation and Islands Evolution. R J. 2017, 9, 187–206. [Google Scholar] [CrossRef] [Green Version]
  39. Scrucca, L. GA: A Package for Genetic Algorithms in R. J. Stat. Softw. 2013, 53, 1–37. [Google Scholar] [CrossRef] [Green Version]
  40. Fritts, H.C. Tree Rings and Climate; Academic Press: London, UK; New York, NY, USA, 1976; ISBN 978-0-12-268450-0. [Google Scholar]
  41. Briffa, K.R.; Jones, P.D.; Wigley, T.M.; Pilcher, J.R.; Baillie, M.G.L. Climate Reconstruction from Tree Rings: Part 1, Basic Methodology and Preliminary Results for England. J. Climatol. 1983, 3, 233–242. [Google Scholar] [CrossRef]
  42. Cook, E.R.; Johnson, A.H.; Blasing, T.J. Forest Decline: Modeling the Effect of Climate in Tree Rings12. Tree Physiol. 1987, 3, 27–40. [Google Scholar] [CrossRef]
  43. Pauling, A.; Luterbacher, J.; Casty, C.; Wanner, H. Five Hundred Years of Gridded High-Resolution Precipitation Reconstructions over Europe and the Connection to Large-Scale Circulation. Clim. Dyn. 2006, 26, 387–405. [Google Scholar] [CrossRef] [Green Version]
  44. Popa, I. Fundamente Metodologice Şi Aplicaţii de Dendrocronologie; Editura Tehnicǎ Silvicǎ, 2004. [Google Scholar]
  45. Hasenauer, H. Dimensional Relationships of Open-Grown Trees in Austria. For. Ecol. Manag. 1997, 96, 197–206. [Google Scholar] [CrossRef]
  46. Kuemmerle, T.; Chaskovskyy, O.; Knorn, J.; Radeloff, V.C.; Kruhlov, I.; Keeton, W.S.; Hostert, P. Forest Cover Change and Illegal Logging in the Ukrainian Carpathians in the Transition Period from 1988 to 2007. Remote Sens. Environ. 2009, 113, 1194–1207. [Google Scholar] [CrossRef]
  47. Sidor, C.G.; Popa, I.; Vlad, R.; Roibu, C.C. Rețeaua Națională de Serii Dendrocronologice Pentru Pinul Silvestru (Pinus Sylvestris) Din România-PIDECRO; Editura Silvica: Voluntari, Romania, 2020; ISBN 978-606-8020-73-0. [Google Scholar]
  48. ICAS. Amenajamentul. O.S. Penteleu; Regia Naţională a pădurilor Romsilva; Institutul de Cercetări şi amenajări Silvice: Bucureşti, Romania, 2012. [Google Scholar]
  49. INCDS. Amenajamentul. O.S. Penteleu; Regia Naţională a pădurilor Romsilva; Institutul de Cercetări şi amenajări Silvice: Bucureşti, Romania, 2021. [Google Scholar]
  50. Cristea, V.C. Dendrometrics and Auxological Identification and Characterization of Forests with Natural and Semi-Natural Structure from Penteleu Massif, Transilvania. Ph.D. Thesis, University of Brașov, Brașov, Romania, 2019. [Google Scholar]
  51. Van Oldenborgh, G.J.; Drijfhout, S.; van Ulden, A.; Haarsma, R.; Sterl, A.; Severijns, C.; Hazeleger, W.; Dijkstra, H. Western Europe Is Warming Much Faster than Expected. Clim. Past 2009, 5, 1–12. [Google Scholar] [CrossRef] [Green Version]
  52. Luterbacher, J. European Seasonal and Annual Temperature Variability, Trends, and Extremes since 1500. Science 2004, 303, 1499–1503. [Google Scholar] [CrossRef] [Green Version]
  53. Zorita, E.; von Storch, H. The Analog Method as a Simple Statistical Downscaling Technique: Comparison with More Complicated Methods. J. Clim. 1999, 12, 2474–2489. [Google Scholar] [CrossRef]
  54. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-Km Spatial Resolution Climate Surfaces for Global Land Areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  55. Giurgiu, V.; Decei, I.; Drăghiciu, D. Metode Şi Tabele Dendrometrice [Methods and Dendrometric Tables]; Ceres Publishing House, Bucharest: Bucharest, Romania, 2004; ISBN 973-40-0639-8. [Google Scholar]
  56. Larson, L.A. Online User Manual for Cybis CooRecorder and CDendro Programs. Available online: http://www.cybis.se/forfun/dendro/ (accessed on 1 July 2020).
  57. Holmes, R.L. Computer-Assisted Quality Control in Tree-Ring Dating and Measurement. Tree-Ring Bull. 1983, 3, 69–78. [Google Scholar]
  58. Cook, B.D.; Krusic, P.J. ARSTAN v. 4.1d: A Tree-Ring Standardization Program Based on Detrending and Autoregressive Time Series Modeling, with Interactive Graphics; Tree-Ring Laboratory, Lamont-Doherty Earth Observatory of Columbia University: Palisades, NY, USA, 2005. [Google Scholar]
  59. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  60. Bunn, A.G. A Dendrochronology Program Library in R (DplR). Dendrochronologia 2008, 26, 115–124. [Google Scholar] [CrossRef]
  61. Zang, C.; Biondi, F. Treeclim: An R Package for the Numerical Calibration of Proxy-Climate Relationships. Ecography 2015, 38, 431–436. [Google Scholar] [CrossRef]
  62. Whitcher, B. Waveslim: Basic Wavelet Routines for One-, Two-, and Three-Dimensional Signal Processing; R Package Version 1.8.2. Available online: https://cran.r-project.org/web/packages/waveslim/ (accessed on 1 January 2021).
  63. Gouhier, T.C.; Grinsted, A.; Simko, V. R Package Biwavelet: Conduct Univariate and Bivariate Wavelet Analyses. Available online: https://cran.r-project.org/web/packages/biwavelet/ (accessed on 1 January 2021).
  64. Duveneck, M.J.; Scheller, R.M.; White, M.A. Effects of Alternative Forest Management on Biomass and Species Diversity in the Face of Climate Change in the Northern Great Lakes Region (USA). Can. J. For. Res. 2014, 44, 700–710. [Google Scholar] [CrossRef]
  65. Scheller, R.M.; Mladenoff, D.J. A Forest Growth and Biomass Module for a Landscape Simulation Model, LANDIS: Design, Validation, and Application. Model. Disturb. Succession For. Landsc. Using LANDIS 2004, 180, 211–229. [Google Scholar] [CrossRef]
  66. Gustafson, E.J.; De Bruijn, A.M.G.; Pangle, R.E.; Limousin, J.-M.; McDowell, N.G.; Pockman, W.T.; Sturtevant, B.R.; Muss, J.D.; Kubiske, M.E. Integrating Ecophysiology and Forest Landscape Models to Improve Projections of Drought Effects under Climate Change. Glob. Chang. Biol. 2015, 21, 843–856. [Google Scholar] [CrossRef]
  67. Gustafson, E.J.; Shifley, S.R.; Mladenoff, D.J.; Nimerfro, K.K.; He, H.S. Spatial Simulation of Forest Succession and Timber Harvesting Using LANDIS. Can. J. For. Res. 2000, 30, 32–43. [Google Scholar] [CrossRef]
  68. Kattge, J.; Díaz, S.; Lavorel, S.; Prentice, I.C.; Leadley, P.; Bönisch, G.; Garnier, E.; Westoby, M.; Reich, P.B.; Wright, I.J.; et al. TRY–A Global Database of Plant Traits. Glob. Chang. Biol. 2011, 17, 2905–2935. [Google Scholar] [CrossRef]
  69. Vittoz, P.; Engler, R. Seed Dispersal Distances: A Typology Based on Dispersal Modes and Plant Traits. Bot. Helv. 2007, 117, 109–124. [Google Scholar] [CrossRef] [Green Version]
  70. Giurgiu, V.; Draghiciu, D. Modele Matematico-Auxologice Si Tabele de Productie Pentru Arborete [Mathematic-Auxologic Models and Yield Tables for Forest Stands]; Ceres Publishing House: Bucharest, Romania, 2004; ISBN 973-40-0637-1. [Google Scholar]
  71. Tanase, M.A.; Villard, L.; Pitar, D.; Apostol, B.; Petrila, M.; Chivulescu, S.; Leca, S.; Borlaf-Mena, I.; Pascu, I.-S.; Dobre, A.-C.; et al. Synthetic Aperture Radar Sensitivity to Forest Changes: A Simulations-Based Study for the Romanian Forests. Sci. Total Environ. 2019, 689, 1104–1114. [Google Scholar] [CrossRef]
  72. Schmitt, L.M. Theory of Genetic Algorithms. Theor. Comput. Sci. 2001, 259, 1–61. [Google Scholar] [CrossRef] [Green Version]
  73. Bivand, R.; Rundel, C. Rgeos: Interface to Geometry Engine-Open Source (‘GEOS’). Available online: https://cran.r-project.org/web/packages/rgeos/ (accessed on 1 January 2021).
  74. Hijmans, R.J. Raster: Geographic Data Analysis and Modeling; R Foundation for Statistical Computing: Vienna, Austria, 2016. [Google Scholar]
  75. Bivand, R.; Keitt, T.; Rowlingson, B. Rgdal: Bindings for the “Geospatial” Data Abstraction Library. Available online: https://cran.r-project.org/web/packages/rgdal/ (accessed on 1 January 2021).
  76. Bivand, R.; Lewin-Koh, N. Maptools: Tools for Reading and Handling Spatial Objects. Available online: https://cran.r-project.org/web/packages/maptools/ (accessed on 1 January 2021).
  77. White, T.; van der Ende, J.; Nichols, T.E. Beyond Bonferroni Revisited: Concerns over Inflated False Positive Research Findings in the Fields of Conservation Genetics, Biology, and Medicine. Conserv. Genet. 2019, 20, 927–937. [Google Scholar] [CrossRef] [Green Version]
  78. Pinheiro, J.; Bates, D.; DebRoy, S.; Sarkar, D.; R Core Team. Nlme: Linear and Nonlinear Mixed Effects Models. Available online: https://cran.r-project.org/web/packages/nlme/ (accessed on 1 January 2021).
  79. Hothorn, T.; Bretz, F.; Westfall, P. Simultaneous Inference in General Parametric Models. Biom. J. 2008, 50, 346–363. [Google Scholar] [CrossRef] [Green Version]
  80. Pretzsch, H.; Dieler, J. The Dependency of the Size-Growth Relationship of Norway Spruce (Picea Abies [L.] Karst.) and European Beech (Fagus Sylvatica [L.]) in Forest Stands on Long-Term Site Conditions, Drought Events, and Ozone Stress. Trees 2011, 25, 355–369. [Google Scholar] [CrossRef] [Green Version]
  81. Bose, A.K.; Weiskittel, A.; Wagner, R.G. Occurrence, Pattern of Change, and Factors Associated with American Beech-Dominance in Stands of the Northeastern USA Forest. For. Ecol. Manag. 2017, 392, 202–212. [Google Scholar] [CrossRef] [Green Version]
  82. Čater, M.; Levanič, T. Beech and Silver Fir’s Response along the Balkan’s Latitudinal Gradient. Sci. Rep. 2019, 9, 16269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Costea, C. Codrul Grădinărit; Editura Agro-silvica: Bucharest, Romania, 1962. [Google Scholar]
  84. Čater, M.; Kobler, A. Light Response of Fagus Sylvatica L. and Abies Alba Mill. in Different Categories of Forest Edge–Vertical Abundance in Two Silvicultural Systems. For. Ecol. Manag. 2017, 391, 417–426. [Google Scholar] [CrossRef]
  85. Töchterle, P.; Yang, F.; Rehschuh, S.; Rehschuh, R.; Ruehr, N.K.; Rennenberg, H.; Dannenmann, M. Hydraulic Water Redistribution by Silver Fir (Abies Alba Mill.) Occurring under Severe Soil Drought. Forests 2020, 11, 162. [Google Scholar] [CrossRef] [Green Version]
  86. Garcia-Duro, J.; Manzoni, L.; Arias, I.; Casal, M.; Cruz, O.; Pesqueira Xose, M.; Munoz, A.; Alvarez, R.; Mariot, L.; Bandini, S.; et al. Hidden Costs of Modelling Post-Fire Plant Community Assembly Using Cellular Automata; Springer: Berlin/Heidelberg, Germany, 2018. [Google Scholar]
  87. García-Duro, J.; Cruz, O.; Casal, M.; Reyes, O. Fire as Driver of the Expansion of Paraserianthes Lophantha (Willd.) I. C. Nielsen in SW Europe. Biol. Invasions 2019, 21, 1427–1438. [Google Scholar] [CrossRef]
  88. García-Duro, J.; Álvarez, R.; Basanta, M.; Casal, M. Aplicación de redes bayesianas ás relacións entre especies vexetais despois de incendio forestal e a súa sensibilidade ó tamaño das unidades de mostraxe. In BIOapps2016. Encontro Galaico-Portugués de Biometría, con Aplicación ás Ciencias da Saúde, á Ecoloxía e ás Ciencias do Medio Ambiente; Ginzo, M.J., Pértega, S., Santiago Pérez, M.I., Eds.; SGAPEIO & SPE: Santiago de Compostela, Spain, 2016. [Google Scholar]
  89. Pausas, J.G.; Fernández-Muñoz, S. Fire Regime Changes in the Western Mediterranean Basin: From Fuel-Limited to Drought-Driven Fire Regime. Clim. Chang. 2012, 110, 215–226. [Google Scholar] [CrossRef] [Green Version]
  90. Hubau, W.; Lewis, S.L.; Phillips, O.L.; Affum-Baffoe, K.; Beeckman, H.; Cuní-Sanchez, A.; Daniels, A.K.; Ewango, C.E.N.; Fauset, S.; Mukinzi, J.M.; et al. Asynchronous Carbon Sink Saturation in African and Amazonian Tropical Forests. Nature 2020, 579, 80–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Study area location in Southern Carpathians.
Figure 1. Study area location in Southern Carpathians.
Forests 12 00885 g001
Figure 2. DBH–age relationship in Penteleu 1 ha mixed–virgin forest plot.
Figure 2. DBH–age relationship in Penteleu 1 ha mixed–virgin forest plot.
Forests 12 00885 g002
Figure 3. Dendrochronological biwavelet between species, after 100 Monte Carlo randomizations. Arrows direction indicate that both dendrochronological series are in phase (pointing to the right) or anti-phase (pointing to the left), the black contour indicates ≤0.05 significance, and the transparency, the cone of influence.
Figure 3. Dendrochronological biwavelet between species, after 100 Monte Carlo randomizations. Arrows direction indicate that both dendrochronological series are in phase (pointing to the right) or anti-phase (pointing to the left), the black contour indicates ≤0.05 significance, and the transparency, the cone of influence.
Forests 12 00885 g003
Figure 4. Saturation of the first and second components of the factorial analysis of beech, silver fir and Norway spruce moving correlations between tree ring growth and monthly climate variables (from April to September, respectively: Apr, May, Jun, Jul, Aug, Sep), here sorted by climate variable (i.e., Tmin, Tmax, BAL, PET, Prec).
Figure 4. Saturation of the first and second components of the factorial analysis of beech, silver fir and Norway spruce moving correlations between tree ring growth and monthly climate variables (from April to September, respectively: Apr, May, Jun, Jul, Aug, Sep), here sorted by climate variable (i.e., Tmin, Tmax, BAL, PET, Prec).
Forests 12 00885 g004
Figure 5. Biomass (tonnes) reconstruction based on tree rings (a) at the stand level (continuous) and (b) the relative abundance of the main species for the steepest (dotted) and flattest (continuous) sectors in the Penteleu 1 ha plot (black), beech (orange), silver fir (green), Norway spruce (blue).
Figure 5. Biomass (tonnes) reconstruction based on tree rings (a) at the stand level (continuous) and (b) the relative abundance of the main species for the steepest (dotted) and flattest (continuous) sectors in the Penteleu 1 ha plot (black), beech (orange), silver fir (green), Norway spruce (blue).
Forests 12 00885 g005
Figure 6. Average quantity of biomass (tonnes/ha) in the 22 ecoregions (different line styles) under the four climate change scenarios: RCP 2.6, RCP 4.5, RCP 6.0 and RCP 8.5.
Figure 6. Average quantity of biomass (tonnes/ha) in the 22 ecoregions (different line styles) under the four climate change scenarios: RCP 2.6, RCP 4.5, RCP 6.0 and RCP 8.5.
Forests 12 00885 g006
Figure 7. Average quantity of biomass (tonnes/ha—width lines) in pure and mixed forests from 1990 to 2100 and confidence intervals (shaded areas).
Figure 7. Average quantity of biomass (tonnes/ha—width lines) in pure and mixed forests from 1990 to 2100 and confidence intervals (shaded areas).
Forests 12 00885 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chivulescu, S.; García-Duro, J.; Pitar, D.; Leca, Ș.; Badea, O. Past and Future of Temperate Forests State under Climate Change Effects in the Romanian Southern Carpathians. Forests 2021, 12, 885. https://doi.org/10.3390/f12070885

AMA Style

Chivulescu S, García-Duro J, Pitar D, Leca Ș, Badea O. Past and Future of Temperate Forests State under Climate Change Effects in the Romanian Southern Carpathians. Forests. 2021; 12(7):885. https://doi.org/10.3390/f12070885

Chicago/Turabian Style

Chivulescu, Serban, Juan García-Duro, Diana Pitar, Ștefan Leca, and Ovidiu Badea. 2021. "Past and Future of Temperate Forests State under Climate Change Effects in the Romanian Southern Carpathians" Forests 12, no. 7: 885. https://doi.org/10.3390/f12070885

APA Style

Chivulescu, S., García-Duro, J., Pitar, D., Leca, Ș., & Badea, O. (2021). Past and Future of Temperate Forests State under Climate Change Effects in the Romanian Southern Carpathians. Forests, 12(7), 885. https://doi.org/10.3390/f12070885

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