- 1Key Laboratory of Tree Breeding and Cultivation, National Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, Beijing, China
- 2Collaborative Innovation Center of Sustainable Forestry in Southern China, Nanjing Forest University, Nanjing, China
- 3Henan Xiaolangdi Earth Critical Zone National Research Station on the Middle Yellow River, Jiyuan, China
Drought-induced forest canopy die-back and tree mortality have been commonly recorded in the lithoid mountainous regions of northern China. However, the capacity of trees to regulate their carbon and water balance in response to drought remains inadequately understood. We measured tree growth, intrinsic water use efficiency (iWUE), vulnerability, and canopy health during drought events using dendrochronology, C isotope measurements, and a tree canopy health survey in a mixed plantation of Quercus variabilis and Robinia pseudoacacia. Resistance (Rt), recovery (Rc), resilience (Rs), and increased amplitude in iWUE compared to the indices 3 years before drought (iWUEr) were calculated for each species across the dominant tree (D), co-dominant tree (CD), and suppressed tree (S). Our results revealed that D and CD showed lower Rt, higher Rc, and higher iWUEr than S. After exposure to multiple sequential drought events, Q. variabilis showed an increasing trend in Rt, and R. pseudoacacia showed a decreasing trend in Rc. R. pseudoacacia exhibited a more conservative strategy towards drought, resulting in a negative SRt-iWUEr (slope of the linear model fitted to capture the trend between Rt and iWUEr) during drought events than Q. variabilis. For individual trees, lower Rc or positive SRt-iWUEr Q. variabilis and negative SRt-iWUEr R. pseudoacacia were more susceptible to canopy die-back. In conclusion, our study offers a new perspective for improved management practices in the design of silvicultural actions for forestry plantations in lithoid mountainous areas with increasing drought risk.
Introduction
Since the 1970s, China has implemented an ecosystem restoration program to repair vulnerable ecosystems and restore trees in previously degraded ecosystems (Cao et al., 2011). However, weak growth rates and tree mortality have arisen in large-scale plantations, even though many tree species have been chosen for drought tolerance (Wu et al., 2004; Sun et al., 2018b; Zhao et al., 2018). For these plantations, despite major progress in understanding how water deficits affect functioning (Wang et al., 2017; Sun et al., 2018a), the relationship between functional indices and species survival under extreme drought remains poorly understood and limits our ability to adequately predict tree mortality risk accompanied by continuing warming and aridification in the future.
Recent syntheses of drought physiology have highlighted the need to consider the capacity of trees to regulate their carbon and water balances under drought conditions (Mitchell et al., 2013; Gentilesca et al., 2021). These frameworks suggest that carbon starvation (Sevanto et al., 2014) or hydraulic (Grote et al., 2016) failure may be responsible for the drought-induced mortality in trees. Also, it has been observed that drought stress changes with selection of tree species or tree species mixtures (Pretzsch et al., 2013; Merlin et al., 2015). Different species growing in the same soil water status may perform differently in a mixed forest, and drought will likely disturb the current equilibrium between the co-existing species (Fekedulegn et al., 2003; Tognetti et al., 2007; Jucker et al., 2014). For example, in mixed plantations on the Loess Plateau, top shoot die-back and mortality are frequent in Populus hopeinsis and Robinia pseudoacacia. However, other species such as Ulmus pumila rarely have similar difficulties (Zhang and Shang, 2002). The reasons for these species-specific differences in response to drought are still not thoroughly understood. In addition to these drivers, the presence of various individual tree social statuses within a stand adds another layer of complexity to the tree’s response to drought (Merlin et al., 2015). This social position suggests a disparity in the accessibility of resources such as water, nutrients, and light, which results in differences in the reactions to severe droughts (Trouvé et al., 2017). Consequently, some studies have comprehensively discussed the effect of water stress on the growth of deciduous trees of various social statuses. However, these findings are inconsistent (Liu and Muller, 1993; He et al., 2005). Therefore, it is critical to incorporate both species and social status into the stand as explanatory factors when evaluating the carbon and water balance responses of mixed plantations to drought.
In particular, tree growth resilience indices and intrinsic water-use efficiency (iWUE) are important indicators for the response analyses of carbon and water balance to drought (Forner et al., 2018; Sun et al., 2018a; Nolan et al., 2021; Zhang et al., 2022). Forest stability reactions during drought events are typically categorized as having either immediate or delayed effects, representing their ability to maintain functions during and after drought events (Trouvé et al., 2017; Zas et al., 2020). iWUE is defined as the net photosynthetic rate to stomatal conductance, which evaluates the water loss through the stomata. The trade-off between carbon capture and water loss is not a serious conflict in favorable climate conditions, resulting in a slight increase in tree growth when iWUE increases (Guerrieri et al., 2019; Lu et al., 2019). Under extreme drought events, the increase in iWUE is mainly induced by closing stomata to save water and maintain their leaf water potential within a safe range at the expense of reduced carbon uptake (Martin-Benito et al., 2017). In this case, the growth resistance of the trees will decrease or remain constant, which has been reported in the recent literature on dead trees under controlled environmental conditions (Sun et al., 2018b; Zadworny et al., 2019). Previous research has provided insights into tree growth and iWUE responses to drought in North China, although the results are inconsistent (Wu et al., 2015; Sun et al., 2018b; Song et al., 2019). However, few empirical studies have evaluated the effectiveness of these indicators for predicting the risk of drought-induced mortality.
In this study, we used dendrochronology and carbon isotopes to understand the growth resilience indices and iWUE for Q. variabilis and R. pseudoacacia across three crown classes. We related growth resilience indices, iWUE, with the current canopy health score, a proxy for vulnerability to canopy die-back, to assess the relationships between these indicators and drought-induced mortality. The following two questions were the focus of our investigation: (1) How did growth resilience indices and iWUE vary between species, crown classes, and drought events? (2) Which post-drought resilience and iWUE indices influence current canopy die-back?
Materials and methods
Study area
This study was conducted in a young mixed forest at a height of approximately 410 m above sea level in the lithoid hilly region of North China (35°02’ N, 112°28’ E) (Figure 1). Warm-temperate continental monsoon weather prevailed at the site over the course of 50 years, with a mean annual temperature of 14.5°C and mean annual precipitation of 609.3 mm. In the study region, there has been a significant increase in the annual mean temperature over the last 30 years (Figure 2A). The geological substrate is composed of limestone. The soil at the site was 40–50 cm deep and classified as brown loam. The sites were covered by 46-year-old mixed deciduous Q. variabilis and R. pseudoacacia trees.
Figure 2 (A) Average annual temperature (T), the dotted line indicates the trends considering T for the 1990–2018 period; (B) summer (June to August) precipitation (P) and SPEI calculated from October of the previous year to September of the current year (SPEI12sep) – dotted line for the 1960–2018 period. The dark line indicates the SPEI threshold of -1.5 (indicative of severe drought).
Climatic data
Monthly mean temperature and total precipitation data were obtained from the Jiyuan County Meteorological Station (35°5’ N, 112°38’ E, 10 km from the site; Figure 1). We obtained the standardized precipitation evapotranspiration index (SPEI) calculated from October of the previous year to September of the current year (SPEI12sep) to evaluate drought severity with the SPEI program (Vicente-Serrano et al., 2010). This methodological approach avoids climatic parameters in the months following stem growth (Conte et al., 2018). Positive SPEI values imply wet conditions, whereas negative values indicate dry conditions. We selected three drought years (years with SPEI12sep< −1.5: 2002, 2009, and 2013) for the period 1960–2018 (Figure 2B).
Field sampling and dendrochronological methods
A forest survey was conducted in August 2019 by establishing 23 inventory plots (30 × 30 m). Geographic coordinates and characteristics of these plots are provided in Supplementary Table S1. In each plot, trees with a diameter at breast height (DBH) ≥ 4.0 cm were permanently marked, measured (DBH, height X–Y coordinates), and recorded (species and crown class). The crown classes were defined as dominant (D), co-dominant (CD), and suppressed (S). In the present study, D trees had crowns extending above the main canopy, CD trees had crowns making up the main canopy, and S trees had crowns extending into the lower portion of the main canopy (Seidel et al., 2011).
Fifty to sixty trees were selected from each plot. In total, 1265 trees were sampled (799 Q. variabilis; 466 R. pseudoacacia). Two complete cores (90° from one another) were collected per sample tree perpendicular to a height of 1.3 m with standard 5.15-mm increment borers (Häglof, Långsele, Sweden). The cores were air-dried and sanded with papers of progressively finer grains until tree rings were clearly visible. Tree rings were dated under a stereomicroscope (Leica, Germany), and tree-ring width was measured using a LINTAB measuring device (Frank Rinn, Heidelberg, Germany, resolution: 0.01 mm). Cross-dating was further validated using COFECHA v11 (Holmes, 1983).
Several statistics commonly used in dendrochronology were calculated to assess the quality of the tree-ring width series using the dplR package (Bunn, 2008): first-order autocorrelation (AC) of raw width data (Fritts, 1976), mean sensitivity (MS) of indexed growth values, mean correlation (Rbt) between trees (Briffa, 1990), and expressed population signal (EPS) (Wigley et al., 1984), which measures the statistical quality of the mean site chronology as compared with a perfect, infinitely replicated chronology (Table 1).
Table 1 Chronology statistics in dominant trees (D), co-dominant trees (CD) and suppressed trees (S) for sampled Q. variabilis and R. pseudoacacia.
The past diameters were extrapolated using the tree-ring widths of the cores. A proportional approach was used to reduce the impact of eccentricity on the growth increment (Bakker, 2005). The tendency caused by the geometrical limitation of adding a volume of wood to a stem with increasing radius was corrected by translating tree-ring widths into basal area increments (BAIs) as follows:
where R is the tree radius and n is the year of tree-ring formation.
Growth resilience indices
To evaluate growth resistance and resilience to each drought event, we calculated three indices for each species and crown class using BAI data: resistance (Rt), recovery (Rc), and resilience (Rs) (Lloret et al., 2011).
where BAID is the mean BAI during a drought event and BAIpre and BAIpost are the mean BAIs in the 3 years before or after the drought event.
Water use efficiency analysis
After cross-dating, several similar-age tree cores of each crown class were manually separated into slivers using a dissecting scope and razor blade. A total of 135 trees were selected (68 Q. variabilis and 67 R. pseudoacacia). To estimate the effects of the iWUE response to drought, we analyzed iWUE in drought events (iWUED) and the average iWUE in the 3 years before the drought event (iWUEB). The wood of each drought event and 3 years before each drought event was homogenized into a sample. All samples were ground into powder using a ball mill instrument. Whole wood samples were subjected to isotopic analysis as recent research has demonstrated that whole wood can be used instead of a-cellulose in short-term eco-physiological and dendrochronological studies to measure stable C isotopes (Au and Tardif, 2009). Approximately 3–4 mg of each sample (dried for 24 h in an oven at 60°C) in a tin cup was completely burnt to gas using a TOC Element Analyzer (Elementar Analysensysteme GmbH. Hanau, Germany). It was then injected into a carbon dioxide isotope analyzer (CCIA-38-EP, Los Gatos Research, USA) to determine δ13C. Three subsamples were randomly selected from each mixed sample, each of which was taken twice. The bi-weight method was then used to average the six values from each mixed sample into a single value. The δ13C values were expressed using the Vienna Pee Dee Belemnite scale.
iWUE (μmol mol−1) was calculated according to standard methodologies and formulae was calculated using (Farquhar and Richards, 1984):
where A is the rate of CO2 assimilation by the leaves, gs is the rate of leaf stomatal conductance, Ci is the leaf intercellular CO2 concentration and Ca is. atmospheric CO2 concentration (Ca). To determine Ci/Ca, we used the equation 6 (Farquhar et al., 1982):
where δ13Ctree represents δ13 measured from tree rings; δ13Catm is the value for the atmospheric value obtained from the Earth System Research Laboratory of the U.S. National Oceanic and Atmospheric Administration. a is the discrimination due to diffusion of 13CO2 through stomata (a = 4.4‰) and b is fractionation discrimination by Rubisco against 13CO2 (b = 27‰).
To evaluate the degree of change in iWUE during drought events, we calculated the ratio of iWUED and iWUEB (iWUEr).
Canopy health
For each sample tree, health condition was evaluated by the canopy health score described previously (Stone et al., 2008; Nolan et al., 2021) based on visual estimates in August 2019, which summed five crown attribute scores (scale of 1 to 5): crown size, crown density, dead branches, crown epicormic growth, and leaf discoloration. Additional details are presented in Supplementary Table S2. Scores ranged from 0 for trees with no leaves left to 25 for healthy trees.
Data analysis
The differences in Rt, Rc, Rs, iWUE, and canopy health score among crown classes of Q. variabilis and R. pseudoacacia were analyzed using the least significant difference (LSD) test. Prior to statistical analysis, we checked all indices for normality and homoscedasticity, and statistical significance was set at p< 0.05.
To test the relationship between radial growth and iWUE during the three drought events, Rt was linearly regressed against iWUEr for each crown class. In addition, a linear model was fitted to capture the trend between Rt and iWUEr for each sample tree, with the slope of this model (SRt-iWUEr) indicating the magnitude of the growth reduction with an increase in iWUE during the three drought events. We evaluated the tree canopy health score as the response variable in linear regressions against the growth resilience indices, iWUEr and SRt-iWUEr.
R 4.1.3 was used to conduct all statistical analyses (Team R C, 2013).
Results
The EPS values were all greater than 0.85 in every instance to suggesting the qualities of chronologies were sufficient good (Table 1). The Rbt, MS, and EPS values were higher for D and CD than for S.
Radial growth analysis
Q. variabilis and R. pseudoacacia for the three canopy classes had a similar pattern in their BAIs: a steady increase during the first 20 years followed by a substantial and significant decrease during the three drought events (Figure 3). A significant negative trend of BAI was observed in suppressed R. pseudoacacia (−1.5% year-1, p<0.01) after the first 20 years, with the trends also being slightly negative but not significant for other types of trees. Q. variabilis and R. pseudoacacia had lower average BAI values in S than in D and CD.
Figure 3 Mean curves of basal area increment (BAI) chronologies of Q. variabilis and R. pseudoacacia for three canopy classes. The black dotted line corresponds to 2002, 2009, and 2013 drought.
Resistance, recovery and resilience to drought
For Q. variabilis and R. pseudoacacia, in terms of drought Rt and Rc, significant differences in responses were observed for each drought event between the three canopy classes. The order of the Rt values was S > CD > D, and the order of the Rc values was D > CD > S (Figure 4). With the increase in drought events experienced by trees, the Rt of Q. variabilis in each canopy class showed an increasing trend, and the Rc of R. pseudoacacia in each canopy class showed a decreasing trend. During the 2002 drought event, D was less resilient than S. On the other hand, S was less resilient than D during the 2009 and 2013 drought events.
Figure 4 Mean resistance (Rt), recovery (Rc), resilience (Rs) calculated for Q. variabilis and R. pseudoacacia of three canopy classes during drought events. Different letters indicate significant difference (p< 0.05) in drought responses between three canopy classes.
iWUE responses to drought
The iWUE of both species significantly decreased in response to drought in relation to the iWUE of the three years preceding drought events (Table 2). For the average iWUEr of the three drought events between the three canopy classes, the order of Q. variabilis was D (1.21) > CD (1.15) > S (1.11) and that of R. pseudoacacia was CD (1.22) > D (1.15) > S (1.10).
Relationships between Rt and iWUEr
For Q. variabilis, the relationships between Rt and iWUEr of D showed significant positive trends (p< 0.05) during the three drought events (Figure 5). In contrast, no significant trends of Rt-iWUEr were observed in the CD and S of Q. variabilis, except for S in the 2002 drought event (p = 0.004). For R. pseudoacacia, the Rt-iWUEr of the three canopy classes showed significant negative trends during the three drought events (p< 0.05). The negative correlation between Rt and iWUEr eased as the slopes of the regression lines increased with an increase in drought events experienced by trees.
Figure 5 Relationship between Rt and iWUEr of Q. variabilis and R. pseudoacacia for three canopy classes during drought events.
Growth resilience indices, iWUEr and canopy health
The significant trends were observed between canopy health scores and Rc and SRt-iWUEr of Q. variabilis, on the other hand, health scores increased with SRt-iWUEr of R. pseudoacacia (Figure 6), indicating that lower Rc and positive SRt-iWUEr Q. variabilis, negative SRt-iWUEr R. pseudoacacia were the more vulnerable to canopy die-back during drought. The canopy health score as a function of crown class was analyzed using ANOVA, and the results showed that S trees for these two species were significantly more susceptible to canopy die-back than D and CD trees.
Figure 6 Relationship between canopy health scores and growth resilience indices, iWUEr, and SRt-iWUEr of Q. variabilis and R. pseudoacacia.
Discussion
Despite receiving far less scholarly attention than natural forests do, plantations play an essential role in rehabilitating endangered ecosystems (Brockerhoff et al., 2008; Rodriguez-Vallejo et al., 2021). By combining resistance indices, iWUE, and canopy health score data, we have provided valuable field data to help predict tree mortality risk in these forests.
Effects of species and crown classes on BAI and iWUE response to drought
Similar to previous observations, the radial growth of both species was significantly reduced under drought-induced water stress (Martín-Benito et al., 2008; Eilmann and Rigling, 2012; Mazza et al., 2021). A decline in the carbon stock available for functions during drought events is associated with reducing photosynthetic products and changing the method of carbon allocation and utilization (Flexas et al., 2006; Körner, 2015). For example, drought promotes more carbohydrates and starches in tree leaves to be decomposed to regulate osmotic pressure (Chaves et al., 2003)e, and some studies have also shown that under drought, trees may accumulate a larger proportion of non-structural carbohydrates through photosynthesis, which is used to maintain respiration in the dormant season and the growth of earlywood in the following year, which leads to the hysteresis of growth (Lacointe, 2000; Richburg, 2005).
D and CD trees showed higher growth reductions than S trees for both species during drought, whereas the recovery of D and CD was enhanced afterward, consistent with some of the existing studies (Martín-Benito et al., 2008; Merlin et al., 2015; Trouvé et al., 2017). One explanation could be that D and CD trees are more exposed to drought as they intercept more resources in forests dominated by asymmetrical competition (Linares et al., 2010). In addition, despite the general drought conditions, S trees may benefit from better growth conditions (slightly cooler, less radiation, or better ventilation) than D and CD trees because they have less leaf area exposed to intense radiation and high temperatures owing to partial shading (Aussenac, 2000; von Arx et al., 2012). Other features of S trees, including smaller crowns, lower root densities (Liu and Muller, 1993; Augspurger and Bartlett, 2003), and lower carbohydrate reserves (Richardson et al., 2013) may hinder their recovery to pre-drought growth levels (Galiano et al., 2011; Grote et al., 2016).
Accumulated stresses from multiple drought events led to an increase in Rt in Q. variabilis and a decrease in Rc in R. pseudoacacia. R. pseudoacacia, as an early successional species, seems to be maladjusted to frequent drought years typical of low altitudes in North China (Tanaka-Oda et al., 2010). Its poor performance in recovering growth after drought might be due to its ineffective control mechanism in dealing with water loss (Eilmann and Rigling, 2012). We anticipated that given the different levels of drought tolerance among the species, the growth releases would come from the death of the less drought-resistant R. pseudoacacia, benefiting Q. variabilis radial growth during the upcoming drought event (Rubio-Cuadrado et al., 2018).
Rs measures the capacity of trees to recover their pre-drought growth rate after drought, with Rs< 1 indicating a persistently detrimental effect caused by drought on subsequent growth. With an increase in drought events experienced by trees, the Rs of D showed increasing trends, and the Rs of CD and S showed decreasing trends. This may indicate that the sensitivity of the drought response changed after withstanding numerous drought occurrences. CD and S will be more vulnerable to drought than D, and are more likely to experience carbon famine (Trifil et al., 2017), increasing the likelihood that they will die.
R. pseudoacacia, a non-native tree species, had slightly higher iWUE than Q. variabilis (a native tree species). Similar results have been reported for other native and nonnative tree species. For example, Song et al. (2019) reported a higher increase in iWUE for non-native trees than for native tree species in forests in Northeast China. When drought occurs, a higher increase in the iWUE for the non-native species compared with the native species seems to be attributable to higher climatic sensitivity or corresponding stronger water stress. In contrast, R. pseudoacacia has more conservative water use (Fan, 2010; Nadal-Sala et al., 2017), which partly reflects the differences in xylem structure. R. pseudoacacia produces more vessels, but its diameter appears to be smaller than that of oaks (Nola et al., 2020; Özden Keleş and Savaci, 2021). Greater resistance to water transportation and cavitation under drought conditions may be achieved with vessels of smaller diameters (Sperry et al., 2006), resulting in lower hydraulic conductivity and higher iWUE. In addition, under drought conditions, R. pseudoacacia closes its stomata at an early stage compared to Q. variabilis, resulting in a rapid decline in stomatal conductance and photosynthesis (Fan, 2010; Jin et al., 2011; Colangelo et al., 2018), at the same time, more apparent morphological changes in adaptation to water limitation were observed in R. pseudoacacia, such as leaf shedding (Wang et al., 2017), which could partially explain the higher iWUE to drought observed in R. pseudoacacia.
S trees were below the main canopy, resulting in less coupling with the atmosphere, which created a higher ratio between boundary layer conductance and canopy stomatal conductance (Wullschleger et al., 2000), as well as assimilated respired CO2 (lower δ13C) within the canopy, and drought induced stomatal closure. Therefore, reduced mesophyll conductance (gm) could increase iWUE. Height alone is supposed to make water transport to the leaves more difficult, resulting in a decrease in gm with height (Cano et al., 2013); the lower the height of S trees compared to D and CD trees, the lower the gm. These features may lead to lower iWUE for S trees.
Drought modifies Rt-iWUEr
In the lithoid hilly areas with shallow soil, the negative trend between Rt and iWUEr could be the result of a more substantial reduction in soil moisture, which causes plants to close their stomata to save water and maintain their leaf water potential within a safe range at the expense of reduced carbon uptake (McDowell, 2011). Unlike R. pseudoacacia, Q. variabilis showed no significant negative trend in Rt-iWUEr, which agrees with previous results on contrastive research of native and non-native tree species (Song et al., 2019). Non-native tree species may have greater stomatal control than native species, consequently enhancing iWUE to the detriment of growth (Gentilesca et al., 2021). In addition, the absence of Rt-iWUEr might be attributed to variable allocation to other tissues, and remobilization of carbohydrate reserves frequently causes a disconnect between tree growth and carbon assimilation (Urrutia-Jalabert et al., 2015).
Canopy health
The difference in relationships between canopy health scores, growth resilience indices, and iWUEr of R. pseudoacacia and Q. variabilis revealed partly different reasons for drought-induced mortality. For R. pseudoacacia, SRt-iWUEr was negatively linked to canopy health scores, suggesting that carbon starvation may cause canopy die-back in trees undergoing severe drought (Wang et al., 2017; Colangelo et al., 2018). The lower SRt-iWUEr indicates that the increased iWUE resulting from multiple drought events prevented carbon uptake by R. pseudoacacia. However, reduced carbon uptake has a negative effect on driving phloem transport, maintaining turgor, and refilling embolized xylem during drought, and mortality occurs as soon as one or more of these processes reaches a threshold (McDowell et al., 2011). R. pseudoacacia showed a growth reduction in drought-induced mortality at a local scale (Wang et al., 2017). According to their findings, no significant difference was observed in the growth rate between dead and healthy R. pseudoacacia before severe drought. After multiple drought events, the growth rate of dying trees fell and remained very low, even resulting in only a row of earlywood conduits without latewood being observed on tree cores during drought events (Wang et al., 2017). Current photosynthetic products are not readily available because of the influence of the main canopy on the photosynthetic capacity of the S trees (Martín-Benito et al., 2008). Structural growth is associated with higher demands for older carbon pools, which are primarily reserves accumulated during favorable growth years. However, replenishing older carbon pools with new sugar inputs may become more challenging after multiple drought events (Richardson et al., 2013). Therefore, for R. pseudoacacia, the lower carbon storage and SRt-iWUEr of S trees compared with those of D and CD trees may explain their higher mortality.
There was a clear link between low Rc capacity and high SRt-iWUEr to drought in Q. variabilis and the increased risk of canopy die-back. Previous studies have revealed a typical pattern for oaks experiencing multiple drought events: oaks usually die a few years after the drought because the time period needed for oak reserves to recover may be too long, resulting in a gradual decline in resilience and increased mortality rates (Galiano et al., 2012; Granda et al., 2013; Granda et al., 2014). For Q. variabilis, a higher percentage of consumption of stored nonstructural carbohydrates during drought events resulted in a positive trend in Rt-iWUEr (higher SRt-iWUEr), further limiting the ability of trees to recover following drought occurrences (Galiano et al., 2012; Rosas et al., 2013). No evidence of such depletion of individuals was found in this study, but inconsistent findings were found between net ecosystem productivity based on eddy covariance and Q. variabilis biomass growth based on biometric methods in a local mixed deciduous plantation dominated by Q. variabilis during drought events, partly because of the consumption of stored non-structural carbohydrates (Jia et al., 2022).
Conclusion
The effect of drought on the growth and functioning ofin plantations has been anone of the important concerns for researchers and forest managers to proactively address climate-related global challenges. Our study confirms the differing responses of three crown classes, R. pseudoacacia, and Q. variabilis, to several specific drought events.: During the drought, D and CD showed a higher reduction in growth and higher increase amplitude in iWUE than S, did but also recovered faster afterward. A more conservative water-saving strategy was reflected in the negative correlation between Rt and iWUEr inof R. pseudoacacia than in Q. variabilis. Lower Rc, and positive SRt-iWUEr Q. variabilis, and negative SRt-iWUEr R. pseudoacacia were more susceptible to canopy die-back during drought due to stress from multiple drought events. This study provides helpful information on species selection and management measures for plantations in lithoid mountainous areas with increasing drought risk.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
HJ and CG contributed to conceptualization, chronology data analysis and wrote the manuscript. CH contributed to stable C isotopes, canopy health score data analysis. JZ, CY, PM contributed to design of the study and funded the study.
Funding
This work was supported by a National Nonprofit Institute Research Grant of the Chinese Academy of Forestry (CAFYBB2018ZA001) and the National Key Research and Development Project (2020YFA0608101). The China Meteorological Data Service Center provided meteorological data.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.1018405/full#supplementary-material
Abbreviations
AC, Autocorrelation; BAI, Basal area increment; CD, Co-dominant; D, Dominant; DBH, Diameter at breast height; EPS, Expressed population signal; iWUE, intrinsic water-use efficiency; MS, Mean sensitivity; S, Suppressed; SPEI, Standardized precipitation evapotranspiration index.
References
Au, R., Tardif, J. C. (2009). Chemical pretreatment of thuja occidentalis tree rings: implications for dendroisotopic studies. Can. J. For. Res. 39 (9), 1777–1784. doi: 10.1139/X09-091
Augspurger, C. K., Bartlett, E. A. (2003). Differences in leaf phenology between juvenile and adult trees in a temperate deciduous forest. Tree Physiol. 23 (8), 517–525. doi: 10.1093/treephys/23.8.517
Aussenac, G. (2000). Interactions between forest stands and microclimate: ecophysiological aspects and consequences for silviculture. Ann. For. Sci. 57 (3), 287–301. doi: 10.1051/forest:2000119
Bakker, J. D. (2005). A new, proportional method for reconstructing historical tree diameters. Can. J. For. Res. 35 (10), 2515–2520. doi: 10.1139/x05-136
Briffa, K. (1990). “Basic chronology statistics and assessment,” In Methods dendochronology. eds. Cook, E. R., Kairiukstis, L. (Dordrecht: Kluwer Academic Publishers), 137–153.
Brockerhoff, E. G., Jactel, H., Parrotta, J. A., Quine, C. P., Sayer, J. (2008). Plantation forests and biodiversity: oxymoron or opportunity? Biodiversity Conserv. 17 (5), 925–951. doi: 10.1006/meth.2000.1142
Bunn, A. G. (2008). A dendrochronology program library in r (dplR). Dendrochronologia 26 (2), 115–124. doi: 10.1016/j.dendro.2008.01.002
Cano, F. J., SÁNCHEZ-GÓMEZ, D., RODRÍGUEZ-CALCERRADA, J., Warren, C. R., Gil, L., Aranda, I. (2013). Effects of drought on mesophyll conductance and photosynthetic limitations at different tree canopy layers. Plant Cell Environ. 36 (11), 1961–1980. doi: 10.1111/pce.12103
Cao, S., Chen, L., Shankman, D., Wang, C., Wang, X., Zhang, H. (2011). Excessive reliance on afforestation in china's arid and semi-arid regions: lessons in ecological restoration. Earth-Sci. Rev. 104 (4), 240–245. doi: 10.1016/j.earscirev.2010.11.002
Chaves, M. M., Maroco, J. P., Pereira, J. S. (2003). Understanding plant responses to drought–from genes to the whole plant. Funct. Plant Biol. 30 (3), 239–264. doi: 10.1071/FP02076
Colangelo, M., Camarero, J. J., Ripullone, F., Gazol, A., Sánchez-Salguero, R., Oliva, J., et al. (2018). Drought decreases growth and increases mortality of coexisting native and introduced tree species in a temperate floodplain forest. Forests 9 (4), 205. doi: 10.3390/f9040205
Conte, E., Lombardi, F., Battipaglia, G., Palombo, C., Altieri, S., La Porta, N., et al. (2018). Growth dynamics, climate sensitivity and water use efficiency in pure vs. mixed pine and beech stands in trentino (Italy). For. Ecol. Manage. 409, 707–718. doi: 10.1016/j.foreco.2017.12.011
Eilmann, B., Rigling, A. (2012). Tree-growth analyses to estimate tree species' drought tolerance. Tree Physiol. 32 (2), 178–187. doi: 10.1093/treephys/tps004
Fan, W. (2010). Water utilization characteristics of 7 typical plant species in hilly region of taihang mountains. Sci. Soil Water Conserv. 8 (5), 61–66. doi: 10.3969/j.issn.1672-3007.2010.05.011
Farquhar, G. D., O'Leary, M. H., Berry, J. A. (1982). On the relationship between carbon isotope discrimination and the intercellular carbon dioxide concentration in leaves. Funct. Plant Biol. 9 (2), 121–137. doi: 10.1071/PP9820121
Farquhar, G., Richards, R. (1984). Isotopic composition of plant carbon correlates with water-use efficiency of wheat genotypes. Funct. Plant Biol. 11 (6), 539–552. doi: 10.1071/pp9840539
Fekedulegn, D., Hicks, R. R., Jr., Colbert, J. (2003). Influence of topographic aspect, precipitation and drought on radial growth of four major tree species in an Appalachian watershed. For. Ecol. Manage. 177 (1-3), 409–425. doi: 10.1016/S0378-1127(02)00446-2
Flexas, J., Bota, J., Galmés, J., Medrano, H., Ribas‐Carbó, M. (2006). Keeping a positive carbon balance under adverse conditions: responses of photosynthesis and respiration to water stress. Physiol Plantarum 127 (3), 343–52. doi: 10.1111/j.1399-3054.2006.00621.x
Forner, A., Valladares, F., Bonal, D., Granier, A., Grossiord, C., Aranda, I. (2018). Extreme droughts affecting Mediterranean tree species’ growth and water-use efficiency: the importance of timing. Tree Physiol. 38 (8), 1127–1137. doi: 10.1093/treephys/tpy022
Galiano, L., Martínez-Vilalta, J., Lloret, F. (2011). Carbon reserves and canopy defoliation determine the recovery of scots pine 4 yr after a drought episode. New Phytol. 190 (3), 750–759. doi: 10.1111/j.1469-8137.2010.03628.x
Galiano, L., Martínez-Vilalta, J., Sabaté, S., Lloret, F. (2012). Determinants of drought effects on crown condition and their relationship with depletion of carbon reserves in a Mediterranean holm oak forest. Tree Physiol. 32 (4), 478–489. doi: 10.1093/treephys/tps025
Gentilesca, T., Battipaglia, G., Borghetti, M., Colangelo, M., Altieri, S., Ferrara, A., et al. (2021). Evaluating growth and intrinsic water-use efficiency in hardwood and conifer mixed plantations. Trees 35 (4), 1329–1340. doi: 10.1007/s00468-021-02120-z
Granda, E., Camarero, J. J., Gimeno, T. E., Martinez-Fernandez, J., Valladares, F. (2013). Intensity and timing of warming and drought differentially affect growth patterns of co-occurring Mediterranean tree species. Eur. J. For. Res. 132 (3), 469–480. doi: 10.1007/s10342-013-0687-0
Granda, E., Scoffoni, C., Rubio-Casal, A. E., Sack, L., Valladares, F. (2014). Leaf and stem physiological responses to summer and winter extremes of woody species across temperate ecosystems. Oikos 123 (11), 1281–1290. doi: 10.1111/oik.01526
Grote, R., Gessler, A., Hommel, R., Poschenrieder, W., Priesack, E. (2016). Importance of tree height and social position for drought-related stress on tree growth and mortality. Trees 30 (5), 1467–1482. doi: 10.1007/s00468-016-1446-x
Guerrieri, R., Belmecheri, S., Ollinger, S. V., Asbjornsen, H., Jennings, K., Xiao, J., et al. (2019). Disentangling the role of photosynthesis and stomatal conductance on rising forest water-use efficiency. Proc. Natl. Acad. Sci. 116 (34), 16909–16914. doi: 10.1073/pnas.1905912116
He, J.-S., Zhang, Q.-B., Bazzaz, F. (2005). Differential drought responses between saplings and adult trees in four co-occurring species of new England. Trees 19 (4), 442–450. doi: 10.1007/s00468-004-0403-2
Holmes, R. L. (1983). Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bull. 44, 69–75. doi: 10.1006/biol.1999.0214
Jia, H., Zhou, Y., Zhang, J., Sun, S., Meng, P. (2022). Early–latewood carbon isotope enhances the understanding of the relationship between tree biomass growth and forest-level carbon fluxes. Agric. For. Meteorol. 315, 108818. doi: 10.1016/j.agrformet.2022.108818
Jin, T., Liu, G., Fu, B., Ding, X., Yang, L. (2011). Assessing adaptability of planted trees using leaf traits: A case study with Robinia pseudoacacia L. in the Loess Plateau, China. Chin. Geogr. Sci. 21 (3), 290. doi: 10.1007/s11769-011-0470-4
Jucker, T., Bouriaud, O., Avacaritei, D., Dănilă, I., Duduman, G., Valladares, F., et al. (2014). Competition for light and water play contrasting roles in driving diversity–productivity relationships in Iberian forests. J. Ecol. 102 (5), 1202–1213. doi: 10.1111/1365-2745.12276
Körner, C. (2015). Paradigm shift in plant growth control. Curr. Opin. Plant Biol. 25, 107–114. doi: 10.1016/j.pbi.2015.05.003
Lacointe, A. (2000). Carbon allocation among tree organs: a review of basic processes and representation in functional-structural tree models. Ann. For. Sci. 57 (5), 521–533. doi: 10.1051/forest:2000139
Linares, J. C., Camarero, J. J., Carreira, J. A. (2010). Competition modulates the adaptation capacity of forests to climatic stress: insights from recent growth decline and death in relict stands of the Mediterranean fir abies pinsapo. J. Ecol. 98 (3), 592–603. doi: 10.1111/j.1365-2745.2010.01645.x
Liu, Y., Muller, R. N. (1993). Effect of drought and frost on radial growth of overstory and undesrstory stems in a deciduous forest. Am. Midland Nat. 129, 19–25. doi: 10.2307/2426431
Lloret, F., Keeling, E. G., Sala, A. (2011). Components of tree resilience: effects of successive low-growth episodes in old ponderosa pine forests. Oikos 120 (12), 1909–1920. doi: 10.1111/j.1600-0706.2011.19372.x
Lu, K., Chen, N., Zhang, C., Dong, X., Zhao, C. (2019). Drought enhances the role of competition in mediating the relationship between tree growth and climate in semi-arid areas of Northwest China. Forests 10 (9), 804. doi: 10.3390/f10090804
Martin-Benito, D., Anchukaitis, K. J., Evans, M. N., Del Río, M., Beeckman, H., Cañellas, I. (2017). Effects of drought on xylem anatomy and water-use efficiency of two co-occurring pine species. Forests 8 (9), 332. doi: 10.3390/f8090332
Martín-Benito, D., Cherubini, P., Del Río, M., Cañellas, I. (2008). Growth response to climate and drought in pinus nigra arn. trees of different crown classes. Trees 22 (3), 363–373. doi: 10.1007/s00468-007-0191-6
Mazza, G., Markou, L., Sarris, D. (2021). Species-specific growth dynamics and vulnerability to drought at the single tree level in a Mediterranean reforestation. Trees 35, 1697–1710. doi: 10.1007/s00468-021-02151-6
McDowell, N. G. (2011). Mechanisms linking drought, hydraulics, carbon metabolism, and vegetation mortality. Plant Physiol. 155 (3), 1051–1059. doi: 10.1104/pp.110.170704
McDowell, N. G., Beerling, D. J., Breshears, D. D., Fisher, R. A., Raffa, K. F., Stitt, M. (2011). The interdependence of mechanisms underlying climate-driven vegetation mortality. Trends Ecol. Evol. 26 (10), 523–532. doi: 10.1016/j.tree.2011.06.003
Merlin, M., Perot, T., Perret, S., Korboulewsky, N., Vallet, P. (2015). Effects of stand composition and tree size on resistance and resilience to drought in sessile oak and scots pine. For. Ecol. Manage. 339, 22–33. doi: 10.1016/j.foreco.2014.11.032
Mitchell, P. J., O'Grady, A. P., Tissue, D. T., White, D. A., Ottenschlaeger, M. L., Pinkard, E. A. (2013). Drought response strategies define the relative contributions of hydraulic dysfunction and carbohydrate depletion during tree mortality. New Phytol. 197 (3), 862–872. doi: 10.1016/j.isprsjprs.2021.11.002
Nadal-Sala, D., Sabaté, S., Sánchez-Costa, E., Poblador, S., Sabater, F., Gracia, C. (2017). Growth and water use performance of four co-occurring riparian tree species in a Mediterranean riparian forest. For. Ecol. Manage. 396, 132–142. doi: 10.1016/j.foreco.2017.04.021
Nola, P., Bracco, F., Assini, S., von Arx, G., Castagneri, D. (2020). Xylem anatomy of robinia pseudoacacia l. and quercus robur l. @ is differently affected by climate in a temperate alluvial forest. Ann. For. Sci. 77 (1), 1–16. doi: 10.1007/s13595-019-0906-z
Nolan, R. H., Gauthey, A., Losso, A., Medlyn, B. E., Smith, R., Chhajed, S. S., et al. (2021). Hydraulic failure and tree size linked with canopy die-back in eucalypt forest during extreme drought. New Phytol. 230 (4), 1354–1365. doi: 10.1111/nph.17298
Özden Keleş, S., Savaci, G. (2021). Seasonal variation of morphological traits, wood formation, and soil properties differs between Quercus robur L. and Robinia pseudoacacia L. saplings. Scandinavian J. For. Res. 5 (36), 344–53. doi: 10.1080/02827581.2021.1941237
Pretzsch, H., Bielak, K., Block, J., Bruchwald, A., Dieler, J., Ehrhart, H.-P., et al. (2013). Productivity of mixed versus pure stands of oak (Quercus petraea (M att.) l iebl. and quercus robur l.) and European beech (Fagus sylvatica l.) along an ecological gradient. Eur. J. For. Res. 132 (2), 263–280. doi: 10.1007/s10342-012-0673-y
Richardson, A. D., Carbone, M. S., Keenan, T. F., Czimczik, C. I., Hollinger, D. Y., Murakami, P., et al. (2013). Seasonal dynamics and age of stemwood nonstructural carbohydrates in temperate forest trees. New Phytologist. 197, 850–861. doi: 10.1111/nph.12042
Richburg, J. A. (2005). Timing Treatments to the Phenology of Root Carbohydrate Reserves to Control Woody Invasive Plants, A Ph.D. dissertation submitted to the University of Massachusetts Amherst.
Rodriguez-Vallejo, C., Navarro-Cerrillo, R. M., Manzanedo, R. D., Rodriguez, G. P., Gazol, A., Camarero, J. J. (2021). High resilience, but low viability, of pine plantations in the face of a shift towards a drier climate. For. Ecol. Manage. 479, 118537. doi: 10.1016/j.foreco.2020.118537
Rosas, T., Galiano, L., Ogaya, R., Peñuelas, J., Martínez-Vilalta, J. (2013). Dynamics of non-structural carbohydrates in three Mediterranean woody species following long-term experimental drought. Front. Plant Sci. 4. doi: 10.3389/fpls.2013.00400
Rubio-Cuadrado, Á., Camarero, J. J., del Río, M., Sánchez-González, M., Ruiz-Peinado, R., Bravo-Oviedo, A., et al. (2018). Long-term impacts of drought on growth and forest dynamics in a temperate beech-oak-birch forest. Agric. For. Meteorology 259, 48–59. doi: 10.1016/j.agrformet.2018.04.015
Seidel, D., Leuschner, C., Müller, A., Krause, B. (2011). Crown plasticity in mixed forests–quantifying asymmetry as a measure of competition using terrestrial laser scanning. For. Ecol. Manage. 261 (11), 2123–2132. doi: 10.1016/j.foreco.2011.03.008
Sevanto, S., Mcdowell, N. G., Dickman, L. T., Pangle, R., Pockman, W. T. (2014). How do trees die? a test of the hydraulic failure and carbon starvation hypotheses. Plant Cell Environ. 37 (1), 153–161. doi: 10.1111/pce.12141
Song, L., Zhu, J., Zhang, J., Wang, K., Lü, L., Wang, F., et al. (2019). Divergent growth responses to warming and drying climates between native and non-native tree species in northeast China. Trees 33 (4), 1143–1155. doi: 10.1007/s00468-019-01848-z
Sperry, J. S., Hacke, U. G., Pittermann, J. (2006). Size and function in conifer tracheids and angiosperm vessels. Am. J. Bot. 93 (10), 1490–1500. doi: 10.3732/ajb.93.10.1490
Stone, C., Kathuria, A., Carney, C., Hunter, J. (2008). Forest canopy health and stand structure associated with bell miners (Manorina melanophrys) on the central coast of new south Wales. Aust. Forestry 71 (4), 294–302. doi: 10.1080/00049158.2008.10675048
Sun, S., Qiu, L., He, C., Li, C., Zhang, J., Meng, P. (2018a). Drought-affected populus simonii carr. show lower growth and long-term increases in intrinsic water-use efficiency prior to tree mortality. Forests 9 (9), 564. doi: 10.3390/f9090564
Sun, S. J., He, C., Qiu, L., Li, C., Zhang, J. S., Meng, P. (2018b). Stable isotope analysis reveals prolonged drought stress in poplar plantation mortality of the three-north shelter forest in northern China. Agric. For. Meteorology 252, 39–48. doi: 10.1016/j.agrformet.2017.12.264
Tanaka-Oda, A., Kenzo, T., Koretsune, S., Sasaki, H., Fukuda, K. (2010). Ontogenetic changes in water-use efficiency (δ13C) and leaf traits differ among tree species growing in a semiarid region of the loess plateau, China. For. Ecol. Manage. 259 (5), 953–957. doi: 10.1016/j.foreco.2009.11.037
Tognetti, R., Cherubini, P., Marchi, S., Raschi, A. (2007). Leaf traits and tree rings suggest different water-use and carbon assimilation strategies by two co-occurring quercus species in a Mediterranean mixed-forest stand in Tuscany, Italy. Tree Physiol. 27 (12), 1741–1751. doi: 10.1093/treephys/27.12.1741
Trifilò, P., Casolo, V., Raimondo, F., Petrussa, E., Boscutti, F., Gullo, M. A. L., et al (2017). Effects of prolonged drought on stem non-structural carbohydrates content and post-drought hydraulic recovery in Laurus nobilis L.: The possible link between carbon starvation and hydraulic failure. Plant Physiol. Biochem. 120, 232–41. doi: 10.1016/j.plaphy.2017.10.003
Trouvé, R., Bontemps, J.-D., Collet, C., Seynave, I., Lebourgeois, F. (2017). Radial growth resilience of sessile oak after drought is affected by site water status, stand density, and social status. Trees 31 (2), 517–529. doi: 10.1007/s00468-016-1479-1
Urrutia-Jalabert, R., Malhi, Y., Barichivich, J., Lara, A., Delgado-Huertas, A., Rodríguez, C. G., et al. (2015). Increased water use efficiency but contrasting tree growth patterns in fitzroya cupressoides forests of southern Chile during recent decades. J. Geophysical Research: Biogeosciences 120 (12), 2505–2524. doi: 10.1002/2015JG003098
Vicente-Serrano, S. M., Beguería, S., López-Moreno, J. I. (2010). A multiscalar drought index sensitive to global warming: the standardized precipitation evapotranspiration index. J. Climate 23 (7), 1696–1718. doi: 10.1175/2009JCLI2909.1
von Arx, G., Dobbertin, M., Rebetez, M. (2012). Spatio-temporal effects of forest canopy on understory microclimate in a long-term experiment in Switzerland. Agric. For. Meteorol. 166, 144–155. doi: 10.1016/j.agrformet.2012.07.018
Wang, L., Dai, Y., Sun, J., Wan, X. (2017). Differential hydric deficit responses of robinia pseudoacacia and platycladus orientalis in pure and mixed stands in northern China and the species interactions under drought. Trees 31 (6), 2011–2021. doi: 10.1007/s00468-017-1605-8
Wigley, T. M., Briffa, K. R., Jones, P. D. (1984). On the average value of correlated time series, with applications in dendroclimatology and hydrometeorology. J. Appl. Meteorology Climatol. 23 (2), 201–213. doi: 10.1175/1520-0450(1984)0232.0.CO;2
Wu, X., Jiang, F., Li, X., Xue, Y., Qiu, S. (2004). Major features of decline of pinus sylvestris var. mongolica plantation on sandy land. J. Appl. Ecol. 15 (12), 2221–2224. doi: CNKI:SUN:YYSB.0.2004-12-004
Wu, G., Liu, X., Chen, T., Xu, G., Wang, W., Zeng, X., et al. (2015). Elevation-dependent variations of tree growth and intrinsic water-use efficiency in schrenk spruce (Picea schrenkiana) in the western tianshan mountains, China. Front. Plant Sci. 6. doi: 10.3389/fpls.2015.00309
Wullschleger, S. D., Wilson, K. B., Hanson, P. J. (2000). Environmental control of whole-plant transpiration, canopy conductance and estimates of the decoupling coefficient for large red maple trees. Agric. For. Meteorol. 104 (2), 157–168. doi: 10.1016/S0168-1923(00)00152-0
Zadworny, M., Jagodziński, A. M., Łakomy, P., Mucha, J., Oleksyn, J., Rodríguez-Calcerrada, J., et al. (2019). Regeneration origin affects radial growth patterns preceding oak decline and death–insights from tree-ring δ13C and δ18O. Agric. For. Meteorology 278, 107685. doi: 10.1016/j.agrformet.2019.107685
Zas, R., Sampedro, L., Solla, A., Vivas, M., Lombardero, M. J., Alía, R., et al. (2020). Dendroecology in common gardens: Population differentiation and plasticity in resistance, recovery and resilience to extreme drought events in pinus pinaster. Agric. For. Meteorol. 291, 108060. doi: 10.1016/j.agrformet.2020.108060
Zhang, L. M., Shang, Z. P. (2002). Relationship between the soil moisture and the vegetation productivity in the loess plateau. Arid Zone Res. (Chinese Version) 19 (4), 59–63. doi: 10.1007/s11769-002-0038-4
Zhang, X., Yu, P., Wang, D., Xu, Z. (2022). Density-and age-dependent influences of droughts and intrinsic water use efficiency on growth in temperate plantations. Agric. For. Meteorol. 325, 109134. doi: 10.1016/j.agrformet.2022.109134
Keywords: drought stress, radial growth, resistance, stable isotopes, canopy die-back
Citation: Jia H, Guan C, Zhang J, He C, Yin C and Meng P (2022) Drought effects on tree growth, water use efficiency, vulnerability and canopy health of Quercus variabilis-Robinia pseudoacacia mixed plantation. Front. Plant Sci. 13:1018405. doi: 10.3389/fpls.2022.1018405
Received: 13 August 2022; Accepted: 22 September 2022;
Published: 12 October 2022.
Edited by:
Fenghui Yuan, University of Minnesota Twin Cities, United StatesReviewed by:
Hongxia Zhang, Northwest Institute of Eco-Environment and Resources (CAS), ChinaKayla Stan, University of Alberta, Canada
Benye Xi, Beijing Forestry University, China
Copyright © 2022 Jia, Guan, Zhang, He, Yin and Meng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: JinSong Zhang, emhhbmdqc0BjYWYuYWMuY24=