Skip to main content
  • Research Paper
  • Published:

CO2 fertilization, transpiration deficit and vegetation period drive the response of mixed broadleaved forests to a changing climate in Wallonia

Abstract

Key message

The change in forest productivity was simulated in six stands in Wallonia (Belgium) following different climate scenarios using a process-based and spatially explicit tree growth model. Simulations revealed a strong and positive impact of the CO 2 fertilization while the negative effect of the transpiration deficit was compensated by longer vegetation periods. The site modulated significantly the forest productivity, mainly through the stand and soil characteristics.

Context

Forest net primary production (NPP) reflects forest vitality and is likely to be affected by climate change.

Aims

Simulating the impact of changing environmental conditions on NPP and two of its main drivers (transpiration deficit and vegetation period) in six Belgian stands and decomposing the site effect.

Methods

Based on the tree growth model HETEROFOR, simulations were performed for each stand between 2011 and 2100 using three climate scenarios and two CO2 modalities (constant vs time dependent). Then, the climate conditions, soils and stands were interchanged to decompose the site effect in these three components.

Results

In a changing climate with constant atmospheric CO2, NPP values remained constant due to a compensation of the negative effect of increased transpiration deficit by a positive impact of longer vegetation periods. With time-dependent atmospheric CO2, NPP substantially increased, especially for the scenarios with higher greenhouse gas (GHG) emissions. For both atmospheric CO2 modalities, the site characteristics modulated the temporal trends and accounted in total for 56 to 73% of the variability.

Conclusion

Long-term changes in NPP were primarily driven by CO2 fertilization, reinforced transpiration deficit, longer vegetation periods and the site characteristics.

1 Introduction

Forests affect the climate in various and complex ways, through biophysical and biogeochemical effects. They can have a warming or a cooling effect depending on the way they modify albedo and evapotranspiration compared with other land covers (Bonan 2008; Stocker et al. 2013). In addition, forest ecosystems decrease the concentration of atmospheric carbon dioxide (CO2) by carbon storage in the soil and in tree biomass (Myhre et al. 2013; Le Quéré et al. 2018). On the other hand, most of the processes involved in forest ecosystem functioning are climate sensitive (Parmesan and Yohe 2003; Kint et al. 2012; Charru et al. 2017). Climate changes could therefore seriously affect forest dynamics and the provision of ecosystem services (Temperli et al. 2012; Mina et al. 2017) that in turn affect climate (Seidl et al. 2014; Thom et al. 2017).

Many environmental manipulation experiments and monitoring studies have highlighted the effects of atmospheric CO2 concentration, air temperature and soil water availability on forest net primary production (NPP). The fertilizing effect of atmospheric CO2 was pointed out from free-air or chamber CO2 enrichment experiments (Ainsworth and Long 2005; Ainsworth and Rogers 2007; Thompson et al. 2017). Norby et al. (2005) found that NPP increased by 20-25% when the CO2 level was elevated to 550 ppm. On the long run, however, this increase was constrained by nutrient availability (Norby et al. 2010; Warren et al. 2015).

Furthermore, global remote sensing and local observational studies have shown that the increase in air temperature tends to advance the budburst date and therefore lengthens the vegetation period (Menzel et al. 2006; Jeong et al. 2011; Richardson et al. 2013; Park et al. 2016; Flynn and Wolkovich 2018) and favours NPP (Baldocchi 2008; Dragoni et al. 2011; Richardson et al. 2013; Fu et al. 2017). Moreover, respiration and photosynthesis, which have an opposite effect on biomass production, are both stimulated by warmer temperatures (Zhang et al. 2017).

Concerning water availability, it is well known that soil water deficit during the vegetation period triggers stomatal closure and reduces carbon assimilation. In addition, it reduces the mineralization, uptake and transport of nutrients (Bréda et al. 2006; Lisar et al. 2012; Osakabe et al. 2014). Syntheses and meta-analyses of rainfall experiments generally report decreases in photosynthesis, NPP, aboveground biomass and soil respiration (Wu et al. 2011; Paschalis et al., 2020). During prolonged water stress periods, tree mortality may occur due to hydraulic failure or carbon starvation (McDowell et al. 2008; Adams et al. 2017). For instance, a particularly drastic rainfall exclusion experiment conducted during a 5-year period in a Mediterranean broadleaf forest (Q. ilex, A. unedo and P. latifolia) led to a 15% decrease in soil water content, lowered the aboveground biomass production by 83% and increased mortality rate by 46% (Ogaya and Peñuelas 2007). Similarly, different studies have highlighted the same trends but with a lower effect during the 2003 drought (Ciais et al. 2005; Bréda et al. 2006; Granier et al. 2007).

The Intergovernmental Panel on Climate Change (IPCC) has adopted different greenhouse gas (GHG) emission scenarios called Representative Concentration Pathway (RCPs). Two commonly used scenarios are RCP4.5 and RCP8.5 and feature CO2 concentrations of 540 and 940 ppm in 2100, respectively. Regarding Atlantic Europe (see Fig. 1 in Jacob et al. 2014), the climate projections for the two abovementioned scenarios at the end of the twenty-first century are characterized by an increase in air temperature (between + 1.3 °C and + 4.2 °C), changes in annual precipitation (− 1% to + 9%), in high-temperature extremes, in drought events and by a decrease in summer precipitation (− 5 to − 25%) compared with the 1971–2000 period (Jacob et al. 2014; Kovats et al. 2014; Jacob et al. 2018).

The projected climate changes are likely to alter the forest dynamics (Lindner et al. 2014). However, the resulting trend is difficult to estimate due to the opposite and interacting effects of climate on the processes regulating forest productivity. To estimate tree growth in unprecedented conditions, process-based models (PBM), which integrate knowledge from in-situ experiments, are recommended since they allow the integration of a wide range of effects (Pretzsch et al. 2015). However, the future forest evolution will not be uniform but is more likely to be site dependent. Indeed, the different components of the site effect—climate, soil and stand—can modify the forest response to changing conditions (Holdridge 1967; Tucker et al. 2008; Steenberg et al. 2015) and therefore, PBM must account for them.

The response of a given tree species to climate change depends on the local climate (first component of the site effect). For sites at the upper limit of the temperature distribution range of a tree species, an increase in temperature would be detrimental for the development of this tree species while it could be positive in colder sites. The same reasoning also holds for a rainfall reduction, which could induce water stress more frequently on drier sites but could be beneficial on wetter ones. In Wallonia, an increase in temperature of 2 °C would confine European beech (Fagus Sylvatica (L.)) in its tolerance limits while sessile oak (Quercus petraea (Matt.) Liebl.) would remain in optimal growing conditions (Petit et al. 2017). The climate change effect can then be modulated by the soil properties. Indeed, greater soil depth and water holding capacity decrease the negative impact of drought events on NPP (Phillips et al. 2016). Finally, the stand characteristics generate differences in the forest response to climate change. Different effects constitute the stand effect: the tree species identity and diversity and the stand structure and density. The tree species identity effect simply depicts that all tree species do not exhibit the same functional traits such as phenological timing (Vitasse et al. 2009; Cole and Sheldon 2017) or drought sensitivity (Leuzinger et al. 2005; Scherrer et al. 2011). Consequently, they display different phenological and stomatal responses to a rise of temperature (Primack et al. 2009) and a change in soil water content or in CO2 levels (Medlyn et al. 2001; Jonard et al. 2011; Raftoyannis and Radoglou 2002). Additionally, when these tree species coexist in the same stand, the response to climate change is even more difficult to predict since it can be affected by complementarity, facilitation and selection effects (Grossiord 2019; Bello et al. 2019a). For instance, Anderegg et al. (2018) observed that mixing species with a high diversity of hydraulic traits decreases their sensitivity to drought. On the other hand, an intensive water use by some tree species mixtures could lead more rapidly to water stress (Pretzsch and Biber 2016). The stand structure modifies the response to climate change since trees respond differently depending on their size, age or social status. Dominant trees seem to be more sensitive to climatic stress than dominated ones, which can be explained by a higher vulnerability to hydraulic stress (McDowell and Allen 2015) and a stronger exposition to atmospheric conditions (Bennett et al. 2015). When soil water availability decreases, small trees keep their stomata open for longer since the position of their crown within the canopy limits the evaporative demand (Carl et al. 2018). Finally, tree density also affects the stand response to climate change as it is positively related to stand leaf area index and, consequently, to evapotranspiration (Bello et al. 2019b). During drought events, thinned broadleaved stands generally exhibit a lower growth decrease due to reduced competition for water among remaining trees (Sohn et al. 2016). Understanding the effect of stand characteristics on forest response to climate change is essential since the forester can adapt it while the adaptability of the soil is much more limited and is practically inapplicable for the climate. In conclusion, sites subject to similar changes in climate conditions can display very different responses due to their climate, soil and stand characteristics. This was already pointed out based on observations by Boisvenue and Running (2006) and can be found in simulation studies (Reyer 2015).

In order to produce realistic projections of NPP for the twenty-first century, one must take into account the future climate under various GHG emission scenarios, the soil properties and the stand characteristics. Yet, most of the currently existing PBM are stand-scale models accounting for the climate and soil effects but only partly integrating the spatial and structural complexity of the stand. However, structurally complex stands (uneven aged and mixed) tend to be favoured by the foresters in order to improve the resilience of their forests (Bauhus et al. 2013; DeRose and Long 2014). To simulate correctly such stands, spatially explicit and individual-based approaches are required (Seidl et al. 2005).

In this study, we will address the question of how and to what extent changing climate and CO2 conditions will impact forest productivity in six contrasting and structurally complex stands in Wallonia and assess how the site components and thinning operations modulate the response. To do so, we will use HETEROFOR (Jonard et al. 2020; de Wergifosse et al. 2020a), a process-based tree growth model running at the individual level with different climate projections based on three GHG emission scenarios.

Wallonia is within a climatic zone that deserves more studies on the impact of climate change on forests since this effect is quite uncertain there for several reasons. Under the current climate conditions, Wallonia is located at the transition between areas where forest growth is limited by temperature, and consequently by the length of the growing period, in the North and by water availability in the South (Bastrup-Birk et al. 2016). In the future, warmer temperatures are expected all over Europe, a rainfall increase in the Northern temperature-limited areas and inversely a precipitation decrease in the Southern water-limited areas (IPCC 2013; Kovats et al. 2014). Therefore, while most simulation studies agree that forest productivity should increase in Boreal forests and decrease or remain constant in Mediterranean area (Reyer et al. 2014; Reyer 2015), to what extent Walloon forest productivity will be constrained by water availability, on average and during extreme events, is not a consensual issue. In addition, an increased water stress could be partly compensated by a longer vegetation period in the future.

More precisely, the objectives of the paper are

  1. I.

    To simulate the temporal changes in the net primary production of six broadleaved stands in Wallonia (Belgium) and in two of its main drivers (transpiration deficit and vegetation period) under various GHG scenarios.

  2. II.

    To differentiate the long-term trend from the inter-annual NPP variations and to evaluate the part of the NPP variability explained by transpiration deficit, vegetation period and atmospheric CO2 concentration.

  3. III.

    To assess to which extent these temporal trends are modulated by the site and how thinning affects transpiration deficit.

  4. IV.

    To decompose the site effect in its three main components (climate, soil and stand) and see how they are accounted for by transpiration deficit and vegetation period.

2 Material and methods

2.1 Site description

For the simulation study, six long-term monitoring plots installed in sessile oak (Quercus petraea (Matt.) Liebl.) and European beech (Fagus sylvatica (L.)) stands were selected as case studies representative of the broadleaved forest in Wallonia (Belgium). Three plots are located in an experimental site in Baileux (50° 01′ N, 4° 24′ E) and have been monitored since 2001. The three other plots are part of the level II network of ICP Forests since 1998 (Ferretti and Fischer 2013) and are located in Louvain-la-Neuve (50° 41′ N, 4° 36′ E), Chimay (50° 07′ N, 5° 34′ E) and Virton (49° 32′ N, 5° 34′ E).

