Skip to main content

Evaporative flux method of leaf hydraulic conductance estimation: sources of uncertainty and reporting format recommendation



The accurate estimation of leaf hydraulic conductance (Kleaf) is important for revealing leaf physiological characteristics and function. However, the Kleaf values are largely incomparable in previous studies for a given species indicating some uncertain influencing factors in Kleaf measurement.


We investigated the potential impacts of plant sampling method, measurement setup, environmental factors, and transpiration steady state identification on Kleaf estimation in Oryza sativa and Cinnamomum camphora using evaporation flux method (EFM). The effects of sampling and rehydration time, the small gravity pressure gradients between water sources and leaves, and water degassing on Kleaf estimation were negligible. As expected, the estimated steady flow rate (E) was significantly affected by multiple environmental factors including airflow around leaf, photosynthetically active radiation (PARa) on leaf surfaces and air temperature. Kleaf decreased by 40% when PARa declined from 1000 to 500 µmol m−2 s−1 and decreased by 15.1% when air temperature increased from 27 to 37 °C. In addition, accurate steady-state flow rate identification and leaf water potential measurement were important for Kleaf estimation.


Based on the analysis of influencing factors, we provided a format for reporting the metadata of EFM-based Kleaf to achieve greater comparability among studies and interpretation of differences.


The influences of measurement setup, environmental factors, and steady state identification on leaf hydraulic conductance measurement were estimated and reporting format of Kleaf were proposed.


Plant hydraulic properties strongly influence photosynthesis and growth. At a given soil water potential, the capacity of leaves to maintain stomata open for photosynthesis mainly depends on plant hydraulic conductance [1]. On average, the hydraulic resistance in leaves accounts for 30% of that in whole-plant, thus leaves constitute an important bottleneck for hydraulic conductance [2]. Due to the important role of leaf hydraulic conductance (Kleaf) in plant, Kleaf has been extensively studied over the last decades [3,4,5,6,7]. The efficiency of water transport from the petiole to the sites of evaporation through the leaf tissues is quantified as Kleaf, and it is generally expressed as water transport efficiency per unit leaf area (mmol m−2 s−1 MPa−1). A number of approaches have been used to estimate Kleaf based on excised leaves, such as the evaporative flux method (EFM), the rehydration kinetics methods (RKM), the high-pressure flowmeter (HPFM), and the vacuum pump method (VPM) etc.

The Kleaf is typically estimated by measuring the ratio of water uptake rate through the leaf to the driving force (e.g., water potential gradient between the petiole to evaporation sites, ΔP). Currently, many methods developed to measure Kleaf, and the EFM has advantage of mimicking the natural transpiration pathways of water movement in the leaf [5]. Furthermore, the EFM method allows other functional traits, such as CO2 assimilation rate and stomatal conductance, to be measured simultaneously. Actually, the EFM has been used to measure stomatal conductance by recording the air humidity, air temperature around the leaf and steady-state water flow rate [8]. As the Kleaf has also been used to explain a range of physiological processes related to photosynthesis, drought tolerance and leaf economical spectrum [9,10,11], simultaneous estimation of Kleaf and other traits will provide more reliable information for understanding plant performance under variable conditions.

While the EFM was used in estimating Kleaf frequently, the Kleaf values estimated by EFM from different groups are largely incomparable even in the same species (Additional file 1: Fig. S1). Great difference in Kleaf has been found in the model species, Arabidopsis thaliana. In some studies [12,13,14], Kleaf of the Columbia (a widely selected ecotype of A. thatliana) was less than 10 mmol m−2 s−1, but it was larger than 50 mmol m−2 s−1 in other studies [15]. Surprisingly, the huge difference in Kleaf values estimated using EFM was even found in a Oryza sativa genotype, Shanyou 63, and the Kleaf values varied greatly from 0.64 to 23 mmol m−2 s−1 in previous studies [16,17,18,19] (more details in Additional file 1: Fig. S1 and Additional file 2: Table S1).

Different Kleaf among studies may be induced by multiple growth environmental factors such as light, temperature, humidity during plant growth. However, it seems unlikely that growth conditions in these studies could have led to such huge differences in Kleaf values [19,20,21]. Scoffoni et al. estimated the Kleaf of six species of lobeliads grown in two irradiances (daily average of 300 vs 800 μmol photons m−2 s−1) and found the largest variation was only 2.5‐fold in Kleaf [21]. Alternately, the differences in Kleaf among studies may be partially due to measurement bias. The environmental irradiance and temperature can influence Kleaf measurement as shown on a range of species, and they must be controlled accurately [16, 22]. However, the more comprehensive sources of uncertainty in using EFM have been less investigated, and, in fact, the details of the Kleaf measurements such as the sampling time, sample selection criteria, temperature, photosynthetically active radiation (PARa) and solutions used for Kleaf measurement were not all available in many studies [16, 17, 20, 23,24,25,26], raising the need for the establishment of transparent and detailed method descriptions and protocols.

