Temporal and spatial variability of phloem structure in Picea abies and Fagus sylvatica and its link to climate.

Using a unique 8-year data set (2010-2017) of phloem data, we studied the effect of temperature and precipitation on the phloem anatomy (conduit area, widths of ring, early and late phloem) and xylem-ring width in two coexisting temperate tree species, Picea abies and Fagus sylvatica, from three contrasting European temperate forest sites. Histometric analyses were performed on microcores taken from tree stems in autumn. We found high interannual variability and sensitivity of phloem anatomy and xylem-ring widths to precipitation and temperature; however, the responses were species- and site-specific. The contrasting response of xylem and phloem-ring widths of the same tree species to weather conditions was found at the two Slovenian sites generally well supplied with precipitation, while at the driest Czech site, the influence of weather factors on xylem and phloem ring widths was synchronised. Since widths of mean annual xylem and phloem increments were narrowest at the Czech site, this site is suggested to be most restrictive for the radial growth of both species. By influencing the seasonal patterns of xylem and phloem development, water availability appears to be the most important determinant of tissue- and species-specific responses to local weather conditions.

A critical thermal transition driving spring phenology of Northern Hemisphere conifers.

Despite growing interest in predicting plant phenological shifts, advanced spring phenology by global climate change remains debated. Evidence documenting either small or large advancement of spring phenology to rising temperature over the spatio-temporal scales implies a potential existence of a thermal threshold in the responses of forests to global warming. We collected a unique data set of xylem cell-wall-thickening onset dates in 20 coniferous species covering a broad mean annual temperature (MAT) gradient (-3.05 to 22.9°C) across the Northern Hemisphere (latitudes 23°-66° N). Along the MAT gradient, we identified a threshold temperature (using segmented regression) of 4.9 ± 1.1°C, above which the response of xylem phenology to rising temperatures significantly decline. This threshold separates the Northern Hemisphere conifers into cold and warm thermal niches, with MAT and spring forcing being the primary drivers for the onset dates (estimated by linear and Bayesian mixed-effect models), respectively. The identified thermal threshold should be integrated into the Earth-System-Models for a better understanding of spring phenology in response to global warming and an improved prediction of global climate-carbon feedbacks.

Photoperiod and temperature as dominant environmental drivers triggering secondary growth resumption in Northern Hemisphere conifers.

Wood formation consumes around 15% of the anthropogenic CO2 emissions per year and plays a critical role in long-term sequestration of carbon on Earth. However, the exogenous factors driving wood formation onset and the underlying cellular mechanisms are still poorly understood and quantified, and this hampers an effective assessment of terrestrial forest productivity and carbon budget under global warming. Here, we used an extensive collection of unique datasets of weekly xylem tissue formation (wood formation) from 21 coniferous species across the Northern Hemisphere (latitudes 23 to 67°N) to present a quantitative demonstration that the onset of wood formation in Northern Hemisphere conifers is primarily driven by photoperiod and mean annual temperature (MAT), and only secondarily by spring forcing, winter chilling, and moisture availability. Photoperiod interacts with MAT and plays the dominant role in regulating the onset of secondary meristem growth, contrary to its as-yet-unquantified role in affecting the springtime phenology of primary meristems. The unique relationships between exogenous factors and wood formation could help to predict how forest ecosystems respond and adapt to climate warming and could provide a better understanding of the feedback occurring between vegetation and climate that is mediated by phenology. Our study quantifies the role of major environmental drivers for incorporation into state-of-the-art Earth system models (ESMs), thereby providing an improved assessment of long-term and high-resolution observations of biogeochemical cycles across terrestrial biomes.

Growth-limiting factors and climate response variability in Norway spruce (Picea abies L.) along an elevation and precipitation gradients in Slovenia.

