Free Access
Issue
Ann. For. Sci.
Volume 66, Number 6, September 2009
Article Number 608
Number of page(s) 7
DOI https://doi.org/10.1051/forest/2009043
Published online 31 July 2009

© INRA, EDP Sciences, 2009

1. INTRODUCTION

Sap flow measurement is an important technique for studying in situ plant-water relations. The method originally devised by Granier (1985; 1987), based on a thermal dissipation probe, is one of the most frequently used. The probe consists of two cylindrical needles inserted radially into the stem. One needle, placed above the other, contains a heating element, supplied at constant power by Joule effect. The heating needle exchanges sensible heat with its environment made of wood and sap. The physical process responsible for the transfer of heat from the needle to the wood involves both conduction in a solid medium (xylem vessels) and convection in a fluid in motion (sap flow). In stationary regime, a relationship linking sap flow density with the temperature difference between the probe needles can be derived and constitutes the basis of Granier’s method of sap flow measurement. This relationship apparently does not depend on tree species when sensors geometry (e.g. effective length) and heating power remain the same, which was confirmed by other studies (Lu et al., 2004). The sap flow for the whole plant is obtained by multiplying the sap flow density by the sapwood area at the point of measurement, which is not always easy to determine accurately. Due to its simplicity, low cost and reliability, Granier’s stationary method has been widely used, particularly in forestry (Andrade et al., 1998; Braun and Schmid, 1999; Lüttschwager and Remus, 2007; Ma et al., 2008; Masmoudi et al., 2004; Meinzer et al., 1999; Oliveras and Llorens, 2001; Sevanto et al., 2008; Wilson et al., 2001; Wullschleger and Hanson, 2008). However, substantial errors on sap flow estimates are induced by heat storage, especially in small stems, and by natural thermal gradients in the wood (Cabibel and Do, 1991; Do and Rocheteau, 2001a). Additionally, the zero flow conditions needed for sap flow density calculation (often observed at predawn) may not occur under certain circumstances, such as during the slow restoration of internal tree water storage following a prolonged drought period (Lu et al., 2004).

Do and Rocheteau (2002b) use a cyclic heating instead of the continuous heating in order to correct sap flow estimates for natural thermal gradients. They measure probe temperature difference at the end of heating periods and subtract from it the probe temperature difference at the end of cooling periods. When cyclic heating is used, the stationary regime does not occur except for high flow densities, which requires a specific calibration to account for the non-stationary regime. The cyclic heating substantially improves the accuracy of sap flow measurement in situations where thermal gradients are important (Lu et al., 2004). Yet, the determination of zero flow conditions is still required. Do and Rocheteau (2002b) pointed out that the measurement principle of the cyclic heating would apply to a single heating needle probe. Our work follows this statement, but it analyses the single probe temperature occurring just after current is switched on or off, in order to develop the basis for a new approach of measurement. Specifically, it explores the pattern of temperature variation during the transient regime and how it is correlated to flow density.

2. MATERIALS AND METHODS

2.1. Theoretical basis

When the thermal dissipation probe is being heated by Joule effect through a given resistance and with a given current intensity, the temperature of the probe reaches a steady-state regime characterized by an equilibrium temperature. The theoretical basis of this stationary regime was first developed by Granier (1985) and then completed by Valancogne and Granier (1991). In this theory the heat flux φ transferred from the probe to the wood is written as

(1)
where T is the temperature of the probe, Tw is the average temperature of the wood around the probe (the temperature field is assumed to be homogeneous) and h is the exchange coefficient, expressed as a function of the flow density u as
(2)
where α and β are two empirical coefficients (β is assumed to be equal to 1 in Granier’s original paper) and h0 is the value of the exchange coefficient at zero density (hence, representing the conductive part of the heat transfer).

When heating current is switched off, the heat lost by the probe is dissipated within the medium and its temperature Tdecreases more or less rapidly down to the temperature of the wood, satisfying the transient heat balance equation written in the form of the following differential equation

