Elsevier

Dendrochronologia

Volume 65, February 2021, 125804
Dendrochronologia

ORIGINAL ARTICLE
The dendroclimatic value of oak stable isotopes

https://doi.org/10.1016/j.dendro.2020.125804Get rights and content

Highlights

  • Oak δ18O and δ13C is not influenced by polyvinyl acetate emulsion.

  • Oak δ18O and δ13C from co-occurring species correlate significantly.

  • Oak δ18O and δ13C from different sites and tree ages correlate significantly.

  • Oak δ18O and δ13C should be annually-resolved for climate reconstruction.

Abstract

Tree-ring stable oxygen and carbon isotope ratios (δ18O and δ13C) are an important archive for climate reconstructions. However, it remains unclear whether the polyvinyl acetate emulsion, often used for the preservation and fixation of wood samples, influences δ18O and δ13C signals. Further uncertainties are associated with the possible effects of geographical origin and cambial age of historical samples. Here, we present annually-resolved and absolutely-dated δ18O and δ13C measurements of 21 living oaks (Quercus robur and Q. petraea) from the Czech Republic. We find that the δ18O and δ13C signals in the extracted alpha-cellulose are not affected by polyvinyl acetate treatment. Covering the entire 20th century and reaching until 2018 CE, our dataset reveals spatial and temporal coherency within and between the individual δ18O and δ13C chronologies of different oak species, sample locations, and tree ages. Highly significant (p < 0.01) Pearson’s correlation coefficients of the site-specific δ13C and δ18O chronologies range from 0.48–0.77 and 0.36–0.56, respectively. The isotopic inter-series correlations of Q. robur and Q. petraea from the same site are 0.75 and 0.43 for the mean δ13C and δ18O values, respectively. Significant (p < 0.01) correlations of 0.49 and 0.84 are found for δ13C and δ18O, respectively, when all measurements from all sampling locations and tree ages are included. Our study shows that non-pooled oak δ18O and δ13C measurements from both species, different locations, and diverse tree ages can be combined into robust isotopic chronologies for climate reconstructions.

Introduction

Changes in tree-ring stable oxygen (δ18O; 18O/16O) and carbon (δ13C; 13C/12C) isotopes are modulated mainly by plant physiological processes associated with photosynthetic CO2 uptake (Farquhar et al., 1989; Gessler et al., 2014), as well as by the transpiration of water (Barbour and Farquhar, 2004). These factors are often influenced by environmental perturbations and/or changes in weather and climate. It has been shown that stable isotopes in tree rings reflect vegetation responses to changes in temperature, relative humidity, precipitation, and/or light (Roden et al., 2000; Leavitt, 2010). However, changes in the isotopic composition of atmospheric CO2 and source water must also be taken into account. Species-specific sensitivity to abiotic factors, as well as the overall constitution of the tree species including the root system, may lead to differences in the δ18O and δ13C ratios (Ponton et al., 2001; Song et al., 2014; Barbeta and Peñuelas, 2017).

Valuable climate information can be extracted from δ18O and δ13C isotopic ratios in tree rings (McCarroll and Loader, 2004; De Boer et al., 2019). Tree-ring stable isotopes can improve high-resolution reconstructions of hydroclimatic variability as compared to more traditional tree-ring parameters including tree-ring width and latewood density (Büntgen et al., 2020; Shestakova and Martínez-Sancho, 2021). Tree-ring stable isotopes have been used to reconstruct various hydroclimatic parameters, such as drought and sunlight duration in the European Alps (Kress et al., 2010); relative humidity, cloud cover, and drought in the eastern Carpathians (Nagavciuc et al., 2019); the Palmer drought severity index (PDSI) in China (Liu et al., 2019); the standardized precipitation evaporation index (SPEI) in southern Central Siberia (Shestakova et al., 2017); and precipitation amounts in Pakistan (Treydte et al., 2009). While there is evidence that δ18O can reflect changes in relative humidity and precipitation in semi-arid regions (Szejner et al., 2018), the links between δ18O and δ13C and climate signals are often stronger under wetter rather than drier conditions (Esper et al., 2018). In humid environments, variations in δ13C reflect the activity of Rubisco, the primary carbon-fixing enzyme in C3 plants (Hafner et al., 2014). In dry environments, δ13C is mainly controlled by stomatal conductance and decreases with increasing drought severity or vapor pressure deficit (McCaroll and Loader et al., 2004; Young et al., 2015). Fractionation of oxygen isotopes prior to and during wood formation is more complex than carbon fractionation (Gessler et al., 2014). Values of δ18O received increasing attention for long-term climate reconstruction in the last decade (Nakatsuka et al., 2020). In addition to providing reliable climate archives, δ18O chronologies have been shown to improve cross-dating compared to conventional tree-ring widths (Loader et al., 2020).