Norway spruce (Picea abies L.) is among the most sensitive coniferous species to ongoing climate change. However, previous studies on its growth response to increasing temperatures have yielded contrasting results (from stimulation to suppression), suggesting highly site-specific responses. Here, we present the first study that applies two independent approaches, i.e. the nonlinear, process-based Vaganov-Shashkin (VS) model and linear daily response functions. Data were collected at twelve sites in Slovenia differing in climate regimes and ranging elevation between 170 and 1300 m a.s.l. VS model results revealed that drier Norway spruce sites at lower elevations are mostly moisture limited, while moist high-elevation sites are generally more temperature limited. Daily response functions match well the pattern of growth-limiting factors from the VS model and further explain the effect of climate on radial growth: prevailing growth-limiting factors correspond to the climate variable with higher correlations. Radial growth correlates negatively with rising summer temperature and positively with higher spring precipitation. The opposite response was observed for the wettest site at the highest elevation, which positively reacts to increased summer temperature and will most likely benefit from a warming climate. For all other sites, the future radial growth of Norway spruce largely depends on the balance between spring precipitation and summer temperature.

Frost and drought: Effects of extreme weather events on stem carbon dynamics in a Mediterranean beech forest.

The effects of short-term extreme events on tree functioning and physiology are still rather elusive. European beech is one of the most sensitive species to late frost and water shortage. We investigated the intra-annual C dynamics in stems under such conditions. Wood formation and stem CO2 efflux were monitored in a Mediterranean beech forest for 3 years (2015-2017), including a late frost (2016) and a summer drought (2017). The late frost reduced radial growth and, consequently, the amount of carbon fixed in the stem biomass by 80%. Stem carbon dioxide efflux in 2016 was reduced by 25%, which can be attributed to the reduction of effluxes due to growth respiration. Counter to our expectations, we found no effects of the 2017 summer drought on radial growth and stem carbon efflux. The studied extreme weather events had various effects on tree growth. Even though late spring frost had a strong impact on beech radial growth in the current year, trees fully recovered in the following growing season, indicating high resilience of beech to this stressful event.

Inter-individual variability in spring phenology of temperate deciduous trees depends on species, tree size and previous year autumn phenology.

We explored the inter-individual variability in bud-burst and its potential drivers, in homogeneous mature stands of temperate deciduous trees. Phenological observations of leaves and wood formation were performed weekly from summer 2017 to summer 2018 for pedunculate oak, European beech and silver birch in Belgium. The variability of bud-burst was correlated to previous' year autumn phenology (i.e. the onset of leaf senescence and the cessation of wood formation) and tree size but with important differences among species. In fact, variability of bud-burst was primarily related to onset of leaf senescence, cessation of wood formation and tree height for oak, beech and birch, respectively. The inter-individual variability of onset of leaf senescence was not related to the tree characteristics considered and was much larger than the inter-individual variability in bud-burst. Multi-species multivariate models could explain up to 66% of the bud-burst variability. These findings represent an important advance in our fundamental understanding and modelling of phenology and tree functioning of deciduous tree species.

Chilling and forcing temperatures interact to predict the onset of wood formation in Northern Hemisphere conifers.

The phenology of wood formation is a critical process to consider for predicting how trees from the temperate and boreal zones may react to climate change. Compared to leaf phenology, however, the determinism of wood phenology is still poorly known. Here, we compared for the first time three alternative ecophysiological model classes (threshold models, heat-sum models and chilling-influenced heat-sum models) and an empirical model in their ability to predict the starting date of xylem cell enlargement in spring, for four major Northern Hemisphere conifers (Larix decidua, Pinus sylvestris, Picea abies and Picea mariana). We fitted models with Bayesian inference to wood phenological data collected for 220 site-years over Europe and Canada. The chilling-influenced heat-sum model received most support for all the four studied species, predicting validation data with a 7.7-day error, which is within one day of the observed data resolution. We conclude that both chilling and forcing temperatures determine the onset of wood formation in Northern Hemisphere conifers. Importantly, the chilling-influenced heat-sum model showed virtually no spatial bias whichever the species, despite the large environmental gradients considered. This suggests that the spring onset of wood formation is far less affected by local adaptation than by environmentally driven plasticity. In a context of climate change, we therefore expect rising winter-spring temperature to exert ambivalent effects on the spring onset of wood formation, tending to hasten it through the accumulation of forcing temperature, but imposing a higher forcing temperature requirement through the lower accumulation of chilling.

