The 2018 European heatwave led to stem dehydration but not to consistent growth reductions in forests.

Heatwaves exert disproportionately strong and sometimes irreversible impacts on forest ecosystems. These impacts remain poorly understood at the tree and species level and across large spatial scales. Here, we investigate the effects of the record-breaking 2018 European heatwave on tree growth and tree water status using a collection of high-temporal resolution dendrometer data from 21 species across 53 sites. Relative to the two preceding years, annual stem growth was not consistently reduced by the 2018 heatwave but stems experienced twice the temporary shrinkage due to depletion of water reserves. Conifer species were less capable of rehydrating overnight than broadleaves across gradients of soil and atmospheric drought, suggesting less resilience toward transient stress. In particular, Norway spruce and Scots pine experienced extensive stem dehydration. Our high-resolution dendrometer network was suitable to disentangle the effects of a severe heatwave on tree growth and desiccation at large-spatial scales in situ, and provided insights on which species may be more vulnerable to climate extremes.

The within-population variability of leaf spring and autumn phenology is influenced by temperature in temperate deciduous trees.

Leaf phenology is a major driver of ecosystem functioning in temperate forests and a robust indicator of climate change. Both the inter-annual and inter-population variability of leaf phenology have received much attention in the literature; in contrast, the within-population variability of leaf phenology has been far less studied. Beyond its impact on individual tree physiological processes, the within-population variability of leaf phenology can affect the estimation of the average budburst or leaf senescence dates at the population scale. Here, we monitored the progress of spring and autumn leaf phenology over 14 tree populations (9 tree species) in six European forests over the period of 2011 to 2018 (yielding 16 site-years of data for spring, 14 for autumn). We monitored 27 to 512 (with a median of 62) individuals per population. We quantified the within-population variability of leaf phenology as the standard deviation of the distribution of individual dates of budburst or leaf senescence (SDBBi and SDLSi, respectively). Given the natural variability of phenological dates occurring in our tree populations, we estimated from the data that a minimum sample size of 28 (resp. 23) individuals, are required to estimate SDBBi (resp. SDLSi) with a precision of 3 (resp. 7) days. The within-population of leaf senescence (average SDLSi = 8.5 days) was on average two times larger than for budburst (average SDBBi = 4.0 days). We evidenced that warmer temperature during the budburst period and a late average budburst date were associated with a lower SDBBi, as a result of a quicker spread of budburst in tree populations, with a strong species effect. Regarding autumn phenology, we observed that later senescence and warm temperatures during the senescence period were linked with a high SDLSi, with a strong species effect. The shares of variance explained by our models were modest suggesting that other factors likely influence the within-population variation in leaf phenology. For instance, a detailed analysis revealed that summer temperatures were negatively correlated with a lower SDLSi.

"Green pointillism": detecting the within-population variability of budburst in temperate deciduous trees with phenological cameras.

Phenological cameras have been used over a decade for identifying plant phenological markers (budburst, leaf senescence) and more generally the greenness dynamics of forest canopies. The analysis is usually carried out over the full camera field of view, with no particular analysis of the variability of phenological markers among trees. Here we show that images produced by phenological cameras can be used to quantify the within-population variability of budburst (WPVbb) in temperate deciduous forests. Using seven site-years of image analyses, we report a strong correlation (r2 = 0.97) between the WPVbb determined with a phenological camera and its quantification through ground observation. We show that WPVbb varies strongly (by a factor of 4) from year to year in a given population and that those variations are linked with temperature conditions during the budburst period, with colder springs associated to a higher differentiation of budburst (higher WPVbb) among trees. Deploying our approach at the continental scale, i.e., throughout phenological cameras networks, would improve the understanding of the spatial (across populations) and temporal (across years) variations of WPVbb, which have strong implications on forest functioning, tree fitness and phenological modelling.

Potential use of the PRI and active fluorescence for the diagnosis of the physiological state of plants under ozone exposure and high atmospheric vapor pressure deficit.

