Cold-season freeze frequency is a pervasive driver of subcontinental forest growth.

As northern latitudes experience rapid winter warming, there is an urgent need to assess the effect of varying winter conditions on tree growth and forest carbon sequestration potential. We examined tree growth responses to variability in cold-season (November­April) frequency of freeze days (FFD) over 1951 to 2018 using tree-ring data from 35,217 trees and 57 species at 4,375 sites distributed across Canada. We found that annual radial growth responses to FFD varied by species, with some commonalities across genera and clades. The growth of gymnosperms with late spring leaf-out strategies was negatively related to FFD; years with high FFD were most detrimental to the annual growth of Pinus banksiana, Pinus contorta, Larix lyalli, Abies amabilis, and Abies lasiocarpa. In contrast, the growth of angiosperms with early leaf-out strategies, namely, Populus tremuloides and Betula papyrifera, was better in the coldest years, and gymnosperms with intermediate leaf-out timing, such as widespread Picea mariana and Picea glauca, had no consistent relationship to FFD. Tree growth responses to FFD were further modulated by tree size, tree age, regional climate (i.e., mean cold-season temperature), and local site conditions. Overall, our results suggest that moderately warming winters may temporarily improve the growth of widespread pines and some high-elevation conifers in western Canada, whereas warming winters may be detrimental to the growth of widespread boreal angiosperms. Our findings also highlight the value of using species-specific climate-growth relationships to refine predictions of forest carbon dynamics.

Asymmetric effects of hydroclimate extremes on eastern US tree growth: Implications on current demographic shifts and climate variability.

Forests around the world are experiencing changes due to climate variability and human land use. How these changes interact and influence the vulnerability of forests are not well understood. In the eastern United States, well-documented anthropogenic disturbances and land-use decisions, such as logging and fire suppression, have influenced forest species assemblages, leading to a demographic shift from forests dominated by xeric species to those dominated by mesic species. Contemporarily, the climate has changed and is expected to continue to warm and produce higher evaporative demand, imposing stronger drought stress on forest communities. Here, we use an extensive network of tree-ring records from common hardwood species across ~100 sites and ~1300 trees in the eastern United States to examine the magnitude of growth response to both wet and dry climate extremes. We find that growth reductions during drought exceed the positive growth response to pluvials. Mesic species such as Liriodendron tulipifera and Acer saccharum, which are becoming more dominant, are more sensitive to drought than more xeric species, such as oaks (Quercus) and hickory (Carya), especially at moderate and extreme drought intensities. Although more extreme droughts produce a larger annual growth reduction, mild droughts resulted in the largest cumulative growth decreases due to their higher frequency. When using global climate model projections, all scenarios show drought frequency increasing substantially (3-9 times more likely) by 2100. Thus, the ongoing demographic shift toward more mesic species in the eastern United States combined with drier conditions results in larger drought-induced growth declines, suggesting that drought will have an even larger impact on aboveground carbon uptake in the future in the eastern United States.

Recent increases in tropical cyclone precipitation extremes over the US east coast.

The impacts of inland flooding caused by tropical cyclones (TCs), including loss of life, infrastructure disruption, and alteration of natural landscapes, have increased over recent decades. While these impacts are well documented, changes in TC precipitation extremes-the proximate cause of such inland flooding-have been more difficult to detect. Here, we present a latewood tree-ring-based record of seasonal (June 1 through October 15) TC precipitation sums (ΣTCP) from the region in North America that receives the most ΣTCP: coastal North and South Carolina. Our 319-y-long ΣTCP reconstruction reveals that ΣTCP extremes (≥0.95 quantile) have increased by 2 to 4 mm/decade since 1700 CE, with most of the increase occurring in the last 60 y. Consistent with the hypothesis that TCs are moving slower under anthropogenic climate change, we show that seasonal ΣTCP along the US East Coast are positively related to seasonal average TC duration and TC translation speed.

