Independence of stem and leaf hydraulic traits in six Euphorbiaceae tree species with contrasting leaf phenology.

Hydraulic traits and hydraulic-related structural properties were examined in three deciduous (Hevea brasiliensis, Macaranga denticulate, and Bischofia javanica) and three evergreen (Drypetes indica, Aleurites moluccana, and Codiaeum variegatum) Euphorbiaceae tree species from a seasonally tropical forest in south-western China. Xylem water potential at 50% loss of stem hydraulic conductivity (P50(stem)) was more negative in the evergreen tree, but leaf water potential at 50% loss of leaf hydraulic conductivity (P50(leaf)) did not function as P50(stem) did. Furthermore, P50(stem) was more negative than P50(leaf) in the evergreen tree; contrarily, this pattern was not observed in the deciduous tree. Leaf hydraulic conductivity overlapped considerably, but stem hydraulic conductivity diverged between the evergreen and deciduous tree. Correspondingly, structural properties of leaves overlapped substantially; however, structural properties of stem diverged markedly. Consequently, leaf and stem hydraulic traits were closely correlated with leaf and stem structural properties, respectively. Additionally, stem hydraulic efficiency was significantly correlated with stem hydraulic resistance to embolism; nevertheless, such a hydraulic pattern was not found in leaf hydraulics. Thus, these results suggest: (1) that the evergreen and deciduous tree mainly diverge in stem hydraulics, but not in leaf hydraulics, (2) that regardless of leaf or stem, their hydraulic traits result primarily from structural properties, and not from leaf phenology, (3) that leaves are more vulnerable to drought-induced embolism than stem in the evergreen tree, but not always in the deciduous tree and (4) that there exists a trade-off between hydraulic efficiency and safety for stem hydraulics, but not for leaf hydraulics.

Introduction

The issue of divergences in functional trait, ecoclimatic adaptation, and global distribution for deciduous and evergreen tree species has long interested physiologists, ecologists, and biogeographer. Many researches have confirmed that seasonal drought can favor deciduous leaves, and that infertile soils can favor long-lived evergreen leaves (Aerts 1995; Salleo et al. 1997; Eamus 1999; Givnish 2002), but the deciduous and evergreen tree species can co-exist in most habitats (Sobrado 1991; Quigley and Platt 2003).

The comparison between the deciduous and evergreen tree species in the same habitats has been widely made. It has been generally accepted that, compared with the evergreen tree species, the deciduous tree species have higher photosynthetic rate, higher relative growth rate, and lower leaf mass per area (Reich et al. 1991; Villar et al. 1995; Cornelissen et al. 1996). On the other hand, lower hydraulic conductivity and lower hydraulic vulnerability have been observed in stems of the evergreen tree species (Choat et al. 2003, 2005; Chen et al. 2009). In contrast, there has also been evidence of the overlap and convergence in stem hydraulic traits for tree species with contrasting phenology (Bucci et al. 2004; Maherali et al. 2004; Choat et al. 2007). Thus, it seems likely that leaf life-span affects photosynthetic characteristics, but stem hydraulic traits may be independent of leaf phenology.

Nowadays, most studies on hydraulics of tree species have concentrated on property of roots and stems (Brodribb et al. 2003; Maherali et al. 2006; Choat et al. 2007), but relatively few have elucidated leaf hydraulics due to methodological barrier. Nevertheless, recent work has shown that resistance to water transport in leaves accounts for 30–80% of the total hydraulic resistance of the whole-plant water transport pathway (Becker et al.1999; Sack et al.2002), even though the water transport pathway in leaves represents a very small fraction of that in the whole plant. Furthermore, compared with the stems, leaves appear to be more sensitive to cavitations, often losing a substantial fraction of their hydraulic conductivity under non-extreme conditions (Brodribb and Holbrook 2004; Woodruff et al.2007; Hao et al. 2008). However, leaf life-span has never been taken into account when the comparison of leaf and stem hydraulic vulnerability is made.

Fortunately, many studies have shown that stem hydraulic traits of tree species are closely correlated with wood structure (Hacke et al. 2001; Rosner et al. 2007; Hao et al. 2008). Accordingly, it has been proposed that the divergence in stem hydraulic property for plants with contrasting phenology may result from the difference in stem functional structure regardless of leaf phenology (Chen et al. 2009). Likewise, Sack and Frole (2006) also found that leaf structural diversity was related to hydraulic capacity in tropical rain forest tree, implying that leaf hydraulic traits may depend on leaf structure, and thus may be independent of stem structure. To our knowledge, the comparison of stem hydraulic traits between the deciduous and evergreen tree species has been made, but the comparison of leaf hydraulic traits is yet unavailable.

