Limits to reproduction and seed size-number trade-offs that shape forest dominance and future recovery.

The relationships that control seed production in trees are fundamental to understanding the evolution of forest species and their capacity to recover from increasing losses to drought, fire, and harvest. A synthesis of fecundity data from 714 species worldwide allowed us to examine hypotheses that are central to quantifying reproduction, a foundation for assessing fitness in forest trees. Four major findings emerged. First, seed production is not constrained by a strict trade-off between seed size and numbers. Instead, seed numbers vary over ten orders of magnitude, with species that invest in large seeds producing more seeds than expected from the 1:1 trade-off. Second, gymnosperms have lower seed production than angiosperms, potentially due to their extra investments in protective woody cones. Third, nutrient-demanding species, indicated by high foliar phosphorus concentrations, have low seed production. Finally, sensitivity of individual species to soil fertility varies widely, limiting the response of community seed production to fertility gradients. In combination, these findings can inform models of forest response that need to incorporate reproductive potential.

Globally, tree fecundity exceeds productivity gradients.

Lack of tree fecundity data across climatic gradients precludes the analysis of how seed supply contributes to global variation in forest regeneration and biotic interactions responsible for biodiversity. A global synthesis of raw seedproduction data shows a 250-fold increase in seed abundance from cold-dry to warm-wet climates, driven primarily by a 100-fold increase in seed production for a given tree size. The modest (threefold) increase in forest productivity across the same climate gradient cannot explain the magnitudes of these trends. The increase in seeds per tree can arise from adaptive evolution driven by intense species interactions or from the direct effects of a warm, moist climate on tree fecundity. Either way, the massive differences in seed supply ramify through food webs potentially explaining a disproportionate role for species interactions in the wet tropics.

Crown defoliation decreases reproduction and wood growth in a marginal European beech population.

BACKGROUND AND AIMS: Abiotic and biotic stresses related to climate change have been associated with increased crown defoliation, decreased growth and a higher risk of mortality in many forest tree species, but the impact of stresses on tree reproduction and forest regeneration remains understudied. At the dry, warm margin of species distributions, flowering, pollination and seed maturation are expected to be affected by drought, late frost and other stresses, eventually resulting in reproduction failure. Moreover, inter-individual variation in reproductive performance versus other performance traits (growth, survival) could have important consequences for population dynamics. This study investigated the relationships among individual crown defoliation, growth and reproduction in a drought-prone population of European beech, Fagus sylvatica. METHODS: We used a spatially explicit mating model and marker-based parentage analyses to estimate effective female and male fecundities of 432 reproductive trees, which were also monitored for basal area increment and crown defoliation over 9 years. KEY RESULTS: Female and male fecundities varied markedly between individuals, more than did growth. Both female fecundity and growth decreased with increasing crown defoliation and competition, and increased with size. Moreover, the negative effect of defoliation on female fecundity was size-dependent, with a slower decline in female fecundity with increasing defoliation for the large individuals. Finally, a trade-off between growth and female fecundity was observed in response to defoliation: some large trees maintained significant female fecundity at the expense of reduced growth in response to defoliation, while some other defoliated trees maintained high growth at the expense of reduced female fecundity. CONCLUSIONS: Our results suggest that, while decreasing their growth, some large defoliated trees still contribute to reproduction through seed production and pollination. This non-coordinated decline of growth and fecundity at individual level in response to stress may compromise the evolution of stress-resistance traits at population level, and increase forest tree vulnerability.

Where is the optimum? Predicting the variation of selection along climatic gradients and the adaptive value of plasticity. A case study on tree phenology.

Many theoretical models predict when genetic evolution and phenotypic plasticity allow adaptation to changing environmental conditions. These models generally assume stabilizing selection around some optimal phenotype. We however often ignore how optimal phenotypes change with the environment, which limit our understanding of the adaptive value of phenotypic plasticity. Here, we propose an approach based on our knowledge of the causal relationships between climate, adaptive traits, and fitness to further these questions. This approach relies on a sensitivity analysis of the process-based model phenofit, which mathematically formalizes these causal relationships, to predict fitness landscapes and optimal budburst dates along elevation gradients in three major European tree species. Variation in the overall shape of the fitness landscape and resulting directional selection gradients were found to be mainly driven by temperature variation. The optimal budburst date was delayed with elevation, while the range of dates allowing high fitness narrowed and the maximal fitness at the optimum decreased. We also found that the plasticity of the budburst date should allow tracking the spatial variation in the optimal date, but with variable mismatch depending on the species, ranging from negligible mismatch in fir, moderate in beech, to large in oak. Phenotypic plasticity would therefore be more adaptive in fir and beech than in oak. In all species, we predicted stronger directional selection for earlier budburst date at higher elevation. The weak selection on budburst date in fir should result in the evolution of negligible genetic divergence, while beech and oak would evolve counter-gradient variation, where genetic and environmental effects are in opposite directions. Our study suggests that theoretical models should consider how whole fitness landscapes change with the environment. The approach introduced here has the potential to be developed for other traits and species to explore how populations will adapt to climate change.

