Long-Term Impacts of Defoliator Outbreaks on Larch Xylem Structure and Tree-Ring Biomass

An open-source machine-learning approach for obtaining high-quality quantitative wood anatomy data from E. grandis and P. radiata xylem

Quantitative wood anatomy is a subfield in dendrochronology that requires effective open-source image analysis tools. In this research, the bioimage analysis software QuPath (v0.4.4) is introduced as a candidate for accurately quantifying the cellular properties of the xylem in an automated manner. Additionally, the potential of QuPath to detect the transition of early- to latewood tracheids over the growing season was evaluated to assess a potential application in dendroecological studies. Various algorithms in QuPath were optimized to quantify different xylem cell types in Eucalyptus grandis and the transition of early- to latewood tracheids in Pinus radiata. These algorithms were coded into cell detection scripts for automatic quantification of stem microsections and compared to a manually curated method to assess the accuracy of the cell detections. The automatic cell detection approach, using QuPath, has been validated to be reproducible with an acceptable error when assessing fibers, vessels, early- and latewood tracheids. However, further optimization for parenchyma is still required. This proposed method developed in QuPath provides a scalable and accurate approach for quantifying anatomical features in stem microsections. With minor amendments to the detection and classification algorithms, this strategy is likely to be viable in other plant species.

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.

The admixture of Quercus sp. in Pinus sylvestris stands influences wood anatomical trait responses to climatic variability and drought events

Introduction

Forests are threatened by increasingly severe and more frequent drought events worldwide. Mono-specific forests, developed as a consequence of widespread management practices established early last century, seem particularly susceptible to global warming and drought compared with mixed-species forests. Although, in several contexts, mixed-species forests display higher species diversity, higher productivity, and higher resilience, previous studies highlighted contrasting findings, with not only many positive but also neutral or negative effects on tree performance that could be related to tree species diversity. Processes underlying this relationship need to be investigated. Wood anatomical traits are informative proxies of tree functioning, and they can potentially provide novel long-term insights in this regard. However, wood anatomical traits are critically understudied in such a context. Here, we assess the role of tree admixture on Pinus sylvestris L. xylem traits such as mean hydraulic diameter, cell wall thickness, and anatomical wood density, and we test the variability of these traits in response to climatic parameters such as temperature, precipitation, and drought event frequency and intensity.

Methods

Three monocultural plots of P. sylvestris and three mixed-stand plots of P. sylvestris and Quercus sp. were identified in Poland and Spain, representing Continental and Mediterranean climate types, respectively. In each plot, we analyzed xylem traits from three P. sylvestris trees, for a total of nine trees in monocultures and nine in mixed stands per study location.

Results

The results highlighted that anatomical wood density was one of the most sensitive traits to detect tree responses to climatic conditions and drought under different climate and forest types. Inter-specific facilitation mechanisms were detected in the admixture between P. sylvestris and Quercus sp., especially during the early growing season and during stressful events such as spring droughts, although they had negligible effects in the late growing season.

Discussion

Our findings suggest that the admixture between P. sylvestris and Quercus sp. increases the resilience of P. sylvestris to extreme droughts. In a global warming scenario, this admixture could represent a useful adaptive management option.

Tree-ring isotopes from the Swiss Alps reveal non-climatic fingerprints of cyclic insect population outbreaks over the past 700 years

Recent experiments have underlined the potential of δ2H in tree-ring cellulose as a physiological indicator of shifts in autotrophic versus heterotrophic processes (i.e., the use of fresh versus stored non-structural carbohydrates). However, the impact of these processes has not yet been quantified under natural conditions. Defoliator outbreaks disrupt tree functioning and carbon assimilation, stimulating remobilization, therefore providing a unique opportunity to improve our understanding of changes in δ2H. By exploring a 700-year tree-ring isotope chronology from Switzerland, we assessed the impact of 79 larch budmoth (LBM, Zeiraphera griseana [Hübner]) outbreaks on the growth of its host tree species, Larix decidua [Mill]. The LBM outbreaks significantly altered the tree-ring isotopic signature, creating a 2H-enrichment and an 18O- and 13C-depletion. Changes in tree physiological functioning in outbreak years are shown by the decoupling of δ2H and δ18O (O-H relationship), in contrast to the positive correlation in non-outbreak years. Across the centuries, the O-H relationship in outbreak years was not significantly affected by temperature, indicating that non-climatic physiological processes dominate over climate in determining δ2H. We conclude that the combination of these isotopic parameters can serve as a metric for assessing changes in physiological mechanisms over time.