The lack of Kleaf estimation and reporting format makes the full and efficient use of Kleaf from other studies difficult. It is essential to explore the estimation and the reporting format of Kleaf for unifying and normalizing Kleaf data from different sources. The study aims to investigate the effects of interference factors during measurement on initial(ψinitial) and final leaf water potential (ψfinal), flow rate, transpiration rate and leaf hydraulic conductance (Kleaf) and to provide the detail recommendations for EFM application and results report.


Kleaf measurements were performed in O. sativa and C. camphora leaves collected at different daily time to test the impacts of sampling time and the sample storage in the lab. The ψinitial of the samples rehydrated in the lab overnight was significantly lower than the one sampled in the morning of the measurement day in O. sativa but not in C. camphora (Fig. 1a). However, the decreased ψinitial of the over nightly hydrated samples had no influence on Kleaf. Actually, no differences in both E and final leaf water potential (ψfinal) were observed between the leaves sampled at the previous night and ones sampled in the morning of the measurement day. The water potentials of nightly sampled leaves were measured three times along the rehydration process. The results showed that ψinitial of leaves sampled at previous night was high and convergent within the first 12 h but decreased after 22 h of hydration. Furthermore, recutting in the morning or storage in sterile water alleviated this decrease (Additional file 1: Fig. S3, Fig. 1). No differences in ψinitial, E, ψfinal, and Kleaf were observed between the O. sativa leaves rehydrated in sterile water and the leaves rehydrated in non-sterile water (Fig. 1).

Fig. 1
figure 1

AD Influences of sampling time on initial water potential (ψinitial), steady flow rate (E), leaf hydraulic conductance (Kleaf), and final water potential (ψfinal), and EH the impacts of storage water sterilization on ψinitial, E, Kleaf, and ψfinal of previous sampled leaves. Green points and yellow points represent O. sativa and C.camphora leaves, respectively. Each point represents one individual leaf. Note, storage water sterilization effects were only investigated in O. sativa. (ns no significance; ***P < 0.001)

We evaluated the impacts of water degassing on Kleaf estimation by comparing degassed water and distilled water used as water source in the cylinder. Our data showed that E, ψfinal, and Kleaf estimations in both species were not impacted by water degassing (Fig. 2). Then, the influences of height gradients between water source and leaf blade were investigated, since the height pressure difference may exist between leaf and water source in cylinder. We found that 2 cm height gradient between leaf and water surface in cylinder exhibited no effect on Kleaf measurement (Fig. 3). In addition, the influences of cylinder water evaporation on E estimation under multiple conditions was quantified. The evaporation rate of water in the cylinder without any intervene on water was 0.115 × 10−3 mmol s−1. The evaporation rate in the cylinder was significantly declined by covering the water surface using liquid wax and/or by maintaining high humidity in weighting chamber (ANOVA, P < 0.001) (Additional file 1: Fig. S2).

Fig. 2
figure 2

Influences of water degassing on A steady flow rate (E), B leaf hydraulic conductance (Kleaf), and C final water potential (ψfinal). Green points and yellow points represent O. sativa and C. camphora leaves, respectively. Each point represents one individual leaf. (ns no significance)

Fig. 3
figure 3

Impacts of height difference between leaf and water surface in cylinder on A steady flow rate (E), and B final water potential (ψfinal), and C leaf hydraulic conductance (Kleaf). 0 indicates same height level between leaf and water surface; − 2 indicates that water surface is 2 cm lower than leaf; 2 indicates that water surface is 2 cm higher than leaf. Green points and yellow points represent O. sativa and C. camphora leaves, respectively. Each point represents one individual leaf. (ns no significance)

The effects of environmental factors including PARa, air temperature, and the airflow through leaf surface on O. sativa Kleaf estimation were investigated. Kleaf at 1000 µmol m−2 s−1 PARa was significantly higher than the values at 500 or 1500 µmol m−2 s−1. The higher Kleaf under 1000 µmol m−2 s−1 PARa was caused by the higher E and higher ψfinal. Low Kleaf were found under high air temperature condition due to the declined ψfinal. Interestingly, although the airflow had the strong effects on transpiration, the Kleaf values estimated under different airflow rates were not significant (Table 1).

Table 1 Influences of ambient photosynthetically active radiation (PARa), temperatureand airflow rate on E, ψfinal and Kleaf in O. sativa

Identifying steady state of flow rate is important because the Kleaf is directly calculated by using the steady water flow rate through leaves. However, our data indicated that the flow rates of some leaves were oscillated and downward, which did not conform to our steady state criteria (Fig. 4). In the current study, about two of three measurements achieved the continuous and steady state of flow rate according to our criteria (Fig. 4 and supplementary raw data of flow rate). The stabilization of flow rate was further confirmed by the estimation of transpiration rate using the gas exchange system (Additional file 1: Figs. S4, S5, S7), and the flow rate values from the two recording systems were only consistent when the leaves were entirely covered by the gas exchange chamber (Fig. 5). Importantly, the variation of ψfinal greatly affected Kleaf. At a given E (e.g., using the E of 8.7 mmol m−2 s−1 in O. sativa under PARa of 1000 µmol m−2 s−1), Kleaf decreased sharply (up to 3-folds) with the decreased ψfinal, especially within the typical O. sativa water potential range of − 0.17 to − 0.45 MPa we observed in the present study, which showed the inverse correlation between Kleaf and ψfinal (Fig. 6). The ψfinal of leaves acclimated at both 0 and 1000 µmol m−2 s−1 irradiance in advance decreased after 60 and 150 min equilibration in zip-lop bag, respectively (Additional file 1: Fig. S8). Conservatively, we showed that 10 min ~ 1 h equilibration time was proper for water potential measurement.