Rapid increases in shrubland and forest intrinsic water-use efficiency during an ongoing megadrought.

Globally, intrinsic water-use efficiency (iWUE) has risen dramatically over the past century in concert with increases in atmospheric CO2 concentration. This increase could be further accelerated by long-term drought events, such as the ongoing multidecadal "megadrought" in the American Southwest. However, direct measurements of iWUE in this region are rare and largely constrained to trees, which may bias estimates of iWUE trends toward more mesic, high elevation areas and neglect the responses of other key plant functional types such as shrubs that are dominant across much of the region. Here, we found evidence that iWUE is increasing in the Southwest at one of the fastest rates documented due to the recent drying trend. These increases were particularly large across three common shrub species, which had a greater iWUE sensitivity to aridity than Pinus ponderosa, a common tree species in the western United States. The sensitivity of both shrub and tree iWUE to variability in atmospheric aridity exceeded their sensitivity to increasing atmospheric [CO2]. The shift to more water-efficient vegetation would be, all else being equal, a net positive for plant health. However, ongoing trends toward lower plant density, diminished growth, and increasing vegetation mortality across the Southwest indicate that this increase in iWUE is unlikely to offset the negative impacts of aridification.

Long-term decrease in Asian monsoon rainfall and abrupt climate change events over the past 6,700 years.

Asian summer monsoon (ASM) variability and its long-term ecological and societal impacts extending back to Neolithic times are poorly understood due to a lack of high-resolution climate proxy data. Here, we present a precisely dated and well-calibrated tree-ring stable isotope chronology from the Tibetan Plateau with 1- to 5-y resolution that reflects high- to low-frequency ASM variability from 4680 BCE to 2011 CE. Superimposed on a persistent drying trend since the mid-Holocene, a rapid decrease in moisture availability between ∼2000 and ∼1500 BCE caused a dry hydroclimatic regime from ∼1675 to ∼1185 BCE, with mean precipitation estimated at 42 ± 4% and 5 ± 2% lower than during the mid-Holocene and the instrumental period, respectively. This second-millennium-BCE megadrought marks the mid-to late Holocene transition, during which regional forests declined and enhanced aeolian activity affected northern Chinese ecosystems. We argue that this abrupt aridification starting ∼2000 BCE contributed to the shift of Neolithic cultures in northern China and likely triggered human migration and societal transformation.

Global tree intrinsic water use efficiency is enhanced by increased atmospheric CO2 and modulated by climate and plant functional types.

We conducted a meta-analysis of carbon and oxygen isotopes from tree ring chronologies representing 34 species across 10 biomes to better understand the environmental drivers and physiological mechanisms leading to historical changes in tree intrinsic water use efficiency (iWUE), or the ratio of net photosynthesis (Anet) to stomatal conductance (gs), over the last century. We show a ∼40% increase in tree iWUE globally since 1901, coinciding with a ∼34% increase in atmospheric CO2 (Ca), although mean iWUE, and the rates of increase, varied across biomes and leaf and wood functional types. While Ca was a dominant environmental driver of iWUE, the effects of increasing Ca were modulated either positively or negatively by climate, including vapor pressure deficit (VPD), temperature, and precipitation, and by leaf and wood functional types. A dual carbon-oxygen isotope approach revealed that increases in Anet dominated the observed increased iWUE in ∼83% of examined cases, supporting recent reports of global increases in Anet, whereas reductions in gs occurred in the remaining ∼17%. This meta-analysis provides a strong process-based framework for predicting changes in tree carbon gain and water loss across biomes and across wood and leaf functional types, and the interactions between Ca and other environmental factors have important implications for the coupled carbon-hydrologic cycles under future climate. Our results furthermore challenge the idea of widespread reductions in gs as the major driver of increasing tree iWUE and will better inform Earth system models regarding the role of trees in the global carbon and water cycles.

Decadal-scale variability and warming affect spring timing and forest growth across the western Great Lakes region.

