Larch Cellulose Shows Significantly Depleted Hydrogen Isotope Values With Respect to Evergreen Conifers in Contrast to Oxygen and Carbon Isotopes

Abstract

The analysis of the stable isotope of the tree-ring cellulose is an important tool for paleo climatic investigations. Long tree-ring chronologies consist predominantly of oaks and conifers in Europe, including larch trees (Larix decidua) and cembran pines (Pinus cembra) that form very long tree ring chronologies in the Alps and grow at the treeline, where tree growth is mainly determined by temperature variations. We analyzed δ¹³C, δ¹⁸O and δ²H isotopes in the cellulose extracted from tree-rings of wood samples collected at high altitude in the Swiss and Tyrol Alps, covering the whole Holocene period. We found that larch cellulose was remarkably more depleted in deuterium than that of cembran pine, with mean δ²H values of −113.4 ± 9.7‰ for larch and of −65.4 ± 11.3‰ for cembran pine. To verify if these depleted values were specific to larch or a property of the deciduous conifers, we extended the analysis to samples from various living conifer species collected at the Bern Botanical Garden. The results showed that not only the larch, but also all the samples of the deciduous larch family had a cellulose composition that was highly depleted in δ²H with regard to the other evergreen conifers including cembran pine, a difference that we attribute to a faster metabolism of the deciduous conifers. The δ¹⁸O values were not statistically different among the species, in agreement with the hypothesis that they are primary signals of the source water. While the δ¹³C values were slightly more depleted for larch than for cembran pine, likely due to metabolic differences of the two species. We conclude that the deciduous larch conifers have specific metabolic hydrogen fractionations and that the larch unique signature of δ²H is useful to recognize it from other conifers in subfossil wood samples collected for paleoclimatic studies. For climate information the absolute δ²H values of larch should be considered carefully and separate from other species.

Introduction

The stable isotope ratios of carbon, oxygen and hydrogen have been studied in various components of plants, including bulk wood samples, lignin, whole leaves and leaf waxes (Borella et al., 1998; Borella et al., 1999; Loader et al., 2003; Kahmen et al., 2011; Kimak and Leuenberger 2015), and those analyzed in tree ring samples were successfully used for the reconstruction of past climate conditions, thus widening the field of dendroclimatology (Borella et al., 1998; Leuenberger 1998; McCarroll and Loader 2004; Treydte et al., 2007; Frank et al., 2015). Tree ring samples are made of cellulose, hemicelluloses and lignin (Khezami et al., 2005) and cellulose is the most used material for isotope ratio studies (Borella et al., 1998; Borella et al., 1999; Loader et al., 2003; McCarroll and Loader 2004). It is a carbohydrate that contains the stable isotope ratios of carbon (δ¹³C), oxygen (δ¹⁸O) and hydrogen (δD alias δ²H) (Cormier et al., 2018) that carry climate information and can act as climate proxies. In dry alpine sites the δ¹³C values are controlled mainly by stomatal conductance, which is linked to summer moisture stress and thus to precipitations (Gagen et al., 2004). Early studies noted strong relationships between cellulose δ¹⁸O and mean annual temperature and humidity (Libby et al., 1976; Burk and Stuiver 1981). More recently in relatively humid locations tree ring δ¹⁸O values have been used to reconstruct precipitation (Brienen et al., 2012; Boysen et al., 2014), whereas in other locations tree-ring δ¹⁸O reflected a combination of vapor pressure deficit, relative humidity and sunshine (Roden and Ehleringer 2007; Xu et al., 2011; Boysen et al., 2014; Hartl-Meier et al., 2014; Labuhn et al., 2014). The δ²H values have been measured in tree rings for paleo climatological investigations for forty years (Schiegl 1974), but their interpretation still remains complex. The interest in δ²H is stimulated by recent developments in isotope ratio mass spectrometry combined with equilibrium methods that allow the measurement of carbon, oxygen and hydrogen isotope ratios at the same time (Filot et al., 2006; Loader et al., 2014). This led to results that indicated that δ²H in plant cellulose is determined by 1) the δ²H value of the water source, 2) the water evaporation in the leaf that enriches heavier water isotopes in the liquid phase and 3) the biosynthetic isotopic fractionation between leaf water and the final organic compounds like cellulose, that includes many complex biochemical processes (Cormier et al., 2019). Altogether, the deuterium fractionation is related to both environmental and physiological factors (Augusti 2007).

The oxygen and hydrogen in the plant cellulose originate mainly from meteoric water and the fractionation processes during pre-photosynthetic, photosynthetic and post-carboxylation reactions affects both δ18O and δ²H, but in a different way (Yakir et al., 1990). Leaf evapotranspiration is expected to enrich heavy isotopes in leaf water (Nabeshima et al., 2018), but some enzymes of glucose biosynthesis are selective for the stable H isotopes and imprint a specific Deuterium (²H) abundance in each H group bound to the C of cellulose (Augusti et al., 2006). It should be noted that the analytical approach of cellulose used here (the cellulose pyrolysis) detects all hydrogens, but the exchangeable hydrogens are equilibrated with a water of known hydrogen isotope ratio (Filot et al., 2006). The non-exchangeable hydrogen isotope ratio is then calculated from the measured and the equilibrated non-exchangeable hydrogens. The specific δ²H in each H bound to the C of cellulose documents a distinct biochemical history (Augusti 2007) that showed minor differences between spruce and oak (Augusti et al., 2008). However, the isotopomer variations seem to have a minor effect on total variation of δ²H in cellulose, and no specific studies on larch isotopomers have been reported so far.

In the last decade various studies have been published on deuterium abundance in cellulose from different tree genus, mainly from Quercus, Fagaceae (Kimak et al., 2015), Picea (Gori et al., 2013), Pinus, Pinales (Nabeshima et al., 2018), Larix (Hafner et al., 2011) and Abeis that all use C3 photosynthesis mechanism (Tang et al., 2000; McCarroll and Loader 2004; Szczepanek et al., 2006).

