Is vulnerability segmentation at the leaf-stem transition a drought resistance mechanism? A theoretical test with a trait-based model for Neotropical canopy tree species

Leaf-stem vulnerability segmentation predicts lower xylem embolism resistance in leaves than stem. However, although it has been intensively investigated these past decades, the extent to which vulnerability segmentation promotes drought resistance is not well understood. Based on a trait-based model, this study theoretically supports that vulnerability segmentation enhances shoot desiccation time across 18 Neotropical tree species. Leaf-stem vulnerability segmentation predicts lower xylem embolism resistance in leaves than stems thereby preserving expensive organs such as branches or the trunk. Although vulnerability segmentation has been intensively investigated these past decades to test its consistency across species, the extent to which vulnerability segmentation promotes drought resistance is not well understood. We investigated the theoretical impact of the degree of vulnerability segmentation on shoot desiccation time estimated with a simple trait-based model. We combined data from 18 tropical rainforest canopy tree species on embolism resistance of stem xylem (flow-centrifugation technique) and leaves (optical visualisation method). Measured water loss under minimum leaf and bark conductance, leaf and stem capacitance, and leaf-to-bark area ratio allowed us to calculate a theoretical shoot desiccation time (tcrit). Large degrees of vulnerability segmentation strongly enhanced the theoretical shoot desiccation time, suggesting vulnerability segmentation to be an efficient drought resistance mechanism for half of the studied species. The difference between leaf and bark area, rather than the minimum leaf and bark conductance, determined the drastic reduction of total transpiration by segmentation during severe drought. Our study strongly suggests that vulnerability segmentation is an important drought resistance mechanism that should be better taken into account when investigating plant drought resistance and modelling vegetation. We discuss future directions for improving model assumptions with empirical measures, such as changes in total shoot transpiration after leaf xylem embolism.


