Understanding the low-temperature limitations to forest growth through calibration of a forest dynamics model with tree-ring data

The sensitivity of altitudinal and latitudinal tree-line ecotones to climate change, particularly that of temperature, has received much attention. To improve our understanding of the factors affecting tree-line position, we used the spatially explicit dynamic forest model TreeMig. Although well-suited because of its landscape dynamics functions, TreeMig features a parabolic temperature growth response curve, which has recently been questioned. and the species parameters are not specifically calibrated for cold temperatures. Our main goals were to improve the theoretical basis of the temperature growth response curve in the model and develop a method for deriving that curve's parameters from tree-ring data. We replaced the parabola with an asymptotic curve, calibrated for the main species at the subalpine (Swiss Alps: Pinus cembra, Larix decidua, Picea abies) and boreal (Fennoscandia: Pinus sylvestris, Betula pubescens, P. abies) tree-lines. After fitting new parameters, the growth curve matched observed tree-ring widths better. For the subalpine species, the minimum degree-day sum allowing, growth (kDDMin) was lowered by around 100 degree-days; in the case of Larix, the maximum potential ring-width was increased to 5.19 mm. At the boreal tree-line, the kDDMin for P. sylvestris was lowered by 210 degree-days and its maximum ring-width increased to 2.943 mm; for Betula (new in the model) kDDMin was set to 325 degree-days and the maximum ring-width to 2.51 mm; the values from the only boreal sample site for Picea were similar to the subalpine ones, so the same parameters were used. However, adjusting the growth response alone did not improve the model's output concerning species' distributions and their relative importance at tree-line. Minimum winter temperature (MinWiT, mean of the coldest winter month), which controls seedling establishment in TreeMig, proved more important for determining distribution. Picea, P. sylvestris and Betula did not previously have minimum winter temperature limits, so these values were set to the 95th percentile of each species' coldest MinWiT site (respectively -7, -11, -13). In a case study for the Alps, the original and newly calibrated versions of TreeMig were compared with biomass data from the National Forest Inventor), (NFI). Both models gave similar, reasonably realistic results. In conclusion, this method of deriving temperature responses from tree-rings works well. However, regeneration and its underlying factors seem more important for controlling species' distributions than previously thought. More research on regeneration ecology, especially at the upper limit of forests. is needed to improve predictions of tree-line responses to climate change further.

Modeling Forest Growth with Management Data: A Matrix Approach for the Italian Alps

Summary

Site- and species-specific responses of forest growth to climate across the European continent

Aim To evaluate the climate sensitivity of model-based forest productivity estimates using a continental-scale tree-ring network. Location Europe and North Africa (30–70° N, 10° W–40° E). Methods We compiled close to 1000 annually resolved records of radial tree growth for all major European tree species and quantified changes in growth as a function of historical climatic variation. Sites were grouped using a neural network clustering technique to isolate spatiotemporal and species-specific climate response patterns. The resulting empirical climate sensitivities were compared with the sensitivities of net primary production (NPP) estimates derived from the ORCHIDEE-FM and LPJ-wsl dynamic global vegetation models (DGVMs). Results We found coherent biogeographic patterns in climate response that depend upon (1) phylogenetic controls and (2) ambient environmental conditions delineated by latitudinal/elevational location. Temperature controls dominate forest productivity in high-elevation and high-latitude areas whereas moisture sensitive sites are widespread at low elevation in central and southern Europe. DGVM simulations broadly reproduce the empirical patterns, but show less temperature sensitivity in the boreal zone and stronger precipitation sensitivity towards the mid-latitudes. Main conclusions Large-scale forest productivity is driven by monthly to seasonal climate controls, but our results emphasize species-specific growth patterns under comparable environmental conditions. Furthermore, we demonstrate that carry-over effects from the previous growing season can significantly influence tree growth, particularly in areas with harsh climatic conditions – an element not considered in most current-state DGVMs. Model–data discrepancies suggest that the simulated climate sensitivity of NPP will need refinement before carbon-cycle climate feedbacks can be accurately quantified.