Pattern of xylem phenology in conifers of cold ecosystems at the Northern Hemisphere.

The interaction between xylem phenology and climate assesses forest growth and productivity and carbon storage across biomes under changing environmental conditions. We tested the hypothesis that patterns of wood formation are maintained unaltered despite the temperature changes across cold ecosystems. Wood microcores were collected weekly or biweekly throughout the growing season for periods varying between 1 and 13 years during 1998-2014 and cut in transverse sections for assessing the onset and ending of the phases of xylem differentiation. The data set represented 1321 trees belonging to 10 conifer species from 39 sites in the Northern Hemisphere and covering an interval of mean annual temperature exceeding 14 K. The phenological events and mean annual temperature of the sites were related linearly, with spring and autumnal events being separated by constant intervals across the range of temperature analysed. At increasing temperature, first enlarging, wall-thickening and mature tracheids appeared earlier, and last enlarging and wall-thickening tracheids occurred later. Overall, the period of wood formation lengthened linearly with the mean annual temperature, from 83.7 days at -2 °C to 178.1 days at 12 °C, at a rate of 6.5 days °C-1 . April-May temperatures produced the best models predicting the dates of wood formation. Our findings demonstrated the uniformity of the process of wood formation and the importance of the environmental conditions occurring at the time of growth resumption. Under warming scenarios, the period of wood formation might lengthen synchronously in the cold biomes of the Northern Hemisphere.

Do variations in leaf phenology affect radial growth variations in Fagus sylvatica?

We used a dendrochronological and leaf phenology network of European beech (Fagus sylvatica) in Slovenia, a transitional area between Mediterranean, Alpine and continental climatic regimes, for the period 1955-2007 to test whether year to year variations in leaf unfolding and canopy duration (i.e. time between leaf unfolding and colouring) influence radial growth (annual xylem production and tree ring widths) and if such influences are more pronounced at higher altitudes. We showed that variability in leaf phenology has no significant effect on variations in radial growth. The results are consistent in the entire region, irrespective of the climatic regime or altitude, although previous studies have shown that leaf phenology and tree ring variation depend on altitude. The lack of relationship between year to year variability in leaf phenology and radial growth may suggest that earlier leaf unfolding--as observed in a previous study--probably does not cause increased tree growth rates in beech in Slovenia.

High preseason temperature variability drives convergence of xylem phenology in the Northern Hemisphere conifers.

Wood growth is key to understanding the feedback of forest ecosystems to the ongoing climate warming. An increase in spatial synchrony (i.e., coincident changes in distant populations) of spring phenology is one of the most prominent climate responses of forest trees. However, whether temperature variability contributes to an increase in the spatial synchrony of spring phenology and its underlying mechanisms remains largely unknown. Here, we analyzed an extensive dataset of xylem phenology observations of 20 conifer species from 75 sites over the Northern Hemisphere. Along the gradient of increase in temperature variability in the 75 sites, we observed a convergence in the onset of cell enlargement roughly toward the 5th of June, with a convergence in the onset of cell wall thickening toward the summer solstice. The increase in rainfall since the 5th of June is favorable for cell division and expansion, and as the most hours of sunlight are received around the summer solstice, it allows the optimization of carbon assimilation for cell wall thickening. Hence, the convergences can be considered as the result of matching xylem phenological activities to favorable conditions in regions with high temperature variability. Yet, forest trees relying on such consistent seasonal cues for xylem growth could constrain their ability to respond to climate warming, with consequences for the potential growing season length and, ultimately, forest productivity and survival in the future.

Incorporating high-resolution climate, remote sensing and topographic data to map annual forest growth in central and eastern Europe.