One study found that in the early annual growing season the cellulose of leaves of deciduous oak and beech was enriched (−20‰), while that of the evergreen taxus was more depleted (−100‰) (Kimak et al., 2015). Similar results were found in tree rings cellulose from oak and beech, with higher δ²H values (−10‰) in the early annual growing season than in the late growing season (−50‰), possibly due to the use of stored carbohydrates (Nabeshima et al., 2018). Thus, δ²H values change during the tree growth phases. The heterotrophic metabolism causes a 2H (D) an enrichment of cellulose higher than that of autotropic metabolism, probably due to a higher recycling of compounds in the Calvin and tricarboxylic acid cycles that is associated with a more complete exchange of C-bound H with the surrounding and enriched foliar water (Ziegler 1995; Nabeshima et al., 2018; Cormier et al., 2019). The effect of metabolism on ²H fractionation in cellulose is supported also by the finding that C4 plants are more ²H -enriched than C3 plants (Ziegler et al., 1976). The different photosynthesis systems and anatomies of C3 and C4 plants influence ²H -biological fractionation via C-bound H exchanges with water in the different anatomical compartments (Zhou et al., 2016).

The data reported above indicate that factors other than climate affect the isotopic fractionation signal of tree ring wood, and they should be identified for their robust use in paleo dendroclimatology. To this aim we have measured δ²H, δ¹⁸O and δ¹³C isotopes in the cellulose of 4399 subfossil samples of 120 conifers cembran pine and larch and of one spruce collected at 12 different sites in Switzerland, Austria and Italy. They cover 8930 years and were collected in the project the Eastern Alpine Conifer Chronology for paleoclimatic studies (Nicolussi et al., 2009). With the aim to recognize the non-climatic signals in the dataset we investigated, in this study, site influences and species differences regarding isotope compositions. These investigations were extended to many living conifers sampled at the Bern Botanical Garden.

Materials and Methods

Sampling Sites in the Alps

The wood samples belong to the Eastern Alpine Conifer Chronology Dataset, the collection and preparation of which have been described earlier (Nicolussi et al., 2009). Part of the dataset includes samples from living trees that cover the last 190 years, from 1820–2010 A.D. These tree ring samples were collected at three upper tree-line sites at the proximity of glacier sites at elevations above 1900 m (Ziehmer et al., 2018). These three sites are: Val Roseg (VRR) where samples of four cembran pines and four larch trees were collected, valley d’Hérens (FPCR) four larches and at Grimselsee (UZAR) four cembran pine trees were sampled (Table 1). The other part of the dataset consists of subfossil wood samples from the dominant tree-line species, the deciduous European larch (Larix decidua Mill) and the evergreen cembran pine (Pinus cembra (L.)) with the inclusion of a single tree of the evergreen spruce (Picea abies (L.) H.Karst). They cover the last 9000 years, the whole Holocene, and were collected at 12 different sites in the European Alps which cover a SW-NE transect, at an elevation range of 1,930 to 2,400 m a.s.l. with different aspects. The site locations are shown in Figure 1 and their characteristics are listed in Table 1. Each species and site are represented by four trees. The samples we analyzed were wood blocks of tree-rings covering five consecutive years (Ziehmer et al., 2018).

FIGURE 1

Wood Samples From Bern Botanical Garden

We sampled tree branches at the Botanical Garden in Bern, some of which were dead. The bark and the cambium were removed from the wood before preparing the samples for analysis. A small wood amount (about 50 mg) was chosen for the extraction of the cellulose and the isotope ratio analysis was done as described below. Each sample was measured at least twice. The samples included 20 different conifer species and one of the broad leaf gymnosperm Ginkgo biloba. The larch family was represented by six different species: Larix decidua, L. decidua subsp. polonica, L. sibirica, L. laricina and one species of Pseudolarix amabilis. The evergreen conifers included eight different species of pinus (Pinus bungeana, P. nigra, P. sylvestris, P. cembra, P. mugo, P. speuce, P. ponderosa, P. wallichiana). The other evergreen conifers were Cephalotaxus harringtonii, Cedrus atlantica, Cupressus nootkatensis, Taxus baccata and Tsuga Canadensis.

Sample Wood Dating

Calendar-dated tree-ring width series of Holocene wood sections are available at the Institute of Geography of the University of Innsbruck, where also the Eastern Alpine Conifer Chronology (EACC) has been established (Nicolussi et al., 2009). For isotope and cellulose analyses, the calendar-dated discs were cut into five-year blocks.

Stable Isotope Analysis

The isotope measurements were carried out using wood blocks of consecutive five tree rings. The separation of wood samples into five-year blocks was performed at the University of Innsbruck, Austria. The procedure of cellulose extraction and the calculation of the cellulose content (cellulose dry weight/ wood dry weight) (Ziehmer et al., 2018) as well as the triple-isotope analysis were described before (Loader et al., 2015). Briefly, we used conventional Isotope Ratio Mass Spectrometry (Isoprime 100) coupled to a pyrolysis unit (HEKAtech GmbH, Germany), which is similar to the previously used TC/EA (for technical details see (Leuenberger 2007)). This approach was extended to measurements of non-exchangeable hydrogen of alpha-cellulose using the on-line equilibration method (Filot et al., 2006; Loader et al., 2015). The results are reported in per mil (‰) relative to the Vienna Pee Dee Belemnite (VPDB) for carbon and to Vienna Standard Mean Ocean Water (VSMOW) for hydrogen and oxygen (Coplen 1994). We applied two-point calibrations using Merck cellulose and IAEA-CH-6 crystalline sugar for carbon and oxygen isotopes as well as Merck cellulose and IAEA-CH-7 PE polyethylene foil for hydrogen isotopes. Wei Ming 101 cellulose was used as target material. Merck and Wei Ming 101 cellulose are well characterized internal laboratory standard materials. The precision of the measurement was ±3.0‰ for hydrogen, ±0.3‰ for oxygen and ±0.15‰ for carbon (Loader et al., 2015).

Statistical Analysis

Statistical analysis was performed using the R software and was aimed at verifying differences between the groups. The type of data distribution was analyzed with the Shapiro-test (Shapiro and Francia 1972). To verify the equality between the groups we used pairwise comparisons by the Wilcoxon test (Wilcoxon et al., 1970). The density of the data is shown by the violin contours in Figure 2 (Hintze and Nelson 1998).

FIGURE 2

Results

Stable Isotopes and Cellulose Content in Living Conifers (1820–2010 A.D)

