Next Article in Journal
Improving Otsu Method Parameters for Accurate and Efficient in LAI Measurement Using Fisheye Lens
Next Article in Special Issue
Predicted Future Changes in the Mean Seasonal Carbon Cycle Due to Climate Change
Previous Article in Journal
Genome-Wide Identification of MYB Gene Family in Peach and Identification of MYBs Involved in Carotenoid Biosynthesis
Previous Article in Special Issue
Environmental Response of Tree Species Distribution in Northeast China with the Joint Species Distribution Model
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 15
Issue 7
10.3390/f15071120
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Stand Age and Climate Change Effects on Carbon Increments and Stock Dynamics

by
Elia Vangi
1,2,*,
Daniela Dalmonech
1,3,
Mauro Morichetti
1,
Elisa Grieco
1,
Francesca Giannetti
2,
Giovanni D’Amico
2,
Mahdi (Andre) Nakhavali
4,
Gherardo Chirici
2,3,5 and
Alessio Collalti
1,3
1
Forest Modelling Laboratory, Institute for Agriculture and Forestry Systems in the Mediterranean, National Research Council of Italy (CNR-ISAFOM), Via Madonna Alta 128, 06128 Perugia, Italy
2
geoLAB-Laboratory of Forest Geomatics, Department of Agriculture, Food, Environment and Forestry, Università degli Studi di Firenze, Via San Bonaventura 13, 50145 Firenze, Italy
3
National Biodiversity Future Centre (NBFC), Piazza Marina 61, 90133 Palermo, Italy
4
International Institute for Applied Systems Analysis (IIASA), 2361 Laxenburg, Austria
5
Fondazione per il Futuro delle Città, 50127 Firenze, Italy
*
Author to whom correspondence should be addressed.
Forests 2024, 15(7), 1120; https://doi.org/10.3390/f15071120
Submission received: 24 May 2024 / Revised: 24 June 2024 / Accepted: 25 June 2024 / Published: 27 June 2024
(This article belongs to the Special Issue Modeling Forest Response to Climate Change)
Browse Figures
Versions Notes

Abstract

:
Carbon assimilation and wood production are influenced by environmental conditions and endogenous factors, such as species auto-ecology, age, and hierarchical position within the forest structure. Disentangling the intricate relationships between those factors is more pressing than ever due to climate change’s pressure. We employed the 3D-CMCC-FEM model to simulate undisturbed forests of different ages under four climate change (plus one no climate change) Representative Concentration Pathways (RCP) scenarios from five Earth system models. In this context, carbon stocks and increment were simulated via total carbon woody stocks and mean annual increment, which depends mainly on climate trends. We find greater differences among different age cohorts under the same scenario than among different climate scenarios under the same age class. Increasing temperature and changes in precipitation patterns led to a decline in above-ground biomass in spruce stands, especially in the older age classes. On the contrary, the results show that beech forests will maintain and even increase C-storage rates under most RCP scenarios. Scots pine forests show an intermediate behavior with a stable stock capacity over time and in different scenarios but with decreasing mean volume annual increment. These results confirm current observations worldwide that indicate a stronger climate-related decline in conifers forests than in broadleaves.

1. Introduction

Assessing the quantity of CO2 equivalent stored in forest ecosystems is one of the main goals for implementing the new European Forest Strategy for 2030, a key component of the European Green Deal, to achieve greenhouse gas emission neutrality by 2050. Within this framework, European forest strategies have been geared towards forest-based mitigation plans [1,2], which makes it essential to estimate the carbon sequestration capacity and potential under future climate conditions.
In the near future, Europe and Mediterranean areas will emerge as focal points (‘hot spots’) of climate change, characterized by heightened temperatures and environmental impacts [3,4]. Carbon assimilation and wood production are influenced by environmental conditions (e.g., precipitation, temperature, atmospheric CO2, etc.) and endogenous factors, such as species auto-ecology, age, and hierarchical position within the forest structure. In the past decades, forest ecosystems proved to be crucial net carbon sinks [5,6], likely due to the positive fertilization effects of rising atmospheric CO2 and temperature [7]. However, whether this effect will remain positive or be compensated by other limiting factors is still a matter of debate [8,9,10]. Some studies suggest that the fertilization effect on carbon storage and biomass production fades with forest aging in temperate forests [11,12] since these positive effects cannot continue indefinitely, complicating the picture of the forest response to climate changes even further. This is already the case in Europe, where forest aging and increased disturbances are causing the saturation and decline of the forest carbon sink [9]. Unfortunately, there is not yet a clear strategy to increase the mitigation potentials of forests, and the factors involved are manifold and entangled together [11,13,14].
The need to disentangle the intricate relationships between those factors is even more pressing under climate change. Our current understanding of how future climate will interact with forests of different age classes is particularly limited, especially since only a few studies have explored the relationship between age and the ecosystem’s carbon balance under changing climate conditions [15].
The climate sensitivity of age cohorts is driven, among all, by different access to environmental resources, such as root depth and, therefore, access to water, as well as height, which affects leaf-level water potential and, thus, stomatal conductance [16]. Rooting depth and height jointly affect the tree’s sensitivity to water scarcity, a key environmental driver of change. Future changes in environmental conditions are expected to impact the age spectrum differently [17,18,19].
Since forest age is determined by management practices and 75% of European forests are even aged [20,21], it is crucial to grasp and pin down the role of age in the sensitivity of forest carbon stocks to climate change to guide and inform adaptative forest management.
Process-based forest models enable the exploration of climate change impacts on various age cohorts within the same area, a task difficult to achieve through direct field measurements, which would require decades or more. In this regard, this study examines the ability of different forest age classes under the same future climate conditions to sustain high productivity and carbon stock capacity. To achieve this goal, we employed the ‘Three Dimensional-Coupled Model Carbon Cycle-Forest Ecosystem Module’ (3D-CMCC-FEM) [22,23], simulating undisturbed forests of different cohorts under four climate change scenarios (and including one ‘no climate change’ scenario), from the moderate one (RCP 2.6) up to the most severe one (RCP 8.5) coming from five Earth system models. In this context, carbon stocks and increment were simulated via total carbon woody stocks (TCWS, i.e., the standing woody biomass in MgC ha−1) and the mean annual increment (MAI, in m3 ha−1year−1), which depend mainly on age and long-term processes, such as climate trends.
The primary aim of this research is (i) to explore the direct effects of climate change on the overall carbon storage capacity across various stands, species, and age classes situated in diverse regions of Europe, and (ii) to elucidate the potential influence of forest age on stand dynamics in adapting to forthcoming climate shifts.

2. Materials and Methods

2.1. Study Sites