Fig. 4
figure 4

Four typical curves and segment of flow rate to time during measurement. Different segments were marked with different color. A, B Non-steady curves identified as the slope of last 15 min segment was not equal to zero. C Non-steady curve identified as measurement time corresponding to the last segment was less than 15 min. D Steady curve, the last segment was longer than 15 min and the slope of last 15 min segment was equal to 0. P represents the P-value of in the last 15 min segment (see details in “Materials and methods” section). ψfinal is the initial leaf water potential

Fig. 5
figure 5

Correlations of stable flow rate measured by balance (E) and by gas exchange system Licor 6800 (Elicor) in O. sativa (green) and C. camphora (yellow). Entire (circle) and partial leaves (triangle) were clamped in gas exchange chamber, respectively. The fitted equations and their p-values for camphor and rice are presented in the upper left of the figure. Blue/red/black fitting line and equation represents the correlation in partly clamped O. sativa leaves/partly clamped C. camphora leaves/entirely clamped C. camphora leaves, respectively

Fig. 6
figure 6

Sensitivity analysis of the influences of final water potential (ψfinal) and steady flow rate (E) on Kleaf. X-axis 0% represents mean value of ψfinal and E in O. sativa measured under 1000 µmol m−2 s−1 PARa. The x-axis ranges of ψfinal and E represent their measured value change ranges, respectively

According to our evaluation, we proposed a reporting format for Kleaf based on EFM to improve reuse and reanalysis valuable data. Not only the parameters used to calculate Kleaf but also the raw data of flow rate relative to time should be conserved and provided to meet the needs for identifying flow rate stabilization. Moreover, the environmental factors including airflow, temperature and PARa also should be provided as the Kleaf estimation was strongly affected by those factors (Additional file 3: Table S2, Additional file 4: Table S3).


In EFM, Kleaf is calculated by the ratio of water flow rate to the water potential gradient driving water movement across the leaf, and we showed that leaf water potential and E were potentially influenced by environmental factors and the criterion for determining physiological stabilization. Our investigation suggested that the remarkable variance in Kleaf among previous studies might be attributed to methodological artifacts. Therefore, providing details relating to Kleaf measurement is important to synthesize and to compare data across different studies.

Influences of sample collection and storage on K leaf estimation

Although the field material sampling at the night before the measurement day and rehydrating overnight in the lab are common practice in previous studies [5, 26,27,28,29], to the best of our knowledges, no study investigated the impacts of sampling and storage time on Kleaf estimation. The tension in xylem is created by transpiration, but the excessive tension in the daytime poses the inherent risk for xylem cavitation, thus substantially reducing the hydraulic conductivity of plants [30, 31]. However, no difference in Kleaf between leaves sampled at the previous night and leaves sampled in the morning was observed in both species. Rehydration overnight was projected to refill cavitated xylem conduits and restore the lost water in tissues caused by transpiration [32]. Interestingly, our data showed that the ψinitial of overnight rehydrated O. sativa leaves was more negative than that leaves sampled in the morning, which wasn’t observed in C. camphora (Fig. 1). The more negative ψinitial of leaves rehydrated overnight may result from xylem being blocked by mechanical wound secreta at the cut, such as callose, suberin, lignin, chitinase, and various phenolics [33,34,35]. Another explanation may be accumulation of sucrose at the cut surface caused by damaged cell wall integrity [36, 37], accelerating microbe attack and water uptake pathway blockage. Indeed, the ψinitial of the overnight rehydrated leaves increased quickly after recutting (Additional file 1: Fig. S3), which might be attributed to the reversible outside-xylem dehydration [38]. Further efforts are still needed to identify the types of microbes and/or the components of the blockage.

Effects of degassing and gravity pressure on K leaf estimation

Water transport in plants occurs under tension, the air-seeding and its expansion in water will create air-vapor embolisms and thus blocked water transportation. To avoid embolisms, the degassed water is used for Kleaf measurement in some previous studies [22, 39]. However, no study compared the Kleaf estimations using degassed water and non-degassed water. Our results showed that degassing had limited influence on Kleaf estimation, which indicates that gas dissolved in water does not necessarily causing xylem embolism. Unfortunately, no direct evidence has been detected in this study and in previous studies. The EFM follows the natural transpiration-driven water movement pathway in leaves, and the water enters the leaf through internal transpiration driving rather than through external water gravity pressure which may be caused by lower position of leaves than the water meniscus in the cylinder. However, the relative position difference between leaves and water source was generally ignored in previous studies [28, 29]. In this study, the estimated Kleaf showed no difference when the leaf was placed below, as high as, or above the meniscus of water in the cylinder. The result indicated that the small height difference between leaves and water source was permissible for Kleaf measurement. Indeed, only 0.196 Pa extra gravity pressure to leaf was changed in the study (Fig. 3).