The Great Lakes region of North America has warmed by 1-2 °C on average since pre-industrial times, with the most pronounced changes observable during winter and spring. Interannual variability in temperatures remains high, however, due to the influence of ocean-atmosphere circulation patterns that modulate the warming trend across years. Variations in spring temperatures determine growing season length and plant phenology, with implications for whole ecosystem function. Studying how both internal climate variability and the "secular" warming trend interact to produce trends in temperature is necessary to estimate potential ecological responses to future warming scenarios. This study examines how external anthropogenic forcing and decadal-scale variability influence spring temperatures across the western Great Lakes region and estimates the sensitivity of regional forests to temperature using long-term growth records from tree-rings and satellite data. Using a modeling approach designed to test for regime shifts in dynamic time series, this work shows that mid-continent spring climatology was strongly influenced by the 1976/1977 phase change in North Pacific atmospheric circulation, and that regional forests show a strengthening response to spring temperatures during the last half-century.

Physiological meaning of bimodal tree growth-climate response patterns.

Correlation coefficients are widely used to identify and quantify climate signals in proxy archives. Significant relationships between tree-ring chronologies and meteorological measurements are typically applied by dendroclimatologists to distinguish between more or less relevant climate variation for ring formation. While insignificant growth-climate correlations are usually found with cold season months, we argue that weak relationships with high summer temperatures not necessarily disprove their importance for xylogenesis. Here, we use maximum latewood density records from ten treeline sites between northern Scandinavia and southern Spain to demonstrate how monthly growth-climate correlations change from narrow unimodal to wide bimodal seasons when vegetation periods become longer and warmer. Statistically meaningful relationships occur when minimum temperatures exceed 'biological zero' at around 5° C. We conclude that the absence of evidence for statistical significance between tree growth and the warmest summer temperatures at Mediterranean sites is no evidence of absence for the physiological importance of high summer temperatures for ring formation. Since correlation should never be confused with causation, statistical values require mechanistic understanding, and different interpretations are needed for insignificant correlations within and outside the growing season.

Responses of stem growth and canopy greenness of temperate conifers to dry spells.

Dry spells strongly influence biomass production in forest ecosystems. Their effects may last several years following a drought event, prolonging growth reduction and therefore restricting carbon sequestration. Yet, our understanding of the impact of dry spells on the vitality of trees' above-ground biomass components (e.g., stems and leaves) at a landscape level remains limited. We analyzed the responses of Pinus sylvestris and Picea abies to the four most severe drought years in topographically complex sites. To represent stem growth and canopy greenness, we used chronologies of tree-ring width and time series of the Normalized Difference Vegetation Index (NDVI). We analyzed the responses of radial tree growth and NDVI to dry spells using superposed epoch analysis and further explored this relationship using mixed-effect models. Our results show a stronger and more persistent response of radial growth to dry spells and faster recovery of canopy greenness. Canopy greenness started to recover the year after the dry spell, whereas radial tree growth remained reduced for the two subsequent years and did not recover the pre-drought level until the fourth year after the event. Stem growth and canopy greenness were influenced by climatic conditions during and after drought events, while the effect of topography was marginal. The opposite responses of stem growth and canopy greenness following drought events suggest a different impact of dry spells on trees´ sink and source compartments. These results underscore the crucial importance of understanding the complexities of tree growth as a major sink of atmospheric carbon.

Tree-ring based summer temperature variability since 1790 CE in the Hindu Kush region of northern Pakistan.

