Next Article in Journal
Distribution Patterns and Environmental Determinants of Invasive Alien Plants on Subtropical Islands (Fujian, China)
Next Article in Special Issue
Research Progress on Ecological Carrying Capacity and Ecological Security, and Its Inspiration on the Forest Ecosystem in the Karst Desertification Control
Previous Article in Journal
Effects and Modification Mechanisms of Different Plasma Treatments on the Surface Wettability of Different Woods
Previous Article in Special Issue
Soil Carbon and Nitrogen Stocks and Their Influencing Factors in Different-Aged Stands of Sand-Fixing Caragana korshinskii in the Mu Us Desert of Northwest China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tree-Ring δ13C and Intrinsic Water-Use Efficiency Reveal Physiological Responses to Climate Change in Semi-Arid Areas of North China

1
Institute of Ecological Conservation and Restoration, Chinese Academy of Forestry, Beijing 100091, China
2
Key Laboratory of Desert Ecosystem and Global Change, State Administration of Forestry and Grassland, Beijing 100091, China
3
Key Laboratory of State Forestry Administration on Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China
4
School of Modern Language, University of Bristol, Bristol BS8 1QU, UK
*
Author to whom correspondence should be addressed.
Forests 2024, 15(7), 1272; https://doi.org/10.3390/f15071272
Submission received: 21 June 2024 / Revised: 9 July 2024 / Accepted: 20 July 2024 / Published: 22 July 2024
(This article belongs to the Special Issue Construction and Maintenance of Desert Forest Plantation)

Abstract

:
Climate change has had a widespread and profound impact on global temperature and precipitation patterns, especially in semi-arid areas. Plant δ13C and iWUE indicate the trade-off between carbon uptake and water loss, which is pivotal for understanding plant responses to climate change. Information about the long-term responses of the physiological and ecological processes of different tree species to climate change is also required. To investigate the impact of different forest stand structures and site conditions on long-term growth and physiological processes of coniferous and broad-leaved trees in the mountainous area of Beijing, we analyzed the tree-ring δ13C variation of four tree species (Platycladus orientalis, Pinus tabuliformis, Quercus variabilis, Robinia pseudoacacia) sampled from 64 plots with varying site and stand conditions. We found that the tree-ring δ13C of the four tree species varied from each other and was mainly affected by density and slope aspect, followed by slope and age. Both tree-ring δ13C and iWUE of the four tree species showed increasing trends over time, mechanistically linked to long-term changes in global CO2 concentration. This indicates the four native tree species have adapted well to climate change, and the risk of decline is relatively low. The increased iWUE translated into different growth patterns which varied with tree species, site, and stand condition. Different tree species have varying sensitivities to environmental factors. The iWUE of coniferous tree species is more sensitive to climate change than that of broad-leaved tree species, especially to temperature (T), the Standardized Precipitation Evapotranspiration Index (SPEI), and vapor pressure deficit (VPD).

1. Introduction

In the past 100 years, the global average air temperature has increased significantly, and the rate of this increase is accelerating, indicating that the rising trend of air temperature will not weaken in the short term [1,2]. Climate change has had a widespread and profound impact on global temperature and precipitation patterns [3]. The “fertilization effect” caused by the increase in CO2 concentration, combined with the opposing effect of drought stress inhibiting tree growth, leads to great uncertainty in growth patterns [4,5]. As the main component of the terrestrial ecosystem, understanding the physiological responses of different tree species in this situation is of great significance for predicting tree growth patterns and the process of tree decline in such climate-sensitive areas. Photosynthesis and transpiration in trees are two important ecological processes that are interconected. The increase in atmospheric CO2 concentration can promote tree growth, but rising temperatures and drought have a negative impact on the forest ecosystem. Under long-term climate change, understanding the response of tree water and carbon processes is of great significance for the integrated management and regulation of regional water and carbon resources. This is especially true in North China, which has become a high-risk area for drought due to the intensification of climate warming [3,6]. Understanding the physiological responses of different tree species to climate change is crucial for understanding and predicting vegetation growth patterns, tree decline processes, and the driving factors in climate sensitive areas.
Carbon isotope effects exist in the process of plants absorbing and utilizing atmospheric CO2, influenced by photosynthetic pathways (C3, C4, and CAM) and environmental factors [7]. Plant δ13C reflects the interactions between various aspects of plant carbohydrate relationships, and can be used more effectively as a comprehensive indicator of whole-plant function compared to gas exchange measurements [8]. It can also reveal how species survive in specific habitats and ensure competitive advantage by regulating their gas exchange processes, resource acquisition, and utilization strategies [9,10]. The dynamic fractionation of 13C during plant photosynthesis is influenced by individual differences and environmental factors [7]. Stable isotopes in tree-rings are a valuable tool for studying the physiological and ecological mechanisms induced by drought, and are often used to reconstruct past changes in temperature and precipitation during growing seasons [11,12]. However, their climate response is still not fully understood, particularly in non-extreme mid-latitude environments with diverse ecological conditions. Tree-ring stable carbon isotope parameters reflect the response of trees to climate change over a long period of time. By tracing the physiological and ecological responses of trees during their growth process, and analyzing the differences between different tree species, we can reveal how different tree species respond to climate change [13]. As an important indicator for evaluating tree growth efficiency, intrinsic water-use efficiency (iWUE) directly reflects the impact of global climate change on the ecosystems [14]. The iWUE can be inferred through δ13C, which simultaneously records the carbon assimilation and stomatal conductance, both of which are susceptible to environmental conditions [15,16]. In the context of climate change, some forests have declined, while others have adapted to changing environmental conditions. However, the underlying mechanisms remain unclear. The information about the various physical and physiological mechanisms of iWUE and tree-ring δ13C change as tree growth progresses is also required.
In order to clarify the long-term physiological and ecological responses of different tree species to climate change, we investigated the impact of different forest stand structures and site conditions on long-term growth and physiological processes of coniferous and broad-leaved trees in semi-arid North China, using tree-ring stable isotopes. Furthermore, we discuss the adaptability of four native tree species to climate change.

2. Methods and Materials

2.1. Study Area