To enhance our understanding of forest carbon sequestration, climate change mitigation and drought impact on forest ecosystems, the availability of high-resolution annual forest growth maps based on tree-ring width (TRW) would provide a significant advancement to the field. Site-specific characteristics, which can be approximated by high-resolution Earth observation by satellites (EOS), emerge as crucial drivers of forest growth, influencing how climate translates into tree growth. EOS provides information on surface reflectance related to forest characteristics and thus can potentially improve the accuracy of forest growth models based on TRW. Through the modelling of TRW using EOS, climate and topography data, we showed that species-specific models can explain up to 52 % of model variance (Quercus petraea), while combining different species results in relatively poor model performance (R2 = 13 %). The integration of EOS into models based solely on climate and elevation data improved the explained variance by 6 % on average. Leveraging these insights, we successfully generated a map of annual TRW for the year 2021. We employed the area of applicability (AOA) approach to delineate the range in which our models are deemed valid. The calculated AOA for the established forest-type models was 73 % of the study region, indicating robust spatial applicability. Notably, unreliable predictions predominantly occurred in the climate margins of our dataset. In conclusion, our large-scale assessment underscores the efficacy of combining climate, EOS and topographic data to develop robust models for mapping annual TRW. This research not only fills a critical void in the current understanding of forest growth dynamics but also highlights the potential of integrated data sources for comprehensive ecosystem assessments.

A meta-analysis of cambium phenology and growth: linear and non-linear patterns in conifers of the northern hemisphere.

BACKGROUND AND AIMS: Ongoing global warming has been implicated in shifting phenological patterns such as the timing and duration of the growing season across a wide variety of ecosystems. Linear models are routinely used to extrapolate these observed shifts in phenology into the future and to estimate changes in associated ecosystem properties such as net primary productivity. Yet, in nature, linear relationships may be special cases. Biological processes frequently follow more complex, non-linear patterns according to limiting factors that generate shifts and discontinuities, or contain thresholds beyond which responses change abruptly. This study investigates to what extent cambium phenology is associated with xylem growth and differentiation across conifer species of the northern hemisphere. METHODS: Xylem cell production is compared with the periods of cambial activity and cell differentiation assessed on a weekly time scale on histological sections of cambium and wood tissue collected from the stems of nine species in Canada and Europe over 1-9 years per site from 1998 to 2011. KEY RESULTS: The dynamics of xylogenesis were surprisingly homogeneous among conifer species, although dispersions from the average were obviously observed. Within the range analysed, the relationships between the phenological timings were linear, with several slopes showing values close to or not statistically different from 1. The relationships between the phenological timings and cell production were distinctly non-linear, and involved an exponential pattern. CONCLUSIONS: The trees adjust their phenological timings according to linear patterns. Thus, shifts of one phenological phase are associated with synchronous and comparable shifts of the successive phases. However, small increases in the duration of xylogenesis could correspond to a substantial increase in cell production. The findings suggest that the length of the growing season and the resulting amount of growth could respond differently to changes in environmental conditions.

Upside down and the game of C allocation.

Non-structural carbohydrates (NSCs) represent the primary carbon (C) reserves and play a crucial role for plant functioning and resilience. Indeed, these compounds are involved in the regulation between C supply and demand, and in the maintenance of hydraulic efficiency. NSCs are stored in parenchyma of woody organs, which is recognized as a proxy for reserve storage capacity of tree. Notwithstanding the importance of NSCs for tree physiology, their long-term regulation and trade-offs against growth were not deeply investigated. This work evaluated the long-term dynamics of mature tree reserves in stem and root, proxied by parenchyma features, and focusing on the trade off and interplay between the resources allocation in radial growth and reserves in stem and coarse root. In a Mediterranean beech forest, NSCs content, stem and root wood anatomy analysis, and eddy covariance data, were combined. The parenchyma fraction (RAP) of beech root and stem was different, due to differences in axial parenchyma (AP) and narrow ray parenchyma (nRP) fractions. However, these parenchyma components and radial growth showed synchronous inter-annual dynamics between the two organs. In beech stem, positive correlations were found among soluble sugars content and nRP, and among starch content and the AP. Positive correlations were found among Net Ecosystem Exchange (NEE) and AP of both organs. In contrast, NEE was negatively correlated to radial growth of root and stem. Our results suggest a different contribution of stem and roots to reserves storage, and a putative partitioning in the functional roles of parenchyma components. Moreover, a long-term trade-off of C allocation between growth and reserve pool was evidenced. Indeed, in case of C source reduction, trees preferentially allocate C towards reserves pool. Conversely, in high productive years, growth represents the major C sink.

Relationships Between Wood-Anatomical Features and Resistance Drilling Density in Norway Spruce and European Beech.

Environmental conditions affect tree-ring width (TRW), wood structure, and, consequently, wood density, which is one of the main wood quality indicators. Although studies on inter- and intra-annual variability in tree-ring features or density exist, studies demonstrating a clear link between wood structure on a cellular level and its effect on wood density on a macroscopic level are rare. Norway spruce with its simple coniferous structure and European beech, a diffuse-porous angiosperm species were selected to analyze these relationships. Increment cores were collected from both species at four sites in Slovenia. In total, 24 European beech and 17 Norway spruce trees were sampled. In addition, resistance drilling measurements were performed just a few centimeters above the increment core sampling. TRW and quantitative wood anatomy measurements were performed on the collected cores. Resistance drilling density values, tree-ring (TRW, earlywood width-EWW, transition-TWW, and latewood width-LWW) and wood-anatomical features (vessel/tracheid area and diameter, cell density, relative conductive area, and cell wall thickness) were then averaged for the first 7 cm of measurements. We observed significant relationships between tree-ring and wood-anatomical features in both spruce and beech. In spruce, the highest correlation values were found between TRW and LWW. In beech, the highest correlations were observed between TRW and cell density. There were no significant relationships between wood-anatomical features and resistance drilling density in beech. However, in spruce, a significant negative correlation was found between resistance drilling density and tangential tracheid diameter, and a positive correlation between resistance drilling density and both TWW + LWW and LWW. Our findings suggest that resistance drilling measurements can be used to evaluate differences in density within and between species, but they should be improved in resolution to be able to detect changes in wood anatomy.

Climatic regulation of leaf and cambial phenology in Quercus pubescens: Their interlinkage and impact on xylem and phloem conduits.

Increased frequency and severity of stressful events affects the growth patterns and functioning of trees which adjust their phenology to given conditions. Here, we analysed environmental effects (temperature, precipitation, VPD and SWC) on the timing of leaf phenology, seasonal stem radial growth patterns, and xylem and phloem anatomy of Quercus pubescens in the sub-Mediterranean in the period 2014-2019, when various adverse weather events occurred, i.e. spring drought in 2015, summer fire in 2016 and summer drought in 2017. Results showed that the timings of leaf and cambium phenology do not occur simultaneously in Q. pubescens, reflecting different environmental and internal constraints. Although year-to-year variability in the timings of leaf and cambial phenology exists, their chronological sequence is fairly fixed. Different effects of weather conditions on different stages of leaf development in spring were observed. Common climatic drivers (i.e., negative effect of hot and dry summers and a positive effect of increasing moisture availability in winter and summer) were found to affect the widths of xylem and phloem increments with more pronounced effect on late formed parts. A legacy effect of the timing of leaf and cambial phenology of the previous growing season on the timing of phenology of the following spring was confirmed. Rarely available phloem data permitted a comprehensive insight into the interlinkage of the timing of cambium and leaf phenology and adjustment strategies of vascular tissues in Mediterranean pubescent oak to various environmental constraints, including frequent extreme events (drought, fire). Our results suggest that predicted changes in autumn/winter and spring climatic conditions for this area could affect the timings of leaf and stem cambial phenology of Q. pubescens in the coming years, which would affect stem xylem and phloem structure and hydraulic properties, and ultimately its performance.

Climate-change-driven growth decline of European beech forests.

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 21st century growth changes across Europe indicate serious ecological and economic consequences that require immediate forest adaptation.

Jet stream position explains regional anomalies in European beech forest productivity and tree growth.

The mechanistic pathways connecting ocean-atmosphere variability and terrestrial productivity are well-established theoretically, but remain challenging to quantify empirically. Such quantification will greatly improve the assessment and prediction of changes in terrestrial carbon sequestration in response to dynamically induced climatic extremes. The jet stream latitude (JSL) over the North Atlantic-European domain provides a synthetic and robust physical framework that integrates climate variability not accounted for by atmospheric circulation patterns alone. Surface climate impacts of north-south summer JSL displacements are not uniform across Europe, but rather create a northwestern-southeastern dipole in forest productivity and radial-growth anomalies. Summer JSL variability over the eastern North Atlantic-European domain (5-40E) exerts the strongest impact on European beech, inducing anomalies of up to 30% in modelled gross primary productivity and 50% in radial tree growth. The net effects of JSL movements on terrestrial carbon fluxes depend on forest density, carbon stocks, and productivity imbalances across biogeographic regions.

Different Wood Anatomical and Growth Responses in European Beech (Fagus sylvatica L.) at Three Forest Sites in Slovenia.

European beech (Fagus sylvatica L.) adapts to local growing conditions to enhance its performance. In response to variations in climatic conditions, beech trees adjust leaf phenology, cambial phenology, and wood formation patterns, which result in different tree-ring widths (TRWs) and wood anatomy. Chronologies of tree ring width and vessel features [i.e., mean vessel area (MVA), vessel density (VD), and relative conductive area (RCTA)] were produced for the 1960-2016 period for three sites that differ in climatic regimes and spring leaf phenology (two early- and one late-flushing populations). These data were used to investigate long-term relationships between climatic conditions and anatomical features of four quarters of tree-rings at annual and intra-annual scales. In addition, we investigated how TRW and vessel features adjust in response to extreme weather events (i.e., summer drought). We found significant differences in TRW, VD, and RCTA among the selected sites. Precipitation and maximum temperature before and during the growing season were the most important climatic factors affecting TRW and vessel characteristics. We confirmed differences in climate-growth relationships between the selected sites, late flushing beech population at Idrija showing the least pronounced response to climate. MVA was the only vessel trait that showed no relationship with TRW or other vessel features. The relationship between MVA and climatic factors evaluated at intra-annual scale indicated that vessel area in the first quarter of tree-ring were mainly influenced by climatic conditions in the previous growing season, while vessel area in the second to fourth quarters of tree ring width was mainly influenced by maximum temperature and precipitation in the current growing season. When comparing wet and dry years, beech from all sites showed a similar response, with reduced TRW and changes in intra-annual variation in vessel area. Our findings suggest that changes in temperature and precipitation regimes as predicted by most climate change scenarios will affect tree-ring increments and wood structure in beech, yet the response between sites or populations may differ.

Timing of spring xylogenesis in temperate deciduous tree species relates to tree growth characteristics and previous autumn phenology.

We explored the timing of spring xylogenesis and its potential drivers in homogeneous mature forest stands in a temperate European region. Three species with contrasting leaf development dynamics and wood anatomy were studied: European beech, silver birch and pedunculate oak. Detailed phenological observations of xylogenesis and leaf phenology were performed from summer 2017 until spring 2018. Cambium reactivation (CR) occurred before the buds of oak and birch were swollen, whereas these two phenological phases were concurrent for beech. On the other hand, initial earlywood vessels were fully differentiated (FDIEV) after leaf unfolding for all three species. Timing of CR was correlated to average ring-width of the last 10 years (2008-17), tree diameter and, partially, with tree age. In addition, the timing of FDIEV was correlated to tree age and previous year's autumn phenology, i.e., timing of wood growth cessation and onset of leaf senescence. Multivariate models could explain up to 68% of the variability of CR and 55% of the variability of FDIEV. In addition to the 'species' factor, the variability could be explained by ca 30% by tree characteristics and previous year's autumn phenology for both CR and FDIEV. These findings are important to better identify which factors (other than environment) can be driving the onset of the growing season, and highlight the influence of tree growth characteristics and previous year's phenology on spring wood phenology, wood formation and, potentially, forest production.