Tree xylem water isotope analysis by Isotope Ratio Mass Spectrometry and laser spectrometry: A dataset to explore tree response to drought.

Water isotopes from plant xylem and surrounding environment are increasingly used in eco-hydrological studies. Carrière et al. [1] analyzed a dataset of water isotopes in (i) the xylem of three different tree species, (ii) the surrounding soil and drainage water and (iii) the underlying karst groundwater, to understand tree water uptake during drought in two different Mediterranean forests on karst setting. The xylem and soil water were extracted by cryogenic distillation. The full dataset was obtained with Isotope Ratio Mass Spectrometry (IRMS) and Isotope Ratio Infrared Spectrometer (IRIS), and included 219 measurements of δ2H and δ18O. Prompted by unexpected isotopic data characterized by a very negative deuterium excess, a subsample of 46 xylem samples and 9 soil water samples were double checked with both analytical techniques. IRMS and IRIS analyses yielded similar data. Therefore, the results reveal that laser spectrometry allows an accurate estimation of xylem and soil water isotopes. The dataset highlights a strong 2H depletion in xylem water for all species. Deuterium does not seem adequate to interpret ecological processes in this dataset given the strong fractionation.

Microgeographic adaptation and the effect of pollen flow on the adaptive potential of a temperate tree species.

In species with long-distance dispersal capacities and inhabiting a large ecological niche, local selection and gene flow are expected to be major evolutionary forces affecting the genetic adaptation of natural populations. Yet, in species such as trees, evidence of microgeographic adaptation and the quantitative assessment of the impact of gene flow on adaptive genetic variation are still limited. Here, we used extensive genetic and phenotypic data from European beech seedlings collected along an elevation gradient, and grown in a common garden, to study the signature of selection on the divergence of eleven potentially adaptive traits, and to assess the role of gene flow in resupplying adaptive genetic variation. We found a significant signal of adaptive differentiation among plots separated by < 1 km, with selection acting on growth and phenological traits. Consistent with theoretical expectations, our results suggest that pollen dispersal contributes to increase genetic diversity for these locally differentiated traits. Our results thus highlight that local selection is an important evolutionary force in natural tree populations and suggest that management interventions to facilitate movement of gametes along short ecological gradients would boost genetic diversity of individual tree populations, and enhance their adaptive potential to rapidly changing environments.

The role of deep vadose zone water in tree transpiration during drought periods in karst settings - Insights from isotopic tracing and leaf water potential.

Karst environments are unusual because their dry, stony and shallow soils seem to be unfavorable to vegetation, and yet they are often covered with forests. How can trees survive in these environments? Where do they find the water that allows them to survive? This study uses midday and predawn water potentials and xylem water isotopes of branches to assess tree water status and the origin of transpired water. Monitoring was conducted during the summers of 2014 and 2015 in two dissimilar plots of Mediterranean forest located in karst environments. The results show that the three monitored tree species (Abies alba Mill, Fagus sylvatica L, and Quercus ilex L.) use deep water resources present in the karst vadose zone (unsaturated zone) more intensively during drier years. Quercus ilex, a species well- adapted to water stress, which grows at the drier site, uses the deep water resource very early in the summer season. Conversely, the two other species exploit the deep water resource only during severe drought. These results open up new perspectives to a better understanding of ecohydrological equilibrium and to improved water balance modeling in karst forest settings.

Research frontiers for improving our understanding of drought-induced tree and forest mortality.

Accumulating evidence highlights increased mortality risks for trees during severe drought, particularly under warmer temperatures and increasing vapour pressure deficit (VPD). Resulting forest die-off events have severe consequences for ecosystem services, biophysical and biogeochemical land-atmosphere processes. Despite advances in monitoring, modelling and experimental studies of the causes and consequences of tree death from individual tree to ecosystem and global scale, a general mechanistic understanding and realistic predictions of drought mortality under future climate conditions are still lacking. We update a global tree mortality map and present a roadmap to a more holistic understanding of forest mortality across scales. We highlight priority research frontiers that promote: (1) new avenues for research on key tree ecophysiological responses to drought; (2) scaling from the tree/plot level to the ecosystem and region; (3) improvements of mortality risk predictions based on both empirical and mechanistic insights; and (4) a global monitoring network of forest mortality. In light of recent and anticipated large forest die-off events such a research agenda is timely and needed to achieve scientific understanding for realistic predictions of drought-induced tree mortality. The implementation of a sustainable network will require support by stakeholders and political authorities at the international level.