We first analyzed a subset of the samples of the project the Eastern Alpine Conifer Chronology for paleoclimatic studies (Nicolussi et al., 2009) represented by living trees grown in the last 190 years. They originate from eight larches and eight cembran pines, collected at three sites (VVR, FPCR and UZAR). The results, shown in Figure 2, are separated for species and for collection sites to illustrate and evaluate the effect of these two factors on isotope signatures.

δ²H Data

The δ²H data (Figure 2A) show a major difference between larch and cembran pine irrespective of the collection site. The mean larch value (4 trees at each site) at VRR site was −127.2 ± 6.6‰ and at FPCR site was −114.9 ± 6.3‰, while the mean cembran pine value at the UAZR site was −69.5 ± 11.0‰ and at VRR site −80.9 ± 11.8‰. The VRR cembran pine group shows a tail in the lower part, which is due to a single tree with a mean value of −96.8‰ that is more negative than the other three trees with mean values of −69.1, −73.7, −76.1‰. This made the distribution of the δ²H non-normal in the VRR cembran pine group, while it was normal in the other three groups. The data indicate that the δ²H values follow mainly species type.

δ¹⁸O Data

The δ¹⁸O data (Figure 2B) show that the density distribution and the means of the larch at the FPCR and at the VRR sites are very similar and are also similar to those of the cembran pine at the VRR site, but significantly different from those of the cembran pine at the UZAR site. A Wilcoxon test indicated a marginally significant difference between two larch groups (p = 0.025), contrarily, the difference between the two cembran pine groups was highly significant (p < 2 × 10^–16). The values of the cembran pine group from VRR were not significantly different from those of the larch from the same VRR site (p = 0.06) but were significantly different from the larch values at the FPCR site (p = 6.4 × 10^–5). The data indicate that the δ¹⁸O values follow mainly the site and not the tree species.

δ¹³C Data

The δ¹³C values of larch tended to be more negative than those of cembran pine (Figure 2C). The dispersion of the larch group values at VRR site was larger than that of larch group at FPCR site, but the means are similar with no significant differences. Also, the cembran pine values at the two regions were not significantly different. The only non-significant difference (Wilcoxon test: p = 0.455), although marginal, between the groups was observed for larch at VRR site and cembran pine at UZAR site. These values showed to be affected only to a minor extent by geography and species.

Cellulose Content

The cellulose content values of the four groups are shown in Figure 2D. The value distribution is more dispersed in the larch groups than in the cembran pine groups, independent from the site origin. A non-significant difference is found with the Wilcoxon test between the two larch groups (p = 0.133). The cellulose did not seem to confer major information and was not considered in the following part of the study.

Stable Isotopes of Subfossil Holocene Conifers

In order to verify whether the differences of living trees from Alpine sites shown above are not only due to the specific sites and restricted to the recent period, we extended the analysis to the complete wooden material of the Eastern Alpine Conifer Chronology (Nicolussi et al., 2009) collected at independent sites of the Alps, with the exclusion of the living trees of Figure 1. It consisted of >4000 wood samples from 120 trees (Table 1) and the sites of sample collection as listed in Figure 1.

δ²H Data

The δ²H values of the samples of the two species collected at the 12 sites are shown in Figure 3A as box-plots. It confirms an evident difference between the more negative values of the deciduous larch and the less negative values of the evergreen cembran pine samples, with no overlap for samples collected at the same site. The total mean for larch is −113.4 ± 9.7‰ and that for cembran pine −61.4 ± 11.4‰, values that are very close to those presented in Figure 2. An exception was the subfossil from spruce tree (Picea abies, PCAB) found at the HIB site the identity of which was initially unclear, since it is not easy to distinguish it from larch specimens by microscopic wood anatomy (Schoch et al., 2004). The δ²H values of these samples were in the range of those of the cembran pine (green in Figure 3A) and overlapped with those of the cembran pine that were collected at the same site, with no statistical difference (p = 0.74) between the two. They are trees of different species and different time periods, altogether the differences between larch and cembran pine are always in the range 40 to 50‰, while that at HIB between cembran pine and spruce was close to 0‰.

FIGURE 3

δ¹⁸O Data

The δ¹⁸O values from each site are shown in Figure 3B. The total mean for larch is 29.6 ± 1.1‰ and that for cembran pine 28.9 ± 0.5‰. A difference between the two species is not evident at sites VRR, HIB, MAZB, UA and ULFI. However, at the site ZER the two species showed a marginally significant difference (difference of the mean values of 0.6‰, p = 0.007), whereas at the sites MIS, AHST, MAZE, MAZF, MM the differences were highly significant (difference of the mean values ok −1 ot 0.6‰ p < 0.001). Boxplots of the 12 sites partially overlap, MIS and TSC are the only sites for which we found a difference of the mean values larger than 1‰. In MAZE, MAZAF, UA, ULFI, VRR sites the mean of larch was smaller than for cembran pine, for the remaining sites it was higher. The differences between the two species are in the range of 0.6 to −1‰ and do not show a common trend. The trees of some sites show values that are rather distant from the mean values and this is probably due to factors linked to the location (elevation, latitude, longitude, aspect) and others that influence the tree metabolism (e.g., soil composition), it can be noted that at sites where δ¹⁸O values are below total mean the corresponding δ²H are above or vice versa, suggesting that the two isotopes are subjected to distinct influences and that soil water is not the primary driver.

δ¹³C Data

The δ¹³C values of the trees are shown in Figure 3C. The total mean for larch is -22.6 ± 0.8‰ and that for cembran pine −22.2 ± 0.6‰. The difference of the values between the two species is significant at all sites with the exception of MAZB site (p = 0.05) and of UA (p = 0.001). The dotted line indicates the total mean value and most of the cembran pine values (blue boxplots) are above or close to the mean value, while the larch values (red boxplots) are below or close to it. In fact, for the sites where both species are present, the larch is always more depleted than the cembran pine. The difference plot shows that the only negative column is at the site HIB, spruce shows significantly less depleted values than cembran pine, while the oxygen values are similar. This might be due to some characteristics of the site, like the aspect, that are known to affect temperature-controlled changes in assimilation regulated by enzymatic isotopic fractionation processes (δ¹³C) (Esper et al., 2018).

Stable Isotopes in Contemporary Conifers