The effects of intrinsic water-use efficiency and climate on wood anatomy

Climate warming may induce growth decline in warm-temperate areas subjected to seasonal soil moisture deficit, whereas increasing atmospheric CO₂ concentration (C_a) is expected to enhance tree growth. An accurate understanding of tree growth and physiological processes responding to climate warming and increasing C_a is critical. Here, we analyzed tree-ring stable carbon isotope and wood anatomical traits of Pinus tabuliformis from Qinling Mountains in China to understand how lumen diameter (LD) determining potential hydraulic conductivity and cell-wall thickness (CWT) determining carbon storage responded to climate and C_a. The effects of climate and C_a on intrinsic water-use efficiency (iWUE) were isolated, and iWUE values due to only-climate (iWUE_Clim) and only-CO₂ effects (iWUE_CO2) were obtained. During a low-iWUE period, the influences of climate on earlywood (EW) LD and latewood (LW) CWT prevailed. During a high-iWUE period, CO₂ fertilization promoted cell enlargement and carbon storage but this was counteracted by a negative influence of climate warming. The limiting direct effects of iWUE_Clim and indirect effects of climate on EW LD were greater than on LW CWT. P. tabuliformis in temperate forests will face a decline of growth and carbon fixation, but will produce embolism-resistant tracheids with narrow lumen responding to future hotter droughts.

Tree-ring anatomy of Pinus cembra trees opens new avenues for climate reconstructions in the European Alps

Tree rings form the backbone of high-resolution palaeoclimatology and represent one of the most frequently used proxy to reconstruct climate variability of the Common Era. In the European Alps, reconstructions were often based on tree-ring width (TRW) and maximum latewood density (MXD) series, with a focus on European larch. By contrast, only a very limited number of dendroclimatic studies exists for long-lived, multi-centennial Pinus cembra, despite the widespread occurrence of the species at treeline sites across the European Alps. This lack of reconstructions can be ascribed to the difficulties encountered in past studies in extracting a robust climate signal from TRW and MXD chronologies. In this study, we tested various wood anatomical parameters from P. cembra as proxies for the reconstruction of past air temperatures. To this end, we measured anatomical cell parameters and TRW of old-growth trees from the God da Tamangur forest stand, known for being the highest pure, and continuous P. cembra forest in Europe. We demonstrate that several wood anatomical parameters allow robust reconstruction of past temperature variability at annual to multidecadal timescales. Best results are obtained with maximum latewood radial cell wall thickness (CWTrad) measured at 40 μm radial band width. Over the 1920-2017 period, the CWTrad chronology explains 62 % and >80 % of interannual and decadal variability of air temperatures during a time window corresponding roughly with the growing season. These values exceed those found in past work on P. cembra and even exceed the values reported for MXD chronologies built with L. decidua and hitherto considered the gold standard for dendroclimatic reconstructions in the European Alps. The wood anatomical analysis of P. cembra records therefore unveils a dormant potential and opens new avenues for a species that has been considered unsuitable for climate reconstructions so far.

Insights into source/sink controls on wood formation and photosynthesis from a stem chilling experiment in mature red maple

Whether sources or sinks control wood growth remains debated with a paucity of evidence from mature trees in natural settings. Here, we altered carbon supply rate in stems of mature red maples (Acer rubrum) within the growing season by restricting phloem transport using stem chilling; thereby increasing carbon supply above and decreasing carbon supply below the restrictions, respectively. Chilling successfully altered nonstructural carbon (NSC) concentrations in the phloem without detectable repercussions on bulk NSC in stems and roots. Ring width responded strongly to local variations in carbon supply with up to seven-fold differences along the stem of chilled trees; however, concurrent changes in the structural carbon were inconclusive at high carbon supply due to large local variability of wood growth. Above chilling-induced bottlenecks, we also observed higher leaf NSC concentrations, reduced photosynthetic capacity, and earlier leaf coloration and fall. Our results indicate that the cambial sink is affected by carbon supply, but within-tree feedbacks can downregulate source activity, when carbon supply exceeds demand. Such feedbacks have only been hypothesized in mature trees. Consequently, these findings constitute an important advance in understanding source-sink dynamics, suggesting that mature red maples operate close to both source- and sink-limitation in the early growing season.