Early-Warning Signals of Individual Tree Mortality Based on Annual Radial Growth.

Tree mortality is a key driver of forest dynamics and its occurrence is projected to increase in the future due to climate change. Despite recent advances in our understanding of the physiological mechanisms leading to death, we still lack robust indicators of mortality risk that could be applied at the individual tree scale. Here, we build on a previous contribution exploring the differences in growth level between trees that died and survived a given mortality event to assess whether changes in temporal autocorrelation, variance, and synchrony in time-series of annual radial growth data can be used as early warning signals of mortality risk. Taking advantage of a unique global ring-width database of 3065 dead trees and 4389 living trees growing together at 198 sites (belonging to 36 gymnosperm and angiosperm species), we analyzed temporal changes in autocorrelation, variance, and synchrony before tree death (diachronic analysis), and also compared these metrics between trees that died and trees that survived a given mortality event (synchronic analysis). Changes in autocorrelation were a poor indicator of mortality risk. However, we found a gradual increase in inter-annual growth variability and a decrease in growth synchrony in the last ∼20 years before mortality of gymnosperms, irrespective of the cause of mortality. These changes could be associated with drought-induced alterations in carbon economy and allocation patterns. In angiosperms, we did not find any consistent changes in any metric. Such lack of any signal might be explained by the relatively high capacity of angiosperms to recover after a stress-induced growth decline. Our analysis provides a robust method for estimating early-warning signals of tree mortality based on annual growth data. In addition to the frequently reported decrease in growth rates, an increase in inter-annual growth variability and a decrease in growth synchrony may be powerful predictors of gymnosperm mortality risk, but not necessarily so for angiosperms.

A synthesis of radial growth patterns preceding tree mortality.

Tree mortality is a key factor influencing forest functions and dynamics, but our understanding of the mechanisms leading to mortality and the associated changes in tree growth rates are still limited. We compiled a new pan-continental tree-ring width database from sites where both dead and living trees were sampled (2970 dead and 4224 living trees from 190 sites, including 36 species), and compared early and recent growth rates between trees that died and those that survived a given mortality event. We observed a decrease in radial growth before death in ca. 84% of the mortality events. The extent and duration of these reductions were highly variable (1-100 years in 96% of events) due to the complex interactions among study species and the source(s) of mortality. Strong and long-lasting declines were found for gymnosperms, shade- and drought-tolerant species, and trees that died from competition. Angiosperms and trees that died due to biotic attacks (especially bark-beetles) typically showed relatively small and short-term growth reductions. Our analysis did not highlight any universal trade-off between early growth and tree longevity within a species, although this result may also reflect high variability in sampling design among sites. The intersite and interspecific variability in growth patterns before mortality provides valuable information on the nature of the mortality process, which is consistent with our understanding of the physiological mechanisms leading to mortality. Abrupt changes in growth immediately before death can be associated with generalized hydraulic failure and/or bark-beetle attack, while long-term decrease in growth may be associated with a gradual decline in hydraulic performance coupled with depletion in carbon reserves. Our results imply that growth-based mortality algorithms may be a powerful tool for predicting gymnosperm mortality induced by chronic stress, but not necessarily so for angiosperms and in case of intense drought or bark-beetle outbreaks.

Towards a common methodology for developing logistic tree mortality models based on ring-width data.

Tree mortality is a key process shaping forest dynamics. Thus, there is a growing need for indicators of the likelihood of tree death. During the last decades, an increasing number of tree-ring based studies have aimed to derive growth-mortality functions, mostly using logistic models. The results of these studies, however, are difficult to compare and synthesize due to the diversity of approaches used for the sampling strategy (number and characteristics of alive and death observations), the type of explanatory growth variables included (level, trend, etc.), and the length of the time window (number of years preceding the alive/death observation) that maximized the discrimination ability of each growth variable. We assess the implications of key methodological decisions when developing tree-ring based growth-mortality relationships using logistic mixed-effects regression models. As examples, we use published tree-ring datasets from Abies alba (13 different sites), Nothofagus dombeyi (one site), and Quercus petraea (one site). Our approach is based on a constant sampling size and aims at (1) assessing the dependency of growth-mortality relationships on the statistical sampling scheme used, (2) determining the type of explanatory growth variables that should be considered, and (3) identifying the best length of the time window used to calculate them. The performance of tree-ring-based mortality models was reasonably high for all three species (area under the receiving operator characteristics curve, AUC > 0.7). Growth level variables were the most important predictors of mortality probability for two species (A. alba, N. dombeyi), while growth-trend variables need to be considered for Q. petraea. In addition, the length of the time window used to calculate each growth variable was highly uncertain and depended on the sampling scheme, as some growth-mortality relationships varied with tree age. The present study accounts for the main sampling-related biases to determine reliable species-specific growth-mortality relationships. Our results highlight the importance of using a sampling strategy that is consistent with the research question. Moving towards a common methodology for developing reliable growth-mortality relationships is an important step towards improving our understanding of tree mortality across species and its representation in dynamic vegetation models.