In the present study, leaf and stem hydraulic traits and their hydraulic-related functional structure were examined in six tree species with contrasting phenology, and the correlation analysis between hydraulic traits and functional structural properties was also made. The six tree species co-occur in a seasonally tropical forest in south-western China. All of the six tree species belong to a single family (Euphorbiaceae), and thus the experimental results might largely rule out the impact of differences in phylogeny. Our working hypothesis was that there may be substantial differences in hydraulic traits between the evergreen and deciduous tree species, but hydraulic traits may be only related to functional structure regardless of leaf phenology, and thus that leaf hydraulic traits may be independent of stem hydraulic traits.

Materials and methods

Study site and plant materials

The study was carried out at Xishuangbanna Tropical Botanical Garden (21°41′N, 101°25′E, and 570 m a.s.l.), the Chinese Academy of Sciences in south-western Yunnan, China. The mean annual air temperature is about 21.7°C and the annual precipitation is about 1,560 mm. However, it has a strong seasonality that more than 85% of precipitation occurs between the months of May and October, thereby resulting in a pronounced rainy season from May to October, and a well-defined dry season from November to April. All measurements in the present study were performed in the middle of rainy season between July and August. Because rainfall frequently occurred in the measured periods, thus stems and leaves collected early in the morning were fully hydrated.

Six tree species of Euphorbiaceae [Hevea brasiliensis (Willd. ex A.Juss.) Muell. Arg., Macaranga denticulate (Bl.) Muell. Arg., Bischofia javanica Bl., Drypetes indica (Muell. Arg.) Pax et Hoffm., Aleurites moluccana (L.) Willd., Codiaeum variegatum (L.) A.Juss.] were selected for our study [see Kunming Institute of Botany, The Chinese Academy of Sciences (2006) Flora of Yunnan for species nomenclature]. H. brasiliensis is an evergreen tree species in its native habitat of South America, but interestingly it is deciduous in the present region and sheds all leaves at the end of February. M. denticulate and B. javanica are native deciduous tree species, and they shed all leaves at the end of January. New leaves of the three deciduous tree species completely flush before the onset of rainy season. Thus, leaf life-span of the three deciduous tree species is less than 12 months. Correspondingly, D. indica, A. moluccana, and C. variegatum are evergreen tree species, and the recorded leaf life-span of A. moluccana, D. indica, and C. variegatum at least exceeds 18, 21, and 24 months, respectively. Leaf life-span was calculated according to the method of Navas et al. (2003). The tree height and the diameter at breast height for C. variegatum is 2–3 m and 3–4 cm, respectively. The tree height for the remaining tree species are 5–7 m and the diameter at breast height are 15–20 cm. All sampled tree species selected in our study are fully exposed to sunlight and three to five mature individuals were located for each of the target species.

Stem hydraulic conductivity

Stem hydraulic conductivity (Kh) was measured on five to seven branches per species, taken from three to five different individuals. Branches of about 1.5 m in length were collected early in the morning, recut immediately under water, and returned to laboratory shortly. Branch segments of more than maximum vessel length and 7–9 mm in diameter were cut underwater, trimmed with a fresh razor blade, and connected to the hydraulic apparatus. Segments were perfused with filtered, 10 mmol oxalic acid solution. The perfusing solution was supplied from a reservoir raised above the flow meter to generate a gravity-induced pressure head of approximately 5 kPa. The segment was perfused with solution until a steady-state flow was attained, at which point the hydraulic conductivity of the segment was calculated as Kh = FL/∆P, where F is the flow rate (kg s−1), ∆P is the pressure drop (MPa) along the stem segment, and L is the length of stem segment (m). The steady-state flow for stem hydraulic measurements was evaluated by the required time for a given volume of water flow through a pipette connected to the water outlet side of the segments with a transparent rubber tube. If the required time for a given volume of water flow was recorded for three times, and it was relatively stable, we considered the water flow having reached a steady-state flow. Upon completion of the measurement, the segment was perfused with 0.01% Safranin dye to determine conductive sapwood area. Leaf surface area distal to the stem segment was measured using a leaf area meter (LI-Cor, Lincoln, NE, USA). Sapwood and leaf area were used to calculate sapwood area-specific hydraulic conductivity (KS, kg m−1 s−1 MPa−1) and leaf area-specific hydraulic conductivity (KL, kg m−1 s−1 MPa−1). Because the stems collected early in the morning were fully hydrated, thus maximum KS (KS-max) and KL (KL-max) was calculated as Kh divided by sapwood and leaf area, respectively.

