Do 2 H and 18 O in leaf water reflect environmental drivers differently?

We compiled hydrogen and oxygen stable isotope compositions (δ2 H and δ18 O) of leaf water from multiple biomes to examine variations with environmental drivers. Leaf water δ2 H was more closely correlated with δ2 H of xylem water or atmospheric vapour, whereas leaf water δ18 O was more closely correlated with air relative humidity. This resulted from the larger proportional range for δ2 H of meteoric waters relative to the extent of leaf water evaporative enrichment compared with δ18 O. We next expressed leaf water as isotopic enrichment above xylem water (Δ2 H and Δ18 O) to remove the impact of xylem water isotopic variation. For Δ2 H, leaf water still correlated with atmospheric vapour, whereas Δ18 O showed no such correlation. This was explained by covariance between air relative humidity and the Δ18 O of atmospheric vapour. This is consistent with a previously observed diurnal correlation between air relative humidity and the deuterium excess of atmospheric vapour across a range of ecosystems. We conclude that 2 H and 18 O in leaf water do indeed reflect the balance of environmental drivers differently; our results have implications for understanding isotopic effects associated with water cycling in terrestrial ecosystems and for inferring environmental change from isotopic biomarkers that act as proxies for leaf water.

Historical changes in the stomatal limitation of photosynthesis: empirical support for an optimality principle.

The ratio of leaf internal (ci ) to ambient (ca ) partial pressure of CO2 , defined here as χ, is an index of adjustments in both leaf stomatal conductance and photosynthetic rate to environmental conditions. Measurements and proxies of this ratio can be used to constrain vegetation model uncertainties for predicting terrestrial carbon uptake and water use. We test a theory based on the least-cost optimality hypothesis for modelling historical changes in χ over the 1951-2014 period, across different tree species and environmental conditions, as reconstructed from stable carbon isotopic measurements across a global network of 103 absolutely dated tree-ring chronologies. The theory predicts optimal χ as a function of air temperature, vapour pressure deficit, ca and atmospheric pressure. The theoretical model predicts 39% of the variance in χ values across sites and years, but underestimates the intersite variability in the reconstructed χ trends, resulting in only 8% of the variance in χ trends across years explained by the model. Overall, our results support theoretical predictions that variations in χ are tightly regulated by the four environmental drivers. They also suggest that explicitly accounting for the effects of plant-available soil water and other site-specific characteristics might improve the predictions.

Millennial-scale tree-ring isotope chronologies from coast redwoods provide insights on controls over California hydroclimate variability.

To understand drivers of hydroclimate variability in north-coastal California, we obtained tree cross-sections from eleven coastal redwoods (mean age of 1232 years old) from the northern half of the species range. Tree rings from eight trees were cross-dated and sampled at sub-annual resolution for carbon isotope discrimination (Δ13C) and oxygen isotope composition (δ18O). Tree-ring Δ13C and δ18O, compared to modern climate data, demonstrate these signals primarily record summertime hydroclimate variability-primarily through variables associated with evaporative conditions and/or precipitation. Our 1100-year stable isotope chronologies showed that north-coastal California did not undergo the megadroughts observed elsewhere in California and the western United States. This result implicates extended periods of low winter precipitation, rather than growing season evaporation, as the primary driver of previous megadroughts across California and neighboring regions. Compared to cool conditions prevailing over the Northern Hemisphere during the Little Ice age (1301-1875 of the common era, CE), the frequency of isotopic events of a certain magnitude was greater during periods with warmer Northern Hemisphere temperatures such as the Medieval Climate Anomaly (900-1300 CE) and the modern period (1876 to present). This association between tree-ring isotopic variability and long-term shifts in temperatures is consistent with the expected patterns in mid-latitude hydroclimate variability expected from arctic amplification (i.e., shifts in equator-to-pole temperature differences that modify jet stream speed and amplitude) or amplified quasi-resonant wave activity (i.e., wave-patterns in high-altitude winds that become "trapped" within a certain pattern, thereby producing a longer-duration periods of drought or wetness) across mid-latitudes during the boreal summer.