The study was conducted in three even-aged, previously managed European forest stands (Figure 1): (i) the Boreal Scots pine (Pinus sylvestris L.) forest of Hyytiälä, Finland (FI-Hyy); (ii) the wet temperate continental Norway spruce (Picea abies (L.) H. Karst) forest of Bílý Krìz in the Czech Republic (CZ-BK1); and (iii) the temperate oceanic European beech (Fagus sylvatica L.) forest of Sorø, Denmark (DK-Sor) where the 3D-CMCC-FEM (in different versions) has been already validated in the past [14,24,25]. An overview of the main site’s characteristics is presented in Table 1 and Table 2.
For each site, daily bias-adjusted downscaled climate data from five Earth system models (i.e., HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, GFDL-ESM2M, and NorESM1-M) driven by four representative concentration pathways, namely RCP 2.6, 4.5, 6.0, and 8.5 were available [26,27] (Figure S1). For more detailed information on the study site characteristics and climate data, see [14,24,25,28]. The chosen sites have been selected due to their long monitoring history and the availability of a wide range of data sources for both carbon fluxes and biometric data for model evaluation, as well as bias-corrected climate scenarios for simulations under climate change scenarios from the ISIMIP-PROFOUND initiatives (https://www.isimip.org/, accessed on 1 January 2024) [25,28]. In addition, these stands (i) represent the most common European tree species; (ii) have a current state that is the result of the legacy of past forest management; (iii) are mainly mono-specific and therefore represent interesting «living labs» to study the effects of climate change on single-species and their productivity, reducing confounding effects which otherwise make models struggle to predict forest growth and carbon dynamics (e.g., [29,30]); and (iv) they have already been investigated in the context of climate-smart-forestry silvicultural scenarios [14].

2.2. The Model

The ‘Three Dimensional-Coupled Model Carbon Cycle-Forest Ecosystem Module’ (3D-CMCC-FEM v 5.6 [12,14,22,23,24,31] is a biogeochemical, biophysical, process-based, stand-level forest model. The model is built to simulate carbon, nitrogen, and water cycles in forest ecosystems, even including forest dynamics, under scenarios of climate change and disturbances (e.g., forest management) and parameterized at the species level. Photosynthesis is modeled through the biogeochemical model of Farquhar von Caemmerer and Berry [32], implemented for sun and shaded leaves [33] and parametrized as in Bernacchi et al. [34,35]. Temperature acclimation of leaf photosynthesis to increasing temperature is accounted for following Kattge and Knorr [36]. Autotrophic respiration (RA) is modeled mechanistically by distinguishing the cost of maintaining already existing tissues (RM) and the cost of synthesizing new ones (RG). Maintenance respiration is controlled by the amount of nitrogen (stoichiometrically fixed fraction of live tissues) and temperature. Temperature effects on enzyme kinetics are modeled through a standard Arrhenius relationship but acclimated for temperature as described in Collalti et al. [24]. The net primary productivity (NPP) is the gross primary productivity (GPP) less RA. Not all the annual NPP goes for biomass production since the model considers the non-structural carbon (NSC) pool, an additional seventh C-pool that includes starch and sugars (undistinguished) used to buffer periods of negative carbon balance (when respiration exceeds assimilation, i.e., RA > GPP). Ultimately, the more trees respire, the more NSC is used to sustain metabolism and NSC pool replenishment, and the less NPP and BP there are (and less carbon is stocked). In the extreme case, when and if all NSCs are depleted because of metabolism without being replenished through current photosynthates, the model predicts stand mortality based on the carbon starvation hypothesis [37,38].
The phenological and allocation schemes are all described extensively in Collalti et al. [22,23,39] and Merganičová et al. [39]. The 3D-CMCC-FEM accounts for the ‘age-effect’ in several ways. Ecological theories of the ‘60s describe [40,41], and past and growing pieces of evidence suggest, that stabilization and a further slight decline follow an initial step-wise increase in forest productivity. The causes of such a decline are debated and include a decline in the GPP because of hydraulic limitation [16,42] as well as an increase in RA because of increased respiring biomass [18,19,43]. The 3D-CMCC-FEM accounts for both by including an age modifier [44], which reduces maximum stomatal conductance (and then also GPP) in the Jarvis model and increases RA because of biomass accumulation during forest development.

2.3. Virtual Stands, Model Runs, and Results Evaluation

The 3D-CMCC-FEM was first evaluated under observed climate and field data for GPP and NPPwoody (i.e., the NPP for woody compound; gC m−2 year−1) and the diameter at breast height (DBH) (see ‘Model validation’ paragraph in Supplementary Materials; [12,14]). The model was forced with the modeled climate under different emission scenarios, corresponding to the RCP atmospheric CO2 concentration values for the period 1997 to 2100, ranging from 421.4 µmol mol−1 in the ‘best-case scenario’ (RCP2.6) to 926.6 µmol mol−1 in the ‘worst-case scenario’ (RCP 8.5) coming from the ISIMIP-PROFOUND initiative. For comparison purposes, we forced the forest model with a detrended and repeated meteorology and atmospheric CO2 concentration from 1996 to 2006. The current climate (i.e., no climate change ‘NoCC’) is considered the baseline to compare against climate change scenarios. At the start of the simulations, we created a composite forest matrix (CFM, composed of both measured stand data and “virtual” stand data), following the approach described in Dalmonech et al. [14], to simulate the potential effect of climate stressors on stands of different ages. The 3D-CMCC-FEM has been run at each site to cover the rotation period of each species (from 1997 to 2099) amid the current climate scenario (fixed atmospheric CO2 concentration at the year 2000 of 368.8 μmol mol−1) consisting of detrended and repeated cycles of the present-day observed meteorology from 1996 to 2006 and the Business-as-Usual (BAU) management practices observed at each site (see [28] for the description of BAU applied at each site). Data required to re-initialize the model at every tenth of the rotation length were retrieved from each simulation. Hence, 10 additional stands were chosen for each age in the composite matrix and added to the CFM. This collection of virtual forest stands was used to set different starting stand ages at the present day (aget0) due, ideally, to the past silvicultural practice and climate. Under this framework, a landscape of eleven different stands (in age and their relative C-pools and forest structure) for each site is created. These new stands were used, each running from 2006 to 2099, to assess the impact of climate forcing, as the model has already been shown to be sensitive to forest stand development and the relative standing biomass.
The 3D-CMCC-FEM was initialized with the structural attributes of the newly created stands from 1997, which was the starting year of all simulations and for all stands. Modeled climate change simulations under different RCP-emissions scenarios started to differentiate in 2006 (up to 2100). The simulation runs from the different stand initial conditions, corresponding to different aget0 classes, were carried out without forest management, as we are interested in the direct climate impact on undisturbed forest stand response, avoiding the confounding effects of forest management on the responses (for forest management effects, see [14]). A total of 825 different simulations were performed, as they combined 5 ESMs × 5 RCPs (4 RCPs + 1 current climate scenario) × 11 aget0 classes × 3 sites. Results are reported for MAI (mean annual increment; m3 ha−1 year−1) and TCWS (total carbon woody stocks; MgC ha−1), respectively, as they are considered some of the most representative and fundamental variables in the carbon cycle and forestry. Following the methodology reported [14] (see Table S1 in Supplementary Materials), we evaluated the model forced with the modeled climate. We compared GPP and NPPwoody against eddy covariance estimates and ancillary data for the years 1997–2005 for DK-Sor and FI-Hyy and 2000–2005 for CZ-BK1. We also compared the diameter at breast height (DBH) in all sites with field measures (see Supplementary Materials).

3. Results

Effect of Age Classes and Climate Change on Total Carbon Woody Stock and Increments

Norway spruce at CZ-BK1 shows mean TCWS values ranging between ~70 and ~140 MgC ha−1 under the NoCC scenario over the century, and from ~70 to ~130 MgC ha−1 with a decreasing pattern across all RCPs (Figure 2). In the Norway spruce stands under some ESMs climate forcing (HadGEM2-ES and GFDL-ESM 2M mostly) and under all climate change scenarios, the 3D-CMCC-FEM simulates mortality events for carbon starvation, which increase across stands under gradually warmer climate scenarios and from the oldest stands to the progressively youngest ones.
Under RCP 8.5, all classes show signs of decay at the end of the century. In the youngest aget0 classes, a sharp decrease in MAI was observed (from 8 to 4 m3 ha−1 year−1), while in the older ones, it holds steady to ~3 m3 ha−1 year−1 with a peak around 2075 (Figure 3). At FI-Hyy, younger aget0 classes (14- to 42-year-old) showed the fastest increase in TCWS (reaching 120–130 MgC ha−1 at the end of the century under all scenarios), also reflected in the pattern of MAI. Older aget0 classes showed a more stable trend throughout the simulation (Figure 2), culminating at ~150 MgC ha−1, with MAI steadily declining from 2.5 to 2 m3 ha−1 year−1. In all scenarios, the Scots pine peaked in the 126 and 56 aget0 in TCWS and MAI, respectively. Minor differences were found in mean TCWS between the NoCC and other RCP scenarios, ranging from −1.6% (140-year-old class under RCP 2.6) to +2.8% (14-year-old class under RCP 6.0). At DK-Sor, the results for TCWS show different patterns to other sites, with the highest values ranging between ~240 MgC ha−1 (under NoCC) and ~255 MgC ha−1 (under RCP 8.5) at the end of the century, with the least TCWS under NoCC. The younger classes showed a shallow increase in TCWS during the simulation period, stabilizing at the end of the century, while the older ones kept growing (Figure 4). DK-Sor was the only site where the tightening of the climate conditions caused a positive effect on the MAI, particularly in the younger classes, reversing the trend from negative to positive at the end of the century.
In summary, a positive growth trend of TCWS over time was found in all sites, with the oldest aget0 classes accounting for the most carbon accumulation. Both conifer stands show a plateau with a reduction in growth at the end of the simulation, which is more pronounced and more severe in the warmest climate scenario. Conversely, the beech stands show a positive growth pattern in all scenarios. Similar results were obtained for MAI, where the conifers showed a decreasing trend over the simulation period despite different magnitudes and patterns among aget0 classes. The beech stands exhibited smaller variations among aget0 than among scenarios concerning other sites. In Table 3, we report the mean value of TCW and MAI over the simulation period for each site and climate scenario.

4. Discussion

Age-Dependent Impacts of Climate Change on Forests’ Increment and C-Stocks

The successional stage, represented by forest age, was the main driver controlling C-storage capacity and biomass accumulation, as already known from previous studies [45,46,47], with differences greater among different age cohorts under the same scenario than among different climate scenarios under the same age class [12,14]. The evidence that the carbon budget is mainly controlled by stand age suggests that the effects of climate change on forest cohorts are generally less significant than the effect of age, mainly in terms of the amount of standing biomass. In this sense, age represents multiple and interacting processes, such as tree size [48,49], forest structural traits (canopy closure and LAI), reduction in stomatal conductance [16], and adaptation to specific environmental conditions which, in turn, make it possible to increases the above-ground biomass (AGB) [50]. The model could reproduce the expected behavior of biomass (and thus carbon) accumulation, simulating rapid growth at a young age and saturation for the oldest age class, but not necessarily at the end of the simulation period. Approaching the physiological optima for the species may benefit the biomass synthesis through an augmented photosynthate supply but may eventually increase the respiratory costs of tissue growth and maintenance despite a strong acclimation capacity [18]. High respiratory costs in warm climates with low precipitation regimes, in the older age classes, lead to C-starvation and mortality phenomena, as modelled for the Norway spruce at the CZ-BK1 site. This indicates that the environment has reached its carrying capacity and that competition for limited resources, such as light and water, is excessively high to sustain more biomass in the oldest age classes.
We found different C-accumulation patterns under climate change between coniferous stands and broadleaves. As expected, increasing temperature and changes in precipitation patterns led to a decline in above-ground biomass in spruce stands, especially in the older age classes. On the contrary, the results show that beech forests at DK-Sor will maintain and even increase C-storage rates under most RCP scenarios. Scots pine forests show an intermediate behavior with a stable stock capacity over time and in different scenarios but with decreasing MAI. These results confirm current observations worldwide that indicate a stronger climate-related decline in conifers forests than in broadleaves [51,52,53]. This contrasting response is explained by the different characteristics of the two phyla, in particular, it is due to the temperature adaptation, with generally lower optimum temperature in conifer in addition to its lower sensitivity to the length of the growing season. Similarly, conifers also show lower efficiency in water management because of the shallower root system, which increases the sensitivity to soil aridity and its vulnerability to drought events [54]. Recent studies confirm that growth decline is more pronounced in conifers than in broadleaf, especially beech forests, in the most northern species distribution [55]. Our results confirm the same growth patterns found by recent studies [47,53,56], where broadleaves outperform conifers in productivity, and climate warming will probably exacerbate these opposite growth patterns.
However, despite some studies suggesting that age modulates different adaptation strategies to some extent, it remains unclear whether younger trees may be more affected by climate change than older ones. Bennett et al. [57], in a global analysis, found that droughts consistently had more severe impacts on larger (older) trees, while Wang et al. [11] observed a more substantial and sharper decline in basal area increment in young Korean pine in China. Hogg et al. [58] found that the percentage decrease in biomass growth was not significantly different for young, productive stands compared to older, less productive ones. Our study suggests that warmer and drier conditions and extended growing seasons will affect younger stands more than older ones, but with different trends among species. In particular, MAI will be positively affected in younger beech forests, while it will remain stable in older stands. On the contrary, climate change will strongly impact the growth rate of young conifers stands more than older ones. Older forests tend to be more stable and resilient than younger ones due to their rugged and stable interaction with climate triggers and better responsiveness to environmental changes. The year-to-year climate variability is buffered by larger carbon pools in sapwood and reservoirs in older trees, leading to higher long-term stability than younger trees [12]. In this sense, ages represent the “memory” of the forest to past climate and disturbance regimes, which align the species-specific traits to the environmental conditions in which they grow, creating the niches in which AGB accumulates [52,59].
Despite numerous efforts to decipher forests’ response to climate change, the intricate methods employed by tree species to withstand extreme climates still need to be fully unveiled. Further research exploiting ecophysiological models explicitly accounting for age, tree-ring experiments, and remote sensing will be critical to understanding forest ecosystems’ adaptation strategies to climate change, particularly in the face of rapid warming and extreme disturbances. A better understanding of the interaction between forests and climate can inform better forest management strategies, ultimately dampening the impacts of climate change on forest ecosystems.

5. Limitations

The presented modeling framework has some limitations that should be considered. Firstly, natural disturbances as consequences of climate change, such as windstorms, forest fires, and insect outbreaks, were not simulated. These disturbances cause changes in carbon stocks, nutrients, and soil conditions and contribute to the global release of CO2 in the atmosphere, ultimately leading to increasing temperature and radiation. In contrast, climate extreme events are considered to be already included in the climate scenarios used to force the model and, thus, already accounted for in the model outputs. Additionally, other indirect alterations due to climate change of key drivers, such as nitrogen deposition, phosphorus, or ozone, which can somewhat amplify or reduce our results, were not assessed. Nonetheless, some studies (e.g., [60]) lend credence to the notion that this phenomenon may not be applicable across the board. They highlight the significant responsiveness of various tree species to CO2 fertilization across a wide range of nutrient availability. Finally, no allowance was made for the possibility of species migration to and from the study areas. However, these dynamics may require longer timescales than those simulated in this study.

6. Conclusions

Forest age is confirmed to be a significant factor in determining the carbon storage capacity and biomass accumulation in forest ecosystems, especially in the context of future climate uncertainty. The effects of species, site location, stand-level characteristics, and development stage vary significantly and are contingent on specific factors. We observed that differences in biomass accumulation were more pronounced among different age cohorts than among different climate scenarios within the same age class, with contrasting carbon accumulation patterns under climate change between coniferous and broadleaf forests. Furthermore, our findings shed light on the differential impacts of climate change on younger versus older forest stands. Warmer and drier conditions are projected to affect younger stands more severely, particularly in coniferous forests. However, older forests will likely exhibit greater stability and resilience due to their accumulated carbon pools and enhanced adaptability to environmental changes. While our study provides valuable insights, it also underscores the need for further research to unravel the complex mechanisms by which forests adapt to climate change. This deeper understanding can inform more effective forest management strategies, helping to mitigate the impacts of climate change on forest ecosystems in the future. The varying responses of different tree species highlight the need for tailored management approaches and conservation efforts to enhance the resilience of our forests.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/f15071120/s1, Figure S1: Evaluation of monthly seasonal GPP (gC m−2 month−1) fluxes (left column) and annual (gC m−2 year−1) fluxes (central column) for the sites of Sorø, Bily Kriz, and Hyytiala (rows). Quality-checked and -filtered GPP values evaluated at the sites by the eddy covariance technique are reported as black dots. The shaded area for seasonal values reports the maximum and minimum monthly values recorded in the time series. The shaded area for annual data represents the relative uncertainty bounds. In the third column, a comparison of the predicted annual DBH increment (cm y−1) with site observations at the three sites is reported. Measured data are shown as black dots. Simulated data are reported as continuous lines. Table S1: Performance statistics (coefficient of determination R2, relative root mean square error RMSE (gC m−2 day−1) and Fractional Mean Bias, FMB) computed from monthly seasonal values and annual series of model gross primary productivity, GPP, against eddy covariance estimated and diametric annual increment data, DBH increment, against measured data. Results are reported for simulations forced with local and modeled climate (i.e., ESM) (ESM1, 2, 3, 4, 5 refer to HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, GFDL-ESM 2M, and NorESM1-M, respectively).

Author Contributions

E.V.: Data curation, Formal analysis, Investigation, Writing—original draft, Writing—review & editing; D.D.: Data curation, Formal analysis, Investigation, Writing—review & editing; M.M.: Writing—review & editing; E.G.: Writing—review & editing; F.G.: Writing—review & editing; G.D.: Writing—review & editing; M.N.: Writing—review & editing; G.C.: Writing—review & editing; A.C.: Formal analysis, Investigation, Writing—original draft, Writing—review & editing, Conceptualization. All authors have read and agreed to the published version of the manuscript.

Funding

E.V. and A.C. acknowledge funding by the project “FORESTNAVIGATOR” Horizon Europe research and innovation program under grant agreement No. 101056875. M.M., E.G., F.G. and A.C. acknowledge funding by the project “OptForEU” Horizon Europe research and innovation program under grant agreement No. 101060554. D.D. and A.C. also acknowledge the project funded under the National Recovery and Resilience Plan (NRRP), Mission 4 Component 2 Investment 1.4-Call for tender No. 3138 of 16 December 2021, rectified by Decree n.3175 of 18 December 2021 of Italian Ministry of University and Research funded by the European Union–NextGenerationEU under award Number: Project code CN_00000033, Concession Decree No. 1034 of 17 June 2022 adopted by the Italian Ministry of University and Research, CUP B83C22002930006, Project title “National Biodiversity Future Centre-NBFC”. E.V. and A.C. also acknowledge funding from the MIUR Project (PRIN 2020) “Unraveling interactions between WATER and carbon cycles during drought and their impact on water resources and forest and grassland ecosySTEMs in the Mediterranean climate (WATERSTEM)” (Project number: 20202WF53Z), “WAFER” at CNR (Consiglio Nazionale delle Ricerche), and by PRIN 2020 (cod 2020E52THS)-Research Projects of National Relevance funded by the Italian Ministry of University and Research entitled: “Multi-scale observations to predict Forest response to pollution and climate change” (MULTIFOR, project number 2020E52THS).

Data Availability Statement

The 3D-CMCC-FEM v.5.6 model code is publicly available and can be found on the GitHub platform at: https://github.com/Forest-Modelling-Lab/3D-CMCC-FEM (accessed on 24 June 2024). The raw data supporting the conclusions of this article will be made available by the authors on request.

Acknowledgments

We thank the ISIMIP project and the COST Action FP1304 “PROFOUND” (Towards Robust Projections of European Forests under Climate Change), supported by COST (European Cooperation in Science and Technology) for providing us the climate historical scenarios and site data used in this work. This work used eddy covariance data acquired and shared by the “FLUXNET” community, including these networks: AmeriFlux, AfriFlux, AsiaFlux, CarboAfrica, CarboEurope-IP, CarboItaly, CarboMont, ChinaFlux, Fluxnet-Canada, GreenGrass, ICOS, KoFlux, LBA, NECC, OzFlux-TERN, TCOS-Siberia, and USCCC. The ERA-Interim reanalysis data are provided by ECMWF and processed by LSCE. The FLUXNET eddy covariance data processing and harmonization was carried out by the European Fluxes Database Cluster, AmeriFlux Management Project, and Fluxdata Project of FLUXNET, with the support of CDIAC and ICOS Ecosystem Thematic Center, and the OzFlux, ChinaFlux, and AsiaFlux offices. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the respective climate modeling groups for producing and making available their model output. The U.S. Department of Energy’s Program for Climate Model Diagnosis and Intercomparison at Lawrence Livermore National Laboratory provides coordinating support for CMIP and led the development of software infrastructure in partnership with the Global Organization for Earth System Science Portals.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Favero, A.; Mendelsohn, R.; Sohngen, B.; Stocker, B. Assessing the long-term interactions of climate change and timber markets on forest land and carbon storage. Environ. Res. Lett. 2021, 16, 014051. [Google Scholar] [CrossRef]
  2. Vangi, E.; D’Amico, G.; Francini, S.; Giannetti, F.; Lasserre, B.; Marchetti, M.; McRoberts, R.E.; Chirici, G. The effect of forest mask quality in the wall-to-wall estimation of growing stock volume. Remote Sens. 2021, 13, 1038. [Google Scholar] [CrossRef]
  3. Noce, S.; Collalti, A.; Valentini, R.; Santini, M. Hot spot maps of forest presence in the Mediterranean basin. iForest Biogeosci. For. 2016, 9, 766. [Google Scholar] [CrossRef]
  4. Lionello, P.; Scarascia, L. The relation between climate change in the Mediterranean region and global warming. Reg. Environ. Chang. 2018, 18, 1481–1493. [Google Scholar] [CrossRef]
  5. 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]
  6. He, Y.; Liu, Y.; Lei, L.; Terrer, C.; Huntingford, C.; Peñuelas, J.; Xu, H.; Piao, S. CO2 fertilization contributed more than half of the observed forest biomass increase in northern extra-tropical land. Glob. Chang. Biol. 2023, 29, 4313–4326. [Google Scholar] [CrossRef]
  7. Roebroek, C.T.J.; Duveiller, G.; Seneviratne, S.I.; Davin, E.L.; Cescatti, A. Releasing global forests from human management: How much more carbon could be stored? Science 2023, 380, 749–753. [Google Scholar] [CrossRef] [PubMed]
  8. Duffy, K.A.; Schwalm, C.R.; Arcus, V.L.; Koch, G.W.; Liang, L.L.; Schipper, L.A. How close are we to the temperature tipping point of the terrestrial biosphere? Sci. Adv. 2021, 7, eaay1052. [Google Scholar] [CrossRef]
  9. Nabuurs, G.J.; Verkerk, P.J.; Schelhaas, M.; González-Olabarria, J.R.; Trasobares, A.; Cienciala, E. Climate-Smart Forestry: Mitigation Impacts in Three European Regions; From Science to Policy 6; European Forest Institute: Joensuu, Finland, 2018. [Google Scholar]
  10. Wang, X.; Pederson, N.; Chen, Z.; Lawton, K.; Zhu, C.; Han, S. Recent rising temperatures drive younger and southern Korean pine growth decline. Sci. Total Environ. 2019, 649, 1105–1116. [Google Scholar] [CrossRef]
  11. Gregor, K.; Krause, A.; Reyer, C.P.O.; Knoke, T.; Meyer, B.F.; Suvanto, S.; Rammig, A. Quantifying the impact of key factors on the carbon mitigation potential of managed temperate forests. Carbon Balance Manag. 2024, 19, 10. [Google Scholar] [CrossRef]
  12. Vangi, E.; Dalmonech, D.; Cioccolo, E.; Marano, G.; Bianchini, L.; Puchi, P.F.; Grieco, E.; Cescatti, A.; Colantoni, A.; Chirici, G.; et al. Stand age diversity and climate change affect forests’ resilience and stability, although unevenly. bioRxiv 2023. [Google Scholar] [CrossRef]
  13. Erb, K.; Haberl, H.; Le Noë, J.; Tappeiner, U.; Tasser, E.; Gingrich, S. Changes in perspective needed to forge ‘no-regret’ forest-based climate change mitigation strategies. GCB Bioenergy 2022, 14, 246–257. [Google Scholar] [CrossRef] [PubMed]
  14. Dalmonech, D.; Marano, G.; Amthor, J.; Cescatti, A.; Lindner, M.; Trotta, C.; Collalti, A. Feasibility of enhancing carbon sequestration and stock capacity in temperate and boreal European forests via changes to forest management. Agric. For. Meteorol. 2022, 327, 109203. [Google Scholar] [CrossRef]
  15. Anderson-Teixeira, K.J.; Miller, A.D.; Mohan, J.E.; Hudiburg, T.W.; Duval, B.D.; Delucia, E.H. Altered dynamics of forest recovery under a changing climate. Glob. Chang. Biol. 2013, 19, 2001–2021. [Google Scholar] [CrossRef] [PubMed]
  16. Ryan, M.G.; Binkley, D.; Fownes, J.H. Age-related decline in forest productivity: Pattern and process. Adv. Ecol. Res. 1997, 27, 213–262. [Google Scholar] [CrossRef]
  17. Goulden, M.L.; McMillan, A.M.; Winston, G.C.; Rocha, A.V.; Manies, K.L.; Harden, J.W.; Bond-Lamberty, B.P. Patterns of NPP, GPP, respiration, and NEP during boreal forest succession. Glob. Chang. Biol. 2011, 17, 855–871. [Google Scholar] [CrossRef]
  18. Collalti, A.; Tjoelker, M.G.; Hoch, G.; Mäkelä, A.; Guidolotti, G.; Heskel, M.; Petit, G.; Ryan, M.G.; Battipaglia, G.; Prentice, I.C. Plant respiration: Controlled by photosynthesis or biomass? Glob. Chang. Biol. 2020, 26, 1739–1753. [Google Scholar] [CrossRef] [PubMed]
  19. Collalti, A.; Ibrom, A.; Stockmarr, A.; Cescatti, A.; Alkama, R.; Fernandez-Martínez, M.; Matteucci, G.; Sitch, S.; Friedlingstein, P.; Ciais, P.; et al. Forest production efficiency increases with growth temperature. Nat. Commun. 2020, 11, 5322. [Google Scholar] [CrossRef] [PubMed]
  20. FOREST EUROPE. State of Europe’s Forests 2018; 2018. [Google Scholar]
  21. FOREST EUROPE. State of Europe’s Forests 2020; 2020; Available online: https://foresteurope.org/wp-content/uploads/2016/08/SoEF_2020.pdf (accessed on 23 May 2024).
  22. Collalti, A.; Marconi, S.; Ibrom, A.; Trotta, C.; Anav, A.; D’Andrea, E.; Matteucci, G.; Montagnani, L.; Gielen, B.; Mammarella, I.; et al. Validation of 3D-CMCC Forest Ecosystem Model (v.5.1) against eddy covariance data for ten European forest sites. Geosci. Model Dev. 2016, 9, 479–504. [Google Scholar] [CrossRef]
  23. Collalti, A.; Dalmonech, D.; Vangi, E.; Marano, G.; Puchi, P.F.; Morichetti, M.; Saponaro, V.; Orrico, M.R.; Grieco, E. Monitoring and Predicting Forest Growth and Dynamics; CNR Edizioni: Rome, Italy, 2024; ISBN 978-88-8080-655-4. [Google Scholar] [CrossRef]
  24. Collalti, A.; Trotta, C.; Keenan, T.F.; Ibrom, A.; Bond-Lamberty, B.; Grote, R.; Vicca, S.; Reyer, C.P.O.; Migliavacca, M.; Veroustraete, F.; et al. Thinning can reduce losses in carbon use efficiency and carbon stocks in managed forests under warmer climate. J. Adv. Model. Earth Syst. 2018, 10, 2427–2452. [Google Scholar] [CrossRef]
  25. Mahnken, M.; Cailleret, M.; Collalti, A.; Trotta, C.; Biondo, C.; D’Andrea, E.; Dalmonech, D.; Marano, G.; Mäkelä, A.; Minunno, F.; et al. Accuracy, realism and general applicability of European forest models. Glob. Chang. Biol. 2022, 28, 6921–6943. [Google Scholar] [CrossRef] [PubMed]
  26. Moss, R.; Edmonds, J.; Hibbard, K. The next generation of scenarios for climate change research and assessment. Nature 2010, 463, 747–756. [Google Scholar] [CrossRef] [PubMed]
  27. 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]
  28. Reyer, C.P.O.; Silveyra Gonzalez, R.; Dolos, K.; Hartig, F.; Hauf, Y.; Noack, M.; Lasch-Born, P.; Rötzer, T.; Pretzsch, H.; Meesenburg, H.; et al. The PROFOUND Database for evaluating vegetation models and simulating climate impacts on European forests. Earth Syst. Sci. Data 2020, 12, 1295–1320. [Google Scholar] [CrossRef]
  29. Vacchiano, G.; Magnani, F.; Collalti, A. Modeling Italian forests: State of the art and future challenges. iForest 2012, 5, 113–120. [Google Scholar] [CrossRef]
  30. Collalti, A.; Perugini, L.; Santini, M.; Chiti, T.; Nolè, A.; Matteucci, G.; Valentini, R. A process-based model to simulate growth in forests with complex structure: Evaluation and use of 3D-CMCC Forest Ecosystem Model in a deciduous forest in Central Italy. Ecol. Model. 2014, 272, 362–378. [Google Scholar] [CrossRef]
  31. Farquhar, G.; von Caemmerer, S.; Berry, J. A biogeochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 1980, 149, 78–90. [Google Scholar] [CrossRef] [PubMed]
  32. de Pury, D.G.G.; Farquhar, G.D. Simple scaling of photosynthesis from leaves to canopies without the errors of big-leaf models. Plant Cell Environ. 1997, 20, 537–557. [Google Scholar] [CrossRef]
  33. Bernacchi, C.J.; Singsaas, E.L.; Pimentel, C.A.R.L.O.S.; Portis Jr, A.R.; Long, S.P. Improved temperature response functions for models of Rubisco-limited photosynthesis. Plant Cell Environ. 2001, 24, 253–259. [Google Scholar] [CrossRef]
  34. Bernacchi, C.J.; Calfapietra, C.A.R.L.; Davey, P.A.; Wittig, V.E.; Scarascia-Mugnozza, G.E.; Raines, C.A.; Long, S.P. Photosynthesis and stomatal conductance responses of poplars to free-air CO2 enrichment (PopFACE) during the first growth cycle and immediately following coppice. New Phytol. 2003, 159, 609–621. [Google Scholar] [CrossRef]
  35. Kattge, J.; Knorr, W. Temperature acclimation in a biochemical model of photosynthesis: A reanalysis of data from 36 species. Plant Cell Environ. 2007, 30, 1176–1190. [Google Scholar] [CrossRef] [PubMed]
  36. McDowell, N. Mechanism linking drought, hydraulics, carbon metabolism, and vegetation mortality. Plant Physiol. 2011, 155, 1051–1059. [Google Scholar] [CrossRef]
  37. Rowland, L.; da Costa, A.C.L.; Galbraith, D.R.; Oliveira, R.S.; Binks, O.J.; Oliveira, A.A.R.; Pullen, A.M.; Doughty, C.E.; Metcalfe, D.B.; Vasconcelos, S.S.; et al. Death from drought in tropical forests is triggered by hydraulics not carbon starvation. Nature 2015, 528, 119–122. [Google Scholar] [CrossRef] [PubMed]
  38. Collalti, A.; Biondo, C.; Buttafuoco, G.; Maesano, M.; Caloiero, T.; Lucà, F.; Pellicone, G.; Ricca, N.; Salvati, R.; Veltri, A.; et al. Simulation, calibration and validation protocols for the model 3D-CMCC-CNR-FEM: A case study in the Bonis’ watershed (Calabria), Italy. For. J. Silvic. For. Ecol. 2017, 14, 247–256. [Google Scholar] [CrossRef]
  39. Merganičová, K.; Merganič, J.; Lehtonen, A.; Vacchiano, G.; Sever, M.Z.O.; Augustynczik, A.L.D.; Grote, R.; Kyselová, I.; Mäkelä, A.; Yousefpour, R.; et al. Forest carbon allocation modelling under climate change. Tree Physiol. 2019, 39, 1937–1960. [Google Scholar] [CrossRef] [PubMed]
  40. Kira, T.; Shidei, T. Primary production and turnover of organic matter in different forest ecosystems of the western Pacific. Jpn. J. Ecol. 1967, 17, 70–87. [Google Scholar]
  41. Odum, E.P. The Strategy of Ecosystem Development: An understanding of ecological succession provides a basis for resolving man’s conflict with nature. Science 1969, 164, 262–270. [Google Scholar] [CrossRef]
  42. Tang, J.; Luyssaert, S.; Richardson, A.D.; Kutsch, W.; Janssens, I.A. Steeper declines in forest photosynthesis than respiration explain age-driven decreases in forest growth. Proc. Natl. Acad. Sci. USA 2014, 111, 8856–8860. [Google Scholar] [CrossRef] [PubMed]
  43. Mäkelä, A.; Valentine, H. The ratio of NPP to GPP: Evidence of change over the course of stand development. Tree Physiol. 2000, 21, 1015–1030. [Google Scholar] [CrossRef]
  44. Landsberg, J.J.; Waring, R.H. A generalised model of forest productivity using simplified concepts of radiation-use efficiency, carbon balance and partitioning. For. Ecol. Manag. 1997, 95, 209–228. [Google Scholar] [CrossRef]
  45. Anderson, K.J.; Allen, A.P.; Gillooly, J.F.; Brown, J.H. Temperature-dependence of biomass accumulation rates during secondary succession. Ecol. Lett. 2006, 9, 673–682. [Google Scholar] [CrossRef]
  46. Cook-Patton, S.C.; Leavitt, S.M.; Gibbs, D.; Harris, N.L.; Lister, K.; Anderson-Teixeira, K.J.; Briggs, R.D.; Chazdon, R.L.; Crowther, T.W.; Ellis, P.W.; et al. Mapping carbon accumulation potential from global natural forest regrowth. Nature 2020, 585, 545–550. [Google Scholar] [CrossRef]
  47. Anderson-Teixeira, K.J.; Herrmann, V.; Morgan, R.B.; Bond-Lamberty, B.; Cook-Patton, S.C.; E Ferson, A.; Muller-Landau, H.C.; Wang, M.M. Carbon cycling in mature and regrowth forests globally. Environ. Res. Lett. 2021, 16, 053009. [Google Scholar] [CrossRef]
  48. Ouyang, S.; Xiang, W.; Wang, X.; Xiao, W.; Chen, L.; Li, S.; Sun, H.; Deng, X.; Forrester, D.I.; Zeng, L.; et al. Effects of stand age, richness and density on productivity in subtropical forests in China. J. Ecol. 2019, 107, 2266–2277. [Google Scholar] [CrossRef]
  49. Ullah, F.; Gilani, H.; Sanaei, A.; Hussain, K.; Ali, A. Stand structure determines aboveground biomass across temperate forest types and species mixture along a local-scale elevational gradient. For. Ecol. Manag. 2021, 486, 118984. [Google Scholar] [CrossRef]
  50. Lee, Y.-J.; Park, G.-E.; Lee, H.-I.; Lee, C.-B. Stand age-driven tree size variation and stand type regulate aboveground biomass in alpine-subalpine forests, South Korea. Sci. Total. Environ. 2024, 915, 170063. [Google Scholar] [CrossRef] [PubMed]
  51. Hlásny, T.; Barka, I.; Kulla, L.; Bucha, T.; Sedmák, R.; Trombik, J. Sustainability of forest management in a Central European mountain forest: The role of climate change. Reg. Environ. Chang. 2015, 17, 65–77. [Google Scholar] [CrossRef]
  52. 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]
  53. Pretzsch, H.; del Río, M.; Arcangeli, C.; Bielak, K.; Dudzinska, M.; Forrester, D.I.; Klädtke, J.; Kohnle, U.; Ledermann, T.; Matthews, R.; et al. Forest growth in Europe shows diverging large regional trends. Sci. Rep. 2023, 13, 15373. [Google Scholar] [CrossRef] [PubMed]
  54. Krejza, J.; Cienciala, E.; Světlík, J.; Bellan, M.; Noyer, E.; Horáček, P.; Štěpánek, P.; Marek, M.V. Evidence of climate-induced stress of Norway spruce along elevation gradient preceding the current dieback in Central Europe. Trees 2021, 35, 103–119. [Google Scholar] [CrossRef]
  55. del Castillo, E.M.; Zang, C.S.; Buras, A.; Hacket-Pain, A.; Esper, J.; Serrano-Notivoli, R.; Hartl, C.; Weigel, R.; Klesse, S.; de Dios, V.R.; et al. Climate-change-driven growth decline of European beech forests. Commun. Biol. 2022, 5, 163. [Google Scholar] [CrossRef] [PubMed]
  56. Gong, X.; Yuan, D.; Zhu, L.; Li, Z.; Wang, X. Long-term changes in radial growth of seven tree species in the mixed broadleaf-Korean pine forest in Northeast China: Are deciduous trees favored by climate change? J. For. Res. 2024, 35, 70. [Google Scholar] [CrossRef]
  57. Bennett, A.C.; McDowell, N.G.; Allen, C.D.; Anderson-Teixeira, K.J. Larger trees suffer most during drought in forests worldwide. Nat. Plants 2015, 1, 15139. [Google Scholar] [CrossRef] [PubMed]
  58. Hogg, E.H.; Michaelian, M.; Hook, T.I.; Undershultz, M.E. Recent climatic drying leads to age-independent growth reductions of white spruce stands in western Canada. Glob. Chang. Biol. 2017, 23, 5297–5308. [Google Scholar] [CrossRef]
  59. Reyer, C.P.O. Forest Productivity Under Environmental Change—A Review of Stand-Scale Modeling Studies. Curr. For. Rep. 2015, 1, 53–68. [Google Scholar] [CrossRef]
  60. Terrer, C.; Vicca, S.; Hungate, B.; Phillips, R.P.; Prentice, I.C. Mycorrhizal association as a primary control of the CO2 fertilization effect. Science 2016, 353, 72–74. [Google Scholar] [CrossRef]
