Next Article in Journal
Comprehensive Air Quality Assessment of the Tobacco Heating System 2.2 under Simulated Indoor Environments
Next Article in Special Issue
A Comparison of the Characteristics of Drought during the Late 20th and Early 21st Centuries over Eastern Europe, Western Russia and Central North America
Previous Article in Journal
Assessing Climate Change Trends and Their Relationships with Alpine Vegetation and Surface Water Dynamics in the Everest Region, Nepal
Previous Article in Special Issue
Spatial and Temporal Analysis of Dry and Wet Spells in the Wadi Cheliff Basin, Algeria
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing the Contrasting Effects of the Exceptional 2015 Drought on the Carbon Dynamics in Two Norway Spruce Forest Ecosystems

1
Department of Matter and Energy and Fluxes, Global Change Research Institute CAS, Bĕlidla 986/4a, 603 00 Brno, Czech Republic
2
Department of Agrosystems and Bioclimatology, Faculty of Agronomy, Mendel University in Brno, Zĕmĕdelská 1, 613 00 Brno, Czech Republic
3
Earth Observation Research and Innovation Centre, University of Energy and Natural Resources, Sunyani P.O. Box 214, Ghana
4
Department of Spatial Planning, Institute of Management, Slovak University of Technology in Bratislava, 811 07 Bratislava, Slovakia
*
Author to whom correspondence should be addressed.
Atmosphere 2021, 12(8), 988; https://doi.org/10.3390/atmos12080988
Submission received: 21 June 2021 / Revised: 23 July 2021 / Accepted: 24 July 2021 / Published: 31 July 2021

Abstract

:
The occurrence of extreme drought poses a severe threat to forest ecosystems and reduces their capability to sequester carbon dioxide. This study analysed the impacts of a central European summer drought in 2015 on gross primary productivity (GPP) at two Norway spruce forest sites representing two contrasting climatic conditions—cold and humid climate at Bílý Kříž (CZ-BK1) vs. moderately warm and dry climate at Rájec (CZ-RAJ). The comparative analyses of GPP was based on a three-year eddy covariance dataset, where 2014 and 2016 represented years with normal conditions, while 2015 was characterized by dry conditions. A significant decline in the forest GPP was found during the dry year of 2015, reaching 14% and 6% at CZ-BK1 and CZ-RAJ, respectively. The reduction in GPP coincided with high ecosystem respiration (Reco) during the dry year period, especially during July and August, when several heat waves hit the region. Additional analyses of GPP decline during the dry year period suggested that a vapour pressure deficit played a more important role than the soil volumetric water content at both investigated sites, highlighting the often neglected importance of considering the species hydraulic strategy (isohydric vs. anisohydric) in drought impact assessments. The study indicates the high vulnerability of the Norway spruce forest to drought stress, especially at sites with precipitation equal or smaller than the atmospheric evaporative demand. Since central Europe is currently experiencing large-scale dieback of Norway spruce forests in lowlands and uplands (such as for CZ-RAJ conditions), the findings of this study may help to quantitatively assess the fate of these widespread cultures under future climate projections, and may help to delimitate the areas of their sustainable production.

1. Introduction

Terrestrial ecosystems, such as forests, play a significant role in the global carbon cycle by sequestering carbon dioxide (CO2) through photosynthesis and the conversion of biomass into stable soil organic compounds [1,2,3,4,5,6]. However, carbon uptake and carbon allocation can be reduced by high vapour pressure deficits (VPDs) and soil water deficits, commonly observed during periods with increased temperature and lack of precipitation [7,8]. These physiological responses to drought in plants vary at both local and global scales, depending on the type of plant species (based on the different levels of resilience to water stress), the local climatic condition, and other additional factors [9,10,11,12,13]. It has also been widely demonstrated that climate change increases the likelihood and severity of such drought events, especially within the temperate regions [14,15,16,17,18]. Thus, despite the resilience of some tree species, the projected increase in the occurrence and severity of extreme drought events will further reduce forest ecosystem photosynthesis in many parts of Europe, especially at unsuitable locations with sub-optimal climatic conditions [19].
Drought does not have a single definition and can be used in different contexts. However, typically, it is used to describe the periods with high VPD (meteorological drought) and low soil moisture (edaphic drought) [20]. The severity of the drought impact on carbon exchange depends on the site characteristics, duration and intensity of the drought periods [18,21,22]. In Europe, the 2003 summer drought period resulted in an estimated loss of about 30% in GPP (0.5 P gC y−1) over both the northern Mediterranean forests and temperate deciduous beech forests [23,24]. This was the result of the prolonged abnormal reduction in the soil moisture and high air temperatures below the wilting point and high air temperatures recorded across Europe during the summer period in 2003. The summer drought of 2015 was considered one of the most severe drought events in Europe after the 2003 summer drought [25]. During this period, much of the European continent (Poland, the Czech Republic, Slovakia, Western Ukraine, and Belarus) was severely affected in June and July 2015, due to the exceptionally high maximum daily air temperatures and the prolonged rainfall shortage by 31% since April [26,27,28,29]. Consequently, this had a significant effect on forest growth, as observed in both spruce and beech forests, within the Czech Republic. Nonetheless, spruce forest showed a higher sensitivity to drought as compared to beech forests [30,31].
The Norway spruce (Picea abies (L.) H. Karst) has significant economic and ecological importance within Europe [32]. However, its shallow root system makes it vulnerable to drought stress, especially at lower altitudes [33]. Since spruce thrives well in cold and humid regions, a significant decrease in soil moisture coupled with high air temperatures could have severe consequences on the spruce forest growth [34]. Therefore, the 2015 summer drought episode provided a good opportunity to further examine the response of spruce forest stands with different climatic conditions to extreme drought stress events.
The eddy covariance (EC) technique provides direct measurements of CO2 exchange between the atmosphere and the underlying ecosystem [35]. It is a convenient and widely used approach for observation of a forest stand carbon uptake and its dynamic response to environmental variables. Obtained CO2 fluxes are integrated over hundreds of square meters and resolved to half-hourly intervals [35,36]. Continuous micrometeorological and EC measurements over multiple years allow to detect drought periods and evaluate their impacts on CO2 exchange.
In this study, we hypothesized that the impact of extreme drought conditions will be more severe on the Norway spruce forest stand located in dry and hot climates, due to higher VPD and the sub-optimal supply of soil moisture in such drier climates. To test this hypothesis, we sought to evaluate the different effects of the 2015 drought during the main growing season in the wet and dry spruce forest ecosystems, compared to normal climatic conditions within the same period of two adjacent years (2014 and 2016). Secondly, we aimed to determine the influence of critical site-specific environmental factors (such as the VPD and soil volumetric water content) on the drought stress response at each of the spruce forest stands located in mountainous (around 900 m a.s.l) and highland (up to 600 m a.s.l) regions within the Czech Republic, using long-term eddy covariance CO2 flux data.

2. Materials and Methods

2.1. Station Description