(3)
where C is the heat capacity of the probe and t is time. Assuming Twto be constant during the integration time, the solution of this first-order linear differential equation is
(4)
where T0 is the temperature of the probe at the initial time of the kinetic (t = 0, when the current is switched off). The temperature of the wood Tw represents the probe temperature at equilibrium (when t → + ∞). According to this theory, probe temperature Tfollows a pseudo-exponential decay with time, which depends upon the heat capacity of the probe C and the coefficient of heat exchange between the probe and the wood h (u). If u is assumed to be constant during the integration time, the coefficient of the exponential reduces to − h (u) t / C and the decay becomes strictly exponential with time. Inverting equation (4) leads to the following expression
(5)
Taking into account equation (2) with β = 1, the integral in equation (5) can be expressed as a function of the mean sap flow density (between 0 and t)
(6)
This means that for a given value of flow density, the function y (t) is a straight line. Its slope varies with flow density: the greater the flow density, the higher the slope of the line. For any intermediate time ti an index I is defined as
(7)
where T (ti) is the temperature of the probe at t = ti. Equations (5) and (6) show that there is a one to one relationship between the mean flow density (between 0 and ti) and the index I (ti). This relationship is expressed as
(8)
where a = C / (αh0) and b = − 1/ αare constant coefficients. Simple linear relationships of the form of equation (8) between mean sap flow density and thermal index I (ti) will be experimentally calibrated. If we consider the heating phase, starting just when the current is switched on, the same reasoning holds, but equation (7) should be written as (see appendix)
(9)
where Te is the equilibrium temperature of the probe when the stationary regime is reached.

thumbnail Figure 1

Schematic view of the testing bench used in the experiment (1): compressed nitrogen tank; (2): pressure regulator; (3): water reservoir; (4): pressure gauge; (5): olive tree branch; (6): TDP; (7): container; (8): balance.

2.2. Experimental setting

An experiment was carried out to test the theory developed above and to calibrate a relationship between flow density and the thermal indices defined above. The basic idea of the method is to record, with a given time step, the temperature of the probe just after the heating current is switched on or switched off. A thermal dissipation probe (TDP) based on the original Granier’s method was designed. This probe consists of a single needle with a heating element and a T-type thermocouple. Cold junction compensation is measured by a temperature sensor, NTC type (Negative Temperature Coefficient) protected by an isothermal shield.

A testing bench (Fig. 1) was used to carry out the experiment. The heating probe was inserted radially into an olive tree (Olea europaea L.) branch of about 5 cm diameter. The bark was removed near the probe so that the effective probe length (2 cm) was entirely located in the sapwood. Olive wood is classified as a diffuse porous xylem where vessels are distributed quite uniformly in the annual growth rings. Water was injected into the olive tree branch from a reservoir maintained under pressure by means of compressed nitrogen, the flow rate being set by means of a pressure regulator. The water flowing out of the branch is collected into a container and weighed by a digital balance. A data logger (Delta-T devices, Cambridge, UK) was used to control the heating and to record the signals from the thermocouples, thermistor and digital balance. A relatively high frequency of temperature measurement is adopted because the temperature variation is fairly rapid for high flow densities. Measurement sequences using 5 s time step are taken during a period of 10 min just after heating current is switched on or off. Measurements were made for twelve levels of flow density from 0 to 9 L h−1 dm−2 covering the whole range of sap flow, as found in field experimentson olive trees.

Sapwood cross area was determined by injecting water containing red-dye safranin through the olive tree branch. As almost all the wood area was colored in fuchsia, we considered that the entire section was conductive.

thumbnail Figure 2

Kinetics of probe temperature for different values of flow density (u) expressed in L h− 1dm− 2: (a) during the heating phase and (b) during the cooling phase. Ambient temperature was not constant during the measurements.

3. RESULTS

3.1. Characteristics of the kinetics