2.1.1 Stand characteristics

The three plots in Baileux were installed to study how species mixtures influence the forest ecosystem functioning (Jonard et al. 2006, 2007, 2008; André et al. 2008a, 2008b, 2010, 2011). Two are located in nearly monospecific stands dominated either by sessile oak (hereafter called Baileux-oak) or by European beech (Baileux-beech) and the third one is in a balanced mixture of both tree species (Baileux-mixed) (Table 1). The stand of Chimay originates from coppice-with-standards and is dominated by mature sessile oaks with a hornbeam understory. In Louvain-la-Neuve and Virton, European beech is the main tree species and is mixed with oaks. In Virton, other broadleaved species (maple, wild cherry, ash, hornbeam) are also present. All these stands cover a wide range in tree size (girth in 1999 or 2001: from 22.6 cm for hornbeam in Baileux-mixed to 190 cm for sessile oak in Virton), in stand density (basal area in 1999 or 2001: from 18.4 to 30.0 m2 ha−1) and in leaf area index (LAI) with values extending from 3.96 to 6.93 m2 m−2. Except for the Louvain-la-Neuve plot installed in an old even-aged beech forest, the study stands display complex structure in terms of species composition and tree age (Table 1).

Table 1 Stand characteristics for the main tree species and for the whole stands

2.1.2 Soil properties

The soils in Baileux and Chimay are cambisols with moder humus (FAO soil taxonomy) developed in a bedrock of sandstone and shales mixed with loess deposits while the soil in Louvain-la-Neuve is an abruptic luvisol with moder humus formed in loamy loess deposits from the interglacial period. Finally, the soil of Virton is a calcaric cambisol with mull humus; it originates from the weathering of a hard limestone bedrock (Table 2).

Table 2 Soil and meteorological (2001–2016 period) characteristics of the different study sites

This large range of soils is reflected in the soil texture. According to the USDA taxonomy, the soil of Baileux is silty clay loam and that of Chimay is intermediate between silty clay loam and clay loam. The soil in Louvain-la-Neuve is silty loam due to its high silt content (62–75%). The highest clay fraction was registered in Virton and it is therefore a clayey soil (Table 2). The soil texture influenced the soil hydraulic properties. In Baileux and Louvain-la-Neuve, the high silt content ensures good drainage. In Chimay, the presence of inflating clay favours the appearance of a perched water table near the surface during winter and drought cracks in the warm and dry period. In Virton, in spite of elevated clay content, the existence of faults in the bedrock enables an efficient drainage.

The stoniness varies a lot among plots. Baileux-beech, Baileux-mixed and Virton have the higher coarse fraction (> 50% in the deep horizons) while the soil in Louvain-la-Neuve contains hardly any large stones. The coarse fraction in the soils of Baileux-oak and Chimay is intermediate (between 30 and 45% at the bottom of the profile). These differences in soil texture and coarse fraction among sites lead to a great diversity of extractable water reserve on a 1.6-m depth. For example, the extractable water reserve in Louvain-la-Neuve is almost three times greater than that of Baileux-mixed and Virton (Table 2). A same soil depth of 1.6 m was retained for all site as it allows to account for most (if not all) the rooting zone and for the sake of comparability. The impact of using a same depth for all sites is limited since water uptake only occurs in the horizon where roots are present. In addition, the lower horizons of some sites are very stony and therefore the associated extractable water is limited.

2.1.3 Climate

Although entire Belgium being characterized by the same climate type (temperate maritime), the four sites of our study are all located in different bioclimatic zones (Van der Perre et al. 2015), which highlights their climate diversity (Table 2). Louvain-la-Neuve is characterized by the warmest (11.0 °C) and driest (818 mm) climate. Despite their proximity, Baileux and Chimay display differences in terms of annual rainfall with 1075 mm in Baileux and 940 mm in Chimay probably due to the more elevated location of Baileux (between 305 and 312 m a.s.l.) than Chimay (260 m a.s.l.). However, yearly average temperatures are similar in both locations (9.8 °C in Baileux and 9.7 °C in Chimay). Finally, Virton features high precipitation values (1060 mm) and an intermediate yearly mean temperature of 9.9 °C (Table 2). Reference crop evapotranspiration was calculated for each site according to the FAO method (Allen et al. 1998) and showed a low variability (between 712 and 745 mm year−1) that can mainly be attributed to air temperature differences.

2.2 Simulation experiments

2.2.1 Model description and performances

For the simulations, we used the individual-based, spatially explicit and process-based model HETEROFOR that has been implemented in the Capsis simulator (Dufour-Kowalski et al. 2012) and is especially convenient to simulate the evolution of structurally complex stands. Hereafter, we present a brief overview of the model functioning limited to the description of the options retained for the simulation experiments carried out in this study. For a more in-depth description, we refer the reader to Jonard et al. (2020) and de Wergifosse et al. (2020a).

To initialize the model, the user must provide different files: a tree species parameter file, an inventory (or stand) file with the tree coordinates and dimensions (tree circumference, total height, height to crown base, height of largest crown extension and crown radius in four cardinal directions), a soil file with the physical properties of each soil horizon (thickness, bulk density, coarse fraction, sand, silt, clay and organic carbon contents and fine root proportion) and a meteorology (or climate) file with hourly data for radiation, air temperature, precipitation, relative humidity and wind speed. After the initialization phase and at the end of each year, HETEROFOR calls the phenology routine that provides daily information concerning the foliage state for the coming year. From a fixed date, a sum of warm temperature is daily accumulated until reaching a threshold, which triggers the budburst and then, similarly, the progressive leaf expansion. From mid-summer, cold temperatures are accumulated until reaching a threshold, which triggers leaf yellowing occurring at a rate proportional to the photoperiod decrease. The leaf falling rate is calculated based on wind speed and accelerated in case of frost events. Phenology can be calculated at the species level or the individual scale to account for the extended vegetation period of understory trees. In this study, we chose to calculate phenology at the species level. Once phenology is fixed, HETEROFOR calculates the proportion of solar radiation intercepted by each tree using the SAMSARALIGHT library based on a ray tracing approach (Courbaud et al. 2003). From the photosynthetically active radiation (PAR) absorbed per unit leaf area and the soil water potential updated hourly with a water balance routine, the gross primary production (GPP) of each tree is estimated hourly with the photosynthesis model of the CASTANEA library (Farquhar et al. 1980; Dufrêne et al., 2005). The individual NPP is then obtained by using a NPP to GPP ratio to take the growth and maintenance respiration into account. NPP is first allocated to foliage and fine roots by ensuring a functional balance and then to structural components using allometric equations, which allows deriving tree dimensional growth. The water balance routine accounts for rainfall partitioning in throughfall, stemflow and interception (André et al. 2008a), for tree transpiration and evaporation from foliage, bark and soil using the Penman-Monteith equation (Monteith 1965), for root water uptake (Couvreur et al. 2012) and for soil water movements based on the Darcy law. This routine can be calculated at the stand or individual scale but calculation time considerably increases when the individual option is selected and consequently, we selected the stand scale in this study.

HETEROFOR was evaluated on the same sites as those used in this study. The evaluation of tree growth conducted in Jonard et al. (2020) demonstrated the model ability to reproduce size-growth relationships and individual radial growth. In the same study, the simulated GPP was related to the NPP reconstructed from tree growth measurement and this relationship displayed high Pearson’s correlation coefficients (0.89 and 0.91 for sessile oak and European beech, respectively). Regarding water balance evaluated in de Wergifosse et al. (2020a), the model satisfactorily simulated the soil water content temporal dynamics (correlation coefficients between simulations and observations ranging from 0.83 to 0.95 according to the site) and the individual transpiration (0.85 and 0.89 for oak and beech, respectively). Finally, the budburst model previously described, a one-phase model based on the Uniforc model (Chuine 2000), which simulates the forcing period, has been chosen among two other budburst models (two-phase models accounting for the forcing and chilling periods) implemented in HETEROFOR as it reproduced best the inter-annual variability. With this budburst model, vegetation period was on average simulated with a RMSE of 6.7 days (unpublished data).

The parameters needed to initialize the model are those described in the two model description papers (Table 2 in Jonard et al. 2020 and Table 1 in de Wergifosse et al. 2020a). Only for Virton, a higher NPP to GPP ratio is used in the species parameter file to account for differences not considered by the model (probably due to a higher site fertility). This ratio was fixed to 0.6 and 0.75 for sessile oak and European beech, respectively. Yet, this parameterization difference is only applied in the first simulation experiment described below. Indeed, one of the objective of this first experiment is to estimate the future productivity of the site, for which as many site characteristics as possible must be integrated while the second simulations aimed at decomposing the site effects in its components without being affected by the parametrization.

2.2.2 First simulation experiment to highlight the climate change impact

One-year simulations were performed for different periods in order to assess the impact of climate change on stand NPP median value and variance, on transpiration deficit and on vegetation period (objectives I and II stated in Section 1) and to evaluate how site and thinning affected these changes (objective III). More specifically, each year, a new simulation was launched starting from the same initial stand. Stand characteristics were therefore reinitialized each year keeping thereby the focus on the climate impact, in contrast to multi-year simulations which could have given rise to diverging stand characteristics with time. Climate projections generated according to three Radiative Concentration Pathway (RCP) scenarios (see below) were used to run simulations on the 2011–2100 period while the 1976–2005 period (called “historical” period) was used as reference for comparisons with RCP scenarios.

One of the major uncertainties when simulating long-term forest productivity is whether or not the positive response of forest to rising CO2 concentration can persist. Indeed, it has been shown that the induced productivity gain may progressively be reduced when other factors such as nutrient availability become limiting (Körner 2006; Norby et al. 2010). In order to cover the range of possible tree responses, from a perfect acclimation to rising atmospheric CO2 to no acclimation at all, the set of simulations was launched considering either a constant atmospheric CO2 concentration (380 ppm) or time-dependent CO2 concentrations corresponding to the RCP scenarios (Reyer et al. 2014). The first case can be seen as a response to increasing CO2 concentration totally constrained by other limiting factors and the second as never constrained. Moreover, holding CO2 concentration constant allowed us to have a better insight into the other climate effects.

The two simulation types were also run for each monitoring plot after applying a virtual thinning to reduce stand basal area by 25% and test the immediate thinning impact. The selection of thinned trees was made among the pool of trees that were effectively cut in each plot during the monitoring period. When the past thinning operations were insufficient to reach 25%, additional trees were randomly selected and removed.

2.2.3 Second simulation experiment to highlight and decompose the site effect

To further investigate the site effect on the NPP variability (objective IV stated in Section 1), a similar set of 1-year simulations was ran for the historical period (1976–2005) by combining the climate, soil and stand input files of all the monitoring plots according to a full factorial design (6 soil types × 6 stands × 4 climates × 30 years). The simulations were performed for a constant atmospheric CO2 concentration in order to limit the number of variation factors. Then, the simulations were repeated for the 2071–2100 period considering the RCP8.5 scenario to test whether the NPP variance decomposition is affected by climate change.

2.2.4 Climate projections

As a basis, the climate projections of the global climate model (GCM) CNRM-CM5 were used here. These global simulations were also included in the Coupled Model Intercomparison Project (CMIP5) on which the IPCC bases most of its conclusions (Fifth Assessment Report: IPCC 2013; IPCC special report on 1.5 C: Masson-Delmotte 2018). However, the horizontal resolution of CNRM-CM5 is 1.4° ( 155 km), which did not allow us to make any distinctions between our study plots. In a first step, the CNRM-CM5 projections were therefore downscaled over the European domain using the Regional Climate Model (RCM) ALARO-0 (Giot et al. 2016; Termonia et al. 2018) following the guidelines of the Coordinated Regional Downscaling Experiment (CORDEX; Giorgi et al. 2009; Jacob et al. 2014). This dynamic downscaling consisted in using ALARO-0 over 50-km resolution areas forced at their boundaries by projections of CNRM-CM5. In a second downscaling step, the simulations over Europe with 50-km resolution were downscaled to a 4-km resolution over Belgium (Rummukainen 2010).