The study used multiyear (2014–2016) measurements from Bílý Kříž and Rájec ecosystem stations that are part of the CzeCOS (Czech Carbon Observation System; http://www.czecos.cz/ (accessed on 6 March 2021)) and FLUXNET (Flux Tower Network; https://fluxnet.fluxdata.org/ (accessed on 6 March 2021)) station network. The FLUXNET station IDs for the wet and dry spruce forest are CZ-BK1 and CZ-RAJ, respectively. CZ-BK1 is also a candidate for ICOS (Integrated Carbon Observation System; https://www.icos-cp.eu/ (accessed on 6 March 2021)) network. The main vegetation cover in both stations is the even-aged stand of Norway spruce. Their main characteristics are presented in Table 1. CZ-BK1 is situated in the Moravian-Silesian Beskids Mountains in the Czech Republic and has a cold and humid climate. CZ-RAJ is situated in the Drahany Highland at a lower altitude and has a moderately warmer and drier climate than CZ-BK1. The reference evapotranspiration amount at each forest stand was computed from in situ data and was additionally used for quantifying the atmospheric evaporative demand from 2014 to 2016 (Table 1).

2.2. Eddy Covariance and Ancillary Measurements

At each station, the EC system consisted of an infrared gas analyser (LI-COR, Lincoln, NE, U.S.A.) and an ultrasonic anemometer (Gill, Hampshire, U.K.) measuring at a 20 Hz frequency (models are specified in Table 2). Each system was mounted on a tower several meters above the forest canopy. The EC measurements were complemented by an extensive set of sensors to collect the required auxiliary meteorological data. Measurements included the air temperature (Tair) and humidity profile with EMS33 temperature and humidity sensors (EMS, Brno, Czech Republic) and precipitation using Precipitation Gauge 386C (Met One Instruments, Grants Pass, OR, U.S.A.). Measurements for the incoming photosynthetic active radiation was made using LI-190R Quantum Sensor (LI-COR, Lincoln, NE, U.S.A.) at Bílý Kříž and EMS12 sensor (EMS, Brno, Czech Republic) at Rájec. Additionally, profiles of soil moisture were measured by using the CS616 (Campbell Scientific, Inc., Logan, UT, U.S.A.) sensors. An overall description of the instrumentation at the study sites is given in Table 2.
The VPD for each site was computed from Tair and relative humidity (RH) as follows [39]:
V P D = 100 R H 100 ( S V P )
where SVP (Pa) is the saturated vapour pressure given as follows:
S V P = 610.7 × 107 . 5 T a i r 237.3 + T a i r

2.3. Data Processing and Analysis

2.3.1. Turbulent Flux Measurements

At both study sites, the EC data from the main growing season (May–September) of 2014–2016 were used for the analysis. Gas analyser measurements of the density of the water vapour in the air and ultrasonic anemometer measurements of three-dimensional wind components and sonic temperature were made [35,36]. The eddy covariance processing was performed, using an open-source software, EddyPro 7.0.6 (Li-COR, Lincoln, NE, U.S.A.). The most recent methods for flux corrections, conversions, and a thorough quality control scheme as proposed by [40] were also applied. This process involves despiking and raw data statistical screening, basic quality checking (QC) of turbulent fluxes (flux stationarity and integral turbulence characteristics tests), coordinate rotation using the planar fit method [41], spectral correction [42,43,44], detecting and compensating for time lags of signals from the ultrasonic anemometer and the gas analyser, footprint estimations and calculating half-hourly fluxes. In EC data post-processing, determining periods with low turbulent mixing is a critical step [35]. Flux measurements over periods with insufficiently developed turbulence, i.e., low friction velocity (u*) need to be detected and filtered out. This filtering procedure assured the exclusion of CO 2 measurements not representative of the biotic flux, i.e., net ecosystem exchange (NEE) [31,45,46].
The R software [47] package ‘REddyProc’ [48] was used to gap-fill the EC data at both sites, using marginal distribution sampling (MDS) [46]. This way, all missing values were replaced by the average value estimated under similar meteorological conditions within a specific time window. If no similar conditions were available within the starting time window of about 7 days, the window was increased. NEE was partitioned into GPP and ecosystem respiration (Reco). The flux partitioning approach by [49] using daytime data was applied to estimate half-hourly GPP ( μ mol m−2 s−1) values. The half-hourly GPP values were aggregated to obtain daily and monthly GPP sums.

2.3.2. Soil Water Content Simulations

Since the soil volumetric water content (SVWC) was not measured throughout the entire period of our analyses (the measurements started at the beginning of 2016), we used a simulated daily SVWC instead. We applied the soil water balance model R-4ET, an R package for Empirical Estimate of Ecosystem EvapoTranspiration as used in [50]. The calibration of the soil water balance model was carried out using Bayesian statistics implemented in the R package BayesianTools [51] with a Differential-Evolution Markov chain Monte Carlo sampler. In Bayesian calibration, the model input parameters are iteratively updated to provide a probability distribution of the calibrated parameters that represents the uncertainty in the measured data and in the model structure. The simulations were repeated 4.8 × 106 times, where the first 1.8 × 106 runs were based on the prior distribution, while the remaining 3 × 106 were constrained by the posterior distribution resulting from the first set of runs. A number of iterations treated as burn-in were set to 0, and the thinning parameter determining the interval in which values are recorded was set to 1. The high number of iterations was necessary to attain a good input parameter convergence with a narrow distribution. For inspecting the convergence, we used Gelman diagnostics with a criterion for the potential scale reduction factor being less than 1.2 for all parameters [52,53]. The final selection of parameters from their probability distributions was conducted by using a maximum a posteriori probability estimate, i.e., the value representing the mode of the posterior distribution.
The main model input variables included meteorological data and the leaf area index (see [50], for more details). The soil was stratified into 12 layers, which increased in thickness with increasing depth down to 3 m. Note that the soil layers were selected in the way that the layers used for optimization matched the depths of all sensors with an extra ±2.5 cm, considering the volume measured by the sensors. At both sites, the available SVWC measurements were available throughout the main part of the root zone. At CZ-BK1, the sensors CS616 (Campbell Scientific, Inc., Logan, UT, U.S.A.) were placed at 0.05, 0.1, 0.22, 0.34, and 0.42 m soil depths and measured since 2016. At CZ-RAJ, the same type of sensors were also placed at 0.05, 0.1, 0.2, 0.5, and 0.8 m depths and measured since 2016. The soil parameters, including SVWC at saturation, field capacity, wilting point, and saturated hydraulic conductivity, were optimized at all 12 depths [50]. Additional single parameters relevant for the SVWC simulations that were optimized included the following: rooting depth with the Beta parameter describing the root profile shape [54], surface resistance and the degree of isohydricity [55], the water interception capacity of leaf and bark area, and curve number representing a runoff parameter [50]. The uniform distribution of priors was used, where the lower and upper limits were set to be within ±50% of the values based on field measurements of the wilting point, field capacity, and saturated water content and ±100% for the remaining parameters that were estimated from the literature or from previous anecdotal analyses. To provide the model with sufficient spin up time to stabilize and ensure reliable and robust parametrization, the simulation was initiated at the beginning of 2010 with initial conditions of SVWC set to the field capacity estimated from soil texture (note that 2010 was one of the wettest years at both sites, according to a computed standardized precipitation-evapotranspiration index over the region). The overall simulation was done for the period 2010–2019 where the observed data for Bayesian calibration spanned the period 2016–2019.
The simulated SVWC averaged over all depths yielded a root mean square error of 0.037 m3 m−3 and 0.017 m3 m−3 for CZ-BK1 and CZ-RAJ, respectively, suggesting realistic SVWC estimates.

2.3.3. Drought Stress Determination

According to [56], drought refers to conditions characterized by low available soil moisture and high atmospheric VPD. As such, in quantifying the dry periods in our study, we used the Standardised Precipitation–Evapotranspiration Index (SPEI), which takes into account both precipitation and potential evapotranspiration for the determination of drought over the main growing season. The SPEI was calculated for various lags (1, 3, 6, 12, and 24 months) from monthly records. We used the period 1981–2010 as the baseline period for the SPEI computation. The computation of SPEI was performed according to [57,58], using the R package SPEI [59]. The SPEI represents an anomaly in the climatological water balance, which is given by the difference between precipitation and reference evapotranspiration. The input data for the SPEI derivation were obtained from climate reanalyses ERA5 monthly averages available at 0.25 spatial resolution, where the reference evapotranspiration was computed from air temperature at 2 m, air dew point temperature in 2 m, wind speed in 2 m converted from the original 10 m height using a logarithmic profile law, and shortwave incoming radiation, following the methodology by [60]. In determining the occurrence of dry and wet events in our study, SPEI classification based on [58] was used (Table 3 and Table 4).

2.3.4. Light Response Curve Fitting

Light response curves (LRC) of daytime GPP were fitted at half-hourly time resolution using the logistic sigmoid approach by [61]:
GPP = 2 × GPP max ( 0.5 1 1 + e x p ( 2 α PAR GPP max ) )
where PAR is the photosynthetically active radiation, α is the apparent quantum yield in mol (CO2) mol−1 (phot.) that describes the effectivity of photosynthesis at low light conditions, and the GPPmax is the asymptotic maximum assimilation rate at the light saturation point in μ mol m2 s−1. The fitting was done separately for half-hourly values from years with normal conditions (2014 and 2016) and the year with drought stress conditions (2015). To eliminate night-time measurements, we used a PAR threshold 10 μ , i.e., light compensation irradiance of 10 μ mol m2 s−1, representing light compensation point of NEE.
Since VPD and soil moisture affect GPP [49,62], the effects of VPD and SVWC on the LRC residuals (GPP values after the removal of the strongest PAR dependence) were analysed for both spruce forest sites within the period under study. A piecewise regression was performed on the LRC residuals versus VPD and SVWC to determine the site-specific environmental stress thresholds at which GPP declined during the growing season periods of normal and drought affected years.

2.3.5. Piecewise Regression Analyses for the Assessment of Drought Effect

Since the LRC model does not account for changes in both VPD and SVWC, piecewise regression of residuals against these environmental variables were analysed. This allowed to determine the response of the studied spruce forest ecosystems with contrasting climates to severe drought conditions over the growing season for normal and drought-affected years. A piecewise regression analysis as part of the R package in ‘segmented’ [63] and the Davies test [64] were used to detect the breakpoints in the regression and also to test for the significant differences in slope parameters of a plot of residuals from the LRC model to drought conditions (high VPD and low SVWC).

3. Results

3.1. Differences in Meteorological Conditions at the Experimental Sites

During May–September from 2014 to 2016, the year 2015 was exceptionally dry compared to 2014 and 2016 (Table 4). Additionally, 2015 was characterized by significantly higher mean VPD values (p < 0.001; Table 5) during the main growing season as compared to the same period in 2014 and 2016 across both sites. Hence, the mean VPD and SVWC values were the best meteorological indicators of drought, while Tair was comparable for years 2015 and 2016. Consequently, the high atmospheric evaporative demand (high VPD values) and low soil moisture (low SVWC values) recorded in 2015 across both sites during the study period further indicated the dryness stress experienced in the main growing season of 2015.
There were observed statistically significant differences (p < 0.05; Mann–Whitney–Wilcoxon ranksum test) in the VPD and SVWC values between the normal and drought regimes across both sites (Figure 1). In addition, VPD was higher and SVWC lower at CZ-RAJ than at CZ-BK1, especially during the drought-affected year (DY).

3.2. Light Response Curves of GPP at Different Climates

The parameters of LRC for normal years (NY) 2014 and 2016 were compared with those for DY for both experimental forest sites (Figure 2; Table 6). Generally, the apparent quantum yield ( α ) at CZ-BK1 was found to be higher than that in CZ-RAJ. During the years with normal conditions, α in CZ-BK1 was observed to be 12% higher than in CZ-RAJ when PAR was less than 500 μ mol m−2 s−1. Moreover, during the DY of 2015, α at CZ-RAJ further declined by 25% as compared to that in CZ-BK1.
The maximum gross primary production at light saturation (GPPmax) for CZ-BK1 was found to be higher than that for CZ-RAJ during the entire study period. During the NY period, GPPmax in CZ-BK1 was 34% higher than in CZ-RAJ. However, during the dry year period, there were significant reductions in the GPPmax values recorded at both spruce forest stations with (18% decline in CZ-BK1 and 17% decline in CZ-RAJ) as compared to those recorded for the years with normal conditions. Additionally, it was further observed that the GPPmax value recorded at CZ-RAJ during the DY period was 33% lower than that in CZ-BK1. Interestingly, the GPPmax value recorded in CZ-BK1 during the DY period was still higher than the GPPmax value recorded at CZ-RAJ during the years with normal conditions.

3.3. Response of Light Response Curve Residuals to VPD and SVWC at Different Climates

The applied LRC model reflects only the light response of GPP and thus, its residuals can be used to assess the influence of additional meteorological parameters (i.e., VPD and SVWC) that are known to affect GPP (Figure 3 and Figure 4). The residual analysis confirmed that the residuals consistently tended to be more negative with increasing VPD and decreasing SVWC (dry conditions). VPD had a significant and stronger effect on GPP than SVWC across both sites, as seen from the piecewise regression analysis, using the residuals from the LRC for the years under study. However, during the DY period, both VPD and SVWC had significant effects on GPP at CZ-RAJ. For CZ-BK1, there was a minimal effect of SVWC on GPP during the years with normal conditions.
Generally, the relationship between the residuals and changes in VPD and SVWC revealed a biphasic response to drought, except for the DY period in CZ-RAJ in Figure 3. All the breakpoints from the piecewise regression were found to be statistically significant (Table 7). However, steeper slopes after the VPD breakpoint values from the initial slope were mostly observed for both DY and NY periods in CZ-RAJ and only for DY period in CZ-BK1. This shows that during all the years under study, GPP decreased at a faster rate after the breakpoint in CZ-RAJ than in CZ-BK1, except for the DY period when there was a significant decline in GPP at a faster rate after the breakpoint across both sites.

3.4. Impact of Drought on Ecosystem Carbon Fluxes at Monthly Timescale

The monthly averages of daily sums of GPP and Reco were analysed separately to assess whether the observed net ecosystem productivity (NEP) reduction during drought was mainly a consequence of increased Reco, reduced GPP, or contribution of both across the forest ecosystems (Figure 5, Figure 6 and Figure 7).
It was observed that there was a significant decline in the total GPP by 14% and 6% during the main growing season period of the dry year for CZ-RAJ and CZ-BK1, respectively. There were also observed significant statistical differences (p < 0.001) in the mean monthly GPP values recorded for July and August during the dry year as compared to the years characterized by normal conditions at both forest sites, especially for CZ-RAJ (Figure 5). Moreover, this significant decline in the mean monthly GPP values at both spruce forest stands during the months of July and August coincided with high daily mean VPD values (>12 hPa) and low SVWC values (<0.16 m3 m−3) as compared to the same period for the adjacent years with normal conditions.
However, an increase in the monthly mean Reco values was observed for July to September of the dry year as compared to the adjacent normal years within CZ-BK1 but with no significant statistical differences, except for the month of August (Figure 6). These three months within 2015 were also characterized by high mean Tair (>20 C) and low mean SVWC (<0.19 m3 m−3) as compared to the two other adjacent years. In addition, during the DY period within CZ-RAJ, the monthly mean Reco values were found to have declined significantly, as compared to the adjacent NY periods. We also observed a 38% significant decline in the mean monthly NEP values during the dry year period for the dry spruce forest in Rájec, compared to a 12% decrease in the mean monthly NEP values for the humid spruce forest in Bílý Kříž (Figure 7). There were also statistically significant differences (p < 0.001) in the mean monthly NEP for the months of July and August across both spruce forest stations. During these months, the observed mean monthly NEP values largely declined during the DY period, as compared to the same period during the NY, especially in CZ-RAJ.
To summarize, the large decline in GPP, Reco and NEP during the dry year period (especially from July to August) for CZ-RAJ showed that the impact of the drought was more severe in CZ-RAJ than in CZ-BK1. However, during July–August of the DY period, Reco at CZ-BK1 significantly increased as compared to that in NY.

4. Discussion

The study sought to assess the different effects of the 2015 summer drought on wet and dry spruce forest ecosystems. Our findings corroborate the earlier published results by [27,28,29] that the year 2015 was characterized by severe drought conditions across Europe (Figure 1 and Table 4). Furthermore, it was found that the months of July–August of 2015 were the most affected by high mean Tair (>20 C), high mean VPD (>10 hPa) and low mean SVWC values (<0.19 m3 m−3) across both spruce forest stands. This shows that the drought effect in 2015 was severe in the months of July–August. Moreover, the results of this study also show that the rate of the forest ecosystem photosynthesis was significantly reduced (likely through the immediate closure of the stomata to protect the tree from desiccation) by the high VPD, Tair and soil water deficit during these two months in 2015. As such, there was an observed strong decline in both forest GPP and NEP across both spruce forest stands during the DY period of 2015, especially in CZ-RAJ compared to CZ-BK1 (M). This strong decline in forest GPP and NEP at CZ-RAJ was mainly due to the high atmospheric evaporative demand coupled with the low SVWC conditions experienced at this forest ecosystem, as compared to that in CZ-BK1 [30]. However, due to the humid climatic conditions at CZ-BK1, an increase in Tair during the months of July–August in 2015 only aided the increase in the kinetics of the enzymes participating in microbial decomposition and root respiration under warm conditions, thereby increasing the overall forest ecosystem respiration [65,66,67,68,69,70,71].
Additionally, results from the piecewise regression analysis using the residuals from the LRC model highlight the effects of both VPD and SVWC on forest ecosystem GPP. This is because both photosynthesis and transpiration are mediated by stomatal conductance, which are affected by these environmental variables. However, a steeper decline in forest ecosystem GPP (Figure 3 and Table 7) was observed with high VPD values across both forest stands, even under non-drought years (when SVWC was non-limiting). This is consistent with recent studies, which show that high VPD values aggravate drought effects in forests, due to the abrupt changes over very short timescales within a day, even without dry soil conditions [7,72,73]. This further explains the strong suppression of GPP by high VPD on even non-drought years (2014 and 2016), as the SPEI (for determining dryness) only captures changes in the contributing environmental drivers (temperature and soil moisture) for longer time scales of weeks or months [74,75]. Thus, we would recommend to include the influence of these abrupt changes in VPD on the rate of photosynthesis with SVWC limitations in future LRC models when analysing the impact of drought on GPP, especially for different forest ecosystems that are exposed to regular strong edaphic droughts [49,62,76,77].

5. Conclusions

The study shows that atmospheric constraints increase the vulnerability of the Norway spruce forest to the severity of drought, especially at sites with a moderately dry climate that are characterized by precipitation that is typically equal or smaller than the atmospheric evaporative demand. It also further highlighted the strong influence of VPD on carbon uptake, which was further worsened by the decline in soil moisture. The effect of SVWC on GPP was especially noticeable during severe drought conditions within the DY period. Consequently, with regards to climate change, our results suggest that elevated temperatures will further exacerbate the drought impacts on forest (Norway spruce) ecosystems at sites with precipitation levels equal or smaller than the atmospheric evaporative demand, such as CZ-RAJ. The decline in GPP and NEP in 2015 found in our study thus questions not only the sustainable productivity, but also the existence of Norway spruce per se in such areas, considering the prolonged period of drought in future climatic conditions. The results of this study may help decision makers to quantitatively assess the performance of Norway spruce in future climatic conditions.

Author Contributions

Conceptualisation, C.M. and L.Š.; methodology, C.M., L.Š., M.F., M.P. and M.V.M.; software, C.M., L.K. and L.Š.; validation, M.F., M.A., G.J., J.D.M. and N.K., E.K.N.; formal analysis, C.M. and L.Š.; investigation, C.M. and L.Š.; resources, M.P. and L.Š.; data curation, L.Š. and G.J.; writing—original draft preparation, C.M.; writing—review and editing, L.Š., M.F., M.A., G.J., L.F., J.D.M., N.K., E.K.N. and M.V.M.; visualisation, L.K.; supervision, L.Š. and M.V.M.; project administration, M.P. and L.Š.; funding acquisition, M.P. and M.V.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC was funded by the Ministry of Education, Youth and Sports of the Czech Republic (CR) within the CzeCOS program, grant number LM2018123.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on FLUXNET and could also be made available on request from the corresponding author.

Acknowledgments

The authors wish to thank Radek Czerný who assisted in the eddy covariance data processing. This work was supported by the Ministry of Education, Youth and Sports of the Czech Republic (CR) within the CzeCOS program, grant number LM2018123. L.Š. was supported by the Ministry of Education, Youth and Sports of CR within Mobility CzechGlobe2 (CZ.02.2.69/0.0/0.0/18 053/0016924). M.F., L.Š. and L.F. acknowledges the support by the Ministry of Education, Youth and Sports of the Czech Republic for SustES—Adaptation strategies for sustainable ecosystem services and food security under adverse environmental conditions (CZ.02.1.01/0.0/0.0/16 019/0000797). M.V.M. acknowledges the support by VEGA project 2/0013/17 (Slovak Agency for Scientific Support): The role of ecosystem services in support of landscape conservation under the global change.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Rossi, G.; Cancelliere, A.; Pereira, L.S. (Eds.) Tools for Drought Mitigation in Mediterranean Regions; Springer Science and Business Media: Berlin, Germany, 2003; Volume 44. [Google Scholar]
  2. Huntington, T.G. Evidence for intensification of the global water cycle: Review and synthesis. J. Hydrol. 2003, 319, 83–95. [Google Scholar] [CrossRef]
  3. Bonan, G.B. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 2008, 320, 1444–1449. [Google Scholar] [CrossRef] [Green Version]
  4. Allen, C.D.; Macalady, A.K.; Chenchouni, H.; Bachelet, D.; McDowell, N.; Vennetier, M.; Gonzalez, P. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For. Ecol. Manag. 2010, 259, 660–684. [Google Scholar] [CrossRef] [Green Version]
  5. Marek, M.V.; Janouš, D.; Taufarová, K.; Havránková, K.; Pavelka, M.; Kaplan, V.; Marková, I. Carbon exchange between ecosystems and atmosphere in the Czech Republic is affected by climate factors. Environ. Pollut. 2011, 159, 1035–1039. [Google Scholar] [CrossRef]
  6. Lal, R.; Smith, P.; Jungkunst, H.F.; Mitsch, W.J.; Lehmann, J.; Nair, P.R.; Skorupa, A.L. The carbon sequestration potential of terrestrial ecosystems. J. Soil Water Conserv. 2018, 73, 145A–152A. [Google Scholar] [CrossRef] [Green Version]
  7. Novick, K.A.; Ficklin, D.L.; Stoy, P.C.; Williams, C.A.; Bohrer, G.; Oishi, A.C.; Papuga, S.A.; Blanken, P.D.; Noormets, A.; Sulman, B.N.; et al. The increasing importance of atmospheric demand for ecosystem water and carbon fluxes. Nat. Clim. Chang. 2016, 6, 1023–1027. [Google Scholar] [CrossRef] [Green Version]
  8. Kowalska, N.; Šigut, L.; Stojanović, M.; Fischer, M.; Kyselova, I.; Pavelka, M. Analysis of floodplain forest sensitivity to drought. Philos. Trans. R. Soc. B 2020, 375, 20190518. [Google Scholar] [CrossRef]
  9. Beer, C.; Reichstein, M.; Tomelleri, E.; Ciais, P.; Jung, M.; Carvalhais, N.; Rödenbeck, C.; Arain, M.A.; Baldocchi, D.; Bonan, G.B.; et al. Terrestrial gross carbon dioxide uptake: Global distribution and covariation with climate. Science 2010, 329, 834–838. [Google Scholar] [CrossRef] [Green Version]
  10. Pan, Y.; Birdsey, R.A.; Fang, J.; Houghton, R.; Kauppi, P.E.; Kurz, W.A.; Phillips, O.L.; Shvidenko, A.; Lewis, S.L.; Canadell, J.G.; et al. A large and persistent carbon sink in the world’s forests. Science 2011, 333, 988–993. [Google Scholar] [CrossRef] [Green Version]
  11. Cox, P.M.; Pearson, D.; Booth, B.B.; Friedlingstein, P.; Huntingford, C.; Jones, C.D.; Luke, C.M. Sensitivity of tropical carbon to climate change constrained by carbon dioxide variability. Nature 2013, 494, 341. [Google Scholar] [CrossRef]
  12. Le Quéré, C.; Andres, R.J.; Boden, T.; Conway, T.; Houghton, R.A.; House, J.I.; Marland, G.; Peters, G.P.; Van der Werf, G.; Ahlström, A.; et al. The global carbon budget 1959–2011. Earth Syst. Sci. Data 2013, 5, 165–185. [Google Scholar] [CrossRef] [Green Version]
  13. Krejza, J.; Cienciala, E.; Světlík, J.; Bellan, M.; Noyer, E.; Horáček, P.; Štěpánek, P.; Marek, M.V. Evidence of climate-induced stress of Norway spruce along elevation gradient preceding the current dieback in Central Europe. Trees 2020, 35, 103–119. [Google Scholar] [CrossRef]
  14. Zeng, N.; Qian, H.; Roedenbeck, C.; Heimann, M. Impact of 1998–2002 midlatitude drought and warming on terrestrial ecosystem and the global carbon cycle. Geophys. Res. Lett. 2005, 32. [Google Scholar] [CrossRef] [Green Version]
  15. Dai, A. Increasing drought under global warming in observations and models. Nat. Clim. Chang. 2013, 3, 52. [Google Scholar] [CrossRef]
  16. Reichstein, M.; Bahn, M.; Ciais, P.; Frank, D.; Mahecha, M.D.; Seneviratne, S.I.; Zscheischler, J.; Beer, C.; Buchmann, N.; Frank, D.C.; et al. Climate extremes and the carbon cycle. Nature 2013, 500, 287. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, L.; Chen, W.; Zhou, W. Assessment of future drought in Southwest China based on CMIP5 multimodel projections. Adv. Atmos. Sci. 2014, 31, 1035–1050. [Google Scholar] [CrossRef]
  18. Xie, Z.; Wang, L.; Jia, B.; Yuan, X. Measuring and modeling the impact of a severe drought on terrestrial ecosystem CO2 and water fluxes in a subtropical forest. J. Geophys. Res. Biogeosci. 2016, 121, 2576–2587. [Google Scholar] [CrossRef] [Green Version]
  19. Zang, C.; Hartl-Meier, C.; Dittmar, C.; Rothe, A.; Menzel, A. Patterns of drought tolerance in major European temperate forest trees: Climatic drivers and levels of variability. Glob. Chang. Biol. 2014, 20, 3767–3779. [Google Scholar] [CrossRef]
  20. De Boeck, H.J.; Verbeeck, H. Drought-associated changes in climate and their relevance for ecosystem experiments and models. Biogeosciences 2011, 8, 1121–1130. [Google Scholar] [CrossRef] [Green Version]
  21. Hui, D.; Deng, Q.; Tian, H.; Luo, Y. Climate change and carbon sequestration in forest ecosystems. In Handbook of Climate Change Mitigation and Adaptation; Chen, W.Y., Suzuki, T., Lackner, M., Eds.; Springer International Publishing: New York, NY, USA, 2017; pp. 555–594. [Google Scholar]
  22. Qie, L.; Lewis, S.L.; Sullivan, M.J.; Lopez-Gonzalez, G.; Pickavance, G.C.; Sunderland, T.; Ashton, P.; Hubau, W.; Salim, K.A.; Aiba, S.I.; et al. Long-term carbon sink in Borneo’s forests halted by drought and vulnerable to edge effects. Nat. Commun. 2017, 8, 1966. [Google Scholar] [CrossRef] [Green Version]
  23. Ciais, P.; Reichstein, M.; Viovy, N.; Granier, A.; Ogée, J.; Allard, V.; Aubinet, M.; Buchmann, N.; Bernhofer, C.; Carrara, A.; et al. Europe-wide reduction in primary productivity caused by the heat and drought in 2003. Nature 2005, 437, 529. [Google Scholar] [CrossRef]
  24. Granier, A.; Reichstein, M.; Breda, N.; Janssens, I.A.; Falge, E.; Ciais, P.; Grünwald, T.; Aubinet, M.; Berbigier, P.; Bernhofer, C.; et al. Evidence for soil water control on carbon and water dynamics in European forests during the extremely dry year: 2003. Agric. For. Meteorol. 2007, 143, 123–145. [Google Scholar] [CrossRef]
  25. Trnka, M.; Hlavinka, P.; Možný, M.; Semerádová, D.; Štěpánek, P.; Balek, J.; Bartošová, L.; Zahradníček, P.; Bláhová, M.; Skalak, P.; et al. Czech Drought Monitor System for Monitoring and Forecasting Agricultural Drought and Drought Impacts. Int. J. Climatol. 2020, 40, 5941–5958. [Google Scholar] [CrossRef]
  26. Orth, R.; Seneviratne, S.I. Introduction of a simple-model-based land surface dataset for Europe. Environ. Res. Lett. 2015, 10, 0442012. [Google Scholar]
  27. Orth, R.; Zscheischler, J.; Seneviratne, S.I. Record dry summer in 2015 challenges precipitation projections in Central Europe. Sci. Rep. 2016, 6, 28334. [Google Scholar] [CrossRef] [PubMed]
  28. Van Lanen, H.A.; Laaha, G.; Kingston, D.G.; Gauster, T.; Ionita, M.; Vidal, J.P.; Vlnas, R.; Tallaksen, L.M.; Stahl, K.; Hannaford, J.; et al. Hydrology needed to manage droughts: The 2015 European case. Hydrol. Process. 2016, 30, 3097–3104. [Google Scholar] [CrossRef] [Green Version]
  29. Ionita, M.; Tallaksen, L.M.; Kingston, D.G.; Stagge, J.H.; Laaha, G.; Van Lanen, H.A.; Scholz, P.; Chelcea, S.M.; Haslinger, K. The European 2015 drought from a climatological perspective. Hydrol. Earth Syst. Sci. 2017, 21, 1397–1419. [Google Scholar] [CrossRef] [Green Version]
  30. Krupková, L.; Havránková, K.; Krejza, J.; Sedlák, P.; Marek, M.V. Impact of water scarcity on spruce and beech forests. J. For. Res. 2018, 30, 899–909. [Google Scholar] [CrossRef]
  31. McGloin, R.; Šigut, L.; Havránková, K.; Dušek, J.; Pavelka, M.; Sedlák, P. Energy balance closure at a variety of ecosystems in Central Europe with contrasting topographies. Agric. For. Meteorol. 2018, 248, 418–431. [Google Scholar] [CrossRef]
  32. Kenderes, K.; Mihók, B.; Standovár, T. Thirty years of gap dynamics in a Central European beech forest reserve. Forestry 2008, 81, 111–123. [Google Scholar] [CrossRef]
  33. Nadezhdina, N.; Urban, J.; Čermák, J.; Nadezhdin, V.; Kantor, P. Comparative study of long-term water uptake of Norway spruce and Douglas-fir in Moravian upland. J. Hydrol. Hydromech. 2014, 62, 1–6. [Google Scholar] [CrossRef]
  34. Kmet, J.; Ditmarová, L.; Kurjak, D.; Priwitzer, T. Physiological response of Norway spruce foliage in the drought vegetation period 2009. Beskydy 2011, 4, 109–118. [Google Scholar]
  35. Aubinet, M.; Vesala, T.; Papale, D. (Eds.) Eddy Covariance: A Practical Huide to Measurement and Data Analysis; Springer Science and Business Media: Berlin, Germany, 2012. [Google Scholar]
  36. Baldocchi, D.D. Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: Past, present and future. Glob. Chang. Biol. 2003, 9, 479–492. [Google Scholar] [CrossRef] [Green Version]
  37. Sedlák, P.; Aubinet, M.; Heinesch, B.; Janouš, D.; Pavelka, M.; Potužníková, K.; Yernaux, M. Night-time airflow in a forest canopy near a mountain crest. Agric. For. Meteorol. 2010, 150, 736–744. [Google Scholar] [CrossRef]
  38. Urban, O.; Klem, K.; Holišová, P.; Šigut, L.; Šprtová, M.; Teslová-Navrátilová, P.; Zitová, M.; Špunda, V.; Marek, M.V.; Grace, J. Impact of elevated CO2 concentration on dynamics of leaf photosynthesis in Fagus sylvatica is modulated by sky conditions. Environ. Pollut. 2014, 185, 271–280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Monteith, J.; Unsworth, M. Principles of Environmental Physics: Plants, Animals, and the Atmosphere; Academic Press: Cambridge, MA, USA, 2013. [Google Scholar]
  40. Foken, T.; Leuning, R.; Oncley, S.R.; Mauder, M.; Aubinet, M. Corrections and data quality control. In Eddy Covariance; Springer: Dordrecht, The Netherlands, 2012; pp. 85–131. [Google Scholar]
  41. Wilczak, J.M.; Oncley, S.P.; Stage, S.A. Sonic anemometer tilt correction algorithms. Bound. Layer Meteorol. 2001, 99, 127–150. [Google Scholar] [CrossRef]
  42. Ibrom, A.; Dellwik, E.; Larsen, S.E.; Pilegaard, K.I.M. On the use of the Webb-Pearman-Leuning theory for closed-path eddy correlation measurements. Tellus B Chem. Phys. Meteorol. 2007, 59, 937–946. [Google Scholar] [CrossRef]
  43. Horst, T.W.; Lenschow, D.H. Attenuation of scalar fluxes measured with spatially-displaced sensors. Bound. Layer Meteorol. 2009, 130, 275–300. [Google Scholar] [CrossRef]
  44. Moncrieff, J.B.; Jarvis, P.G.; Valentini, R. Canopy fluxes. In Methods in Ecosystem Science; Springer: New York, NY, USA, 2000; pp. 161–180. [Google Scholar]
  45. Papale, D.; Reichstein, M.; Aubinet, M.; Canfora, E.; Bernhofer, C.; Kutsch, W.; Longdoz, B.; Rambal, S.; Valentini, R.; Vesala, T.; et al. Towards a standardized processing of Net Ecosystem Exchange measured with eddy covariance technique: Algorithms and uncertainty estimation. Biogeosciences 2006, 3, 571–583. [Google Scholar] [CrossRef] [Green Version]
  46. Reichstein, M.; Falge, E.; Baldocchi, D.; Papale, D.; Aubinet, M.; Berbigier, P.; Bernhofer, C.; Buchmann, N.; Gilmanov, T.; Granier, A.; et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: Review and improved algorithm. Glob. Chang. Biol. 2005, 11, 1424–1439. [Google Scholar] [CrossRef]
  47. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  48. Wutzler, T.; Lucas-Moffat, A.; Migliavacca, M.; Knauer, J.; Sickel, K.; Šigut, L.; Menzer, O.; Reichstein, M. Basic and extensible post-processing of eddy covariance flux data with REddyProc. Biogeosciences 2018, 15, 5015–5030. [Google Scholar] [CrossRef] [Green Version]
  49. Lasslop, G.; Reichstein, M.; Papale, D.; Richardson, A.D.; Arneth, A.; Barr, A.; Stoy, P.; Wohlfahrt, G. Separation of net ecosystem exchange into assimilation and respiration using a light response curve approach: Critical issues and global evaluation. Glob. Chang. Biol. 2010, 16, 187–208. [Google Scholar] [CrossRef] [Green Version]
  50. Fischer, M.; Zenone, T.; Trnka, M.; Orság, M.; Montagnani, L.; Ward, E.J.; Tripathi, A.M.; Hlavinka, P.; Seufert, G.; Žalud, Z.; et al. Water requirements of short rotation poplar coppice: Experimental and modelling analyses across Europe. Agric. For. Meteorol. 2018, 250, 343–360. [Google Scholar] [CrossRef]
  51. Hartig, F.; Minunno, F.; Paul, S. BayesianTools: General-Purpose MCMC and SMC Samplers and Tools for Bayesian Statistics. R Package Version 0.1.6. 2019. Available online: https://CRAN.R-project.org/package=BayesianTools (accessed on 20 March 2021).
  52. Gelman, A.; Rubin, D.B. Inference from iterative simulation using multiple sequences. Stat. Sci. 1992, 7, 457–472. [Google Scholar] [CrossRef]
  53. Brooks, S.P.; Gelman, A. General methods for monitoring convergence of iterative simulations. J. Comput. Graph. Stat. 1998, 7, 434–455. [Google Scholar]
  54. Jackson, R.B.; Canadell, J.; Ehleringer, J.R.; Mooney, H.A.; Sala, O.E.; Schulze, E.D. A global analysis of root distributions for terrestrial biomes. Oecologia 1996, 108, 389–411. [Google Scholar] [CrossRef] [PubMed]
  55. Oren, R.; Sperry, J.S.; Katul, G.G.; Pataki, D.E.; Ewers, B.E.; Phillips, N.; Scháfer, K.V.R. Survey and synthesis of intra-and interspecific variation in stomatal sensitivity to vapour pressure deficit. Plant Cell Environ. 1999, 22, 1515–1526. [Google Scholar] [CrossRef] [Green Version]
  56. Mishra, A.K.; Singh, V.P. A review of drought concepts. J. Hydrol. 2010, 391, 202–216. [Google Scholar] [CrossRef]
  57. Beguería, S.; Vicente-Serrano, S.M.; Reig, F.; Latorre, B. Standardized precipitation evapotranspiration index (SPEI) revisited: Parameter fitting, evapotranspiration models, tools, datasets and drought monitoring. Int. J. Climatol. 2014, 34, 3001–3023. [Google Scholar] [CrossRef] [Green Version]
  58. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A Multiscalar Drought Index Sensitive to Global Warming: The Standardized Precipitation Evapotranspiration Index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef] [Green Version]
  59. Beguería, S.; Vicente-Serrano, S.M. SPEI: Calculation of the Standardised Precipitation-Evapotranspiration Index. R Package Version 1.7. 2017. Available online: https://CRAN.R-project.org/package=SPEI (accessed on 13 January 2021).
  60. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop evapotranspiration—Guidelines for computing crop water requirements. FAO Irrig. Drain. Pap. 1998, 56, 290. [Google Scholar]
  61. Moffat, A.M. A New Methodology to Interpret High Resolution Measurements of Net Carbon Fluxes between Terrestrial Ecosystems and the Atmosphere. Ph.D. Thesis, Friedrich-Schiller-Universität Jena, Fakultät für Mathematik und Informatik, Jena, Germany, 2012. [Google Scholar]
  62. Migliavacca, M.; Reichstein, M.; Richardson, A.D.; Colombo, R.; Sutton, M.A.; Lasslop, G.; Tomelleri, E.; Wohlfahrt, G.; Carvalhais, N.; Cescatti, A.; et al. Semiempirical modeling of abiotic and biotic factors controlling ecosystem respiration across eddy covariance sites. Glob. Chang. Biol. 2011, 17, 390–409. [Google Scholar] [CrossRef]
  63. Muggeo, V.M. Estimating regression models with unknown break-points. Stat. Med. 2003, 22, 3055–3071. [Google Scholar] [CrossRef]
  64. Davies, R.B. Hypothesis testing when a nuisance parameter is present only under the alternative: Linear model case. Biometrika 2002, 89, 484–489. [Google Scholar] [CrossRef]
  65. Reichstein, M.; Ciais, P.; Papale, D.; Valentini, R.; Running, S.; Viovy, N.; Cramer, W.; Granier, A.; Ogee, J.; Allard, V.; et al. Reduction of ecosystem productivity and respiration during the European summer 2003 climate anomaly: A joint flux tower, remote sensing and modelling analysis. Glob. Chang. Biol. 2007, 13, 634–651. [Google Scholar] [CrossRef]
  66. Jassal, R.S.; Black, T.A.; Chen, B.; Roy, R.; Nesic, Z.; Spittlehouse, D.L.; Trofymow, J.A. N2O emissions and carbon sequestration in a nitrogen-fertilized Douglas fir stand. J. Geophys. Res. Biogeosci. 2008, 113. [Google Scholar] [CrossRef] [Green Version]
  67. Mahecha, M.D.; Reichstein, M.; Carvalhais, N.; Lasslop, G.; Lange, H.; Seneviratne, S.I.; Vargas, R.; Ammann, C.; Arain, M.A.; Cescatti, A.; et al. Global convergence in the temperature sensitivity of respiration at ecosystem level. Science 2010, 329, 838–840. [Google Scholar] [CrossRef] [Green Version]
  68. Wu, S.H.; Jansson, P.E.; Kolari, P. Modeling seasonal course of carbon fluxes and evapotranspiration in response to low temperature and moisture in a boreal Scots pine ecosystem. Ecol. Model. 2011, 222, 3103–3119. [Google Scholar] [CrossRef]
  69. Frank, D.; Reichstein, M.; Bahn, M.; Thonicke, K.; Frank, D.; Mahecha, M.D.; Smith, P.; Van der Velde, M.; Vicca, S.; Babst, F.; et al. Effects of climate extremes on the terrestrial carbon cycle: Concepts, processes and potential future impacts. Glob. Chang. Biol. 2015, 21, 2861–2880. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Van Gorsel, E.; Wolf, S.; Cleverly, J.; Isaac, P.; Haverd, V.; Ewenz, C.; Arndt, S.; Beringer, J.; Resco de Dios, V.; Evans, B.J.; et al. Carbon uptake and water use in woodlands and forests in southern Australia during an extreme heat wave event in the “angry Summer” of 2012/2013. Biogeosciences 2016, 13, 5947–5964. [Google Scholar] [CrossRef] [Green Version]
  71. von Buttlar, J.; Zscheischler, J.; Rammig, A.; Sippel, S.; Reichstein, M.; Knohl, A.; Jung, M.; Menzer, O.; Arain, M.A.; Buchmann, N.; et al. Impacts of droughts and extreme-temperature events on gross primary production and ecosystem respiration: A systematic assessment across ecosystems and climate zones. Biogeosciences 2018, 15, 1293–1318. [Google Scholar] [CrossRef] [Green Version]
  72. Ruehr, N.K.; Law, B.E.; Quandt, D.; Williams, M. Effects of heat and drought on carbon and water dynamics in a regenerating semi-arid pine forest: A combined experimental and modeling approach. Biogeosciences 2014, 11, 4139–4156. [Google Scholar] [CrossRef] [Green Version]
  73. Sulman, B.N.; Roman, D.T.; Yi, K.; Wang, L.; Phillips, R.P.; Novick, K.A. High atmospheric demand for water can limit forest carbon uptake and transpiration as severely as dry soil. Geophys. Res. Lett. 2016, 43, 9686–9695. [Google Scholar] [CrossRef]
  74. Sheffield, J.; Wood, E.F.; Roderick, M.L. Little change in global drought over the past 60 years. Nature 2012, 491, 435–438. [Google Scholar] [CrossRef] [PubMed]
  75. Potopová, V.; Boroneant, C.; Možný, M.; Soukup, J. Driving role of snow cover on soil moisture and drought development during the growing season in the Czech Republic. Int. J. Climatol. 2016, 36, 3741–3758. [Google Scholar] [CrossRef] [Green Version]
  76. Laaha, G.; Gauster, T.; Tallaksen, L.M.; Vidal, J.P.; Stahl, K.; Prudhomme, C.; Heudorfer, B.; Vlnas, R.; Ionita, M.; Van Lanen, H.A.; et al. The European 2015 drought from a hydrological perspective. Hydrol. Earth Syst. Sci. 2017, 21, 3001–3024. [Google Scholar] [CrossRef] [Green Version]
  77. Pretzsch, H.; Schütze, G.; Biber, P. Drought can favour the growth of small in relation to tall trees in mature stands of Norway spruce and European beech. For. Ecosyst. 2018, 5, 20. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Histogram showing the frequency of occurrence of (a) vapour pressure deficit (VPD) and (b) soil volumetric water content (SVWC) conditions during May–September between years with normal conditions (dark colour) and drought stress (grey colour) in CZ-BK1 and CZ-RAJ.