Simulating local adaptation to climate of forest trees with a Physio-Demo-Genetics model.

One challenge of evolutionary ecology is to predict the rate and mechanisms of population adaptation to environmental variations. The variations in most life history traits are shaped both by individual genotypic and by environmental variation. Forest trees exhibit high levels of genetic diversity, large population sizes, and gene flow, and they also show a high level of plasticity for life history traits. We developed a new Physio-Demo-Genetics model (denoted PDG) coupling (i) a physiological module simulating individual tree responses to the environment; (ii) a demographic module simulating tree survival, reproduction, and pollen and seed dispersal; and (iii) a quantitative genetics module controlling the heritability of key life history traits. We used this model to investigate the plastic and genetic components of the variations in the timing of budburst (TBB) along an elevational gradient of Fagus sylvatica (the European beech). We used a repeated 5 years climatic sequence to show that five generations of natural selection were sufficient to develop nonmonotonic genetic differentiation in the TBB along the local climatic gradient but also that plastic variation among different elevations and years was higher than genetic variation. PDG complements theoretical models and provides testable predictions to understand the adaptive potential of tree populations.

Modeling climate change effects on the potential production of French plains forests at the sub-regional level.

We modeled the effects of climate change and two forest management scenarios on wood production and forest carbon balance in French forests using process-based models of forest growth. We combined data from the national forest inventory and soil network survey, which were aggregated over a 50 x 50-km grid, i.e., the spatial resolution of the climate scenario data. We predicted and analyzed the climate impact on potential forest production over the period 1960-2100. All models predicted a slight increase in potential forest yield until 2030-2050, followed by a plateau or a decline around 2070-2100, with overall, a greater increase in yield in northern France than in the south. Gross and net primary productivities were more negatively affected by soil water and atmospheric water vapor saturation deficits in western France because of a more pronounced shift in seasonal rainfall from summer to winter. The rotation-averaged values of carbon flux and production for different forest management options were estimated during four years (1980, 2015, 2045 and 2080). Predictions were made using a two-dimensional matrix covering the range of local soil and climate conditions. The changes in ecosystem fluxes and forest production were explained by the counterbalancing effect of rising CO2 concentration and increasing water deficit. The effect of climate change decreased with rotation length from short rotations with high production rates and low standing biomasses to long rotations with low productivities and greater standing biomasses. Climate effects on productivity, both negative and positive, were greatest on high fertility sites. Forest productivity in northern France was enhanced by climate change, increasingly from west to east, whereas in the southwestern Atlantic region, productivity was reduced by climate change to an increasing degree from west to east.

Modeling annual production and carbon fluxes of a large managed temperate forest using forest inventories, satellite data and field measurements.

We evaluated annual productivity and carbon fluxes over the Fontainebleau forest, a large heterogeneous forest region of 17,000 ha, in terms of species composition, canopy structure, stand age, soil type and water and mineral resources. The model is a physiological process-based forest ecosystem model coupled with an allocation model and a soil model. The simulations were done stand by stand, i.e., 2992 forest management units of simulation. Some input parameters that are spatially variable and to which the model is sensitive were calculated for each stand from forest inventory attributes, a network of 8800 soil pits, satellite data and field measurements. These parameters are: (1) vegetation attributes: species, age, height, maximal leaf area index of the year, aboveground biomass and foliar nitrogen content; and (2) soil attributes: available soil water capacity, soil depth and soil carbon content. Main outputs of the simulations are wood production and carbon fluxes on a daily to yearly basis. Results showed that the forest is a carbon sink, with a net ecosystem exchange of 371 g C m(-2) year(-1). Net primary productivity is estimated at 630 g C m(-2) year(-1) over the entire forest. Reasonably good agreement was found between simulated trunk relative growth rate (2.74%) and regional production estimated from the National Forest Inventory (IFN) (2.52%), as well as between simulated and measured annual wood production at the forest scale (about 71,000 and 68,000 m(3) year(-1), respectively). Results are discussed species by species.