The influence of sampling design on tree-ring-based quantification of forest growth

Tree-rings offer one of the few possibilities to empirically quantify and reconstruct forest growth dynamics over years to millennia. Contemporaneously with the growing scientific community employing tree-ring parameters, recent research has suggested that commonly applied sampling designs (i.e. how and which trees are selected for dendrochronological sampling) may introduce considerable biases in quantifications of forest responses to environmental change. To date, a systematic assessment of the consequences of sampling design on dendroecological and-climatological conclusions has not yet been performed. Here, we investigate potential biases by sampling a large population of trees and replicating diverse sampling designs. This is achieved by retroactively subsetting the population and specifically testing for biases emerging for climate reconstruction, growth response to climate variability, long-term growth trends, and quantification of forest productivity. We find that commonly applied sampling designs can impart systematic biases of varying magnitude to any type of tree-ring-based investigations, independent of the total number of samples considered. Quantifications of forest growth and productivity are particularly susceptible to biases, whereas growth responses to short-term climate variability are less affected by the choice of sampling design. The world's most frequently applied sampling design, focusing on dominant trees only, can bias absolute growth rates by up to 459% and trends in excess of 200%. Our findings challenge paradigms, where a subset of samples is typically considered to be representative for the entire population. The only two sampling strategies meeting the requirements for all types of investigations are the (i) sampling of all individuals within a fixed area; and (ii) fully randomized selection of trees. This result advertises the consistent implementation of a widely applicable sampling design to simultaneously reduce uncertainties in tree-ring-based quantifications of forest growth and increase the comparability of datasets beyond individual studies, investigators, laboratories, and geographical boundaries.

Growth rate of a terra firme rain forest in Brazilian Amazonia over an eight-year period in response to logging

This paper deals with growth rates of trees > 5cm dbh over an eight-year period from 257 species at the Tapajós National Forest. The discussion is centred on the behaviour of the forest after logging. Permanent sample plots were established in 1981 and measured at the first time. The area was logged in 1982. Measurements after logging occurred in 1983, 1987 and 1989. Considering all species together, diameter increment was similar for both intensities of logging until five years after logging. Light-demanding species showed significantly higher growth rates than shade-tolerant species in the logged forest, with greater increment in the heavier treatment intensity. Commercial species also had higher growth rates in the heavier logged area, although those were significantly different only in the period from one to five years after logging. In the undisturbed forest, growth rates increased with increasing dbh size. At species level, growth rate varied between and within treatments, as well as between trees within species, depending mainly on degree of canopy opening. The logging favoured the growth of commercial species, chiefly the light-demanders. Therefore, if the same growth conditions continue being given, for example by silvicultural treatments, to those species of commercial interest, the forest would reach a stock available for harvesting around year 30 after logging. However, the high variation in increment rates indicates that an eight-year period is not sufficient to allow predictions on cutting cycles or polycyclic management systems for the study forest.

Biomass partitioning and growth efficiency in four naturally regenerated forest tree species

Current forest growth models and yield tables are almost exclusively based on data from mature trees, reducing their applicability to young and developing stands. To address this gap, young European beech, sessile oak, Scots pine and Norway spruce trees approximately 0 to 10 years old were destructively sampled in a range of naturally regenerated forest stands in Central Europe. Diameter at base and height were first measured in situ for up to 175 individuals per species. Subsequently, the trees were excavated and dry biomass of foliage, branches, stems and roots was measured. Allometric relations were then used to calculate biomass allocation coefficients (BAC) and growth efficiency (GE) patterns in young trees. We found large differences in BAC and GE between broadleaves and conifers, but also between species within these categories. Both BAC and GE are strongly age-specific in young trees, their rapidly changing values reflecting different growth strategies in the earliest stages of growth. We show that linear relationships describing biomass allocation in older trees are not applicable in young trees. To accurately predict forest biomass and carbon stocks, forest growth models need to include species and age specific parameters of biomass allocation patterns.

Response of nutrient cycles of an old-growth montane forest in Ecuador to experimental low-level nutrient amendments