Figure 2 shows the variation of probe temperature during the heating phase and the cooling phase for different values of flow density. For u= 3 L dm−2 h−1 and u= 7 L dm−2 h−1 and for both heating and cooling periods, an equilibrium state is reached within 10 min and a plateau is clearly visible. The time needed for the kinetic to be completed decreases when flow density increases. For low densities no clear plateau can be noticed within the 10 min. It has been found that the difference between probe temperature measured after 10 and 20 min of cooling does not exceed 2% of the temperature recorded at 10 min for a wide range of flow density. Once the kinetic is completed, temperature remains approximately constant, if there is no change in the environmental conditions. This temperature depends on the level of flow density. Equation (A.2) in the appendix shows that for a given heating power, the higher the flow density (and thus the heat exchange coefficient h (u)), the lower the equilibrium temperature Te. During the cooling phase, the stationary temperature Tw adjusts to ambient temperature variations (Fig. 2b).

In fact, the probe temperature was found not to fit the strict exponential variation with time that could be expected from the theory developed above assuming flow density to remain constant. In similar experiments involving non-stationary regimes, Benet et al. (1977) reported that the exponential decrease occurs only after a short period of time, during which the variation of temperature is faster than the exponential variation predicted by the theory, for a reason not completely elucidated yet. This behavior was verified in our experiment: a rapid cooling or heating is observed during a short period, and then the temperature variation slows down, the experimental points approximately fitting an exponential curve. The delay needed for the exponential variation to appear depends on flow density level. The duration of the exponential phase also varies with flow density. The theory developed above (Eq. (5)) demonstrates that under a constant flow density, and thus a constant heat exchange coefficient, y (t) = ln [(T0Tw) / (T (t) − Tw)] should be a linear function of time during the cooling phase (Tw being approximated by the probe temperature 10 min after switch off). As shown in Figure 3, where y (t) is plotted against time, the strict exponential phase, which corresponds here to the linear part of the curve, is in fact delayed with respect to the switch off time. The delay before the linear phase is reached and the duration of this linear phase decrease with increasing flow density. It appears also from Figure 3 that the greater the flow density, the higher the slope of the line.

thumbnail Figure 3

Linear part of function y (t) defined by equation (5) versus time during the cooling phase of the probe for different levels of flow density u (L h− 1 dm− 2).

3.2. Calibration and validation

The experimental data set corresponding to several repetitions of measurement sequences of heating and cooling was split into two subsets, one (A) for calibration and the other (B) for validation. The consecutive cycles were alternatively put in subset A and B. Thermal indices I (ti) were calculated for tiranging from 10 to 60 s every 10 s, Tw and Te being taken at t = 600s. The temperature recorded after 10 min was used as an approximation of the probe temperature at equilibrium for both heating and cooling periods. Correlations between thermal indices I (ti) and gravimetrically measured water flow densities (ume) were determined for both heating and cooling periods with data subset A. The coefficients of the regression lines ume = aI (ti) + bwith the corresponding coefficient of determination (R2 ) are given in Table I. First, it appears that the proposed thermal indices I (ti) are generally well correlated to measured water flow density with R2 values in the range 0.83–0.96. Second, the linear regressions are better for cooling mode than for heating mode. Third, the highest R2 values are obtained for thermal indices I (ti) calculated at ti = 10 s, 20 s and 30 s during the cooling phase, the corresponding regression lines are shown in Figure 4. Beyond ti = 60s, the R2 values decrease substantially.

Table I

Calibration equations u = aI (ti) + b between the flow density u and the thermal index I for ti ranging from10 s to 60 s. R2 is the coefficient of determination of the regression line (17 and 13 observations respectively for heating and cooling modes).

thumbnail Figure 4

Regression lines of the gravimetrically measured water flow density ume versus I (ti) for ti = 10, 20 and 30 s during the cooling phase (data set A).

thumbnail Figure 5