Flow rate stabilization identification and water potential measurement

In previous studies, the flow rate stabilization was generally claimed a coefficient of variation < 5% over 3–10 min [22, 28, 40], but no study provided the flow rate variation over time. Apparently, this criterion was not suitable in our Kleaf measurement. Stricter criteria (P > 0.001, 15 min) for identifying steady state of flow rate were applied in this study. Unexpectedly, we found that the flow rate dramatically oscillated or even reduced after a short-term increase in the measurement of some leaves (Fig. 4). Stomatal oscillations have frequently been observed in previous studies due to the hydraulic mismatch between leaf water transpiration and water supply by xylem [41]. The simultaneous estimation of liquid flow rate and transpiration rate supported this mismatching hypothesis (Additional file 1: Figs. S5, S7), and the oscillation was frequent in high transpiration condition (such as 37 °C condition in Additional file 1: Fig. S6). The declined flow rate after reaching a peak might result from disruption of ionic homeostasis, which might be further attributed to the wide use of deionized water in EFM. Bundle sheath cells permeability have been reported to be related to xylem pH, ionic concentration and component [42,43,44]. Since the flow rate is highly dynamic over the measurement, and the identification of steady-state flow rate influences Kleaf estimation, reporting the details of flow rate estimation as well as the raw data of flow rate will be helpful for researchers to interpret the results and syntheses data for meta-analysis in the future.

Besides E, the Kleaf was also strongly affected by ψfinal. In this study, the ψfinal was generally higher than − 0.5 MPa, and, for a leaf with such a high leaf potential, a small error in ψfinal estimation would result a large change in Kleaf, emphasizing the importance of accuracy ψfinal measurement (Fig. 6). Water potential is typically measured using pressure chamber technique, and the methodological artifacts have been discussed for decades [45,46,47,48,49]. For instance, Levin [46] reported that the contrasting results were obtained by different operators. In addition, it is suggested that leaves need to be equilibrated in bags before the water potential measurement. However, the effects of equilibration time have rarely been reported in literature [46]. In this study, the O. sativa leaves sampled under dark and light had different ψfinal, but leaves sampled under both conditions rapidly achieved water potential equilibrium in bags (Additional file 1: Fig. S8). Further works on improving the accuracy of water potential estimation are needed.

Interestingly, the liquid flow rate and transpiration rate were equal when the entire leaves were clamped in the gas exchange chamber (Fig. 5). The shifted correlations of two-phase water flow rates for partly clamped leaves may be caused by the heterogeneity of transpiration along leaf blade [50]. Our result reminds us of the cautious use of in situ Kleaf measurement with photosynthetic instruments, another widely used method [27, 51].


We investigated the potential methodological artifacts in Kleaf estimation using EFM and showed that environmental factors, such as PARa, air temperature and airflow around leaf, identifications of steady-state flow rate, and ψfinal significantly affected Kleaf estimation. It is important to consider the environmental settings, the flow rate stabilization and precise water potential measurement when estimating Kleaf. In parallel, providing the details of the measurements is also necessary with greater expectations for data preservation, reproducible and open research [52, 53]. We recommend a table-like format (Additional file 3: Table S2) convenient for Kleaf data measurement and storage.


Plant materials

All the plant materials grew outdoors in Huazhong Agricultural University, Wuhan, China (114o22′E, 30o29′N). A monocot species, Oryza sativa L., cv Huanghuazhan (HHZ), was selected and the O. sativa plants were growing in paddy field for 50–70 days before sampling. O. sativa plants were well watered and fertilized, free of diseases, pests, and weeds. Meanwhile, a dicot species, Cinnamomum camphora L. was selected on campus of Huazhong Agricultural University.

Harvest time and sample storage

Oryza sativa tillers were cut off under water in the early morning (between 5:30 and 6:00 am) or the previous night (between 18:30 and 19:00 pm) of the measurement day. The fresh cuts of tillers were soaked in ultra-pure water, and tillers were covered by double black plastic bags. The samples collected at night before the measurement day were conserved in ordinary ultra-pure water and sterile ultra-pure water, and half the tillers were randomly selected and recut in the morning of the measurement day to test the effect of blockage at the cut surface (Additional file 1: Fig. S3). Longer than 0.5 m C. camphora branches were also sampled in the early morning or previous night and conserved in ordinary ultra-pure water. It took 10–15 min to transfer the samples to the laboratory. Branches and tillers were recut under ultra-pure water in the laboratory. Their cut ends were soaked in water, and other parts were covered with double black plastic bags at least 1 h.

Equipment settings