Atmospheric nitrogen (N) and phosphorus (P) depositions are expected to increase in the tropicsrnas a consequence of increasing human activities in the next decades. Furthermore, a possiblernshortened El Niño Southern Oscillation cycle might come along with more frequent calcium (Ca)rndepositions on the eastern slope of the Ecuadorian Andes originating from Saharan dust. It isrncrucial to understand the response of the old-growth montane forest in Ecuador to increasedrnnutrient deposition to predict the further development of this megadiverse ecosystem.rnI studied experimental additions of N, P, N+P and Ca to the forest and an untreatedrncontrol, all in a fourfold replicated randomized block design. These experiments were conductedrnin the framework of a collaborative research effort, the NUtrient Manipulation EXperimentrn(NUMEX). I collected litter leachate, mineral soil solution (0.15 and 0.30 m depths), throughfallrnand fine litterfall samples and determined N, P and Ca concentrations and fluxes. This approachrnalso allowed me to assess whether N, P and/or Ca are limiting nutrients for forest growth.rnFurthermore, I evaluated the response of fine root biomass, leaf area index, leaf area and specificrnleaf area, tree diameter growth and basal area increment contributed from a cooperating group inrnthe Ca applied and control treatments.rnDuring the observation period of 16 months after the first fertilizer application, less thanrn10, 1 and 5% of the applied N, P and Ca, respectively, leached below the organic layer whichrncontained almost all roots but no significant leaching losses occurred to the deeper mineral soil.rnDeposited N, P and Ca from the atmosphere in dry and wet form were, on balance, retained in therncanopy in the control treatment. Retention of N, P and Ca in the canopy in their respectiverntreatments was reduced resulting in higher concentrations and fluxes of N, P and Ca inrnthroughfall and litterfall. Up to 2.5% of the applied N and 2% of the applied P and Ca werernrecycled to the soil with throughfall. Fluxes of N, P and Ca in throughfall+litterfall were higher inrnthe fertilized treatments than in the control; up to 20, 5 and 25% of the applied N, P and Ca,rnrespectively, were recycled to the soil with throughfall+litterfall.rnIn the Ca-applied plots, fine root biomass decreased significantly. Also the leaf area of thernfour most common tree species tended to decrease and the specific leaf area increasedrnsignificantly in Graffenrieda emarginata Triana, the most common tree species in the study area.rnThese changes are known plant responses to reduced nutrient stress. Reduced aluminium (Al)rntoxicity as an explanation of the Ca effect was unlikely, because of almost complete organocomplexationrnof Al and molar Ca:Al concentration ratios in solution above the toxicity threshold.rnThe results suggest that N, P and Ca co-limit the forest ecosystem functioning in thernnorthern Andean montane forests in line with recent assumptions in which different ecosystemrncompartments and even different phenological stages may show different nutrient limitationsrn(Kaspari et al. 2008). I conclude that (1) the expected elevated N and P deposition will bernretained in the ecosystem, at least in the short term and hence, quality of river water will not bernendangered and (2) increased Ca input will reduce nutrient stress of the forest.

Coexistence of forest and savanna in an Amazonian area from a geological perspective