In the present study, the segment length of more than the maximum vessel length was selected to reflect intrinsic efficiency of water transport for each species. The maximum vessel length was evaluated by air method measurements as described by Cohen et al. (2003) and Thorne et al. (2006). In brief, an approximately 40 cm unbranched stem or branch segment was cut directly from the plant. Segment diameters were all between 7 and 9 mm. The segment was connected to the hydraulic apparatus, and filtered; pressurized air at a maximum of 80 kPa, which was low enough to avoid inducing cavitations in xylem vessels, was applied to the cut end of segments. Then, the segments were placed under water and incisions were made at 1–2 mm intervals at the air outlet side. The first appearance of a stream of bubbles was evidence of an open xylem vessel, because any water in the open vessels would be pushed out in the first few seconds, and the stable bubbles were then observed. Thus, the distance between two sided of this segment was considered just as the maximum vessel length. Here, it should be pointed out that when low-pressure air is forced through xylem vessels of fresh tissue, air will only pass through open vessels because wet pit membranes will block air flow, and if vessels are already embolized, air will also be able to flow through those vessels, which could compromise the accuracy of continuous vessel length measurements. The measurements of the maximum vessel length were conducted on 7–9 segments from 3–5 individuals of each species.

Stem hydraulic vulnerability curves

Stem hydraulic vulnerability curve was determined by bench drying method on detached branches (Sperry et al. 1988). Large (about 1.5 m) branches from three to five individuals of each species were collected early in the morning, immediately sealed in plastic bags and transported to the laboratory. In the laboratory, these branches were allowed to dehydrate for different periods of time to reach a large range of water potential. The branches were then sealed into double layers of black plastic bags with wet paper towels for at least 2 h to equilibrate between different parts of the stems. A series of measurements were then made on branch tissue using a tubing apparatus to assess the percentage loss of hydraulic conductivity (PLC) as a result of embolism. The xylem water potential of branch segments was assessed by measuring the water potential of three to five leaves distal to the stem segment. During the measurements, we should carefully distinguish latex from water in some study species. A stem segment was then cut under water and connected to the hydraulic apparatus. For a fixed pressure head, Kh is proportional to volumetric flow rate (Jv) of water through stem segments, and thus in practice the calculations of PLC were based on Jv rather than Kh. After the initial measurement of flow rate (Ji), the segments were flushed using filtered, 10 mmol oxalic acid solution under a constant pressure of 0.15–0.18 MPa for 10–30 min to remove embolisms until stable readings of flow rate were reached. The maximum flow rate (Jmax) was then measured using the same pressure head as before. PLC was then calculated as: PLC = 100(Jmax − Ji)/Jmax. For each branch, three to five segments were selected for measurement, and the total number of measurement for each species was between 50 and 100. The vulnerability curves were generated for each species by fitting curves for PLC against xylem water potential (Pammenter and Vander Willigen 1998): PLC = 100/(1 + exp (a (Ψ − b))), where Ψ is xylem water potential, a is the slope of the line and b is the xylem water potential at which 50% loss of stem hydraulic conductivity occurs (P50stem). It should be stressed that water potential at 50% loss of hydraulic conductivity (P50) is commonly used to characterize and interpret interspecific variation in cavitation resistance (Tyree and Ewers 1991).

Leaf hydraulic conductivity

The branches from three to five individuals of each species were collected early in the morning and leaf hydraulic conductivity (Kleaf) was then measured according to the method as described by Franks (2006). In brief, Kleaf = Δν/(10 × Aleaf × (Ψ1 − Ψ2)), where Δv was the volume of sap expressed in the first 10 s following the step increase in pressure from Ψ1 to Ψ2, and Aleaf was leaf area (m2). The chamber pressure was increased to the balance pressure (Ψ1) and allowed to equilibrate for about 2 min. The chamber pressure was then rapidly increased to Ψ2, and the sap expressed in the first 10 s was collected and weight on analytical balance to obtain Δv. Ψ2 was more 0.5 MPa than Ψ1 (accurately measured to ±0.01 MPa). Leaf area was measured with a leaf area meter (LI-Cor, Lincoln, NE, USA). In our experiment, maximum leaf hydraulic conductivity (Kleaf-max, mmol m−2 S−1 MPa−1) was obtained due to full hydration of the sampled leaves collected in the early morning.