The enigma of effective path length for (18) O enrichment in leaf water of conifers.

The Péclet correction is often used to predict leaf evaporative enrichment and requires an estimate of effective path length (L). Studies to estimate L in conifer needles have produced unexpected patterns based on Péclet theory and leaf anatomy. We exposed seedlings of six conifer species to different vapour pressure deficits (VPD) in controlled climate chambers to produce steady-state leaf water enrichment (in (18) O). We measured leaf gas exchange, stable oxygen isotopic composition (δ(18) O) of input and plant waters as well as leaf anatomical characteristics. Variation in bulk needle water δ(18) O was strongly related to VPD. Conifer needles had large amounts of water within the vascular strand that was potentially unenriched (up to 40%). Both standard Craig-Gordon and Péclet models failed to accurately predict conifer leaf water δ(18) O without taking into consideration the unenriched water in the vascular strand and variable L. Although L was linearly related to mesophyll thickness, large within-species variation prevented the development of generalizations that could allow a broader use of the Péclet effect in predictive models. Our results point to the importance of within needle water pools and isolating mechanisms that need further investigation in order to integrate Péclet corrections with 'two compartment' leaf water concepts.

Predicting leaf wax n-alkane 2H/1H ratios: controlled water source and humidity experiments with hydroponically grown trees confirm predictions of Craig-Gordon model.

The extent to which both water source and atmospheric humidity affect δ(2)H values of terrestrial plant leaf waxes will affect the interpretations of δ(2)H variation of leaf waxes as a proxy for hydrological conditions. To elucidate the effects of these parameters, we conducted a long-term experiment in which we grew two tree species, Populus fremontii and Betula occidentalis, hydroponically under combinations of six isotopically distinct waters and two different atmospheric humidities. We observed that leaf n-alkane δ(2)H values of both species were linearly related to source water δ(2)H values, but with slope differences associated with differing humidities. When a modified version of the Craig-Gordon model incorporating plant factors was used to predict the δ(2)H values of leaf water, all modelled leaf water values fit the same linear relationship with n-alkane δ(2)H values. These observations suggested a relatively constant biosynthetic fractionation factor between leaf water and n-alkanes. However, our calculations indicated a small difference in the biosynthetic fractionation factor between the two species, consistent with small differences calculated for species in other studies. At present, it remains unclear if these apparent interspecies differences in biosynthetic fractionation reflect species-specific biochemistry or a common biosynthetic fractionation factor with insufficient model parameterization.

Reconstructing relative humidity from plant delta18O and deltaD as deuterium deviations from the global meteoric water line.

Cellulose delta18O and deltaD can provide insights on climates and hydrological cycling in the distant past and how these factors differ spatially. However, most studies of plant cellulose have used only one isotope, most commonly delta18O, resulting in difficulties partitioning variation in delta18O of precipitation vs. evaporative conditions that affect leaf water isotopic enrichment. Moreover, observations of pronounced diurnal differences from conventional steady-state model predictions of leaf water isotopic fractionation have cast some doubt on single isotope modeling approaches for separating precipitation and evaporation drivers of cellulose delta18O or deltaD. We explore a dual isotope approach akin to the concept of deuterium-excess (d), to establish deuterium deviations from the global meteoric water line in leaf water (deltad(l)) as driven by relative humidity (RH). To demonstrate this concept, we survey studies of leaf water delta18O and deltaD in hardwood vs. conifer trees. We then apply the concept to cellulose delta18O and deltaD using a mechanistic model of cellulose delta18O and deltaD to reconstruct deuterium deviations from the global meteoric water line (deltad(c)) in Quercus macrocarpa, Q. robur, and Pseudotsuga menziesii. For each species, deltad(c) showed strong correlations with RH across sites. deltad(c) agreed well with steady-state predictions for Q. macrocarpa, while for Q. robur, the relationship with RH was steeper than expected. The slope of deltad(c) vs. RH of P. menziesii was also close to steady-state predictions, but deltad(c) were more enriched than predicted. This is in agreement with our leaf water survey showing conifer deltad(l) was more enriched than predicted. Our data reveal that applications of this method should be appropriate for reconstructing RH from cellulose delta18O and deltaD after accounting for differences between hardwoods and conifers. Hence, deltad(c) should be useful for understanding variability in RH associated with past climatic cycles, across regional climates, or across complex terrain where climate modeling is challenging. Furthermore, deltad(c) and inferred RH values should help in constraining variation in source water delta18O.

