Wood anatomical traits in black spruce reveal latent water constraints on the boreal forest

Future carbon sequestration potential in a widespread transcontinental boreal tree species: Standing genetic variation matters!

Climate change (CC) necessitates reforestation/afforestation programs to mitigate its impacts and maximize carbon sequestration. But comprehending how tree growth, a proxy for fitness and resilience, responds to CC is critical to maximize these programs' effectiveness. Variability in tree response to CC across populations can notably be influenced by the standing genetic variation encompassing both neutral and adaptive genetic diversity. Here, a framework is proposed to assess tree growth potential at the population scale while accounting for standing genetic variation. We applied this framework to black spruce (BS, Picea mariana [Mill] B.S.P.), with the objectives to (1) determine the key climate variables having impacted BS growth response from 1974 to 2019, (2) examine the relative roles of local adaptation and the phylogeographic structure in this response, and (3) project BS growth under two Shared Socioeconomic Pathways while taking standing genetic variation into account. We modeled growth using a machine learning algorithm trained with dendroecological and genetic data obtained from over 2600 trees (62 populations divided in three genetic clusters) in four 48-year-old common gardens, and simulated growth until year 2100 at the common garden locations. Our study revealed that high summer and autumn temperatures negatively impacted BS growth. As a consequence of warming, this species is projected to experience a decline in growth by the end of the century, suggesting maladaptation to anticipated CC and a potential threat to its carbon sequestration capacity. This being said, we observed a clear difference in response to CC within and among genetic clusters, with the western cluster being more impacted than the central and eastern clusters. Our results show that intraspecific genetic variation, notably associated with the phylogeographic structure, must be considered when estimating the response of widespread species to CC.

Intra-annual stem radial growth of Qinghai spruce and its environmental drivers in the Qilian Mountains, northwestern China

Tree stem radial growth could be used to estimate forest productivity, which plays a dominant role in the carbon sink of terrestrial ecosystems. However, it is still obscure how intra-annual stem radial growth is regulated by environmental variables. Here, we monitored Qinghai spruce stem radial growth over seven years and analyzed the environmental drivers of the intra-annual stem radial changes in the Qilian Mountains at low (2700 m) and high altitudes (3200 m). We found that stem radial growth initiated when the daily mean minimum air temperature reached 1.6oC, while the cessation of stem growth was unrelated to temperatures and water conditions. Initiations of stem growth at 2700 m were significantly earlier than that at 3200 m. Maximum growth rates were observed before the summer solstice at low altitude, whereas at high altitude, the majority of them occurred after the summer solstice. Most variability in annual stem increment (AI) can be explained by the rate (Rm) than by the duration of stem growth (∆t), and 78.9 % and 69.6 % of the variability in AI were attributable to Rm for the lower and upper site, respectively. Structural equation modeling revealed that precipitation (P) could both directly positively influence stem radial increment (SRI) and indirectly positively influence SRI through influencing relative humidity (RH), but the positive effect of P on SRI was higher at low altitude than at high altitude. Daily minimum air temperature (Tmin) was also the main direct diver of SRI, and the positive effect of Tmin on SRI was higher at high altitude than at low altitude. Considering the trends in climate warming and humidification over the past decades, climate changes would result in earlier initiation of Qinghai spruce stem growth and promote the growth through positive response to increased precipitation in low altitude and through elevated temperature in high altitude, respectively.

Impact of intra-annual wood density fluctuation on tree hydraulic function: insights from a continuous monitoring approach

Climate change significantly impacts global forests, leading to tree decline and dieback. To cope with climate change, trees develop several functional traits, such as intra-annual density fluctuations (IADFs) in tree rings. The formation of these traits facilitates trees to optimize resource allocation, allowing them to withstand periods of stress and eventually recover when the conditions become favourable again. This study focuses on a Pinus pinaster Aiton forest in a warm, drought-prone Mediterranean area, comparing two growing seasons with different weather patterns. The innovative continuous monitoring approach used in this study combines high-resolution monitoring of sap flow (SF), analysis of xylogenesis and quantitative wood anatomy. Our results revealed the high plasticity of P. pinaster in water use and wood formation, shedding light on the link between IADFs and tree conductance. Indeed, the capacity to form large cells in autumn (as IADFs) improves the total xylem hydraulic conductivity of this species. For the first time, a continuous SF measurement system captured the dynamics of bimodal SF during the 2022 growing season in conjunction with the bimodal growth pattern observed through xylogenesis monitoring. These results highlight the intricate interplay between environmental conditions, water use, wood formation and tree physiology, providing valuable insights into the acclimation mechanisms employed by P. pinaster to cope with weather fluctuations.