Leaf hydraulic vulnerability curves

In brief, the branches from three to five individuals of each species were cut early in the morning while leaf water potential (Ψleaf) was fully hydrated, and most leaves were removed except for terminal clusters of six to eight leaves. These branches were then allowed to desiccate very slowly ensuring all leaves remained at similar Ψleaf. Periodically, branches were bagged and placed in dark for 30 min to ensure that stomata were closed and Ψleaf was homogenous among leaves. Two leaves were then removed to gauge the Ψleaf remaining on the branch, after which the remaining leaves were detached to measure Kleaf according to the method described as above: Kleaf = Δv/(10 × Aleaf × (Ψ1 − Ψ2)). Ψ1 was defined as Ψleaf, and leaf hydraulic vulnerability curves were generated by plotting Kleaf against Ψleaf with an exponential sigmoid equation: Kleaf = a/(1 + exp (−b (Ψleaf − c))), where a is maximum Kleaf, b is the slope of the line, and c is Ψleaf at which 50% loss of maximum Kleaf occurs (P50leaf). It should be emphasized that this method used to measure Kleaf involves applying a positive pressure to squeeze water out of the leaf, and this may refill cavitations in the leaf xylem when the water is pushed through the leaf at substantial pressure. Consequently, this might lead to an overestimation of leaf hydraulic conductivity for dehydrated leaves, and thus to an underestimation of P50. Therefore, cautions should be used when comparing values of Kleaf and P50 given here with values determined using different methods in other studies.

Leaf and stem structural property

Leaves in branches sampled for the measurement of leaf hydraulic traits were also used to determine leaf mass per area (LMA, g cm−2) and leaf anatomical structure. The thickness of the leaf blade, upper epidermis, lower epidermis, palisade mesophyll, and spongy mesophyll were determined from a micrograph at 200× magnification. Vessel anatomical traits were determined from transverse sections of xylem tissue. Transverse sections of the xylem were made, at least 50–70 vessels diameter was measured in each of the samples to determine mean vessels diameter (Dm–v), and the vessels number (VN) per square millimeter was counted. Sapwood density (g cm−3) was determined for each sample that was tested hydraulically. Sapwood with both bark and pith removed was immersed in tap water overnight to saturate the samples. After the surface was wiped dry, the volume of sapwood was measured immediately using the water displacement method. Sapwood density = Md/Vf, where Md is the oven dry mass (dried at 80°C for 72 h), and Vf is the fresh volume.

Statistical analysis

All analyses were performed with the SPSS software package (Chicage, IL), and the variables were given as mean ± SD. The comparison of hydraulic traits between evergreen and deciduous tree was confirmed by t tests with species nested within leaf habit type. Pearson’s correlations were made among hydraulic traits and structural properties. In addition, only the P values of less than 0.05 could be accepted as a significant level.

Results

The stem hydraulic vulnerability curves indicated that values of P50stem for the three deciduous tree species of H. brasiliensis, M. denticulate, and B. javantica were −1.27, −1.14, and −1.27 MPa, respectively; correspondingly, values of P50stem for the three evergreen tree species were more negative than −2.0 MPa (Fig. 1). However, values of P50leaf in the six tree species overlapped considerably when considering leaf habit types (Fig. 2). Thus, compared with the deciduous tree species, P50stem was significantly negative in the evergreen tree species (P < 0.05), but P50leaf did not show a significant difference between the evergreen and deciduous tree species (Fig. 3). In addition, P50stem was more negative than P50leaf in the evergreen tree species (P < 0.05), but such a pattern was not observed in the deciduous tree species (Fig. 3).

Fig. 1

Vulnerability to water stress-induced embolism in six Euphorbiaceae tree species: H. brasiliensis (HB), M. denticulate (MD), B. javanica (BJ), D. indica (DI), A. moluccana (AM), and C. variegatum (CV). Each point represents one measurements and percentage loss of conductivity (PLC) was plotted against xylem tension for each tree species. The vertical dashed lines indicated the xylem tension at which 50% loss of hydraulic conductivity occurs (P50)

Fig. 2

Leaf hydraulic conductivity (Kleaf) as a function of leaf water potential for six Euphorbiaceae tree species: H. brasiliensis, M. denticulate, B. javanica, D. indica, A. moluccana, and C. variegatum. Values for each point are mean ± SD of more than ten measurements on the slowly dark-desiccated leaves from 3 to 5 individual trees (n = 3–5), and a sigmoid function was fitted to the data. The vertical dashed lines indicated the leaf water potential at which 50% loss of maximum Kleaf occurs (P50)