Mask, Train, Repeat! Artificial Intelligence for Quantitative Wood Anatomy

The recent developments in artificial intelligence have the potential to facilitate new research methods in ecology. Especially Deep Convolutional Neural Networks (DCNNs) have been shown to outperform other approaches in automatic image analyses. Here we apply a DCNN to facilitate quantitative wood anatomical (QWA) analyses, where the main challenges reside in the detection of a high number of cells, in the intrinsic variability of wood anatomical features, and in the sample quality. To properly classify and interpret features within the images, DCNNs need to undergo a training stage. We performed the training with images from transversal wood anatomical sections, together with manually created optimal outputs of the target cell areas. The target species included an example for the most common wood anatomical structures: four conifer species; a diffuse-porous species, black alder (Alnus glutinosa L.); a diffuse to semi-diffuse-porous species, European beech (Fagus sylvatica L.); and a ring-porous species, sessile oak (Quercus petraea Liebl.). The DCNN was created in Python with Pytorch, and relies on a Mask-RCNN architecture. The developed algorithm detects and segments cells, and provides information on the measurement accuracy. To evaluate the performance of this tool we compared our Mask-RCNN outputs with U-Net, a model architecture employed in a similar study, and with ROXAS, a program based on traditional image analysis techniques. First, we evaluated how many target cells were correctly recognized. Next, we assessed the cell measurement accuracy by evaluating the number of pixels that were correctly assigned to each target cell. Overall, the "learning process" defining artificial intelligence plays a key role in overcoming the issues that are usually manually solved in QWA analyses. Mask-RCNN is the model that better detects which are the features characterizing a target cell when these issues occur. In general, U-Net did not attain the other algorithms' performance, while ROXAS performed best for conifers, and Mask-RCNN showed the highest accuracy in detecting target cells and segmenting lumen areas of angiosperms. Our research demonstrates that future software tools for QWA analyses would greatly benefit from using DCNNs, saving time during the analysis phase, and providing a flexible approach that allows model retraining.

Manipulating phloem transport affects wood formation but not local nonstructural carbon reserves in an evergreen conifer

How variations in carbon supply affect wood formation remains poorly understood in particular in mature forest trees. To elucidate how carbon supply affects carbon allocation and wood formation, we attempted to manipulate carbon supply to the cambial region by phloem girdling and compression during the mid- and late-growing season and measured effects on structural development, CO₂ efflux and nonstructural carbon reserves in stems of mature white pines. Wood formation and stem CO₂ efflux varied with a location relative to treatment (i.e., above or below the restriction). We observed up to twice as many tracheids formed above versus below the treatment after the phloem transport manipulation, whereas the cell-wall area decreased only slightly below the treatments, and cell size did not change relative to the control. Nonstructural carbon reserves in the xylem, needles and roots were largely unaffected by the treatments. Our results suggest that low and high carbon supply affects wood formation, primarily through a strong effect on cell proliferation, and respiration, but local nonstructural carbon concentrations appear to be maintained homeostatically. This contrasts with reports of decoupling of source activity and wood formation at the whole-tree or ecosystem level, highlighting the need to better understand organ-specific responses, within-tree feedbacks, as well as phenological and ontogenetic effects on sink-source dynamics.

Triggering Bimodal Radial Stem Growth in Pinus sylvestris at a Drought-Prone Site by Manipulating Stem Carbon Availability