A tree-ring perspective on the terrestrial carbon cycle.

Tree-ring records can provide valuable information to advance our understanding of contemporary terrestrial carbon cycling and to reconstruct key metrics in the decades preceding monitoring data. The growing use of tree rings in carbon-cycle research is being facilitated by increasing recognition of reciprocal benefits among research communities. Yet, basic questions persist regarding what tree rings represent at the ecosystem level, how to optimally integrate them with other data streams, and what related challenges need to be overcome. It is also apparent that considerable unexplored potential exists for tree rings to refine assessments of terrestrial carbon cycling across a range of temporal and spatial domains. Here, we summarize recent advances and highlight promising paths of investigation with respect to (1) growth phenology, (2) forest productivity trends and variability, (3) CO2 fertilization and water-use efficiency, (4) forest disturbances, and (5) comparisons between observational and computational forest productivity estimates. We encourage the integration of tree-ring data: with eddy-covariance measurements to investigate carbon allocation patterns and water-use efficiency; with remotely sensed observations to distinguish the timing of cambial growth and leaf phenology; and with forest inventories to develop continuous, annually-resolved and long-term carbon budgets. In addition, we note the potential of tree-ring records and derivatives thereof to help evaluate the performance of earth system models regarding the simulated magnitude and dynamics of forest carbon uptake, and inform these models about growth responses to (non-)climatic drivers. Such efforts are expected to improve our understanding of forest carbon cycling and place current developments into a long-term perspective.

Tree-ring isotopes adjacent to Lake Superior reveal cold winter anomalies for the Great Lakes region of North America.

Tree-ring carbon isotope discrimination (Δ13C) and oxygen isotopes (δ18O) collected from white pine (Pinus strobus) trees adjacent to Lake Superior show potential to produce the first winter-specific paleoclimate reconstruction with inter-annual resolution for this region. Isotopic signatures from 1976 to 2015 were strongly linked to antecedent winter minimum temperatures (Tmin), Lake Superior peak ice cover, and regional to continental-scale atmospheric winter pressure variability including the North American Dipole. The immense thermal inertia of Lake Superior underlies the unique connection between winter conditions and tree-ring Δ13C and δ18O signals from the following growing season in trees located near the lake. By combining these signals, we demonstrate feasibility to reconstruct variability in Tmin, ice cover, and continental-scale atmospheric circulation patterns (r ≥ 0.65, P < 0.001).

Summer precipitation influences the stable oxygen and carbon isotopic composition of tree-ring cellulose in Pinus ponderosa.

The carbon and oxygen isotopic composition of tree-ring cellulose was examined in ponderosa pine (Pinus ponderosa Dougl.) trees in the western USA to study seasonal patterns of precipitation inputs. Two sites (California and Oregon) had minimal summer rainfall inputs, whereas a third site (Arizona) received as much as 70% of its annual precipitation during the summer months (North American monsoon). For the Arizona site, both the delta(18)O and delta(13)C values of latewood cellulose increased as the fraction of annual precipitation occurring in the summer (July through September) increased. There were no trends in latewood cellulose delta(18)O with the absolute amount of summer rain at any site. The delta(13)C composition of latewood cellulose declined with increasing total water year precipitation for all sites. Years with below-average total precipitation tended to have a higher proportion of their annual water inputs during the summer months. Relative humidity was negatively correlated with latewood cellulose delta(13)C at all sites. Trees at the Arizona site produced latewood cellulose that was significantly more enriched in (18)O compared with trees at the Oregon or California site, implying a greater reliance on an (18)O-enriched water source. Thus, tree-ring records of cellulose delta(18)O and delta(13)C may provide useful proxy information about seasonal precipitation inputs and the variability and intensity of the North American monsoon.