Fig. 3

Water potential at 50% loss of hydraulic conductivity (P50) for stem (shaded rectangle) and leaf (striped rectangle) in the deciduous and evergreen Euphorbiaceae tree species. The mean ± SD is shown for the values of P50 (n = 3). Within the graph, different lower case letters represent significant differences (P < 0.05, t tests)

The values of Kleaf-max fluctuated between 8 and 20 mmol m−2 s−1 MPa−1 in the six tree species, but Kleaf-max did not exhibit a significant difference between plants with contrasting phenology (Fig. 4a). Of the six tree species, the minimum value of Kleaf-max was observed in the deciduous tree species of H. brasiliensis, followed by the evergreen tree species of D. indica. On the other hand, the values of KL-max and KS-max were generally lower in the three evergreen tree species except for the evergreen tree species of A. moluccana, the KL-max and KS-max of which was as high as those of the deciduous tree species (Fig. 4b, c). Especially, the value of KS-max for the deciduous tree species was significantly higher than that for the evergreen tree species (P < 0.05; Fig. 4c).

Fig. 4

Maximum leaf hydraulic conductivity (Kleaf-max, a), maximum leaf area-specific hydraulic conductivity (KL-max, b), and maximum sapwood area-specific hydraulic conductivity (KS-max, c) in six Euphorbiaceae tree species with contrasting leaf phenology. HB, H. brasiliensis; MD, M. denticulate; BJ, B. javanica; DI, D. indica; AM, A. moluccana; CV, C. variegatum; DT, deciduous tree, and ET, evergreen tree. For individual hydraulic traits of six tree species, values are mean ± SD of more than ten measurements on leaves or stem from 3 to 5 individual trees (n = 3–5). For hydraulic traits of deciduous and evergreen tree species, values are mean ± SD (n = 3) and different lower case letters represent significant differences (P < 0.05, t tests) within the graph

Structural properties of leaves, such as the thickness of upper epidermis, lower epidermis, palisade mesophyll, spongy mesophyll, and leaf blade and LMA overlapped considerably between the deciduous and evergreen tree species (Table 1). In contrast, structural properties of stems were markedly different between plants with contrasting phenology (Table 1). In comparison with the deciduous tree species, the evergreen tree species generally exhibited higher sapwood density, lower Dm–v, and higher VN per square millimeter.

Table 1

Leaf and stem structural traits in six Euphorbiaceae tree species: H. brasiliensis, M. denticulate, B. javanica, D. indica, A.moluccana, C. variegatum

Tree species	Leaf phenology	UE (μm)	LE (μm)	PM (μm)	SM (μm)	LT (μm)	LMA (g/m2)	ρ (g/cm−3)	Dm–v (μm)	VN (mm−2)
H. brasiliensis	Deciduous	25.0 ± 0.9	28.6 ± 0.8	100.2 ± 8.3	98.0 ± 4.2	251.8 ± 4.9	74.50 ± 8.9	0.48 ± 0.04	80.9 ± 12.6	27 ± 3.9
M. denticulate	Deciduous	26.0 ± 1.2	20.5 ± 1.3	114.1 ± 7.5	67.0 ± 3.2	227.7 ± 7.3	54.08 ± 6.8	0.41 ± 0.03	97.3 ± 10.9	21 ± 4.3
B. javanica	Deciduous	67.6 ± 1.9	34.8 ± 1.8	185.4 ± 11.0	276.3 ± 15.2	564.1 ± 6.5	157.51 ± 18.0	0.43 ± 0.04	85.4 ± 13.9	38 ± 5.5
D. indica	Evergreen	30.4 ± 0.9	23.1 ± 0.5	80.6 ± 5.2	155.6 ± 13.1	289.7 ± 10.1	105.20 ± 11.3	0.67 ± 0.04	33.6 ± 4.5	135 ± 9.1
A. moluccana	Evergreen	45.4 ± 3.6	31.9 ± 3.2	129.6 ± 12.3	154.7 ± 9.1	361.6 ± 8.1	82.17 ± 9.8	0.52 ± 0.03	49.2 ± 3.2	69 ± 13.1
C. variegatum	Evergreen	47.1 ± 2.5	61.5 ± 5.2	162.2 ± 10.8	251.5 ± 14.6	522.2 ± 13.5	104.87 ± 14.5	0.49 ± 0.05	38.0 ± 5.6	92 ± 10.1