The meteorological variables that served as input for HETEROFOR include hourly values of the longwave and shortwave radiations, air temperature, surface temperature, rainfall, specific humidity and wind speed. Finally, relative humidity was calculated based on temperature, specific humidity and atmospheric pressure. All values were taken at the grid points closest to the four sites for the historical period (1976–2005) and for the 2011–2100 period according to three RCP scenarios: RCP2.6, RCP4.5, RCP8.5. The scenario names depict the increase in radiative forcing in 2100 relative to pre-industrial levels (+ 2.6 W m−2, + 4.5 W m−2, + 8.5 W m−2). The climate projections should be considered as sensitivity experiments. In other words, the climate changes rather than the absolute climate values are of importance as the model climatology (during the historical period) is known to differ from the observed one (Maraun and Widmann 2018). However, for our case, there were important positive model biases in rainfall ranging from 7 to 35% when compared with the observed values at the considered sites. A bias correction was therefore performed (Maraun and Widmann 2018). More specifically, for air and soil temperatures, data were corrected according to

$$ {x}_{\mathrm{corr}\_\mathrm{t}}={x}_{\mathrm{simul}\_\mathrm{t}}+\left(\overline{x_{\mathrm{obs}}}-\overline{x_{\mathrm{simul}}}\right) $$
(1)

with

xcorr_t, the corrected value of a variable at time t

xsimul_t, the variable simulated by the regional model at time t

\( \overline{x_{obs}} \) and \( \overline{x_{simul}} \), the average observed and simulated values for the period 2001–2016.

This method is, however, not suitable for variables which cannot take negative values. For these variables (radiation, rainfall, relative humidity and wind speeds), data were corrected using a multiplicative scaling

$$ {x}_{\mathrm{corr}\_\mathrm{t}}={x}_{\mathrm{simul}\_\mathrm{t}}\ \frac{\overline{x_{\mathrm{obs}}}}{\overline{x_{\mathrm{simul}}}} $$
(2)

The same correction was applied to the three RCP scenarios using \( \overline{x_{\mathrm{obs}}} \) and \( \overline{x_{\mathrm{simul}}} \) based on the period 2001–2016. The average corrected mean air temperature, rainfall and reference evapotranspiration are presented in Table 3 for the three RCP scenarios during the 2071–2100 period and also during the historical period.

Table 3 Mean annual rainfall, air temperature and reference crop evapotranspiration for the historical period (1976–2005) and the various RCP scenarios (2071–2100). Mean values for the vegetation period (fixed from the 1st of May to the 30th of September) are indicated between parentheses. For each site and climatic variable, same letters indicate that the scenarios are not statistically different (Student’s t test)

2.2.5 Model simulation analysis

The HETEROFOR model generates many fluxes and stocks of carbon and water as outputs. For this study, we focus on the actual and potential tree transpiration (obtained without considering any limitation from soil water) to determine the transpiration deficit, the daily foliage status of each tree species to calculate the vegetation period length and the yearly NPP to characterize forest productivity (de Wergifosse et al. 2020b).

From the daily foliage status, the yearly vegetation period length was defined as the number of days between the day the green leaf proportion reaches 50% (budburst period) until the day it drops below this threshold (yellowing and then falling periods). Annual stand NPP values (gC m−2) were simply the sum of individual tree NPP (gC) divided by the stand area (m2), with NPP derived from GPP after accounting for the growth and maintenance respirations (see Jonard et al. 2020 for details). The annual transpiration deficit was calculated for each tree as the difference between actual and potential transpiration (in L). Then, individual transpiration deficits were summed and divided by the stand area to obtain a transpiration deficit in mm. As described in de Wergifosse et al. (2020a), the stomatal conductance is considered as decreasing exponentially with the soil water potential in the model. Therefore, this difference depicts the transpiration deficit induced by the soil water limitation. Using the vegetation period and the transpiration deficit is interesting since these variables are sufficiently integrative to summarize the model functioning but not too general so that we can disentangle the effect of two main mechanisms through which climate change affects forest ecosystems functioning (phenology and water).

In order to compare stand NPP, transpiration deficit and vegetation period (objective I) among RCP scenarios and time periods (1976–2005 for past climate and 2011–2040, 2041–2070 and 2071-2100 for the future climate), we used two different statistical tests to assess whether the distributions were significantly different. An unpaired Mann-Whitney test (Wilcoxon 1945) was performed when the two periods were not related (e.g. for the comparison of the 1976–2005 and 2041–2070), while a Wilcoxon signed-rank test (Wilcoxon 1945) was used when comparing RCP scenarios for a same period. These tests were chosen given the non-normality of the investigated variables. The Wilcoxon signed-rank test was also adopted to assess the effect of thinning on transpiration deficit. In order to test the equality of variance among distributions, a Levene test was performed as it is less sensitive to non-normality than other commonly used tests.

To differentiate the long-term effect of climate change on NPP from that of the inter-annual climate variability (objective II) while taking the site effect into account (objective III), a linear mixed model was fitted on the simulated stand NPPs of the first simulation experiment including both thinning modalities. For a same location, the thinned and unthinned stands were considered as two different sites in the linear mixed models We have chosen to use linear mixed models to account for the correlation structure of the simulated dataset and to avoid an overparameterization of the model (for parsimony reasons). Some factors were important to consider to estimate their relative importance in explaining the NPP variability and to represent correctly the correlation structure of the data. However, we did not need to know accurately the value taken by each level of these factors, which, therefore, were considered as random and characterized with a limited number of parameters (one parameter per factor instead of one per factor level). In contrast, we wanted to accurately quantify the effects of other factors which were considered as fixed. For this reason, in the first mixed model, the temporal trend of each RCP scenario (time × scenario) is considered as a fixed effect and the site (site) and its effect on the temporal trend (site × time × scenario) as random factor effects. This model was applied for both atmospheric CO2 modalities.

$$ \mathrm{NPP}\sim {\left[\mathrm{time}\times \mathrm{scenario}\right]}_{\mathrm{fixed}}+{\left[\mathrm{site}+\mathrm{site}\times \mathrm{time}\times \mathrm{scenario}+\upepsilon \right]}_{\mathrm{random}} $$
(3)

The continuous variable characterizing the time effect is the number of years since 2011. In this way, no effect of the RCP scenario is considered in 2011, which allowed us to avoid considering the scenario as a main effect.

Besides, we adjusted another linear mixed model containing yearly vegetation period (VP), transpiration deficit (TD) and atmospheric CO2 concentration (CO2) as fixed effects in addition to the effects already considered in the previous model (Eq. (3)) in order to assess the extent to which these three variables accounted for the site effect, the long-term trend and the inter-annual variability:

$$ \mathrm{NPP}\sim {\left[\mathrm{time}\times \mathrm{scenario}+\mathrm{VP}+\mathrm{TD}+\mathrm{CO}2\right]}_{\mathrm{fixed}}+{\left[\mathrm{site}+\mathrm{site}\times \mathrm{time}\ \mathrm{x}\ \mathrm{scenario}+\upepsilon \right]}_{\mathrm{random}} $$
(4)

This model was applied for both atmospheric CO2 modalities, except that the atmospheric CO2 concentration was logically not considered for the modality with constant atmospheric CO2 concentration.

Using the outputs of the second simulation experiment, a linear mixed model was applied to decompose the site effect in its climate, soil and stand components (objective IV). Three random factors were used to characterize the site components.

$$ \mathrm{NPP}\sim {\left[\mathrm{climate}+\mathrm{soil}+\mathrm{stand}+\upepsilon \right]}_{\mathrm{random}} $$
(5)

Finally, to estimate how transpiration deficit and vegetation period accounted for the three components of the site effect, we fitted a linear mixed model containing these two drivers as fixed effects in addition to the effects considered in Eq. (5).

$$ \mathrm{NPP}\sim {\left[\mathrm{VP}+\mathrm{TD}\right]}_{\mathrm{fixed}}+{\left[\mathrm{climate}+\mathrm{soil}+\mathrm{stand}+\upepsilon \right]}_{\mathrm{random}} $$
(6)

For all the effects in the various models, the partial R2 was calculated as the difference between the R2 of the model with and without the considered effect. This method assumes that the effects are independent. As, in reality, this not always the case, the sum of the partial R2 can be lower than the R2 of the full model. All the figures, statistical tests and linear mixed models were realized using R Studio software (RStudio Team 2015).

3 Results

3.1 Objective I: Climate change impact on NPP, transpiration deficit and vegetation period

Net primary production

The differences in NPP between the RCP scenarios were generally non-significant when the atmospheric CO2 concentration was kept constant. For this modality, the only significant difference with the historical period occurred between 2041 and 2070 for the RCP4.5 and 8.5 scenarios and remained limited: an increase of 3 and 5%, respectively (Fig. 1a). The site-by-site examination of the NPP projections revealed that the only sites with a significant positive effect of the RCP scenarios were Baileux-oak (RCP4.5 between 2071 and 2100) and Chimay (RCP4.5 and 8.5 for the 2041–2070 period and RCP8.5 during the 2071–2100 period) (Fig. 4).

Fig. 1
figure 1

NPP comparisons among RCP scenarios and with the historical period (1976–2005) for three 30-year periods considering all the sites together, with constant (a) and time-dependent (b) atmospheric CO2 concentrations. The horizontal line corresponds to the median, the box ends indicate the upper and lower quartiles and the whiskers show the values above and below these quartiles within 1.5 interquartile. Common letters indicate that the distributions are non-significantly different according to a paired Wilcoxon signed-rank (between scenarios of the same period) or an unpaired Mann-Whitney (between scenarios of different periods) tests

For the simulations accounting for the time-dependent atmospheric CO2, NPP increased significantly over time, especially for the scenarios with the higher CO2 emission levels in 2100. Upon comparison with the historical period, NPP in 2071–2100 increased by 9%, 20% and 34% for RCP2.6, 4.5 and 8.5, respectively (Fig. 1b). The trends were similar for the different sites taken individually (Fig. 5).

The impact of climate change on NPP cannot be based solely on the change in its median value. The variability is a key component of the temporal evolution as well. However, as depicted by the boxplot width and whisker length of Fig. 1 and confirmed by the Levene tests, no consistent increase in NPP variability was observed in our simulations.

Transpiration deficit

The results obtained for transpiration deficit were identical for the constant and time-dependent atmospheric CO2 concentration since the way stomatal conductance for water was calculated does not account for the atmospheric CO2 concentration effect. Therefore, only those of the constant atmospheric CO2 were presented (Fig. 2a). Transpiration deficit did not change in comparison with the historical period during the 2011–2040 period. During the next period, the RCP2.6 and 8.5 scenarios displayed a significant increase in transpiration deficit of 24% and 19%, respectively. For this period, the only inter-scenario difference occurred between RCP2.6 and 4.5. During the last period, all scenarios were different from each other and the RCP4.5 and 8.5 scenarios showed respectively a 21% and 42% increase in transpiration deficit compared with the historical period (Fig. 2a).

Fig. 2
figure 2

Transpiration deficit (a), vegetation period of sessile oak (b) and European beech (c) comparisons among RCP scenarios and with the historical period (1976–2005) for three 30-year periods considering all the sites together, with constant atmospheric CO2 concentrations. The horizontal line corresponds to the median, the box ends indicate the upper and lower quartiles and the whiskers show the values above and below these quartiles within 1.5 interquartile. Common letters indicate that the distributions are non-significantly different according to a paired Wilcoxon signed-rank (between scenarios of the same period) or an unpaired Mann-Whitney (between scenarios of different periods) tests

