Temperature-sensitive biochemical 18 O-fractionation and humidity-dependent attenuation factor are needed to predict δ18 O of cellulose from leaf water in a grassland ecosystem

Half of the 18O enrichment of leaf sucrose is conserved in leaf cellulose of a C3 grass across atmospheric humidity and CO2 levels

The 18O enrichment (Δ18O) of cellulose (Δ18OCel) is recognized as a unique archive of past climate and plant function. However, there is still uncertainty regarding the proportion of oxygen in cellulose (pex) that exchanges post-photosynthetically with medium water of cellulose synthesis. Particularly, recent research with C3 grasses demonstrated that the Δ18O of leaf sucrose (Δ18OSuc, the parent substrate for cellulose synthesis) can be much higher than predicted from daytime Δ18O of leaf water (Δ18OLW), which could alter conclusions on photosynthetic versus post-photosynthetic effects on Δ18OCel via pex. Here, we assessed pex in leaves of perennial ryegrass (Lolium perenne) grown at different atmospheric relative humidity (RH) and CO2 levels, by determinations of Δ18OCel in leaves, Δ18OLGDZW (the Δ18O of water in the leaf growth-and-differentiation zone) and both Δ18OSuc and Δ18OLW (adjusted for εbio, the biosynthetic fractionation between water and carbohydrates) as alternative proxies for the substrate for cellulose synthesis. Δ18OLGDZW was always close to irrigation water, and pex was similar (0.53 ± 0.02 SE) across environments when determinations were based on Δ18OSuc. Conversely, pex was erroneously and variably underestimated (range 0.02-0.44) when based on Δ18OLW. The photosynthetic signal fraction in Δ18OCel is much more constant than hitherto assumed, encouraging leaf physiological reconstructions.

Quantifying isotope parameters associated with carbonyl-water oxygen exchange during sucrose translocation in tree phloem

Stable oxygen isotope ratio of tree-ring α-cellulose (δ18Ocel) yields valuable information on many aspects of tree-climate interactions. However, our current understanding of the mechanistic controls on δ18Ocel is incomplete, with a knowledge gap existent regarding the fractionation effect characterizing carbonyl-water oxygen exchange during sucrose translocation from leaf to phloem. To address this insufficiency, we set up an experimental system integrating a vapor 18O-labeling feature to manipulate leaf-level isotopic signatures in tree saplings enclosed within whole-canopy gas-exchange cuvettes. We applied this experimental system to three different tree species to determine their respective relationships between 18O enrichment of sucrose in leaf lamina (Δ18Ol_suc) and petiole phloem (Δ18Ophl_suc) under environmentally/physiologically stable conditions. Based on the determined Δ18Ophl_suc-Δ18Ol_suc relationships, we estimated that on average, at least 25% of the oxygen atoms in sucrose undergo isotopic exchange with water along the leaf-to-phloem translocation path and that the biochemical fractionation factor accounting for such exchange is c. 34‰, markedly higher than the conventionally assumed value of 27‰. Our study represents a significant step toward quantitative elucidation of the oxygen isotope dynamics during sucrose translocation in trees. This has important implications with respect to improving the δ18Ocel model and its related applications in paleoclimatic and ecophysiological contexts.

Modeling the response of Norway spruce tree-ring carbon and oxygen isotopes to selection harvest on a drained peatland forest

Continuous cover forestry (CCF) has gained interest as an alternative to even-aged management particularly on drained peatland forests. However, relatively little is known about the physiological response of suppressed trees when larger trees are removed as a part of CCF practices. Consequently, studies concentrating on process-level modeling of the response of trees to selection harvesting are also rare. Here, we compared, modeled and measured harvest response of previously suppressed Norway spruce (Picea abies) trees to a selection harvest. We quantified the harvest response by collecting Norway spruce tree-ring samples in a drained peatland forest site and measuring the change in stable carbon and oxygen isotopic ratios of wood formed during 2010-20, including five post-harvest years. The measured isotopic ratios were compared with ecosystem-level process model predictions for ${\kern0em }^{13}$C discrimination and ${\kern0em }^{18}$O leaf water enrichment. We found that the model predicted similar but lower harvest response than the measurements. Furthermore, accounting for mesophyll conductance was important for capturing the variation in ${\kern0em }^{13}$C discrimination. In addition, we performed sensitivity analysis on the model, which suggests that the modeled ${\kern0em }^{13}$C discrimination is sensitive to parameters related to CO2 transport through stomata to the mesophyll.