The study area is located in Beijing, China (39°28′–41°05′ N, 115°25′–117°30′ E). This area experiences a temperate continental climate. The average annual precipitation is 500–600 mm, while the average annual evaporation is about 1000 mm. The majority of precipitation occurs in summer, accounting for 64%–69% of the annual total. The average annual temperature is 10–12 °C. The geomorphology of the study area is complex, with vegetation covering more than half of the total area, playing an important role in ecological protection. There are various tree species in the study area. The distribution areas of each species are approximately 86,000 hm2 of Platycladus orientalis, 77,000 hm2 of Pinus tabuliformis, 53,000 hm2 of Quercus variabilis, and 15,000 hm2 of Robinia pseudoacacia.

2.2. Sample Plot and Tree-Ring Collection

Two coniferous species, P. orientalis and P. tabulaeformis, and two broad-leaved species, Q. variabilis and R. pseudoacacia, were selected for the study, which are widely distributed naturally in North China and have a long history of growth. These trees have strong adaptability to soil with no strict requirements for texture, and they are tolerant of drought and barren conditions. The environmental factors affecting δ13C are categorized into stand factors (including age and density), site factors (including slope and aspect), and climate factors (including precipitation, air temperature, and relative humidity). Due to the thin soil layer in the Beijing mountainous area, and the small altitude variation (400–1000 m) in the main distribution areas of the studied tree species, further subdivision was not performed. Changes in soil water content are mainly influenced by precipitation, slope, and aspect, with precipitation having the most significant effect. Therefore, the response to soil water content is encompassed within the response to slope, aspect, and precipitation, and is thus no longer discussed separately.
Sixteen plot types for each tree species were selected in this study, resulting in a total of 64 sample plots (Figure 1). The specific classifications are shown in Table 1. Plot surveys and log scaling were carried out on all plots, recording the basic information of each 20 m × 20 m plot and the data of each standing tree. Within each plot, four trees were selected for sampling. Two cores were taken from each tree at breast height (1.3 m) along the east–west and north–south directions, respectively, using a 5 mm increment borer. Except for the damaged tree cores, a total of 512 samples were collected (133 cores for P. orientalis, 137 cores for P. tabulaeformis, 123 cores for Q. variabilis, and 119 cores for R. pseudoacacia). To prevent contamination from other carbon sources, the collected samples were stored in glass tubes.

2.3. Tree-Ring Stable Carbon Isotope Analysis

The core samples were sanded with grain paper of varying frits from 100 to 800 mesh to make the tree-rings more clearly visible for cross-dating. COFECHA (version XP 2007) was then used to verify dating, identify errors in the individual series, and improve the cross-dating accuracy [17]. Prior to the isotope analysis, the first 10 years of each core were removed to avoid “juvenile age effects” on the tree-ring isotopes. Therefore, the tree-rings from 1960 to 2020 were ultimately preserved for isotope analysis. All tree cores from each plot were used for stable carbon isotope analysis after the establishing the tree-ring chronology. The cores were separated into individual rings using a scalpel to cut along the ring lines under a stereomicroscope (40× magnification). Rings from the same year were pooled to form a composite sample representing the plot-aggregated isotopic composition changes [18,19]. Whole wood was used for the δ13C analysis in this study to trace temporal trends and responses to climate variables, based on the strong correlation between the carbon sources of whole wood and cellulose [20,21]. Each pooled sample was ground into a powder using a ball mill, sieved through an 80-mesh sieve, and then the milligram wood powder was placed into tin capsules for δ13C analyses using mass spectrometry (DELTA V Advantage; Thermo Fisher Scientific Inc., Waltham, MA, USA; precision 0.1‰).

2.4. Carbon Isotope Discrimination and iWUE

The carbon isotope effect in plant photosynthesis can usually be expressed by the following formula [22]:
δ 13 C p = δ 13 C a a b a × C i C a
where δ13Ca and δ13Cp represent the ratios in atmospheric CO2 and tree-rings, respectively; a and b represent the isotope fractionation coefficient of the CO2 diffusion process (4.4‰) [23] and the carbon isotope fractionation coefficient of the photosynthetic carboxylase carboxylation process (27‰) [24]; Ci and Ca are the intercellular and ambient CO2 concentrations, respectively. The isotope effect can also be expressed by isotope discrimination value (Δ):
Δ = δ 13 C a δ 13 C p 1 + δ 13 C p
Therefore, Δ is also related to the ratio of the Ci to Ca mole fraction [24] as:
Δ = a + ( b a ) × C i C a
iWUE = A G s = C a C i 1.6
iWUE (μmol mol−1) is defined as the ratio of net carbon assimilation (A) to stomatal conductance to water vapor (gw). The linear relationship between Ci/Ca and Δ can be used to calculate the iWUE [25]:
iWUE = C a b Δ 1.6 b a

2.5. Climate Data