The analysis of the temporal change in actual and potential transpiration enabled us to get a better insight in the origin of the transpiration deficit. It appeared that both variables were characterized by an increasing trend but with a more pronounced one for potential transpiration. The augmentation in actual transpiration ranged from 2.7 to 10.2% during the last period (2071–2100) while the rise in potential transpiration varied from 3.4 to 19.8% (Fig. 6).

Vegetation period

As phenology is not CO2 dependent in the HETEROFOR model, no distinction was made between atmospheric CO2 modalities. Even though the length of the vegetation period differed between oak and beech (206 days for oak and 209 days for beech), its temporal change due to climate change was very similar for both tree species. The vegetation period length increased significantly with time for RCP4.5 and 8.5 while, for RCP2.6, it peaked between 2041 and 2070 before returning to the 2011–2040 level in 2071–2100. The increase of the vegetation period length (calculated with regards to the last period) amounted to 0.6, 1.8 and 5.3 days per decade for RCP2.6, 4.5 and 8.5, respectively (Fig. 2 b and c).

3.2 Objective II: Long-term trend and inter-annual variations of NPP

When a constant atmospheric CO2 concentration was considered for the model described by Eq. (3), no mean temporal trend in NPP was detected (time × scenario) though the random effect associated with this trend (site × time × scenario) explained 34% of the variability. This means that the temporal trend oscillated around 0, being slightly positive in some sites and slightly negative in others. The rest of the site effect accounted for 39% of the NPP variability and the unexplained variability amounted to 27% (Table 4 (a)).

Table 4 Estimate, standard error, P value and partial R2 of the parameters of the mixed linear model (Eq. (3): NPP~time*scenario + site + site*time*scenario) used to differentiate the temporal trend from the inter-annual variability in NPP (gC m−2) for all sites, periods and RCP scenarios considering constant (a) and time-dependent (b) atmospheric CO2 concentration

For the simulations with changing atmospheric CO2 concentrations, we observed a significant positive temporal trend dependent on both the RCP scenario (time × scenario) and on the site (site × time × scenario). These two effects accounted for 22% and 24% of the NPP variability, respectively. The remaining site effect explained 32% of the NPP variability while the unexplained variability was slightly lower than the simulations with constant CO2, accounting for 23% of the NPP variability (Table 4 (b)).

3.3 Objective III: Influence of atmospheric CO2 concentration, transpiration deficit, vegetation period, site and thinning on NPP

Relative importance of atmospheric CO2 concentration, transpiration deficit, vegetation period and site in explaining NPP variability

For both atmospheric CO2 modalities, when comparing the linear mixed models containing only the temporal trend and the site effect (Eq. (3)) with those including also transpiration deficit, vegetation period and atmospheric CO2 concentration (Eq. (4)), one can observe that these variables explained the temporal trend (time × scenario) and its modulation by the site (site × time × scenario) and all the inter-annual variability (residuals). In addition, their inclusion slightly lowered the variability associated with the rest of the site effect. When atmospheric CO2 concentration was held constant, the transpiration deficit effect was negative with a partial R2 of 0.58 while the vegetation period effect was positive with a partial R2 of 0.025 (Table 5 (a)). For the time-dependent atmospheric CO2 simulations, the transpiration deficit and vegetation period effects had the same sign as for the constant atmospheric CO2 modality but their partial R2 was lower (0.43 and 0.016 respectively) since the atmospheric CO2 concentration also explained a significant part of the variability (partial R2 = 0.22, Table 5 (b)).

Table 5 Estimate, standard error, P value and partial R2 of the parameters of the mixed linear model (Eq. (4): NPP~time*scenario + VP + TD + CO2 + site + site*time*scenario) used to explain the NPP (gC m−2) variability based on transpiration deficit, vegetation period and CO2 (only for time-dependent CO2 modality) for all the sites, periods and RCP scenarios considering constant (a) and time-dependent (b) atmospheric CO2 concentration

Thinning effect on transpiration deficit

Transpiration deficit was much higher in beech dominated (between 205 and 267 mm) than in oak-dominated stands (between 75 and 104 mm) and this was observed after thinning as well (Fig. 3). Thinning significantly decreased transpiration deficit in all situations (P value < 0.001). For the historical period, RCP2.6, 4.5 and 8.5, the transpiration deficit in the thinned oak-dominated stands was 30, 34, 32 and 39 mm lower than in the unthinned ones (Fig. 3a) while in beech-dominated stands, this decrease amounted to 65, 82, 76 and 91 mm, respectively (Fig. 3b). The relative decrease was, however, more similar between the two tree species with a decrease between 35 and 40% for oak and between 32 and 34% for beech.

Fig. 3
figure 3

Thinning effect on transpiration deficit for the historical period (1976–2005) and the three RCP scenarios (2011–2100) considering a constant atmospheric CO2 concentration. Results are shown separately for oak-dominated stands (Baileux-oak and Chimay) (a) and beech-dominated stands (Baileux-beech, Louvain-la-Neuve and Virton) (b). Significance of the Wilcoxon test evaluating the thinning effect by climate scenario is represented as NS (non-significant), *p < 0.05, **p < 0.01 and ***p < 0.001

3.4 Objective IV: Decomposition of the site effect in its climate, soil and stand components

The second simulation experiment, aiming at decomposing the site effect in its climate, soil and stand components, was performed for two contrasted periods and RCP scenarios (i.e. the 1977–2005 historical period vs the 2071–2100 period for the RCP8.5 scenario). As the results were quite similar among periods, only the results for the historical period are presented. The stand and soil partial R2 were close and amounted to 0.321 and 0.264, respectively, and were much higher than that of the climate effect (0.016) (Eq. (5)). Introducing the transpiration deficit and the vegetation period in the model (Eq. (6)) accounted for nearly the entire climate (94%) and soil (98%) effects but only partly for the stand effect (7%). The consideration of these drivers also strongly reduced the unexplained part of the variability (from 40 to 7%). The transpiration deficit had a negative effect while the impact of vegetation period was positive, with a much higher partial R2 for transpiration deficit (0.536) than for vegetation period (0.057) (Table 6).

Table 6 Estimate, standard error, P value and partial R2 of the parameters of (a) the mixed linear model (Eq. (5): NPP~stand + soil + climate) and (b) the mixed linear model (Eq. (6): NPP~VP + TD + stand + soil + climate) used to decompose the site effect in its climate, soil and stand components in order to explain the NPP (gC m−2) variability for all the sites during the historical period, using simulations with constant atmospheric CO2 concentration

4 Discussion

4.1 What are the possible evolutions of broadleaved forest NPP according to different climate projections?

At first glance, our results appear quite clear and easy to interpret. When the atmospheric CO2 level was held constant, no long-term changes were observed but upon changing atmospheric CO2, NPP increased up to 34% for the 2071–2100 period. However, using a process-based model with many different outputs gave us the opportunity to understand more deeply the underlying mechanisms.

Based on the bias-corrected downscaled climate projections of Alaro-0 for our four sites in Wallonia, mean air temperature is expected to be significantly higher during the 2071–2100 period compared with the historical period (from 0.7 to 3.4 °C) while mean annual rainfall would remain stable or increase a bit. During the vegetation period, rainfall would slightly decrease but this decrease would be significant only for Virton under RCP4.5 scenario (Table 3). According to our simulations, these changing climate conditions would substantially increase the vegetation period (Fig. 2 b and c) since the budburst is triggered earlier when late winter temperatures are warmer while leaf yellowing occurs later under warm conditions. The decrease of the vegetation period length between 2041–2070 and 2071–2100 for RCP2.6 scenario might look surprising but it is simply generated by the scenario that projects a temperature peak around the middle of the twenty-first century and a subsequent progressive decrease. Some limitations concerning the phenology date simulations are discussed hereafter. First, the budburst model is a one-phase model that only accounts for the accumulation of warm temperature to trigger budburst while it is commonly accepted that a chilling period is a prerequisite for the start of forcing period (corresponding to the endodormancy break) and then budburst. However, when no observations of endodormancy break are available and the species considered are not located at the margin of their species distribution area, one-phase models are often preferred to more complex ones (Chuine et al. 2016). In addition, our approach does not account for the impact of photoperiod on budburst, which can become significant when chilling requirements are not met (Vitasse and Basler 2013; Pletsers et al. 2015), the influence of late frost or water stress on the leaf development and senescence (Sanz-Perez and Castro-Diez 2010; Fu et al. 2014; Morin and Chuine 2014; Xie et al. 2018). However, as shown in the review of Piao et al. (2019), the modelling of these second order processes is extremely difficult and inaccurate because the interactions between these factors are still poorly known, observations are available only for a few phenological stages and environmental modifying experiments have not been conducted to disentangle such a complexity.

In addition to the strong impact of climate change on the vegetation period, the model simulates that the transpiration deficit could be moderately reinforced through an increased evapotranspiration (Fig. 2a). As the transpiration deficit was a stronger NPP driver than the vegetation period (Tables 5 and 6), their opposite effects on the long-term trends in NPP had more or less the same magnitude and were compensating. However, the long-term NPP evolution was slightly affected by the site: it tended to be slightly negative (positive) in sites with a low (high) soil holding capacity and for drought-sensitive (-tolerant) tree species such as beech (oak) (Fig. 4). Furthermore, by reducing the transpiration deficit, thinning contributed to make the climate change effect on NPP more positive. As the oak stomatal sensitivity to soil drying is lower than that of beech, the simulated impact of climate change was rather positive on oak and null on beech-dominated stands, which is consistent with simulation studies estimating that oak competitiveness could exceed that of beech under the projected future climate conditions in temperate European forests (Kint et al. 2012; Mette et al. 2013; Zimmermann et al. 2013; Rubio-Cuadrado et al. 2018).

With our simulations, we were unable to distinguish between the transpiration deficit under constant and time-dependent atmospheric CO2 concentration since the water balance is calculated before photosynthesis in HETEROFOR and the stomatal conductance for water does not depend on atmospheric CO2 concentration contrary to that for CO2 (Dufrêne et al., 2005; Jonard et al. 2020; de Wergifosse et al. 2020a). Therefore, our transpiration simulations are more reliable for the constant atmospheric CO2 modality. For the time-dependent CO2 modality, tree transpiration is probably slightly overestimated so that the increase in the transpiration deficit could still be lower than that simulated with HETEROFOR which is already limited. According to our simulations, water stress on Walloon broadleaved forests should not drastically increase in the future.

As stated in the description of the first simulation experiment (Section 2.2.2), the objective of running two similar sets of simulations with constant and changing atmospheric CO2 concentrations was to define a possible range of NPP change depending on the way the fertilizing effect of CO2 is constrained by other limiting factors. As soil water limitation was taken into account in the model, the more obvious constraint would come from the soil nutrient supply. The increased atmospheric CO2 could potentially have no effect in sites where trees suffer from severe nutrient deficiencies. The CO2 fertilization, on the other hand, would be manifested in full for stands with an optimal mineral nutrition (Oren et al. 2001; Fernandez-Martinez et al. 2014a). On average, this CO2 fertilizing effect seems already constrained by nutrient availability since many European tree species, especially European beech and oak, are experiencing significant deterioration of their foliar nutrition (Jonard et al. 2015). Therefore, among our study sites, the plots in Chimay, Virton and Louvain-la-Neuve that already display some latent deficiency regarding P concentrations (Titeux et al. 2018) should behave closer to the simulations with constant CO2 than the Baileux plots which present better foliar nutrition. On the other hand, one could also consider that a decreased soil water content would reduce nutrient availability and consequently NPP but this aspect is not taken into account in this study. However, as the simulated increase in transpiration deficit is limited, one can consider that the impact on nutrient availability would remain very low.