Introduction
The increasing frequency of extreme drought events due to climate change induces global tree mortality and forest dieback (Allen et al. 2010;Aleixo et al. 2019;Brodribb et al. 2020). Current climate predictions suggest drought and mortality phenomena to increase in space and time (Duffy et al. 2015), probably leading to significant impacts on forest species composition, their functioning, and productivity (Brienen et al. 2015(Brienen et al. , 2020Esquivel-Muelbert et al. 2019). This situation urges a better understanding of the mechanisms underlying drought-induced tree mortality, and also the possibility to quantify and predict tree drought resistance (McDowell et al. 2018;Brodribb et al. 2020). The field of plant hydraulics can contribute to this discussion because water transport and storage in plants are tightly coupled with photosynthesis and growth. In particular, hydraulic failure of the water transport system due to high water stress is now recognised as one of main causes of drought-induced tree mortality (Urli et al. 2013;Martin-StPaul et al. 2017;Blackman et al. 2019b;Brodribb et al. 2020). Other hydraulic traits related to water absorption and rooting depth, transpiration control, and water storage have been identified in shaping drought resistance of plants (Blackman et al. 2016;Martin-StPaul et al. 2017;Brum et al. 2019). However, we still lack a clear understanding of how various drought resistance traits determine drought resistance in a complementary way.

Quantifying plant drought resistance
Plant drought resistance is usually defined as the desiccation time between the onset of water stress with stomatal closure, various stages of tissue dehydration, and eventually mortality (Blackman et al. 2016;Volaire 2018). This definition emphasises on the maximisation of survival during one particular water shortage, and can be simply quantified as the time between the onset of a water stress and plant death by desiccation, i.e. desiccation time. All the traits enabling the expansion of this desiccation time can be considered as drought resistance mechanisms. The plant water status is controlled by stomata, and stomatal closure is generally the first plant response to water stress (Martin-StPaul et al. 2017). Comparison of stomatal closure point across species is allowed by measuring the leaf water potential inducing ca. 90-95% loss of maximum stomatal conductance (Ψ closure ; MPa). The earlier the stomatal closure, the slower the reduction in xylem water potential (Martin-StPaul et al. 2017). An excessive water potential drop of the xylem can cause embolism (i.e. bubble formation within conduits). Embolism blocks water transport through conduits, and accumulation of embolism drives the loss of xylem hydraulic conductivity (Tyree and Sperry 1989), finally leading to plant death by desiccation (Urli et al. 2013). Drought-induced xylem embolism resistance is quantified by measuring the xylem pressure at which a certain percentage of xylem hydraulic conductivity is lost. For comparison across plant species or individuals, the xylem water potential inducing 50% loss of xylem hydraulic conductivity (P 50 ; MPa) is generally used as a reference (Choat et al. 2012). Desiccation time is theoretically maximised by more negative P 50 values, and a large difference between Ψ closure and P 50 , i.e. the stomatalhydraulic safety margin (SHSM; MPa). Earlier hydraulic studies have largely focused on xylem embolism resistance, and allowed to demonstrate that P 50 is tightly related to drought mortality and drought resistance in controlled experiments (Resco et al. 2009;Brodribb and Cochard 2009;Urli et al. 2013;Blackman et al. 2019a).
The desiccation time between Ψ closure and P 50 is influenced by water pools stored in plant tissues, and by residual water losses (through leaves and bark) which determine the rate of depletion of these plant water reservoirs. These processes are analogical to a source-sink model. Indeed, when stomata close, leaf transpiration is not null. The leaf has a certain permeability, and closed stomata could still lose water via transpiration, which determines minimum leaf conductance (g min ; mmol m −2 s −1 ; Duursma et al. 2019). The bark is equally prone to transpiration, because of an incomplete permeability of the different tissue layers that constitute the bark (Roth 1981), or lenticels, which affect drought resistance (Wolfe 2020). However, quantification of bark conductance (g bark ; mmol m −2 s −1 ) across species remains scarce, particularly in the perspective of drought resistance (Blackman et al., 2019b; but see Wolfe, 2020). Overall, the lower the minimum conductance of the bark and leaves, the longer the desiccation time. The plant water reservoirs, on the other hand, are determined by tissue capacitance (W; kg m −3 MPa −1 ), which measures the capacity of a tissue to release water per unit of tissue volume and per unit of pressure drop (Tyree and Ewers 1991). The capacitance is determined by the empty space per unit of tissue volume, thus W is strongly related to tissue density (Scholz et al. 2011;Jupa et al. 2016;Santiago et al. 2018). Since the product of g min and total leaf area determines the whole-plant minimum transpiration (E min ; mmol s −1 ) after stomatal closure, and the product of W and the volume of wood or stem determines the total capacity (C; kg MPa −1 ), it can be hypothesised that the ratio between total leaf area and stem volume could have a strong functional significance regarding desiccation time.
Desiccation time can be experimentally quantified as the time between the onset of water stress and death of an entire plant or a plant organ (e.g. a branch, a leaf). This approach is relatively easy to implement in controlled experiments such as in greenhouses, where desiccation time has been demonstrated to be high for drought-resistant species (Urli et al., 2013;Li et al., 2016;Blackman et al., 2019a,b). However, observation of desiccation time is hardly feasible in natural conditions, especially for trees. To better predict desiccation time in the field, Blackman et al. (2016) have recently developed a trait-based model, integrating the drought resistance traits described above (P 50 , Ψ closure , g min , W). The model calculates desiccation time (t crit ; hrs) as the ratio between (1) the total amount of water releasable by tissues across the water potential drop from Ψ closure to critical levels of embolism (P crit ; MPa) and (2) the rate of transpiration after stomatal closure, such as: where t crit (h) is the desiccation time, Ψ closure is the leaf water potential at stomatal closure (MPa), W is the tissue capacitance (kg m −3 MPa −1 ), g min is the minimum leaf conductance (mmol m −2 s −1 ), M H20 is the water molar mass of 18.015 × 10 −6 kg mmol −1 , V plant (m 3 ) is the volume of plant tissue releasing water, A L (m 2 ) is the total leaf area, and VPD is the vapour pressure deficit (mol mol −1 ). The factor 3600 converts seconds in hours. The total leaf area is multiplied by two because leaf minimum transpiration-determined by g min -occurs at both abaxial and adaxial leaf sides. Blackman et al. (2019b) found a strong relationship between observed and calculated t crit across eight Eucalyptus species, in an experiment based on potted individuals, strengthening this trait-based model approach. This framework is promising for exploring the relative importance of mechanisms behind drought resistance, as well as the diversity of drought resistance strategies.

Vulnerability segmentation and drought resistance
Vulnerability segmentation has long been hypothesised to have a role in plant drought resistance (Tyree and Zimmermann 2002). The vulnerability segmentation hypothesis (VSH) predicts distal organs such as leaves and roots to be more vulnerable to xylem embolism relative to perennial organs such as the trunk or branches of a tree. One implicit implication of the preservation of perennial organs is that survival time is enhanced. The a posteriori metaphor of "safety valves" or "safety fuse" has more explicitly presented the VSH as a potential mechanism related to drought resistance (Chen et al. 2009(Chen et al. , 2010. At the leaf-stem transition, the transpiration flow is directly proportional to leaf area. Then, a decrease in the evaporative leaf area following the rupture of the leaf water column through drought-induced leaf xylem embolism can alleviate the water potential drop in the stem during drought. In other words, the relative contribution of minimum leaf transpiration decreases with the accumulation of leaf xylem embolism, thereby strengthening the relative contribution of the bark to total transpiration. Therefore, vulnerability segmentation should enhance the time from leaf P crit to stem P crit . In the last 20 years, the VSH has gained renewed interest and great efforts have been made to demonstrate the existence or the lack of vulnerability segmentation across species (Johnson et al. 2016;Hochberg et al. 2016;Charrier et al. 2016;Skelton et al. 2017Skelton et al. , 2018; Rodriguez-Dominguez Klepsch et al. 2018;Levionnois et al. 2020;Smith-Martin et al. 2020;Li et al. 2020). However, even if it is clear that vulnerability segmentation at the leaf-stem transition can be common, and is species dependant, the functional significance of vulnerability segmentation in the perspective of drought resistance remains unclear and poorly investigated. Recently, Blackman et al. (2019b) showed that leaf-stem vulnerability segmentation explains t crit for Eucalyptus species. Such approach needs to be replicated to identify to what extent we need to recognise vulnerability segmentation as an important drought resistance mechanism in most biomes and lineages across the world, which would enable better predictions of tree and vegetation dynamics under climate change (Martin-StPaul et al. 2017;Brodribb et al. 2020).
The strength of the effect of vulnerability segmentation on desiccation time can theoretically be affected by several factors. For instance, vulnerability segmentation implies that leaf xylem embolism reduces the total leaf evaporative surface. But the difference between leaf and bark area (or the leaf area-to-bark area ratio) at a given level (e.g. the shoot, the branch, the tree) can vary across species, potentially affecting the functional significance of vulnerability segmentation on desiccation time across species. Moreover, to our knowledge, a comparison between g min and g bark has not been analysed yet. It can be hypothesised that the minimum conductance is different between leaf and bark, further affecting desiccation time after vulnerability segmentation across species. Finally, the ratio between the total water content and the evaporative surface is another parameter behind vulnerability segmentation, potentially varying across species and then affecting desiccation time. However, these parameters have never been investigated in the perspective of vulnerability segmentation controlling for desiccation time across species.
Here, we aim to test the theoretical effect of leaf-stem vulnerability segmentation on desiccation time of shoots across Neotropical rainforest canopy tree species. Rainforest species have previously been suggested to lack vulnerability segmentation (Zhu et al. 2016), but Levionnois et al. (2020) recently showed that vulnerability segmentation was common across nearly half of the tropical tree species they studied. Based on data from Levionnois et al. (2020), and by incorporating novel measurements, we calculated desiccation time at the shoot level, with a trait-based model taking or not vulnerability segmentation into account. We also added bark conductance and surface-area measurements to fully characterise the theoretical effect of vulnerability segmentation. We specifically addressed the following questions: • Does leaf-stem vulnerability segmentation have a theoretical effect on desiccation time of tropical tree shoots? After answering this question, we also discuss future directions to improve estimations of desiccation time to better take into account vulnerability segmentation. • What are the mechanisms (i.e. P crit,leaf -P crit,stem ; g ming bark ; A leaf -A bark ) allowing to expand shoot desiccation time through vulnerability segmentation?

Study site, species, and design
The experiment was conducted in French Guiana at the Paracou experimental station (https:// parac ou. cirad. fr/ websi te; 5° 16′ 26″ N, 52° 55′ 26″ W), which represents a lowland tropical rainforest (Gourlet-Fleury et al. 2004). The warm and wet tropical climate of French Guiana is highly seasonal due to the North-South movement of the Intertropical Convergence Zone. Mean (± SE) annual air temperature is 25.7 °C ± 0.1 °C and the mean annual precipitation is 3,102 mm ± 70 mm (data between 2004 and 2014; Aguilos et al. 2019). There is a dry season lasting from mid-August to mid-November, during which rainfall is < 100 mm month −1 . We sampled only dominant canopy trees, and only growing on terra firme. A total of 50 trees belonging to 18 species and 50 trees were sampled, with three trees per species for 14 species, and two trees per species for four species (Bocoa prouacensis, Iryanthera sagotiana, Moronobea coccinea, Qualea rosea). These 50 trees represent a subset of the sampling in Ziegler et al. (2019), for which stem embolism resistance was initially measured. The same dataset was used by Levionnois et al. (2020), who added data on leaf embolism resistance. Measurements of stem or leaf embolism resistance were particularly challenging for species for which data on only two individuals were available. Although several samples have been collected for a given species, only two vulnerability curves (VC) per species could be obtained (Ziegler et al. 2019;Levionnois et al. 2020). The studied species covered a broad phylogenetic diversity such that the main clades of the flowering plants were represented, i.e. magnoliids, rosids, and asterids. Branches were sampled by professional tree climbers.
We calculated desiccation time at the shoot level, with a shoot defined in the context of this study as an unbranched and leafy stem, with the largest possible number of leaves and with all nodes bearing a leaf. The traits necessary for calculating shoot desiccation time were measured as much as possible at the shoot level, except for stem embolism resistance and stem capacitance, which were measured further upstream, on older branch sections (~2 cm in diameter).

Stem embolism resistance
Stem embolism resistance was measured as the stem water potential inducing 50% and 88% loss of conductivity (P 50,stem and P 88,stem , respectively; MPa; Table 1). Sample collection was held between January and July 2017, with one 3-m-long canopy branch sampled per tree. Only sun-exposed branches were sampled, except for the understory species Gustavia hexapetala. Vulnerability curves were measured with the flow-centrifugation technique using a Cavi1000 system (DGMeca, BIOGECO lab, Gradignan, France). A detailed description of the methods and the data for each species were provided and discussed by Ziegler et al. (2019).

Leaf embolism resistance
Leaf embolism resistance was measured as the leaf water potential inducing 50% and 88% loss of conductivity, respectively (P 50,leaf and P 88,leaf respectively; MPa; Table 1). The field sampling was held between November 2018 and March 2019, on the same trees that were sampled for the measurements of stem embolism resistance. We generally sampled three trees per sampling day, during the morning and before solar midday, in order to avoid too negative leaf water potentials and potential native leaf embolism. Only sun-exposed branches were sampled as much as possible. One 1-m-long canopy branch was sampled per tree, with ~50 leaves or leaflets for the monitoring of water potentials, as described below. A second small branch was also sampled for bark conductance and morphological measurements, as described below. To measure leaf xylem embolism resistance, we relied on the optical light transmission method (Brodribb et al., 2016a(Brodribb et al., , 2016b; https:// www. opens ource ov. org). For details about the method and data for each species of this study, we refer the reader to Levionnois et al. (2020).

Leaf pressure-volume curves
Stomatal closure is generally measured from the stomatal response curve to leaf water potential based on stomatal conductance measurements. However, Ψ closure and leaf water potential at leaf turgor loss point (Ψ tlp ; MPa; Table 1) are strongly related (Brodribb et al. 2003;Rodriguez-Dominguez et al. 2016;Bartlett et al. 2016), leading to the possibility to use Ψ tlp as a surrogate of Ψ closure (Martin-StPaul et al. 2017). Ψ tlp was measured by pressure-volume analysis using the bench drying technique (Koide et al., 1989). Pressure-volume measurements were based on the measurement of one leaf per tree. Before the branch used for the measurement of leaf embolism resistance was removed from water and plastic bags, one fully hydrated and healthy leaf was sampled. At regular intervals, leaf water potential was measured with a Scholander pressure chamber (Model 600, PMS), directly followed by the measurement of leaf mass with a precision balance (AB 204-S, Mettler Toledo). The different water potential and mass measurements were made to cover the range from 0.00 to −3.00 MPa (or more if necessary as the dynamic of water potential drop is generally clearly noticeable during measurement), with ca. 0.30 MPa intervals. Then, the leaf was scanned to measure leaf area, and dried for 3 days at 70 °C to measure leaf dry mass Ψ tlp being the leaf water potential at the phase transition (Koide et al. 1989). The pressure-volume curve was also used to derive leaf capacitance after turgor loss W tlp,leaf (mmol m −2 MPa −1 ; Table 1), since the plant water stock used for the calculation of t crit is the water available in plant tissues after stomatal closure. W tlp,leaf was calculated as the slope of the loss of leaf water mass according to leaf water potential for water potentials below Ψ tlp , and divided by the leaf area.
We used the estimation of wood capacitance as a proxy for whole stem capacitance. We are aware that the different stem tissues (pith, xylem, phloem, cortex) can exhibit different capacitances, but we assume capacitances between stem tissues to be coordinated across species. During field work for the measurements of stem embolism resistance, we sampled short stem segments (~1.5 cm in length; ~1 cm in diameter) for WD measurements. One wood sample was sampled per branch, on the distal extremity of the sample used for the Cavi1000 measurement. The fresh volume of the sample was calculated using an inverse Archimedes principle and a precision balance (CP224S, Sartorius), as the difference between fresh mass and immersed mass of the sample. Dry mass was derived after drying at 103 °C for 3 days. WD was calculated as the ratio between dry mass and fresh volume. In this present study, WD varied from 0.49 to 0.85, whereas WD in the study of Santiago et al. (2018) varied from 0.33 to 0.80.

Leaf and bark minimum conductance
Leaf minimum conductance g min (mmol m −2 MPa −1 ; Table 1) was measured between July and November 2017, with the method of weight loss of detached leaves (Blackman et al. 2019b;Duursma et al. 2019). During field work for the measurements of stem embolism resistance, while sampling the branches, strings were tied to sunlit shoots to facilitate later canopy sampling without climbing. g min was Stem capacity to bark surface W ft,stem V stem /A bark measured on one leaf per tree. Immediately after leaf sampling, leaves were placed in small sealed plastic bags with wet paper towels to avoid dehydration. At the laboratory, the entire leaf surface was gently padded dry with dry paper towels. The cut petiole was sealed with nail polish such that the total leaf transpiration was only mediated by the lamina. Leaves were placed 1 m above air fans and let to dehydrate in a closed air-conditioned room. Air temperature and relative humidity were maintained at ca. 25 °C and ca. 50%, respectively. At regular intervals, leaf mass was recorded with a precision balance (AB 204-S, Mettler Toledo). Time, air temperature, and relative air humidity were noted for each measurement. g min was calculated as the slope of decreasing leaf mass according to time, in conjunction with VPD-calculated from the known air temperature and relative humidity and divided by the doubled leaf area. Particular attention was paid to the kinetics of water loss, by estimating g min only during the solely linear phase of leaf water mass decrease with time after stomata had closed. The bark conductance g bark (mmol m −2 MPa −1 ; Table 1) was measured during the measurement of leaf embolism resistance, for which a second small branch was sampled. From this branch, three 20-cm-long and ~1-cm-diameter stem segments were sampled, by paying attention to have linear and unbranched stem segments. Then, the segments were placed in small sealed plastic bags with wet paper towels to avoid dehydration and eventual bark shrinkage. At the laboratory, the entire bark surface was dried with paper towels. The two extremities of each stem segment were sealed with melted wax such that the total stem segment transpiration was only mediated by the bark. g bark was measured exactly in the same manner as described for leaves, but the slope of the decreasing water mass of the stem segment over time was divided by the bark surface area. Particular attention was paid to the kinetics of water loss only during the solely linear phase of leaf water mass decrease with time. This also allowed to avoid biases due to the potential release of water that would be just recently absorbed by the bark. The bark surface area was measured by considering the stem segment as a cylinder, based on measurements of stem length and diameter.

Morphological parameters
Shoot leaf area (A leaf ; m 2 ) and shoot bark area (A bark ; m 2 ) were quantified during leaf embolism resistance measurements, with a second small branch used for bark conductance measurements only. On this branch, a shoot was sampled, defined as an unbranched stem supporting as many leaves as possible, but with no missing or damaged leaves between the youngest and the oldest leaf. Once sampled from the tree, shoots were placed in sealed plastic bags with wet paper towels to avoid dehydration. At the laboratory, all the surfaces were dried with paper towels. All the leaves were cut and scanned to measure total shoot leaf area (Table 1). The shoot bark area was measured by considering the shoot stem as a cylinder, from length and diameter measurements ( Table 1). The stem volume (V stem ; m 3 ) was calculated as a cylinder.

The desiccation time model
To test the effect of vulnerability segmentation on t crit , we calculated t crit for species exhibiting a positive vulnerability segmentation (i.e. P 88,leaf > P 88,stem ; Levionnois et al., 2020), representing 12 species. Then, we calculated for the same species a theoretical t crit with no segmentation, such as P 88,leaf = P 88,stem . To take into account the vulnerability segmentation in our calculation of t crit , we used the framework of Blackman et al. (2016) and refined the model. We added two successive desiccation times: (i) leaf desiccation time, from stomatal closure to leaf xylem hydraulic failure, and (ii) stem desiccation time, from leaf xylem to stem xylem hydraulic failure. The model assumes that during phase (i), shoot transpiration is driven by both g min and g bark ; but that during phase (ii), shoot transpiration is driven by g bark only. Contrary to Blackman et al. (2016) we also added the effect of the leaf capacitance after leaf turgor loss point, but before leaf xylem hydraulic failure. We used P 88,leaf and P 88,stem as parameters for critical leaf and stem hydraulic failure, respectively (Blackman et al. 2016), as P 88,stem has been shown to trigger plant mortality for angiosperm species (Urli et al. 2013). Therefore, t crit_seg for species exhibiting a positive vulnerability segmentation was calculated as: where Ψ tlp is the water potential (MPa) at turgor loss point used as a surrogate of stomatal closure; P 88,leaf is the leaf water potential (MPa) at 88% loss of leaf xylem conductivity; P 88,stem is the water potential (MPa) at 88% loss of stem xylem conductivity; W ft,stem is the stem capacitance (kg m −3 MPa −1 , but divided by the water molar mass of 18.015 × 10 −6 kg mmol −1 ) at full turgor; W tlp,leaf is the leaf capacitance (mmol m −2 MPa −1 ) after leaf turgor loss point; V stem is the stem volume (m 3 ); A leaf is the shoot leaf area (m 2 ); g min is the leaf minimum conductance (mmol m −2 s −1 ); g bark is the bark conductance (mmol m −2 s −1 ); and VPD is the vapour pressure deficit (mol mol −1 ). The factor 3,600 converts seconds in hours. Then, for species exhibiting a positive vulnerability segmentation, we artificially t crit_seg = Ψ tlp − P 88,leaf * W ft,stem V stem + W tlp,leaf A leaf 2A leaf g min + A bark g bark * 3600 * VPD + P 88,leaf − P 88,stem * W ft,stem V stem A bark g bark * 3600 * VPD cancelled the segmentation by calculating a t crit_noseg , with P 88,leaf = P 88,stem , and therefore: All the parameters of the models and per species are presented in the Table 2. To better characterise how vulnerability segmentation expands t crit , we calculated additional traits at the shoot and stem levels by combining couples of aforementioned traits ( Table 1) that we compared across shoot and stem levels. We calculated leaf minimum transpiration rate (g min *2A leaf ; mmol s −1 ) as the product of g min and shoot leaf area multiplied by 2 (because transpiration occurs at both abaxial and adaxial sides). We calculated bark transpiration rate (g bark *A bark ; mmol s −1 ) at the shoot level. For the A bark considered, we come back to our definition of the leafy shoot, as a single unbranched stem supporting the most leaves, with no missing or damaged leaves between the youngest and the oldest leaf. Then, A bark was directly related to the number of nodes and leaves that constitute the leafy shoot. We calculated shoot capacity to total evaporative surface (C shoot A leaf+bark −1 ; mmol MPa −1 m −2 ), by multiplying A leaf by 2. We calculated stem capacity to bark surface (C stem A bark −1 ; mmol MPa −1 m −2 ). Capacity to evaporative surface translates the available stored water which will evaporate during dehydration as a function of total evaporative surface, with a higher capacity to evaporative surface amplifying the positive effect of stored water on drought resistance.

Statistical analyses
All statistical analyses were performed in R (R Core Team 2018). Data were tested for normality (Shapiro-Wilk test; α = 0.05). For correlations between traits, we used Pearson or Spearman correlation analyses depending on the normality. Comparison tests were conducted with Student's, Welch's, or Mann-Whitney-Wilcoxon's test, depending on the parameters of the samples (sample size, normality of distribution, variance).
To test for the effect of a given trait on t crit , it would be appropriate to apply a variance decomposition, which requires a multiple linear model based on the sum of explanatory variables. However, the t crit model is rather based on products. Switching from a product-based to a sum-based model is feasible with a log transformation. But in this case, this rather leads to a more complex model where the application of variance decomposition would lose its relevance. Therefore, as R 2 in bivariate relationships conveys information on explained variance, we calculated R 2 for relationships between t crit (response variable) and a given constitutive trait of the model (predictive variable). This was realised t crit_noseg = Ψ tlp − P 88,stem * W ft,stem V stem + W tlp,leaf A leaf 2A leaf g min + A bark g bark * 3600 * VPD for all traits constitutive of the t crit model and ordered the relationships according to R 2 . Then, we assumed that traits with higher R 2 have the strongest relative contribution to t crit .

Results
For the 12 species exhibiting a positive vulnerability segmentation (P 88,leaf > P 88,stem ), there was a significant difference between t crit_seg and t crit_noseg (Fig. 1), with a mean of 235.9 ± 71.3 h (±SE) and 9.0 ± 2.1 h, respectively. When considering all species, including species with both positive (P 88,leaf > P 88,stem ) and negative (P 88,leaf < P 88,stem ) vulnerability segmentation, t crit was positively related to the degree of P 88 -segmentation (Fig. 2a) as well as the degree of P 50 -segmentation (Fig. 2b), such that larger positive segmentation degree was associated with longer shoot desiccation time. t crit × VPD was also positively related to the stomatal-hydraulic safety margin (SHSM; MPa) when considering both P 88,stem (Fig. 2c) and P 50,stem (Fig. 2d), such that larger safety margin was associated with longer shoot desiccation time.
From all the constitutive traits of the t crit model, P 88 -segmentation and SHSM were the traits with the highest R 2 , both when considering species altogether or when considering only species exhibiting a positive vulnerability segmentation (P 88,leaf > P 88,stem ) ( Table 3). This suggests that P 88 -segmentation and SHSM had the highest relative contribution to t crit .
For the 12 species exhibiting a positive vulnerability segmentation (P 88,leaf > P 88,stem ), there was a significant difference between P 88,leaf and P 88,stem across species (Fig. 3a), with a mean of −4.10 ± 0.29 MPa and −6.08 ± 0.62 MPa, respectively. There was no significant difference between g min and g bark (Fig. 3b), with a mean of 5.529 ± 0.993 mmol m −2 s −1 and 4.285 ± 0.386 mmol m −2 s −1 , respectively. There was a significant difference in evaporative surfaces, between 2*A leaf and A bark (Fig. 3c), with a mean of 780.0 ± 171.6 cm 2 and 13.40 ± 1.91 cm 2 , respectively. A significant difference was also found between minimum leaf minimum transpiration (g min *2A leaf ) and bark transpiration (g bark *A bark ; Fig. 3d), with a mean of 0.513 ± 0.251 mmol s −1 and 0.006 ± 0.001 mmol s −1 , respectively. There was also a significant difference between shoot capacity per evaporative surface (C shoot A leaf+bark −1 ) and stem capacity per bark area (C stem A bark −1 ; Fig. 3e), with a mean of 640.8 ± 78.1 mmol MPa −1 m −2 and 15.226 ± 2.406 mmol MPa −1 m −2 , respectively, such that the total capacity per evaporative surface was higher following leaf embolism.
For the 12 species exhibiting a positive vulnerability segmentation (P 88,leaf > P 88,stem ), t crit_seg was not related to the difference between g min and g bark (Fig. 4a), nor to the  (Fig. 4b), nor to the difference between g min *2A leaf and g bark *A bark (Fig. 4c). t crit_seg was also not related to the difference between C shoot A leaf+bark −1 and C stem A bark −1 (Fig. 4d), even if a tendency exists at P < 0.1.

Discussion
Our study demonstrates that leaf-stem vulnerability segmentation theoretically enhances shoot desiccation time, supporting that vulnerability segmentation is an efficient drought resistance mechanism in tropical rainforest tree species. The effect of vulnerability segmentation on t crit is mainly due to the difference in leaf vs. bark area at the shoot level, but not to any difference in area-based leaf vs. bark conductance. The results are based on theoretical estimations and the model assumes that no water flows from stem to leaf xylem after embolism has occurred in leaf xylem. Further empirical tests of this assumption are needed, but we argue that not incorporating vulnerability segmentation can result in a considerable underestimation of drought resistance.

Leaf-stem vulnerability segmentation is a drought resistance mechanism
Not accounting for vulnerability segmentation into our desiccation time model for species with a positive vulnerability segmentation (P 88,leaf > P 88,stem ) dramatically underestimated shoot desiccation time, which lost one order of magnitude, shifting from 235.9 ± 71.3 h to 9.0 ± 2.1 h (Fig. 1). Moreover, the segmentation degree was the strongest predictor of t crit (Fig. 2a; Table 3). This finding highlights that leaf-stem vulnerability segmentation enhances-at least theoretically-shoot desiccation time, and is potentially a strong determinant of drought resistance at the leafy shoot level. This finding is in agreement with Blackman et al. (2019b), who found a relationship between the leaf-stem segmentation degree and observed t crit across eight Eucalyptus species in a controlled experiment. Furthermore, the degree of vulnerability segmentation seems to be as strong as the stomatal-hydraulic safety margin in determining t crit (SHSM; Ψ tlp -P crit,stem ; Table 3) in the context of the present study, bearing in mind that P crit,stem can vary across species (Li et al. 2016). This finding is also in agreement with Blackman et al. (2019b), who reported a higher coefficient of determination for P 50,leaf -P 50,stem than SHSM, when predicting their observed t crit . SHSM is an important drought resistance parameter, as it defines the complete sequence of drought stress, from stomatal closure to plant death, reflecting the water potential drop, the plant is able to sustain before die-back is likely. SHSM has already been shown to directly determine desiccation time in other studies, whether based on the integrative SurEau model (Martin-StPaul et al. 2017), or based on observations of desiccation time in controlled conditions (Blackman et al. 2019a). All these elements suggest vulnerability segmentation to be a drought resistance mechanism comparable in effectiveness to SHSM. Consequently, future trait-based assessment of plant drought resistance should incorporate measurements of both P 50,leaf and P 50,stem . Positive leaf-stem vulnerability segmentation is species dependant and common (Johnson et al. 2016;Hochberg et al. 2016;Charrier et al. 2016;Skelton et al. 2017Skelton et al. , 2018Rodriguez-Dominguez et al. 2018;Klepsch et al. 2018;Levionnois et al. 2020;Smith-Martin et al. 2020;Li et al. 2020). Then, omitting P 50,leaf and assuming no vulnerability segmentation for each species can lead to a strong underestimation of plant drought resistance. Such underestimation can be an issue if hydraulic trait data are used in vegetation prediction models (Brodribb et al. 2020). However, in a long-term perspective, more data are needed to fully characterise vulnerability segmentation, desiccation time (measured or estimated), and drought resistance together. For instance, the way we consider drought resistance is a maximisation of survival during one particular water shortage, such as an exceptional long dry season. However, maximising of desiccation time through vulnerability segmentation comes at the cost of the functional crown, and carbon uptake. Comparison for shoot desiccation time between the 12 positively segmented species (P 88,leaf > P 88,stem ) and the same species but artificially not segmented (P 88,leaf is set as same value as P 88,stem for each species). All abbreviations are presented in Table 1 Repeatedly relying on vulnerability segmentation as a consequence of exceptional drought events could deplete carbon reserves, and finally impede growth. Overall, the resilience of trees and forests under repetitive periods of drought stress is largely unknown (DeSoto et al. 2020).

Improving desiccation time estimations
Future directions should also alleviate the assumptions on which the model is based. First, our model assumes that following leaf xylem embolism, there is no longer water flow from the stem to the leaf. This assumption directly and simply defines the abrupt loss of total evaporative surface following leaf embolism, shifting from shoot leaf area plus bark area, to bark area alone. This assumption directly determines the effectiveness of the model to expand t crit for positively segmented species. Therefore, if this assumption is invalid, t crit may be critically overestimated. For instance, Wolfe et al. (2016) have shown for one out of six species that leaf shedding was not associated with a protection against a water potential drop in a tropical dry forest.
To our knowledge, no study has monitored the water mass loss of a shoot or a branch preceding and following leaf xylem embolism or leaf shedding. It is probable that some remaining water losses occur at the leaf level, even if leaf minimum transpiration is strongly reduced after leaf embolism. It would be beneficial if future studies would monitor total transpiration flow of whole shoots or branches, over desiccation time, to empirically assess if the hypothesised mechanism behind segmentation (reduction of the evaporative surface due to leaf embolism) is true.
Second, stem desiccation may start before leaf hydraulic failure, which may impact the estimation of t crit . For instance, a recent study has suggested that pressure-driven embolism spreading is not the exclusive mechanism of embolism spreading, but that pre-existing embolism may also be determinant (Guan et al. 2020). Such pre-existing embolism has been shown to be more likely in large and long vessels characterising stems rather than leaves.
Third, the loss of evaporative surface, either driven by leaf embolism or leaf shedding, can be progressive and gradual in space and time, rather than sudden as assumed Fig. 2 Relationships between shoot desiccation time and vulnerability segmentation for all species. a Shoot desiccation time according to P 88 -segmentation (P 88,leaf -P 88,stem ). b Shoot desiccation time according to P 50 -segmentation (P 50,leaf -P 50,stem ). c Shoot desiccation time according to the stomatalhydraulic safety margin (Ψ tlp -P 88,stem ). d Shoot desiccation time according to the stomatalhydraulic safety margin (Ψ tlp -P 50,stem ). All abbreviations are presented in Table 1 Annals of Forest Science (2021) 78: 87 by the model. This has been shown for leaf shedding by Blackman et al. (2019b), and for leaf shedding mediated by vulnerability segmentation by Hochberg et al. (2017). In case of progressive leaf shedding mediated by vulnerability segmentation of grapevine (Hochberg et al. 2017), it implies that P crit,leaf exhibits a certain degree of variation within a species, a crown, or a shoot, which has been properly demonstrated for Persea americana (Cardoso et al. 2020). However, our current model assumes a homogenous P crit,leaf across leaves. Not taking into account such progressive loss of evaporative surface in our model may overestimate t crit . Therefore, future investigations should (i) try to quantify intra-branch or intra-crown P 50,leaf variability, and (ii) to monitor the kinetics of leaf embolism or shedding over time, and its impact on desiccation time across species.

How does leaf-stem vulnerability segmentation expand shoot desiccation time?
Positively segmented species are species for which leaf xylem is more vulnerable to embolism than stem xylem ( Fig. 3a), which increases shoot desiccation time (Fig. 2s,  b). Across those species, the shoot leaf area is 30-fold higher than the stem bark area (Fig. 3c), but there is no difference between leaf minimum conductance and bark conductance (Fig. 3b). At the shoot level, the difference between leaf and bark surface area is sufficient to drive a cross-species difference in transpiration rate between the leaf compartment and the bark compartment (Fig. 3d). These findings suggest that the effect of vulnerability segmentation on t crit is largely driven by the difference in evaporative surface between leaf and bark, and not by any difference in conductance. Bark conductance is not particularly low, as it is not different from g min , with a mean of 4.29 mmol m −2 s −1 . This is comparable to the g min mean of 4.90 mmol m −2 s −1 found by Duursma et al. (2019) across 221 species based on a meta-analysis. This is also in agreement with g bark measurements of Wolfe (2020) across eight tropical tree species (0.86-12.98 mmol m −2 s −1 ). On the one hand, at the shoot level, bark surface represents a small fraction of the total evaporative surface. However, for larger levels (i.e. branch, trunk), the leaf area-to-bark area ratio may decrease, as the leaf area to stem volume ratio generally decreases with tree size (Scholz et al. 2011). At these levels, bark area in combination with g bark may drive a larger fraction of the total plant transpiration, further playing a more important role on t crit . On the other hand, g bark may be prone to dramatic ontogenetic effects by secondary growth (Lendzian 2006;Paine et al. 2010). In the present study, g bark measurements were based on relatively young stems, of ca. 0.5-to 1-cm diameter. At this stage, the bark is relatively thin, even green and poorly lignified for some species, and may be a photosynthetic tissue (Santiago et al. 2016). As the tree grows in diameter, the cortex and the periderm become thicker. Therefore, for large trees, lower g bark may alleviate the effect of bark transpiration on t crit despite a larger bark surface. Future work should incorporate g bark measurements, potential ontogenetic effects, and assess the sensitivity of the model to these parameters (g bark ; A leaf -to-A bark ratio) when scaling up, from shoot to branch to tree levels.
The difference between shoot and stem capacity per evaporative surface is not related to t crit (Fig. 4d), but there is a tendency of positive relationship if we consider P < 0.1. This tendency along with the dramatic shift following leaf embolism (Fig. 3e) suggests that the capacity per unit of evaporative surface is a strong determinant of t crit at the shoot level. However, this may be partly explained by the fact that the capacity per evaporative surface is explicit in the t crit formula. At the shoot level, stem volume and thus the amount stored water are relatively small. It is well known that stem volume and total stem water content to leaf area ratio increase with plant size (Scholz et al. 2011). Therefore, the role of capacitance and water reservoirs in controlling t crit should be more important for adult canopy trees for which wood volume represents several cubic meters. Moreover, here we considered the shoot as hydraulically independent and separated from the rest of the tree, whereas capacitive effects of upstream organs such as branches and the trunk may have a role in supplying water to the shoot, which would then expand the shoot desiccation time. Hence, it is critical to scale-up t crit estimations at higher levels, such as the whole-tree level. However, this is challenging as it will require data on crown area, volume of capacitive wood, and even root traits (Martin-StPaul et al. 2017), which is particularly difficult to measure or estimate for tall tropical rainforest trees.
Overall, we found consistent cross-species differences in key parameters (xylem embolism resistance, surface area, total transpiration, capacity per evaporative area) between leaves and stems ( Fig. 3c-e). However, we did not find tight relationships between the leaf-stem difference in these parameters and t crit (Fig. 4), except for embolism resistance. First, this supports the degree of vulnerability segmentation to be the leading parameter in determining t crit . Second, this does not necessary suggest that other parameters are not significant determinants, but it could be that the relative contribution of these parameters may vary depending on the species. In other words, the effect Fig. 3 Leaf-stem comparison in traits underlying vulnerability segmentation for the 12 species with a positive vulnerability segmentation (P 88,leaf > P 88,stem ). a Comparison in P 88 between leaf and stem across species. b Comparison in minimum conductance between leaf (g min ) and bark (g bark ) across species. c Comparison in surface area between leaves and bark across species. d Comparison between total transpiration between leaves and bark across species. e Comparison in capacity per evaporative area between the entire shoot (leaf + stem) and stem only. Leaf surface area is multiplied by two as both abaxial and adaxial sides are evaporative. All abbreviations are presented in Table 1 of interactions between parameters into the model may blur a direct relationship between a given trait and t crit . Third, the relative contribution of the different parameters may also be sensitive to the focused level of organisation (i.e. shoot, branch, whole plant). Therefore, more experimental (such as Blackman et al., 2019b) and theoretical (such as this present study) measurements of t crit are required to better identify (i) the relative contribution of parameters behind the effectiveness of vulnerability segmentation on t crit , and (ii) the relative contribution of traits determining t crit and drought resistance in general, which is critical for parsimonious but powerful parameterisations of models to predict plant drought resistance and future vegetation dynamics under climate change (McDowell et al. 2018;Brodribb et al. 2020).

Conclusion
In this present study, we investigated the effect of leaf-stem vulnerability segmentation on shoot t crit , as also the potential mechanisms by which vulnerability segmentation allows to expand t crit .
This study reports two important findings. (1) Vulnerability segmentation theoretically expands shoot t crit , suggesting that vulnerability segmentation is an important drought resistance mechanism.
(2) Increasing t crit is mostly driven by the difference in evaporative surface between leaf and bark-not by the difference between g min and g bark -and potentially the difference between shoot and stem capacity per evaporative surface. Consequently, we suggest that future trait-based assessments of plant drought resistance should incorporate vulnerability segmentation, to avoid underestimations of plant drought resistance. This is particularly required in the perspective of individual-based modelling aiming at predicting vegetation dynamics under different climate change scenarios. As plant drought resistance remains difficult to quantify, the trait-based model used in the present study is promising, but requires (i) further tests on the model assumptions, especially regarding the incorporation of vulnerability segmentation, (ii) comparisons with experimental t crit measurements and more mechanistic models (e.g. SurEau; Cochard et al. 2020)), and (iii) investigations and improvements to allow whole-tree level upscaling, which is an important issue for tropical rainforest canopy trees. Fig. 4 Relationships between shoot desiccation time and drought resistance traits for the 12 species with a positive vulnerability segmentation (P 88,leaf > P 88,stem ). a Shoot desiccation time according to the leaf-bark difference in minimum conductance. b Shoot desiccation time according to the leaf-bark difference in surface area. c Shoot desiccation time according to the leaf-bark difference in evapotranspiration. d Shoot desiccation time according to the total shoot vs. stem capacity per total evaporative surface area. Leaf surface area is multiplied by two as both abaxial and adaxial sides are evaporative. All abbreviations are presented in Table 1 Author contribution S.L, S.C, and P.H. designed the study; S.L, C.Z., S.C., C.S., E.C., J.-Y.G., and P.H. collected field samples; S.L., C.Z., S. C., D.B., J.-Y.G, and E.C. performed measurements or provided data; S.L. and C.Z. performed data analysis; S.L. wrote the manuscript with contribution of C.Z., S.J., P.H., and S.C; all authors discussed the results and contributed valuable comments on the manuscript.
Funding This study was funded by the GFclim project (FEDER 20142020, Project GY0006894). This work has benefited from an "Investissement d'Avenir" grant managed by Agence Nationale de la Recherche (CEBA, ref. ANR-10-LABX-25-01). S.L. was supported by a doctoral fellowship from CEBA.

Data availability
The datasets generated during the current study are available from the corresponding author on reasonable request.

Declarations
Ethical approval All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

Conflict of interest The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.