A multi-millennial standard oak ring-width chronology of ∼4000 core and disc samples from living trees, historical timbers, archaeological remains, and subfossil wood has been developed for the Czech Republic (Dobrovolný et al., 2015; Prokop et al., 2017; Dobrovolný et al., 2018). However, the suitability of stable isotopes from such a heterogeneous compilation has to be demonstrated before it can be used in paleoclimatic studies. Putative uncertainties related to species identification, geographical origin, cambial age, and the presence of polyvinyl may possibly limit the exploitation of all relict samples for climate reconstructions. Here, we present annually-resolved and absolutely-dated δ18O and δ13C ratios in 21 living oaks from different locations along elevational, longitudinal and ecological gradients in the Czech Republic. We use this dataset to elucidate the natural variability of oak stable isotopes. Moreover, we test whether polyvinyl acetate, an emulsion often used for the preservation of relict samples, has any influence on isotopic signals. More specifically, we aim to address the following issues: (1) The δ18O and δ13C ratios in the extracted alpha cellulose should be tested to ensure they are not influenced by the polyvinyl acetate emulsion. (2) Since the dataset includes two widespread oak species—English oak (Quercus robur L.) and sessile oak (Quercus petraea [Matt.] Liebl.)—which are anatomically indistinguishable from one another (Schoch et al., 2004), the effect of any species differences on δ18O and δ13C ratios must be investigated. Previous studies have revealed contradictory results (Ponton et al., 2001; Rybníček et al., 2016). (3) Uncertainties regarding the extent to which the geographical origin of the wood influences δ18O and δ13C signals must be clarified (Roden et al., 2005; Shestakova et al., 2017; Zalloni et al., 2019). (4) Age-related changes in physiological activity and tree-ring morphology, including tree‐ring width and maximum latewood density (Genet et al., 2010), might have an effect on tree-ring stable isotopes (Etien et al., 2008; Labuhn et al., 2013). If this is the case, it makes using historical wood samples for isotopic-based climate reconstruction particularly problematic.

Section snippets

Tree-ring network