The dominant effect of the CO2 fertilization in the long-term trend of forest productivity is consistent with similar studies. In his global review mainly focused on temperate and boreal forests in Europe and North America, Reyer (2015) showed that, for the simulation studies in which the atmospheric CO2 levels were held constant, the simulated NPP change with regards to reference conditions varied between − 20 and + 33% for a median value of + 5%. When both climate change and atmospheric CO2 rise were taken into account, most of the simulated biomass production increased relative to the historical period. In our simulations, when the atmospheric CO2 was maintained constant, the NPP increase ranged from 0.1 to 5.0% (Fig. 1a), which is close to the value pointed out in the abovementioned review. When the atmospheric CO2 concentration changed, the NPP increase was between 7.8 and 34.2% (Fig. 1b), which is again in good agreement with Reyer (2015) that displayed a median value 20% higher than that of the historical period.

Fernandez-Martinez et al. (2014b) analysed an extensive global dataset and showed that NPP was mainly determined by water availability, warm period length and nitrogen deposition. These results, which come from temporally averaged measurements and therefore reflect spatial rather than temporal patterns, can be considered as a spatial corollary of our results that highlight the significance of the same variables or related ones.

We have mainly discussed here the average impact of climate change. However, many papers highlight the importance of extreme heat and drought waves on tree growth and mortality (Fuhrer et al. 2006; Lindner et al. 2010; Allen et al. 2010; Teskey et al. 2015). These extreme events are important because tree functioning is a complex set of non-linear mechanisms where threshold exceedance can generate feedbacks and totally deregulate their functioning (Thompson 2011; Reyer et al. 2015; D’Orangeville et al. 2018). In HETEROFOR, the leaf-level processes (photosynthesis, respiration and transpiration) are climate dependent and take the impact of heat or drought waves into account, especially on tree growth. In addition, as most processes in HETEROFOR are calculated at the hourly time scale, temperature peaks are not smoothed as in models working at the daily or monthly time scale. The impact of these extreme climate events is, however, only partly accounted for since mortality by hydraulic failure and leaf shedding is not considered. Furthermore, tree mortality driven by pests and diseases, which is often promoted by a succession of drought and heat waves (Allen et al. 2010) is not yet included in HETEROFOR. This could partly explain why we did not observe any changes in NPP variability in our simulations between the different scenarios and time periods. Moreover, as these elements are mostly detrimental, our simulation results should be seen as the upper estimates for NPP.

4.2 What can be learnt from the decomposition of the NPP variability?

The main originality of our study is the decomposition of the NPP variability, which allows estimating the relative importance of the temporal trend compared with the inter-annual variations and evaluating the extent to which site components (climate, soil and stand) could modulate the impact of climate change on NPP. The part of the NPP variability explained by the stand effect gives an idea of the leeway left to forest managers to adapt the forests to changing conditions. Furthermore, this is of primary concern for forest managers because NPP can be considered as closely related to the timber volume under the hypothesis of allometry conservation and as a good proxy for forest health and the provision of most of the ecosystem services (e.g. Costanza et al., 1998; Dobbertin, 2005; Costanza et al., 2007; Van Oudenhoven et al., 2012; Vargas et al., 2019).

When the atmospheric CO2 was kept constant, no overall temporal trend was observed but considering a trend randomly changing with the site explained 34% of the variability (Eq. (3) and Table 4 (a)). This random effect was ascribed to differences in soil water properties and in tree species sensitivity to drought and phenology since it almost totally disappeared when the transpiration deficit and the vegetation period were included in the linear mixed model (Eq. (4) and Table 5 (a)). These factors accounted also for some of the site effect that decreased from 39.2 to 36.5%. The remaining “unexplained” variability (27%) was mainly due to inter-annual climate variations since it totally vanished when the transpiration deficit and the vegetation period were added in the model (Tables 4 (a) and 5 (a)). Among these two drivers, the transpiration deficit had a much greater explanatory power than the vegetation period (partial R2 of 0.582 vs 0.025).

Using a time-dependent atmospheric CO2 concentration generated a strong temporal trend. This temporal trend accounted for 22% of the NPP variability while its modulation by the site explained another 24% (Eq. (3) and Table 4 (b)). The integration of the atmospheric CO2 concentration, the transpiration deficit and the vegetation period in the linear mixed model made disappear the entire temporal trend and most of its variation among sites (3% remaining after the inclusion). The part of the variability explained by the atmospheric CO2 concentration corresponded exactly to that associated with the temporal trend (22%) (Eq. (4) and Table 5 (b)). The site-dependent component of the trend was mainly ascribed to differences in transpiration deficit among sites with a minor role also played by the vegetation period. The rest of the NPP variability associated with the site amounted to 30%. As for the constant atmospheric CO2 modality, we considered that the remaining “unexplained” variability (23%) was mainly due to inter-annual climate variations since it disappeared when the transpiration deficit and the vegetation period were added in the linear mixed model (Tables 4 (b) and 5 (b)). In this case, the transpiration deficit had also a much higher explanatory power than vegetation period (43% vs 2%).

Interchanging the climate, soil and stand files for the historical period under constant CO2 concentration allowed us to get a deeper understanding of the site effect. According to Table 6 (a) (Eq. 5), when they were the only variables included, the stand and soil components had a similar contribution in explaining the site effect (32 and 26%, respectively) while the part explained by differences in climate among sites was very low (1.6%). This was however not surprising as we examined a broad range of soils and stands but a much narrower range of climates. Anyway, integrating the transpiration deficit and the vegetation period in the linear mixed model explained most of the soil effect (which dropped to 0.6%) but only a little part of the stand effect, which remained at 30% (Eq. (6) and Table 6 (b)). These 30% represent the amount of freedom the forest managers have to influence the forest productivity under the current climate conditions. For the future, one must also consider the interaction between the site effect and the climate change (including atmospheric CO2) whose relative importance in explaining the NPP variability is of the same order of magnitude than the site effect (Table 4). Foresters can also act on the transpiration deficit through tree species selection and thinning even if transpiration deficit is also strongly determined by climate conditions and soil water properties.

In all stands, thinning significantly decreased the transpiration deficit and this decrease was much more pronounced in beech-dominated stands. Still, the transpiration deficit levels before and after thinning were lower in oak-dominated than in beech-dominated stands (Fig. 3). This positive short-term effect of thinning on a stand response to drought was pointed out by various studies. For example, a Douglas-fir stand showed a decrease in evapotranspiration of 30 mm (17%) during the first year after thinning. This effect was progressively reduced during the next 4 years before evapotranspiration returned to its original level (Aussenac and Granier 1988). For broadleaved species, a meta-analysis highlighted the potential of thinning to mitigate growth reduction during drought events by increasing soil water availability (Sohn et al. 2016). As a result, thinning seems to be an interesting practice to reduce the projected increase in transpiration deficit and its detrimental effect on tree growth, especially for drought-sensitive tree species. However, one must keep in mind that thinning has a transitory effect and that its impact on drought resistance progressively decreases (Guillemot et al. 2015; Sohn et al. 2016). Thinning abruptly modifies stand characteristics and forest functioning due to tree removal. Then, the remaining trees react by expanding their crown and increasing their growth rate benefiting from the higher availability of resources per tree. With time, the openings in the canopy close and the effect of thinning decreases. In this study, we simulate the first year after the thinning when its effect is maximal and when the tree dimensions (especially the crown extension) are still those that characterize a stand with a higher density. Finally, in addition to thinning, some efforts to promote oak regeneration could also be recommended to increase the stand resistance to drought as this tree species is more drought tolerant.

5 Conclusion

Understanding how NPP is going to be affected in the future due to environmental changes is crucial in order to create consistent climate change adaptation strategies and preserve the forest ecosystem services. This paper aimed at assessing, for six Belgian stands, the temporal change in NPP and in two of its main drivers: transpiration deficit and vegetation period length. Concomitantly, the influence of the CO2 fertilization effect and the impact of thinning operations were evaluated. We did not detect any trend under the three contrasted GHG emission scenarios (RCP2.6, 4.5 and 8.5) when atmospheric CO2 concentration was held constant but NPP showed a significant increase ranging from 9.4 to 34.2% for the time-dependent atmospheric CO2 concentration. Behind the apparent lack of temporal trend in NPP for simulations with constant atmospheric CO2 lies a compensatory effect of the transpiration deficit that slightly increased with time and had a pronounced negative effect on NPP and the vegetation period that became substantially longer but with a less marked impact on NPP. The site effect modulated these temporal trends and accounted for a substantial part of the NPP variability, which is encouraging for forest managers who have still a lot of possibilities to adapt their forest to changing conditions. Among others, thinning appeared very effective to decrease transpiration deficit, especially in beech-dominated stands. Forest practitioners could regularly decrease stand density or promote oak regeneration to limit the negative effect of the transpiration deficit. In the future, we plan to extend our methodology at the European scale to expand the validity of our results.

Data availability

The datasets generated during and/or analysed during the current study are available from the Zenodo repository https://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.3744345 and citable (de Wergifosse et al. 2020b).