Figure 1. Histogram showing the frequency of occurrence of (a) vapour pressure deficit (VPD) and (b) soil volumetric water content (SVWC) conditions during May–September between years with normal conditions (dark colour) and drought stress (grey colour) in CZ-BK1 and CZ-RAJ.
Atmosphere 12 00988 g001
Figure 2. Response of gross primary productions to photosynthetically active radiation during the years with normal conditions (black) and affected by drought stress (red) in CZ-BK1 and CZ-RAJ. The half-hourly GPP values (points) were fitted using the logistic sigmoid light response curves (lines).
Figure 2. Response of gross primary productions to photosynthetically active radiation during the years with normal conditions (black) and affected by drought stress (red) in CZ-BK1 and CZ-RAJ. The half-hourly GPP values (points) were fitted using the logistic sigmoid light response curves (lines).
Atmosphere 12 00988 g002
Figure 3. Relationship between the light response curve (LRC) residuals of the light response curve of gross primary production and the vapour pressure deficit (VPD) at the spruce forest sites in CZ-BK1 and CZ-RAJ for the normal (NY) and dry years (DY). Red represents low SVWC conditions and blue shows high SVWC conditions. The grey dashed lines represent the piecewise regression model slope, whereas the black dashed lines show the breakpoint values.
Figure 3. Relationship between the light response curve (LRC) residuals of the light response curve of gross primary production and the vapour pressure deficit (VPD) at the spruce forest sites in CZ-BK1 and CZ-RAJ for the normal (NY) and dry years (DY). Red represents low SVWC conditions and blue shows high SVWC conditions. The grey dashed lines represent the piecewise regression model slope, whereas the black dashed lines show the breakpoint values.
Atmosphere 12 00988 g003
Figure 4. Relationship between the light response curve (LRC) residuals of the light response curve of gross primary production and the soil volumetric water content (SVWC) at the spruce forest sites in CZ-BK1 and CZ-RAJ for the normal (NY) and dry years (DY). Red represents high VPD conditions and blue shows low VPD conditions. The grey dashed lines represent the piecewise regression model slope, whereas the black dashed lines show the breakpoint values.
Figure 4. Relationship between the light response curve (LRC) residuals of the light response curve of gross primary production and the soil volumetric water content (SVWC) at the spruce forest sites in CZ-BK1 and CZ-RAJ for the normal (NY) and dry years (DY). Red represents high VPD conditions and blue shows low VPD conditions. The grey dashed lines represent the piecewise regression model slope, whereas the black dashed lines show the breakpoint values.
Atmosphere 12 00988 g004
Figure 5. Monthly averages of daily sums of gross primary productivity (GPP) for May–September of the normal years (NY) and dry year (DY) in CZ-BK1 and CZ-RAJ. The tables within the figure represent the mean monthly vapour pressure deficit (VPD) and soil volumetric water content (SVWC) values from May to September for each forest station.
Figure 5. Monthly averages of daily sums of gross primary productivity (GPP) for May–September of the normal years (NY) and dry year (DY) in CZ-BK1 and CZ-RAJ. The tables within the figure represent the mean monthly vapour pressure deficit (VPD) and soil volumetric water content (SVWC) values from May to September for each forest station.
Atmosphere 12 00988 g005
Figure 6. Monthly averages of daily sums of ecosystem respiration (Reco) for May–September of the normal years (NY) and dry year (DY) in CZ-BK1 and CZ-RAJ. The tables within the figure represent the mean monthly air temperature (Tair) and soil volumetric water content (SVWC) values from May to September for each forest station.
Figure 6. Monthly averages of daily sums of ecosystem respiration (Reco) for May–September of the normal years (NY) and dry year (DY) in CZ-BK1 and CZ-RAJ. The tables within the figure represent the mean monthly air temperature (Tair) and soil volumetric water content (SVWC) values from May to September for each forest station.
Atmosphere 12 00988 g006
Figure 7. Monthly averages of daily sums of NEP for May–September of the normal years (NY) and dry year (DY) in CZ-BK1 and CZ-RAJ.
Figure 7. Monthly averages of daily sums of NEP for May–September of the normal years (NY) and dry year (DY) in CZ-BK1 and CZ-RAJ.
Atmosphere 12 00988 g007
Table 1. Characteristics of the study sites.
Table 1. Characteristics of the study sites.
Site NameCZ-BK1CZ-RAJ
LocationMoravian-Silesian Beskids MountainsDrahany Highland
Coordinates49 30 08 N, 18 32 13 E49 26 37 N, 16 41 48 E
Elevation (in m a.s.l)875625
TopographyMountainous (13 slope with SSW exposure, located close to a mountain ridge)Hilly (5 slope with NNE exposure)
Ecosystem TypeConiferous evergreen forestConiferous evergreen forest
Prevailing speciesNorway Spruce (Picea abies (L.) H.Karst.)Norway Spruce (Picea abies (L.) H.Karst.)
Canopy height (m)16 (mean, as of 2015)33 (mean, as of 2015)
Stand age (years)35 (as of 2016)113 (as of 2016)
Mean annual air temperature (May–September; C)7.2 (2014–2016)8.3 (2014–2016)
Mean annual precipitation (May–September; mm)1143 (2014–2016)610 (2014–2016)
Mean annual reference evapotranspiration (May–September; mm)569 (2014–2016)649 (2014–2016)
Soil typeHaplic and Entic PodzolModal Cambisol oligotrophic
Max and min fetch (m)717 (WNW) and 115 (ENE)697 (SSW) and 96 (NNW)
References[31,37,38][31]
Table 2. Description of the eddy covariance systems at the investigated sites.
Table 2. Description of the eddy covariance systems at the investigated sites.
Site NameCZ-BK1CZ-RAJ
Ultrasonic Anemometer
InstrumentGill HS-50 ultrasonic anemometer, Gill Instruments, Hampshire, U.K.Gill R3-100 ultrasonic anemometer, Gill Instruments, Hampshire, U.K. but later changed to Gill HS-50 on 5 June 2015
Gas Analyser
InstrumentLI-7200 enclosed gas analyser, LI-COR, Lincoln, U.S.A.Initially LI-7000 (IRG-0226) closed-path gas analyser, LI-COR, Lincoln, U.S.A. but later changed to LI-7200 on 5 June 2015
Measurement Height for the Eddy covariance Set-up (m)Initially 20.5, but later changed to 25 m on 7 June 201641
Air Temperature and Humidity Profile
InstrumentEMS33 temperature and humidity sensor (EMS, Brno, CZ)EMS33 temperature and humidity sensor (EMS, Brno, CZ)
Measurement Height (m)2.0, 7.6, 12.6, 13.5, 14.3, 14.8, 15.4, 16.5, 18.72.0, 11.0, 23.0, 29.0, 35.0, 42.0
Net Radiation
InstrumentCNR1 Net RadiometerCNR1 Net Radiometer
Measurement Height (m)2242
Soil Moisture
InstrumentCS616 (Campbell Scientific, Inc., Logan, UT, U.S.A.)CS616 (Campbell Scientific, Inc., Logan, UT, U.S.A.)
Measurement Depths (m)0.05, 0.1, 0.22, 0.34, 0.420.05, 0.1, 0.2, 0.5, 0.8
Table 3. The Standardised Precipitation-Evapotranspiration Index (SPEI) categories based on the classification of SPEI values by [58].
Table 3. The Standardised Precipitation-Evapotranspiration Index (SPEI) categories based on the classification of SPEI values by [58].
Drought/ Wet SeveritySPEI Value
Extremely Wet≥2.00
Severely Wet1.50–1.99
Moderately Wet1.00–1.49
Near Normal−0.99–(0.99)
Moderate Drought−1.00–(−1.49)
Severe Drought−1.50–(−1.99)
Extreme Drought≤−2.00
Table 4. Categorization of dryness/wetness using Standardised Precipitation-Evapotranspiration Index (SPEI) indices for CZ-BK1 and CZ-RAJ stations in years (2014–2016).
Table 4. Categorization of dryness/wetness using Standardised Precipitation-Evapotranspiration Index (SPEI) indices for CZ-BK1 and CZ-RAJ stations in years (2014–2016).
YEARSCZ-BK1CZ-RAJ
SPEI VALUECLASSSPEI VALUECLASS
20140.94Near Normal0.56Near Normal
2015−1.75Severe Drought−1.55Severe Drought
2016−0.20Near Normal−0.87Near Normal
Table 5. Mean values of vapour pressure deficit (VPD), air temperature (Tair) and soil volumetric water content (SVWC) during the main growing season in years for both spruce forest stands.
Table 5. Mean values of vapour pressure deficit (VPD), air temperature (Tair) and soil volumetric water content (SVWC) during the main growing season in years for both spruce forest stands.
YEARSCZ-BK1CZ-RAJ
VPD (hPa)Tair ( C)SVWC (m3 −3)VPD (hPa)Tair ( C)SVWC (m3 −3)
20147.613.40.288.714.90.21
201510.214.60.1911.215.80.19
20167.714.40.259.316.00.22
Table 6. Light response curve parameters for the normal (2014 and 2016) and dry (2015) years for the wet and dry climates (CZ-BK1 and CZ-RAJ respectively) within the main growing season period of May–September. The apparent quantum yield ( α ), the maximum gross primary production at light saturation (GPPmax) and the coefficient of determination (R2) are also shown.
Table 6. Light response curve parameters for the normal (2014 and 2016) and dry (2015) years for the wet and dry climates (CZ-BK1 and CZ-RAJ respectively) within the main growing season period of May–September. The apparent quantum yield ( α ), the maximum gross primary production at light saturation (GPPmax) and the coefficient of determination (R2) are also shown.
VariantsCZ-BK1CZ-RAJ
Years with Normal Conditions (2014 & 2016)Dry Year (2015)Years with Normal Conditions (2014 & 2016)Dry Year (2015)
α [mol (CO2) mol−1 (phot.)]0.0383 ± 0.00030.0446 ± 0.00060.0338 ± 0.00030.0312 ± 0.0006
GPPmax[ μ mol m−2 s−1]26.91 ± 0.1422.03 ± 0.1317.77 ± 0.0814.75 ± 0.12
R20.880.780.780.73
Table 7. Breakpoints and Slopes from the piecewise regression of the light response curve residuals to vapour pressure deficit (VPD) and soil volumetric water content (SVWC) for the period of May–September.
Table 7. Breakpoints and Slopes from the piecewise regression of the light response curve residuals to vapour pressure deficit (VPD) and soil volumetric water content (SVWC) for the period of May–September.
VariantsCZ-BK1CZ-RAJ
Years with Normal Conditions (2014 & 2016)Dry Year (2015)Years with Normal Conditions (2014 & 2016)Dry Year (2015)
VPD [hPa]2.6 * * * 5.3 * * * 6.8 * * * 23.5 * * *
Slope before breakpoint in VPD1.12 ± 0.300.02 ± 0.120.08 ± 0.05−0.55 ± 0.01
Slope after breakpoint in VPD−0.39 ± 0.02−0.70 ± 0.02−0.54 ± 0.02−0.22 ± 0.06
SVWC [m3 m−3]0.20 * * * 0.16 * * * 0.18 * * * 0.19 * * *
Slope before breakpoint in SVWC−85.49 ± 68.37305.61 ± 23.43164.61 ± 25.14187.22 ± 8.78
Slope after breakpoint in SVWC14.75 ± 2.605.03 ± 2.51−1.27 ± 2.8452.16 ± 6.24
Signif. code for the breakpoint values: p < 0.001 ‘ * * * ’.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mensah, C.; Šigut, L.; Fischer, M.; Foltýnová, L.; Jocher, G.; Acosta, M.; Kowalska, N.; Kokrda, L.; Pavelka, M.; Marshall, J.D.; et al. Assessing the Contrasting Effects of the Exceptional 2015 Drought on the Carbon Dynamics in Two Norway Spruce Forest Ecosystems. Atmosphere 2021, 12, 988. https://doi.org/10.3390/atmos12080988

AMA Style

Mensah C, Šigut L, Fischer M, Foltýnová L, Jocher G, Acosta M, Kowalska N, Kokrda L, Pavelka M, Marshall JD, et al. Assessing the Contrasting Effects of the Exceptional 2015 Drought on the Carbon Dynamics in Two Norway Spruce Forest Ecosystems. Atmosphere. 2021; 12(8):988. https://doi.org/10.3390/atmos12080988

Chicago/Turabian Style

Mensah, Caleb, Ladislav Šigut, Milan Fischer, Lenka Foltýnová, Georg Jocher, Manuel Acosta, Natalia Kowalska, Lukáš Kokrda, Marian Pavelka, John David Marshall, and et al. 2021. "Assessing the Contrasting Effects of the Exceptional 2015 Drought on the Carbon Dynamics in Two Norway Spruce Forest Ecosystems" Atmosphere 12, no. 8: 988. https://doi.org/10.3390/atmos12080988

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

Article Metrics

Back to TopTop