The Dual Method Approach (DMA) Resolves Measurement Range Limitations of Heat Pulse Velocity Sap Flow Sensors

: Sap ﬂow, the movement of ﬂuid in the xylem of plants, is commonly measured with the heat pulse velocity ( V h ) family of methods. The observable range of V h in plants is ~ − 10 to ~+270 cm/h. However, most V h methods only measure a limited portion of this range, which restricts their utility. Previous research attempted to extend the range of V h methods, yet these approaches were analytically intensive or impractical to implement. The Dual Method Approach (DMA), which is derived from the optimal measurement ranges of two V h methods, the Tmax and the heat ratio method (HRM), also known as the “slow rates of ﬂow” method (SRFM), is proposed to measure the full range of sap ﬂow observable in plants. The DMA adopts an algorithm to dynamically choose the optimal V h measurement via the Tmax or HRM/SRFM. The DMA was tested by measuring sap ﬂux density ( J s ) on Tecoma capensis (Thunb.) Lindl., stems and comparing the results against J s measured gravimetrically. The DMA successfully measured the entire range of V h observed in the experiment from 0.020 to 168.578 cm/h, whereas the HRM/SRFM range was between 0.020 and 45.063 cm/h, and the Tmax range was between 2.049 cm/h and 168.578 cm/h. A linear regression of DMA J s against gravimetric J s found an R 2 of 0.918 and error of 1.2%, whereas the HRM had an R 2 of 0.458 and an error of 49.1%, and the Tmax had an R 2 of 0.826 and an error of 0.5%. Different methods to calculate sapwood thermal diffusivity ( k ) were also compared with the k Vand method showing better accuracy. This study demonstrates that the DMA can measure the entire range of V h in plants and improve the accuracy of sap ﬂow measurements.