Photosynthetic gas exchange response of poplars to steady-state and dynamic light environments.

The steady-state and dynamic photosynthetic response of two poplar species (Populus tremuloides and P. fremontii) to variations in photon flux density (PFD) were observed with a field portable gas exchange system. These poplars were shown to be very shade intolerant with high light saturation (800 to 1300 µmol photons m-2 s-1) and light compensation (70 to 100 µmol m-2 s-1) points. Understory poplar leaves showed no physiological acclimation to understory light environments. These plants become photosynthetically induced quickly (10 min). Activation of Rubisco was the primary limitation for induction, with stomatal opening playing only a minor role. Leaves maintained high stomatal conductances and stomata were unresponsive to variations in PFD. Leaves were very efficient at utilizing rapidly fluctuating light environments similar to those naturally occurring in canopies. Post-illumination CO2 fixation contributed proportionally more to the carbon gain of leaves during short frequent lightflecks than longer less frequent ones. The benefits of a more dynamic understory light environment for the carbon economy of these species are discussed.

Effect of leaf flutter on the light environment of poplars.

The dynamics of the canopy light environment for two poplar species (Populus tremuloides Michx., and P. fremontii Wats.) were characterized with an array of photocells in fixed positions within the canopy or attached directly to leaves and using a data logger that recorded photon flux density (PFD) at frequencies from 1 to 20 Hz. The majority of sunflecks were short in duration (<1 s) with a similar short interval between sunflecks. Sunflecks contribute as much as 90% of the total daily PFD in the lower canopy. Leaf flutter may cause high frequency (3 to 5 Hz) variations of PFD in poplar canopies. The amount of light intercepted by a fluttering leaf at the top of the canopy decreased with increasing flutter, whereas a fluttering lower canopy leaf showed no such trend. When leaves fluttered at the top of the canopy the understory light environment showed an increased number of shorter sunflecks. Leaf flutter may increase mean PFD for understory leaves. It also creates a canopy light environment that is more dynamic temporally and more evenly distributed spatially. The potential benefits of these changes in light dynamics are discussed.

A controlled test of the dual-isotope approach for the interpretation of stable carbon and oxygen isotope ratio variation in tree rings.

Seedlings of a conifer (Pinus radiata D. Don) and a broad leaf angiosperm (Eucalyptus globulus Labill.) were grown for 100 days in two growth cabinets at 45 or 65% relative humidity. The seedlings were exposed to treatments designed to modify carbon assimilation rates and capacities, stomatal conductance and transpiration to test conceptual models that attempt to clarify the interpretation of carbon isotope discrimination (Δ(13)C) by using oxygen isotope enrichment (Δ(18)O). Differences in relative humidity and within-cabinet treatments (including lower irradiance, lower nitrogen inputs, higher leaf temperature and lower moisture status than control seedlings) produced significant differences in assimilation rates, photosynthetic capacities, stomatal conductance, leaf transpiration rates and leaf evaporative enrichment. The dual-isotope approach accurately interpreted the cause of variation in wood cellulose Δ(13)C for some of the treatments, but not for others. We also tested whether we could use Δ(13)C variation to constrain the interpretation of δ(18)O variation. Carbon isotope discrimination appears to be influenced by transpiration (providing information on leaf evaporative enrichment), but the results did not provide a clear way to interpret such variation. The dual-isotope approach appears to be valid conceptually, but more work is needed to make it operational under different scenarios.

Frost tolerance and ice formation in Pinus radiata needles: ice management by the endodermis and transfusion tissues.