The Hindu Kush high-altitude regions of Pakistan are currently experiencing severe consequences as a result of global warming. In this sense, increasing soil erosion and the quick melting of glaciers are two particularly evident effects. In such a scenario, understanding long-term temperature changes is crucial for making accurate forecasts about how the Hindu Kush region may experience regional temperature changes in the future. In this study, the climate tree-ring width (TRW) analysis designated a positive and significant correlation (r = 0.622, p < 0.001) between the TRW chronology and the June to September (summer) mean maximum temperature (MMT). Using the tree-ring width of Pinus wallichiana A. B. Jackson, we reconstructed summer temperatures in the Hindu Kush region from 1790 CE. Statistical analysis showed that the reconstruction model has explained 38.7% of the climate variance during the instrumental period of 1967 to 2018 CE. Five extremely warm summer periods (≥ 4 years; before the instrumental period 1967-2018 CE) of 1804-1830, 1839-1862, 1876-1879, 1905-1910, 1923-1935 CE, and six cold summer periods of 1790-1803, 1832-1838, 1863-1875, 1880-1904, 1911-1922, and 1936-1945 CE have been observed during the past 229 years. Individually, the year 1856 CE experienced severe warmth (31.85 °C), whereas 1794 CE was relatively cooler (29.60 °C). The spectral multi-taper method (MTM) shows significant (p < 0.05) cycles, which take place about every 9.3, 5.7, 4.2, and 3.6 years. In particular, the 9.3-year cycle, which closely aligns with the 11-year solar activity cycle, suggests a potential correlation between solar activity and local temperature fluctuations. Moreover, our reconstruction demonstrates a significant degree of consistency when compared to actual climate data and regional temperature reconstruction series, reporting a strong logic of trust in the reliability and accuracy of our findings. This evidence reaffirms that our reconstruction shows significant and dependable regional temperature signals, notably being representative for the Hindu Kush region.

Hydrogen isotope fractionation in carbohydrates of leaves and xylem tissues follows distinct phylogenetic patterns: a common garden experiment with 73 tree and shrub species.

Recent methodological advancements in determining the nonexchangeable hydrogen isotopic composition (δ2 Hne ) of plant carbohydrates make it possible to disentangle the drivers of hydrogen isotope (2 H) fractionation processes in plants. Here, we investigated the influence of phylogeny on the δ2 Hne of twig xylem cellulose and xylem water, as well as leaf sugars and leaf water, across 73 Northern Hemisphere tree and shrub species growing in a common garden. 2 H fractionation in plant carbohydrates followed distinct phylogenetic patterns, with phylogeny reflected more in the δ2 Hne of leaf sugars than in that of twig xylem cellulose. Phylogeny had no detectable influence on the δ2 Hne of twig or leaf water, showing that biochemistry, not isotopic differences in plant water, caused the observed phylogenetic pattern in carbohydrates. Angiosperms were more 2 H-enriched than gymnosperms, but substantial δ2 Hne variations also occurred at the order, family, and species levels within both clades. Differences in the strength of the phylogenetic signals in δ2 Hne of leaf sugars and twig xylem cellulose suggest that the original phylogenetic signal of autotrophic processes was altered by subsequent species-specific metabolism. Our results will help improve 2 H fractionation models for plant carbohydrates and have important consequences for dendrochronological and ecophysiological studies.

Unenriched xylem water contribution during cellulose synthesis influenced by atmospheric demand governs the intra-annual tree-ring δ18 O signature.

The oxygen isotope composition (δ18 O) of tree-ring cellulose is used to evaluate tree physiological responses to climate, but their interpretation is still limited due to the complexity of the isotope fractionation pathways. We assessed the relative contribution of seasonal needle and xylem water δ18 O variations to the intra-annual tree-ring cellulose δ18 O signature of larch trees at two sites with contrasting soil water availability in the Swiss Alps. We combined biweekly δ18 O measurements of soil water, needle water, and twig xylem water with intra-annual δ18 O measurements of tree-ring cellulose, xylogenesis analysis, and mechanistic and structural equation modeling. Intra-annual cellulose δ18 O values resembled source water δ18 O mean levels better than needle water δ18 O. Large parts of the rings were formed under high proportional exchange with unenriched xylem water (pex ). Maximum pex values were achieved in August and imprinted on sections at 50-75% of the ring. High pex values were associated with periods of high atmospheric evaporative demand (VPD). While VPD governed needle water δ18 O variability, we estimated a limited Péclet effect at both sites. Due to a variable pex , source water has a strong influence over large parts of the intra-annual tree-ring cellulose δ18 O variations, potentially masking signals coming from needle-level processes.