A bimodal radial growth (RG) pattern, i.e., growth peaks in spring and autumn, was repeatedly found in trees in the Mediterranean regions, where summer drought causes reduction or cessation of cambial activity. In a dry inner Alpine valley of the Eastern Alps (Tyrol, Austria, 750 m asl), Pinus sylvestris shows unimodal RG with onset and cessation of cambial activity in early April and late June, respectively. A resumption of cambial activity after intense summer rainfall was not observed in this region. In a field experiment, we tested the hypothesis that early cessation of cambial activity at this drought-prone site is an adaptation to limited water availability leading to an early and irreversible switch of carbon (C) allocation to belowground. To accomplish this, the C status of young P. sylvestris trees was manipulated by physical blockage of phloem transport (girdling) 6 weeks after cessation of cambial cell division. Influence of manipulated C availability on RG was recorded by stem dendrometers, which were mounted above the girdling zone. In response to blockage of phloem flow, resumption of cambial activity was detected above girdling after about 2 weeks. Although the experimentally induced second growth surge lasted for the same period as in spring (c. 2 months), the increment was more than twice as large due to doubling of daily maximum RG rate. After girdling, wood anatomical traits above girdling no longer showed any significant differences between earlywood and latewood tracheids indicating pronounced effects of C availability on cell differentiation. Below girdling, no reactivation of cambial activity occurred, but cell wall thickness of last formed latewood cell was reduced due to lack of C supply after girdling. Intense RG resumption after girdling indicates that cessation of cambial activity can be reversed by manipulating C status of the stem. Hence, our girdling study yielded strong support for the hypothesis that belowground organs exert high C sink strengths on the drought-prone study site. Furthermore, this work highlights the need of in-depth experimental studies in order to understand the interactions between endogenous and exogenous factors on cambial activity and xylem cell differentiation more clearly.

Phenotypic plasticity of European larch radial growth and wood density along a-1,000 m elevational gradient

Phenotypic plasticity is a key mechanism for sedentary long-living species to adjust to changing environment. Here, we use mature Larix decidua tree-ring variables collected along an elevational transect in the French Alps to characterize the range of individual plastic responses to temperature. Stem cores from 821 mature Larix decidua trees have been collected from four plots distributed along a 1,000-m elevational gradient in a natural forest to build up individual linear reaction norms of tree-ring microdensity traits to temperature. The sign, magnitude and spread of variations of the slopes of the individual reaction norms were used to characterize variation of phenotypic plasticity among plots and traits. Results showed a large range of phenotypic plasticity (with positive and negative slopes) at each elevational plot and for each tree-ring variable. Overall, phenotypic plasticity tends to be larger but positive at higher elevation, negative at the warmer lower sites, and more variable in the center of the elevation distribution. Individual inter-ring reaction norm is a valuable tool to retrospectively characterize phenotypic plasticity of mature forest trees. This approach applied to Larix decidua tree-ring micro-density traits along an elevation gradient showed the existence of large inter-individual variations that could support local adaptation to a fast-changing climate.

Tree physiological monitoring of the 2018 larch budmoth outbreak: preference for leaf recovery and carbon storage over stem wood formation in Larix decidua

Insect defoliation impacts forest productivity worldwide, highlighting the relevance of plant-insect interactions. The larch budmoth (Zeiraphera griseana Hübner) is one of the most extensively studied defoliators, where numerous tree ring-based analyses on its host (Larix decidua Mill.) have aided in identifying outbreak dynamics over the past millennia. Yet, outbreaks have been widely absent after the early 1980s, and little is known about the in situ tree physiological responses and the allocation of carbon resources during and after defoliation. In summer 2018, we tracked an ongoing larch budmoth outbreak in a well-studied larch forest in the Swiss Alps. We performed biweekly monitoring on an affected and unaffected site using a unique combination of xylogenesis observations, measurements of non-structural carbohydrates, isotopic analysis of needle assimilates and ground-based and remote-sensed leaf trait observations. The budmoth induced a defoliation that lasted 40 days and could be detected by satellite observations. Soluble sugars significantly decreased in needles and stem phloem of the defoliated trees, while starch levels remained stable in the stem and root xylem compared to the control. Carbon and oxygen isotope ratios in needle assimilates indicated that neither photosynthetic assimilation rates nor stomatal conductance was different between sites before, during and after the outbreak. Defoliated trees ceased cell wall thickening 17 days earlier than unaffected trees, showing the earliest halt of ring formation recorded from 2007 untill 2013 and causing significant thinner cell walls, particularly in the latewood. No significant differences were found for cell enlargement rates and ring width. Our study revealed that an outbreak causes a downregulation of cell wall thickening first, while no starch is mobilized or leaf physiology is adjusted to compensate for the reduced carbon source due to defoliation. Our observations suggest that affected larch trees prioritize leaf recovery and carbon storage over wood biomass development.