Conifers are among the most frost tolerant tree species. Cryo-scanning electron microscopy (cryo-SEM) was used to visualise ice formation in pine needles to better understand how conifer leaves manage extracellular ice. Acclimated and unacclimated needles of Pinus radiata (D.Don) were subjected to freezing treatments (at a rate of 2°C h-1), tested for electrolyte leakage and sampled for cryo-SEM analysis. Half maximal electrolyte leakage occurred at -4 and -12°C for unacclimated and acclimated needles, respectively. Ice nucleation occurred at similar temperatures (-3°C) in both acclimated and unacclimated pine needles, indicating that frost tolerance did not increase supercooling. During freezing and thawing, the tissues outside and inside the endodermis shrank and swelled independently, with little or no transfer of water between the two regions. During freezing, mesophyll cells shrank, exhibiting cytorrhysis, and extracellular ice accumulated in gas spaces of the mesophyll tissue. Mesophyll cells from acclimated needles recovered their structure after thawing, and unacclimated mesophyll showed significant damage. In the vascular cylinder, ice accumulated in transfusion tracheids which expanded to occupy areas made vacant by shrinkage of transfusion parenchyma, Strasburger cells and the endodermis. This behaviour was reversible in acclimated tissue, and may play an important role in the management of ice during freeze/thaw events.

Comparative analysis of cellulose preparation techniques for use with 13C, 14C, AND 18O isotopic measurements.

A number of operationally defined methods exist for pretreating plant tissues in order to measure C, N, and O isotopes. Because these isotope measurements are used to infer information about environmental conditions that existed at the time of tissue growth, it is important that these pretreatments remove compounds that may have exchanged isotopes or have been synthesized after the original formation of these tissues. In stable isotope studies, many pretreatment methods focus on isolating "cellulose" from the bulk tissue sample because cellulose does not exchange C and O isotopes after original synthesis. We investigated the efficacy of three commonly applied pretreatment methods, the Brendel method and two variants of the Brendel method, the Jayme-Wise method and successive acid/base/acid washes, for use on three tissue types (wood, leaves, roots). We then compared the effect of each method on C and O isotope composition (13C, 14C, 18O), C and N content, and chemical composition of the residue produced (using 13C nuclear magnetic resonance (NMR)). Our results raised concerns over use of the Brendel method as published, as it both added C and N to the sample and left a residue that contains remnant lipids and waxes. Furthermore, this method resulted in 18O values that are enriched relative to the other methods. Modifying the Brendel method by adding a NaOH step (wash) solved many of these problems. We also found that processed residues vary by tissue type. For wood and root tissues, the 13C NMR spectra and the 18O and 13C data showed only small differences between residues for the Jayme-Wise and modified Brendel methods. However, for leaf tissue, 13C NMR data showed that Jayme-Wise pretreatments produced residues that are more chemically similar to cellulose than the other methods. The acid/base/acid washing method generated 13C NMR spectra with incomplete removal of lignin for all tissues tested and both isotopic, and 13C NMR results confirmed that this method should not be used if purified cellulose is desired.

A rapid and precise method for sampling and determining the oxygen isotope ratio of atmospheric water vapor.

A quantitative method for cryogenically sampling atmospheric water vapor on the temporal scale of 10 to 15 min in the field or laboratory is described. The sample apparatus is lightweight, affordable, and easy to assemble. The method allows for H2O:CO2 equilibration within the same sampling tubes and hence increases turnaround time for delta18O analysis. Quantitative analysis in the laboratory showed recovery of a vaporized, known, 18O water standard to 0.2 per thousand precision.

Expressing leaf water and cellulose oxygen isotope ratios as enrichment above source water reveals evidence of a Péclet effect.

There is an increasing ecological interest in understanding the gradients in H(2)(18)O enrichment in leaf water (i.e. a Péclet effect), because an appreciation of the significance of the Péclet effect is important for improving our understanding of the mechanistic processes affecting the (18)O composition of leaf water and plant organic material. In data sets where both source water and leaf water (18)O data are available, we can evaluate the potential contribution of a Péclet effect. As an example, we recalculate data published earlier by Roden and Ehleringer (1999, Oecologia 121:467-477) as enrichments in leaf water (Delta(L)) and cellulose (Delta(cell)) above source water. Based on these recalculations, we present support for the relevance of a Péclet effect in leaves. Further, we demonstrate that the subtle variations in Delta(L) and Delta(cell) caused by a Péclet effect may be masked in experimental systems in which variation in the source water oxygen isotope ratio is considerable.