The variables are given as mean ± SD (n = 3–5)

UE upper epidermis, LE lower epidermis, PM palisade mesophyll, SM spongy mesophyll, LT leaf thickness, LMA leaf mass per area, ρ sapwood density, Dm–v mean vessel diameter, VN vessels number per square millimeter

Person’s correlations between hydraulic traits and structural properties were made across the deciduous and evergreen tree species (Table 2). LMA was independent of stem and leaf hydraulic traits. However, negative correlations of P50leaf with the thickness of upper epidermis, lower epidermis, palisade mesophyll, and leaf blade were observed; in contrast, positive correlations of Kleaf-max with the thickness of upper epidermis, palisade mesophyll, and leaf blade were found across the six tree species with contrasting phenology.

Table 2

Pearson’s correlation coefficients between the traits of leaf and stem in six Euphorbiaceae tree species: H. brasiliensis, M. denticulate, B. javanica, D. indica, A.moluccana, C. variegatum

	P50leaf	Kleaf-max	UE	LE	SM	PM	LT	LMA	P50stem	KL-max	KS-max	ρ	Dm-v	VN
P50leaf		−0.37	−0.70*	−0.69*	−0.54	−0.91***	−0.73**	−0.45	0.29	−0.66*	−0.50	0.65*	−0.36	0.30
Kleaf-max			0.76**	0.31	0.47	0.71*	0.62*	0.32	0.17	0.23	0.14	−0.30	0.04	0.15
P50stem			0.21	0.40	0.47	−0.14	0.29	0.27		0.66*	0.80**	0.79**	0.98***	−0.96***
KL-max			0.21	0.04	−0.17	0.53	0.07	−0.20			0.93***	−0.83**	0.79**	−0.78**
KS-max			0.11	−0.25	−0.28	0.36	−0.08	−0.17				−0.78**	0.91***	−0.85**

To facilitate log10-transformation for analysis, P50leaf and P50stem values were converted from negatives to positives. *P < 0.10, **P < 0.05, ***P < 0.01

P50 leaf water potential at 50% loss of leaf hydraulic conductivity, K maximum leaf hydraulic conductivity, P50 xylem water potential at 50% loss of stem hydraulic conductivity, K maximum leaf area-specific hydraulic conductivity, K maximum sapwood area-specific hydraulic conductivity, UE upper epidermis, LE lower epidermis, PM palisade mesophyll, SM spongy mesophyll, LT leaf thickness, LMA leaf mass per area, ρ sapwood density, D mean vessel diameter, VN vessels number per square millimeter

However, no correlations between stem hydraulic traits and leaf structural properties were observed, but stem hydraulic traits were significantly correlated with stem structural properties across the six tree species. Furthermore, no correlations of leaf hydraulic traits with stem anatomical properties, such as Dm–v and VN, were found across the six tree species. On the other hand, P50stem was significantly correlated with KS-max, but P50leaf was uncorrelated with Kleaf-max (Table 2). In addition, no correlation of leaf hydraulic traits with stem hydraulic traits was observed across the six tree species, except that P50leaf was weakly correlated with KL-max.

Discussion

The value of Kleaf-max obtained in our experiments is within the range of Kleaf-max reported by Sack and Holbrook (2006), who found that the average value of Kleaf-max in tropical wood angiosperms was about 15 mmol m−2 s−1 MPa−1. Beyond our expectations that the value of Kleaf-max may be higher in short-lived deciduous leaf, the divergence in Kleaf-max was not observed between the short-lived deciduous leaf and long-lived evergreen leaf. Nevertheless, it has been reported that the value of Kleaf-max is highly variable even within a life form, varying by tenfold among coexisting tree species (Sack et al. 2005). Furthermore, the focused comparison of Kleaf-max between the deciduous and evergreen tree species has not been made so far. We analyzed the data of the study four tree species of Brodribb and Holbrook (2003b), and also found that there was not a significant difference in Kleaf-max between two deciduous and two evergreen tree species. Thus, it seems likely that leaf life-span may not have an effect on Kleaf-max.