Question: How can the coexistence of savanna and forest in Amazonian areas with relatively uniform climates be explained? Location: Eastern Marajo Island, northeast Amazonia, Brazil. Methods: The study integrated floristic analysis, terrain morphology, sedimentology and delta(13)C of soil organic matter. Floristic analysis involved rapid ecological assessment of 33 sites, determination of occurrence, specific richness, hierarchical distribution and matrix of floristic similarity between paired vegetation types. Terrain characterization was based on analysis of Landsat images using 4(R), 5(G) and 7(B) composition and digital elevation model (DEM). Sedimentology involved field descriptions of surface and core sediments. Finally, radiocarbon dating and analysis of delta(13)C of soil profile organic matter and natural ecotone forest-savanna was undertaken. Results: Slight tectonic subsidence in eastern Marajo Island favours seasonal flooding, making it unsuitable for forest growth. However, this area displays slightly convex-up, sinuous morphologies related to paleochannels, covered by forest. Terra-firme lowland forests are expanding from west to east, preferentially occupying paleochannels and replacing savanna. Slack, running water during channel abandonment leads to disappearance of varzea/gallery forest at channel margins. Long-abandoned channels sustain continuous terra-firme forests, because of longer times for more species to establish. Recently abandoned channels have had less time to become sites for widespread tree development, and are either not vegetated or covered by savanna. Conclusion: Landforms in eastern Marajo Island reflect changes in the physical environment due to reactivation of tectonic faults during the latest Quaternary. This promoted a dynamic history of channel abandonment, which controlled a set of interrelated parameters (soil type, topography, hydrology) that determined species location. Inclusion of a geological perspective for paleoenvironmental reconstruction can increase understanding of plant distribution in Amazonia.

Improved understanding of drought controls on seasonal variation in Mediterranean forest canopy CO2 and water fluxes through combined in situ measurements and ecosystem modelling

Water stress is a defining characteristic of Mediterranean ecosystems, and is likely to become more severe in the coming decades. Simulation models are key tools for making predictions, but our current understanding of how soil moisture controls ecosystem functioning is not sufficient to adequately constrain parameterisations. Canopy-scale flux data from four forest ecosystems with Mediterranean-type climates were used in order to analyse the physiological controls on carbon and water flues through the year. Significant non-stomatal limitations on photosynthesis were detected, along with lesser changes in the conductance-assimilation relationship. New model parameterisations were derived and implemented in two contrasting modelling approaches. The effectiveness of two models, one a dynamic global vegetation model ('ORCHIDEE'), and the other a forest growth model particularly developed for Mediterranean simulations ('GOTILWA+'), was assessed and modelled canopy responses to seasonal changes in soil moisture were analysed in comparison with in situ flux measurements. In contrast to commonly held assumptions, we find that changing the ratio of conductance to assimilation under natural, seasonally-developing, soil moisture stress is not sufficient to reproduce forest canopy CO2 and water fluxes. However, accurate predictions of both CO2 and water fluxes under all soil moisture levels encountered in the field are obtained if photosynthetic capacity is assumed to vary with soil moisture. This new parameterisation has important consequences for simulated responses of carbon and water fluxes to seasonal soil moisture stress, and should greatly improve our ability to anticipate future impacts of climate changes on the functioning of ecosystems in Mediterranean-type climates.

The impact of fragmentation and density regulation on forest succession in the Atlantic rain forest

The Atlantic Rain Forest, an important biodiversity hot spot, has faced severe habitat loss since the last century which has resulted in a highly fragmented landscape with a large number of small forest patches (<100 ha). For conservation planning it is essential to understand how current and future forest regeneration depends on ecological processes, fragment size and the connection to the regional seed pool. We have investigated the following questions by applying the forest growth simulation model FORMIND to the situation of the Atlantic Forest in the state of Sao Paulo, SE Brazil: (1) which set of parameters describing the local regeneration and level of density regulation can reproduce the biomass distribution and stem density of an old growth forest in a reserve? (2) Which additional processes apart from those describing the dynamics of an old growth forest, drive forest succession of small isolated fragments? (3) Which role does external seed input play during succession? Therefore, more than 300 tree species have been classified into nine plant functional types (PFTs), which are characterized by maximum potential height and shade tolerance. We differentiate between two seed dispersal modes: (i) local dispersal, i.e. all seedlings originated from fertile trees within the simulated area and (ii) external seed rain. Local seed dispersal has been parameterized following the pattern oriented approach, using biomass estimates of old growth forest. We have found that moderate density regulation is essential to achieve coexistence for a broad range of regeneration parameters. Considering the expected uncertainty and variability in the regeneration processes it is important that the forest dynamics are robust to variations in the regeneration parameters. Furthermore, edge effects such as increased mortality at the border and external seed rain have been necessary to reproduce the patterns for small isolated fragments. Overall, simulated biomass is much lower in the fragments compared to the continuous forest, whereas shade tolerant species are affected most strongly by fragmentation. Our simulations can supplement empirical studies by extrapolating local knowledge on edge effects of fragments to larger temporal and spatial scales. In particular our results show the importance of external seed rain and therefore highlight the importance of structural connectivity between regenerating fragments and mature forest stands. (C) 2009 Elsevier B.V. All rights reserved.