References

  • Adams HD, Zeppel MJ, Anderegg WR, et al. (59 more authors) (2017) A multi-species synthesis of physiological mechanisms in drought-induced tree mortality. Nat. Ecol. Evol. 1(9): 1285, 1291.

  • Ainsworth EA, Long SP (2005) What have we learned from 15 years of free-air CO2 enrichment (FACE)? A meta-analytic review of the responses of photosynthesis, canopy properties and plant production to rising CO2. New Phytol. 165(2):351–372

    PubMed  Google Scholar 

  • Ainsworth EA, Rogers A (2007) The response of photosynthesis and stomatal conductance to rising [CO2]: mechanisms and environmental interactions. Plant, Cell Environ. 30(3):258–270

    CAS  Google Scholar 

  • Allen RG, Pereira LS, Raes D, Smith M (1998) Crop evapotranspiration-guidelines for computing crop water requirements-FAO irrigation and drainage paper 56. FAO, Rome.

  • Allen CD, Macalady AK, Chenchouni H, Bachelet D, McDowell N, Vennetier M, Kitzberger T, Rigling A, Breshears DD, Hogg EH, Gonzalez P, Fensham R, Zhang Z, Castro J, Demidova N, Lim JH, Allard G, Running SW, Semerci A, Cobb N (2010) A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. Forest Ecol. Manag. 259(4):660–684

    Google Scholar 

  • Anderegg WR, Konings AG, Trugman AT, Yu K, Bowling DR, Gabbitas R, Zenes N (2018) Hydraulic diversity of forests regulates ecosystem resilience during drought. Nature 561(7724):538–541

    CAS  PubMed  Google Scholar 

  • André F, Jonard M, Ponette Q (2008a) Influence of species and rain event characteristics on stemflow volume in a temperate mixed oak–beech stand. Hydrol. Process. 22(22):4455–4466

    Google Scholar 

  • André F, Jonard M, Ponette Q (2008b) Precipitation water storage capacity in a temperate mixed oak–beech canopy. Hydrol. Process. 22(20):4130–4141

    Google Scholar 

  • André F, Jonard M, Ponette Q (2010) Biomass and nutrient content of sessile oak (Quercus petraea (Matt.) Liebl.) and beech (Fagus sylvatica L.) stem and branches in a mixed stand in southern Belgium. Sci.Total Environ 408(11):2285–2294

    PubMed  Google Scholar 

  • André F, Jonard M, Jonard F, Ponette Q (2011) Spatial and temporal patterns of throughfall volume in a deciduous mixed-species stand. J. Hydrol. 400(1-2):244–254

    Google Scholar 

  • Aussenac G, Granier A (1988) Effects of thinning on water stress and growth in Douglas-fir. Can. J. For. Res. 18(1):100–105

    Google Scholar 

  • Baldocchi D (2008) ‘Breathing’ of the terrestrial biosphere: lessons learned from a global network of carbon dioxide flux measurement systems. Aust. J. Bot. 56(1):1–26

    CAS  Google Scholar 

  • Bastrup-Birk A, Reker J, Zal N (2016) European Forest Ecosystems: State and Trends. EEA Report, Copenhagen, Denmark

    Google Scholar 

  • Bauhus J, Puettmann KJ, Kühne C (2013) Close-to-nature forest management in Europe: does it support complexity and adaptability of forest ecosystems. In: Messier C, Puettmann KJ, Coates KD (eds) Managing forests as complex adaptive systems: building resilience to the challenge of global change, Oxford, pp. 187-213.

  • Bello J, Hasselquist NJ, Vallet P, Kahmen A, Perot T, Korboulewsky N (2019a) Complementary water uptake depth of Quercus petraea and Pinus sylvestris in mixed stands during an extreme drought. Plant Soil 437(1-2):93–115

    CAS  Google Scholar 

  • Bello J, Vallet P, Perot T, Balandier P, Seigner V, Perret S, Couteau C, Korboulewsky N (2019b) How do mixing tree species and stand density affect seasonal radial growth during drought events? Forest Ecol. Manag. 432:436–445

    Google Scholar 

  • Bennett AC, McDowell NG, Allen CD, Anderson-Teixeira KJ (2015) Larger trees suffer most during drought in forests worldwide. Nat. Plants 1(10):15139

    PubMed  Google Scholar 

  • Boisvenue C, Running SW (2006) Impacts of climate change on natural forest productivity–evidence since the middle of the 20th century. Global Change Biol. 12(5):862–882

    Google Scholar 

  • Bonan GB (2008) Forests and climate change: forcings, feedbacks, and the climate benefits of forests. Science 320(5882):1444–1449

    CAS  PubMed  Google Scholar 

  • Bréda N, Huc R, Granier A, Dreyer E (2006) Temperate forest trees and stands under severe drought: a review of ecophysiological responses adaptation processes and long-term consequences. Ann. For. Sci. 63(6):625–644

    Google Scholar 

  • Carl C, Biber P, Veste M, Landgraf D, Pretzsch H (2018) Key drivers of competition and growth partitioning among Robinia pseudoacacia L. trees. Forest Ecol. Manag 430:86–93

    Google Scholar 

  • Charru M, Seynave I, Hervé JC, Bertrand R, Bontemps JD (2017) Recent growth changes in Western European forests are driven by climate warming and structured across tree species climatic habitats. Ann. For. Sci. 74(2):33

    Google Scholar 

  • Chuine I (2000) A united model for budburst of trees. J Theor Biol 2007:337–347

    Google Scholar 

  • Chuine I, Bonhomme M, Legave JM, García de CortázarAtauri I, Charrier G, Lacointe A, Améglio T (2016) Can phenological models predict tree phenology accurately in the future? The unrevealed hurdle of endodormancy break. Glob. Change Biol 22:3444–3460

    Google Scholar 

  • Ciais P, Reichstein M, Viovy N, Granier A, Ogée J, Allard V, Aubinet M, Buchmann N, Bernhofer C, Carrara A, Chevallier F, De Noblet N, Friend AD, Friedlingstein P, Grünwald T, Heinesch B, Keronen P, Knohl A, Krinner G, Loustau D, Manca G, Matteucci G, Miglietta F, Ourcival JM, Papale D, Pilegaard K, Rambal S, Seufert G, Soussana JF, Sanz MJ, Schulze ED, Vesala T, Valentini R (2005) Europe-wide reduction in primary productivity caused by the heat and drought in 2003. Nature 437(7058):529–533

    CAS  PubMed  Google Scholar 

  • Cole EF, Sheldon BC (2017) The shifting phenological landscape: within-and between-species variation in leaf emergence in a mixed-deciduous woodland. Ecol. Evol. 7(4):1135–1147

    PubMed  PubMed Central  Google Scholar 

  • Costanza, R, d’Arge R, De Groot R, Farber S, Grasso M, Hannon B, Limburg K, Naeem S, O'Neill RV, Paruelo J, Raskin RG, Sutton P, van den Belt M (1998) The value of ecosystem services: putting the issues in perspective. Ecol. Econ. 25(1): 67–72

  • Costanza R, Fisher B, Mulder K, Liu S, Christopher T (2007) Biodiversity and ecosystem services: A multi-scale empirical study of the relationship between species richness and net primary production. Ecol. Econ. 61(2-3): 478-491

  • Courbaud B, de Coligny F, Cordonnier T (2003) Simulating radiation distribution in a heterogeneous Norway spruce forest on a slope. Agr. Forest Meteorol. 116:1–18

    Google Scholar 

  • Couvreur V, Vanderborght J, Javaux M (2012) A simple threedimensional macroscopic root water uptake model based on the hydraulic architecture approach. Hydrol. Earth Syst. Sci. 16:2957–2971

    Google Scholar 

  • D’Orangeville L, Houle D, Duchesne L, Phillips RP, Bergeron Y, Kneeshaw D (2018) Beneficial effects of climate warming on boreal tree growth may be transitory. Nat. Commun. 9(1):3213

    PubMed  PubMed Central  Google Scholar 

  • de Wergifosse L, André F, Beudez N, de Coligny F, Goosse H, Jonard F, Ponette Q, Titeux H, Vincke C, Jonard M (2020a) HETEROFOR 1.0: a spatially explicit model for exploring the response of structurally complex forests to uncertain future conditions. II. Phenology and water cycle. Geosci. Model Dev. 13:1459–1498

    Google Scholar 

  • de Wergifosse L, André F, Goosse H, Caluwaerts S, de Cruz L, de Troch R, van Schaeybroeck B and Jonard M (2020b) CO2 fertilization, transpiration deficit and vegetation period drive the response of mixed broadleaved forests to a changing climate in Wallonia. V1. Zenodo. [Dataset]. https://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.3744345

  • DeRose RJ, Long JN (2014) Resistance and resilience: a conceptual framework for silviculture. For. Sci. 60(6):1205–1212

    Google Scholar 

  • Dragoni D, Schmid HP, Wayson CA, Potter H, Grimmond CSB, Randolph JC (2011) Evidence of increased net ecosystem productivity associated with a longer vegetated season in a deciduous forest in south-central Indiana, USA. Global Change Biol. 17(2):886–897

    Google Scholar 

  • Dobbertin M (2005) Tree growth as indicator of tree vitality and of tree reaction to environmental stress: a review. Eur. J. For. Res. 124(4): 319–333

  • Dufrêne E, Davi H, François C, Le Maire G, Le Dantec V, Granier A (2005) Modelling carbon and water cycles in a beech forest: Part I: Model description and uncertainty analysis on modelled NEE. Ecol. Model. 185: 407–436

  • Dufour-Kowalski S, Courbaud B, Dreyfus P, Meredieu C, De Coligny F (2012) Capsis: an open software framework and community for forest growth modelling. Ann. For. Sci. 69:221–233

    Google Scholar 

  • Farquhar GD, von Caemmerer S, Berry JA (1980) A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 149:78–80

    CAS  PubMed  Google Scholar 

  • Fernandez-Martinez M, Vicca S, Janssens IA et al (2014a) Nutrient availability as the key regulator of global forest carbon balance. Nature Climate Change 4:471–476

  • Fernández-Martínez M, Vicca S, Janssens IA, Luyssaert S, Campioli M, Sardans J, Estiarte M, Peñuelas J (2014b) Spatial variability and controls over biomass stocks, carbon fluxes, and resource-use efficiencies across forest ecosystems. Trees 28(2):597–611

  • Ferretti M, Fischer R (2013) Forest monitoring: methods for terrestrial investigations in Europe with an overview of North America and Asia. Elsevier, Oxford, United Kingdom

    Google Scholar 

  • Flynn DFB, Wolkovich EM (2018) Temperature and photoperiod drive spring phenology across all species in a temperate forest community. New Phytol. 219(4):1353–1362

    CAS  PubMed  Google Scholar 

  • Fu YH, Piao S, Zhao H, Jeong SJ, Wang X, Vitasse Y, Ciais P, Janssens IA (2014) Unexpected role of winter precipitation in determining heat requirement for spring vegetation green-up at northern middle and high latitudes. Glob. Change Biol. 20(12):3743–3755

    Google Scholar 

  • Fu Z, Stoy PC, Luo Y, Chen J, Sun J, Montagnani L, Wohlfahrt G, Rahman AF, Rambal S, Bernhofer C, Wang J, Shirkey G, Niu S (2017) Climate controls over the net carbon uptake period and amplitude of net ecosystem production in temperate and boreal ecosystems. Agr. Forest Meteorol. 243:9–18

    Google Scholar 

  • Fuhrer J, Beniston M, Fischlin A, Frei C, Goyette S, Jasper K, Pfister C (2006) Climate risks and their impact on agriculture and forests in Switzerland. Clim. Change 79:79–102

    CAS  Google Scholar 

  • Giorgi F, Jones C, Asrar GR (2009) Addressing climate information needs at the regional level: the CORDEX framework. WMO Bulletin 58(3):175

    Google Scholar 

  • Giot O, Termonia P, Degrauwe D, De Troch R, Caluwaerts S, Smet G, Berckmans J, Deckmyn A, De Cruz L, De Meutter P, Duerinckx A, Gerard L, Hamdi R, Van den Bergh J, Van Ginderachter M, Van Schaeybroeck B (2016) Validation of the ALARO-0 model within the EURO-CORDEX framework. Geosci. Model Dev. 9(3):1143–1152

    Google Scholar 

  • Granier A, Reichstein M, Bréda N, Janssens IA, Falge E, Ciais P, Grünwald T, Aubinet M, Berbigier P, Bernhofer C, Buchmann N, Facini O, Grassi G, Heinesch B, Ilvesniemi H, Keronen P, Knohl A, Köstner B, Lagergren F, Lindroth A, Longdoz B, Loustau D, Mateus J, Montagnani L, Nys C, Moors E, Papale D, Peiffer M, Pilegaard K, Pita G, Pumpanen P, Rambal S, Rebmann C, Rodrigues A, Seufert G, Tenhunen J, Vesala T, Wang Q (2007) Evidence for soil water control on carbon and water dynamics in European forests during the extremely dry year: 2003. Agr. Forest Meteorol. 143(1-2):123–145

    Google Scholar 

  • Grossiord C (2019) Having the right neighbors: how tree species diversity modulates drought impacts on forests. New Phytol. doi.org/10.1111/nph.15667.

  • Guillemot J, Klein EK, Davi H, Courbet F (2015) The effects of thinning intensity and tree size on the growth response to annual climate in Cedrus atlantica: a linear mixed modeling approach. Ann. For. Sci. 72(5):651–663

    Google Scholar 

  • Holdridge LR (1967) Life zone ecology. Tropical Science Center, San Jose, Costa Rica

    Google Scholar 

  • IPCC (2013) Summary for Policymakers. In: Stocker TF, Qin D, Plattner GK, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Midgley PM (eds) Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, pp 3–29

    Google Scholar 

  • Jacob D, Petersen J, Eggert B, Alias A., Christensen O.B., Bouwer L.M., Braun A., Colette A., Déqué M., Georgievski G., Georgopoulou E., Gobiet A., Menut L., Nikulin G., Haensler A., Hempelmann N., Jones C., Keuler K., Kovats S., Kröner N., Kotlarski S., Kriegsmann A., Martin E., van Meijgaard E., Moseley C., Pfeifer S., Preuschmann S., Radermacher C., Radtke K., Rechid D., Rounsevell M., Samuelsson P., Somot S., Soussana J.F., Teichmann C., Valentini R., Vautard R., Weber B., Yiou P. (36 more authors) (2014). EURO-CORDEX: new high-resolution climate change projections for European impact research. Reg. Environ. Change 14(2): 563-578.

  • Jacob D, Kotova L, Teichmann C, Sobolowski SP, Vautard R, Donnelly C, Koutroulis AG, Grillakis MG, Tsanis IK, Damm A, Sakalli A, van Vliet M (2018) Climate impacts in Europe under +1.5°C global warming. Earths Future 6(2):264–285

    Google Scholar 

  • Jeong SJ, Ho CH, Gim HJ, Brown ME (2011) Phenology shifts at start vs. end of growing season in temperate vegetation over the Northern Hemisphere for the period 1982–2008. Global Change Biol. 17(7):2385–2399

    Google Scholar 

  • Jonard M, Andre F, Ponette Q (2006) Modeling leaf dispersal in mixed hardwood forests using a ballistic approach. Ecology 87(9):2306–2318

    PubMed  Google Scholar 

  • Jonard M, André F, Jonard F, Mouton N, Procès P, Ponette Q (2007) Soil carbon dioxide efflux in pure and mixed stands of oak and beech. Ann. For. Sci. 64(2):141–150

    CAS  Google Scholar 

  • Jonard M, André F, Ponette Q (2008) Tree species mediated effects on leaf litter dynamics in pure and mixed stands of oak and beech. Can. J. For. Res. 38(3):528–538

    Google Scholar 

  • Jonard F, André F, Ponette Q, Vincke C, Jonard M (2011) Sap flux density and stomatal conductance of European beech and common oak trees in pure and mixed stands during the summer drought of 2003. J. Hydrol. 409(1-2):371–381

    Google Scholar 

  • Jonard M, Fürst A, Verstraeten A, Thimonier A, Timmermann V, Potočić N, Waldner P, Benham S, Hansen K, Merilä P, Ponette Q, de la Cruz AC, Roskams P, Nicolas M, Croisé L, Ingerslev M, Matteucci G, Decinti B, Bascietto M, Rautio P (2015) Tree mineral nutrition is deteriorating in Europe. Global Change Biol. 21(1):418–430

    Google Scholar 

  • Jonard M, André F, de Coligny F, de Wergifosse L, Beudez N, Davi H, Ligot G, Ponette Q, Vincke C (2020) HETEROFOR 1.0: a spatially explicit model for exploring the response of structurally complex forests to uncertain future conditions. I. Carbon fluxes and tree dimensional growth. Geosci. Model Dev 13:905–935

    Google Scholar 

  • Kint V, Aertsen W, Campioli M, Vansteenkiste D, Delcloo A, Muys B (2012) Radial growth change of temperate tree species in response to altered regional climate and air quality in the period 1901–2008. Clim. Change 115(2):343–363

    CAS  Google Scholar 

  • Körner C (2006) Plant CO2 responses: an issue of definition, time and resource supply. New Phytol. 172(3):393–411

    PubMed  Google Scholar 

  • Kovats R, Valentini RS, Bouwer M, Georgopoulou E, Jacob D, Martin E, Rounsevell M, Soussana JF (2014) Chapter 23 - Europe In: Barros VR, Field CB, Dokken DJ, Mastrandrea MD, Mach KJ, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL (eds) Climate Change 2014: Impacts Adaptation and Vulnerability Part B: Regional Aspects, Cambridge University Press, Cambridge, pp 1267-1326.

  • Le Quéré C, Andrew RM, Friedlingstein P, Sitch S, Hauck J, Pongratz J, Arneth A (2018) Global carbon budget 2018. Earth Syst. Sci. Data 10:2141–2194

    Google Scholar 

  • Leuzinger S, Zotz G, Asshoff R, Körner C (2005) Responses of deciduous forest trees to severe drought in Central Europe. Tree Physiol. 25(6):641–650

    PubMed  Google Scholar 

  • Lindner M, Maroschek M, Netherer S, Kremer A, Barbati A, Garcia-Gonzalo J, Seidl R, Delzon S, Corona P, Kolström M, Lexer MJ, Marchetti M (2010) Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems. Forest Ecol. Manag. 259(4):698–709

    Google Scholar 

  • Lindner M, Fitzgerald JB, Zimmermann NE, Reyer C, Delzon S, van der Maaten E, Schelhaas MJ, Lasch P, Eggers J, van der Maaten-Theunissen M, Suckow F, Psomas A, Poulter B, Hanewinkel M (2014) Climate change and European forests: what do we know, what are the uncertainties, and what are the implications for forest management? J. Environ. Manage. 146:69–83

    PubMed  Google Scholar 

  • Lisar SY, Motafakkerazad R, Hossain MM, Rahman IM (2012) Water stress in plants: causes effects and responses. In: Rahman I, Hasegawa H (eds) Water stress. InTech, Rijeka, Croatia, pp 15–42

    Google Scholar 

  • Maraun D, Widmann M (2018) Statistical downscaling and bias correction for climate research. Cambridge University Press, Cambridge, United Kingdom

    Google Scholar 

  • Masson-Delmotte V (2018) Global Warming of 1.5°C: An IPCC Special Report on the Impacts of Global Warming of 1.5° C Above Pre-industrial Levels and Related Global Greenhouse Gas Emission Pathways, in the Context of Strengthening the Global Response to the Threat of Climate Change, Sustainable Development, and Efforts to Eradicate Poverty. Cambridge University Press, Cambridge, United Kingdom.

  • McDowell NG, Allen CD (2015) Darcy’s law predicts widespread forest mortality under climate warming. Nat. Clim. Change 5(7):669

    Google Scholar 

  • McDowell NG, Pockman WT, Allen CD, Breshears DD, Cobb N, Kolb T, Yepez EA (2008) Mechanisms of plant survival and mortality during drought: why do some plants survive while others succumb to drought? New Phytol. 178(4):719–739

    PubMed  Google Scholar 

  • Medlyn BE, Barton CVM, Broadmeadow MSJ, Ceulemans R, De Angelis P, Forstreuter M, Rey A (2001) Stomatal conductance of forest species after long-term exposure to elevated CO2 concentration: a synthesis. New Phytol. 149(2):247–264

    Google Scholar 

  • Menzel A, Sparks TH, Estrella N, Koch E, Aasa A, Ahas R, Alm-Kübler K, Bissolli P, Braslavska O, Briede A, Chmielewski FM, Crepinsek Z, Curnel Y, Dahl A, Defila C, Donnelly A, Filella Y, Jatczak K, Mage F, Mestre A, Nordli O, Peñuelas J, Pirinen P, Remisova V, Scheifinger H, Striz M, Susnik A, Van Vliet A, Wielgolaski FE, Zach S, Zust A (2006) European phenological response to climate change matches the warming pattern. Global Change Biol. 12(10):1969–1976

    Google Scholar 

  • Mette T, Dolos K, Meinardus C, Bräuning A, Reineking B, Blaschke M, Wellstein C (2013) Climatic turning point for beech and oak under climate change in Central Europe. Ecosphere 4(12):1–19

    Google Scholar 

  • Mina M, Bugmann H, Cordonnier T, Irauschek F, Klopcic M, Pardos M, Cailleret M (2017) Future ecosystem services from European mountain forests under climate change. J. Appl. Ecol. 54(2):389–401

    Google Scholar 

  • Monteith JL (1965) Evaporation and environment. Symposia of the Society for Experimental Biology 19:205–234

    CAS  PubMed  Google Scholar 

  • Morin X, Chuine I (2014) Will tree species experience increased frost damage due to climate change because of changes in leaf phenology? Can. J. For. Res. 44(12):1555–1565

    Google Scholar 

  • Myhre G, Shindell D, Breon FM, Collins W, Fuglestvedt J, Huang J, Koch D, Lamarque JF, Lee D, Mendoza B, Nakajima T, Robock A, Stephens G, Takemura T, Zhang H (2013) Anthropogenic and Natural Radiative Forcing. In: Stocker TF, Qin D, Plattner GK, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Migley PM (eds) Climate Change 2013: The Physical Science Basis. Cambridge University Press, Cambridge, United Kingdom, pp 659–740

    Google Scholar 

  • Norby RJ, DeLucia EH, Gielen B, Calfapietra C, Giardina CP, King JS, Ledford J, McCarthy HR, Moore D, Ceulemans R, De Angelis P, Finzi AC, Karnosky DF, Kubiske ME, Lukac M, Pregitzer KS, Scarascia-Mugnozza G, Schlesinger WH, Oren R (2005) Forest response to elevated CO2 is conserved across a broad range of productivity. Proc. Natl. Acad. Sci. 102(50):18052–18056

    CAS  PubMed  Google Scholar 

  • Norby RJ, Warren JM, Iversen CM, Medlyn BE, McMurtrie RE (2010) CO2 enhancement of forest productivity constrained by limited nitrogen availability. Proc. Natl. Acad. Sci. 107(45):19368–19373

    CAS  PubMed  Google Scholar 

  • Ogaya R, Peñuelas J (2007) Tree growth mortality and above-ground biomass accumulation in a holm oak forest under a five-year experimental field drought. Plant Ecol. 189(2):291–299

    Google Scholar 

  • Oren R, Ellsworth DS, Johnsen KH, Phillips N, Ewers BE, Maier C, Schäfer KVR, McCarthy H, Hendrey G, McNulty SG, Katul GG (2001) Soil fertility limits carbon sequestration by forest ecosystems in CO2-enriched atmosphere. Nature 411:469–472

    CAS  PubMed  Google Scholar 

  • Osakabe Y, Osakabe K, Shinozaki K, Tran LSP (2014) Response of plants to water stress. Front. Plant Sci. 5:86

    PubMed  PubMed Central  Google Scholar 

  • Park T, Ganguly S, Tømmervik H, Euskirchen ES, Høgda KA, Karlsen SR, Brovkin V, Nemani RR, Myneni RB (2016) Changes in growing season duration and productivity of northern vegetation inferred from long-term remote sensing data. Environ. Res. 11(8):084001

    Google Scholar 

  • Parmesan C, Yohe G (2003) A globally coherent fingerprint of climate change impacts across natural systems. Nature 421(6918):37–42

    CAS  PubMed  Google Scholar 

  • Petit S, Claessens H, Vincke C, Ponette Q, Marchal D (2017) Le Fichier écologique des essences, version 2.0. Forêt. Nature 143:12–19

    Google Scholar 

  • Phillips RP, Ibáñez I, D’Orangeville L, Hanson PJ, Ryan MG, McDowell NG (2016) A belowground perspective on the drought sensitivity of forests: towards improved understanding and simulation. Forest Ecol. Manag. 380:309–320

    Google Scholar 

  • Paschalis A, Fatichi S, Zscheischler J, Ciais P, Bahn M, Boysen L, Chang J, De Kauwe M, Estiarte M, Goll D, Hanson PJ, Harper AB, Hou E, Kigel J, Knapp AK, Larsen KS, Li W, Lienert S, Luo Y, Meir P, Nabel J, Ogaya R, Parolari AJ, Peng C, Peñuelas J, Pongratz J, Rambal S, Schmidt IK, Shi H, Sternberg M, Tian H, Tschumi E, Ukkola A, Vicca S, Viovy N, Wang YP, Wang Z, Williams K, Wu D, Zhu Q (2020) Rainfall manipulation experiments as simulated by terrestrial biosphere models: Where do we stand?. Global Change Biol. 26(6): 3336–3355

  • Piao S, Liu Q, Chen A, Janssens IA, Fu Y, Dai J, Liu L, Shen M, Zhu X (2019) Plant phenology and global climate change: current progresses and challenges. Glob. Change Biol. 25(6):1922–1940

    Google Scholar 

  • Pletsers A, Caffarra A, Kelleher CT, Donnelly A (2015) Chilling temperature and photoperiod influence the timing of bud burst in juvenile Betula pubescens Ehrh. and Populus tremulaL.trees. Ann. For. Sci. 72:941–953

    Google Scholar 

  • Pretzsch H, Biber P (2016) Tree species mixing can increase maximum stand density. Can. J. For. Res. 46(10):1179–1193

    Google Scholar 

  • Pretzsch H, Forrester DI, Rötzer T (2015) Representation of species mixing in forest growth models. A review and perspective. Ecol. Model. 313:276–292

    Google Scholar 

  • Primack RB, Ibáñez I, Higuchi H, Lee SD, Miller-Rushing AJ, Wilson AM, Silander JA Jr (2009) Spatial and interspecific variability in phenological responses to warming temperatures. Biol. Conserv. 142(11):2569–2577

    Google Scholar 

  • Raftoyannis Y, Radoglou K (2002) Physiological responses of beech and sessile oak in a natural mixed stand during a dry summer. Ann. Bot. 89(6):723–730

    PubMed  PubMed Central  Google Scholar 

  • Reyer C (2015) Forest productivity under environmental change—a review of stand-scale modelling studies. Curr. For. Rep. 1(2):53–68

    Google Scholar 

  • Reyer C, Lasch-Born P, Suckow F, Gutsch M, Murawski A, Pilz T (2014) Projections of regional changes in forest net primary productivity for different tree species in Europe driven by climate change and carbon dioxide. Ann. For. Sci. 71(2):211–225

    Google Scholar 

  • Reyer CP, Brouwers N, Rammig A, Brook BW, Epila J, Grant RF, Holmgren M, Langerswich F, Leuzinger S, Lucht W, Medlyn B, Pfeifer M, Steinkamp J, Vanderwel MC, Verbeeck H, Villela DM (2015) Forest resilience and tipping points at different spatio-temporal scales: approaches and challenges. J. Ecol. 103(1):5–15

    Google Scholar 

  • Richardson AD, Keenan TF, Migliavacca M, Ryu Y, Sonnentag O, Toomey M (2013) Climate change, phenology, and phenological control of vegetation feedbacks to the climate system. Agr. Forest Meteorol. 169:156–173

    Google Scholar 

  • RStudio Team (2015) RStudio: Integrated Development for R. RStudio. Boston, MA (Computer Software v0.98.1074).

  • Rubio-Cuadrado Á, Camarero JJ, Del Rio M, Sánchez-González M, Ruiz-Peinado R, Bravo-Oviedo A, Montes F (2018) Drought modifies tree competitiveness in an oak-beech temperate forest. Forest Ecol. Manag. 429:7–17

    Google Scholar 

  • Rummukainen M (2010) State-of-the-art with regional climate models. WIRES Clim. Change 1(1):82–96

    Google Scholar 

  • Sanz-Perez V, Castro-Diez P (2010) Summer water stress and shade alter bud size and budburst date in three mediterranean Quercus species. Trees 24:89–97

    Google Scholar 

  • Scherrer D, Bader MK, Körner C (2011) Drought-sensitivity ranking of deciduous tree species based on thermal imaging of forest canopies. Agr. Forest Meteorol. 151(12):1632–1640

    Google Scholar 

  • Seidl R, Lexer MJ, Jäger D, Hönninger K (2005) Evaluating the accuracy and generality of a hybrid patch model. Tree Physiol. 25(7):939–951

    PubMed  Google Scholar 

  • Seidl R, Schelhaas MJ, Rammer W, Verkerk PJ (2014) Increasing forest disturbances in Europe and their impact on carbon storage. Nat. Clim. Change 4(9):806

    CAS  Google Scholar 

  • Sohn JA, Saha S, Bauhus J (2016) Potential of forest thinning to mitigate drought stress: a meta-analysis. Forest Ecol. Manag. 380:261–273

    Google Scholar 

  • Steenberg JW, Millward AA, Duinker PN, Nowak DJ, Robinson PJ (2015) Neighbourhood-scale urban forest ecosystem classification. J. Environ. Manage. 163:134–145

    PubMed  Google Scholar 

  • Stocker TF, Qin D, Plattner GK, Tignor M, Allen SK, Boschung J, Midgley PM (2013) Climate change 2013: The physical science basis. Cambridge Univ Press, Cambridge, United Kingdom

    Google Scholar 

  • Temperli C, Bugmann H, Elkin C (2012) Adaptive management for competing forest goods and services under climate change. Ecol. Appl. 22(8):2065–2077

    PubMed  Google Scholar 

  • Termonia P, Van Schaeybroeck B, De Cruz L, De Troch R, Caluwaerts S, Giot O, Hamdi R, Vannitsem S, Duchêne F, Willems P, Tabari H, Van Uytven E, Hosseinzadehtalaei P, Van Lipzig N, Wouters H, Vanden Broucke S, van Ypersele JP, Marbaix P, Villanueva-Birriel C, Fettweis X, Wyard C, Scholzen C, Doutreloup S, De Ridder K, Gobbin A, Lauwaet D, Stavrakou T, Bauwens M, Müller JF, Luyten P, Ponsar S, Van den Eynde D, Pottiaux E (2018) The CORDEX. be initiative as a foundation for climate services in Belgium. Clim. Serv. 11:49–61

    Google Scholar 

  • Teskey R, Wertin T, Bauweraerts I, Ameye M, McGuire MA, Steppe K (2015) Responses of tree species to heat waves and extreme heat events. Plant Cell Environ. 38(9):1699–1712

    PubMed  Google Scholar 

  • Thom D, Rammer W, Seidl R (2017) The impact of future forest dynamics on climate: interactive effects of changing vegetation and disturbance regimes. Ecol. Monogr. 87(4):665–684

    PubMed  PubMed Central  Google Scholar 

  • Thompson I (2011) Biodiversity, ecosystem thresholds, resilience and forest degradation. Unasylva 238(62):25–30

    Google Scholar 

  • Thompson M, Gamage D, Hirotsu N, Martin A, Seneweera S (2017) Effects of elevated carbon dioxide on photosynthesis and carbon partitioning: a perspective on root sugar sensing and hormonal crosstalk. Front. Physiol. 8:578

    PubMed  PubMed Central  Google Scholar 

  • Titeux H, Jonard M, Elodie B, Christian L, Patrick M, Ponette Q (2018) L’évolution de la santé des forêts wallonnes. Principaux enseignements des réseaux de monitoring ICP-forests. Forêt.Nature 147:48–60

    Google Scholar 

  • Tucker C, Randolph J, Evans T, Andersson K, Persha L, Green G (2008) An approach to assess relative degradation in dissimilar forests: toward a comparative assessment of institutional outcomes. Ecol. Soc. 13(1):4

    Google Scholar 

  • Van der Perre R, Bythell S, Bogaert P, Claessens H, Ridremont F, Tricot C, Vincke C, Ponette Q (2015) La carte bioclimatique de Wallonie: un nouveau découpage écologique du territoire pour le choix des essences forestières. Forêt-Nature 135:47–58

    Google Scholar 

  • Van Oudenhoven AP, Petz K, Alkemade R, Hein L, de Groot RS (2012) Framework for systematic indicator selection to assess effects of land management on ecosystem services. Ecol. Indic. (21): 110–122

  • Vitasse Y, Basler D (2013) What role for photoperiod in the bud burst phenology of European beech. Eur. J. Forest Res 132:1–8

    Google Scholar 

  • Vargas L, Willemen L, Hein L (2019) Assessing the Capacity of Ecosystems to Supply Ecosystem Services Using Remote Sensing and An Ecosystem Accounting Approach. Environ. Manage. 63(1): 1–15

  • Vitasse Y, Delzon S, Dufrêne E, Pontailler JY, Louvet JM, Kremer A, Michalet R (2009) Leaf phenology sensitivity to temperature in European trees: do within-species populations exhibit similar responses? Agr. Forest Meteorol. 149(5):735–744

    Google Scholar 

  • Warren JM, Jensen AM, Medlyn BE, Norby RJ, Tissue DT (2015) Carbon dioxide stimulation of photosynthesis in Liquidambar styraciflua is not sustained during a 12-year field experiment. AoB Plants 7.

  • Wilcoxon F (1945) Individual comparisons by ranking methods. Biometrics bulletin 1(6):80–83

    Google Scholar 

  • Wu Z, Dijkstra P, Koch GW, Peñuelas J, Hungate BA (2011) Responses of terrestrial ecosystems to temperature and precipitation change: A meta‐analysis of experimental manipulation. Global Change Biol. 17(2): 927–942

  • Xie Y, Wang X, Wilson AM, Silander JA Jr (2018) Predicting autumn phenology: how deciduous tree species respond to weather stressors. Agr. Forest Meteorol. 250:127–137

    Google Scholar 

  • Zhang Z, Zhang R, Cescatti A, Wohlfahrt G, Buchmann N, Zhu J, Chen G, Moyano F, Pumpanen J, Hirano T, Takagi K, Merbold L (2017) Effect of climate warming on the annual terrestrial net ecosystem CO2 exchange globally in the boreal and temperate regions. Sci. Rep. 7(1):3108

    CAS  PubMed  PubMed Central  Google Scholar 

  • Zimmermann NE, Normand S, Pearman PB, Psomas A (2013) Chapter III: Future ranges in European tree species. In: Fitzergald J, Lindner M (eds) Adapting to climate change in European forests–Results of the MOTIVE project. Pensoft Publishers, Sofia, Bulgaria, pp 15–21

    Google Scholar 