Climate response of a glacial relict conifer across its distribution range is invariant in space but not in time

Climate change impacts on forest trees will be particularly severe for relict species endemic to the subalpine forest, such as Pinus cembra in the Alps and Carpathians. Most current knowledge about the response of this species to climate comes from tree-ring width analysis. However, this approach cannot perform in-depth and highly time-resolved analysis on the climate influence on specific growth processes and xylem functions. We analyzed xylem anatomical traits from six sites covering most of the longitudinal range of this species. Associations between climate and cell number, lumen area and cell wall thickness were computed for the 1920-2010 period using climate records aligned to degree-day temperature sum thresholds. The anatomical chronologies were clearly distinct between the Alps and Carpathians. However, climate responses were similar for all sites, suggesting common species-specific response mechanisms. Temperature showed a positive correlation with both cell number and cell wall thickness. Cell lumen size exhibited an early positive association, followed by strong negative association with temperature and a positive one with precipitation. This highlights that the cell enlargement process was negatively related to high temperature at high elevation, where meristematic processes are rather supposed to be constrained by low temperatures. Therefore, long-term climate warming can have negative consequences on the xylem potential to transport water at all investigated sites. Moreover, in the last 30 years, we observed a slight anticipation of some responses and a decrease in climate sensitivity of some xylem parameters. Our findings provide evidence of temporally unstable but spatially consistent climate response of Pinus cembra from the Alps to the Carpathians. The low diversity in xylem phenotypic responses to climate suggests that future warming could extensively and evenly affect the species throughout its entire distribution.

Climatic responses and variability in bark anatomical traits of 23 Picea species

In woody plants, bark is an important protective tissue which can participate in photosynthesis, manage water loss, and transport assimilates. Studying the bark anatomical traits can provide insight into plant environmental adaptation strategies. However, a systematic understanding of the variability in bark anatomical traits and their drivers is lacking in woody plants. In this study, the bark anatomical traits of 23 Picea species were determined in a common garden experiment. We analyzed interspecific differences and interpreted the patterns in bark anatomical traits in relation to phylogenetic relationships and climatic factors of each species according to its global distribution. The results showed that there were interspecific differences in bark anatomical traits of Picea species. Phloem thickness was positively correlated with parenchyma cell size, possibly related to the roles of parenchyma cells in the radial transport of assimilates. Sieve cell size was negatively correlated with the radial diameter of resin ducts, and differences in sieve cells were possibly related to the formation and expansion of resin ducts. There were no significant phylogenetic signals for any bark anatomical trait, except the tangential diameter of resin ducts. Phloem thickness and parenchyma cell size were affected by temperature-related factors of their native range, while sieve cell size was influenced by precipitation-related factors. Bark anatomical traits were not significantly different under wet and dry climates. This study makes an important contribution to our understanding of variability in bark anatomical traits among Picea species and their ecological adaptations.

Influence of monsoon anomalies on intra-annual density fluctuations of Chinese pine in the Loess Plateau