A new probabilistic canopy dynamics model (SLCD) that is suitable for evergreen and deciduous forest ecosystems

There are strong uncertainties regarding LAI dynamics in forest ecosystems in response to climate change. While empirical growth & yield models (G&YMs) provide good estimations of tree growth at the stand level on a yearly to decennial scale, process-based models (PBMs) use LAI dynamics as a key variable for enabling the accurate prediction of tree growth over short time scales. Bridging the gap between PBMs and G&YMs could improve the prediction of forest growth and, therefore, carbon, water and nutrient fluxes by combining modeling approaches at the stand level.Our study aimed to estimate monthly changes of leaf area in response to climate variations from sparse measurements of foliage area and biomass. A leaf population probabilistic model (SLCD) was designed to simulate foliage renewal. The leaf population was distributed in monthly cohorts, and the total population size was limited depending on forest age and productivity. Foliage dynamics were driven by a foliation function and the probabilities ruling leaf aging or fall. Their formulation depends on the forest environment.The model was applied to three tree species growing under contrasting climates and soil types. In tropical Brazilian evergreen broadleaf eucalypt plantations, the phenology was described using 8 parameters. A multi-objective evolutionary algorithm method (MOEA) was used to fit the model parameters on litterfall and LAI data over an entire stand rotation. Field measurements from a second eucalypt stand were used to validate the model. Seasonal LAI changes were accurately rendered for both sites (R-2 = 0.898 adjustment, R-2 = 0.698 validation). Litterfall production was correctly simulated (R-2 = 0.562, R-2 = 0.4018 validation) and may be improved by using additional validation data in future work. In two French temperate deciduous forests (beech and oak), we adapted phenological sub-modules of the CASTANEA model to simulate canopy dynamics, and SLCD was validated using LAI measurements. The phenological patterns were simulated with good accuracy in the two cases studied. However, IA/max was not accurately simulated in the beech forest, and further improvement is required.Our probabilistic approach is expected to contribute to improving predictions of LAI dynamics. The model formalism is general and suitable to broadleaf forests for a large range of ecological conditions. (C) 2014 Elsevier B.V. All rights reserved.

Modelling the effects of nitrogen deposition and carbon dioxide enrichment on forest carbon balance

Atmospheric CO2 concentration ([CO2]) has increased over the last 250 years, mainly due to human activities. Of total anthropogenic emissions, almost 31% has been sequestered by the terrestrial biosphere. A considerable contribution to this sink comes from temperate and boreal forest ecosystems of the northern hemisphere, which contain a large amount of carbon (C) stored as biomass and soil organic matter. Several potential drivers for this forest C sequestration have been proposed, including increasing atmospheric [CO2], temperature, nitrogen (N) deposition and changes in management practices. However, it is not known which of these drivers are most important. The overall aim of this thesis project was to develop a simple ecosystem model which explicitly incorporates our best understanding of the mechanisms by which these drivers affect forest C storage, and to use this model to investigate the sensitivity of the forest ecosystem to these drivers. I firstly developed a version of the Generic Decomposition and Yield (G’DAY) model to explicitly investigate the mechanisms leading to forest C sequestration following N deposition. Specifically, I modified the G’DAY model to include advances in understanding of C allocation, canopy N uptake, and leaf trait relationships. I also incorporated a simple forest management practice subroutine. Secondly, I investigated the effect of CO2 fertilization on forest productivity with relation to the soil N availability feedback. I modified the model to allow it to simulate short-term responses of deciduous forests to environmental drivers, and applied it to data from a large-scale forest Free-Air CO2 Enrichment (FACE) experiment. Finally, I used the model to investigate the combined effects of recent observed changes in atmospheric [CO2], N deposition, and climate on a European forest stand. The model developed in my thesis project was an effective tool for analysis of effects of environmental drivers on forest ecosystem C storage. Key results from model simulations include: (i) N availability has a major role in forest ecosystem C sequestration; (ii) atmospheric N deposition is an important driver of N availability on short and long time-scales; (iii) rising temperature increases C storage by enhancing soil N availability and (iv) increasing [CO2] significantly affects forest growth and C storage only when N availability is not limiting.