Covariation between oxygen and hydrogen stable isotopes declines along the path from xylem water to wood cellulose across an aridity gradient

Oxygen and hydrogen isotopes of cellulose in plant biology are commonly used to infer environmental conditions, often from time series measurements of tree rings. However, the covariation (or the lack thereof) between δ18 O and δ2 H in plant cellulose is still poorly understood. We compared plant water, and leaf and branch cellulose from dominant tree species across an aridity gradient in Northern Australia, to examine how δ18 O and δ2 H relate to each other and to mean annual precipitation (MAP). We identified a decline in covariation from xylem to leaf water, and onwards from leaf to branch wood cellulose. Covariation in leaf water isotopic enrichment (Δ) was partially preserved in leaf cellulose but not branch wood cellulose. Furthermore, whilst δ2 H was well-correlated between leaf and branch, there was an offset in δ18 O between organs that increased with decreasing MAP. Our findings strongly suggest that postphotosynthetic isotope exchange with water is more apparent for oxygen isotopes, whereas variable kinetic and nonequilibrium isotope effects add complexity to interpreting metabolic-induced δ2 H patterns. Varying oxygen isotope exchange in wood and leaf cellulose must be accounted for when δ18 O is used to reconstruct climatic scenarios. Conversely, comparing δ2 H and δ18 O patterns may reveal environmentally induced shifts in metabolism.

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.

18 O enrichment of sucrose and photosynthetic and nonphotosynthetic leaf water in a C3 grass-atmospheric drivers and physiological relations

The ¹⁸ O enrichment (Δ¹⁸ O) of leaf water affects the Δ¹⁸ O of photosynthetic products such as sucrose, generating an isotopic archive of plant function and past climate. However, uncertainty remains as to whether leaf water compartmentation between photosynthetic and nonphotosynthetic tissue affects the relationship between Δ¹⁸ O of bulk leaf water (Δ¹⁸ O_LW ) and leaf sucrose (Δ¹⁸ O_Sucrose ). We grew Lolium perenne (a C₃ grass) in mesocosm-scale, replicated experiments with daytime relative humidity (50% or 75%) and CO₂ level (200, 400 or 800 μmol mol^-1 ) as factors, and determined Δ¹⁸ O_LW , Δ¹⁸ O_Sucrose and morphophysiological leaf parameters, including transpiration (E_leaf ), stomatal conductance (g_s ) and mesophyll conductance to CO₂ (g_m ). The Δ¹⁸ O of photosynthetic medium water (Δ¹⁸ O_SSW ) was estimated from Δ¹⁸ O_Sucrose and the equilibrium fractionation between water and carbonyl groups (ε_bio ). Δ¹⁸ O_SSW was well predicted by theoretical estimates of leaf water at the evaporative site (Δ¹⁸ O_e ) with adjustments that correlated with gas exchange parameters (g_s or total conductance to CO₂ ). Isotopic mass balance and published work indicated that nonphotosynthetic tissue water was a large fraction (~0.53) of bulk leaf water. Δ¹⁸ O_LW was a poor proxy for Δ¹⁸ O_Sucrose , mainly due to opposite Δ¹⁸ O responses of nonphotosynthetic tissue water (Δ¹⁸ O_non-SSW ) relative to Δ¹⁸ O_SSW , driven by atmospheric conditions.

Capping the hydroxyl groups (-OH) of α-cellulose to reduce Hy-groscopicity for accurate 18O/16O measurement by EA/Py/IRMS

Obtaining an accurate measurement of ¹⁸O/¹⁶O at natural abundance level for land plants-derived α-cellulose with the currently popular EA/Py/IRMS (elemental analysis/pyrolysis/isotope ratio mass spectrometry) method is a challenge due to the hygroscopic nature of the exposed hydroxyl groups, as the ¹⁸O/¹⁶O of adsorbed moisture is usually different from that of the α-cellulose and the relative amount of adsorbed moisture is sample- and relative humidity-dependent. To minimize the hygroscopicity-related measurement error, we capped the hydroxyl groups of α-cellulose by benzylation to various degrees and found that the ¹⁸O/¹⁶O ratio of α-cellulose increased with the degree of benzyl substitution (DS), consistent with the theoretical prediction that a reduced presence of exposed hydroxyl groups should lead to a more accurate (and therefore more reliable) α-cellulose ¹⁸O/¹⁶O measurement. We propose the establishment of a moisture adsorption-degree of substitution or percentage of oxygen-¹⁸O/¹⁶O ratio equation, based on the measurement of C%, O% and δ¹⁸O of variably capped α-cellulose, so that a robust correction can be made in a plant species- and laboratory conditions-specific manner. Failure to do so will lead to an average underestimate of α-cellulose δ¹⁸O by 3.5 mUr under "average" laboratory conditions.