Kleaf was determined using evaporative flux method (EFM) reported previously [5, 22]. To minimize estimation errors, the water evaporation in the graduated cylinder without leaf under multiple conditions was quantified. The system transpiration was measured under the following four conditions: water without intervene, water surface covered by liquid wax, maintaining high humidity in the balance chamber by putting wet tissue papers, and the combination of liquid wax cover and putting wet tissue papers in the balance chamber. Based on our results (Additional file 1: Fig. S2), the combination of liquid wax cover and putting wet tissue papers in the balance chamber was adopted in the subsequent experiments due to its superior capacity to prevent water loss. In order to avoid ion deposition in leaf, ultra-pure water rather than ionic solution was adopted for leaf uptake[54, 55]. Ultra-pure water was vacuumed for 8 h to remove bubbles or directly stored overnight for Kleaf measuring of two species.

One end of low-resistance transparent tube (inner diameter = 2 mm, Oupli campany, Shanghai, China) filled with water was connected to the graduated cylinder containing with water on a balance (± 0.01 mg; Mettler MS205DU, Mettler-Toledo GmbH, Greifensee, Switzerland) (see the equipment diagram in Additional file 1: Fig. S9). Finally, water volume in cylinder was adjusted to ensure that leaves were placed 2 cm below the meniscus of the water in the cylinder for Kleaf measurement of two species.

Changes in environmental factors

To investigate the influences of environmental factors on Kleaf estimation, three environmental factors—air temperature, ambient photosynthetically active radiation (PARa), and airflow around O. sativa leaf—were individually changed. The environment conditions were as follows: the temperature was set as 37 °C\(\pm 1\) or 27 °C\(\pm 1\); PARa at leaf surface was set as 500, 1000, or 1500 µmol m−2 s−1; the airflow around leaves was set as 1 and 0 m s−1.

K leaf measurement

The newly- and fully-expanded leaves with 2 cm sheathes or petioles were cut from tiller under distilled water. The petioles of C. camphora leaves were connected to the water pipe using a hose tape. A hose tape and a cork were used to achieve seamless connection between O. sativa leaf sheathes and tube. Leaf was lifted higher than water surface to detect whether bubbles occurred in the connection. After leaf was placed on fish line net, leaf surface was wiped with tissue paper and irradiated by a lamp (600 W, Weichuang Company, Wuhan, China). A box fan (Comfort Zone 20 Inch Box Fan, the factory Depot Advantages, Inc, USA) was used to minimize the boundary resistance. At the same time, water weight in the graduated cylinder and the slope between weight and time were recorded every 3 s.

The water loss rate into the leaves was recorded until it was stable for a period of time (> 15 min). The detail identification of steady state was described in “Statistical analysis” section below. The temperature of the blade middle was determined as leaf temperature using a thermocouple (XimaAS877, Wanchuang electronic products Co., Ltd., Dongguan, China). Afterwards, leaf area and the final leaf water potential (ψfinal) were measured. Kleaf was calculated according to the following formula:

$${K}_{leaf}=E/(0-{\Psi }_{\mathrm{final}})$$

All the Kleaf values were normalized to those at 25 °C considering that water viscosity varied with temperature [56]. The measurements were performed from 8:00 am to 18:00 pm since there was no correlation between measuring time and Kleaf or E (data not shown).

Leaf water potentials

Upper and lower leaves adjacent to the target leaf used for hydraulic conductance measurement were cut from the tiller before Kleaf measurement, quickly put in an exhaled double-layer zip-lock bag, and placed in a foam box for water potential equilibration. Subsequently, leaf initial water potential (ψinitial) was detected in pressure chamber (PMS Instrument Company, Albany, OR, USA). Constant slow pressurization rate (< 0.05 MP s−1) was maintained during measurement. After flow measurement, the final leaf water potential (ψfinal) was determined as described above. To investigate the influences of equilibration time on leaf water potential estimation, water potential measurements were conducted on leaves in foam box for 10, 20, 30, 60, 90, 120, 150, and 180 min. Since there was no difference in ψfinal under 10–60 min equilibration, a 30 min of equilibration to leaves were applied in this study.

Liquid flow and gas flow comparison

In order to test the consistency of liquid flow rate and gas flow rate, the balance based liquid flow rate and Licor 6800 (LI-COR Inc., Lincoln, NE, USA) based gas flow rate were simultaneously measured. Leaves were quickly clamped into a 6 × 6 cm transparent chamber (Li-6800-13, LI-COR Inc., Lincoln, NE, USA) in the middle of O. sativa leaf or the entire leaf blade of the C. camphora after being put on net. In addition, a 3 × 3 cm transparent chamber (Li-6800-02) was used to clamp part C. camphora leaf blades. The chamber environment was set as coincident with ambient environment as possible. Auto-log was conducted with 30-s interval until a steady state was reached.

Statistical analysis