Introduction
Sap flow is synonymous with water movement in plants, and its accurate measurement is critical for physiologists, hydrologists, modelers, and growers [1,2]. Sap flow is widely measured with thermometric sensors based on a variety of theoretical or empirical approaches. The heat pulse velocity (HPV) family of methods is a popular approach [3] which is derived from theory on thermal conduction and convection through porous materials [4]. Numerous methods have been derived from the conduction/convection equation, including the popular Tmax method [5,6] and heat ratio method (HRM) [4,7].
Although the accuracy of each method is debatable [8], it is universally accepted that the Tmax method is limited at slow flows, and the HRM is limited at fast flows [3,6,[8][9][10]. The measurement range limitations are caused by sensor design, thermal properties of xylem, and data logging or electronic noise [8]. The minimum heat velocity (V h ) typically detected by the Tmax is approximately 5 to 10 cm/h; whereas the maximum V h detected by the HRM is 15 to 45 cm/h [7][8][9][10][11][12][13].

Theoretical Background
The movement of heat through a porous medium can be determined via the thermal conductance and convection equation. Marshall [4] provided the thermal conductance and convection equation as: where ∆T is the change in temperature before and after the application of heat, q is heat input, k is thermal diffusivity, t is time (s), and x and y are the distance between the heat source and temperature sensor in the axial and tangential directions, respectively, and V h is heat velocity. The Sapflow+ method explicitly handles y; however, the other HPV methods ignore heat movement in the tangential direction by assigning y as zero [17].

Tmax Method
All HPV methods are derivatives of Equation (1). For example, the Tmax method, following Cohen et al. [5], is given as: where t m is the time taken (s) to reach maximum temperature following the start of a heat pulse. Equation (2) assumes an instantaneous heat input, which is not possible in practice. Kluitenberg and Ham [18] proposed a modified version of the Tmax method which explicitly accounts for the duration of heat input: where t 0 is the length (s) of the heat pulse.

Heat Ratio Method (HRM) or Slow Rates of Flow Method (SRFM)
Marshall [4] derived an additional method from Equation (1) where it was desired to measure zero and slow rates of sap flow. Marshall [4] proposed that two temperature sensors be installed, one above and one below the heater source at equal distance in a three-probe configuration, with the following equation: where ∆T d and ∆T u are temperature changes following a heat pulse in the downstream and upstream temperature sensors from the heater source. Marshall's SRFM was later renamed the "heat ratio method" (HRM) by Burgess et al. [7].

Thermal Diffusivity (k)
Thermal diffusivity (k) is an important parameter in heat pulse velocity calculations, and can be measured in situ or from sapwood samples [19]. Uncertainties in the calculation of k arise because of variation in sapwood moisture content and of the density of dry wood [20], which means it must be carefully measured. Cohen et al. [5] and Kluitenberg and Ham [18] proposed methods to measure k derived from the Tmax equation. However, this approach is not recommended as it requires deriving k when there is zero sap flow, and the Tmax method has poor resolution at slow to zero sap flow (see Table 1 and [6]). Table 1. The theoretical heat velocity (V h ) measurement range for the Tmax and heat ratio method (HRM), also known as the "slow rates of flow" method (SRFM), at varying sapwood thermal diffusivity (k). The V h values are based on a probe spacing (x) of 0.6 cm and a 3 s heat pulse (t 0 ). The minimum data from Tmax is based on a one-second measurement resolution and Equation (3). The theoretical SRFM calculation is based on a ∆T d /∆T u ratio of 20 and Equation (4). V m_critical was calculated via Equation (10). Instead, k is determined via [4,19,20]: where K is the thermal conductivity of sapwood (W/m/ • C), ρ is the basic density of fresh sapwood (kg/m 3 ), and c is the specific heat capacity of fresh sapwood (J/kg/ • C). Edwards and Warwick [21] defined c as: where w d and w f are the sapwood's dry and fresh weight (kg), respectively; c d and c w are the specific heat capacity of the dry wood matrix and sap solution, respectively. The parameters c d and c w are assumed to be constants with assigned values of 1200 and 4182 (J/kg/ • C), respectively [22]. There are two similar methods to determine K. The first method (k Hogg ) follows Hogg et al. [23], which was also labelled k_Burg by Looker et al. [20] who assigned the method to the later publication of Burgess et al. [7], and is determined by: where K w is the thermal conductivity of water (0.5984 W/m/ • C), K d is the thermal conductivity of dry sapwood, m c is sapwood moisture content (kg/kg), and ρ d and ρ w are the basic density (kg m −3 ) of dry sapwood and water, respectively. The ρ d value is found by dividing the sapwood's dry weight by fresh volume, and ρ w is a constant with a value of 1000. The second method (k Vand ) follows Vandegehuchte and Steppe [19], and is similar to Equation 8, though with the inclusion of a fiber saturation point (FSP) parameter: The parameter m c_FSP , or sapwood moisture content at the fiber saturation point, can be quantified via Barkas [24] or given the nominal value of 0.26 (26%) following Kollmann and Cote [25]. The parameter F v_FSP is calculated as Vandegehuchte and Steppe [19]: where ρ cw is cell wall density and is assumed to be equal to 1530 (kg/m 3 ).

The Dual Method Approach (DMA)
The DMA combines the optimal measurement ranges of Tmax and the SRFM to output a single V h value. The DMA works on an algorithmic basis, where a decision-making process is implemented to decide whether to use V h measured via Tmax or SRFM. Figure 1 is the algorithmic flow chart describing the DMA process. The first step is to complete a concurrent measurement cycle for Tmax and SRFM. The measurement cycle outputs V m , the measured heat velocity. If V m <SRFM max (the maximum V h via SRFM), then the measurement from the SRFM is used for V h ; if V m is >Tmax min (the minimum V h via Tmax), then the measurement from Tmax is adopted for V h .
Determining SRFM max and Tmax min is difficult, as it can vary with k and, possibly, other factors. A V m_critical value can be statistically determined via a regression of V h , as determined by the SRFM, against an independent measure of sap flow, such as a gravimetric weighing lysimeter. When correlated against an independent measure of sap flow, V h via the SRFM will approach a plateau. A piecewise linear regression analysis, with V h via the SRFM on the Y-axis and an independent measure of sap flow on the X-axis, can statistically determine a V m_critical value via the breakpoint. Values below the breakpoint are assigned to the SRFM and values greater than the breakpoint are assigned to the Tmax. However, this statistical approach is practical if a sap flow sensor can be calibrated, which may not be possible for large trees. Additionally, this approach is only applicable once data has been collected and a V m_critical value cannot be assigned prior to the commencement of a data campaign. Therefore, this statistical approach is called a posteriori analysis-an analysis derived from observational data.
An a priori analysis, or an analysis derived from theoretical deduction, can also be conducted to determine V m_critical . The theoretical minimum and maximum V h for Tmax and SRFM can be determined via Equations (1), (3), and (4). Tmax is typically limited by a 1 s cycle speed of contemporary data-loggers, and the SRFM is limited by the maximum value of ∆T d /∆T u of 20 [4,7]. Furthermore, the minimum and maximum V h for Tmax and SRFM, respectively, vary with k. The observable range in plants for k is approximately 0.0015 to 0.004 cm 2 /s [4,20]. Table 1 displays the minimum and maximum V h for Tmax and SRFM, respectively, for varying k with a probe spacing of 0.6 cm and heat pulse of 3 s duration. The values in Table 1 will vary with different probe spacings and heat pulse durations.
An a priori analysis may choose with SRFM max or Tmax min as the V m_critical . However, the values listed in Table 1 are theoretical and, in practice, observable SRFM max or Tmax min may be unknown. Therefore, it is proposed that the mid-point between SRFM max and Tmax min is adopted for V m_critical via the following equation: As an example, calculation based on data presented in Table 1, for a x of 0.6 cm, a t 0 of 3 s, and k of 0.0025 cm 2 /s, yields a V m_critical value of 26.3 cm/h. An a priori analysis may choose with SRFMmax or Tmaxmin as the Vm_critical. However, the values listed in Table 1 are theoretical and, in practice, observable SRFMmax or Tmaxmin may be unknown. Therefore, it is proposed that the mid-point between SRFMmax and Tmaxmin is adopted for Vm_critical via the following equation: As an example, calculation based on data presented in Table 1, for a x of 0.6 cm, a t0 of 3 s, and k of 0.0025 cm 2 /s, yields a Vm_critical value of 26.3 cm/h. Figure 1. The Dual Method Approach (DMA) algorithm to decide which heat velocity measurement (Vm) should be taken from HRM (heat ratio method)/SRFM ("slow rates of flow" method) of Tmax methods. The Vm_critical is a theoretically (a priori) or statistically (posteriori) determined value.

Converting Heat Velocity to Sap Flux Density
Once Vh has been determined via the SRFM or Tmax method for the DMA, it is necessary to correct for wounding, then convert corrected heat velocity (Vc) to sap flux density (Js).

Converting Heat Velocity to Sap Flux Density
Once V h has been determined via the SRFM or Tmax method for the DMA, it is necessary to correct for wounding, then convert corrected heat velocity (V c ) to sap flux density (J s ).

The wound effect from insertion of probes into sapwood is corrected via V c (cm/h):
where a, b, and c are coefficients that are calculated via finite-difference numerical modelling and vary with wound width, probe size, and spacing [7,26].
The V c is then converted to J s (cm 3 /cm 2 /h), which is an equation that incorporates wood and sap parameters of the measured plant [27]:

Study Species and Plant Material
Plant material was sourced from Tecoma capensis (Thunb.) Lindl., a fast-growing southern Africa species which is an exotic weed in Australia [28]. Stems were harvested from a shrub growth form of Tecoma, 2.3 m height, in the Kingston Heath Reserve, an urban bushland in southeast Melbourne, Victoria, Australia. Ten stems were used in this study, between 1.1 and 1.3 cm in diameter and 20 cm in length, cut from eight branches from a single Tecoma shrub. Samples were harvested across two days, with four branches harvested on each day. Five stems were measured on the first day, and an additional five stems were measured on the second day. Following harvesting from the shrub, the cut end of the branch was wrapped in a moistened tissue and bagged. Samples were then taken to a nearby Implexx Sense laboratory (approx. 5 min drive) for immediate testing. In the laboratory, stems were recut underwater just prior to commencement of measurements to minimize xylem embolism.

Heat Pulse Velocity Sensor & Data Logger
A commercially available heat pulse velocity sensor (HPV-06, Implexx Sense, Melbourne, Australia) was used for measurements. The HPV-06 sensor consisted of three probes, 30 mm in length, 1.3 mm in diameter, and 6 mm axial distance apart ( Figure 2). The downstream and upstream probes were temperature probes, and the central probe was the line heater. The heater probe produced 4 W to 4.9 W of power, depending on battery voltages typical of a charge cycle. Temperature measurements were made via two 10 kΩ negative temperature co-efficient (NTC) thermistors in each probe spaced at 10 mm and 20 mm consecutively from the tip of the probe. Each thermistor was wired as a Half bridge in series with a precision 24.9 kΩ resistor mounted on the circuit board in the sensor body. The accuracy of the thermistors was ±0.2 • C, and the resolution was 0.001 • C. The HPV-06 was connected to the CR1000X data-logger (Campbell Scientific, Logan, UT, USA) and individual heat-pulse events were at 10 min intervals. thermistors in each probe spaced at 10 mm and 20 mm consecutively from the tip of the probe. Each thermistor was wired as a Half bridge in series with a precision 24.9 kΩ resistor mounted on the circuit board in the sensor body. The accuracy of the thermistors was ±0.2 °C, and the resolution was 0.001 °C. The HPV-06 was connected to the CR1000X data-logger (Campbell Scientific, Logan, UT, USA) and individual heat-pulse events were at 10 min intervals. indicate the outer and inner measurement zone, respectively, with an effective measurement radial radius of 5 mm. T d and T u represent the downstream and upstream temperature probes, respectively. Each black enclosed circle represents a thermistor measurement point. The central probe is the heater. The distance between T d or T u and the Heater was 6 mm.

Sensor Installation
The HPV-06 sensor was installed using a drill guide into the Tecoma stem. Sections of stem internode between 1.1 and 1.3 cm diameter were chosen as sites of installation. Prior to insertion, the HPV-06 probes were lubricated with electrical grease (Inox MX6 Grease) to improve probe thermal contact with sapwood. The outer thermistor of the HPV-06 sensor was carefully positioned in the middle of the stem, as the xylem radius was~0.5 cm and the zone of detection of the thermistor also had a~0.5 cm radius. Therefore, positioning the thermistor in the middle of the stem captured the entire area of sapwood. Following insertion of the HPV-06 sensor in the Tecoma stem, a segment of the stem 20 cm in length, without leaves or branches, was cut underwater. This segment, with the sensor installed, remained underwater until the experimental apparatus was ready, and was immediately transferred at the commencement of measurements.

Wood Properties and Thermal Diffusivity
An additional five stem segments,~3.0 cm in length and~1.0 cm in diameter, were cut underwater from the sample branches for measurement of wood properties, including sapwood fresh weight, sapwood dry weight (oven-dried at 60 • C for 72 h), and sapwood fresh volume, and the subsequent calculation of m c , ρ d , k Hogg , and k Vand , using Equations (5)-(9) (see above). The average of the five stem segments was subsequently used in converting heat velocity to sap flow and data analyses.

Heat Velocity and Sap Flux Density Measurements
During the experiment, t m and ∆T d /∆T u were logged and later used for the calculation of V h . Six methods to calculate V h were used in this study: V h_Mar with k Hogg (Equations (4) and (7)), V h_Mar with k Vand (Equations (4) and (8)), V h_Coh with k Hogg (Equations (2) and (7)), V h_Coh with k Vand (Equations (2) and (8)), V h_Klu with k Hogg (Equations (3) and (7)), and V h_Klu with k Vand (Equations (3) and (8)). For each method, V h was then converted to V c , where the wound diameter was determined visually with a digital caliper, following Equation (11) and Burgess et al. [7]. Then, for each method, V c was converted to J s following Equation (12).

Gravimetric Testing of Sap Flow Measurements
The accuracy of sap flux density, derived from the DMA, HRM/SRFM, and Tmax methods, was tested via a water pressure system maintaining a constant head, based on the experiments of Steppe et al. [29] and López-Bernal et al. [30]. Teflon tubing with an inner diameter of 1.3 cm was attached to the upstream end of the stem segment, and was made water-and gas-tight via pipe clamps and polytetrafluoroethylene (PTFE) thread seal tape (Kinetic Supply Pty Ltd., Broadmeadow, Australia) and was thoroughly tested for leaks. Various length tubing was connected to a water reservoir that was changed in height to create variable pressure gradients. The heights were (cm) : 15, 30, 45, 60, 75,  90, 105, 120, 135, 150, 170, 200, 235, 255, and 315. At each height, flow was left to settle for at least 20 min prior to duplicate measurements at 10 min intervals. Sap passing through the downstream cut segment was collected and measured every 10 min using a three-point precision balance (WTC 200, Radwag, Radom, Poland). Concurrent sap flow data, estimated via the sap flow method and measured from the electronic balance (gravimetric sap flow), were then compared.

Statistical Analysis
Data measured on the ten stems were collated into a linear regression analysis that compared J s estimated from the V h method against J s measured gravimetrically. Linear regression slopes were fitted through the intercept. A slope of 1, from the linear regression curve, indicated zero error between estimated and measured J s , and the R 2 was a measure of precision. For the posteriori determination of V m_critical , a piecewise linear regression model, with gravimetric J s on the X-axis and estimated J s from the sensor on the Y-axis, was used to determine the breakpoint.

Measurement Range Limitation of SRFM and Tmax
The results confirmed the measurement range limitation of the HRM/SRFM and Tmax methods ( Data measured on the ten stems were collated into a linear regression analysis that compared Js estimated from the Vh method against Js measured gravimetrically. Linear regression slopes were fitted through the intercept. A slope of 1, from the linear regression curve, indicated zero error between estimated and measured Js, and the R 2 was a measure of precision. For the posteriori determination of Vm_critical, a piecewise linear regression model, with gravimetric Js on the X-axis and estimated Js from the sensor on the Y-axis, was used to determine the breakpoint.

Measurement Range Limitation of SRFM and Tmax
The results confirmed the measurement range limitation of the HRM/SRFM and Tmax methods ( Sap flux density values ranged from approximately 0 to 170 cm 3 /cm 2 /h ( Figure 3). However, the piecewise linear regression analysis determined that there was a breakpoint in the HRM/SRFM data at 32.540 cm 3 /cm 2 /h (Figure 4), which was considered the maximum reliable sap flux density measurement derived from the HRM/SRFM.    Table 2.

Accuracy of HRM/SRFM, Tmax and DMA
Across the full range of measured flow, the HRM/SRFM had the lowest levels of accuracy and precision, followed by the Tmax and DMA methods ( Table 2 and Supplementary Materials). The HRM/SRFM had an error between 49.1% and 63.4%, and underestimated gravimetric Js measurements. The Tmax had an error between 0.5% and 35.5% but overestimated gravimetric Js measurements. The DMA had an error between 1.2% and 29.0% and either closely estimated, or overestimated, Js measurements. The linear regression statistical analyses showed an R 2 of 0.458 for HRM/SRFM, 0.807 for the Tmax and 0.888 to 0.918 for the DMA ( Table 2). The results from the root mean square error (RMSE) of the linear regression analyses was consistent with this result, showing the DMA posteriori analysis as having the lowest RMSE, indicating this was the most accurate model.
An additional problem was observed with the Tmax measurements at low range. There were data points outside of the minimum measurement range (Table 1) Table 2. Sap flux density values ranged from approximately 0 to 170 cm 3 /cm 2 /h ( Figure 3). However, the piecewise linear regression analysis determined that there was a breakpoint in the HRM/SRFM data at 32.540 cm 3 /cm 2 /h (Figure 4), which was considered the maximum reliable sap flux density measurement derived from the HRM/SRFM.

Accuracy of HRM/SRFM, Tmax and DMA
Across the full range of measured flow, the HRM/SRFM had the lowest levels of accuracy and precision, followed by the Tmax and DMA methods ( Table 2 and Supplementary Materials). The HRM/SRFM had an error between 49.1% and 63.4%, and underestimated gravimetric J s measurements. The Tmax had an error between 0.5% and 35.5% but overestimated gravimetric J s measurements. The DMA had an error between 1.2% and 29.0% and either closely estimated, or overestimated, J s measurements. The linear regression statistical analyses showed an R 2 of 0.458 for HRM/SRFM, 0.807 for the Tmax and 0.888 to 0.918 for the DMA ( Table 2). The results from the root mean square error (RMSE) of the linear regression analyses was consistent with this result, showing the DMA posteriori analysis as having the lowest RMSE, indicating this was the most accurate model. An additional problem was observed with the Tmax measurements at low range. There were data points outside of the minimum measurement range (Table 1) causing errors in the calculation of Tmax at slow flows. Therefore, there was a lower sample size (n) for Tmax measurements shown in Table 2 and the Supplementary Materials.

Determination of V m_critical and the DMA
The V m_critical level, determined via the a priori and the posteriori methods ( Figure 4) yielded similar results ( Table 2). The slope of the linear regression curve of J s measured via the DMA against the gravimetric technique showed a value of 0.984 and 0.988 for the a priori and posteriori methods, respectively. The precision of the data was also similar, with an R 2 of 0.888 and 0.918 for the a priori and posteriori methods, respectively. The RMSE analysis was also similar, but indicated that the posteriori approach was more accurate than the a priori approach ( Table 2).

Accuracy of Cohen's versus Kluitenberg's Tmax Calculation
The different methods for calculating Tmax (V h_Coh and V h_Klu ; Equations (2) and (3)) showed varying levels of accuracy and precision. The V h_Coh equation showed better accuracy but slightly lower precision, with a linear regression curve slope of 1.005, R 2 of 0.807, and RMSE of 12.816 (Supplementary Materials). The V h_Klu equation, in contrast, resulted in a linear regression curve slope of 1.178, R 2 of 0.885, and RMSE of 14.300 (Supplementary Materials). Tmax, calculated via V h_Coh , had a smaller RMSE, indicating this model was more accurate than the V h_Klu model.

Accuracy of Hogg's versus Vandegehuchte's Thermal Diffusivity (k) Calculations
The k value differed according to calculations via Hogg's (k Hogg ) or Vandegehuchte's (k Vand ) method. The mean k Hogg value was 0.001648 cm 2 /s (±SD 9 × 10 −5 ) and the mean k Vand value was 0.002293 cm 2 /s (±SD 5 × 10 −5 ). The difference in k resulted in different levels of accuracy when sap flux density was calculated using k Hogg or k Vand (Supplementary Materials). Generally, J s calculated with k Hogg yielded lower accuracy and precision than J s calculated with k Vand .

Discussion
The DMA successfully resolved the measurement range limitation of the Tmax and HRM/SRFM heat pulse velocity methods. In this study, the DMA measured near zero to extremely fast rates of sap flow. Although reverse flows were not explicitly tested in this study, the probe design and theoretical calculations suggest the DMA can also resolve upstream movement of xylem sap.
The measurement range limitation of nearly all sap flow methods has restricted their utility for research and applied applications. For example, Forster [31] demonstrated that, on average, 12% of sap flow occurs during the night, or during periods of slow sap flow, which cannot be quantified by many sap flow methods. Conversely, economically important crops, such as grapevines, are known to have fast daytime flows, which cannot be quantified by the HRM [8,10]. The DMA is a practical method to implement that can successfully quantify sap flow under these previously limiting applications.
This study advances previous attempts to combine heat pulse velocity methods to overcome the measurement range limitation. For example, Cohen et al. [32,33] combined the Tmax method and compensation heat pulse method (CHPM) for measuring slow to very high rates of flow, and Pearsall et al. [10] combined the CHPM and HRM to measure reverse, zero, slow, and high rates of flow in grapevines. The CHPM cannot measure zero or reverse flows, with one study reporting an observed minimum flow of 1.8 cm/h [34]. Therefore, the combination of CHPM with Tmax does not resolve the entire measurement range limitation observed in plants. The combination of CHPM and HRM/SRFM by Pearsall et al. [10] did overcome the measurement range limitation, as CHPM can measure high flows and HRM/SRFMcan measure slow, zero, and reverse flows. However, combining CHPM with HRM/SRFM requires either two installation sites or up to four probes [10], which can lead to errors caused by probe misalignment and wounding. The DMA advocated in this study reduces the number of probes to three, which minimizes errors caused by probe misalignment and wounding. Additionally, Pearsall et al. [10] relied on a visual assessment to determine when to switch from HRM/SRFM to CHPM, whereas a statistical and theoretical approach was demonstrated in this study to determine the V m_critical to switch from one heat pulse velocity method to another.

Accuracy of the DMA, HRM/SRFM and Tmax
When compared with a gravimetric measure of J s , the DMA had a slope close to 1, or an error of 1.2% and an R 2 >0.842. This compares with an average error of 16.9% and R 2 of 0.916 for the HRM/SRFM, and an error of 35.560% and R 2 of 0.859 for the Tmax reported in a meta-analysis [8]. Therefore, the DMA significantly improved the accuracy of J s estimations with similar precision to existing sap flow methods.
In this study, the HRM/SRFM poorly estimated J s , with an error between 49% and 63% ( Table 2, Supplementary Materials), compared with an average error of 16.9% from a meta-analysis of 11 studies [8]. However, Vandegehuchte and Steppe [17] and Wang et al. [35] found an error of 39% and 43%, respectively, when the HRM was compared against a gravimetric measure of sap flow. Other studies testing the accuracy of the HRM found better accuracy at slow flows (e.g., [36,37]) than fast flows (e.g., [10,35]). Therefore, within its limited measurement range, the HRM may be an accurate method to estimate sap flow, but it should not be deployed on stems with fast sap fluxes.
In contrast, this study found an error as low as 0.5% for the Tmax, similar to a 0% error reported by Green [38], but contrasted with errors between 19% and 65% reported by Intrigliolo et al. [39]. The high levels of error reported by Intrigliolo et al. [39] were measured on Vitis vinifera stems, a reportedly difficult plant to measure accurately with sap flow sensors [10], and compared against a canopy gas exchange instrument. It is possible that the canopy gas exchange instrument did not accurately reflect sap flow dynamics, as it measures different plant physiological processes. The 0% error reported by Green [38] was a study also conducted on Vitis vinifera, but against a gravimetric measure of sap flow.
Heat pulse velocity sap flow sensors have an additional limitation of only measuring a small area of the potential total sapwood area in plants. In this study, this potential source of error was removed because a small stem was measured that had the same area which could be detected with the HPV-06 sap flow sensor. In other conditions, such as when measuring large trees in a forest or orchard, a heat pulse velocity sap flow sensor will potentially sample only a portion of the total cross-sectional sapwood area. Previous studies, for example Eliades et al. [40], deployed sensors at multiple sapwood radial depths and azimuthal locations to increase the sampled sapwood area of the measured tree. The DMA can overcome potential sampling errors by also measuring at multiple radial and azimuthal locations. The DMA can better address sampling issues than other heat pulse velocity methods, as it is the only method capable of measuring the entire observable range of sap flow in plants of all sizes.

Determination of V m_critical
A statistical (posteriori) and theoretical (a priori) approach were used to determine V m_critical , with almost identical results. The piecewise linear regression used to determine a breakpoint is a robust method to determine when the HRM/SRFM reaches its maximum measurement range. It is recommended to adopt a statistical approach, where possible, to determine V m_critical to avoid bias and error in the estimation of the maximum HRM/SRFM measurement range. Currently, the maximum range of HRM/SRFM cannot be predicted, because the causes are varying and poorly understood [8]. It is hypothesized that at high flows, the ∆T u parameter in Equation (4) is constant whereas ∆T d decreases, leading to errors in calculating V h [3,17]. However, it is not certain how or when this occurs. Forster [8] noted that, under extremely high flow conditions, the heat pulse travels outside the zone of measurement of the downstream temperature probe, which is possibly related to a decrease in the q parameter, or heat input, of Equation (1). The maximum range of the HRM/SRFM is also related to the plant species and its hydraulic properties. These hypotheses have not been systematically tested; therefore, it cannot currently be predicted when the HRM/SRFM will reach its maximum limit.
In many applications, it will not be feasible to wait until the end of a study period to determine V m_critical . In such applications, the user requires reliable sap flow data in, or near, real time. The a priori method to determine V m_critical can be adopted for such applications which, in this study, yielded almost identical results to the statistical, posteriori approach. The a priori approach used a maximum ∆T d /∆T u ratio of 20, which was hypothesized by Marshall [4], and later used by Burgess et al. [7], to determine the maximum measurement limit of the HRM/SRFM. The importance of the ∆T d /∆T u ratio of 20 is unknown, as Marshall nor Burgess et al. did not provide an explanation for why the value of 20 is important. In the author's experience, ∆T d /∆T u ratios rarely exceed 5 (personal observation); therefore, a value of 20 may be an overestimation. Values less than 20 may be equally valid in the determination of V m_critical . A more robust examination of the optimal ∆T d /∆T u ratio for the calculation of V m_critical , which was outside the scope of this study, is required.

Methods to Estimate Thermal diffusivity
The k parameter in HPV calculations is critical, yet some studies adopt Marshall's suggestion of 0.0025 cm 2 /s as a default value (e.g., [41]). Many studies measure k at the start or end of a measurement campaign through a stem core and the use of Equation (7) (k Hogg ; e.g., [42,43]). Looker et al. [20] and Vandegehuchte and Steppe [19] suggested that k Hogg may underestimate k, and instead recommended the use of Equation (8) (k Vand ), as this method explicitly accounts for the bound and unbound moisture content of sapwood. The results of this study support the conclusion that k Hogg is smaller than k Vand . Furthermore, in this study J s calculated via k Vand was more accurate and precise than J s calculated with k Hogg . Therefore, the use of k Vand may also improve the overall accuracy of sap flow measurements in general. For example, Forster [8] found that sap flow estimates from HPV methods underestimated true sap flow by an average of 34.706%. Most of the studies collated as part of Forster's [8] meta-analysis used k Hogg in calculating sap flow estimates. The use of k Vand , instead, may result in improved accuracy of some of these studies, as suggested by the results of this and other studies (e.g., [44]).

Conclusions
The DMA successfully resolved the measurement range limitation of heat pulse velocity methods and reduced errors from as high as 63.4% to 1.2% without loss of precision. It is recommended that the DMA is adopted for sap flow measurements, rather than the HRM/SRFM or Tmax, as it supersedes the utility of these latter methods. The appropriate calculation of thermal diffusivity via using the k Vand method also improved the accuracy of sap flow measurements.
Supplementary Materials: The following are available online at http://www.mdpi.com/1999-4907/10/1/46/s1, Supplementary Material S1: Detailed results from the piecewise linear regression analysis of the various sap flow methods against gravimetric measurements of sap flow.

Conflicts of Interest:
The author is an owner of Implexx Sense-a company that manufactures and distributes sap flow sensors commercially.