In the past few decades, the East Asian summer monsoon (EASM) has experienced an unprecedented weakening, exacerbating drought in northern China, especially in the monsoon margin area. Improving our understanding of monsoon variability will benefit agricultural production, ecological construction, and disaster management. Tree-ring is widely used as proxy data for extending the monsoon history. However, in the East Asian monsoon margin, the tree-ring width were mostly formed before the rainy season, thus may have limited ability to indicate the monsoon variability. Intra-annual density fluctuations (IADFs) can provide higher resolution information on tree growth as well as evidence of short-term climate events. Here, we used Chinese pine (Pinus tabuliformis Carr.) samples from the eastern edge of the Chinese Loess Plateau (CLP), where the climate is deeply affected by monsoon, to investigate the response of tree growth and IADFs frequency to climate variation. We show that tree-ring width and IADFs record significantly different climatic signals. The former was mainly affected by moisture conditions at the end of the previous growing season and the current spring. While the latter was common in years when severe droughts occurred in June and July, especially in June. This period coincides with the onset of the EASM, so we further analyzed the relationship between IADFs frequency and the rainy season. Both correlation analysis and the GAM model suggest that the frequent occurrence of IADFs may be related to the late start of the monsoon rainy season, meaning that we have found a new indicator in tree-ring records that can capture monsoon anomalies. Our results provide further insight into drought variation in the eastern CLP, which also implicates the Asian summer monsoon dynamic.

Climate-driven sapwood-specific hydraulic conductivity and the Huber value but not leaf-specific hydraulic conductivity on a global scale

Efficient water transport is crucial for plant growth and survival. Plant hydraulic conductivity varies between functional groups and biomes and is strongly influenced by changing environmental conditions. However, correlations of conductivity-related hydraulic traits with climatic variables are not fully understood, preventing clarification of plant form and function under climate change scenarios. By compiling leaf-specific hydraulic conductivity (K_L), sapwood-specific hydraulic conductivity (K_s), and Huber values (H_v, sapwood area to leaf area ratio) along with climatic variables including mean annual temperature (MAT), mean annual precipitation (MAP) and aridity index (AI) for 428 species across a wide range of plant functional types (PFTs) and biomes at a global scale, we found greater variability of K_L within PFTs and biomes than across PFTs and biomes. Interaction effects between PFTs and biomes on K_L and K_s were found. The interaction between MAT and MAP played a significant role in K_s and H_v (t = 3.89, P < 0.001 for K_s and t = -5.77, P < 0.001 for H_v). With increasing AI, K_s increased and H_v decreased. K_L was not influenced by the investigated climatic variables. Our study provides a better understanding of the dynamics of hydraulic structure and function across functional groups and biomes and of the abiotic drivers of their large-scale variations.

Retrospective analysis of wood anatomical traits and tree-ring isotopes suggests site-specific mechanisms triggering Araucaria araucana drought-induced dieback

In 2010-2018, Northern Patagonia featured the longest severe drought of the last millennium. This extreme dry spell triggered widespread growth decline and forest dieback. Nonetheless, the roles played by the two major mechanisms driving dieback, hydraulic failure and carbon starvation, are still not clear and understudied in this seasonally dry region. Here, for the 1800-2017 period, we apply a retrospective analysis of radial growth, wood anatomical traits (lumen area, cell-wall thickness) and δ¹³ C and δ¹⁸ O stable isotopes to assess dieback causes of the iconic conifer Araucaria araucana. We selected three stands where declining (defoliated) and nondeclining (not defoliated) trees coexisted along a precipitation gradient from the warm-dry Coastal Range to the cool-wet Andes. At all sites declining trees showed lower radial growth and lower theoretical hydraulic conductivity, suggesting a long-lasting process of hydraulic deterioration in their water transport system compared to nondeclining, coexisting trees. Wood anatomical traits evidenced that this divergence between declining and nondeclining trees started at least seven decades before canopy dieback. In the drier stands, declining trees showed higher water-use efficiency (WUE) throughout the whole period, which we attributed to early stomatal closure, suggesting a greater carbon starvation risk consistent with thinner cell walls. In the wettest stand, we found the opposite pattern. Here, a reduction in WUE coupled with thicker cell walls suggested increased carbon assimilation rates and exposure to drought-induced hydraulic failure. The δ¹⁸ O values indicated different strategies of gas exchange between sites, which are likely a consequence of microsite conditions and water sources. Multiproxy, retrospective quantifications of xylem anatomical traits and tree-ring isotopes provide a robust tool to identify and forecast, which stands or trees will show dieback or, on the contrary, which will likely withstand and be more resilient to future hotter droughts.

Origin of Intra-annual Density Fluctuations in a Semi-arid Area of Northwestern China

Intra-annual density fluctuation (IADF) is a structural modification of the tree ring in response to fluctuations in the weather. The expected changes in monsoon flow would lead to heterogeneous moisture conditions during the growing season and increase the occurrence of IADF in trees of the arid ecosystems of continental Asia. To reveal the timings and physiological mechanisms behind IADF formation, we monitored cambial activity and wood formation in Chinese pine (Pinus tabuliformis) during 2017-2019 at three sites in semi-arid China. We compared the dynamics of xylem formation under a drought event, testing the hypothesis that drought affects the process of cell enlargement and thus induces the production of IADF. Wood microcores collected weekly from April to October were used for anatomical analyses to estimate the timings of cambial activity, and the phases of enlargement, wall thickening, and lignification of the xylem. The first cells started enlargement from late April to early May. The last latewood cells completed differentiation in mid-September. Trees produced IADF in 2018. During that year, a drought in June limited cell production in the cambium, only 36% of the xylem cells being formed in IADF trees, compared to 68% in normal tree rings. IADF cells enlarged under drought in early July and started wall thickening during the rainfall events of late July. The drought restricted cell enlargement and affected wall thickening, resulting in narrow cells with wide walls. Cambium and cell enlargement recovered from the abundant rainfall, producing a new layer with large earlywood tracheids. IADF is a specific adaptation of trees to cope with water deficit events occurring during xylem formation. Our findings confirmed the hypothesis that the June-July drought induces latewood-like IADFs by limiting the process of cell enlargement in the xylem. Our finding suggests a higher occurrence of IADF in trees of arid and semi-arid climates of continental Asia if the changes to monsoon flows result in more frequent drought events during the earlywood formation in June.

Direct and Indirect Effects of Environmental Limitations on White Spruce Xylem Anatomy at Treeline

Treeline ecosystems are of great scientific interest to study the effects of limiting environmental conditions on tree growth. However, tree growth is multidimensional, with complex interactions between height and radial growth. In this study, we aimed to disentangle effects of height and climate on xylem anatomy of white spruce [Picea glauca (Moench) Voss] at three treeline sites in Alaska; i.e., one warm and drought-limited, and two cold, temperature-limited. To analyze general growth differences between trees from different sites, we used data on annual ring width, diameter at breast height (DBH), and tree height. A representative subset of the samples was used to investigate xylem anatomical traits. We then used linear mixed-effects models to estimate the effects of height and climatic variables on our study traits. Our study showed that xylem anatomical traits in white spruce can be directly and indirectly controlled by environmental conditions: hydraulic-related traits seem to be mainly influenced by tree height, especially in the earlywood. Thus, they are indirectly driven by environmental conditions, through the environment's effects on tree height. Traits related to mechanical support show a direct response to environmental conditions, mainly temperature, especially in the latewood. These results highlight the importance of assessing tree growth in a multidimensional way by considering both direct and indirect effects of environmental forcing to better understand the complexity of tree growth responses to the environment.

Growth of 19 conifer species is highly sensitive to winter warming, spring frost and summer drought

BACKGROUND AND AIMS

Conifers are key components of many temperate and boreal forests and are important for forestry, but species differences in stem growth responses to climate are still poorly understood and may hinder effective management of these forests in a warmer and drier future.

METHODS

We studied 19 Northern Hemisphere conifer species planted in a 50-year-old common garden experiment in the Netherlands to (1) assess the effect of temporal dynamics in climate on stem growth, (2) test for a possible positive relationship between the growth potential and climatic growth sensitivity across species, and (3) evaluate the extent to which stem growth is controlled by phylogeny.

KEY RESULTS

Eighty-nine per cent of the species showed a significant reduction in stem growth to summer drought, 37 % responded negatively to spring frost and 32 % responded positively to higher winter temperatures. Species differed largely in their growth sensitivity to climatic variation and showed, for example, a four-fold difference in growth reduction to summer drought. Remarkably, we did not find a positive relationship between productivity and climatic sensitivity, but instead observed that some species combined a low growth sensitivity to summer drought with high growth potential. Both growth sensitivity to climate and growth potential were partly phylogenetically controlled.

CONCLUSIONS

A warmer and drier future climate is likely to reduce the productivity of most conifer species. We did not find a relationship between growth potential and growth sensitivity to climate; instead, some species combined high growth potential with low sensitivity to summer drought. This may help forest managers to select productive species that are able to cope with a warmer and drier future.

Environmental and developmental factors driving xylem anatomy and micro-density in black spruce

Wood density is the product of carbon allocation for structural growth and reflects the trade-off between mechanical support and water conductivity. We tested a conceptual framework based on the assumption that micro-density depends on direct and indirect relationships with endogenous and exogenous factors. The dynamics of wood formation, including timings and rates of cell division, cell enlargement, and secondary wall deposition, were assessed from microcores collected weekly between 2002 and 2016 from five black spruce stands located along a latitudinal gradient in Quebec, Canada. Cell anatomy and micro-density were recorded by anatomical analyses and X-ray measurements. Our structural equation model explained 80% of micro-density variation within the tree-ring with direct effects of wall thickness (σ = 0.61), cell diameter (σ = -0.51), and photoperiod (σ = -0.26). Wood formation dynamics had an indirect effect on micro-density. Micro-density increased under longer periods of cell-wall deposition and shorter durations of enlargement. Our results fill a critical gap in understanding the relationships underlying micro-density variation in conifers. We demonstrated that short-term responses to environmental variations could be overridden by plastic responses that modulate cell differentiation. Our results point to wood formation dynamics as a reliable predictor of carbon allocation in trees.

Role of Mixed-Species Stands in Attenuating the Vulnerability of Boreal Forests to Climate Change and Insect Epidemics

We investigated whether stand species mixture can attenuate the vulnerability of eastern Canada's boreal forests to climate change and insect epidemics. For this, we focused on two dominant boreal species, black spruce [Picea mariana (Mill.) BSP] and trembling aspen (Populus tremuloides Michx.), in stands dominated by black spruce or trembling aspen ("pure stands"), and mixed stands (M) composed of both species within a 36 km² study area in the Nord-du-Québec region. For each species in each stand composition type, we tested climate-growth relations and assessed the impacts on growth by recorded insect epidemics of a black spruce defoliator, the spruce budworm (SBW) [Choristoneura fumiferana (Clem.)], and a trembling aspen defoliator, the forest tent caterpillar (FTC; Malacosoma disstria Hübn.). We implemented linear models in a Bayesian framework to explain baseline and long-term trends in tree growth for each species according to stand composition type and to differentiate the influences of climate and insect epidemics on tree growth. Overall, we found climate vulnerability was lower for black spruce in mixed stands than in pure stands, while trembling aspen was less sensitive to climate than spruce, and aspen did not present differences in responses based on stand mixture. We did not find any reduction of vulnerability for mixed stands to insect epidemics in the host species, but the non-host species in mixed stands could respond positively to epidemics affecting the host species, thus contributing to stabilize ecosystem-scale growth over time. Our findings partially support boreal forest management strategies including stand species mixture to foster forests that are resilient to climate change and insect epidemics.

Contribution of Xylem Anatomy to Tree-Ring Width of Two Larch Species in Permafrost and Non-Permafrost Zones of Siberia

Plants exhibit morphological and anatomical adaptations to cope the environmental constraints of their habitat. How can mechanisms for adapting to contrasting environmental conditions change the patterns of tree rings formation? In this study, we explored differences in climatic conditions of permafrost and non-permafrost zones and assessed their influence on radial growth and wood traits of Larix gmelinii Rupr (Rupr) and Larix sibirica L., respectively. We quantified the contribution of xylem cell anatomy to the tree-ring width variability. Comparison of the anatomical tree-ring parameters over the period 1963–2011 was tested based on non-parametric Mann-Whitney U test. The generalized linear modeling shows the common dependence between TRW and the cell structure characteristics in contrasting environments, which can be defined as non-specific to external conditions. Thus, the relationship between the tree-ring width and the cell production in early- and latewood are assessed as linear, whereas the dependence between the radial cell size in early- and latewood and the tree-ring width becomes significantly non-linear for both habitats. Moreover, contribution of earlywood (EW) and latewood (LW) cells to the variation of TRW (in average 56.8% and 24.4% respectively) was significantly higher than the effect of cell diameters (3.3% (EW) and 17.4% (LW)) for the environments. The results show that different larch species from sites with diverging climatic conditions converge towards similar xylem cell structures and relationships between xylem production and cell traits. The work makes a link between climate and tree-ring structure, and promotes a better understanding the anatomical adaptation of larch species to local environment conditions.

Xylem Anatomical Variability in White Spruce at Treeline Is Largely Driven by Spatial Clustering

The ecological function of boreal forests is challenged by drastically changing climate conditions. Although an increasing number of studies are investigating how climate change is influencing growth and distribution of boreal tree species, there is a lack of studies examining the potential of these species to genetically adapt or phenotypically adjust. Here, we sampled clonally and non-clonally growing white spruce trees (Picea glauca [Moench] Voss) to investigate spatial and genetic effects on tree ring width and on six xylem anatomical traits representing growth, water transport, mechanical support, and wood density. We compared different methods for estimating broad sense heritability (H2) of each trait and we evaluated the effects of spatial grouping and genetic grouping on the xylem anatomical traits with linear models. We found that the three different methods used to estimate H2 were quite robust, showing overall consistent patterns, while our analyses were unsuccessful at fully separating genetic from spatial effects. By evaluating the effect size, we found a significant effect of genetic grouping in latewood density and earlywood hydraulic diameter. However, evaluating model performances showed that spatial grouping was a better predictor than genetic grouping for variance in earlywood density, earlywood hydraulic diameter and growth. For cell wall thickness neither spatial nor genetic grouping was significant. Our findings imply that (1) the variance in the investigated xylem anatomical traits and growth is mainly influenced by spatial clustering (most probably caused by microhabitat conditions), which (2) makes it rather difficult to estimate the heritability of these traits in naturally grown trees in situ. Yet, (3) latewood density and earlywood hydraulic diameter qualified for further analysis on the genetic background of xylem traits and (4) cell wall thickness seems a useful trait to investigate large-scale climatic effects, decoupled from microclimatic, edaphic and genetic influences.

Comparing the Cell Dynamics of Tree-Ring Formation Observed in Microcores and as Predicted by the Vaganov-Shashkin Model

New insights into the intra-annual dynamics of tree-ring formation can improve our understanding of tree-growth response to environmental conditions at high-resolution time scales. Obtaining this information requires, however, a weekly monitoring of wood formation, sampling that is extremely time-intensive and scarcely feasible over vast areas. Estimating the timing of cambial and xylem differentiation by modeling thus represents an interesting alternative for obtaining this important information by other means. Temporal dynamics of cambial divisions can be extracted from the daily tree-ring growth rate computed by the Vaganov-Shashkin (VS) simulation model, assuming that cell production is tightly linked to tree-ring growth. Nonetheless, these predictions have yet to be compared with direct observations of wood development, i.e., via microcoring, over a long time span. We tested the performance of the VS model by comparing the observed and predicted timing of wood formation in black spruce [Picea mariana (Mill.)]. We obtained microcores over 15 years at 5 sites along a latitudinal gradient in Quebec (Canada). The measured variables included cell size and the timing of cell production and differentiation. We calibrated the VS model using daily temperature and precipitation recorded by weather stations located on each site. The predicted and observed timing of cambial and enlarging cells were highly correlated (R ² = 0.8); nonetheless, we detected a systematic overestimation in the predicted timing of cambial cells, with predictions delayed by 1-20 days compared with observations. The growth rate of cell diameter was correlated with the predicted growth rate assigned to each cambial cell, confirming that cell diameter developmental dynamics have the potential to be inferred by the tree-ring growth curve of the VS model. Model performances decrease substantially in estimating the end of wood formation. The systematic errors suggest that the actual relationships implemented in the model are unable to explain the phenological events in autumn. The mismatch between the observed and predicted timing of wood formation in black spruce within our study area can be reduced by better adapting the VS model to wet sites, a context for which this model has been rarely used.