As the flow data was typically dynamic with time, the judgment on whether the flow rate has stabilized is challenging. An effective method to restrict segment lengths of given flow data is to explicitly allow high variance of segments, and segment length restriction was achieved via the break-point penalty parameter P in ‘dpseg’ package. In our analysis, high P value will allow high variance of the individual segments to produce long segments. The flow rate—time curve was segmented according to P-value, and the obtained segments were marked with different color. Different segments within 1 min separated by a few outliers were deemed invalid (Fig. 4). The time of last segment was required to be longer than 15 min, and the t-test P-value of the curve slope corresponding to the last 15 min (about 300 points) was required to be larger than 0.001.

One-way analysis of variation (ANOVA) and multiple comparisons (least significant ranges, ‘agricolae’ package) were conducted to test the significance of different treatments. The correlation in Fig. 5 was fitted using ‘ggpmisc’ package. All figures were plotted using ‘tidyverse’ package. All of the statistics and plotting were performed in R version 3.6.1 (

Availability of data and materials

The datasets supporting the findings of this study are available within the paper and within its additional files published online.


  1. Scoffoni C, Chatelet DS, Pasquet-kok J, et al. Hydraulic basis for the evolution of photosynthetic productivity. Nat Plants. 2016;2:16072.

    Article  CAS  PubMed  Google Scholar 

  2. Sack L, Holbrook NM. Leaf hydraulics. Annu Rev Plant Biol. 2006;57:361–81.

    Article  CAS  PubMed  Google Scholar 

  3. Brodribb TJ, Holbrook NM. Stomatal closure during leaf dehydration, correlation with other leaf physiological traits. Plant Physiol. 2003;132:2166–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Blackman CJ, Brodribb TJ. Two measures of leaf capacitance: insights into the water transport pathway and hydraulic conductance in leaves. Funct Plant Biol. 2011;38:118–26.

    Article  PubMed  Google Scholar 

  5. Sack L, Melcher PJ, Zwieniecki MA, Holbrook NM. The hydraulic conductance of the angiosperm leaf lamina: a comparison of three measurement methods. J Exp Bot. 2002;53:2177–84.

    Article  CAS  PubMed  Google Scholar 

  6. Nardini A, Tyree MT, Salleo S. Xylem cavitation in the leaf of Prunus laurocerasus and its impact on leaf hydraulics. Plant Physiol. 2001;125:1700–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Cochard H, Venisse JS, Barigah TS, et al. Putative role of aquaporins in variable hydraulic conductance of leaves in response to light. Plant Physiol. 2007;143:122–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Guyot G, Scoffoni C, Sack L. Combined impacts of irradiance and dehydration on leaf hydraulic conductance: insights into vulnerability and stomatal control. Plant, Cell Environ. 2012;35:857–71.

    Article  CAS  Google Scholar 

  9. Xiong D, Flexas J. Leaf economics spectrum in rice: Leaf anatomical, biochemical and physiological trait trade-offs. J Exp Bot. 2018;69:5599–609.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Wright IJ, Reich PB, Westoby M, et al. The worldwide leaf economics spectrum. Nature. 2004;428:821–7.

    Article  CAS  PubMed  Google Scholar 

  11. Sack L, Scoffoni C, John GP, et al. How do leaf veins influence the worldwide leaf economic spectrum? Review and synthesis. J Exp Bot. 2013;64:4053–80.

    Article  CAS  PubMed  Google Scholar 

  12. Sade N, Shatil-Cohen A, Attia Z, et al. The role of plasma membrane aquaporins in regulating the bundle sheath-mesophyll continuum and leaf hydraulics. Plant Physiol. 2014;166:1609–20.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Pantin F, Monnet F, Jannaud D, et al. The dual effect of abscisic acid on stomata. New Phytol. 2013;197:65–72.

    Article  CAS  PubMed  Google Scholar 

  14. Scoffoni C, Albuquerque C, Cochard H, et al. The causes of leaf hydraulic vulnerability and its influence on gas exchange in Arabidopsis thaliana. Plant Physiol. 2018;178:1584–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Attia Z, Dalla A, Moshelion M. Vascular bundle sheath and mesophyll cells modulate leaf water balance in response to chitin. Plant J. 2020;101:1368–77.

    Article  CAS  PubMed  Google Scholar 

  16. Yang Y, Zhang Q, Huang G, Peng S, Li Y. Temperature responses of photosynthesis and leaf hydraulic conductance in rice and wheat. Plant Cell Environ. 2020;43:1437–51.

    Article  CAS  PubMed  Google Scholar 

  17. Huang G, Zhang Q, Wei X, Peng S, Li Y. Nitrogen can alleviate the inhibition of photosynthesis caused by high temperature stress under both steady-state and flecked irradiance. Front Plant Sci. 2017;8:945.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Xiong D, Yu T, Zhang T, et al. Leaf hydraulic conductance is coordinated with leaf morpho-anatomical traits and nitrogen status in the genus Oryza. J Exp Bot. 2015;66:741–8.

    Article  CAS  PubMed  Google Scholar 

  19. Ding L, Li Y, Gao L, et al. Aquaporin expression and water transport pathways inside leaves are affected by nitrogen supply through transpiration in rice plants. Int J Mol Sci. 2018;19:256.

    Article  PubMed Central  CAS  Google Scholar 

  20. Locke AM, Ort DR. Leaf hydraulic conductance declines in coordination with photosynthesis, transpiration and leaf water status as soybean leaves age regardless of soil moisture. J Exp Bot. 2014;65:6617–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Scoffoni C, Kunkle J, Pasquet-Kok J, et al. Light-induced plasticity in leaf hydraulics, venation, anatomy, and gas exchange in ecologically diverse Hawaiian lobeliads. New Phytol. 2015;207:43–58.

    Article  CAS  PubMed  Google Scholar 

  22. Sack L, Scoffoni C. Measurement of leaf hydraulic conductance and stomatal conductance and their responses to irradiance and dehydration using the evaporative flux method (EFM). J Vis Exp. 2012;70:e4179.

    Google Scholar 

  23. Ounapuu-Pikas E, Sellin A. Plasticity and light sensitivity of leaf hydraulic conductance to fast changes in irradiance in common hazel (Corylus avellana L.). Plant Sci. 2020;290:110299.

    Article  CAS  PubMed  Google Scholar 

  24. Scoffoni C, Rawls M, McKown A, Cochard H, Sack L. Decline of leaf hydraulic conductance with dehydration: relationship to leaf size and venation architecture. Plant Physiol. 2011;156:832–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Trifilo P, Petruzzellis F, Abate E, Nardini A. The extra-vascular water pathway regulates dynamic leaf hydraulic decline and recovery in Populus nigra. Physiol Plant. 2021;172:29–40.

    Article  CAS  PubMed  Google Scholar 

  26. Xiong D, Douthe C, Flexas J. Differential coordination of stomatal conductance, mesophyll conductance, and leaf hydraulic conductance in response to changing light across species. Plant Cell Environ. 2018;41:436–50.

    Article  CAS  PubMed  Google Scholar 

  27. Sack L, Cowan PD, Holbrook NM. The major veins of mesomorphic leaves revisited: tests for conductive overload in Acer saccharum (Aceraceae) and Quercus rubra (Fagaceae). Am J Bot. 2003;90:32–9.

    Article  PubMed  Google Scholar 

  28. Brocious CA, Hacke UG. Stomatal conductance scales with petiole xylem traits in Populus genotypes. Funct Plant Biol. 2016;43:553–62.

    Article  CAS  PubMed  Google Scholar 

  29. Wu T, Tissue DT, Li X, et al. Long-term effects of 7-year warming experiment in the field on leaf hydraulic and economic traits of subtropical tree species. Globe Chang Biol. 2020;26:7144–57.

    Article  Google Scholar 

  30. Stiller V, Lafitte HR, Sperry JS. Hydraulic properties of rice and the response of gas exchange to water stress. Plant Physiol. 2003;132:1698–706.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Locke AM, Ort DR. Diurnal depression in leaf hydraulic conductance at ambient and elevated [CO2] reveals anisohydric water management in field-grown soybean and possible involvement of aquaporins. Environ Exp Bot. 2015;116:39–46.

    Article  CAS  Google Scholar 

  32. Stiller V, Sperry JS, Lafitte R. Embolized conduits of rice (Oryza sativa, Poaceae) refill despite negative xylem pressure. Am J Bot. 2005;92:1970–4.

    Article  PubMed  Google Scholar 

  33. Savatin DV, Gramegna G, Modesti V, Cervone F. Wounding in the plant tissue: the defense of a dangerous passage. Front Plant Sci. 2014;5:470.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Jacobs AK, Lipka V, Burton RA, et al. An Arabidopsis callose synthase, GSL5, is required for wound and papillary callose formation. Plant Cell. 2003;15:2503–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Clarke HRG, Lawrence SD, Flaskerud J, Korhnak TE, Gordon MP, Davis JM. Chitinase accumulates systemically in wounded poplar trees. Physiol Plant. 1998;103:154–61.

    Article  CAS  Google Scholar 

  36. Heil M. Nightshade wound secretion: the world’s simplest extrafloral nectar? Trends Plant Sci. 2016;21:637–8.

    Article  CAS  PubMed  Google Scholar 

  37. Bellincampi D, Cervone F, Lionetti V. Plant cell wall dynamics and wall-related susceptibility in plant–pathogen interactions. Front Plant Sci. 2014; 5.

  38. Scoffoni C, Albuquerque C, Brodersen C, et al. Outside-xylem vulnerability, not xylem embolism, controls leaf hydraulic decline during dehydration. Plant Physiol. 2017;173:1197–210.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Fu X, Meinzer FC, Woodruff DR, et al. Coordination and trade-offs between leaf and stem hydraulic traits and stomatal regulation along a spectrum of isohydry to anisohydry. Plant Cell Environ. 2019;42:2245–58.

    Article  CAS  PubMed  Google Scholar 

  40. Brodribb TJ, Feild TS, Jordan GJ. Leaf maximum photosynthetic rate and venation are linked by hydraulics. Plant Physiol. 2007;144:1890–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Farquhar GD, Cowan IR. Oscillations in stomatal conductance: the influence of environmental gain. Plant Physiol. 1974;54:769–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Geilfus CM. The pH of the apoplast: dynamic factor with functional impact under stress. Mol Plant. 2017;10:1371–86.

    Article  CAS  PubMed  Google Scholar 

  43. Grunwald Y, Wigoda N, Sade N, et al. Arabidopsis leaf hydraulic conductance is regulated by xylem sap pH, controlled, in turn, by a P-type H(+) -ATPase of vascular bundle sheath cells. Plant J. 2021;106:301–13.

    Article  CAS  PubMed  Google Scholar 

  44. Secchi F, Pagliarani C, Cavalletto S, et al. Chemical inhibition of xylem cellular activity impedes the removal of drought-induced embolisms in poplar stems—new insights from micro-CT analysis. New Phytol. 2021;229:820–30.

    Article  CAS  PubMed  Google Scholar 

  45. Turner NC. Measurement of plant water status by the pressure chamber technique. Irrig Sci. 1988;9:289–308.

    Article  Google Scholar 

  46. Levin AD. Re-evaluating pressure chamber methods of water status determination in field-grown grapevine (Vitis spp.). Agric Water Manag. 2019;221:422–9.

    Article  Google Scholar 

  47. Hochberg U. Facilitating protocols while maintaining accuracy in grapevine pressure chamber measurements—comments on Levin 2019. Agric Water Manag. 2020;227: 105836.

    Article  Google Scholar 

  48. Boyer JS. Isopiestic technique: measurement of accurate leaf water potentials. Science. 1966;154:1459–60.

    Article  CAS  PubMed  Google Scholar 

  49. Scholander PF, Hammel HT, Bradstreet ED, Hemmingsen EA. Sap pressure in vascular plants. Science. 1965;148:339–46.

    Article  CAS  PubMed  Google Scholar 

  50. Xiong D, Yu T, Liu X, et al. Heterogeneity of photosynthesis within leaves is associated with alteration of leaf structural features and leaf N content per leaf area in rice. Funct Plant Biol. 2015;42:687–96.

    Article  CAS  PubMed  Google Scholar 

  51. Attia Z, Domec JC, Oren R, Way DA, Moshelion M. Growth and physiological responses of isohydric and anisohydric poplars to drought. J Exp Bot. 2015;66:4373–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Ely KS, Rogers A, Agarwal DA, et al. A reporting format for leaf-level gas exchange data and metadata. Ecol Informatics. 2021;61:101232.

    Article  Google Scholar 

  53. Kattge J, Bönisch G, Díaz S, et al. TRY plant trait database—enhanced coverage and open access. Glob Change Biol. 2019;26:119–88.

    Article  Google Scholar 

  54. Gasco A, Nardini A, Gortan E, Salleo S. Ion-mediated increase in the hydraulic conductivity of Laurel stems: role of pits and consequences for the impact of cavitation on water transport. Plant Cell Environ. 2006;29:1946–55.

    Article  CAS  PubMed  Google Scholar 

  55. Hu Y, Schmidhalter U. Spatial distributions and net deposition rates of mineral elements in the elongating wheat (Triticum aestivum L.) leaf under saline soil conditions. Planta. 1998;204:212–9.

    Article  CAS  Google Scholar 

  56. Yang S, Tyree MT. Hydraulic resistance in Acer saccharum shoots and its influence on leaf water potential and transpiration. Tree Physiol. 1993;12:231–42.

    Article  CAS  PubMed  Google Scholar 

Download references


Not applicable.


This study was supported by the National Natural Science Foundation of China (No. 32022060).

Author information

Authors and Affiliations



XW and JZ contributed to data collection; WX conducted data analysis and writing; JH and SP read and approved the manuscript; DX contributed to conceptualization, data analysis and writing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Dongliang Xiong.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

A literature survey on Kleaf measuring methods. Figure S2. The water loss from the cylinder without leaf under different preventions. Figure S3. Influences of sample storage time and recutting on initial leaf water potential. Figure S4. Dynamic water flow rate measured by a Licor 6800 and a balance in the same leaf. Figure S5. Oscillation of water flow rate, transpiration rate, and stomatal conductance. Figure S6. Temperature effects on water flow rate estimation. Figure S7. Modeled final leaf water potential and Kleaf changes over the time. Figure S8. Effects of equilibration time on leaf water potential estimation. Figure S9. Experimental setup used in determining Kleaf. Figure S10. Meteorological data of rice growing season.

Additional file 2: Table S1.

The meta-data sheet of reported Kleaf.

Additional file 3: Table S2.

All tidy data used in the main figures.

Additional file 4: Table S3.

All raw flow data used in the in the main figures.

Rights and permissions

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, X., Zhao, J., Huang, J. et al. Evaporative flux method of leaf hydraulic conductance estimation: sources of uncertainty and reporting format recommendation. Plant Methods 18, 63 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: