•Sensitivity analysis of LWF-Brook90 revealed a distinct parameter importance ranking•Bayesian calibration indicated high water use for a poplar short rotation forest•The LWF-Brook90 SVAT model ...excellently captures temporal variations of soil moisture•The LWFBrook90R package facilitates complex statistical analysis and parallelization•The modelling case study with code examples demonstrates the utility of the R package
Soil vegetation atmosphere transport (SVAT) models are important for the quantification of water fluxes, soil water availability, drought stress and their uncertainties under climate change. We present LWFBrook90R, an enhanced implementation of the well-established, process-based SVAT model LWF-Brook90 for the R environment for statistical computing. The package provides new functions and sub-models for model parameterization, and facilitates parallel computing, sensitivity analysis and inverse calibration of the model. A case study comprising i) basic forward water balance simulations for temperate grassland vegetation, deciduous and evergreen forest, ii) a parallelized sensitivity analysis, and iii) Bayesian calibrations based on soil water storage observed in a poplar (Populus nigra × P. maximowiczii) Short Rotation Forest (SRF) demonstrates the utility of the R package. The sensitivity analysis revealed parameters affecting plant-available soil water storage capacity and the vegetation's timing and level of water demand to be most important for the annual course of simulated soil water storage, with seasonal and interannual differences in parameter importance rankings. The subsequent calibration yielded a very high agreement between daily simulated and observed soil water storage (0-200 cm soil depth) for the calibration and validation datasets, with Nash-Sutcliffe efficiencies of 0.97 and 0.95, respectively. The final model predicted high though realistic rates of annual evapotranspiration (2011: 844 ± 3.8 mm y-1, 2012: 733 ± 4.5 mm y-1) for the poplar SRF, regularly exceeding grass reference evaporation (ET0) by 20-47% during the months of the growing season. However, basing calibrations solely on observed soil water storage probably resulted in biased partitioning of evapotranspiration towards interception losses. The integration of the LWF-Brook90 hydrological model into R with its wide variety of extensions was successfully tested and may provide efficient, reliable and reproducible water balance predictions by facilitating complex statistical analyses and large-scale applications of the model.
► We analysed the age structure of a virgin beech forest based on 164 increment cores. ► Beech can reach ages of up to 500
years, although old trees are prone to stem rot. ► Beech can survive long ...suppression periods of more than 100
years. ► Small-scale disturbances lead to continuous tree recruitment and structural diversity.
The Carpathian mountains harbour the largest areas of virgin European beech (
Fagus sylvatica L
.) forest. Understanding the growth dynamics of European beech under natural conditions without human intervention is crucial for close-to-nature management of beech forests in Europe. In this study we give an insight into the natural disturbance dynamics of the virgin beech forest Uholka, based on a structural analysis and a dendroecological reconstruction of its history. On four circular plots of 0.1
ha each, DBH and tree height of all living trees ⩾6
cm DBH were measured. Increment cores of all the 164 trees were taken for age estimation and growth pattern analysis. To identify significant growth releases and reconstruct the disturbance history a dendroecological approach, referred to as the boundary line method, was employed. The density of the living trees per plot ranged from 270 to 590 stems per ha and the volume from 525 to 1237
m
3 per ha. The longest tree-ring series was 451
years long, however age estimation methods showed that beech could reach an age of up to 550
years. All four plots covered an age span of at least 300
years and can be characterised as uneven-aged with continuous tree establishment. Beech can survive long suppression periods of over 100
years and shows regularly distributed growth releases over the analysed period (1870–1999). The results suggest that stand dynamics in this forest are driven by periodic small disturbances, while larger events occur only rarely. The high percentage of rotten trees in the upper canopy indicates that individual trees are prone to windbreakage, which promotes these small-scale dynamics.
Determining the drivers of shifting forest disturbance rates remains a pressing global change issue. Large‐scale forest dynamics are commonly assumed to be climate driven, but appropriately scaled ...disturbance histories are rarely available to assess how disturbance legacies alter subsequent disturbance rates and the climate sensitivity of disturbance. We compiled multiple tree ring‐based disturbance histories from primary Picea abies forest fragments distributed throughout five European landscapes spanning the Bohemian Forest and the Carpathian Mountains. The regional chronology includes 11,595 tree cores, with ring dates spanning the years 1750–2000, collected from 560 inventory plots in 37 stands distributed across a 1,000 km geographic gradient, amounting to the largest disturbance chronology yet constructed in Europe. Decadal disturbance rates varied significantly through time and declined after 1920, resulting in widespread increases in canopy tree age. Approximately 75% of current canopy area recruited prior to 1900. Long‐term disturbance patterns were compared to an historical drought reconstruction, and further linked to spatial variation in stand structure and contemporary disturbance patterns derived from LANDSAT imagery. Historically, decadal Palmer drought severity index minima corresponded to higher rates of canopy removal. The severity of contemporary disturbances increased with each stand's estimated time since last major disturbance, increased with mean diameter, and declined with increasing within‐stand structural variability. Reconstructed spatial patterns suggest that high small‐scale structural variability has historically acted to reduce large‐scale susceptibility and climate sensitivity of disturbance. Reduced disturbance rates since 1920, a potential legacy of high 19th century disturbance rates, have contributed to a recent region‐wide increase in disturbance susceptibility. Increasingly common high‐severity disturbances throughout primary Picea forests of Central Europe should be reinterpreted in light of both legacy effects (resulting in increased susceptibility) and climate change (resulting in increased exposure to extreme events).
Climate change has been linked to increasing forest disturbance rates and large‐scale disturbance histories are needed to assess how disturbance legacies modulate climate–disturbance dynamics. Historical disturbance rates (c. 1750–2000), reconstructed with tree cores collected from primary Picea abies forests throughout Central‐Eastern Europe, have exhibited moderate sensitivity to drought extremes; developmental feedbacks are important determinants of a forest's responsiveness to climate‐induced events. Widespread increases in forest age have occurred throughout the 20th century and these developmental patterns, in addition to climate change, are responsible for the severity of recent disturbances.
•Radial-growth averaging and absolute-increase had fewer type I and II errors.•Radial-growth averaging require the least a priori data.•Time-series produces information on magnitude and duration of ...release events.•False positives were more common in forests with infrequent and low severity disturbance.•Yearly binning performs better to identify local maxima indicating disturbance.
The retrospective study of abrupt and sustained increases in the radial growth of trees (hereinafter ‘releases’) by tree-ring analysis is an approach widely used for reconstructing past forest disturbances. Despite the range of dendrochronological methods used for release-detection, a lack of in-depth comparison between them can lead researchers to question which method to use and, potentially, increases the uncertainties of disturbance histories derived with different methods.
Here, we investigate the efficacy and sensitivity of four widely used release detection methods using tree-ring width series and complete long-term inventories of forest stands with known disturbances. We used support vector machine (SVM) analysis trained on long-term forest census data to estimate the likelihood that Acer rubrum trees experiencing reductions in competition show releases in their tree-ring widths. We compare methods performance at the tree and stand level, followed by evaluation of method sensitivity to changes in their parameters and settings.
Disturbance detection methods agreed with 60–76% of the SVM-identified growth releases under high canopy disturbance and 80–94% in a forest with canopy disturbance of low severity and frequency. The median competition index change (CIC) of trees identified as being released differed more than two-fold between methods, from −0.33 (radial-growth averaging) to −0.68 (time-series). False positives (type I error) were more common in forests with low severity disturbance, whereas false negatives (type II error) occurred more often in forests with high severity disturbance. Sensitivity analysis indicated that reductions of the detection threshold and the length of the time window significantly increased detected stand-level disturbance severity across all methods.
Radial-growth averaging and absolute-increase methods had lower levels of type I and II error in detecting disturbance events with our datasets. Parameter settings play a key role in the accuracy of reconstructing disturbance history regardless of the method. Time-series and radial-growth averaging methods require the least amount of a priori information, but only the time-series method quantified the subsequent growth increment related to a reduction in competition. Finally, we recommend yearly binning of releases using a kernel density estimation function to identify local maxima indicating disturbance. Kernel density estimation improves reconstructions of forest history and, thus, will further our understanding of past forest dynamics.
The growth of past, present, and future forests was, is and will be affected by climate variability. This multifaceted relationship has been assessed in several regional studies, but spatially ...resolved, large-scale analyses are largely missing so far. Here we estimate recent changes in growth of 5800 beech trees (Fagus sylvatica L.) from 324 sites, representing the full geographic and climatic range of species. Future growth trends were predicted considering state-of-the-art climate scenarios. The validated models indicate growth declines across large region of the distribution in recent decades, and project severe future growth declines ranging from -20% to more than -50% by 2090, depending on the region and climate change scenario (i.e. CMIP6 SSP1-2.6 and SSP5-8.5). Forecasted forest productivity losses are most striking towards the southern distribution limit of Fagus sylvatica, in regions where persisting atmospheric high-pressure systems are expected to increase drought severity. The projected 21
century growth changes across Europe indicate serious ecological and economic consequences that require immediate forest adaptation.
Process-based forest models are important tools for predicting forest growth and their vulnerability to factors such as climate change or responses to management. One of the most widely used ...stand-level process-based models is the 3-PG model (Physiological Processes Predicting Growth), which is used for applications including estimating wood production, carbon budgets, water balance and susceptibility to climate change. Few 3-PG parameter sets are available for central European species and even fewer are appropriate for mixed-species forests. Here we estimated 3-PG parameters for twelve major central European tree species using 1418 long-term permanent forest monitoring plots from managed forests, 297 from un-managed forest reserves and 784 Swiss National Forest Inventory plots. A literature review of tree physiological characteristics, as well as regression analyses and Bayesian inference, were used to calculate the 3-PG parameters.
The Swiss-wide calibration, based on monospecific plots, showed a robust performance in predicting forest stocks such as stem, foliage and root biomass. The plots used to inform the Bayesian calibration resulted in posterior ranges of the calibrated parameters that were, on average, 69% of the prior range. The bias of stem, foliage and root biomass predictions was generally less than 20%, and less than 10% for several species. The parameter sets also provided reliable predictions of biomass and mean tree sizes in mixed-species forests. Given that the information sources used to develop the parameters included a wide range of climatic, edaphic and management conditions and long time spans (from 1930 to present), these species parameters for 3-PG are likely to be appropriate for most central European forests and conditions.
The expected future intensification of forest disturbance as a consequence of ongoing anthropogenic climate change highlights the urgent need to more robustly quantify associated biotic responses. ...Saproxylic beetles are a diverse group of forest invertebrates representing a major component of biodiversity that is associated with the decomposition and cycling of wood nutrients and carbon in forest ecosystems. Disturbance-induced declines or shifts in their diversity indicate the loss of key ecological and/or morphological species traits that could change ecosystem functioning. Functional and phylogenetic diversity of biological communities is commonly used to link species communities to ecosystem functions. However, our knowledge on how disturbance intensity alters functional and phylogenetic diversity of saproxylic beetles is incomplete. Here, we analyzed the main drivers of saproxylic beetle abundance and diversity using a comprehensive dataset from montane primary forests in Europe. We investigated cascading relationships between 250 years of historical disturbance mechanisms, forest structural attributes and the taxonomic, phylogenetic and functional diversity of present-day beetle communities. Our analyses revealed that historical disturbances have significant effects on current beetle communities. Contrary to our expectations, different aspects of beetle communities, that is, abundance, taxonomic, phylogenetic and functional diversity, responded to different disturbance regime components. Past disturbance frequency was the most important component influencing saproxylic beetle communities and habitat via multiple temporal and spatial pathways. The quantity of deadwood and its diameter positively influenced saproxylic beetle abundance and functional diversity, whereas phylogenetic diversity was positively influenced by canopy openness. Analyzing historical disturbances, we observed that current beetle diversity is far from static, such that the importance of various drivers might change during further successional development. Only forest landscapes that are large enough to allow for the full range of temporal and spatial patterns of disturbances and post-disturbance development will enable long-term species coexistence and their associated ecosystem functions.
Accurately capturing medium- to low-frequency trends in tree-ring data is vital to assessing climatic response and developing robust reconstructions of past climate. Non-climatic disturbance can ...affect growth trends in tree-ring-width (RW) series and bias climate information obtained from such records. It is important to develop suitable strategies to ensure the development of chronologies that minimize these medium- to low-frequency biases. By performing high density sampling (760 trees) over a ~40-ha natural high-elevation Norway spruce (Picea abies) stand in the Romanian Carpathians, this study assessed the suitability of several sampling strategies for developing chronologies with an optimal climate signal for dendroclimatic purposes. There was a roughly equal probability for chronologies (40 samples each) to express a reasonable (r = 0.3–0.5) to non-existent climate signal. While showing a strong high-frequency response, older/larger trees expressed the weakest overall temperature signal. Although random sampling yielded the most consistent climate signal in all sub-chronologies, the outcome was still sub-optimal. Alternative strategies to optimize the climate signal, including very high replication and principal components analysis, were also unable to minimize this disturbance bias and produce chronologies adequately representing climatic trends, indicating that larger scale disturbances can produce synchronous pervasive disturbance trends that affect a large part of a sampled population. The Curve Intervention Detection (CID) method, used to identify and reduce the influence of disturbance trends in the RW chronologies, considerably improved climate signal representation (from r = 0.28 before correction to r = 0.41 after correction for the full 760 sample chronology over 1909–2009) and represents a potentially important new approach for assessing disturbance impacts on RW chronologies. Blue intensity (BI) also shows promise as a climatically more sensitive variable which, unlike RW, does not appear significantly affected by disturbance. We recommend that studies utilizing RW chronologies to investigate medium- to long-term climatic trends also assess disturbance impact on those series.
Display omitted
•High summer ozone and foliar symptoms were frequent in European forests.•They had stagnating, slightly declining time trends and different spatial pattern.•Ozone was higher in the ...Mediterranean and the Alpine biogeographic regions.•Symptoms tend to be driven by functional leaf traits.•Ozone had a significant effect on symptoms on the most sensitive species.
Tropospheric ozone (O3) increased globally in the 20th century, contributes to climate change and can have adverse effects on terrestrial ecosystems. The response of forest vegetation to ozone is modulated by species- and site-specific factors and visible foliar symptoms (VFS) are the only direct evidence of ozone effects on vegetation. VFS have been observed and reproduced under (semi-) controlled conditions and their field assessment has been largely harmonized in Europe. We analyzed ozone concentration and VFS data as measured at (respectively) 118 and 91 intensive monitoring sites of the International Co-Operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests (ICP Forests) spanning over five European biogeographic regions from 2005 to 2018. Average values for VFS were calculated accounting for the number of species present and their observed frequency. Spatial and temporal variation of ozone concentrations, VFS, and their relationships across Europe were then investigated by applying Generalized Linear Mixed Models (GLMMs) and combined GLMMs. Ozone concentrations exceeded 40 ppb on 37.3 % of the sites and were significantly higher (p < 0.05) in the Alpine and the Mediterranean regions. Over the 2005–2018 period there was a substantial stagnation of ozone concentrations with a tendency towards decreasing values in the Alpine-Boreal sites and increasing values in the Atlantic sites. Ozone left a “fingerprint” in terms of VFS on 38 % of the observed broadleaved woody species across Europe, with no significant difference among biogeographic regions. Overall, and again with the exception of an increase at the Atlantic sites, the frequency of VFS remained unchanged or has been slightly declining over the investigated period. We found positive relationship between ozone concentrations and VFS across Europe (p < 0.05), while their temporal trends (both insignificant) were not related. The species with the highest frequency of VFS were those classified as sensitive species under controlled/semi-controlled experimental conditions. Frequency of VFS tends to be modulated by vegetation traits such as specific leaf area and leaf thickness (p < 0.10). Our results showed that, although ozone levels suggested a North-to-South gradient of increasing potential risk to vegetation with hot spots in the Alps and in the Mediterranean, VFS observed on the actual species assemblage at the sites modifies this picture. According to frequency of VFS, ozone risk for vegetation may be higher in parts of the Alpine and Continental Europe than in the Mediterranean region.
The TreeNet research and monitoring network has been continuously collecting data from point dendrometers and air and soil microclimate using an automated system since 2011. The goal of TreeNet is to ...generate high temporal resolution datasets of tree growth and tree water dynamics for research and to provide near real-time indicators of forest growth performance and drought stress to a wide audience. This paper explains the key working steps from the installation of sensors in the field to data acquisition, data transmission, data processing, and online visualization. Moreover, we discuss the underlying premises to convert dynamic stem size changes into relevant biological information. Every 10 min, the stem radii of about 420 trees from 13 species at 61 sites in Switzerland are measured electronically with micrometer precision, in parallel with the environmental conditions above and below ground. The data are automatically transmitted, processed and stored on a central server. Automated data processing (R-based functions) includes screening of outliers, interpolation of data gaps, and extraction of radial stem growth and water deficit for each tree. These long-term data are used for scientific investigations as well as to calculate and display daily indicators of growth trends and drought levels in Switzerland based on historical and current data. The current collection of over 100 million data points forms the basis for identifying dynamics of tree-, site- and species-specific processes along environmental gradients. TreeNet is one of the few forest networks capable of tracking the diurnal and seasonal cycles of tree physiology in near real-time, covering a wide range of temperate forest species and their respective environmental conditions.