Assessing photosynthesis rates with remote sensing is important for tracking the physiological state of plants. The photochemical reflectance index (PRI) is a good estimator of short-term light-use efficiency (LUE) at the leaf scale but its responses to environmental factors are poorly understood. In this study, we assessed changes in the responses of the PRI to ozone exposure and to an increase in atmospheric drought (separately and combined) in oak (Quercus robur) and holm oak (Quercus ilex) that were planted in climatic cells under controlled conditions. The aim was to evaluate the ability of the PRI as a relevant indicator to assess the impact of abiotic factors on photosynthesis. Leaf-scale measurements of biochemical, physiological and spectral properties, including the PRI in dim light on dark-adapted leaves (PRI0), kinetics of PRI responses to PAR variations (photosynthetically active radiation), and leaf chlorophyll fluorescence parameters were performed. The results show that PRI0 is a good proxy of the leaf chlorophyll content, and is correlated to chlorophyll fluorescence parameters on dark adapted leaves (Fo, Fm). The correction of the PRI from the leaf chlorophyll content variations (PRIc) significantly improves correlations between the PRI and NPQ (non-photochemical quenching). The variability of PARsat (estimated PAR value at PRI saturation using PRI vs. PAR relationships) depends on ozone exposure and on the increase in atmospheric vapor pressure deficit. For Quercus robur, results highlight that PARsat is linked to abiotic stress indicating that the PRI may be used as a relevant indicator of abiotic factors limiting the photosynthesis. Quercus ilex did not show significant variability in PRI0 and PARsat, which suggest that it is a more drought resistant species than Q. robur.

Environmental control of carbon allocation matters for modelling forest growth.

We aimed to evaluate the importance of modulations of within-tree carbon (C) allocation by water and low-temperature stress for the prediction of annual forest growth with a process-based model. A new C allocation scheme was implemented in the CASTANEA model that accounts for lagged and direct environmental controls of C allocation. Different approaches (static vs dynamic) to modelling C allocation were then compared in a model-data fusion procedure, using satellite-derived leaf production estimates and biometric measurements at c. 104 sites. The modelling of the environmental control of C allocation significantly improved the ability of CASTANEA to predict the spatial and year-to-year variability of aboveground forest growth along regional gradients. A significant effect of the previous year's water stress on the C allocation to leaves and wood was reported. Our results also are consistent with a prominent role of the environmental modulation of sink demand in the wood growth of the studied species. Data available at large scales can inform forest models about the processes driving annual and seasonal C allocation. Our results call for a greater consideration of C allocation drivers, especially sink-demand fluctuations, for the simulations of current and future forest productivity with process-based models.

Wood phenology, not carbon input, controls the interannual variability of wood growth in a temperate oak forest.

Although the analysis of flux data has increased our understanding of the interannual variability of carbon inputs into forest ecosystems, we still know little about the determinants of wood growth. Here, we aimed to identify which drivers control the interannual variability of wood growth in a mesic temperate deciduous forest. We analysed a 9-yr time series of carbon fluxes and aboveground wood growth (AWG), reconstructed at a weekly time-scale through the combination of dendrometer and wood density data. Carbon inputs and AWG anomalies appeared to be uncorrelated from the seasonal to interannual scales. More than 90% of the interannual variability of AWG was explained by a combination of the growth intensity during a first 'critical period' of the wood growing season, occurring close to the seasonal maximum, and the timing of the first summer growth halt. Both atmospheric and soil water stress exerted a strong control on the interannual variability of AWG at the study site, despite its mesic conditions, whilst not affecting carbon inputs. Carbon sink activity, not carbon inputs, determined the interannual variations in wood growth at the study site. Our results provide a functional understanding of the dependence of radial growth on precipitation observed in dendrological studies.

Explaining the variability of the photochemical reflectance index (PRI) at the canopy-scale: Disentangling the effects of phenological and physiological changes.

Assessing photosynthesis rates at the ecosystem scale and over large regions is important for tracking the global carbon cycle and remote sensing has provided new and useful approaches for performing this assessment. The photochemical reflectance index (PRI) is a good estimator of short-term light-use efficiency (LUE) at the leaf scale; however, confounding factors appear at larger temporal and spatial scales. In this study, canopy-scale PRI variability was investigated for three species (Fagus sylvatica L., Quercus robur L. and Pinus sylvestris L.) growing under contrasting soil moisture conditions. Throughout the growing season, no significant differences in chlorophyll content and in violaxanthin, antheraxanthin and zeaxanthin were found between species or treatments. The daily PRI vs PAR (photosynthetically active radiation) relationships were determined using continuous measurements obtained at high frequency throughout the entire growing season, from early spring budburst to later autumn senescence, and were used to deconvolute the physiological PRI variability related to LUE variations due to phenological variability and related to temporal changes in the biochemical and structural canopy attributes. The PRI vs PAR relationship is used to show that the canopy-scale PRI measured at low radiation depends on the chlorophyll content of the canopy. The range of PRI variations at an intra-daily scale and the dynamics of the xanthophyll pool do not vary between days, which suggests that the PRI responds to a xanthophyll ratio. The PAR values at PRI saturation are mainly related to the canopy chlorophyll content during budburst and senescence and to the soil moisture content when the chlorophyll content is no longer a limiting factor. This parameter is significantly lower in the oak species that experience less stress from variations in soil moisture and is species dependant. These results provide new insights regarding the analysis and the meaning of PRI variability as a proxy for LUE at the canopy scale.

