- 1College of Water Resources and Architectural Engineering, Northwest A&F University, Yangling, China
- 2College of Science, Northwest A&F University, Yangling, China
- 3Department of Soil Science, University of Saskatchewan, Saskatoon, SK, Canada
Introduction: Parameterizing the process of trees from the comfort zone to mortality during progressive drought is important for, but is not well represented in, vegetation models, given the lack of appropriate indices to gauge the response of trees to droughts. The objective of this study was to determine reliable and readily available tree drought stressindices and the thresholds at which droughts activate important physiological responses.
Methods: We analyzed the changes in the transpiration (T), stomatal conductance, xylem conductance, and leaf health status due to a decrease in soil water availability (SWA), predawn xylem water potential (ψpd), and midday xylem water potential (ψmd) in Robinia pseudoacacia seedlings during progressive drought.
Results: The results showed that ψmd was a better indicator of drought stress than SWA and ψpd, because ψmd was more closely related to the physiological response (defoliation and xylem embolization) during severe drought and could be measured more conveniently. We derived the following five stress levels from the observed responses to decreasing ψmd: comfort zone (ψmd > -0.9 MPa), wherein transpiration and stomatal conductance are not limited by SWA; moderate drought stress (-0.9 to -1.75 MPa), wherein transpiration and stomatal conductance are limited by drought; high drought stress (-1.75 to -2.59 MPa), wherein transpiration decreases significantly (T< 10%) and stomata closes completely; severe drought stress (-2.59 to -4.02 MPa), wherein transpiration ceases (T< 0.1%) and leaf shedding orwilting is > 50%; and extreme drought stress (< -4.02 MPa), leading to tree mortality due to xylem hydraulic failure.
Discussion: To our knowledge, our scheme is the first to outline the quantitative thresholds for the downregulation of physiological processes in R. pseudoacacia during drought, therefore, can be used to synthesize valuable information for process-based vegetation models.
1 Introduction
Drought is one of the principal stresses that affect terrestrial ecosystems carbon balance with the potential to cause forest decline and mortality at regional and continental scales (Huxman et al., 2004; Zhao and Running, 2010; Choat et al., 2012; Reichstein et al., 2013; Huang et al., 2016). Global climate change is expected to intensify the regional-scale droughts (Dai, 2012; Barkhordarian et al., 2019; Yuan et al., 2019), posing a serious threat to forest productivity and growth (Anderegg et al., 2019; Taccoen et al., 2019). Accurate predictions of future drought-induced forest dynamics require reliable data on temporal and spatial variations in environmental factors as well as key physiological processes. The latter determines species productivity and mortality during drought and represents crucial plant characteristics or appropriate thresholds in vegetation models (Choat et al., 2018). It is urgent to understand and quantify the decline in the physiological functions of trees and the corresponding thresholds during drought (Venturas et al., 2018; Nadal-Sala et al., 2021), especially the processes related to productivity and mortality. However, there is a lack of appropriate indices for quantifying the response of trees to drought.
Soil water availability (SWA) and plant water potential have been used as indicators to determine threshold values for various physiological responses of plants to drought, including stomatal closure, tissue damage, hydraulic failure, and other physiological processes (Bolte et al., 2016; Yan et al., 2017b; Trueba et al., 2019; Walthert et al., 2021). However, because of the spatial heterogeneity of soil hydraulic properties (Bormann and Klaassen, 2008; Lu et al., 2019) and the large spatial and temporal variabilities in root distribution and water absorption characteristics (Christina et al., 2017; Li et al., 2019), accurate measurement of SWA in the root zone of trees may require huge financial and human resources, especially for trees with deep roots. The existing conceptual framework states that plant water potential is the most straightforward indicator of plant water status as it integrates the effects of soil, plant, and atmospheric conditions (Jones, 2004; Suter et al., 2019). Walthert et al. (2021) used predawn leaf water potential to quantify the critical thresholds of stem growth ceases, xylem embolism, and crown dieback in the European beech during drought. However, predawn water potential only reflects soil water status of plants under conditions of none-to-very little atmospheric evaporation demand rather than that of maximum evaporation demand (Scholander et al., 1965; Suter et al., 2019; Tian and Schreiner, 2021).
Tree growth status is not only affected by SWA, but also by the atmospheric environment (Weithmann et al., 2021; Tumajer et al., 2022). Therefore, it is more accurate to use an index that includes comprehensive impact information considering the plants, soil, and atmosphere, when quantifying the relationship between plant physiological response and drought degree. Midday plant water potential (including the xylem and leaf water potentials) reflects water stress experienced by plants at the maximum evaporation demand (Pilar et al., 2007). However, the leaf water potential only represents water potential of a single leaf and may not always reflect the true water condition of the whole plant, whereas the xylem water potential indicates the water condition of all of the leaves (Jones, 2007). Although leaf water potential has been successfully used by some researchers, the advantages of xylem water potential have been demonstrated in recent studies (Moriana et al., 2012; Cole and Pagay, 2015; Santesteban et al., 2019; Suter et al., 2019). For example, the midday stomatal conductance and transpiration are positively correlated with midday stem water potential rather than midday leaf water potential (Choné, 2001; Zhang et al., 2013). The midday xylem water potential (ψmd) has also been used in a few studies to explore the hydraulic failure threshold of some species (Urli et al., 2013; Hammond et al., 2019).
Robinia pseudoacacia is a widespread and ecologically important deciduous tree (Schwarzel et al., 2019; Puchalka et al., 2021), but its production thresholds (stomatal closure and leaf loss) and survival threshold (hydraulic failure) during drought are still difficult to be determined. This gap has hampered our understanding of the survival of R. pseudoacacia during drought and the development of ecological models in predicting tree growth status. It may take years or even decades to complete the measurements of adult trees that go from no water stress to drought-caused death in the wild under non-extreme climatic conditions. Although there may be some variation (e.g., embolism resistance) in the physiological response of trees to drought with changes in tree age (Domec and Gartner, 2003; Rodriguez-Zaccaro et al., 2019), the age-related variation may be much smaller than differences between species. Therefore, we studied the physiological processes of R. pseudoacacia seedlings during drought grown in a greenhouse to collect important information for predicting their growth state under different growing environments and climatic conditions.
We measured the changes in stomatal conductance (gs), transpiration rate (T), percentage loss of xylem hydraulic conductivity (PLC), and total leaf area (AL) in the seedlings due to a decrease in SWA, ψpd, and ψmd during droughts. The purpose of this study was to explore: 1) which index of SWA, ψpd or ψmd could better reflect the drought stress of trees; and 2) the thresholds corresponding to production (stomatal closure and leaf shedding) and survival (hydraulic failure) during the entire drought process.
2 Materials and methods
2.1 Plant materials and experimental design
We obtained 2-year-old Robinia pseudoacacia saplings (n =26, height: 30 - 50 cm) from a local nursery (Changwu, Shaanxi, China), and then transplanted each to a 18.5 L pot filled with 60% of clay loam (particle size 0 - 2 mm), 20% of coarse river sand (particle size 2 - 5 mm), and 20% of peat moss mixed with 2.5 kg m−3 slow-release fertilizer (Osmocote, NPK 15:10:12 + 2MgO + micro elements) in April 2021. The soil surfaces were covered with 3-cm-thick perlite to reduce evaporation. All of the trees were watered well after transplantation, and grown for two months to achieve acclimation to the glasshouse before setting different treatments and conducting the measurements. None of the trees died during the two months after transplantation.
All of the saplings were exposed to a 2-week pre-experimental drought treatment, by withholding water until stomatal closure approached 80%, to avoid subjecting the tree saplings that had never experienced water stress to our experiments. No trees died during the acclimation to drought. The trees were then watered well for two weeks before starting the measurements. At the beginning of the drought experiment in July 2021, the sapling heights were 1.8 ± 0.12 m (mean ± standard deviation), and the stem diameters at atmosphere-soil interface were 5.00 ± 0.24 cm. Two different water regimes were applied after three months of growth to allow acclimation. Drought was induced by withholding irrigation for 20 saplings in the water-stressed treatment group, whereas the remaining six control saplings were watered well to maintain the soil water status at field capacity.
2.2 Index measurements
The xylem water potential, stomatal conductance, transpiration rate, and conductivity from the rhizosphere to xylem were measured on sunny leaves from the upper crown of the selected plants at least twice a week in the early drought period and once a day in the late drought period. These indices were measured on a total of 26 trees, including six control trees. The above eco-physiological measurements were conducted until leaf desiccation (no more leaves, complete necrosis, or withering).
2.2.1 Xylem water potential and stomatal conductance
The xylem water potential was measured using a Scholander pressure bomb (model 1505D-Exp; PMS Instruments, Albany, OR, USA). Predawn and midday xylem water potentials (ψpd and ψmd) were measured at 4:00 - 6:00 and 13:00 - 15:00, respectively. The leaves that were used to measure ψmd were wrapped with tin foil paper for 1 h to equilibrate the water potential between the xylem and leaf.
The stomatal conductance was measured on at least four selected leaves per plant at 09:00 - 11:00 using a porometer (AP-4; Delta-T Devices Ltd., Cambridge, UK). Fully expanded mature leaves from the upper crowns were selected for stomatal conductance measurements during progressive dehydration, and the adjacent leaves were used for water potential measurements.
2.2.2 Transpiration and conductivity from rhizosphere to xylem
The plants were weighed to an accuracy of ± 0.1 g (Langge R30-A01, Langge Technology, Beijing, China) between 11:00 and 13:00. For each plant, transpiration rate (T, g m−2 h−1) was calculated by the loss of weight of each plant, as follows (Urli et al., 2013):
where Δw/Δt (g h−1) is the loss of water over a given time interval (60 - 120 min) and Total leaf area (AL, m2) is the product of the average area of a single leaf and the total number of leaves per plant. Leaf area was measured on at least 10 leaf samples using a flatbed scanner (Epson model V700).
The hydraulic conductivity from the rhizosphere to xylem at midday (KRX, mmol m−2 s-1 MPa-1) was calculated based on the assumption that soil water potential equals to ψpd. Under these circumstances (Venturas et al., 2018):
where T (g m−2 h−1) is the transpiration, ψmd (MPa) and ψpd (MPa) are predawn and midday xylem water potentials, respectively. The ψmd was measured immediately after the transpiration was measured.
2.2.3 Percentage loss of xylem hydraulic conductivity and vulnerability curve
The saplings were transported to a nearby laboratory for the measurement of hydraulic conductivity. In the laboratory, stem segments over 30 cm long were cut from the saplings and their bases were immediately placed in filtered KCl solution (0.2 μm, 0.01 mol L-1). The initial excision segments were recut approximately 2 cm from the cut base and placed under a filtered solution for at least 30 min to release the tension. The segments, approximately 28 cm in length, were excised and mounted on a conductivity system filled with deionized, filtered (0.2 μm), and degassed 0.01 mol L-1 KCl solution. The initial conductivity (Ki) was measured by flowing a KCl solution from a reservoir through the segment and onto a computer-interfaced balance with a pressure difference of approximately 3 kPa. The stems were then flushed at a pressure of 150 kPa for 20 min to remove air bubbles. The hydraulic conductivity was determined again, and flushing was repeated until reaching the maximum conductivity (Kmax). The percentage loss of conductivity (PLC, %) was computed as follows:
The xylem vulnerability curve measurements were conducted according to the bench-top dehydration method. Sample shoots were collected early morning before sunrise and dehydrated on a bench in the laboratory to obtain the PLC and corresponding xylem tension. The measured shoot was wrapped in a black plastic bag to equilibrate the water potential between the xylem and leaf when the shoot had reached the desired tension. After a minimum of 1 h of equilibration, three leaves were excised from the wrapped shoots. The equilibrated leaf water potential was measured using a pressure chamber (model 1505D; PMS Instruments, OR, USA). Soon after excising the leaves, the PLC due to embolism was measured in three 10 cm long stem segments. The xylem vulnerability curve was obtained by fitting a single Weibull function to the measured PLC value at each pressure (Pammenter and Vander Willigen, 1998).
2.2.4 Soil water availability
Soil water availability (SWA, %) was expressed as (Granier et al., 2007):
where θA is the actual water content (g g-1), θFC and θPWP are the soil water contents (g g-1) at field capacity (FC, − 0.06 MPa soil water potential) and the permanent wilting point (PWP, − 1.5 MPa soil water potential), respectively. θFC and θPWP were obtained from the soil water characteristic curve of the measured soil water potential against soil water content. Soil water potential was measured using a water potential probe buried 10 cm deep in the pot. Soil water content was obtained using the oven-dry method. To determine the SWA, each pot was weighed three times per week until seedlings were rewatered.
2.3 Data analyses
To eliminate the error caused by environmental conditions during the measurement days, we normalized the stomatal conductance and transpiration data by dividing the data for each drought plant by the means of the control plants measured on the same day. The lethal threshold of hydraulic failure was defined as the percentage of xylem hydraulic conductivity loss associated with a 50% mortality rate in the population. The interpretation degree of drought indices (SWA, ψpd, and ψmd) to related physiological traits were expressed by determination coefficient. When analyzing the correlation between drought indices, the Ψpd was calculated by the fitting model between SWA and Ψpd, and the Ψmd was replaced by the Ψ of vulnerability curve measured by the table dehydration method. Because the leaves had completely dropped off, the Ψpd and Ψmd could not be measured when the xylem conductance loss reached >50%. All of the statistical analyses were performed using the SPSS software (IBM SPSS Statistics for Windows, Version 22.0; IBM Corp., Armonk, NY, USA), and p<0.05 was considered significant.
3 Results
3.1 Midday xylem water potential is better than soil water availability and predawn water potential in quantifying drought stress degree
Figure 1 shows changes of soil water availability (SWA) with days after withholding irrigation and the relationship between SWA and Ψpd. The SWA decreased rapidly within the first 15 days of drought treatment from 100% to 14%, and then continued to decrease gradually to the lowest 6.6% at 24th day (Figure 1A). The SWA was closely correlated with the predawn water potentials Ψpd, and the Ψpd decreased exponentially with a decrease in the SWA (R2 = 0.92, p< 0.01, Figure 1B). Ψpd was −0.47 MPa at 100% SWA and decreased to −2.59 MPa at 6.6% SWA.
Figure 1 Variation in soil water availability (SWA) on drought days (A) and relationship between predawn xylem water potential (Ψpd) and SWA (B). The SWA values were the mean values of 20 individuals, and the bars indicate the standard error on the mean value. The fitted curve is an exponential function with statistical significance (p< 0.01, **). Each value corresponds to a single measurement.
The standardized stomatal conductance (gs) (10:00–11:00) and the transpiration (T) (11:00–13:00) loss percentage increased with the decline of SWA, Ψpd and Ψmd. The gs and T were significantly correlated with SWA, Ψpd and Ψmd (p< 0.01). The coefficients of determination of gs and T for SWA, Ψpd and Ψmd ranged from 0.85 to 0.92, and the difference between them was very small (Figure 2).
Figure 2 Response of stomatal conductance (gs) and transpiration rate (T) to soil water availability (100 - SWA) (A, B), predawn xylem water potential (Ψpd) (C, D), and midday xylem water potential (Ψmd) (E, F). These data were fitted with a single Weibull function. The curve depicts modeled percent loss of gs, and T, and the orange and red dotted vertical lines depict estimates of water potentials inducing 50% (gs50, T50) and 90% (gs90, T90) loss of gs and T. The blue dotted vertical line depict estimate of water potential inducing 10% (T10) loss of T. Significance codes: **(p< 0.01).
The R. pseudoacacia leaves began to wilt and shed when drought stress reached a certain degree, and the proportion of wilted and shed leaves gradually increased as the drought worsened. The proportion of withered leaves was negatively correlated with Ψpd and Ψmd (p< 0.05, r = 0.52; p< 0.01, r = 0.66) and was not significantly correlated with the SWA (p > 0.05) (Figure 3). The total leaf area decreased by 50% when the Ψmd reached −2.59 MPa (Figure 3C).
Figure 3 Percentage of fallen leaves (relative to that before the drought) as a function of soil water availability (SWA) (A), predawn xylem water potential (Ψpd) (B), and midday xylem water potential (Ψmd) (C). Points depict pooled data from individuals (n = 20). The xylem water potential at which the total leaf area reached 50% was calculated as the point at which the regression line intersects 50% of fallen leaves number (dotted line). Significance codes: asterisks indicate the level of significance (**, P< 0.01; ns, nonsignificant relationships).
The percentage of xylem conductance loss (PLC) had no significant correlation with the SWA or ΨPd (p > 0.05) (Figures 4A, B). In other words, the PLC did not increase significantly with decreasing SWA and ΨPd as expected. However, the PLC increased as the Ψmd decreased and was close to 90% at −6.16 MPa and almost 10% at −1.2 MPa (Figure 4C).
Figure 4 Percentage xylem conductance loss (PLC) as a function of soil water availability (SWA) (A), predawn xylem water potential (Ψpd) (B), and midday xylem water potential (Ψmd) (C). The blue, orange, and red vertical dash dot lines depict estimates of water potentials that represent the water potential thresholds for xylem embolization initiation (Pe), xylem conductance loss of 50% (P50), and hydraulic failure (Phf), respectively. Significance codes: asterisks indicate the level of significance (**, P< 0.01; ns, nonsignificant relationships). (The SWA was measured by oven-dry method when the xylem conductance loss was measured. The Ψpd was calculated by fitting the model between SWA and Ψpd (Figure 1B). The Ψmd was replaced by the Ψ of vulnerability curve measured by the table dehydration method. Because the leaves had completely dropped off, the Ψpd and Ψmd could not be measured when the xylem conductance loss reached >50%. The PLC was associated with the lowest xylem water potential experienced by the plant, and the Ψmd is generally considered to be the lowest water potential a plant experiences throughout the day.).
3.2 Thresholds of key physiological response in the shift from no stress to hydraulic failure
When stomatal conductance decreased by 50% (gs50) and 90% (gs90), the midday xylem water potential was −1.35 ± 0.32 MPa and −1.75 ± 0.32 MPa (mean ± standard deviation), respectively (Figure 2E). The difference between the values of the two nodes was only 0.4 MPa. Compared with the control group, the Ψmd of T decreased by 10% (T10), 50% (T50) and 90% (T90 due to drought were −0.89 ± 0.33 MPa, −1.40 ± 0.33 MPa and −1.87 ± 0.33 MPa, respectively (Figure 2F). The Ψmd was −2.59 ± 0.72 MPa and −4.10 ± 0.72 MPa when the leaf wilting or shedding rate reached 50% and 90%, respectively (Figure 4C). The Ψmd for the onset of xylem embolization (Pe), loss of hydraulic conductivity 50% (P50), and hydraulic failure (Phf) was −1.29 ± 0.34 MPa, −2.95 ± 0.34 MPa and −4.02 ± 0.34 MPa, respectively (Figure 4C).
The rhizosphere-to-xylem conductivity (KRX) decreased significantly with an increase in the Ψmd, which was closer to the R shape curve than to the xylem vulnerability curve. Compared with the control group, the Ψmd with 50% and 90% KRX loss was −1.35 ± 0.34 MPa and −2.58 ± 0.34 MPa (mean ± standard deviation), respectively (Figure 5). Notablely, the xylem water potential with a 50% decrease in stomatal conductance was the same as that with a 50% decrease in conductivity from the rhizosphere to xylem. In addition, the KRX decreased with the decrease in Ψmd much faster than the xylem conductivity.
Figure 5 Response of normalized rhizosphere-to-xylem conductivity (KRX) loss to midday xylem water potential (Ψmd). Points depict pooled data from individuals (n = 20). The curve depicts modeled percent loss of gs KRX, and the orange and red dotted vertical lines depict estimates of water potentials inducing 50% (KRX50) and 90% (KRX90) loss of KRX. Significance codes: **(p<.01).
Based on the observed responses to decreasing Ψmd, we derived the following five stress levels (Figure 6): comfort zone (Ψmd > −0.9 MPa), loss of stomatal conductance and transpiration rate was less than 10%, that is, transpiration and stomatal conductance were not limited by water availability; moderate drought stress (−0.9 to 1.75 MPa), loss of stomatal conductance and transpiration rate reached 90%, that is, stomatal conductance and transpiration were limited by drought; high drought stress (−1.75 to −2.59 MPa), transpiration rate decreased significantly (T< 10%) and stomata closed completely; severe drought stress (−2.59 to −4.02 MPa), transpiration ceased (T< 0.1%) and leaf shedding or wilting reached > 50%; extreme drought stress (< −4.02 MPa), tree mortality (p > 50%) due to xylem hydraulic failure.
Figure 6 Schematic synthesis of physiological responses of Robinia pseudoacacia seedlings to progressive drought. The derivation of the five stress levels and positioning of the thresholds were determined based on important physiological responses under decreasing midday xylem water potential (Ψmd).
4 Discussion
Our results showed that soil water availability (SWA), predawn xylem water potential (Ψpd), and midday xylem water potential (Ψmd) of R. pseudoacacia were significantly correlated with the maximum stomatal conductance and T during the early drought period (Figure 2). In other words, these three indicators can indicate the degree of drought stress faced by R. pseudoacacia in the early stage, and the performance of Ψpd and Ψmd were better than SWA. However, the correlation between leaf wilting or shedding rate and Ψmd was stronger during the late drought period. In addition, there was no significant correlation between the xylem hydraulic conductivity loss (PLC), Ψpd, and SWA, but the correlation between the PLC and the lowest xylem water potential was very high. It is worth noting that the midday water potential of plants is closer to the lowest water potential of the day than the water potentials at other times (Urli et al., 2013; Hammond et al., 2019). Midday xylem water potential has also been used in past studies on the hydraulic failure threshold of trees, which is an important hydraulic characteristic (Urli et al., 2013; Hammond et al., 2019). Therefore, according to the variation in stomatal behavior, transpiration rate, leaf health status and xylem embolization with the changes in the SWA, Ψpd, and Ψmd during the whole drought process, it can be concluded that Ψmd is the best indicator of the degree of drought stress faced by R. pseudoacacia during the whole drought process.
Reduction in the transpiration (T), root water uptake rate and stem contraction are the first manifestations of trees under drought stress, and these processes occur almost simultaneously (Walthert et al., 2021). Therefore, we defined the starting point of drought stress as the point at which the T began to decrease significantly (T10). Stomatal closure is the primary mechanism that limits excessive water loss and xylem tension accumulation in trees under drought stress (Choat et al., 2018). In many studies, a 50% and 90% decrease in stomatal conductance was defined as the beginning and complete closure of the stomata, respectively (Martin-StPaul et al., 2017; Li et al., 2018). The leaf photochemical apparatus and related enzyme activities are not significantly damaged before stomatal closure in the process of dehydration caused by drought, and stomatal conductance is the main limiting factor of photosynthetic carbon assimilation during the early stage of drought (Yan et al., 2017b; Trueba et al., 2019; Zhu et al., 2021). Thus, the water potential at complete stomatal closure can be considered the first node at which tree productivity is limited by drought stress.
Leaf wilting or shedding can effectively reduce water vapor leakage in the cuticle (Nadal-Sala et al., 2021), protect the hydraulic security of the core organ xylem, and delay the time for trees to reach hydraulic failure during drought (Blackman et al., 2019). This is considered another important drought-avoidance hydraulic strategy of deciduous and semi-deciduous trees at the expense of growth and productivity, in addition to stomatal closure because rebuilding the canopy structure requires additional carbon investment, whether from non-structural carbohydrate reserves or the assimilation of newly grown leaves after drought stress relief (Yan et al., 2017a). Therefore, the water potential at leaf wilting or shedding can be considered the node where the photosynthetic production capacity of trees is completely lost during drought. The threshold of hydraulic failure, defined as the percentage of xylem conductance loss when the community mortality exceeds 50% (Hammond et al., 2019), is an important parameter for predicting the probability of tree mortality caused by drought (Choat et al., 2018). The vulnerability curve is correlated to the loss of xylem hydraulic conductivity loss due to the decline in xylem water potential (Cochard et al., 2013), so the death node of trees in drought can be defined by the water potential corresponding to the wood hydraulic failure threshold. Therefore, the midday xylem water potentials when the reduced transpiration rate reduced by 10%, stomatal closure by 90%, leaf wilting or shedding by 50%, and hydraulic failure were selected as the nodes for starting drought, moderate drought, severe drought, and extreme drought stress in R. pseudoacacia, respectively.
The transpiration rate and leaf stomatal conductance of R. pseudoacacia showed similar trends with xylem water potential, and the Ψmd of the beginning of stomatal closure (gs50) and complete closure (gs90) were very close to the Ψmd of 50% (T50) and 90% (T90), respectively (Figure 6). This indicated that stomata controlled the main transpiration of R. pseudoacacia seedlings. stomatal closure (gs50) was very close to the water potential at the beginning of xylem embolization (Pe) and is quite different from the water potential at hydraulic failure (Pfh), suggesting that stomatal closure may play a protective role in inhibiting rapid reduction of water potential and resisting catastrophic embolism. However, stomatal closure slows down tree water consumption and xylem water potential drop but leads to a rapid decrease in photosynthesis (Venturas et al., 2018). Low photosynthetic rates associated with stomatal closure can lead to the depletion of non-structural carbohydrate pools, which interfere with the translocation of sugars through the phloem and production of chemical defense compounds needed to prevent herbivory and disease over long time periods (McDowell et al., 2008; Choat et al., 2018). This behavior of R pseudoacacia greatly reduces the risk of xylem hydraulic failure during short-term drought but significantly increases the threat of pests and diseases caused by non-structural carbon depletion during long-term drought.
Our results showed that the water potential of R. pseudoacacia at a 50% decrease in the transpiration rate was very close to that at a 50% loss in the rhizosphere-to-xylem conductivity. The rhizosphere-to-xylem conductivity decreased more rapidly than the xylem conductivity with xylem water potential, indicating that the rhizosphere was more vulnerable than xylem during drought. Rhizosphere conductivity includes the conductivity of the root system, root-soil interface, and soil. In many angiosperms, the roots are not more susceptible to embolism than the xylem (Rodriguez-Dominguez et al., 2018; Wu et al., 2020). The water potential of the roots will not be lower than that of the xylem during drought. Therefore, the conductivity loss of roots should be less than that of xylem, and the root should not be a limiting factor in the rhizosphere to xylem conductivity during drought. The difference between the hydraulic conductivity of various soils and xylem conductance of trees is more than 5 orders of magnitude, with the soil water potential ranging from −0.06 to −1.5 MPa (Carminati and Javaux, 2020), that is, xylem conductance is much greater than soil hydraulic conductivity. Therefore, it is reasonable to speculate that the transpiration and productivity of R. pseudoacacia were first affected by the conductance loss of the root-soil interface or soil rather than that of the root and xylem during early drought (Carminati and Javaux, 2020; Rodriguez-Dominguez and Brodribb, 2020). Calculation based on Sperry’s model showed that the restriction of transpiration by the rhizosphere or xylem depended on the root-shoot ratio and soil texture, that is, if the root-shoot ratio is very large, the restriction is mainly caused by the xylem, and if the root-shoot ratio is not very large (< 40) and under the condition of fine soil texture, the restriction of tree transpiration is mainly caused by the rhizosphere (Sperry et al., 1998).
Plant transpiration was limited by rhizosphere depending on both soil textures and root hydraulic phenotypes in drying soils, a root phenotype with low root hydraulic conductance, long and dense fine root postpones rhizosphere limitation, and coarse textured soils exhibit an earlier limitation to root water uptake due to their sudden loss in conductivity for decreasing soil water potentials (Cai et al., 2022). Additionally, if roots are in dry soil for a long time, the suberisation degree will increase and a large number of nests will be formed, which will affect the roots water absorbing capacity (Moura et al., 2010; Cuneo et al., 2016; Zhang et al., 2020). Even when the soil is rewet, it is difficult for the rhizosphere conductance to quickly recover to the original level, because the structural change of the root is irreversible (Steudle, 2000; Cuneo et al., 2016). For deep-rooted plant with two water sources (deep soil water and precipitation), the ratio of fine roots in deep soil to shallow soil also affects whether transpiration is limited by the rhizosphere during drought (Oliveira et al., 2005; Li et al., 2019). However, the influence of drought on transpiration may be delayed if there is hydraulic redistribution in the root zone (Bleby et al., 2010; Domec et al., 2010). That’s because trees respond to drought by changing water uptake among existing roots, rather than growing new ones (Mackay et al., 2020). Therefore, further field experiments are needed to demonstrate the specific process of rhizosphere influence on tree transpiration in natural arid environment.
5 Conclusions
The results of this study indicate that ψmd can be used as a reliable quantitative index of the degree of drought stress in trees. The R. pseudoacacia seedlings experienced rhizosphere conductivity decreased, stomatal closure, transpiration rate decreased, xylem conductivity decreased, leaf wilted, and hydraulic failure with the ψmd decrease from comfort zone to drought death. The transpiration rate is mainly affected by stomatal conductance, and stomatal conductance is primarily affected by the decrease in rhizosphere conductivity rather than xylem conductivity. Therefore, we should not focus only on the aboveground parts of plants, but consider the hydraulic characteristics of the aboveground and belowground parts, especially the soil and the connection between roots and soil, when using ecological models to study the response of tree transpiration or photosynthetic production to drought.
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 authors.
Author contributions
XW: experiment design, data collection and analysis, visualization, and original manuscript writing and editing. YF: data collection, visualization, and manuscript editing. CZ: data collection and manuscript editing. YZ: provided substantial inspiration and revision. GD: date validation, statistical analysis, and manuscript editing and revision. ML: date validation, statistical analysis, funding acquisition, and original manuscript editing and revision. BS: experiment design, funding acquisition, data validation and analysis, and manuscript editing and revision. All authors contributed to the article and approved the submitted version.
Funding
This research was partially supported by the National Key R&D Program of China (2022YFE0100300), the National Natural Science Foundation of China (41630860, 41877017), and the 111Project (B12007).
Acknowledgments
We thank Siyi Chen, and Xiaoyun Ding for their invaluable assistance in collecting root and soil field data.We appreciate the technical help from Jingjing Jin, Institute of Water-Saving Agriculture in Arid Areas of China, Northwest A&F University.
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.
References
Anderegg, W. R. L., Anderegg, L. D. L., Kerr, K. L., Trugman, A. T. (2019). Widespread drought-induced tree mortality at dry range edges indicates that climate stress exceeds species' compensating mechanisms. Glob Chang Biol. 25 (11), 3793–3802. doi: 10.1111/gcb.14771
Barkhordarian, A., Saatchi, S. S., Behrangi, A., Loikith, P. C., Mechoso, C. R. (2019). A recent systematic increase in vapor pressure deficit over tropical south America. Sci. Rep. 9 (1), 15331. doi: 10.1038/s41598-019-51857-8
Blackman, C. J., Li, X, Choat, B., Rymer, P. D., De Kauwe, M. G., Duursma, R. A., et al. (2019). Desiccation time during drought is highly predictable across species of eucalyptus from contrasting climates. New Phytol. 224 (2), 632–643. doi: 10.1111/nph.16042
Bleby, T. M., McElrone, A. J., Jackson, R. B. (2010). Water uptake and hydraulic redistribution across large woody root systems to 20 m depth. Plant Cell Environ. 33 (12), 2132–2148. doi: 10.1111/j.1365-3040.2010.02212.x
Bolte, A., Czajkowski, T., Cocozza, C., Tognetti, R., de Miguel, M., Psidova, E., et al. (2016). Desiccation and mortality dynamics in seedlings of different European beech (Fagus sylvatica l.) populations under extreme drought conditions. Front. Plant Sci. 7, 751. doi: 10.3389/fpls.2016.00751
Bormann, H., Klaassen, K. (2008). Seasonal and land use dependent variability of soil hydraulic and soil hydrological properties of two northern German soils. Geoderma 145 (3-4), 295–302. doi: 10.1016/j.geoderma.2008.03.017
Cai, G. C., Ahmed, M. A., Abdalla, M., Carminati, A. (2022). Root hydraulic phenotypes impacting water uptake in drying soils. Plant Cell Environ 45 (3), 650–663. doi: 10.1111/pce.14259
Carminati, A., Javaux, M. (2020). Soil rather than xylem vulnerability controls stomatal response to drought. Trends Plant Sci. 25 (9), 868–880. doi: 10.1016/j.tplants.2020.04.003
Choat, B., Brodribb, T. J., Brodersen, C. R., Duursma, R. A., Lopez, R., Medlyn, B. E., et al. (2012). Global convergence in the vulnerability of forests to drought. Nature 491 (7426), 752–755. doi: 10.1038/s41586-018-0240-x
Choat, B., Jansen, S., Brodribb, T. J., Cochard, H., Delzon, S., Bhaskar, R., et al. (2018). Triggers of tree mortality under drought. Nature 558 (7711), 531–539. doi: 10.1038/nature11688
Choné, X. (2001). Stem water potential is a sensitive indicator of grapevine water status. Ann. Bot. 87 (4), 477–483. doi: 10.1006/anbo.2000.1361
Christina, M., Nouvellon, Y., Laclau, J. P., Stape, J. L., Bouillet, J. P., Lambais, G. R., et al. (2017). Importance of deep water uptake in tropical eucalypt forest. Funct. Ecol. 31 (2), 509–519. doi: 10.1111/1365-2435.12727
Cochard, H., Badel, E., Herbette, S., Delzon, S., Choat, B., Jansen, S. (2013). Methods for measuring plant vulnerability to cavitation: A critical review. J. Exp. Bot. 64 (15), 4779–4791. doi: 10.1093/jxb/ert193
Cole, J., Pagay, V. (2015). Usefulness of early morning stem water potential as a sensitive indicator of water status of deficit-irrigated grapevines (Vitis vinifera l.). Sci. Hortic. 191, 10–14. doi: 10.1016/j.scienta.2015.04.034
Cuneo, I. F., Knipfer, T., Brodersen, C. R., McElrone, A. J. (2016). Mechanical failure of fine root cortical cells initiates plant hydraulic decline during drought. Plant Physiol. 172 (3), 1669–1678. doi: 10.1104/pp.16.00923
Dai, A. (2012). Increasing drought under global warming in observations and models. Nat. Clim Change 3 (1), 52–58. doi: 10.1038/nclimate1633
Domec, J. C., King, J. S., Noormets, A., Treasure, E., Gavazzi, M. J., Sun, G., et al. (2010). Hydraulic redistribution of soil water by roots affects whole-stand evapotranspiration and net ecosystem carbon exchange. New Phytol. 187 (1), 171–183. doi: 10.1111/j.1469-8137.2010.03245.x
Domec, J. C., Gartner, B. L. (2003). Relationship between growth rates and xylem hydraulic characteristics in young, mature and old-growth ponderosa pine trees. Plant Cell Environ. 26 (3), 471–483. doi: 10.1046/j.1365-3040.2003.00978.x
Granier, A., Reichstein, M., Bréda, N., Janssens, I. A., Falge, E., Ciais, P., et al. (2007). Evidence for soil water control on carbon and water dynamics in European forests during the extremely dry year: 2003. Agr For. Meteorol 143 (1-2), 123–145. doi: 10.1016/j.agrformet.2006.12.004
Hammond, W. M., Yu, K. L., Wilson, L. A., Will, R. E., Anderegg, W. R. L, Adams, H. D. (2019). Dead or dying? Quantifying the point of no return from hydraulic failure in drought-induced tree mortality. New Phytol. 223 (4), 1834–1843. doi: 10.1111/nph.15922
Huang, L., He, B., Chen, A. F., Wang, H. Y., Liu, J. J., Lu, A. F., et al. (2016). Drought dominates the interannual variability in global terrestrial net primary production by controlling semi-arid ecosystems. Sci. Rep. 6, 6. doi: 10.1038/srep24639
Huxman, T. E., Smith, M. D., Fay, P. A., Knapp, A. K., Shaw, M. R., Loik, M. E., et al. (2004). Convergence across biomes to a common rain-use efficiency. Nature 429 (6992), 651–654. doi: 10.1038/nature02561
Jones, H. G. (2004). Irrigation scheduling: Advantages and pitfalls of plant-based methods. J. Exp. Bot. 55 (407), 2427–2436. doi: 10.1093/jxb/erh213
Jones, H. G. (2007). Monitoring plant and soil water status: established and novel methods revisited and their relevance to studies of drought tolerance. J. Exp. Bot. 58 (2), 119–130. doi: 10.1093/jxb/erl118
Li, X., Blackman, C. J., Choat, B., Duursma, R. A., Rymer, P. D., Medlyn, B. E., et al. (2018). Tree hydraulic traits are coordinated and strongly linked to climate-of-origin across a rainfall gradient. Plant Cell Environ. 41 (3), 646–660. doi: 10.1111/pce.13129
Li, H. J., Si, B. C., Wu, P. T., McDonnell, J. J. (2019). Water mining from the deep critical zone by apple trees growing on loess. Hydrol Process 33 (2), 320–327. doi: 10.1002/hyp.13346
Lu, Y. W, Si, B. C., Li, H. J., Biswas, A. (2019). Elucidating controls of the variability of deep soil bulk density. Geoderma 348. doi: 10.1016/j.geoderma.2019.04.033
Mackay, D. S., Savoy, P. R., Grossiord, C., Tai, X., Pleban, J. R., Wang, D. R., et al. (2020). Conifers depend on established roots during drought: Results from a coupled model of carbon allocation and hydraulics. New Phytol. 225 (2), 679–692. doi: 10.1111/nph.16043
Martin-StPaul, N., Delzon, S., Cochard, H. (2017). Plant resistance to drought depends on timely stomatal closure. Ecol. Lett. 20 (11), 1437–1447. doi: 10.1111/ele.12851
McDowell, N., Pockman, W. T., Allen, C. D., Breshears, D. D., Cobb, N., Kolb, T., et al. (2008). Mechanisms of plant survival and mortality during drought: Why do some plants survive while others succumb to drought? New Phytol. 178 (4), 719–739. doi: 10.1111/j.1469-8137.2008.02436.x
Moriana, A., Pérez-López, D., Prieto, M. H., Ramírez-Santa-Pau, M., Pérez-Rodriguez, J. M. (2012). Midday stem water potential as a useful tool for estimating irrigation requirements in olive trees. Agr Water Manage 112, 43–54. doi: 10.1016/j.agwat.2012.06.003
Moura, J., Bonine, C. A. V., Viana, J. D. F., Dornelas, M. C., Mazzafera, P. (2010). Abiotic and biotic stresses and changes in the lignin content and composition in plants. J. Integr. Plant Biol. 52 (4), 360–376. doi: 10.1111/j.1744-7909.2010.00892.x
Nadal-Sala, D., Grote, R., Birami, B., Knuver, T., Rehschuh, R., Schwarz, S., et al. (2021). Leaf shedding and non-stomatal limitations of photosynthesis mitigate hydraulic conductance losses in scots pine saplings during severe drought stress. Front. Plant Sci. 12, 715127. doi: 10.3389/fpls.2021.715127
Oliveira, R. S., Bezerra, L., Davidson, E. A., Pinto, F., Klink, C. A., Nepstad, D. C., et al. (2005). Deep root function in soil water dynamics in cerrado savannas of central Brazil. Funct. Ecol. 19 (4), 574–581. doi: 10.1111/j.1365-2435.2005.01003.x
Pammenter, N. W., Vander Willigen, C. (1998). A mathematical and statistical analysis of the curves illustrating vulnerability of xylem to cavitation. Tree Physiol. 18 (8-9), 589–593. doi: 10.1093/treephys/18.8-9.589
Pilar, B., Sanchez-De-Miguel, P., Centeno, A., Junquera, P., Linares, R., Lissarrague, J. R., et al. (2007). Water relations between leaf water potential, photosynthesis and agronomic vine response as a tool for establishing thresholds in irrigation scheduling. Sci. Hortic. 114 (3), 151–158. doi: 10.1016/j.scienta.2007.06.012
Puchalka, R., Dyderski, M. K., Vitkova, M., Sadlo, J., Klisz, M., Netsvetov, M., et al. (2021). Black locust (Robinia pseudoacacia l.) range contraction and expansion in Europe under changing climate. Glob Chang Biol. 27 (8), 1587–1600. doi: 10.1111/gcb.15486
Reichstein, M., Bahn, M., Ciais, P., Frank, D., Mahecha, M. D., Seneviratne, S. I., et al. (2013). Climate extremes and the carbon cycle. Nature 500 (7462), 287–295. doi: 10.1038/nature12350
Rodriguez-Dominguez, C. M., Brodribb, T. J. (2020). Declining root water transport drives stomatal closure in olive under moderate water stress. New Phytol. 225 (1), 126–134. doi: 10.1111/nph.16177
Rodriguez-Dominguez, C. M., Carins Murphy, M. R., Lucani, C., Brodribb, T. J. (2018). Mapping xylem failure in disparate organs of whole plants reveals extreme resistance in olive roots. New Phytol. 218 (3), 1025–1035. doi: 10.1111/nph.15079
Rodriguez-Zaccaro, F. D., Valdovinos-Ayala, J., Percolla, M. I., Venturas, M. D., Pratt, R. B., Jacobsen, A. L., et al. (2019). Wood structure and function change with maturity: Age of the vascular cambium is associated with xylem changes in current-year growthKey-words. Plant Cell Environ. 42 (6), 1816–1831. doi: 10.1111/pce.13528
Santesteban, L. G., Miranda, C., Marín, D., Sesma, B., Intrigliolo, D. S., Mirás-Avalos, J. M., et al. (2019). Discrimination ability of leaf and stem water potential at different times of the day through a meta-analysis in grapevine (Vitis vinifera l.). Agr Water Manage 221, 202–210. doi: 10.1016/j.agwat.2019.04.020
Scholander, P. F., Hammel, H. T., Bradstreet, E. D., Hemmingsen, E. A. (1965). Sap pressure in vascular plants - negative hydrostatic pressure can be measured in plants. Science 148 (3668), 339–33+. doi: 10.1126/science.148.3668.339
Schwarzel, K., Zhang, L. l., Montanarella, L., Wang, Y. H., Sun, G. (2019). How afforestation affects the water cycle in drylands: A process-based comparative analysis. Glob Chang Biol 26 (2), 944–959. doi: 10.1111/gcb.14875
Sperry, J. S., Adler, F. R., Campbell, G. S., Comstock, J. P. (1998). Limitation of plant water use by rhizosphere and xylem conductance: Results from a model. Plant Cell Environ. 21 (4), 347–359. doi: 10.1046/j.1365-3040.1998.00287.x
Steudle, E. (2000). Water uptake by roots: Effects of water deficit. J. Exp. Bot. 51 (350), 1531–1542. doi: 10.1093/jexbot/51.350.1531
Suter, B., Triolo, R., Pernet, D., Dai, Z., Van Leeuwen, C. (2019). Modeling stem water potential by separating the effects of soil water availability and climatic conditions on water status in grapevine (Vitis vinifera l.). Front. Plant Sci. 10, 1485. doi: 10.3389/fpls.2019.01485
Taccoen, A., Piedallu, C., Seynave, I., Perez, V., Gegout-Petit, A., Nageleisen, L. M., et al. (2019). Background mortality drivers of European tree species: Climate change matters. Proc. R. Soc. B-Biological Sci. 286 (1900), 10. doi: 10.1098/rspb.2019.0386
Tian, T., Schreiner, R. P. (2021). Appropriate time to measure leaf and stem water potential in north-south oriented, vertically shoot-positioned vineyards. Am. J. Enology Viticulture 72 (1), 64–72. doi: 10.5344/ajev.2020.20020
Trueba, S., Pan, R. H., Scoffoni, C., John, G. P., Davis, S. D., Sack, L., et al. (2019). Thresholds for leaf damage due to dehydration: Declines of hydraulic function, stomatal conductance and cellular integrity precede those for photochemistry. New Phytol. 223 (1), 134–149. doi: 10.1111/nph.15779
Tumajer, J., Scharnweber, T., Smiljanic, M., Wilmking, M. (2022). Limitation by vapour pressure deficit shapes different intra-annual growth patterns of diffuse- and ring-porous temperate broadleaves. New Phytol. 233 (6), 2429–2441. doi: 10.1111/nph.17952
Urli, M., Porte, A. J., Cochard, H., Guengant, Y., Burlett, R., Delzon, S., et al. (2013). Xylem embolism threshold for catastrophic hydraulic failure in angiosperm trees. Tree Physiol. 33 (7), 672–683. doi: 10.1093/treephys/tpt030
Venturas, M. D., Sperry, J. S., Love, D. M., Frehner, E. H., Allred, M. G., Wang, Y., et al. (2018). A stomatal control model based on optimization of carbon gain versus hydraulic risk predicts aspen sapling responses to drought. New Phytol. 220 (3), 836–850. doi: 10.1111/nph.15333
Walthert, L., Ganthaler, A., Mayr, S., Saurer, M., Waldner, P., Walser, M., et al. (2021). From the comfort zone to crown dieback: Sequence of physiological stress thresholds in mature European beech trees across progressive drought. Sci. Total Environ. 753 (991), 753. doi: 10.1016/j.scitotenv.2020.141792
Weithmann, G., Schuldt, B., Link, R. M., Heil, D., Hoeber, S., John, H., et al. (2021). Leaf trait modification in European beech trees in response to climatic and edaphic drought. Plant Biol. (Stuttg) 24 (7), 1272–1286. doi: 10.1111/plb.13366
Wu, M., Zhang, Y., Oya, T., Marcati, C.R., Pereira, L., Jansen, S., et al. (2020). Root xylem in three woody angiosperm species is not more vulnerable to embolism than stem xylem. Plant Soil 450 (1-2), 479–495. doi: 10.1007/s11104-020-04525-0
Yan, W., Zhong, Y., Shangguan, Z. (2017a). Rapid response of the carbon balance strategy in robinia pseudoacacia and amorpha fruticosa to recurrent drought. Environ. Exp. Bot. 138, 46–56. doi: 10.1016/j.envexpbot.2017.03.009
Yan, W., Zhong, Y., Shangguan, Z. (2017b). Responses of different physiological parameter thresholds to soil water availability in four plant species during prolonged drought. Agr For. Meteorol 247, 311–319. doi: 10.1016/j.agrformet.2017.08.017
Yuan, X., Wang, L., Wu, P., Ji, P., Sheffield, J., Zhang, M., et al. (2019). Anthropogenic shift towards higher risk of flash drought over China. Nat. Commun. 10 (1), 4661. doi: 10.1038/s41467-019-12692-7
Zhang, L., Merlin, I., Pascal, S., Bert, P. F., Domergue, F., Gambetta, G. A., et al. (2020). Drought activates MYB41 orthologs and induces suberization of grapevine fine roots. Plant Direct 4, 17. doi: 10.1002/pld3.278
Zhang, Y. J., Meinzer, F. C., Qi, J. H., Goldstein, G., Cao, K. F. (2013). Midday stomatal conductance is more related to stem rather than leaf water status in subtropical deciduous and evergreen broadleaf trees. Plant Cell Environ. 36 (1), 149–158. doi: 10.1111/j.1365-3040.2012.02563.x
Zhao, M., Running, S. W. (2010). Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. Science 329 (5994), 940–943. doi: 10.1126/science.1192666
Keywords: trees response to drought, xylem water potential, soil water availability, drought stress level, transpiration, stomatal conductance, xylem to rhizosphere conductance
Citation: Wang X, Fan Y, Zhang C, Zhao Y, Du G, Li M and Si B (2023) From comfort zone to mortality: Sequence of physiological stress thresholds in Robinia pseudoacacia seedlings during progressive drought. Front. Plant Sci. 14:1149760. doi: 10.3389/fpls.2023.1149760
Received: 23 January 2023; Accepted: 27 February 2023;
Published: 16 March 2023.
Edited by:
Santiago Signorelli, Universidad de la República, UruguayReviewed by:
Boris Lazarević, University of Zagreb, CroatiaSebastián Ricardo Simondi, Universidad Nacional de Cuyo, Argentina
Copyright © 2023 Wang, Fan, Zhang, Zhao, Du, Li and Si. 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: Bingcheng Si, QmluZy5TaUB1c2Fzay5jYQ==; Min Li, bGltaW4yMDE2QG53YWZ1LmVkdS5jbg==