Root mean square errors (RMSE) and bias (B) between the gravimetrically measured water flow densities and values estimated by equation (8) using I (ti) in the cooling phase for ti ranging from 10 s to 60 s.

thumbnail Figure 6

Gravimetrically measured water flow densities versus estimated values using I (ti) in the cooling phase with the intermediate time ti = 20 s (data set B).

In the validation work, the coefficients of the linear regression u = aI (ti) + bgiven in Table I for the cooling phase were used to estimate the values of flow density (ues) with the independent data subset B. In order to select the most appropriate value of ti, the root mean square errors (RMSE) and the bias (B) between ues (estimated flow density) and ume

(gravimetrically measured flow density) were calculated for ti varying form 10 to 60 s (Fig. 5). The lowest RMSE is observed for ti = 20s with an acceptable bias, suggesting that it is the most appropriate time, at least for our experiment. Figure 6 compares the measured values of water flow density with those estimated using the regression equation corresponding to ti = 20 s. With a flow density in the range 0–7 L h−1 dm−2 the RMSE value is 0.46 L h−1 dm−2 and the bias is 0.14 L h−1 dm−2 which is rather satisfactory. For ti = 10 s and ti = 30s, RMSE values are also acceptable with respectively 0.50 L h−1 dm−2 and 0.52 L h−1 dm−2 . These results confirm the robustness of this new approach based upon the transient regime of a single probe, which can be a reliable method of sap flow measurement.

4. DISCUSSION

The techniques of sap flow measurement by thermal dissipation probe with constant heating (Granier, 1985) and cyclic heating (Do and Rocheteau, 2002b) have been revisited. It has been shown that the transient regime of the single probe occurring just after the heating current is switched off allows the definition of a thermal index well correlated to flow density. This index is expressed as a combination of three temperatures of the probe measured at three different times of the kinetic: initial temperature T0, final temperature Tw and an intermediate temperature T (ti) taken between ti = 10 s and 30 s. This thermal index is the basis of a transient approach of sap flow measurement by a single probe.

When compared to the previous methods, this approach appears to have several advantages. The first practical advantage is its simplicity, since only one probe is needed with the same complexity in electronic control as the cyclic method of Do and Rocheteau (2002b). Second, the problems linked to thermal gradients within the sapwood, when constant heating is used, are eliminated since the temperature is measured only at one location. Additionally, the method does not need, in principle, the determination of zero flow conditions, which are not easy to identify practically, and which are required in Granier’s method (to eliminate the heating power RI2 in the equations) as well as in Do and Rocheteau’s one. Nevertheless, if the thermal index in zero flow conditions I0 (ti) is known, the calibration relationship (represented by equation (8)) simply transforms into u = − bK, where K = [I (ti) / I0 (ti)] − 1 is a dimensionless index similar to the one used in Granier’s method. However, this new relationship, in which the heat exchange coefficient at zero flow (h0) is eliminated, is not fundamentally different from equation (8) and apparently not more accurate. An inherent problem with the method of Do and Rocheteau (2001b) is the non-steady state nature of the sap flow during the measurement period. Indeed, in periods of rapid changing meteorological conditions (early morning, late afternoon or cloudy conditions), steady flow rarely happens. It is not clear how this difficulty is resolved in their method. This problem apparently does not exist in our approach. The theory shows that it is the mean flow between t = 0and the intermediate time ti, which is taken into account (a period of time lower than 30 s in principle). The period after t = ti, is necessary to allow for the probe to reach the wood temperature Tw which should remain constant.

In our experiment, measurement periods were arbitrarily set to 10 min in order to approximate stationary wood temperature. For both cooling and heating modes, additional 10 min were used respectively for heating and cooling between consecutive measurement sequences. However, as the initial temperature T0 involved in the thermal index for cooling mode is not necessarily at equilibrium, a shorter heating period could be adopted and would result in a faster cooling. The number of sap flow measurements per hour could therefore be increased. For the selection of cooling duration, a compromise should be found: if the cooling phase is too short, the probe will not approach the wood equilibrium temperature (Tw); if it is too long, flow density and external conditions may change and Tw will no longer be related to the initial phase (ti). Additional experiments will be carried out to confirm the possibility of shortening the measurement sequences and further field work will test the single probe transient method for measuring tree sap flow. It will be also relevant to check if the calibration relationships obtained with an olive tree branch are independent of the medium used, as for Granier’s relationship.

Appendix: Theoretical basis of the transient method in the heating phase

When the heating phase is considered, starting when the Joule effect (through a resistance R and a current intensity I) is switched on, the transient heat balance equation is written as

(A1)
with the same symbols as those used in the main text. When the steady-state regime is reached, assuming all the heat produced by the probe to be exchanged with the wood, the equilibrium temperature of the probe is denoted by Te and equation (A.1) transforms into
(A2)
This means that equation (A.1) can be rewritten as
(A3)
The solution of which is
(A4)
where T0 is the temperature of the probe at t = 0and Te is the equilibrium temperature of the probe when the stationary regime is reached. Consequently, an equation very similar to equation (5) can be written as
(A5)
which leads to equation (9).

List of symbols

C: heat capacity of the probe

h: exchange coefficient of the probe

h0: exchange coefficient of the probe at zero flow

I: thermal index defined by equation (7) or (9)

T (t) : temperature of the probe at time t

Te: temperature of the probe at equilibrium with constant heating

ti: intermediate time of the kinetic

T0: initial temperature of the probe when current is switched off (or on)

T w: temperature of the wood

u: sap flow density

α: parameter involved in the exchange coefficient hof the probe(Eq. (2))

β: parameter involved in the exchange coefficient hof the probe (Eq. (2))

References

  1. Andrade J.A., Meinzer F.C., Goldstein G., Holbrook N.M., Cavelier J., Jackson P., and Silvera K., 1998. Regulation of water flux through trunks, branches and leaves in trees of a lowland tropical forest. Oecologia 115: 463–471 [CrossRef].
  2. Benet J.C., Pomade P., and Hinzelin P., 1977. Mesure de la conductivité thermique des sols par sonde cylindrique en régime transitoire. Etude de l'influence du mode de mise en place de la sonde. Matér. constr. 10: 385–392.
  3. Braun P. and Schmid J., 1999. Sap flow measurements in grapevines (Vitis vinifera L.). 2. Granier measurements. Plant Soil 215: 47–55.
  4. Cabibel B. and Do F., 1991. Mesures thermiques des flux de sève dans les troncs et les racines et fonctionnement hydrique des arbres. I. Analyse théorique des erreurs sur la mesure des flux et validation des mesures en présence de gradients thermiques extérieurs. Agronomie 11: 669–678.
  5. Do F. and Rocheteau A., 2002a. A cyclic thermal dissipation system for measuring sap flow under high natural wood temperature gradients. I. Gradients and signals in the field. Tree Physiol. 22: 641–648.
  6. Do F. and Rocheteau A., 2002b. A cyclic thermal dissipation system for measuring sap flow under high natural wood temperature gradients. II. Analysis in an artificial flow system. Tree Physiol. 22: 649–654.
  7. Granier A., 1985. Une nouvelle méthode pour la mesure du flux de sève brute dans le tronc des arbres. Ann. Sci. For. 42: 193–200 [CrossRef].
  8. Granier A., 1987. Evaluation of transpiration in a Douglas-fir stand by means of sap flow measurements. Tree Physiol. 3: 309–320 [PubMed].
  9. Lu P., Urban L., and Zaho P., 2004. Granier's Thermal dissipation probe (TDP) method for measuring sap flow in trees: theory and practice. Acta Bot. Sin. 46: 631–646.
  10. Lüttschwager D. and Remus R., 2007. Radial distribution of sap flux density in trunks of a mature beech stand. Ann. For. Sci. 46: 431–438 [CrossRef].
  11. Ma L., Lu P., Zhao P., Rao X., Cai X., and Zeng X., 2008. Diurnal, daily, seasonal and annual patterns of sap-flux-scaled transpiration from an Acacia mangium plantation in southern China. Ann. For. Sci. 65: 402 [EDP Sciences].
  12. Marshall D.C., 1958. Measurement of sap flow in conifers by heat transport. Plant Physiol. 6: 385–396 [CrossRef].
  13. Masmoudi M.M., Mahjoub I., Charfi-Masmoudi C., Abid-Karray J., and Ben Mechlia N., 2004. Mise au point d'un dispositif de mesure du flux de sève xylémique chez l'olivier. Rev. Rég. Arides ns: 242–251.
  14. Meinzer F.C., Goldstein G., Franco A.C., Bustamente M., Igler E., Jackson P., Caldas L., and Rundel P.W., 1999. Atmospheric and hydraulic limitations on transpiration in Brazilian cerrado woody species. Funct. Ecol. 13: 273–282 [CrossRef].
  15. Oliveras I. and Llorens P., 2001. Medium-term sap flux monitoring in a Scots pine stand: analysis of the operability of the heat dissipation method for hydrological purposes. Tree Physiol. 21: 473–480 [PubMed].
  16. Sevanto S., Nikinmaa E., Riikonen A., Daley M., Pettyjohn J.C., Mikkelsen T.N., Phillips N., and Holbrook N.M., 2008. Linking xylem diameter variation with sap flow measurements. Plant soil 305: 77–90 [CrossRef].
  17. Valancogne C. and Granier A., 1991. Intérêt des méthodes thermiques de mesure du flux de sève pour l'étude du bilan hydrique des savannes. In Soil Water Balance in the Sudano-Sahelian Zone, Proceedings of the Niamey Workshop, IAHS Publ. 199: 387–400.
  18. Wilson K.B., Hanson P.J., Mulholland P.J., Baldocch D.D., and Wullschleger S.D., 2001. A comparison of methods for determining forest evapotranspiration and its components: sap-flow, soil water budget, eddy covariance and catchment water balance. Agric. For. Meteorol. 106: 153–168 [CrossRef].
  19. Wullschleger S.D. and Hanson P.J., 2008. Sensitivity of canopy transpiration in an upland oak forest: evidence from a long-term field manipulation study. Glob. Chang. Biol. 12: 97–109 [CrossRef].

All Tables

Table I

Calibration equations u = aI (ti) + b between the flow density u and the thermal index I for ti ranging from10 s to 60 s. R2 is the coefficient of determination of the regression line (17 and 13 observations respectively for heating and cooling modes).

All Figures

thumbnail Figure 1

Schematic view of the testing bench used in the experiment (1): compressed nitrogen tank; (2): pressure regulator; (3): water reservoir; (4): pressure gauge; (5): olive tree branch; (6): TDP; (7): container; (8): balance.

In the text
thumbnail Figure 2

Kinetics of probe temperature for different values of flow density (u) expressed in L h− 1dm− 2: (a) during the heating phase and (b) during the cooling phase. Ambient temperature was not constant during the measurements.

In the text
thumbnail Figure 3

Linear part of function y (t) defined by equation (5) versus time during the cooling phase of the probe for different levels of flow density u (L h− 1 dm− 2).

In the text
thumbnail Figure 4

Regression lines of the gravimetrically measured water flow density ume versus I (ti) for ti = 10, 20 and 30 s during the cooling phase (data set A).

In the text
thumbnail Figure 5

Root mean square errors (RMSE) and bias (B) between the gravimetrically measured water flow densities and values estimated by equation (8) using I (ti) in the cooling phase for ti ranging from 10 s to 60 s.

In the text
thumbnail Figure 6

Gravimetrically measured water flow densities versus estimated values using I (ti) in the cooling phase with the intermediate time ti = 20 s (data set B).

In the text