In the six tree species, leaf hydraulic resistance to embolism was also highly variable, but long-lived evergreen leaf did not have more negative P50 (Fig. 3), indicating that the capacity of a leaf to resist drought-induced losses of hydraulic function may not affect a leaf’s life-span. Originally, we had expected that short-lived deciduous leaf is more vulnerable to drought-induced cavitations, because many researches have shown that, because of the failure to protect from drought-induced losses of hydraulic function, the deciduous tree species avoids seasonal drought by shedding leaves (Gill and Mahall 1986; Holbrook et al. 1995; Brodribb and Holbrook 2003a). However, the adaptation of short-lived deciduous leaf to seasonal drought was not reflected in the capacity of a leaf to resist embolisms estimated by P50. Based on the previously reported leaf hydraulic vulnerability curves of four tree species with contrasting phenology (Brodribb and Holbrook 2003b), it has not also been found that short-lived deciduous leaf was more vulnerable to water stress-induced cavitations. The overlaps in leaf hydraulic traits for the deciduous and evergreen leaf may result from leaf structural properties, because leaf structural properties overlapped considerably in the six tree species (Table 1).

Evidently, we found that the value of Kleaf-max was related to the thickness of upper epidermis, palisade mesophyll, and leaf blade (Table 2). Similarly, there has been evidence that the value of Kleaf-max correlates with palisade thickness and palisade/spongy mesophyll ratio for tropical rainforest tree species (Sack and Frole 2006), and with total leaf thickness and water storage capacitance per area for temperate woody species (Sack et al. 2003). Meanwhile, some research has shown that leaves with high Kleaf-max tend to have wider xylem conduits in the midrib and higher venation densities (Aasamaa et al. 2001; Sack et al. 2004). Recently, Brodribb et al. (2007) have also found that the post-vein traverse as determined by characters such as vein density, leaf thickness, and cell shape was strongly correlated with the hydraulic conductivity of foliage. The most likely basis for the relationships is that hydraulic traits and mesophyll traits converge for effective whole-leaf design in given habitats (Sack and Holbrook 2006). Species adapted to establishment in high irradiance would benefit from high vein density, conferring higher hydraulic capacity to supply these species higher maximum gas exchange rates, and also providing effective phloem translocation. Interestingly, we also found that leaf hydraulic resistance to embolism was significantly correlated with the thickness of palisade mesophyll (Table 2).

Dissimilarly, stem hydraulic traits displayed a pronounced difference between the plants with contrasting phenology. Obviously, stem hydraulic resistance to drought-induced embolism was weaker in the deciduous tree species (Fig. 3). It is consistent with the findings of Choat et al. (2003) that the deciduous tree species was more susceptible to water stress-induced embolism than the co-existing evergreen tree species. In contrast, the overlap in P50stem was also found in the comparisons between 62 deciduous tree species and 66 evergreen tree species (Maherali et al. 2004). Thus, it has been proposed that the structure of xylem vessels should be seen as an important factor in determining the occurrence of water stress-induced embolism (Zimmermann 1983), and that the vulnerability of xylem to water stress-induced embolism is more closely related to the porosity of inter-vessel pit membranes (Choat et al. 2003). More recently, there have been evidences of the strongly positive correlation between P50stem and mean vessel diameter (Maherali et al. 2006; Chen et al. 2009). In addition, our results also confirmed that stem hydraulic resistance to cavitations is significantly related to xylem functional structure, such as Dm–v and VN (Table 2).

Indeed, the differences in vessels anatomical properties cannot fully explain the divergence in the resistance of xylem to cavitations. Wood density in many tree species from different ecosystems tend to be correlated with a variety of hydraulic-related traits such as the capacity of stem water storage, the efficiency of xylem water transport, the regulation of leaf water status, and the avoidance of turgor loss (Meinzer 2003; Bucci et al. 2004; Gartner and Meinzer 2005). Thus, not surprisingly, xylem becomes more resistant to water stress-induced embolism with increasing wood density (Hacke et al. 2001; Hao et al. 2008). Similarly, sapwood density was significantly correlated with stem hydraulic resistance to cavitations in our study six tree species. Therefore, we think that, regardless of leaf phenology, a series of sapwood characteristics such as diameter and length of vessel, porosity, strength, and durability of pit membrane, and wood density affects stem hydraulic resistance to drought-induced embolisms.