The above data showed that larch has δ²H values of cellulose that are much lower than those of cembran pine and spruce (Figures 2,3). In order to verify if this difference is attributed to the metabolism of the evergreen conifers versus that of deciduous larch, we looked for a larger number of conifer species. The Bern Botanic Garden provided us with samples of six different species of larch (deciduous), one pseudo-larch (deciduous), and 13 different species of evergreen conifers including pinus, taxus and cedrus, and the gymnosperm, broadleaf, deciduous Ginkgo Biloba. We processed the samples and analysed their isotope content as describe above. Figure 4 shows the results, in red are the deciduous conifers (larch and pseudo-larch), in green the evergreen conifers and in blue the values of Ginkgo Biloba, the horizontal lines are the mean value of each group and the dashed lines shows the interval of the standard deviation.

FIGURE 4

δ²H

The δ²H values of the deciduous larch group (mean −91.8 ± 7.6‰) is well separated from that of evergreen conifers (mean −51.1 ± 12.3‰) (Figure 4). This confirms that the deciduous conifer trees are strongly depleted in comparison with the evergreen ones with no overlap between the two groups and a highly significant difference (p = 1.92 e − 05). The values of Pseudolarix amabilis (tree #7 in Figure 4) are in full agreement with those of the other larches (#1–6). While the values of the Cedrus atlantica (#8), Cupressus nootkatensis (#10) and Taxus baccata (#19) are in agreement with those of the Pinus (#11–18), which are all evergreen conifers. Ginkgo biloba (#21) does not belong to the conifer order and is deciduous, but its values are close to those of the evergreen conifers.

δ¹⁸O

The difference between the δ¹⁸O value of the deciduous conifers (29.8 ± 0.8‰) and that of the evergreen conifers (29.1 ± 1.2‰) was not statistically significant (Figure 4). The distance between the two means was barely significant (p = 0.084). The absolute δ¹⁸O and δ²H values of samples from the botanic garden conifers are less depleted from the ones of the of the Alpine trees, this was attributed to higher depletion of the heavy isotopes of hydrogen and oxygen of source water at high elevations (Kern et al., 2014).

δ¹³C

Mean δ¹³C values show no statistically significant difference between the deciduous (−25.1 ± 1‰) and the evergreen groups (−25.4 ± 1.5‰) (Figure 4).

Discussion

A major finding of present work is the much higher depletion of deuterium in the cellulose of the deciduous larches compared to the evergreen conifer cembran pine growing at the same site. The δ²H values of larch that we found are in good agreement with the only previous measurement of the hydrogen isotope in larch cellulose (Hafner et al., 2011). The difference between δ²H of larch and cembran pine has been observed in eight trees covering the last 190 years, grown in Val Roseg valley, Grimselstausee and Val d’Herens valley (Figure 2) and the 120 trees covering the whole Holocene (Figure 3)The difference is large, with a δ²H mean of −127‰ for larch and −81‰ for the cembran pine and no overlapping values. This is a non-climatic signal and was not observed in the other two stable isotopes we analyzed, δ¹³C and δ¹⁸O, the values of which were similar in the two conifer species.

All but one of the trees of the Eastern Alpine Conifer Chronology were larch and cembran pine, and thus it was unclear whether the observed differences were specific to the two species or caused by different physiological processes of evergreen and deciduous conifers. The absence of difference between cembran pine and Picea abies in the site HIB is surprising as they are of different species and originate from different time periods, the cembran pine covered the period 5695–5516 B2K, some 1500 years before the spruce tree that covered the 4000–3766 B2K period. Thus, we measured conifer samples from the Bern Botanical Garden and the results showed that the cellulose of all the larch species including the Pseudo-larix amabilis had more depleted δ²H values than that of the evergreen conifers (Figure 4). The deciduous non-conifer Ginkgo biloba resulted to have δ²H value in the range of the evergreen ones. The range of δ²H values of the evergreen conifers (−30 to −70‰) is similar to the ones reported in deciduous broadleaf trees such as oak (−75 and −45‰) (Szczepanek et al., 2006; Etien et al., 2009), and in the conifers Pinus sylvestris (Gori et al., 2013) and spruce (−20 and −50‰) in three regions that cover a vast geographical area with different conditions (Loader et al., 2008). These values are quite distinct from the δ²H values we found in the larch family with a mean of −120‰ and range −150 to −86‰. Thus, all the deciduous conifers (larch and pseudo larches), present δ²H values that are very different from those of the evergreen conifers. This can be used for the identification of larch from spruce wood that have similar wood anatomic characteristics and are therefore very difficult to distinguish in subfossil specimens (Schoch et al., 2004). Moreover, our data show that for paleoclimatic use of isotopic dendrochronology the δ²H values of larch cellulose must not mixed with those of evergreen conifers, but can be used separately or after applying suitable corrections or scaling.

The major δ²H difference indicates that deciduous and evergreen conifer handle hydrogen in a different way, while the similar δ¹⁸O values indicate that they handle oxygen in the same way. The soil water is the source of the two elements and their isotopes, and its uptake by the roots does not cause fractionation of both isotope compositions (White et al., 1985; Ehleringer and Dawson 1992; Dawson and Ehleringer 1993; Lin and Sternberg 1993). Evapotranspiration in the leaves causes enrichment in similar extent of the two heavy isotopes (Flanagan et al., 1991; Yakir 1992; McCarroll and Loader 2004; Studer et al., 2015), thus we can conclude that this does not cause the difference of δ²H between the two species but leaving the δ¹⁸O in the two species unchanged (Figure 2). The difference will occur in some subsequent biological steps, likely at carboxylation or post-carboxylation levels that might differ in deciduous and evergreen conifers.

Deciduous larches are characterized by a short photosynthesis period that involves a higher photosynthesis capacity to compete with the surrounding evergreen trees (Tranquillini 1962; Tranquillini 1979). Consequently, they have faster rates of photosynthesis and net assimilation and a more efficient carboxylation system than evergreen conifers (Sweet and Wareing 1968; Larcher 1969; Fry and Phillips 1977; Gowin et al., 1980). The biochemical steps that amplify δ²H depletion when photosynthetic activity is high have been identified (Cormier et al., 2018): when the supply of photosynthetic carbohydrates is low, the produced organic compounds are relatively ²H-enriched (Cormier et al., 2018), when the photosynthetic activity is high the produced organic compounds are relatively ²H-depleted (Yakir 1992). Moreover, larch is known to have a lower heterotrophic metabolism than the surrounding evergreen conifers (Vaganov et al., 2009; Rinne et al., 2015).

In cembran pine heterotrophism increase with the altitude (Wieser et al., 2019). The finding that the δ²H difference between larches and evergreen conifers is lower in Bern (∼500 m elevation) than at >2000 m of tree-line (40‰ in Figure 4 vs 52‰ in Figure 3) support the hypothesis that heterotrophism affect δ²H but the difference between the two metabolisms is maintained.

Regarding the δ¹³C values they were all within the range of −23‰ to −28‰, as previously described in C3 plants (Sternberg and DeNiro 1983; Leuenberger 1998). This value in leaf sugar is controlled mainly by the ratio c_i/c_a (c_i being the CO₂ concentration in the intercellular space of leaves and c_a the CO₂ concentration of the atmosphere), which is related to photosynthetic water-use efficiency, WUE, (also called transpiration efficiency) (Farquhar and Sharkey 1982). In our results δ¹³C values of larch were consistently slightly more negative than those of cembran pine (Figure 3), in agreement with previous data showing that larch species have more negative carbon values than the evergreen conifers because of their lower water-use efficiency (WUE) (Gower and Richards, 1990; Kloeppel et al., 1998). Altogether the results indicate that neither the species nor geography have major effect on δ¹³C signature.

Regarding the δ¹⁸O values we found no evident difference between the two species, in accordance with the interpretation that these values are more sensitive to the source water signal than to other external variables, such as relative humidity and vapor pressure deficit. We found a δ¹⁸O difference between the sites in the Swiss and the Tyrolian Alps, likely due to the source water, the δ¹⁸O of which is known to be linked to the altitude and to a continental effect (Rozanski et al., 1993). This indicates that δ¹⁸O values of different conifers species can be used together for paleoclimatic studies, after taking into account the geographical effects.

Conclusion

We show that the deciduous larch family of trees presents unique highly negative δ²H values compared to evergreen conifers. We attributed the difference to a higher photosynthetic and metabolic activity of larch that drive a reduced exchange of δ²H sucrose with the surrounding enriched leaf water leading to δ²H depletion of the photosynthetic products, in comparison with cembran pine and more in general with the evergreen conifers. No significant differences of δ¹⁸O values were found between trees grown at the same site and under similar conditions, but differences were found among trees grown at different sites. Minor more negative δ¹³C values of larch were consistently found than those of the cembran pine, but no evident driver of the difference could be identified. We also conclude that for paleoclimatic studies, for δ²H values the evergreen and deciduous conifers should be keep separately or normalized for the tree species, while δ¹⁸O does not seem to be affected by tree species, but by the source water of each site. Regarding δ¹³C we do not have yet a clear driver of the difference in carbon isotope. We suggest the determination of the δ²H values as tool for distinguishing larch samples from those of other conifers, such as spruce, which are difficult to identify in subfossil wood.

Funding

The project is funded by the Swiss National Science Foundation (SNSF, 2000212_144255, 200020_172550) as well as by the Austrian Science Fund (FWF, grant I-1183-N19) and is supported by the Oeschger Center for Climate Change Research, University of Bern, Bern, Switzerland (OCCR).

References

• 1

  AugustiA.BetsonT. R.SchleucherJ. (2006). Hydrogen exchange during cellulose synthesis distinguishes climatic and biochemical isotope fractionations in tree rings. New Phytol.172, 490–499. 10.1111/j.1469-8137.2006.01843.x

  • CrossRef
  • Google Scholar

• 2

  AugustiA. (2007). Monitoring climate and plant physiology using deuterium isotopomers of carbohydrates Fysiologisk botanik. PhD thesis. Umea (Sweden): Umea University.

  • Google Scholar

• 3

  AugustiA.BetsonT. R.SchleucherJ. (2008). Deriving correlated climate and physiological signals from deuterium isotopomers in tree rings. Chemical Geology252 (1–2), 1–8. 10.1016/j.chemgeo.2008.01.011

  • CrossRef
  • Google Scholar

• 4

  BorellaS.LeuenbergerM.SaurerM. (1999). Analysis of δ¹⁸O in tree rings: wood-cellulose comparison and method dependent sensitivity. J. Geophys. Res.104, 19267–19273. 10.1029/1999jd900298

  • CrossRef
  • Google Scholar

• 5

  BorellaS.LeuenbergerM.SaurerM.SiegwolfR. (1998). Reducing uncertainties in δ¹³C analysis of tree rings: pooling, milling, and cellulose extraction. J. Geophys. Res.103, 19519–19526. 10.1029/98jd01169

  • CrossRef
  • Google Scholar

• 6

  BoysenB. M.EvansM. N.BakerP. J. (2014). δ¹⁸O in the tropical conifer Agathis robusta records ENSO-related precipitation variations. PLoS One9, e102336. 10.1371/journal.pone.0102336

  • CrossRef
  • Google Scholar

• 7

  BrienenR. J. W.HelleG.PonsT. L.GuyotJ.-L.GloorM. (2012). Oxygen isotopes in tree rings are a good proxy for amazon precipitation and El Nino-Southern oscillation variability. Proc. Natl. Acad. Sci.109, 16957–16962. 10.1073/pnas.1205977109

  • CrossRef
  • Google Scholar

• 8

  BurkR. L.StuiverM. (1981). Oxygen isotope ratios in trees reflect mean annual temperature and humidity. Science211, 1417–1419. 10.1126/science.211.4489.1417

  • CrossRef
  • Google Scholar

• 9

  CoplenT. B. (1994). Reporting of stable hydrogen, carbon, and oxygen isotopic abundances (technical report). Pure Appl. Chem.66, 273–276. 10.1351/pac199466020273

  • CrossRef
  • Google Scholar

• 10

  CormierM.-A.WernerR. A.LeuenbergerM. C.KahmenA. (2019). 2H-enrichment of cellulose and n-alkanes in heterotrophic plants. Oecologia189, 365–373. 10.1007/s00442-019-04338-8

  • CrossRef
  • Google Scholar

• 11

  CormierM.-A.WernerR. A.SauerP. E.GröckeD. R.LeuenbergerM. C.WielochT.et al (2018). 2H-fractionations during the biosynthesis of carbohydrates and lipids imprint a metabolic signal on the δ₂H values of plant organic compounds. New Phytol.218, 479–491. 10.1111/nph.15016

  • CrossRef
  • Google Scholar

• 12

  DawsonT. E.EhleringerJ. R. (1993). Isotopic enrichment of water in the “woody” tissues of plants: implications for plant water source, water uptake, and other studies which use the stable isotopic composition of cellulose. Geochimica et Cosmochimica Acta57, 3487–3492. 10.1016/0016-7037(93)90554-a

  • CrossRef
  • Google Scholar

• 13

  EhleringerJ. R.DawsonT. E. (1992). Water uptake by plants: perspectives from stable isotope composition. Plant Cell Environ.15, 1073–1082. 10.1111/j.1365-3040.1992.tb01657.x

  • CrossRef
  • Google Scholar

• 14

  EsperJ.HolzkämperS.BüntgenU.SchöneB.KepplerF.HartlC.et al (2018). Site-specific climatic signals in stable isotope records from Swedish pine forests. Trees32, 855–869. 10.1007/s00468-018-1678-z

  • CrossRef
  • Google Scholar

• 15

  EtienN.DauxV.Masson-DelmotteV.MestreO.StievenardM.GuilleminM. T.et al (2009). Summer maximum temperature in northern France over the past century: instrumental data versus multiple proxies (tree-ring isotopes, grape harvest dates and forest fires). Climatic Change94, 429–456. 10.1007/s10584-008-9516-8

  • CrossRef
  • Google Scholar

• 16

  FarquharG. D.SharkeyT. D. (1982). Stomatal conductance and photosynthesis. Annu. Rev. Plant. Physiol.33 (1), 317–345. 10.1146/annurev.pp.33.060182.001533

  • CrossRef
  • Google Scholar

• 17

  FilotM. S.LeuenbergerM.PazdurA.BoettgerT. (2006). Rapid online equilibration method to determine the D/H ratios of non-exchangeable hydrogen in cellulose. Rapid Commun. Mass Spectrom.20, 3337–3344. 10.1002/rcm.2743

  • CrossRef
  • Google Scholar

• 18

  FlanaganL. B.ComstockJ. P.EhleringerJ. R. (1991). Comparison of modeled and observed environmental influences on the stable oxygen and hydrogen isotope composition of leaf water in Phaseolus vulgaris L. Plant Physiol.96, 588–596. 10.1104/pp.96.2.588

  • CrossRef
  • Google Scholar

• 19

  FrankD.PoulterB.SaurerM.EsperJ.HuntingfordC.HelleG.et al (2015). Water-use efficiency and transpiration across European forests during the Anthropocene. Nat. Clim. Change5, 579–583. 10.1038/nclimate2614

  • CrossRef
  • Google Scholar

• 20

  FryD.PhillipsI. (1977). Photosynthesis of conifers in relation to annual growth cycles and dry matter production: II. Seasonal photosynthetic capacity and mesophyll ultrastructure in Abies grandis, Picea sitchensis, Tsuga heterophylla and Larix leptolepis growing in SW England. Physiol. Plant.40, 300–306. 10.1111/j.1399-3054.1977.tb04077.x

  • CrossRef
  • Google Scholar

• 21

  GagenM.McCarrollD.EdouardJ.-L. (2004). Latewood width, maximum density, and stable carbon isotope ratios of pine as climate indicators in a dry subalpine environment, French Alps. Alpine Research Arctic, Antarctic36 (2), 166–171. 10.1657/1523-0430(2004)036[0166:lwmdas]2.0.co;2

  • CrossRef
  • Google Scholar

• 22

  GoriY.WehrensR.GreuleM.KepplerF.ZillerL.La PortaN.CaminF. (2013). Carbon, hydrogen and oxygen stable isotope ratios of whole wood, cellulose and lignin methoxyl groups of Picea abies as climate proxies. Rapid Commun. Mass Spectro.27 (1), 265–275. 10.1002/rcm.6446

  • CrossRef
  • Google Scholar

• 23

  GowerS. T.RichardsJ. H. (1990). Larches: deciduous conifers in an evergreen world. Bioscience40, 818–826. 10.2307/1311484

  • CrossRef
  • Google Scholar

• 24

  GowinT.LourtiouxA.MousseauM. (1980). Influence of constant growth temperature upon the productivity and gas exchange of seedlings of Scots pine and European larch. Forest Science26, 301–309. 10.1093/forestscience/26.2.301

  • CrossRef
  • Google Scholar

• 25

  HafnerP.RobertsonI.McCarrollD.LoaderN. J.GagenM.BaleR. J.et al (2011). Climate signals in the ring widths and stable carbon, hydrogen and oxygen isotopic composition of Larix decidua growing at the forest limit in the southeastern European Alps. Trees25 (6), 1141–1154. 10.1007/s00468-011-0589-z

  • CrossRef
  • Google Scholar

• 26

  HilasvuoriC.ZangC.BüntgenU.EsperJ.RotheA.GöttleinA.et al (2014). Uniform climate sensitivity in tree-ring stable isotopes across species and sites in a mid-latitude temperate forest. Tree Physiol.35, 4–15. 10.1093/treephys/tpu096

  • CrossRef
  • Google Scholar

• 27

  HintzeJ. L.NelsonR. D. (1998). Violin plots: a box plot-density trace synergism. Am. Stat.52, 181–184. 10.1080/00031305.1998.10480559

  • CrossRef
  • Google Scholar

• 28

  KahmenA.SachseD.ArndtS. K.TuK. P.FarringtonH.VitousekP. M.et al (2011). Cellulose δ¹⁸O is an index of leaf-to-air vapor pressure difference (VPD) in tropical plants. Proc. Natl. Acad. Sci. USA108, 1981–1986. 10.1073/pnas.1018906108

  • CrossRef
  • Google Scholar

• 29

  KernZ.KohánB.LeuenbergerM. (2014). Precipitation isoscape of high reliefs: interpolation scheme designed and tested for monthly resolved precipitation oxygen isotope records of an Alpine domain. Atmos. Chem. Phys.14, 1897–1907. 10.5194/acp-14-1897-2014

  • CrossRef
  • Google Scholar

• 30

  KhezamiL.ChetouaniA.TaoukB.CapartR. (2005). Production and characterisation of activated carbon from wood components in powder: cellulose, lignin, xylan. Powder Technology157, 48–56. 10.1016/j.powtec.2005.05.009

  • CrossRef
  • Google Scholar

• 31

  KimakA.KernZ.LeuenbergerM. (2015). Qualitative distinction of autotrophic and heterotrophic processes at the leaf level by means of triple stable isotope (C–O–H) patterns. Front. Plant Sci.6, 1008. 10.3389/fpls.2015.01008

  • CrossRef
  • Google Scholar

• 32

  KimakA.LeuenbergerM. (2015). Are carbohydrate storage strategies of trees traceable by early-latewood carbon isotope differences?. Trees29 (3), 859–870. 10.1007/s00468-015-1167-6

  • CrossRef
  • Google Scholar

• 33

  KloeppelB. D.TreichelI. W.KharukS.GowerS. T. (1998). Foliar carbon isotope discrimination in Larix species and sympatric evergreen conifers: a global comparison, Oecologia. 114, 153–159. 10.1007/s004420050431

  • CrossRef
  • Google Scholar

• 34

  LabuhnI.DauxV.PierreM.StievenardM.GirardclosO.FéronA.et al (2014). Tree age, site and climate controls on tree ring cellulose δ¹⁸O: a case study on oak trees from south-western France. Dendrochronologia32, 78–89. 10.1016/j.dendro.2013.11.001

  • CrossRef
  • Google Scholar

• 35

  LarcherW. (1969). Effect of environmental and physiological variables on the carbon dioxide gas exchange of trees. Photosynthetica3, 167–198.

  • Google Scholar

• 36

  LeuenbergerM. (1998). Stable isotopes in tree rings as climate and stress indicators. Zürich, Switzerland: Vdf Hochschulverlag AG.

  • Google Scholar

• 37

  LeuenbergerM. (2007). To what extent can ice core data contribute to the understanding of plant ecological developments of the past?. Terres. Eco.1, 211–233. 10.1016/s1936-7961(07)01014-7

  • CrossRef
  • Google Scholar

• 38

  LibbyL. M.PandolfiL. J.PaytonP. H.MarshallJ.IIIBeckerB.Giertz-SienbenlistV. (1976). Isotopic tree thermometers. Nature261, 284. 10.1038/261284a0

  • CrossRef
  • Google Scholar

• 39

  LinG.SternbergL. S. L. (1993). “Hydrogen isotopic fractionation by plant roots during water uptake in coastal wetland plants,” in Stable isotopes and plant carbon-water relations (Cambridge, MA: Academic Press), 497–510. 10.1016/b978-0-08-091801-3.50041-6

  • CrossRef
  • Google Scholar

• 40

  LoaderN. J.RobertsonI.McCarrollD. (2003). Comparison of stable carbon isotope ratios in the whole wood, cellulose and lignin of oak tree-rings. Palaeogeogr. Palaeocl.196, 395–407. 10.1016/s0031-0182(03)00466-8

  • CrossRef
  • Google Scholar

• 41

  LoaderN. J.SantilloP. M.Woodman-RalphJ. P.RolfeJ. E.HallM. A.GagenM.et al (2008). Multiple stable isotopes from oak trees in southwestern Scotland and the potential for stable isotope dendroclimatology in maritime climatic regions. Chem. Geol.252, 62–71. 10.1016/j.chemgeo.2008.01.006

  • CrossRef
  • Google Scholar

• 42

  LoaderN. J.Street-PerrottF. A.DaleyT. J.HughesP. D. M.KimakA.LevaničT.et al (2014). Simultaneous determination of stable carbon, oxygen, and hydrogen isotopes in cellulose. Anal. Chem.87, 376–380. 10.1021/ac502557x

  • CrossRef
  • Google Scholar

• 43

  LoaderN. J.Street-PerrottF. A.DaleyT. J.HughesP. D. M.KimakA.LevaničT.et al (2015). Simultaneous determination of stable carbon, oxygen, and hydrogen isotopes in cellulose. Anal. Chem.87, 376–380. 10.1021/ac502557x

  • CrossRef
  • Google Scholar

• 44

  McCarrollD.LoaderN. J. (2004). Stable isotopes in tree rings. Quaternary Science Reviews23, 771–801. 10.1016/j.quascirev.2003.06.017

  • CrossRef
  • Google Scholar

• 45

  NabeshimaE.NakatsukaT.KagawaA.HiuraT.FunadaR. (2018). Seasonal changes of δD and δ¹⁸O in tree-ring cellulose of Quercus crispulasuggest a change in post-photosynthetic processes during earlywood growth. Tree Physiol.38, 1829–1840. 10.1093/treephys/tpy068

  • CrossRef
  • Google Scholar

• 46

  NicolussiK.KaufmannM.MelvinT. M.Van Der PlichtJ.SchießlingP.ThurnerA. (2009). A 9111 year long conifer tree-ring chronology for the European Alps: a base for environmental and climatic investigations. The Holocene19, 909–920. 10.1177/0959683609336565

  • CrossRef
  • Google Scholar

• 47

  RinneK. T.SaurerM.KirdyanovA. V.LoaderN. J.BryukhanovaM. V.WernerR. A.et al (2015). The relationship between needle sugar carbon isotope ratios and tree rings of larch in Siberia. Tree Physiol.35 (11), 1192–205. 10.1093/treephys/tpv096

  • CrossRef
  • Google Scholar

• 48

  RodenJ. S.EhleringerJ. R. (2007). Summer precipitation influences the stable oxygen and carbon isotopic composition of tree-ring cellulose in Pinus ponderosa. Tree Physiol.27, 491–501. 10.1093/treephys/27.4.491

  • CrossRef
  • Google Scholar

• 49

  RozanskiK.Araguás‐AraguásL.GonfiantiniR. (1993). Isotopic patterns in modern global precipitation. Clim. Change Contin. Isoto. Rec.78, 1–36.

  • Google Scholar

• 50

  ScheideggerY.SaurerM.BahnM.SiegwolfR. (2000). Linking stable oxygen and carbon isotopes with stomatal conductance and photosynthetic capacity: a conceptual model. Oecologia125, 350–357. 10.1007/s004420000466

  • CrossRef
  • Google Scholar

• 51

  SchieglW. E. (1974). Climatic significance of deuterium abundance in growth rings of Picea. Nature251, 582. 10.1038/251582a0

  • CrossRef
  • Google Scholar

• 52

  SchochW.HellerI.SchweingruberF. H.KienastF. (2004). Wood anatomy of central European species. Birmensdorf, Switzerland: Swiss Federal Institute for Forest.

  • Google Scholar

• 53

  ShapiroS. S.FranciaR. S. (1972). An approximate analysis of variance test for normality. J. Am. Stat. Assoc.67, 215–216. 10.1080/01621459.1972.10481232

  • CrossRef
  • Google Scholar

• 54

  SternbergL.DeNiroM. J. (1983). Isotopic composition of cellulose from C3, C4, and CAM plants growing near one another. Science220, 947–949. 10.1126/science.220.4600.947

  • CrossRef
  • Google Scholar

• 55

  StuderM. S.SiegwolfR. T. W.LeuenbergerM.AbivenS. (2015). Multi-isotope labelling of organic matter by diffusion of. Biogeosciences12, 1865–1879. 10.5194/bg-12-1865-2015

  • CrossRef
  • Google Scholar

• 56

  SweetG. B.WareingP. (1968). A comparison of the seasonal rates of dry matter production of three coniferous species with contrasting patterns of growth. Ann. Bot.32, 721–734. 10.1093/oxfordjournals.aob.a084245

  • CrossRef
  • Google Scholar

• 57

  SzczepanekM.PazdurA.PawełCzykS.BöttgerT.HauptM.HałasS.et al (2006). Hydrogen, carbon and oxygen isotopes in pine and oak tree rings from southern Poland as climatic indicators in years 1900–2003. Geochronometria: Journal on Methods Applications of Absolute Chronology 25, 67–73.

  • Google Scholar

• 58

  TangK.FengX.EttlG. J. (2000). The variations in δD of tree rings and the implications for climatic reconstruction. Geochem. Cosmochim. Acta64, 1663–1673. 10.1016/s0016-7037(00)00348-3

  • CrossRef
  • Google Scholar

• 59

  TranquilliniW. (1962). Beitrag zur Kausalanalyse des Wettbewerbs ökologisch verschiedener Holzarten. Ber. Deut. Bot. Ges.75, 356–364.

  • Google Scholar

• 60

  TranquilliniW.FrankD.EsperJ.AndreuL.BednarzZ.BerningerF.et al (1979). Physiological ecology of the alpine timberline. Berlin, Germany: Springer-Verlag.

  • Google Scholar

• 61

  TreydteK. (2007). Signal strength and climate calibration of a European tree‐ring isotope network. Geophys. Res. Lett.34 (24). 10.1029/2007GL031106

  • CrossRef
  • Google Scholar

• 62

  VaganovE. A.SchulzeE. D.SkomarkovaM. V.KnohlA.BrandW. A.RoscherC. (2009). Intra-annual variability of anatomical structure and δ 13 C values within tree rings of spruce and pine in alpine, temperate and boreal Europe. Oecologia161 (4), 729–745. 10.1007/s00442-009-1421-y

  • CrossRef
  • Google Scholar

• 63

  WhiteJ. W.CookE. R.LawrenceJ. R. (1985). The DH ratios of sap in trees: implications for water sources and tree ring DH ratios. Geochem. Cosmochim. Acta49, 237–246. 10.1016/0016-7037(85)90207-8

  • CrossRef
  • Google Scholar

• 64

  WieserG.OberhuberW.GruberA. (2019). Effects of climate change at treeline: Lessons from space-for-time studies, manipulative experiments, and long-term observational records in the Central Austrian Alps. Forests10 (6), 508. 10.3390/f10060508

  • CrossRef
  • Google Scholar

• 65

  WilcoxonF.KattiS.WilcoxR. A. (1970). “Critical values and probability levels for the Wilcoxon rank sum test and the Wilcoxon signed rank test,” in Selected tables in mathematical statistics. (Providence, RI, American Mathematical Soc.), 1, 171–259.

  • Google Scholar

• 66

  XuC.SanoM.NakatsukaT. (2011). Tree ring cellulose δ¹⁸O of Fokienia hodginsii in northern Laos: a promising proxy to reconstruct ENSO?. J. Geophys. Res.: Atmosphere116. 10.1029/2011jd016694

  • CrossRef
  • Google Scholar

• 67

  YakirD.DeNiroM.EphrathJ. (1990). Effects of water stress on oxygen, hydrogen and carbon isotope ratios in two species of cotton plants. Plant Cell Environ.13, 949–955. 10.1111/j.1365-3040.1990.tb01985.x

  • CrossRef
  • Google Scholar

• 68

  YakirD. (1992). Variations in the natural abundance of oxygen‐18 and deuterium in plant carbohydrates. Plant Cell Environ.15, 1005–1020. 10.1111/j.1365-3040.1992.tb01652.x

  • CrossRef
  • Google Scholar

• 69

  ZhouY.GriceK.Stuart‐WilliamsH.HocartC. H.GesslerA.FarquharG. D. (2016). Hydrogen isotopic differences between C3 and C4 land plant lipids: consequences of compartmentation in C4 photosynthetic chemistry and C3 photorespiration. Plant Cell Environ.39, 2676–2690. 10.1111/pce.12821

  • CrossRef
  • Google Scholar

• 70

  ZieglerH. (1995). “Deuterium content in organic material of hosts and their parasites.” in Ecophysiology of photosynthesis. (Berlin, Germany: Springer), 393–408.

  • Google Scholar

• 71

  ZieglerH.OsmondC.StichlerW.TrimbornP. (1976). Hydrogen isotope discrimination in higher plants: correlations with photosynthetic pathway and environment. Planta128, 85–92. 10.1007/bf00397183

  • CrossRef
  • Google Scholar

• 72

  ZiehmerM. M.NicolussiK.SchlüchterC.LeuenbergerM. (2018). Preliminary evaluation of the potential of tree-ring cellulose content as a novel supplementary proxy in dendroclimatology. Biogeosciences15, 1047–1064. 10.5194/bg-15-1047-2018

  • CrossRef
  • Google Scholar