On the Chemical Purity and Oxygen Isotopic Composition of α-Cellulose Extractable from Higher Plants and the Implications for Climate, Metabolic, and Physiological Studies

The ¹⁸O/¹⁶O ratio of α-cellulose in land plants has proved of interest for climate, environmental, physiological, and metabolic studies. Reliable application of such a ratio may be compromised by the presence of hemicellulose impurities in the α-cellulose product obtainable with current extraction methods, as the impurities are known to be isotopically different from that of the α-cellulose. We first compared the quality of hydrolysates of "α-cellulose products" obtained with four representative extraction methods (Jayme and Wise; Brendel; Zhou; Loader) and quantified the hemicellulose-derived non-glucose sugars in the α-cellulose products from 40 land grass species using gas chromatography-mass spectrometry (GC/MS). Second, we performed compound-specific isotope analysis of the hydrolysates using GC/Pyrolysis/IRMS. These results were then compared with the bulk isotope analysis using EA/Pyrolysis/IRMS of the α-cellulose products. We found that overall, the Zhou method afforded the highest purity α-cellulose as judged by the minimal presence of lignin and the second-lowest presence of non-glucose sugars. Isotopic analysis then showed that the O-2-O-6 of the α-cellulose glucosyl units were all depleted in ¹⁸O by 0.0-4.3 mUr (average, 1.9 mUr) in a species-dependent manner relative to the α-cellulose products. The positive isotopic bias of using the α-cellulose product instead of the glucosyl units stems mainly from the fact that the pentoses that dominate hemicellulose contamination in the α-cellulose product are relatively enriched in ¹⁸O (compared to hexoses) as they inherit only the relatively ¹⁸O-enriched O-2-O-5 moiety of sucrose, the common precursor of pentoses and hexoses in cellulose, and are further enriched in ¹⁸O by the (incomplete) hydrolysis.

Explicitly accounting for needle sugar pool size crucial for predicting intra-seasonal dynamics of needle carbohydrates δ18 O and δ13 C

We explore needle sugar isotopic compositions (δ¹⁸ O and δ¹³ C) in boreal Scots pine (Pinus sylvestris) over two growing seasons. A leaf-level dynamic model driven by environmental conditions and based on current understanding of isotope fractionation processes was built to predict δ¹⁸ O and δ¹³ C of two hierarchical needle carbohydrate pools, accounting for the needle sugar pool size and the presence of an invariant pinitol pool. Model results agreed well with observed needle water δ¹⁸ O, δ¹⁸ O and δ¹³ C of needle water-soluble carbohydrates (sugars + pinitol), and needle sugar δ¹³ C (R²  = 0.95, 0.84, 0.60, 0.73, respectively). Relative humidity (RH) and intercellular to ambient CO₂ concentration ratio (C_i /C_a ) were the dominant drivers of δ¹⁸ O and δ¹³ C variability, respectively. However, the variability of needle sugar δ¹⁸ O and δ¹³ C was reduced on diel and intra-seasonal timescales, compared to predictions based on instantaneous RH and C_i /C_a , due to the large needle sugar pool, which caused the signal formation period to vary seasonally from 2 d to more than 5 d. Furthermore, accounting for a temperature-sensitive biochemical ¹⁸ O-fractionation factor and mesophyll resistance in ¹³ C-discrimination were critical. Interpreting leaf-level isotopic signals requires understanding on time integration caused by mixing in the needle sugar pool.

Species variation in the hydrogen isotope composition of leaf cellulose is mostly driven by isotopic variation in leaf sucrose