Placing unprecedented recent fir growth in a European-wide and Holocene-long context

Forest decline played a pivotal role in motivating Europe's political focus on sustainability around 35 years ago. Silver fir (Abies alba) exhibited a particularly severe dieback in the mid-1970s, but disentangling biotic from abiotic drivers remained challenging because both spatial and temporal data were lacking. Here, we analyze 14 136 samples from living trees and historical timbers, together with 356 pollen records, to evaluate recent fir growth from a continent-wide and Holocene-long perspective. Land use and climate change influenced forest growth over the past millennium, whereas anthropogenic emissions of acidic sulfates and nitrates became important after about 1850. Pollution control since the 1980s, together with a warmer but not drier climate, has facilitated an unprecedented surge in productivity across Central European fir stands. Restricted fir distribution prior to the Mesolithic and again in the Modern Era, separated by a peak in abundance during the Bronze Age, is indicative of the long-term interplay of changing temperatures, shifts in the hydrological cycle, and human impacts that have shaped forest structure and productivity.

Unusual Southern Hemisphere tree growth patterns induced by changes in the Southern Annular Mode

Recent changes in the summer climate of the Southern Hemisphere extra-tropics are primarily related to the dominance of the positive phase of the Southern Annular Mode1, 2. This shift in the behaviour of the Southern Annular Mode—essentially a measure of the pressure gradient between Southern Hemisphere mid and high latitudes—has been predominantly induced by polar stratospheric ozone depletion2, 3, 4. The concomitant southward expansion of the dry subtropical belts5, 6 could have consequences for forest growth. Here, we use tree-ring records from over 3,000 trees in South America, Tasmania and New Zealand to identify dominant patterns of tree growth in recent centuries. We show that the foremost patterns of growth between 1950 and 2000 differed significantly from those in the previous 250 years. Specifically, growth was higher than the long-term average in the subalpine forests of Tasmania and New Zealand, but lower in the dry-mesic forests of Patagonia. We further demonstrate that variations in the Southern Annular Mode can explain 12–48% of the tree growth anomalies in the latter half of the twentieth century. Tree-ring-based reconstructions of summer Southern Annular Mode indices suggest that the high frequency of the positive phase since the 1950s is unprecedented in the past 600 years. We propose that changes in the Southern Annular Mode have significantly altered tree growth patterns in the Southern Hemisphere.

An empirical mixed model to quantify climate influence on the growth of Pinus halepensis Mill. stands in South-Eastern Spain

The influence of climate on forest stand composition, development and growth is undeniable. Many studies have tried to quantify the effect of climatic variables on forest growth and yield. These works become especially important because there is a need to predict the effects of climate change on the development of forest ecosystems. One of the ways of facing this problem is the inclusion of climatic variables into the classic empirical growth models. The work has a double objective: (i) to identify the indicators which best describe the effect of climate on Pinus halepensis growth and (ii) to quantify such effect in several scenarios of rainfall decrease which are likely to occur in the Mediterranean area. A growth mixed model for P. halepensis including climatic variables is presented in this work. Growth estimates are based on data from the Spanish National Forest Inventory (SNFI). The best results are obtained for the indices including rainfall, or rainfall and temperature together, with annual precipitation, precipitation effectiveness, Emberger?s index or free bioclimatic intensity standing out among them. The final model includes Emberger?s index, free bioclimatic intensity and interactions between competition and climate indices. The results obtained show that a rainfall decrease about 5% leads to a decrease in volume growth of 5.5?7.5% depending on site quality.