Seasonal changes in carbon and nitrogen compound concentrations in a Quercus petraea chronosequence.

Forest productivity declines with tree age. This decline may be due to changes in metabolic functions, resource availability and/or changes in resource allocation (between growth, reproduction and storage) with tree age. Carbon and nitrogen remobilization/storage processes are key to tree growth and survival. However, studies of the effects of tree age on these processes are scarce and have not yet considered seasonal carbon and nitrogen variations in situ. This study was carried out in a chronosequence of sessile oak (Quercus petraea Liebl.) for 1 year to survey the effects of tree age on the seasonal changes of carbon and nitrogen compounds in several tree compartments, focusing on key phenological stages. Our results highlight a general pattern of carbon and nitrogen function at all tree ages, with carbon reserve remobilization at budburst for growth, followed by carbon reserve formation during the leafy season and carbon reserve use during winter for maintenance. The variation in concentrations of nitrogen compounds shows less amplitude than that of carbon compounds. Storage as proteins occurs later, and mainly depends on leaf nitrogen remobilization and root uptake in autumn. We highlight several differences between tree age groups, in particular the loss of carbon storage function of fine and medium-sized roots with tree ageing. Moreover, the pattern of carbon compound accumulation in branches supports the hypothesis of a preferential allocation of carbon towards growth until the end of wood formation in juvenile trees, at the expense of the replenishment of carbon stores, while mature trees start allocating carbon to storage right after budburst. Our results demonstrate that at key phenological stages, physiological and developmental functions differ with tree age, and together with environmental conditions, influence the carbon and nitrogen concentration variations in sessile oaks.

Assessing the effects of management on forest growth across France: insights from a new functional-structural model.

BACKGROUND AND AIMS: The structure of a forest stand, i.e. the distribution of tree size features, has strong effects on its functioning. The management of the structure is therefore an important tool in mitigating the impact of predicted changes in climate on forests, especially with respect to drought. Here, a new functional-structural model is presented and is used to assess the effects of management on forest functioning at a national scale. METHODS: The stand process-based model (PBM) CASTANEA was coupled to a stand structure module (SSM) based on empirical tree-to-tree competition rules. The calibration of the SSM was based on a thorough analysis of intersite and interannual variability of competition asymmetry. The coupled CASTANEA-SSM model was evaluated across France using forest inventory data, and used to compare the effect of contrasted silvicultural practices on simulated stand carbon fluxes and growth. KEY RESULTS: The asymmetry of competition varied consistently with stand productivity at both spatial and temporal scales. The modelling of the competition rules enabled efficient prediction of changes in stand structure within the CASTANEA PBM. The coupled model predicted an increase in net primary productivity (NPP) with management intensity, resulting in higher growth. This positive effect of management was found to vary at a national scale across France: the highest increases in NPP were attained in forests facing moderate to high water stress; however, the absolute effect of management on simulated stand growth remained moderate to low because stand thinning involved changes in carbon allocation at the tree scale. CONCLUSIONS: This modelling approach helps to identify the areas where management efforts should be concentrated in order to mitigate near-future drought impact on national forest productivity. Around a quarter of the French temperate oak and beech forests are currently in zones of high vulnerability, where management could thus mitigate the influence of climate change on forest yield.

Comparing the intra-annual wood formation of three European species (Fagus sylvatica, Quercus petraea and Pinus sylvestris) as related to leaf phenology and non-structural carbohydrate dynamics.

