Uncertainty, Emergence, and Statistics in Dendrochronology

Climate-driven tree growth and mortality in the Black Forest, Germany-Long-term observations

Episodic tree mortality can be caused by various reasons. This study describes climate-driven tree mortality and tree growth in the Black Forest mountain range in Germany. It is based on a 68-year consistent data series describing the annual mortality of all trees growing in a forest area of almost 250 thousand ha. The study excludes mortality caused by storm, snow and ice, and fire. The sequence of the remaining mortality, the so-called "desiccated trees," is analyzed and compared with the sequence of the climatic water balance during the growing season and the annual radial growth of Norway spruce in the Black Forest. The annual radial growth series covers 121 years and the climatic water balance series 140 years. These unique time series enable a quantitative assessment of multidecadal drought and heat impacts on growth and mortality of forest trees on a regional spatial scale. Data compiled here suggest that the mortality of desiccated trees in the Black Forest during the last 68 years is driven by the climatic water balance. Decreasing climatic water balance coincided with an increase in tree mortality and growth decline. Consecutive hot and dry summers enhance mortality and growth decline as a consequence of drought legacies lasting several years. The sensitivity of tree growth and mortality to changes in the climatic water balance increases with the decreasing trend of the climatic water balance. The findings identify the climatic water balance as the main driver of mortality and growth variation during the 68-year observation period on a landscape-scale including a variety of different sites. They suggest that bark beetle population dynamics modify mortality rates. They as well provide evidence that the mortality during the last 140 years never was as high as in the most recent years.

Engineered bacterial swarm patterns as spatial records of environmental inputs

A diverse array of bacteria species naturally self-organize into durable macroscale patterns on solid surfaces via swarming motility-a highly coordinated and rapid movement of bacteria powered by flagella. Engineering swarming is an untapped opportunity to increase the scale and robustness of coordinated synthetic microbial systems. Here we engineer Proteus mirabilis, which natively forms centimeter-scale bullseye swarm patterns, to 'write' external inputs into visible spatial records. Specifically, we engineer tunable expression of swarming-related genes that modify pattern features, and we develop quantitative approaches to decoding. Next, we develop a dual-input system that modulates two swarming-related genes simultaneously, and we separately show that growing colonies can record dynamic environmental changes. We decode the resulting multicondition patterns with deep classification and segmentation models. Finally, we engineer a strain that records the presence of aqueous copper. This work creates an approach for building macroscale bacterial recorders, expanding the framework for engineering emergent microbial behaviors.

Sensitivity and threshold dynamics of Pinus strobus and Quercus spp. in response to experimental and naturally occurring severe droughts

Increased drought frequency and severity are a pervasive global threat, yet the capacity of mesic temperate forests to maintain resilience in response to drought remains poorly understood. We deployed a throughfall removal experiment to simulate a once in a century drought in New Hampshire, USA, which coupled with the region-wide 2016 drought, intensified moisture stress beyond that experienced in the lifetimes of our study trees. To assess the sensitivity and threshold dynamics of two dominant northeastern tree genera (Quercus and Pinus), we monitored sap flux density (Js), leaf water potential and gas exchange, growth and intrinsic water-use efficiency (iWUE) for one pretreatment year (2015) and two treatment years (2016-17). Results showed that Js in pine (Pinus strobus L.) declined abruptly at a soil moisture threshold of 0.15 m3 m-3, whereas oak's (Quercus rubra L. and Quercus velutina Lam.) threshold was 0.11 m3 m-3-a finding consistent with pine's more isohydric strategy. Nevertheless, once oaks' moisture threshold was surpassed, Js declined abruptly, suggesting that while oaks are well adapted to moderate drought, they are highly susceptible to extreme drought. The radial growth reduction in response to the 2016 drought was more than twice as great for pine as for oaks (50 vs 18%, respectively). Despite relatively high precipitation in 2017, the oaks' growth continued to decline (low recovery), whereas pine showed neutral (treatment) or improved (control) growth. The iWUE increased in 2016 for both treatment and control pines, but only in treatment oaks. Notably, pines exhibited a significant linear relationship between iWUE and precipitation across years, whereas the oaks only showed a response during the driest conditions, further underscoring the different sensitivity thresholds for these species. Our results provide new insights into how interactions between temperate forest tree species' contrasting physiologies and soil moisture thresholds influence their responses and resilience to extreme drought.

Global assessment of relationships between climate and tree growth

Tree-ring records provide global high-resolution information on tree-species responses to global change, forest carbon and water dynamics, and past climate variability and extremes. The underlying assumption is a stationary (time-stable), quasi-linear relationship between tree growth and environment, which however conflicts with basic ecological and evolutionary theory. Indeed, our global assessment of the relevant tree-ring literature demonstrates non-stationarity in the majority of tested cases, not limited to specific proxies, environmental parameters, regions or species. Non-stationarity likely represents the general nature of the relationship between tree-growth proxies and environment. Studies assuming stationarity however score two times more citations influencing other fields of science and the science-policy interface. To reconcile ecological reality with the application of tree-ring proxies for climate or environmental estimates, we provide a clarification of the stationarity concept, propose a simple confidence framework for the re-evaluation of existing studies and recommend the use of a new statistical tool to detect non-stationarity in tree-ring proxies. Our contribution is meant to stimulate and facilitate discussion in light of our results to help increase confidence in tree-ring-based climate and environmental estimates for science, the public and policymakers.

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.

Multi-century tree-ring precipitation record reveals increasing frequency of extreme dry events in the upper Blue Nile River catchment

Climate-related environmental and humanitarian crisis are important challenges in the Great Horn of Africa (GHA). In the absence of long-term past climate records in the region, tree-rings are valuable climate proxies, reflecting past climate variations and complementing climate records prior to the instrumental era. We established annually resolved multi-century tree-ring chronology from Juniperus procera trees in northern Ethiopia, the longest series yet for the GHA. The chronology correlates significantly with wet-season (r = .64, p < .01) and annual (r = .68, p < .01) regional rainfall. Reconstructed rainfall since A.D. 1811 revealed significant interannual variations between 2.2 and 3.8 year periodicity, with significant decadal and multidecadal variations during 1855-1900 and 1960-1990. The duration of negative and positive rainfall anomalies varied between 1-7 years and 1-8 years. Approximately 78.4% (95%) of reconstructed dry (extreme dry) and 85.4% (95%) of wet (extreme wet) events lasted for 1 year only and corresponded to historical records of famine and flooding, suggesting that future climate change studies should be both trend and extreme event focused. The average return periods for dry (extreme dry) and wet (extreme wet) events were 4.1 (8.8) years and 4.1 (9.5) years. Extreme-dry conditions during the 19th century were concurrent with drought episodes in equatorial eastern Africa that occurred at the end of the Little Ice Age. El Niño and La Niña events matched with 38.5% and 50% of extreme-dry and extreme-wet events. Equivalent matches for positive and negative Indian Ocean Dipole events were weaker, reaching 23.1 and 25%, respectively. Spatial correlations revealed that reconstructed rainfall represents wet-season rainfall variations over northern Ethiopia and large parts of the Sahel belt. The data presented are useful for backcasting climate and hydrological models and for developing regional strategic plans to manage scarce and contested water resources. Historical perspectives on long-term regional rainfall variability improve the interpretation of recent climate trends.

Assessing Conifer Ray Parenchyma for Ecological Studies: Pitfalls and Guidelines

Ray parenchyma is an essential tissue for tree functioning and survival. This living tissue plays a major role for storage and transport of water, nutrients, and non-structural carbohydrates (NSC), thus regulating xylem hydraulics and growth. However, despite the importance of rays for tree carbon and water relations, methodological challenges hamper knowledge about ray intra- and inter-tree variability and its ecological meaning. In this study we provide a methodological toolbox for soundly quantifying spatial and temporal variability of different ray features. Anatomical ray features were surveyed in different cutting planes (cross-sectional, tangential, and radial) using quantitative image analysis on stem-wood micro-sections sampled from 41 mature Scots pines (Pinus sylvestris). The percentage of ray surface (PERPAR), a proxy for ray volume, was compared among cutting planes and between early- and latewood to assess measurement-induced variability. Different tangential ray metrics were correlated to assess their similarities. The accuracy of cross-sectional and tangential measurements for PERPAR estimates as a function of number of samples and the measured wood surface was assessed using bootstrapping statistical technique. Tangential sections offered the best 3D insight of ray integration into the xylem and provided the most accurate estimates of PERPAR, with 10 samples of 4 mm(2) showing an estimate within ±6.0% of the true mean PERPAR (relative 95% confidence interval, CI95), and 20 samples of 4 mm(2) showing a CI95 of ±4.3%. Cross-sections were most efficient for establishment of time series, and facilitated comparisons with other widely used xylem anatomical features. Earlywood had significantly lower PERPAR (5.77 vs. 6.18%) and marginally fewer initiating rays than latewood. In comparison to tangential sections, PERPAR was systematically overestimated (6.50 vs. 4.92%) and required approximately twice the sample area for similar accuracy. Radial cuttings provided the least accurate PERPAR estimates. This evaluation of ray parenchyma in conifers and the presented guidelines regarding data accuracy as a function of measured wood surface and number of samples represent an important methodological reference for ray quantification, which will ultimately improve the understanding of the fundamental role of ray parenchyma tissue for the performance and survival of trees growing in stressed environments.

Filling the Eastern European gap in millennium-long temperature reconstructions

Tree ring-based temperature reconstructions form the scientific backbone of the current global change debate. Although some European records extend into medieval times, high-resolution, long-term, regional-scale paleoclimatic evidence is missing for the eastern part of the continent. Here we compile 545 samples of living trees and historical timbers from the greater Tatra region to reconstruct interannual to centennial-long variations in Eastern European May-June temperature back to 1040 AD. Recent anthropogenic warming exceeds the range of past natural climate variability. Increased plague outbreaks and political conflicts, as well as decreased settlement activities, coincided with temperature depressions. The Black Death in the mid-14th century, the Thirty Years War in the early 17th century, and the French Invasion of Russia in the early 19th century all occurred during the coldest episodes of the last millennium. A comparison with summer temperature reconstructions from Scandinavia, the Alps, and the Pyrenees emphasizes the seasonal and spatial specificity of our results, questioning those large-scale reconstructions that simply average individual sites.