The North American Monsoon buffers forests against the ongoing megadrought in the Southwestern United States.

The US Southwest has been entrenched in a two-decade-long megadrought (MD), the most severe since 800 CE, which threatens the long-term vitality and persistence of regional montane forests. Here, we report that in the face of record low winter precipitation and increasing atmospheric aridity, seasonal activity of the North American Monsoon (NAM) climate system brings sufficient precipitation during the height of the summer to alleviate extreme tree water stress. We studied seasonally resolved, tree-ring stable carbon isotope ratios across a 57-year time series (1960-2017) in 17 Ponderosa pine forests distributed across the NAM geographic domain. Our study focused on the isotope dynamics of latewood (LW), which is produced in association with NAM rains. During the MD, populations growing within the core region of the NAM operated at lower intrinsic and higher evaporative water-use efficiencies (WUEi and WUEE , respectively), compared to populations growing in the periphery of the NAM domain, indicating less physiological water stress in those populations with access to NAM moisture. The disparities in water-use efficiencies in periphery populations are due to a higher atmospheric vapor pressure deficit (VPD) and reduced access to summer soil moisture. The buffering advantage of the NAM, however, is weakening. We observed that since the MD, the relationship between WUEi and WUEE in forests within the core NAM domain is shifting toward a drought response similar to forests on the periphery of the NAM. After correcting for past increases in the atmospheric CO2 concentration, we were able to isolate the LW time-series responses to climate alone. This showed that the shift in the relation between WUEi and WUEE was driven by the extreme increases in MD-associated VPD, with little advantageous influence on stomatal conductance from increases in atmospheric CO2 concentration.

The contribution of Indigenous stewardship to an historical mixed-severity fire regime in British Columbia, Canada.

Indigenous land stewardship and mixed-severity fire regimes both promote landscape heterogeneity, and the relationship between them is an emerging area of research. In our study, we reconstructed the historical fire regime of Ne Sextsine, a 5900-ha dry, Douglas fir-dominated forest in the traditional territory of the T'exelc (Williams Lake First Nation) in British Columbia, Canada. Between 1550 and 1982 CE, we found median fire intervals of 18 years at the plot level and 4 years at the study-site level. Ne Sextsine was characterized by an historical mixed-severity fire regime, dominated by frequent, low-severity fires as indicated by fire scars, with infrequent, mixed-severity fires indicated by cohorts. Differentiating low- from mixed-severity plots over time was key to understanding the drivers of the fire regime at Ne Sextsine. Low-severity plots were coincident with areas of highest use by the T'exelc, including winter village sites, summer fishing camps, and travel corridors. The high fire frequency in low-severity plots ceased in the 1870s, following the smallpox epidemic, the forced relocation of Indigenous peoples into small reserves, and the prohibition of Indigenous burning. In contrast, the mixed-severity plots were coincident with areas where forest resources, such as deer or certain berry species, were important. The high fire frequency in the mixed-severity plots continued until the 1920s when industrial-scale grazing and logging began, facilitated by the establishment of a nearby railway. T'exelc oral histories and archeological evidence at Ne Sextsine speak to varied land stewardship, reflected in the spatiotemporal complexity of low- and mixed-severity fire plots. Across Ne Sextsine, 63% of cohorts established and persisted in the absence of fire after colonial impacts beginning in the 1860s, resulting in a dense, homogeneous landscape that no longer supports T'exelc values and is more likely to burn at uncharacteristic high severities. This nuanced understanding of the Indigenous contribution to a mixed-severity fire regime is critical for advancing proactive fire mitigation that is ecoculturally relevant and guided by Indigenous expertise.