On the other hand, the deciduous tree species generally had higher stem hydraulic efficiency, especially KS-max (Fig. 4c). These results are similar to those reported for the evergreen and deciduous tree species from a seasonally dry forest in Venezuela, in which KL-max was 2–4 times lower and KS-max was 2–6 times lower in the evergreen tree species (Sobrado 1993). In contrast, Goldstein et al. (1989) early found that two evergreen tree species from a Venezuelan tropical savanna had higher KL-max and KS-max than two co-existing deciduous tree species. More recently, Brodribb et al. (2002) have observed that the values of KS-max overlapped substantially for a range of the evergreen and deciduous tree species from a dry forest in Costa Rica, although the highest values recorded were for two deciduous species. Similarly, in our study, KL-max for the evergreen tree species of A. moluccana was higher than the deciduous tree species of H. brasiliensis and M. denticulate, and KS-max for the former also approached that for the latter (Fig. 4b, c). In addition, from the Hagen-Poiseuille law, it is predicted that the wider vessels of tree species would result in higher flow rates for a given pressure gradient. This is confirmed by our results that KL-max and KS-max were positively correlated with Dm–v (Table 2). Therefore, the differences in hydraulic efficiency between the evergreen and deciduous tree species possibly rely mainly on the divergence in xylem structure.

Surprisingly, we found that leaf hydraulic traits were uncorrelated with stem hydraulic traits except for the weak correlation of P50leaf with KL-max, and that leaf and stem hydraulic traits was only correlated with leaf and stem structure, respectively (Table 2). This indicates that leaf and stem each have relatively independent hydraulic system. In addition, P50stem was more negative than P50leaf in the evergreen tree species, but this did not always hold true in the deciduous tree species (Fig. 3). Thus, our result is not fully consistent with the findings of Hao et al. (2008), who found that, because of leaf hydraulic traits being under stronger selective pressure than being stem hydraulic traits in savanna and adjacent forest habitats with high evaporative demand and potentially low soil water availability, leaves were more vulnerable to drought-induced cavitations than terminal braches regardless of genus. However, our findings indicate that leaves function as a safety valve to a greater extent in the evergreen than deciduous tree species.

The correlation of hydraulic vulnerability with hydraulic efficiency has been observed within (Kavanagh et al. 1999; Domec and Gartner 2002; Rosner et al. 2006) and across species (Piñol and Sala 2000; Maherali et al. 2004). Similarly, P50stem was significantly related to KS-max within our study six tree species with contrasting phenology. Understandably, tree with higher hydraulic conductivity have wider vessel diameter, and such vessel structure is more susceptible to water stress-induced embolism. This line of reasoning is also confirmed by our results that stem hydraulic efficiency (KL-max and KS-max) was positively related to Dm–v, while P50stem was negatively correlated with Dm–v (Table 2). Actually, this indicates that tree employs different strategies of acquiring water resource based on a trade-off between hydraulic safety and efficiency. Higher hydraulic efficiency might guarantee sufficient acquisition of water resource for tree with higher hydraulic vulnerability before embolism occurs, while low hydraulic vulnerability would ensure tree with low hydraulic conductivity sufficient acquisition of water resource by delaying the occurrence of embolism.

Unexpectedly, one of the most commonly investigated patterns of a trade-off between hydraulic efficiency and resistance to drought-induced cavitations in stem hydraulics was not found in leaf hydraulics of the study six tree species, because Kleaf-max was not correlated with P50leaf (Table 2). In the analysis of the available data for 13 species of ferns, trees, and herbs varying 42-fold in Kleaf-max, Sack and Holbrook (2006) found that P50leaf ranged from −1.3 to −3.0 MPa, and was uncorrelated with Kleaf-max. It is widely accepted that stomata closure facilitates preventing the formation of embolisms (Pockman and Sperry 2000; Martínez-Vilalta et al. 2002; Froux et al. 2002). However, we consider that stomata closure may regulate losses of leaf hydraulic function to a larger extent, and losses of stem hydraulic function to a lesser extent. Then, the reason for the absence of a trade-off between leaf hydraulic efficiency and safety may be that the existence of stomata regulation for losses of hydraulic function disrupts such a trade-off.

In conclusion, the evident differences between the evergreen and deciduous tree species were found in stem hydraulic traits, but not in leaf hydraulic traits. However, regardless of leaf or stem, their hydraulic traits might result primarily from functional structure, and not from leaf phenology. On the other hand, the findings of this comparative study was that leaves are more vulnerable to drought-induced loss of hydraulic function than stem among the evergreen tree species, but this was not always the case among the deciduous tree species. In addition, our results also demonstrated that there exists a trade-off between hydraulic efficiency and safety for stem hydraulics, but not for leaf hydraulics, and the reason may be that stomata regulation for drought-induced losses of hydraulic function disrupts a trade-off between leaf hydraulic efficiency and safety.