Download references

Acknowledgements

We thank Rafiq Hamdi and Piet Termonia from Royal Meteorological Institute of Belgium for enabling the generation of the climate data.

Funding

This study was funded by FRIA (grant no. 1.E005.18).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Louis de Wergifosse.

Ethics declarations

Conflict of interest

The authors declare that they have no conflict of interest.

Additional information

Handling Editor: Céline Meredieu

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributions of the co-authors

Conceptualization: Louis DE WERGIFOSSE, Mathieu JONARD; Methodology: Louis DE WERGIFOSSE, Mathieu JONARD; Resources: Steven CALUWAERTS, Lesley DE CRUZ, Rozemien DE TROCH, Bert VAN SCHAEYBROECK; Writing – original draft: Louis DE WERGIFOSSE; Writing – review & editing: Louis DE WERGIFOSSE, Frédéric ANDRE, Hugues GOOSSE, Bert VAN SCHAEYBROECK, Mathieu JONARD; Supervision: Hugues GOOSSE, Mathieu JONARD

This article is part of the topical collection on Mensuration and modelling for forestry in a changing environment

Annexes

Annexes

Fig. 4
figure 4

NPP comparisons among RCP scenarios and with the historical period (1976–2005) for three 30-year periods in the different study sites (a-f) and with constant CO2 atmospheric concentrations. The horizontal line corresponds to the median, the box ends indicate the upper and lower quartiles and the whiskers show the values above and below these quartiles within 1.5 interquartile. Common letters indicate that the distributions are non-significantly different according to a paired Wilcoxon signed-rank (between scenarios of the same period) or an unpaired Mann-Whitney (between scenarios of different periods) tests