Securing timelines in the ancient Mediterranean using multiproxy annual tree-ring data.

Calendar-dated tree-ring sequences offer an unparalleled resource for high-resolution paleoenvironmental reconstruction. Where such records exist for a few limited geographic regions over the last 8,000 to 12,000 years, they have proved invaluable for creating precise and accurate timelines for past human and environmental interactions. To expand such records across new geographic territory or extend data for certain regions further backward in time, new applications must be developed to secure "floating" (not yet absolutely dated) tree-ring sequences, which cannot be assigned single-calendar year dates by standard dendrochronological techniques. This study develops two approaches to this problem for a critical floating tree-ring chronology from the East Mediterranean Bronze-Iron Age. The chronology is more closely fixed in time using annually resolved patterns of 14C, modulated by cosmic radiation, between 1700 and 1480 BC. This placement is then tested using an anticorrelation between calendar-dated tree-ring growth responses to climatically effective volcanism in North American bristlecone pine and the Mediterranean trees. Examination of the newly dated Mediterranean tree-ring sequence between 1630 and 1500 BC using X-ray fluorescence revealed an unusual calcium anomaly around 1560 BC. While requiring further replication and analysis, this anomaly merits exploration as a potential marker for the eruption of Thera.

Recent anthropogenic curtailing of Yellow River runoff and sediment load is unprecedented over the past 500 y.

The Yellow River (YR) is the fifth-longest and the most sediment-laden river in the world. Frequent historical YR flooding events, however, have resulted in tremendous loss of life and property, whereas in recent decades YR runoff and sediment load have fallen sharply. To put these recent changes in a longer-term context, we reconstructed natural runoff for the middle reach of the YR back to 1492 CE using a network of 31 moisture-sensitive tree-ring width chronologies. Prior to anthropogenic interference that started in the 1960s, the lowest natural runoff over the past 500 y occurred during 1926 to 1932 CE, a drought period that can serve as a benchmark for future planning of YR water allocation. Since the late 1980s, the low observed YR runoff has exceeded the natural range of runoff variability, a consequence of the combination of decreasing precipitation and increasing water consumption by direct and indirect human activities, particularly agricultural irrigation. This reduced runoff has resulted in an estimated 58% reduction of the sediment load in the upper reach of the YR and 29% reduction in the middle reach.

The Western Himalayan fir tree ring record of soil moisture in Pakistan since 1855.

Data on historical soil moisture is crucial for assessing and responding to droughts that commonly occur in climate change-affected countries. The Himalayan temperate forests in Pakistan are particularly at risk of climate change. Developing nations lack the means to gather surface soil moisture (SSM) information. Tree rings are one way to bridge this gap. Here, we employed dendrochronological methods on climate-sensitive tree rings from Abies pindrow to reconstruct the SSM in the Western Himalayan mountain region of Pakistan from 1855 to 2020. December (r = 0.41), May (r = 0.40), and June (r = 0.65) SSMs were found to be the limiting factors for A. pindrow growth. However, only the June SSM showed reconstruction possibility (coefficient of efficiency = 0.201 and reduction of error = 0.325). Over the studied period, we found 6 years (wet year) when June SSM was above the threshold of 32.04 (mean + 2 δ) and 1 year (dry year) when June SSM was below the threshold of 21.28 (mean - 2 δ). It was revealed that 1921 and 1917 were the driest and wettest SSM of all time, with means of 19.34 and 36.49, respectively. Our study shows that winter soil moisture is critical for the growing season in the context of climate change. Climate change has broad impacts on tree growth in the Western Himalayas. This study will assist various stakeholders in understanding and managing local and regional climate change.

Distribution and pollution characteristics of organophosphate esters: reflected by tree rings of arbor species.

Organophosphate esters (OPEs) are emerging pollutants. Currently, research on OPEs in tree rings is still limited. In this study, tree rings of five arbor species from Sichuan Province, China, were sampled to study the occurrence and distribution of six OPEs, which were quantitatively analyzed by gas chromatography-mass spectrometry (GC-MS). The total concentrations of OPEs in all samples ranged from 189.79 (Fir species) to 341.23 ng/g (Toona sinensis), with average concentration of 284.77 ± 46.66 ng/g. So, arbor could be used as good passive samplers for OPEs. The levels of OPEs among five arbor species showed no significant difference (p = 0.668 > 0.05), suggesting that the pollution status of OPEs in a region or country could be roughly assessed by any arbor tree species. In this study area, tris(2-butoxyethyl) phosphate (TBEP) was the dominant OPEs followed by tri(2-chloroethyl) phosphate (TCEP). Tris(2-ethylhexyl) phosphate (TEHP) and tri-n-butyl phosphate (TnBP) showed relatively stable concentrations in each arbor species, while the other four OPEs including TBEP, triphenyl phosphate (TPhP), tri(chloropropyl) phosphate (TCPP) and TCEP had significantly different concentrations. Interestingly, the absorption and accumulation of OPEs by tree rings of arbor species were quite different from that of inorganic elements reported by other studies.

Local hydrological gradients structure high intraspecific variability in plant hydraulic traits in two dominant central Amazonian tree species.

Addressing the intraspecific variability of functional traits helps understand how climate change might influence the distribution of organismal traits across environments, but this is notably understudied in the Amazon, especially for plant hydraulic traits commonly used to project drought responses. We quantified the intraspecific trait variability of leaf mass per area, wood density, and xylem embolism resistance for two dominant central Amazonian tree species, along gradients of water and light availability, while accounting for tree age and height. Intraspecific variability in hydraulic traits was high, with within-species variability comparable to the whole-community variation. Hydraulic trait variation was modulated mostly by the hydrological environment, with higher embolism resistance of trees growing on deep-water-table plateaus compared with shallow-water-table valleys. Intraspecific variability of leaf mass per area and wood density was mostly modulated by intrinsic factors and light. The different environmental and intrinsic drivers of variation among and within individuals lead to an uncoupled coordination among carbon acquisition/conservation and water-use traits. Our findings suggest multivariate ecological strategies driving tropical tree distributions even within species, and reflect differential within-population sensitivities along environmental gradients. Therefore, intraspecific trait variability must be considered for accurate predictions of the responses of tropical forests to climate change.

isocalcR: An R package to streamline and standardize stable isotope calculations in ecological research.

The use of stable isotopes to characterize ecosystem dynamics and infer leaf gas exchange processes has become increasingly prevalent over the last few decades within the ecological community. While advancements in theory and our understanding of the physiological processes controlling isotopic signatures in plants has been well-documented, no standardized tool currently exists to facilitate the computation of common isotope-derived plant physiological indices. Here, we present isocalcR, an R package intended to facilitate the use of stable isotope data from plant tissues by providing an integrated collection of functions and recommended reference data. The isocalcR R package contains a suite of functions that compute leaf carbon isotope discrimination (∆13 C), leaf intercellular [CO2 ], the ratio of leaf intercellular to atmospheric [CO2 ], the difference between atmospheric and leaf intercellular [CO2 ], and intrinsic water use efficiency from carbon isotope signatures in leaf or wood tissue with minimal inputs from the user. isocalcR also implements and provides recommended input atmospheric [CO2 ] (ppm) and atmospheric δ13 CO2 (‰) data for the period 0-2021 C.E. A major goal of isocalcR is to provide a standardized, open-source tool to streamline the calculation of reproducible physiological indices from stable isotope signatures in plant tissues, incorporating the most up-to-date theory, while simultaneously eliminating potential errors associated with complex calculations. isocalcR can be used for any location globally as long as the user provides information regarding temperature and elevation to the main workhorse functions.