Our sampling sites represent the natural, spatial and altitudinal distribution of oak (Q. robur and Q. petraea) in the Czech Republic (Fig. 1). Moreover, young (<50-year-old) and old (>100-year-old) trees were sampled to investigate the impact of tree age on stable isotopes. In total, 21 trees from seven localities in western Bohemia (sites #1, 2, and 6) and eastern Moravia (sites #3–5, and 7) were sampled (Table 1). This network includes trees from low (170–250 m a.s.l.; sites #3, and 5–7) and

Site effect

The TRW, δ18O, and δ13C time series (1920–2018 CE) of our sampled trees (Table 1) are shown in Fig. 3. Both the uncorrected δ13C time series and the time series corrected for the increase in atmospheric CO2 due to the burning of fossil fuels are shown (Figs. 3C, D). Hereafter, all presented δ13C time series refer to atmospherically corrected data.

Inter-site correlations of δ13C series, expressed as a Pearson’s correlation coefficient (r), range between 0.48 and 0.77 (p < 0.01) at most

Conclusions

Our annually-resolved and absolutely-dated δ18O and δ13C measurements of 21 living oak trees from the Czech Republic exhibit significant coherency across species, locations, and tree ages. Moreover, the applied extraction method removes any biases resulting from the polyvinyl acetate emulsions that are commonly used to preserve wood samples. The pooling of samples is often motivated by lower costs and shorter analysis times (Leavitt, 2010); however, non-pooled data enable the important

Declaration of Competing Interest

The authors report no declarations of interest.

Acknowledgements

This research was supported by the Czech Science Foundation (grant 17-22102S). The work by O.U. and A.A. on the development of the methods was supported by SustES – Adaptation strategies for sustainable ecosystem services and food security under adverse environmental conditions (CZ.02.1.01/0.0/0.0/16_019/0000797). The authors would also like to thank Inna Roshka and Josef Čáslavský for their work in the laboratory, and two anonymous referees who kindly comment on previous versions of this

References (54)

  • T.A. Shestakova et al.

    Stories hidden in tree rings: a review on the application of stable carbon isotopes to dendrosciences

    Dendrochronologia

    (2021)
  • K.S. Treydte et al.

    Impact of climate and CO2 on a millennium-long tree-ring carbon isotope record

    Geochim. Cosmochim. Acta

    (2009)
  • R. Alley

    The Two-Mile Time Machine

    (2000)
  • A. Barbeta et al.

    Increasing carbon discrimination rates and depth of water uptake favor the growth of Mediterranean evergreen trees in the ecotone with temperate deciduous forests

    Glob. Change Biol.

    (2017)
  • M.M. Barbour et al.

    Do pathways of water movement and leaf anatomical dimensions allow development of gradients in H218O between veins and the sites of evaporation within leaves?

    Plant Cell Environ.

    (2004)
  • T. Boettger et al.

    Wood cellulose preparation methods and mass spectrometric analyses of δ13C, δ18O, and nonexchangeable δ2H values in cellulose, sugar, and starch: an interlaboratory comparison

    Anal. Chem.

    (2007)
  • U. Büntgen et al.

    No age trends in oak stable isotopes

    Paleoceanogr. Palaeoclimatol.

    (2020)
  • H.J. de Boer et al.

    Tree-ring isotopes suggest atmospheric drying limits temperature–growth responses of treeline bristlecone pine

    Tree Physiol.

    (2019)
  • P. Dobrovolný et al.

    A tree-ring perspective on temporal changes in the frequency and intensity of hydroclimatic extremes in the territory of the Czech Republic since 761 AD

    Clim. Past

    (2015)
  • P. Dobrovolný et al.

    Recent growth coherence in long-term oak (Quercus spp.) ring width chronologies in the Czech Republic

    Clim. Res.

    (2016)
  • P. Dobrovolný et al.

    May–July precipitation reconstruction from oak tree-rings for Bohemia (Czech Republic) since AD 1040

    Int. J. Climatol.

    (2018)
  • J. Dombrosky

    A∼ 1000-year 13C Suess correction model for the study of past ecosystems

    Holocene

    (2020)
  • J. Esper et al.

    Long-term summer temperature variations in the Pyrenees from detrended stable carbon isotopes

    Geochronometria

    (2015)
  • J. Esper et al.

    Site-specific climatic signals in stable isotope records from Swedish pine forests

    Trees

    (2018)
  • N. Etien et al.

    A bi-proxy reconstruction of Fontainebleau (France) growing season temperature from A.D. 1596 to 2000

    Clim. Past

    (2008)
  • G.D. Farquhar et al.

    Carbon isotope discrimination and photosynthesis

    Annu. Rev. Plant Biol.

    (1989)
  • H. Genet et al.

    Age-related variation in carbon allocation at tree and stand scales in beech (Fagus sylvatica L.) and sessile oak (Quercus petraea (Matt.) Liebl.) using a chronosequence approach

    Tree Physiol.

    (2010)
  • Cited by (0)

    View full text