Monitoring cambial phenology and intra-annual growth dynamics is a useful approach for characterizing the tree growth response to climate change. However, there have been few reports concerning intra-annual wood formation in lowland temperate forests with high time resolution, especially for the comparison between deciduous and coniferous species. The main objective of this study was to determine how the timing, duration and rate of radial growth change between species as related to leaf phenology and the dynamics of non-structural carbohydrates (NSC) under the same climatic conditions. We studied two deciduous species, Fagus sylvatica L. and Quercus petraea (Matt.) Liebl., and an evergreen conifer, Pinus sylvestris L. During the 2009 growing season, we weekly monitored (i) the stem radial increment using dendrometers, (ii) the xylem growth using microcoring and (iii) the leaf phenology from direct observations of the tree crowns. The NSC content was also measured in the eight last rings of the stem cores in April, June and August 2009. The leaf phenology, NSC storage and intra-annual growth were clearly different between species, highlighting their contrasting carbon allocation. Beech growth began just after budburst, with a maximal growth rate when the leaves were mature and variations in the NSC content were low. Thus, beech radial growth seemed highly dependent on leaf photosynthesis. For oak, earlywood quickly developed before budburst, which probably led to the starch decrease quantified in the stem from April to June. For pine, growth began before the needles unfolding and the lack of NSC decrease during the growing season suggested that the substrates for radial growth were new assimilates of the needles from the previous year. Only for oak, the pattern determined from the intra-annual growth measured using microcoring differed from the pattern determined from dendrometer data. For all species, the ring width was significantly influenced by growth duration and not by growth rate, which differs from previous studies. The observed between-species difference at the intra-annual scale is key information for anticipating suitability of future species in temperate forests.

Climate change impacts on tree ranges: model intercomparison facilitates understanding and quantification of uncertainty.

Model-based projections of shifts in tree species range due to climate change are becoming an important decision support tool for forest management. However, poorly evaluated sources of uncertainty require more scrutiny before relying heavily on models for decision-making. We evaluated uncertainty arising from differences in model formulations of tree response to climate change based on a rigorous intercomparison of projections of tree distributions in France. We compared eight models ranging from niche-based to process-based models. On average, models project large range contractions of temperate tree species in lowlands due to climate change. There was substantial disagreement between models for temperate broadleaf deciduous tree species, but differences in the capacity of models to account for rising CO(2) impacts explained much of the disagreement. There was good quantitative agreement among models concerning the range contractions for Scots pine. For the dominant Mediterranean tree species, Holm oak, all models foresee substantial range expansion.

Comparison of seasonal variations in water-use efficiency calculated from the carbon isotope composition of tree rings and flux data in a temperate forest.

Tree-ring δ(13) C is often interpreted in terms of intrinsic water-use efficiency (WUE) using a carbon isotope discrimination model established at the leaf level. We examined whether intra-ring δ(13) C could be used to assess variations in intrinsic WUE (W(g), the ratio of carbon assimilation and stomatal conductance to water) and variations in ecosystem WUE (W(t) , the ratio of C assimilation and transpiration) at a seasonal scale. Intra-ring δ(13) C was measured in 30- to 60-µm-thick slices in eight oak trees (Quercus petraea). Canopy W(g) was simulated using a physiologically process-based model. High between-tree variability was observed in the seasonal variations of intra-ring δ(13) C. Six trees showed significant positive correlations between W(g) calculated from intra-ring δ(13) C and canopy W(g) averaged over several days during latewood formation. These results suggest that latewood is a seasonal recorder of W(g) trends, with a temporal lag corresponding to the mixing time of sugars in the phloem. These six trees also showed significant negative correlations between photosynthetic discrimination Δ calculated from intra-ring δ(13) C, and ecosystem W(t), during latewood formation. Despite the observed between-tree variability, these results indicate that intra-ring δ(13) C can be used to access seasonal variations in past W(t).

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.

Estimation of deciduous forest leaf area index using direct and indirect methods.

This study evaluated one semi-direct and three indirect methods for estimating leaf area index (LAI) by comparing these estimates with direct estimates derived from litter collection. The semi-direct method uses a thin metallic needle to count a number of contacts across fresh litter layers. One indirect method is based on the penetration of diffuse global radiation measured over the course of a day. The second indirect method uses the LAI-2000 plant canopy analyser (PCA) which measures diffuse light penetration from five different sky sectors simultaneously. The third indirect method uses the "Demon" portable light sensor to measure the penetration of direct beam sunlight at different zenith angles over the course of half a day. The Poisson model of gap frequency was applied to estimate plant area index (PAI) from observed transmittances using the second and third methods. Litter collection from 11 temperate decidous forests gave values of LAI ranging from 1.7 to 7.5. Estimates based on the needle method showed a significant linear relationship with LAI values obtained from litter collections but were systematically lower (by 6-37%). PAI estimates using all three indirect techniques (fixed light sensor system, LAI-2000 and Demon) showed a strong linear relationship with LAI derived from litter collection. Differences, averaged over all forest stands, between PAI estimates from each of the three indirect methods and LAI from litter collections were below 2%. If we consider that LAI=PAI-WAI (wood area index) then, all three indirect methods underestimated LAI by an additional factor close to the value of WAI. One reason could be a local clumping of architectural canopy components: in particular, the spatial dispositions of branchlets and leaves are not independent, leading to a non-random relationship between the distributions of these two canopy components.