Annual meteorological data (precipitation, temperature, and relative humidity) of the study area from 1960 to 2020 were downloaded from the National Meteorological Science Data Center (http://data.cma.cn/, accessed on 12 June 2023). Kendall tests were conducted to confirm the data validity and representativeness of the local climatic trends. The Standardized Precipitation Evapotranspiration Index (SPEI) was calculated using the SPEI package in R 3.5.1 (http://sac.csic.es/spei/, accessed on 5 September 2023) and was used for drought assessment and analysis. Vapor pressure deficit (VPD), an important index reflecting the degree of atmospheric drought, was calculated according to the general formula [16].

2.6. Data Analysis

One-way analysis of variance (ANOVA) was used to test for significant differences in tree-ring δ13C between different plots and different tree species. Random forest analysis was used to identify the effects of site conditions on tree-ring δ13C by constructing a random forest model incorporating tree-ring δ13C and various factors. Regression analysis was used to explore the trends in tree-ring δ13C and iWUE. Pearson’s correlations were used to indicate the influence of precipitation, temperature, relative humidity, SPEI, and VPD on tree-ring δ13C and iWUE.

3. Results

3.1. Long-Term Tree-Ring δ13C of Different Tree Species

The average value of tree-ring δ13C for the four tree species ranges from −25.2‰ to −23.2‰, with significant differences among species. The order is as follows: δ13CP.o (−23.963‰~−23.054‰) > δ13CP.t (−24.340‰~−23.737‰) > δ13CQ.v (−24.718‰~−24.203‰) > δ13CR.p (−25.187‰~−24.794‰) (Figure 2). The annual variations in tree-ring δ13C of the four tree species between 1960 and 2020, shown in Figure 3, indicate that the average tree-ring δ13C of coniferous trees is higher than that of broad-leaved tree species during the study period. Although the growth trend of δ13CP.o has weakened in the past decade, it still shows a significant gradual increase, while the trends for other tree species are relatively gentle. Additionally, δ13CP.o exhibits the largest fluctuation range, with a larger standard error than other tree species, indicating more dispersed average values.

3.2. Comparison of Tree-Ring δ13C among Plots and Species

Differences in tree-ring δ13C among plots are shown in Figure 4, revealing both similarities and differences among the tree species. The density of plots P1~P8 is relatively high, while the density of plots P9~P16 is relatively low. According to the numerical orders of δ13C among different plots for four tree species shown in Table 2, we found that the δ13C of high-density sample plots (75%~87.5%) is higher than that of low-density plots, suggesting that density has a positive effect on tree-ring δ13C. Plots P1~P4 and P9~P12 are relatively older stands, while plots P5~P8 and P13~P16 are relatively younger stands. The average ranking shows that tree-ring δ13CP.o and δ13CP.t do not exhibit a clear pattern with age, whereas 62.8% of δ13CQ.v and δ13CR.p values are higher in older trees compared to younger ones. Plots P1, P2, P5, P6, P11, P12, P15, and P16 are steep slopes, and the rest are on gentle slopes. There is no obvious regularity shown in the rank of different slopes. Plots P1, P3, P5, P7, P9, P11, P13, and P15 are on sunny slopes, while the rest are on shady slopes. Tree-ring δ13C values for 62.5%~75% of the sample plots on sunny slopes are higher than those on shady slopes.
Random forest analysis was used to determine the influence of density, aspect, slope, and age on tree-ring δ13C of P. orientalis, P. tabulaeformis, Q. variabilis, and R. pseudoacacia (Figure 5). The results show that the tree-ring δ13C of four tree species was mainly affected by density and aspect, followed by slope and age. There are also differences in the importance ranking of these factors for different tree species. For example, tree-ring δ13C of P. orientalis is most affected by aspect, while others are most affected by density. In addition, tree age has a negative effect on tree-ring δ13C of P. orientalis, P. tabulaeformis, and Q. variabilis, but a positive effect on R. pseudoacacia.

3.3. Relations of Tree-Ring δ13C and iWUE with Climate Index

There are significant differences in Δ among the four tree species, with ΔR.p > ΔQ.v > ΔP.t > ΔP.o (Figure 6). Decreasing trends in Δ were observed over time for stands of all tree species, with approximately downward-slopes of the curves. The intrinsic water-use efficiency (iWUE) of the four tree species, derived from tree-ring δ13C measurements (Figure 6), showed similar increasing trends. Both Ci and iWUE show similar increasing trends, with the ranges of Ci being CiR.p > CiQ.v > CiP.t > CiP.o and iWUE being iWUEP.o > iWUEP.t > iWUEQ.v > iWUER.p. The growth rates of iWUE for the four tree species are similar, with iWUEP.o having the largest growth rate at 0.88 μmol mol−1 per year.
Meteorological data displayed significant fluctuations in precipitations and rising air temperatures between 1960 and 2020. No significant increasing or decreasing trend was found for VPD, but the VPD anomaly exhibited increasing occurrences of higher VPD. In addition, the frequency of negative SPEI increased. Correlation analysis of the averaged tree-ring δ13C from 16 sample plots of four tree species with meteorological factors and drought shows that tree-ring δ13C has negative correlations with P and RH, and positive correlations with T, SPEI, and VPD (Figure 7). Overall, the tree-ring δ13C of coniferous tree species is more sensitive to climate change than that of broad-leaved tree species. Comparison among the four tree species reveals that δ13CP.o is the most sensitive to climate change, showing a significant positive correlation with T, SPEI, and VPD (p < 0.05), and a significant negative correlation with P (p < 0.05). The sensitivity of the tree-ring δ13C responses of the four tree species to P is significantly different, while the sensitivity to T is slightly different. There are no significant differences in the correlations between tree-ring δ13C and RH.
Similar results are found in the correlation of the averaged iWUE of the four tree species with meteorological factors and drought, showing that iWUE has negative correlations with P and RH, and positive correlations with T, SPEI, and VPD (Figure 8). It is also found that the iWUE of coniferous tree species is more sensitive to climate change than broad-leaved tree species, especially to T, SPEI, and VPD. The responses of iWUE to P are relatively weak compared to those of tree-ring δ13C. There are no significant differences in the correlations between iWUE and P, or between iWUE and RH.

4. Discussion

4.1. Effects of Site Conditions, Stand Factors, Climate Change, and Drought on Tree-Ring δ13C

Tree-ring stable isotopes can exhibit robust climate signals [26,27]. However, understanding the influence of various parameters on tree-ring δ13C is challenging due to interactions among environmental and stand parameters [28]. The analysis results indicate that tree-ring δ13C in four tree species is primarily affected by the density and slope aspect, followed by slope and age. Light and water availability fluctuate in the dynamic field environment. Understory trees face more growth limitations due to drought [29]. The sampling sites differ in soil moisture and light, primarily controlled by slope differences and aspect [11]. This site appeared to be the driest due to the shallow soil, steep slope, and southern exposure. Because drought and shade have contrasting effects on isotope discrimination, it is necessary to distinguish between water and light limitations on tree growth. Our study shows no significant differences in tree-ring δ13C between trees growing on steep slopes and those on gentle slopes, suggesting that within the sampling slope range, the physiological processes of the four tree species were not affected by different slopes. There may be two possible reasons for this phenomenon. First, the sampled trees did not distinguish between uphill and downhill slopes. Second, the gentle slope may not have an impact on soil moisture. Compared to slope degree, the impact of slope direction is more significant in this study. CO2 assimilation is greater on sunny slopes compared to shady slopes [30]. Trees growing on sunny slopes have significantly higher tree-ring δ13C than those growing on shady slopes. Similarly, Chen et al. [31] found that plants grown on the sunny slopes were enriched with more 13C compared to those on the shady slopes, attributed to the drier and warmer climate on the sunny slopes relative to the shady slopes.
Stand density influences the intensity of interspecific and intraspecific competition in forests, thereby affecting water-use strategies and tree growth [32,33]. However, there is limited understanding of the relationship between stand density, age, droughts, stable carbon discrimination, and tree growth in plantations [34,35]. Tree-ring δ13C records signatures influenced by neighborhood effects, where the structure of tree canopies is shaped by stand density and functional traits of neighboring trees [36,37]. In this study, we found that 75% to 87.5% of high-density plots had higher tree-ring δ13C than low-density plots. Higher tree densities increase competition for resources (light, water, and CO2) in the ecosystem. Trees in this forest type enhance water-use efficiency, leading to higher tree-ring δ13C values and increased iWUE. Gouveia and Freitas [32] found stand density-dependent differences in leaf carbon isotope discrimination in a water-limited oak forest. They identified an optimal stand density by the comparing of δ13C signatures among stands with varying stem densities. Nevertheless, it is still unclear if tree-ring stable isotopes contain age-related trends. In addition, it is important to understand the changes of climate response in tree-ring data of different tree ages, yet there is limited knowledge on the impact of tree age on the climate response of tree-ring stable isotopes. While Gagen et al. [38] observed signs of the juvenile growth effect of δ13C in Finnish pine trees, they found no long-term age-related trend in tree-ring stable carbon isotopes. Likewise, Daux et al. [39] did not observe age-related low-frequency behavior in the carbon isotope characteristics of Larix decidua from the French Alps, consistent with findings from live pine and residual pine in northern Sweden [40]. In addition, Young et al. [41] reported that there was no long-term age-related trend in the stable carbon and oxygen isotopes of living pine trees grown in the marine conditions in Northwest Norway. Xu et al. [42] measured tree-ring δ13C in three age groups of Schrenk spruce (Picea schrenkiana) trees in northwestern China to determine if age differences in tree-ring δ13C could lead to varying climate responses. They found that young trees exhibited a stronger climate response in δ13C compared to old trees.

4.2. Tree-Ring δ13C and iWUE Reflect Long-Term Tree Responses to Environment

The precipitation pattern and air temperature in North China are changing, manifested by significant fluctuations in precipitation and rising air temperatures. The increased frequency of negative SPEI values indicates the intensification of drought severity. iWUE reflects the balance between the carbon uptake and water loss, influenced by atmospheric CO2 concentration and water availability. When exposed to long-term water constraints, trees improve their iWUE to reduce transpiration water loss. This is achieved through strict stomatal regulation to maintain the hydraulic function of the xylem, or through stomatal closure due to the loss of hydraulic conductivity. This understanding is crucial for interpreting plant responses to climate change [43,44]. It is important to assess the effects of climate change, such as increasing atmospheric CO2 concentration and water stress, on carbon sequestration and iWUE [6,45]. A meta-analysis revealed a ~40% increase in tree iWUE globally since 1901, paralleling a ∼34% increase in atmospheric CO2. However, mean iWUE, and its rate of increase varied across biomes and leaf and wood functional types [46]. We found clear trends over time in tree-ring δ13C and iWUE for four native tree species, mechanistically linked to long-term changes in global CO2 concentration, which is consistent with previous studies [18,47,48,49]. The long-term growth of iWUE indicates an increase in net photosynthesis and/or a decrease in stomatal conductance (and consequently transpiration). The prolonged rise in atmospheric CO2 concentration is expected to promote both photosynthesis and reduce transpiration [50]. Drought often enhances iWUE, improving the adaptability of trees by enhancing carbon assimilation and reducing water loss. However, worsening drought stress does not always lead to increased iWUE. In the American Southwest, a severe, multidecadal drought has significantly increased plant water-use efficiency and has reached a critical threshold [43]. Currently, a section of the Three-North Shelterbelt Forest in China has experienced decline, particularly in poplar plantations [19,51]. Despite extensive research on the decline of plantations in arid and semi-arid areas, the understanding of the underlying mechanisms remains unclear [52,53]. Integrating tree-ring data and iWUE to elucidate the growth process of trees is crucial for understanding both the mechanisms of decline and the adaptability of trees to future climate change.
Increased iWUE has significantly positive impacts on tree recovery (17%) and resilience (15%), indicating enhanced restoration capabilities [54]. Our study shows that tree-ring δ13C and iWUE are still increasing, suggesting that the physiological parameters of the four tree species have not yet peaked. This performance indicates that these native tree species have adapted well to local environmental changes, with a relatively low risk of decline. The increased iWUE has resulted in varying growth patterns depending on tree species, site, and stand conditions [46,55]. Tree species exhibit varying sensitivities to environmental factors. Our study also revealed that coniferous tree species are more sensitive to climate change in terms of iWUE, particularly to T, SPEI, and VPD, compared to broad-leaved species. Both shrub and tree iWUE showed greater sensitivity to variability in atmospheric aridity than to increasing atmospheric CO2 concentration [44]. Liu et al. [56] determined that temperature and precipitation were the primary factors influencing water-use efficiency in trees growing in seasonal drought areas, based on correlation analyses between needle δ13C values and meteorological data. They found that temperature had a stronger impact on water-use efficiency in middle-aged forests. Identifying and selecting tree species with high water-use efficiency is crucial for sustaining forest benefits under changing climate conditions and limited water supply [57].

5. Conclusions

Our study reveals distinct differences in the tree-ring δ13C among four native tree species, primarily influenced by density and slope aspect, followed by slope and age. Over time, both tree-ring δ13C and iWUE have shown consistent increases, correlating with global CO2 concentration changes. Different tree species exhibit varying sensitivities to environmental factors. Coniferous species, for instance, demonstrate higher sensitivity to temperature, the Standardized Precipitation Evapotranspiration Index (SPEI), and vapor pressure deficit (VPD) compared to broad-leaved species. Tree-ring δ13C and iWUE serve as integrated indicators of environmental variability and the impacts of stand and site conditions on tree physiology and growth. Our findings underscore the adaptive capacity of these native tree species in response to changing environmental conditions. The increasing trends in iWUE highlight their ability to optimize water-use efficiency amidst rising atmospheric CO2 levels.

Author Contributions

Conceptualization, W.L. and X.Y.; methodology, G.J.; software, Y.G.; validation, W.L. and L.W.; formal analysis, W.L.; investigation, W.L.; resources, B.W.; data curation, W.L.; writing—original draft preparation, W.L.; writing—review and editing, A.L.; supervision, A.L.; project administration, B.W.; funding acquisition, B.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key R&D Program of China [2023YFF1305304] and National Natural Science Foundation of China [32001372].

Data Availability Statement

The data that support the findings of this study are available on request from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lucarini, V.; Ragone, F.; Lunkeit, F. Predicting climate change using response theory: Global averages and spatial patterns. J. Stat. Phys. 2017, 166, 1036–1064. [Google Scholar] [CrossRef]
  2. Devi, N.M.; Kukarskih, V.V.; Galimova, A.A.; Mazepa, V.S.; Grigoriev, A.A. Climate change evidence in tree growth and stand productivity at the upper treeline ecotone in the Polar Ural Mountains. For. Ecosyst. 2020, 7, 7. [Google Scholar] [CrossRef]
  3. IPCC. Climate Change 2021: The Physical Science Basis; IPCC: Geneva, Switzerland, 2021. [Google Scholar]
  4. Zhao, M.; Running, S.W. Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. Science 2010, 329, 940–943. [Google Scholar] [CrossRef] [PubMed]
  5. Boysen, L.R.; Lucht, W.; Gerten, D.; Heck, V.; Lenton, T.M.; Schellnhuber, H.J. The limits to global-warming mitigation by terrestrial carbon removal. Earth’s Future 2017, 5, 463–474. [Google Scholar] [CrossRef]
  6. Lu, W.; Yu, X.; Jia, G.; Li, H.; Liu, Z. Responses of Intrinsic Water-use Efficiency and Tree Growth to Climate Change in Semi-Arid Areas of North China. Sci. Rep. 2018, 8, 308. [Google Scholar]
  7. Farquhar, G.D.; Ehleringer, J.R.; Hubick, K.T. Carbon isotope discrimination and photosynthesis. Annu. Rev. Plant Biol. 1989, 40, 503–537. [Google Scholar] [CrossRef]
  8. Dawson, T.E.; Mambelli, S.; Plamboeck, A.H.; Templer, P.H.; Tu, K.P. Stable isotopes in plant ecology. Annu. Rev. Ecol. Syst. 2002, 33, 507–559. [Google Scholar] [CrossRef]
  9. Ehleringer, J.R. Variation in leaf carbon isotope discrimination in Encelia farinosa: Implications for growth, competition, and drought survival. Oecologia 1993, 95, 340–346. [Google Scholar] [CrossRef]
  10. Ehleringer, J.R.; Cooper, T.A. Correlations between carbon isotope ratio and microhabitat in desert plants. Oecologia 1988, 76, 562–566. [Google Scholar] [CrossRef]
  11. Hartl-Meier, C.; Zang, C.; Büntgen, U.; Esper, J.; Rothe, A.; Göttlein, A.; Dirnböck, T.; Treydte, K. Uniform climate sensitivity in tree-ring stable isotopes across species and sites in a mid-latitude temperate forest. Tree Physiol. 2015, 35, 4–15. [Google Scholar] [CrossRef]
  12. Eilmann, B.; Buchmann, N.; Siegwolf, R.; Saurer, M.; Cherubini, P.; Rigling, A. Fast response of Scots pine to improved water availability reflected in tree-ring width and δ13C. Plant Cell Environ. 2010, 33, 1351–1360. [Google Scholar] [CrossRef] [PubMed]
  13. Marchand, W.; Girardin, M.P.; Hartmann, H.; Depardieu, C.; Bergeron, Y. Strong overestimation of water-use efficiency responses to rising CO2 in tree-ring studies. Glob. Chang. Biol. 2020, 26, 4538–4558. [Google Scholar] [CrossRef] [PubMed]
  14. Sidorova, O.V.C.; Fonti, M.V.; Siegwolf, R.T.W.; Saurer, M.; Myglan, V.S. Impact of Recent Climate Change on Water-Use Efficiency Strategies of Larix sibirica in the Altai-Sayan Mountain Range. Forests 2020, 11, 1103. [Google Scholar] [CrossRef]
  15. Timofeeva, G.; Treydte, K.; Bugmann, H.; Rigling, A.; Schaub, M.; Siegwolf, R.; Saurer, M. Long-term effects of drought on tree-ring growth and carbon isotope variability in Scots pine in a dry environment. Tree Physiol. 2017, 37, 1028–1041. [Google Scholar] [CrossRef]
  16. Jia, G.; Chen, L.; Yu, X.; Liu, Z. Soil water stress overrides the benefit of water-use efficiency from rising CO2 and temperature in a cold semi-arid poplar plantation. Plant Cell Environ. 2022, 45, 1172–1186. [Google Scholar] [CrossRef]
  17. Nock, C.A.; Baker, P.J.; Wanek, W.; Leis, A.; Grabner, M.; Bunyavejchewin, S.; Hietz, P. Long-term increases in intrinsic water-use efficiency do not lead to increased stem growth in a tropical monsoon forest in western Thailand. Glob. Change Biol. 2011, 17, 1049–1063. [Google Scholar] [CrossRef]
  18. Liu, X.; Wang, W.; Xu, G.; Zeng, X.; Wu, G.; Zhang, X.; Qin, D. Tree growth and intrinsic water-use efficiency of inland riparian forests in northwestern China: Evaluation via δ13C and δ18O analysis of tree rings. Tree Physiol. 2014, 34, 966–980. [Google Scholar] [CrossRef]
  19. Lu, W.; Yu, X.; Jia, G. Retrospective analysis of tree decline based on intrinsic water-use efficiency in semi-arid areas of north China. Atmosphere 2020, 11, 577. [Google Scholar] [CrossRef]
  20. Loader, N.; Robertson, I.; McCarroll, D. Comparison of stable carbon isotope ratios in the whole wood, cellulose and lignin of oak tree-rings. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2003, 196, 395–407. [Google Scholar] [CrossRef]
  21. Sun, S.; He, C.; Qiu, L.; Li, C.; Zhang, J.; Meng, P. Stable isotope analysis reveals prolonged drought stress in poplar plantation mortality of the Three-North Shelter Forest in Northern China. Agric. For. Meteorol. 2018, 252, 39–48. [Google Scholar] [CrossRef]
  22. Farquhar, G.D.; Sharkey, T.D. Stomatal conductance and photosynthesis. Annu. Rev. Plant Physiol. 1982, 33, 317–345. [Google Scholar] [CrossRef]
  23. O’Leary, M.H. Carbon isotope fractionation in plants. Phytochemistry 1981, 20, 553–567. [Google Scholar] [CrossRef]
  24. Farquhar, G.; Richards, R. Isotopic composition of plant carbon correlates with water-use efficiency of wheat genotypes. Funct. Plant Biol. 1984, 11, 539–552. [Google Scholar] [CrossRef]
  25. McCarroll, D.; Loader, N.J. Stable isotopes in tree rings. Quat. Sci. Rev. 2004, 23, 771–801. [Google Scholar] [CrossRef]
  26. Cherubini, P.; Battipaglia, G.; Innes, J.L. Tree vitality and forest health: Can tree-ring stable isotopes be used as indicators? Curr. For. Rep. 2021, 7, 69–80. [Google Scholar] [CrossRef]
  27. Panek, J.A.; Waring, R.H. Carbon isotope variation in Douglas-fir foliage: Improving the δ13C–climate relationship. Tree Physiol. 1995, 15, 657–663. [Google Scholar] [CrossRef]
  28. Li, X.; Chen, G.; Qin, W.; Wang, X.; Wang, P. Differences in responses of tree-ring δ13C in angiosperms and gymnosperms to climate change on a global scale. For. Ecol. Manag. 2021, 492, 119247. [Google Scholar] [CrossRef]
  29. Brienen, R.; Helle, G.; Pons, T.; Boom, A.; Gloor, M.; Groenendijk, P.; Clerici, S.; Leng, M.; Jones, C. Paired analysis of tree ring width and carbon isotopes indicates when controls on tropical tree growth change from light to water limitations. Tree Physiol. 2022, 42, 1131–1148. [Google Scholar] [CrossRef]
  30. Lu, W.; Yu, X.; Jia, G. Instantaneous and long-term CO2 assimilation of Platycladus orientalis estimated from 13C discrimination. Ecol. Indic. 2019, 104, 237–247. [Google Scholar] [CrossRef]
  31. Chen, Z.; Wang, G.; Jia, Y. Foliar δ13C showed no altitudinal trend in an arid region and atmospheric pressure exerted a negative effect on plant δ13C. Front. Plant Sci. 2017, 8, 250867. [Google Scholar] [CrossRef]
  32. Gouveia, A.C.; Freitas, H. Intraspecific competition and water use efficiency in Quercus suber: Evidence of an optimum tree density? Trees 2008, 22, 521–530. [Google Scholar] [CrossRef]
  33. Jack, S.B.; Long, J.N. Analysis of stand density effects on canopy structure: A conceptual approach. Trees 1991, 5, 44–49. [Google Scholar] [CrossRef]
  34. Toillon, J.; Fichot, R.; Dallé, E.; Berthelot, A.; Brignolas, F.; Marron, N. Planting density affects growth and water-use efficiency depending on site in Populus deltoides × P. nigra. For. Ecol. Manag. 2013, 304, 345–354. [Google Scholar] [CrossRef]
  35. Zhang, X.; Yu, P.; Wang, D.; Xu, Z. Density-and age-dependent influences of droughts and intrinsic water use efficiency on growth in temperate plantations. Agric. For. Meteorol. 2022, 325, 109134. [Google Scholar] [CrossRef]
  36. Mölder, I.; Leuschner, C.; Leuschner, H.H. δ13C signature of tree rings and radial increment of Fagus sylvatica trees as dependent on tree neighborhood and climate. Trees 2011, 25, 215–229. [Google Scholar] [CrossRef]
  37. McDowell, N.G.; Allen, C.D.; Marshall, L. Growth, carbon-isotope discrimination, and drought-associated mortality across a Pinus ponderosa elevational transect. Glob. Change Biol. 2010, 16, 399–415. [Google Scholar] [CrossRef]
  38. Gagen, M.; McCarroll, D.; Robertson, I.; Loader, N.J.; Jalkanen, R. Do tree ring δ13C series from Pinus sylvestris in northern Fennoscandia contain long-term non-climatic trends? Chem. Geol. 2008, 252, 42–51. [Google Scholar] [CrossRef]
  39. Daux, V.; Édouard, J.-L.; Masson-Delmotte, V.; Stievenard, M.; Hoffmann, G.; Pierre, M.; Mestre, O.; Danis, P.; Guibal, F. Can climate variations be inferred from tree-ring parameters and stable isotopes from Larix decidua? Juvenile effects, budmoth outbreaks, and divergence issue. Earth Planet. Sci. Lett. 2011, 309, 221–233. [Google Scholar] [CrossRef]
  40. Loader, N.; Young, G.; Grudd, H.; McCarroll, D. Stable carbon isotopes from Torneträsk, northern Sweden provide a millennial length reconstruction of summer sunshine and its relationship to Arctic circulation. Quat. Sci. Rev. 2013, 62, 97–113. [Google Scholar] [CrossRef]
  41. Young, G.H.; Demmler, J.C.; Gunnarson, B.E.; Kirchhefer, A.J.; Loader, N.J.; McCarroll, D. Age trends in tree ring growth and isotopic archives: A case study of Pinus sylvestris L. from northwestern Norway. Glob. Biogeochem. Cycles 2011, 25, GB2020. [Google Scholar] [CrossRef]
  42. Xu, G.; Wu, G.; Liu, X.; Chen, T.; Wang, B.; Hudson, A.; Trouet, V. Age-related climate response of tree-ring δ13C and δ18O from spruce in northwestern China, with implications for relative humidity reconstructions. J. Geophys. Res. Biogeosci. 2020, 125, e2019JG005513. [Google Scholar] [CrossRef]
  43. Kannenberg, S.A.; Driscoll, A.W.; Szejner, P.; Anderegg, W.R.; Ehleringer, J.R. Rapid increases in shrubland and forest intrinsic water-use efficiency during an ongoing megadrought. Proc. Natl. Acad. Sci. USA 2021, 118, e2118052118. [Google Scholar] [CrossRef] [PubMed]
  44. Lu, H.; Qin, Z.; Lin, S.; Chen, X.; Chen, B.; He, B.; Wei, J.; Yuan, W. Large influence of atmospheric vapor pressure deficit on ecosystem production efficiency. Nat. Commun. 2022, 13, 1653. [Google Scholar] [CrossRef] [PubMed]
  45. Metz, J.; Annighöfer, P.; Schall, P.; Zimmermann, J.; Kahl, T.; Schulze, E.D.; Ammer, C. Site-adapted admixed tree species reduce drought susceptibility of mature European beech. Glob. Change Biol. 2016, 22, 903–920. [Google Scholar] [CrossRef] [PubMed]
  46. Mathias, J.M.; Thomas, R.B. Global tree intrinsic water use efficiency is enhanced by increased atmospheric CO2 and modulated by climate and plant functional types. Proc. Natl. Acad. Sci. USA 2021, 118, e2014286118. [Google Scholar] [CrossRef] [PubMed]
  47. Tang, Y.; Sahlstedt, E.; Young, G.; Schiestl-Aalto, P.; Saurer, M.; Kolari, P.; Jyske, T.; Bäck, J.; Rinne-Garmston, K.T. Estimating intraseasonal intrinsic water-use efficiency from high-resolution tree-ring δ13C data in boreal Scots pine forests. New Phytol. 2023, 237, 1606–1619. [Google Scholar] [CrossRef] [PubMed]
  48. Rezaie, N.; D’Andrea, E.; Bräuning, A.; Matteucci, G.; Bombi, P.; Lauteri, M. Do atmospheric CO2 concentration increase, climate and forest management affect iWUE of common beech? Evidences from carbon isotope analyses in tree rings. Tree Physiol. 2018, 38, 1110–1126. [Google Scholar] [CrossRef] [PubMed]
  49. Liu, Z.; Ji, X.; Ye, L.; Jiang, J. Inherent water-use efficiency of different forest ecosystems and its relations to climatic variables. Forests 2022, 13, 775. [Google Scholar] [CrossRef]
  50. van der Sleen, J.P.; Groenendijk, P.; Vlam, M.; Anten, N.P.R.; Boom, A.; Bongers, F.; Pons, T.L.; Terburg, G.; Zuidema, P.A. No growth stimulation of tropical trees by 150 years of CO2 fertilization but water-use efficiency increased. Nat. Geosci. 2015, 8, 24–28. [Google Scholar] [CrossRef]
  51. Sun, S.; Qiu, L.; He, C.; Li, C.; Zhang, J.; Meng, P. Drought-Affected Populus simonii Carr. Show Lower Growth and Long-Term Increases in Intrinsic Water-Use Efficiency Prior to Tree Mortality. Forests 2018, 9, 564. [Google Scholar] [CrossRef]
  52. Allen, C.D.; Macalady, A.K.; Chenchouni, H.; Bachelet, D.; Mcdowell, N.; Vennetier, M.; Kitzberger, T.; Rigling, A.; Breshears, D.D.; Hogg, E.H.T. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For. Ecol. Manag. 2010, 259, 660–684. [Google Scholar] [CrossRef]
  53. Mcdowell, N.G. Mechanisms Linking Drought, Hydraulics, Carbon Metabolism, and Vegetation Mortality. Plant Physiol. 2011, 155, 1051–1059. [Google Scholar] [CrossRef]
  54. Liu, X.; Wu, G.; Xu, G.; Chen, T.; Wang, B. The positive contribution of iWUE to the resilience of Schrenk spruce (Picea schrenkiana) to extreme drought in the western Tianshan Mountains, China. Acta Physiol. Plant 2020, 42, 168. [Google Scholar]
  55. Martínez-Sancho, E.; Dorado-Lián, I.; Gutiérrez-Merino, E.; Matiu, M.; Helle, G.; Heinrich, I.; Menzel, A. Increased water-use efficiency translates into contrasting growth patterns of Scots pine and sessile oak at their southern distribution limits. Glob. Change Biol. 2017, 24, 1012–1028. [Google Scholar] [CrossRef]
  56. Liu, Y.; Wu, J.; Wu, D.; Li, S.; Wang, L. Seasonal variation in δ13C of Pinus. yunnanensis and Pinus. armandii at different stand ages. Sci. Rep. 2023, 13, 7938. [Google Scholar] [CrossRef] [PubMed]
  57. Jongrungklang, N.; Toomsan, B.; Vorasoot, N.; Jogloy, S.; Patanothai, A. Identification of Peanut Genotypes with High Water Use Efficiency under Drought Stress Conditions from Peanut Germplasm of Diverse Origins. Asian J. Plant Sci. 2008, 7, 628–638. [Google Scholar] [CrossRef]
Figure 1. Geographic location of the study area and distribution of the sampling plots. “Forests 15 01272 i001” indicates P. orientalis stands; “Forests 15 01272 i002” indicates P. tabulaeformis stands; “Forests 15 01272 i003” indicates Q. variabilis stands; “Forests 15 01272 i004” indicates R. pseudoacacia stands.
Figure 1. Geographic location of the study area and distribution of the sampling plots. “Forests 15 01272 i001” indicates P. orientalis stands; “Forests 15 01272 i002” indicates P. tabulaeformis stands; “Forests 15 01272 i003” indicates Q. variabilis stands; “Forests 15 01272 i004” indicates R. pseudoacacia stands.
Forests 15 01272 g001
Figure 2. Average tree-ring δ13C of four tree species in each sample plot.
Figure 2. Average tree-ring δ13C of four tree species in each sample plot.
Forests 15 01272 g002
Figure 3. Annual variations in tree-ring δ13C of four tree species between 1960 and 2020. The distributions are given in the box plots.
Figure 3. Annual variations in tree-ring δ13C of four tree species between 1960 and 2020. The distributions are given in the box plots.
Forests 15 01272 g003
Figure 4. Averaged tree-ring δ13C of each plot. The green bar indicates the average tree-ring δ13C of all the 16 types of this tree species. Independent-sample t-test was used to test the differences between the average tree-ring δ13C of all the 16 types and each type. * p < 0.05, ** p < 0.01.
Figure 4. Averaged tree-ring δ13C of each plot. The green bar indicates the average tree-ring δ13C of all the 16 types of this tree species. Independent-sample t-test was used to test the differences between the average tree-ring δ13C of all the 16 types and each type. * p < 0.05, ** p < 0.01.
Forests 15 01272 g004
Figure 5. Random forest analysis indicating the effects of density, aspect, slope, and age on tree-ring δ13C. %IncMSE, percentage of increase in mean square error; IncNodePurity, increase in node purity. * p < 0.05.
Figure 5. Random forest analysis indicating the effects of density, aspect, slope, and age on tree-ring δ13C. %IncMSE, percentage of increase in mean square error; IncNodePurity, increase in node purity. * p < 0.05.
Forests 15 01272 g005
Figure 6. Temporal series of tree-ring Δ13C, Ci and iWUE of four tree species. The solid line represents the observed values, while the dotted line indicates the linear trend.
Figure 6. Temporal series of tree-ring Δ13C, Ci and iWUE of four tree species. The solid line represents the observed values, while the dotted line indicates the linear trend.
Forests 15 01272 g006
Figure 7. Relationships between tree-ring δ13C of four tree species and the variables P, T, RH, SPEI, and VPD. The blocks are color-coded to represent different Pearson’s correlation coefficients.
Figure 7. Relationships between tree-ring δ13C of four tree species and the variables P, T, RH, SPEI, and VPD. The blocks are color-coded to represent different Pearson’s correlation coefficients.
Forests 15 01272 g007
Figure 8. Relationships between iWUE of four tree species and the variables P, T, RH, SPEI, and VPD. The blocks are color-coded to indicate different Pearson’s correlation coefficients.
Figure 8. Relationships between iWUE of four tree species and the variables P, T, RH, SPEI, and VPD. The blocks are color-coded to indicate different Pearson’s correlation coefficients.
Forests 15 01272 g008
Table 1. Basic information on different plot types in the mountainous area of Beijing. Based on stand density, the plots were classified into low-density plots (density < 1000 plants/ha) and high-density plots (density > 1000 plants/ha). Based on stand age, coniferous forest plots were divided into young forest plots (tree age < 50 a) and middle-aged forest plots (tree age > 50 a). Broad-leaved forest plots were divided into middle-aged forest plots (30 a < tree age < 50 a) and near-mature forest plots (tree age > 50 a). Based on slope, the plots were classified into gentle slope plots (slope < 15°) and steep slope plots (slope > 15°). Based on aspect, the plots were classified into sunny slope plots and shady slope plots.
Table 1. Basic information on different plot types in the mountainous area of Beijing. Based on stand density, the plots were classified into low-density plots (density < 1000 plants/ha) and high-density plots (density > 1000 plants/ha). Based on stand age, coniferous forest plots were divided into young forest plots (tree age < 50 a) and middle-aged forest plots (tree age > 50 a). Broad-leaved forest plots were divided into middle-aged forest plots (30 a < tree age < 50 a) and near-mature forest plots (tree age > 50 a). Based on slope, the plots were classified into gentle slope plots (slope < 15°) and steep slope plots (slope > 15°). Based on aspect, the plots were classified into sunny slope plots and shady slope plots.
DensityAgeSlopeAspect
HighLowMatureMiddleSteepGentleSunnyShadow
S1
S2
S3
S4
S5
S6
S7
S8
S9
S10
S11
S12
S13
S14
S15
S16
Table 2. The numerical orders of δ13C among different plots for four tree species.
Table 2. The numerical orders of δ13C among different plots for four tree species.
Tree SpeciesNumerical Order of Tree-Ring δ13C among Plots
P. orientalisP1 > P3 > P6 > P5 > P13 > P2 > P9 > P7 > P10 > P11 > P15 > P4 > P8 > P14 > P16 > P12
P. tabulaeformisP1 > P3 > P6 > P2 > P5 > P9 > P13 > P7 > P11 > P8 > P4 > P10 > P12 > P14 > P15 > P16
Q. variabilisP1 > P3 > P6 > P5 > P2 > P9 > P7 > P4 > P15 > P10 > P11 > P13 > P8 > P14 > P12 > P16
R. pseudoacaciaP1 > P3 > P6 > P5 > P2 > P9 > P7 > P4 > P15 > P10 > P11 > P13 > P8 > P14 > P12 > P16
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lu, W.; Wu, B.; Yu, X.; Jia, G.; Gao, Y.; Wang, L.; Lu, A. Tree-Ring δ13C and Intrinsic Water-Use Efficiency Reveal Physiological Responses to Climate Change in Semi-Arid Areas of North China. Forests 2024, 15, 1272. https://doi.org/10.3390/f15071272

AMA Style

Lu W, Wu B, Yu X, Jia G, Gao Y, Wang L, Lu A. Tree-Ring δ13C and Intrinsic Water-Use Efficiency Reveal Physiological Responses to Climate Change in Semi-Arid Areas of North China. Forests. 2024; 15(7):1272. https://doi.org/10.3390/f15071272

Chicago/Turabian Style

Lu, Weiwei, Bo Wu, Xinxiao Yu, Guodong Jia, Ying Gao, Lili Wang, and Anran Lu. 2024. "Tree-Ring δ13C and Intrinsic Water-Use Efficiency Reveal Physiological Responses to Climate Change in Semi-Arid Areas of North China" Forests 15, no. 7: 1272. https://doi.org/10.3390/f15071272

APA Style

Lu, W., Wu, B., Yu, X., Jia, G., Gao, Y., Wang, L., & Lu, A. (2024). Tree-Ring δ13C and Intrinsic Water-Use Efficiency Reveal Physiological Responses to Climate Change in Semi-Arid Areas of North China. Forests, 15(7), 1272. https://doi.org/10.3390/f15071272

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop