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.

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.

A synthesis of radial growth patterns preceding tree mortality.

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

Separating water-potential induced swelling and shrinking from measured radial stem variations reveals a cambial growth and osmotic concentration signal.

The quantification of cambial growth over short time periods has been hampered by problems to discern between growth and the swelling and shrinking of a tree stem. This paper presents a model, which separates cambial growth and reversible water-potential induced diurnal changes from simultaneously measured whole stem and xylem radial variations, from field-measured Scots pine trees in Finland. The modelled growth, which includes osmotic concentration changes, was compared with (direct) dendrometer measurements and microcore samples. In addition, the relationship of modelled growth and dendrometer measurements to environmental factors was analysed. The results showed that the water-potential induced changes of tree radius were successfully separated from stem growth. Daily growth predicted by the model exhibited a high correlation with the modelled daily changes of osmotic concentration in phloem, and a temperature dependency in early summer. Late-summer growth saw higher dependency on water availability and temperature. Evaluation of the model against dendrometer measurements showed that the latter masked a true environmental signal in stem growth due to water-potential induced changes. The model provides better understanding of radial growth physiology and offers potential to examine growth dynamics and changes due to osmotic concentration, and how the environment affects growth.

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.

CASSIA--a dynamic model for predicting intra-annual sink demand and interannual growth variation in Scots pine.

The control of tree growth vs environment by carbon sources or sinks remains unresolved although it is widely studied. This study investigates growth of tree components and carbon sink-source dynamics at different temporal scales. We constructed a dynamic growth model 'carbon allocation sink source interaction' (CASSIA) that calculates tree-level carbon balance from photosynthesis, respiration, phenology and temperature-driven potential structural growth of tree organs and dynamics of stored nonstructural carbon (NSC) and their modifying influence on growth. With the model, we tested hypotheses that sink demand explains the intra-annual growth dynamics of the meristems, and that the source supply is further needed to explain year-to-year growth variation. The predicted intra-annual dimensional growth of shoots and needles and the number of cells in xylogenesis phases corresponded with measurements, whereas NSC hardly limited the growth, supporting the first hypothesis. Delayed GPP influence on potential growth was necessary for simulating the yearly growth variation, indicating also at least an indirect source limitation. CASSIA combines seasonal growth and carbon balance dynamics with long-term source dynamics affecting growth and thus provides a first step to understanding the complex processes regulating intra- and interannual growth and sink-source dynamics.

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.

Partial asynchrony of coniferous forest carbon sources and sinks at the intra-annual time scale.

As major terrestrial carbon sinks, forests play an important role in mitigating climate change. The relationship between the seasonal uptake of carbon and its allocation to woody biomass remains poorly understood, leaving a significant gap in our capacity to predict carbon sequestration by forests. Here, we compare the intra-annual dynamics of carbon fluxes and wood formation across the Northern hemisphere, from carbon assimilation and the formation of non-structural carbon compounds to their incorporation in woody tissues. We show temporally coupled seasonal peaks of carbon assimilation (GPP) and wood cell differentiation, while the two processes are substantially decoupled during off-peak periods. Peaks of cambial activity occur substantially earlier compared to GPP, suggesting the buffer role of non-structural carbohydrates between the processes of carbon assimilation and allocation to wood. Our findings suggest that high-resolution seasonal data of ecosystem carbon fluxes, wood formation and the associated physiological processes may reduce uncertainties in carbon source-sink relationships at different spatial scales, from stand to ecosystem levels.

A physiological model of softwood cambial growth.

Cambial growth was modelled as a function of detailed levelled physiological processes for cell enlargement and water and sugar transport to the cambium. Cambial growth was described at the cell level where local sugar concentration and turgor pressure induce irreversible cell expansion and cell wall synthesis. It was demonstrated how transpiration and photosynthesis rates, metabolic and physiological processes and structural features of a tree mediate their effects directly on the local water and sugar status and influence cambial growth. Large trees were predicted to be less sensitive to changes in the transient water and sugar status, compared with smaller ones, as they have more water and sugar storage and were, therefore, less coupled to short-term changes in the environment. Modelling the cambial dynamics at the individual cell level turned out to be a complex task as the radial short-distance transport of water and sugars and control signals determining cell division and cessation of cell enlargement and cell wall synthesis had to be described simultaneously.

The effect of artificially induced drought on radial increment and wood properties of Norway spruce.

We studied experimentally the effects of water availability on height and radial increment as well as wood density and tracheid properties of Norway spruce (Picea abies (L.) Karst.). The study was carried out in two long-term N-fertilization experiments in Southern Finland (Heinola and Sahalahti). At each site, one fertilized and one control plot was covered with an under-canopy roof preventing rainwater from reaching the soil. Two uncovered plots were monitored at each site. The drought treatment was initiated in the beginning of growing season and lasted for 60-75 days each year. The treatment was repeated for four to five consecutive years depending on the site. Altogether, 40 sample trees were harvested and discs sampled at breast height. From the discs, ring width and wood density were measured by X-ray densitometry. Tracheid properties were analysed by reflected-light microscopy and image analysis. Reduced soil water potential during the growing season decreased annual radial and height increment and had a small influence on tracheid properties and wood density. No statistically significant differences were found in the average tracheid diameter between the drought-treated and control trees. The average cell wall thickness was somewhat higher (7-10%) for the drought treatment than for the control, but the difference was statistically significant only in Sahalahti. An increased cell wall thickness was found in both early- and latewood tracheids, but the increase was much greater in latewood. In drought-treated trees, cell wall proportion within an annual ring increased, consequently increasing wood density. No interaction between the N fertilization and drought treatment was found in wood density. After the termination of the drought treatment, trees rapidly recovered from the drought stress. According to our results, severe drought due to the predicted climate change may reduce Norway spruce growth but is unlikely to result in large changes in wood properties.

Low growth resilience to drought is related to future mortality risk in trees.

Severe droughts have the potential to reduce forest productivity and trigger tree mortality. Most trees face several drought events during their life and therefore resilience to dry conditions may be crucial to long-term survival. We assessed how growth resilience to severe droughts, including its components resistance and recovery, is related to the ability to survive future droughts by using a tree-ring database of surviving and now-dead trees from 118 sites (22 species, >3,500 trees). We found that, across the variety of regions and species sampled, trees that died during water shortages were less resilient to previous non-lethal droughts, relative to coexisting surviving trees of the same species. In angiosperms, drought-related mortality risk is associated with lower resistance (low capacity to reduce impact of the initial drought), while it is related to reduced recovery (low capacity to attain pre-drought growth rates) in gymnosperms. The different resilience strategies in these two taxonomic groups open new avenues to improve our understanding and prediction of drought-induced mortality.

Volcanic dust veils from sixth century tree-ring isotopes linked to reduced irradiance, primary production and human health.

The large volcanic eruptions of AD 536 and 540 led to climate cooling and contributed to hardships of Late Antiquity societies throughout Eurasia, and triggered a major environmental event in the historical Roman Empire. Our set of stable carbon isotope records from subfossil tree rings demonstrates a strong negative excursion in AD 536 and 541-544. Modern data from these sites show that carbon isotope variations are driven by solar radiation. A model based on sixth century isotopes reconstruct an irradiance anomaly for AD 536 and 541-544 of nearly three standard deviations below the mean value based on modern data. This anomaly can be explained by a volcanic dust veil reducing solar radiation and thus primary production threatening food security over a multitude of years. We offer a hypothesis that persistently low irradiance contributed to remarkably simultaneous outbreaks of famine and Justinianic plague in the eastern Roman Empire with adverse effects on crop production and photosynthesis of the vitamin D in human skin and thus, collectively, human health. Our results provide a hitherto unstudied proxy for exploring the mechanisms of 'volcanic summers' to demonstrate the post-eruption deficiencies in sunlight and to explain the human consequences during such calamity years.

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

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

Reliability of temperature signal in various climate indicators from northern Europe.

We collected relevant observational and measured annual-resolution time series dealing with climate in northern Europe, focusing in Finland. We analysed these series for the reliability of their temperature signal at annual and seasonal resolutions. Importantly, we analysed all of the indicators within the same statistical framework, which allows for their meaningful comparison. In this framework, we employed a cross-validation procedure designed to reduce the adverse effects of estimation bias that may inflate the reliability of various temperature indicators, especially when several indicators are used in a multiple regression model. In our data sets, timing of phenological observations and ice break-up were connected with spring, tree ring characteristics (width, density, carbon isotopic composition) with summer and ice formation with autumn temperatures. Baltic Sea ice extent and the duration of ice cover in different watercourses were good indicators of winter temperatures. Using combinations of various temperature indicator series resulted in reliable temperature signals for each of the four seasons, as well as a reliable annual temperature signal. The results hence demonstrated that we can obtain reliable temperature information over different seasons, using a careful selection of indicators, combining the results with regression analysis, and by determining the reliability of the obtained indicator.

Predicting the decomposition of Scots pine, Norway spruce, and birch stems in Finland.

Models were developed for predicting the decomposition of dead wood for the main tree species in Finland, based on data collected from long-term thinning experiments in southern and central Finland. The decomposition rates were strongly related to the number of years after tree death. In contrast to previous studies, which have used the first-order exponential model, we found that the decomposition rate was not constant. Therefore, the Gompertz and Chapman-Richard's functions were fitted to the data. The slow initial decomposition period was mainly due to the fact that most dead trees remained standing as snags after their death. The initial period was followed by a period of rapid decomposition and, finally, by a period of moderately slow decomposition. Birch stems decomposed more rapidly than Scots pine and Norway spruce stems. Decomposition rates of Norway spruce stems were somewhat lower than those of Scots pine. Because the carbon concentration of decaying boles was relatively stable (about 50%) the rate of carbon loss follows that of mass loss. Models were also developed for the probability that a dead tree remains standing as a snag. During the first years after death, the probability was high. Thereafter, it decreased rapidly, the decrease being faster for birch stems than for Scots pine and Norway spruce stems. Almost all stems had fallen down within 40 years after their death. In Scots pine and Norway spruce, most snags remained hard and belonged to decay class 1. In birch, a higher proportion of snags belonged to the more advanced decay classes. The models provide a framework for predicting dead wood dynamics in managed as well as dense unthinned stands. The models can be incorporated into forest management planning systems, thereby facilitating estimates of carbon dynamics.

Woody biomass production lags stem-girth increase by over one month in coniferous forests.

Wood is the main terrestrial biotic reservoir for long-term carbon sequestration(1), and its formation in trees consumes around 15% of anthropogenic carbon dioxide emissions each year(2). However, the seasonal dynamics of woody biomass production cannot be quantified from eddy covariance or satellite observations. As such, our understanding of this key carbon cycle component, and its sensitivity to climate, remains limited. Here, we present high-resolution cellular based measurements of wood formation dynamics in three coniferous forest sites in northeastern France, performed over a period of 3 years. We show that stem woody biomass production lags behind stem-girth increase by over 1 month. We also analyse more general phenological observations of xylem tissue formation in Northern Hemisphere forests and find similar time lags in boreal, temperate, subalpine and Mediterranean forests. These time lags question the extension of the equivalence between stem size increase and woody biomass production to intra-annual time scales(3, 4, 5, 6). They also suggest that these two growth processes exhibit differential sensitivities to local environmental conditions. Indeed, in the well-watered French sites the seasonal dynamics of stem-girth increase matched the photoperiod cycle, whereas those of woody biomass production closely followed the seasonal course of temperature. We suggest that forecasted changes in the annual cycle of climatic factors(7) may shift the phase timing of stem size increase and woody biomass production in the future.

Seasonal changes in stem radius and production of new tracheids in Norway spruce.

The progress of xylem formation in Norway spruce (Picea abies (L.) Karst.) was measured during one growing season in southern Finland. Stem radius was monitored continuously with band dendrometers, and the formation of new tracheids was determined by examination of small increment cores taken twice weekly. Tracheid production started in June and ceased in August. Xylem formation was fastest in early July, when 0.75-1.25 new tracheids were formed per day. The rate of xylem formation was significantly correlated with mean daily temperature. Synchronous fluctuations in tracheid and lumen diameters were observed at the same relative positions within each annual ring, but no relationship existed between the diameters and weather variables. The timing of changes in stem radius differed from the timing of actual xylem formation. Stem radius increased in April and May, and the fastest daily increments were recorded in June. Increases in stem radius slowed in July, but small increases were measured more than a month after xylem formation had ceased. Daily changes in stem radius were correlated with daily precipitation, reflecting changes in stem water content. Therefore, dendrometers are of dubious value for measuring the timing of actual xylem formation. Small increment cores proved to be useful in assessing actual xylem formation, but the method is laborious.

Effects of nutrient optimization on intra-annual wood formation in Norway spruce.

In the Nordic countries, growth of Norway spruce (Picea abies (L.) Karst.) is generally limited by low availability of nutrients, especially nitrogen. Optimizing forest management requires better insight on how growth responds to the environmental conditions and their manipulation. The aim of this study was to analyse the effects of nutrient optimization on timing and the rate of tracheid formation of Norway spruce and to follow the differentiation of newly formed tracheids. The study was performed during two growing seasons in a long-term nutrient optimization experiment in northern Sweden, where all essential macro- and micronutrients were supplied in irrigation water every second day from mid-June to mid-August. The control plots were without additional nutrients and water. Tracheid formation in the stem was monitored throughout the growing season by weekly sampling of microcores at breast height. The onset of xylogenesis occurred in early June, but in early summer there were no significant between-treatment differences in the onset and relative rate of tracheid formation. In both treatments, the onset of secondary cell wall formation occurred in mid-June. The maximum rate of tracheid formation occurred close to the summer solstice and 50% of the tracheids had been accumulated in early July. Optimized nutrition resulted in the formation of ∼50% more tracheids and delayed the cessation of tracheid formation, which extended the tracheid formation period by 20-50%, compared with control trees. The increased growth was mainly an effect of enhanced tracheid formation rate during the mid- and later-part of the growing season. In the second year, the increased growth rate also resulted in 11% wider tracheids. We conclude that the onset and rate of tracheid formation and differentiation during summer is primarily controlled by photoperiod, temperature and availability of nutrients, rather than supply of carbohydrates.