Fig. 5
figure 5

NPP comparisons among RCP scenarios and with the historical period (1976–2005) for three 30-year periods in the different study sites (a-f) and with time-dependent CO2 atmospheric concentrations. The horizontal line corresponds to the median, the box ends indicate the upper and lower quartiles and the whiskers show the values above and below these quartiles within 1.5 interquartile. Common letters indicate that the distributions are non-significantly different according to a paired Wilcoxon signed-rank (between scenarios of the same period) or an unpaired Mann-Whitney (between scenarios of different periods) tests

Fig. 6
figure 6

Actual transpiration (left) and potential transpiration (right) comparisons among RCP scenarios and with the historical period (1976–2005) for three 30-year periods in the different study sites and with constant CO2 atmospheric concentrations. The horizontal line corresponds to the median, the box ends indicate the upper and lower quartiles and the whiskers show the values above and below these quartiles within 1.5 interquartile. Common letters indicate that the distributions are non-significantly different according to a paired Wilcoxon signed-rank (between scenarios of the same period) or an unpaired Mann-Whitney (between scenarios of different periods) tests

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

de Wergifosse, L., André, F., Goosse, H. et al. CO2 fertilization, transpiration deficit and vegetation period drive the response of mixed broadleaved forests to a changing climate in Wallonia. Annals of Forest Science 77, 70 (2020). https://0-doi-org.brum.beds.ac.uk/10.1007/s13595-020-00966-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1007/s13595-020-00966-w

Keywords