Forests 15 01120 g001
Figure 1. Test site locations in Europe. The red “+” in the overwiew panel represent the sitre locations.
Figure 1. Test site locations in Europe. The red “+” in the overwiew panel represent the sitre locations.
Forests 15 01120 g001
Forests 15 01120 g002
Figure 2. Modeled total carbon woody stock (TCWS) (MgC ha−1) for age classes at the three sites in all scenarios along the simulation period (2006–2099). Lines represent the moving average of 10 years. The solid line corresponds to the real stand, while the dotted lines correspond to the virtual ones. The shaded area represents two standard deviations from the mean predictions with the results from the five ESMs’ climate change scenarios.
Figure 2. Modeled total carbon woody stock (TCWS) (MgC ha−1) for age classes at the three sites in all scenarios along the simulation period (2006–2099). Lines represent the moving average of 10 years. The solid line corresponds to the real stand, while the dotted lines correspond to the virtual ones. The shaded area represents two standard deviations from the mean predictions with the results from the five ESMs’ climate change scenarios.
Forests 15 01120 g002
Forests 15 01120 g003
Figure 3. Modeled mean annual increment (MAI) (m3 ha−1 year−1) for age classes at the three sites in all scenarios along the simulation period (2006–2099). Lines represent the moving average of 10 years. The solid line corresponds to the real stand, while the dotted lines correspond to the virtual ones. The shaded area represents two standard deviations from the mean predictions with the results from the five ESMs climate change scenarios.
Figure 3. Modeled mean annual increment (MAI) (m3 ha−1 year−1) for age classes at the three sites in all scenarios along the simulation period (2006–2099). Lines represent the moving average of 10 years. The solid line corresponds to the real stand, while the dotted lines correspond to the virtual ones. The shaded area represents two standard deviations from the mean predictions with the results from the five ESMs climate change scenarios.
Forests 15 01120 g003
Forests 15 01120 g004
Figure 4. Boxplot of modeled total carbon woody stock (TCWS) (left, MgC ha−1) and mean annual increment (MAI) (right, m3 ha−1 year−1) for age classes at the three sites in the four RCPs scenarios compared to the NoCC (no climate change) scenario. Boxplots with thick borders correspond to the real stand. Lines are fitted throughout the median of the values of the variables using a generalized additive model.
Figure 4. Boxplot of modeled total carbon woody stock (TCWS) (left, MgC ha−1) and mean annual increment (MAI) (right, m3 ha−1 year−1) for age classes at the three sites in the four RCPs scenarios compared to the NoCC (no climate change) scenario. Boxplots with thick borders correspond to the real stand. Lines are fitted throughout the median of the values of the variables using a generalized additive model.
Forests 15 01120 g004
Table 1. Overview of the main site characteristics provided for each forest site. Years of obs. refers to the first and last year of measurement; the temporal resolution of measurement is annual. For the stand values (DBH, height, BA, age, and tree density), the range corresponds to the first and last field measurement according to the years of obs. Column. DBH = diameter at breast height; BA = basal area.
Table 1. Overview of the main site characteristics provided for each forest site. Years of obs. refers to the first and last year of measurement; the temporal resolution of measurement is annual. For the stand values (DBH, height, BA, age, and tree density), the range corresponds to the first and last field measurement according to the years of obs. Column. DBH = diameter at breast height; BA = basal area.
NameSpeciesLatLongAspectElevationSlopeYears of Obs.DBH
(cm)		Height (m)BA (m2 ha−1)AgeTree Density (ha−1)
Bílý KrìzPicea abies49.318.3218087512.51997–20158.16–20.476.26–15.2610.33–36.9616–342408–1252
HyytiäläPinus sylvestris61.8524.2918018521995–201115.89–20.5812.61–18.6212.64–18.3334–50870–684
SorøFagus sylvatica55.4911.64-4001994–201728.99–48.2524.23–31.1518.50–29.7662–87407–199
Table 2. Yearly averages of the daily maximum temperature (Tmax), daily minimum temperature (Tmin), daily mean temperature (Tmean), annual precipitation sum (P), daily mean relative humidity (RH), daily mean air pressure (AP), and annual sum of global radiation (R, direct + diffuse shortwave radiation) for each of the sites. The column “Years” indicates the data’s acquisition year and the period the average values refer to.
Table 2. Yearly averages of the daily maximum temperature (Tmax), daily minimum temperature (Tmin), daily mean temperature (Tmean), annual precipitation sum (P), daily mean relative humidity (RH), daily mean air pressure (AP), and annual sum of global radiation (R, direct + diffuse shortwave radiation) for each of the sites. The column “Years” indicates the data’s acquisition year and the period the average values refer to.
SiteSourceYearsTmax (°C)Tmean (°C)Tmin (°C)P (mm)RH (%)R (J cm−2)
Bílý KrìzLocal2000–200811.57.363.81434.5681.99378 774.86
HyytiäläLocal1996–20147.44.361.13604.0177.95309 628.86
SorøLocal1996–201210.668.265.91760.5282.95360 687.83
Table 3. Mean values of total carbon woody stock (TCWS) and mean annual increment (MAI) over the simulation period (2006–2099) for each scenario and age class. CZ-BK1 = Bílý Krìz; FI-Hyy = Hyytiälä; DK-Sor = Sorø.
Table 3. Mean values of total carbon woody stock (TCWS) and mean annual increment (MAI) over the simulation period (2006–2099) for each scenario and age class. CZ-BK1 = Bílý Krìz; FI-Hyy = Hyytiälä; DK-Sor = Sorø.
ScenarioScenario
NoCCRCP 2.6RCP 4.5RCP 6.0RCP 8.5NoCCRCP 2.6RCP 4.5RCP 6.0RCP 8.5
AgeTCWS (Mg Cha−1)MAI (m3 ha−1)
CZ-BK11272.2269.9265.0470.4265.934.894.754.574.804.63
1675.7173.6073.2074.1071.885.034.864.674.914.72
2486.5384.4784.0685.0282.704.734.604.414.634.45
3697.4695.3795.0096.1393.644.304.194.014.224.07
48101.7699.6999.3391.0099.583.813.713.713.743.68
60110.19108.17107.8899.34107.003.583.513.503.533.47
72118.19116.00115.64117.06114.553.423.353.343.373.31
84121.56119.38119.04120.58117.503.183.123.113.143.08
96120.19118.12117.98119.70116.902.882.832.832.852.80
108126.78121.26124.28117.57112.332.742.702.702.712.66
120145.50124.43142.41135.02127.842.722.692.682.692.65
FI-Hyy1466.4266.6267.3968.4966.981.911.931.951.991.94
2879.8379.0079.7381.0579.372.212.192.212.252.20
36104.34102.42103.17104.75102.863.203.143.163.213.15
4289.8288.5589.3590.7989.082.352.312.332.372.33
56118.76116.62117.57119.20117.363.203.143.173.213.17
70119.29117.43118.61120.25118.652.912.862.892.922.89
84125.92123.95125.19126.85125.272.762.722.742.772.74
98133.33131.06132.29134.12132.332.662.622.642.672.64
112141.29138.87140.00141.92140.192.572.532.552.572.55
126148.26145.43146.69148.63146.912.492.452.472.492.47
140155.77153.11154.50156.52154.792.472.432.452.472.45
DK-Sor1498.53107.89110.48109.94113.322.462.722.792.772.86
28110.94120.06122.68122.18125.532.602.862.922.902.99
42113.56122.68125.31124.82128.162.542.782.852.832.92
56115.97125.07127.73127.23130.582.522.762.822.802.89
70116.15125.14127.79127.30130.622.472.692.752.732.81
76120.65130.45133.16132.55136.052.763.003.063.043.13
84124.68133.99136.65136.14139.532.562.782.842.822.90
98135.87145.41148.06147.49150.932.823.023.073.063.13
112154.13164.32166.58165.71168.773.003.213.253.233.30
126166.84176.33178.38177.63180.323.063.243.283.263.31
140170.49180.08181.84181.26184.052.913.083.113.093.14
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vangi, E.; Dalmonech, D.; Morichetti, M.; Grieco, E.; Giannetti, F.; D’Amico, G.; Nakhavali, M.; Chirici, G.; Collalti, A. Stand Age and Climate Change Effects on Carbon Increments and Stock Dynamics. Forests 2024, 15, 1120. https://doi.org/10.3390/f15071120

AMA Style

Vangi E, Dalmonech D, Morichetti M, Grieco E, Giannetti F, D’Amico G, Nakhavali M, Chirici G, Collalti A. Stand Age and Climate Change Effects on Carbon Increments and Stock Dynamics. Forests. 2024; 15(7):1120. https://doi.org/10.3390/f15071120

Chicago/Turabian Style

Vangi, Elia, Daniela Dalmonech, Mauro Morichetti, Elisa Grieco, Francesca Giannetti, Giovanni D’Amico, Mahdi (Andre) Nakhavali, Gherardo Chirici, and Alessio Collalti. 2024. "Stand Age and Climate Change Effects on Carbon Increments and Stock Dynamics" Forests 15, no. 7: 1120. https://doi.org/10.3390/f15071120

APA Style

Vangi, E., Dalmonech, D., Morichetti, M., Grieco, E., Giannetti, F., D’Amico, G., Nakhavali, M., Chirici, G., & Collalti, A. (2024). Stand Age and Climate Change Effects on Carbon Increments and Stock Dynamics. Forests, 15(7), 1120. https://doi.org/10.3390/f15071120

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