Experimental approaches to isolate drivers of variation in the carbon-bound hydrogen isotope composition (δ² H) of plant cellulose are rare and current models are limited in their application. This is in part due to a lack in understanding of how ² H-fractionations in carbohydrates differ between species. We analysed, for the first time, the δ² H of leaf sucrose along with the δ² H and δ¹⁸ O of leaf cellulose and leaf and xylem water across seven herbaceous species and a starchless mutant of tobacco. The δ² H of sucrose explained 66% of the δ² H variation in cellulose (R²  = 0.66), which was associated with species differences in the ² H enrichment of sucrose above leaf water ( ε sucrose <math altimg="urn:x-wiley:01407791:media:pce14362:pce14362-math-0001" wiley:location="equation/pce14362-math-0001.png" xmlns="http://www.w3.org/1998/Math/MathML"><mrow><msub><mtext>\unicode{x003B5}</mtext><mtext>sucrose</mtext></msub></mrow></math> : -126% to -192‰) rather than by variation in leaf water δ² H itself. ε sucrose <math altimg="urn:x-wiley:01407791:media:pce14362:pce14362-math-0002" wiley:location="equation/pce14362-math-0002.png" xmlns="http://www.w3.org/1998/Math/MathML"><mrow><msub><mtext>\unicode{x003B5}</mtext><mtext>sucrose</mtext></msub></mrow></math> was positively related to dark respiration (R²  = 0.27), and isotopic exchange of hydrogen in sugars was positively related to the turnover time of carbohydrates (R²  = 0.38), but only when ε sucrose <math altimg="urn:x-wiley:01407791:media:pce14362:pce14362-math-0003" wiley:location="equation/pce14362-math-0003.png" xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mrow><msub><mi mathvariant="normal">\unicode{x003B5}</mi><mtext>sucrose</mtext></msub></mrow></mrow></math> was fixed to the literature accepted value of - 171 <math altimg="urn:x-wiley:01407791:media:pce14362:pce14362-math-0004" wiley:location="equation/pce14362-math-0004.png" xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mrow><mo>\unicode{x02212}</mo><mn>171</mn></mrow></mrow></math> ‰. No relation was found between isotopic exchange of hydrogen and oxygen, suggesting large differences in the processes shaping post-photosynthetic fractionation between elements. Our results strongly advocate that for robust applications of the leaf cellulose hydrogen isotope model, parameterization utilizing δ² H of sugars is needed.

Stomatal conductance limited the CO2 response of grassland in the last century

BACKGROUND

The anthropogenic increase of atmospheric CO2 concentration (ca) is impacting carbon (C), water, and nitrogen (N) cycles in grassland and other terrestrial biomes. Plant canopy stomatal conductance is a key player in these coupled cycles: it is a physiological control of vegetation water use efficiency (the ratio of C gain by photosynthesis to water loss by transpiration), and it responds to photosynthetic activity, which is influenced by vegetation N status. It is unknown if the ca-increase and climate change over the last century have already affected canopy stomatal conductance and its links with C and N processes in grassland.

RESULTS

Here, we assessed two independent proxies of (growing season-integrating canopy-scale) stomatal conductance changes over the last century: trends of δ18O in cellulose (δ18Ocellulose) in archived herbage from a wide range of grassland communities on the Park Grass Experiment at Rothamsted (U.K.) and changes of the ratio of yields to the CO2 concentration gradient between the atmosphere and the leaf internal gas space (ca - ci). The two proxies correlated closely (R2 = 0.70), in agreement with the hypothesis. In addition, the sensitivity of δ18Ocellulose changes to estimated stomatal conductance changes agreed broadly with published sensitivities across a range of contemporary field and controlled environment studies, further supporting the utility of δ18Ocellulose changes for historical reconstruction of stomatal conductance changes at Park Grass. Trends of δ18Ocellulose differed strongly between plots and indicated much greater reductions of stomatal conductance in grass-rich than dicot-rich communities. Reductions of stomatal conductance were connected with reductions of yield trends, nitrogen acquisition, and nitrogen nutrition index. Although all plots were nitrogen-limited or phosphorus- and nitrogen-co-limited to different degrees, long-term reductions of stomatal conductance were largely independent of fertilizer regimes and soil pH, except for nitrogen fertilizer supply which promoted the abundance of grasses.

CONCLUSIONS

Our data indicate that some types of temperate grassland may have attained saturation of C sink activity more than one century ago. Increasing N fertilizer supply may not be an effective climate change mitigation strategy in many grasslands, as it promotes the expansion of grasses at the disadvantage of the more CO2 responsive forbs and N-fixing legumes.