On the Use of the Eddy Covariance Latent Heat Flux and Sap Flow Transpiration for the Validation of a Surface Energy Balance Model

Actual evapotranspiration is assessed via surface energy balance at an hourly rate. However, a robust estimation of daily evapotranspiration from hourly values is required. Outcomes of surface energy balance are frequently determined via measures of eddy covariance latent heat flux. Surface energy balance can be applied on images acquired at different times and spatial resolutions. In addition, hourly actual evapotranspiration needs to be integrated at a daily rate for operational uses. Questions arise whether the validation of surface energy balance models can benefit from complementary in situ measures of latent heat flux and sap flow transpiration. Here, validation was driven by image acquisition time, spatial resolution, and temporal integration. Thermal and optical images were collected with a proximity-sensing platform on an olive orchard at different acquisition times. Actual latent heat fluxes from canopy and sap flux at tree trunks were measured with a flux tower and heat dissipation probes. The latent heat fluxes were then further analyzed. A surface energy balance was applied over proximity sensing images re-sampled at different spatial resolutions with resulting latent heat fluxes compared to in situ ones. A time lag was observed and quantified between actual latent heat fluxes from canopy and sap flux at the tree trunk. Results also indicate that a pixel resolution comparable to the average canopy size was suitable for estimating the actual evapotranspiration via a single source surface energy balance model. Images should not be acquired at the beginning or the end of the diurnal period. Findings imply that sap flow transpiration can be used to measure surface energy balance at a daily rate or when images are found at an hourly rate near noon, and the existing time lag between the latent heat flux at the canopy and the sap flow at the trunk does not need to be taken into account.


Introduction
An accurate estimation of agricultural water is essential since it accounts for 70% of total global freshwater withdrawals and is expected to increase a further 10% by 2050 [1].Thus, implementing wise agricultural water practices is an important task for global hydrologic cycle management.For precision irrigation, both the amount and timing of water supply are important.These can be obtained using criteria of plant physiological status such as the plant water potential needed to manage the water stress of the crop.An inverse relationship exists between predawn leaf water potential and the ratio between actual and reference evapotranspiration [2].On the other hand, remote-sensing-based energy balance models allow for characterization of the spatial distribution of actual hourly evapotranspiration, ET (mm•d −1 ).This allows for the extension of point-based to spatially distributed water potential.
The best known and widely accepted methodology for validating surface energy balance models is the eddy covariance (EC), which allows for the measurement of both latent heat flux, λET (W•m −2 ), and sensible heat flux, H (W•m −2 ) [3].Some authors [4] have demonstrated that the joint use of eddy covariance and sap flow (SF) techniques is effective for evaluating the components of evapotranspiration in orchards (i.e., olive trees).
Surface energy balance estimations can be validated by comparing micro-meteorological measures (e.g., acquired by flux towers, FT) and eco-physiological sensors (e.g., sap flow).More specifically, sap flow is based on a heat dissipation technique [5] suited to measure canopy transpiration.Flux towers, on the other hand, are based on the eddy covariance technique [6], which provides an estimation of the sum of surface evaporation (coincident with soil evaporation over dry canopy) and canopy transpiration.
Surface energy balance models can estimate instantaneous flux values using images acquired at different times and spatial resolutions.Instantaneous values are assumed to represent hourly fluxes.For operational purposes, actual evapotranspiration needs to be estimated over longer time periods.Several methods have been proposed to time-scale instantaneous λET into evapotranspiration daily values, ET D (mm•d −1 ) (from now on, the subscript "D" indicates daily variables) [7].Several authors (e.g., [8,9]) use a method based on two main hypotheses: (i) the evaporative fraction at the acquisition time, Λ(-), is assumed equal to its average diurnal value; in other words, the fraction of available energy used for the evapotranspiration process is almost constant during diurnal hours (self-preservation hypothesis); (ii) the daily soil heat flux at ground level can be neglected when compared to the daily net radiation, R n,D , since, during a day cycle, daytime heat flux almost balances the nighttime one.
These two hypotheses are uncertain for the estimation of daily evapotranspiration, so the application of this time-scaling approach might be limited under certain meteorological conditions.Changes in temperature, humidity, and wind speed may play a prominent role when partitioning the energy budget, and intermittent cloudiness may cause a significant change in Λ [9,10].However, Crago [10] showed that the self-preservation hypothesis is surprisingly robust on clear-sky days.
An alternative and simple upscaling approach hypothesizes a proportionality between evapotranspiration and incoming shortwave radiation, R s (W•m −2 ), during the whole day [11], which can be applied in partially clear-sky conditions as well.
The main objective of this research was to assess whether or not evaporation and transpiration fluxes can be determined at hourly and daily scales using flux towers and sap flow measures.The research also aimed to point out the advantages and limitations of each technique.The analysis was carried out by applying two different time-scale methods to data obtained at varying spatial resolutions and different times.
To this goal, latent heat flux and sap flux were monitored at an olive orchard by deploying a flux tower and heat dissipation probes for fruit development and berry growth stages.At the end of the investigation period, an experimental platform, equipped with thermal and optical radiometric cameras, was used to estimate the spatial distribution of latent heat flux via surface energy balance.

Methods
Latent heat flux estimated via surface energy balance models can be validated via eddy covariance and sap flow measurements.Other measurements are based on scintillometry, lisimetry, and surface renewal (e.g., [12][13][14]).

The Surface Energy Balance Model
Latent heat flux can be determined as a residual term of the surface energy balance equation: where R n , G 0 , and H are, respectively, the net radiation, the ground heat, and sensible heat fluxes at the acquisition time (W•m −2 ); the difference (R n − G 0 ) is the available energy [15].SEBAL is a single-source surface energy balance model schematizing soil-vegetation as a sole resistance layer ([8,16,17]).ET D is extrapolated from hourly values using Λ as an integration parameter.
The diurnal average evaporative fraction, Λ d , is obtained by increasing Λ by 10% [15] to account for the Λ temporal behavior during cloudless days.Indeed, hourly Λ usually exhibits a slightly concave-up shape [18] characterized by a minimum near solar noon and two maxima near sunrise and sunset (with high variability before and after sunrise and sunset, respectively).
The use of Λ d implies that nighttime λET is small and thus negligible compared to daily λET [9].Accordingly, Malek [19] found that nighttime λET (from sunset to sunrise) was 1.7% of λET D during a complete growing cycle of alfalfa.However, this ratio increased up to 14% when the nighttime wind speed was high.
The instantaneous daily integration process assumes that during a daily cycle the sum of incoming and outgoing ground heat fluxes is negligible compared to the daily net radiation [8].It also assumes that Λ can be supposed as a constant if limited to only diurnal hours of cloudless days (evaporative fraction self-preservation).
Although this latter assumption has been documented (e.g., [17,20]), it could lead to the underestimation of ET D in arid climates where afternoon advection (or increased afternoon wind speed) may increase λET in proportion to available energy [21].Diurnal course of Λ responds to several factors including the relative amplitude of turbulent heat fluxes.The validation of λET by using the difference between the G 0 and R n (the available energy) explains the Λ daily shape.In addition, the daytime behavior of Λ depends on environmental factors including a sharp rise in soil water content and leaf area index, LAI (m 2 •m −2 ) [22].
Due to these factors, the inherent concave-up shape of Λ causes underestimation of daily evapotranspiration if the instantaneous evaporative fraction value is assumed to represent the diurnal average value.Experimental work (e.g., [23,24]) has also demonstrated that the self-preservation hypothesis holds true for environmental conditions where soil water content does not significantly change and advection does not occur.
According to Chàvez et al. [7], the upscaling of ET to ET D (mm•d −1 ), under the hypothesis of self-preservation of Λ, is obtained as follows: where 86,400 is a time unit conversion factor (s•d −1 ) for ET D , λ is the latent heat of vaporization (MJ•Kg −1 ), and ρ w is the water density (kg•m −3 ).
The second approach estimates the ET D using R s as a reference integration variable [25] under the hypothesis of proportionality between ET and R s at a daily scale.Even though this approach requires the measurement of solar irradiation as additional input, it is sometimes preferred to Equation (3), which can have inaccuracies associated with the variability of cloud presence (Equation ( 4)):

The Flux Time-Scaling via Micro-Meteorological Measurements
Integration based on Λ self-preservation hypothesis was tested using EC tower data.R n is obtained by balancing short-wave and long-wave radiations measured through a four components net radiometer.G 0 is measured through a self-calibrating soil heat plate.λET is calculated based on the covariance between the vertical velocity fluctuation of wind and the water vapor density.
H and λET were derived from high frequency sonic air temperature and water vapor concentration.Both H and λET were determined in phase with three-dimensional wind speed measurements, according to the eddy covariance [6].H and λET were adjusted by introducing several corrections to raw data, as described in Cammalleri et al. [9].
However, failure in the closing of the surface energy balance often occurs due to several factors, including the averaging and coordination of rotation periods, a difference in the fluxes source areas, and an influence on the low-frequency part of the turbulence spectra [26].Thus, surface energy balance components were adjusted by including the energy closure with the goal of preserving the Bowen ratio, β R = H•λET −1 [27].Measures that have a surface energy balance closure error within ±15% were considered reliable [28].According to Schuepp et al. [29], the cumulative normalized contribution to the flux measurements allowed for estimation of the upwind distance, which gave a significant fraction of fluxes (e.g., 70% of the fluxes, FP 0.70 (m)), as well as the distance for the maximum flux contribution, X max (m).

Tree Transpiration Measurement: The Heat Dissipation Technique
Tree transpiration can be estimated via the heat dissipation technique [30] by measuring the sap velocity of the woody stems (5): The plant transpiration T p (mm) is obtained by measuring the temperature difference, ∆T (K), between a heated upper probe and an unheated lower probe radially inserted into the trunk to intercept the sapwood.The faster the sap fluxes, q (m 3 •s −1 ), the more the heat dissipates and the lower ∆T.When sap flux rate is negligible (during night time) ∆T is the maximal, ∆T max .This approach allows for estimation of the q of a single plant by multiplying the sap flow velocity for the cross-sectional area of the sapwood, S (m 2 ); q is then upscaled through a plant pertinence area, A p (m 2 ).The coefficients α and β are assumed to be equal to 0.714 ms −1 and 1.213, respectively [5].The daily volume of water consumed by a single plant is determined via daily integration every half hour, called the time step, t, assuming the effect of the tree capacity at the daily time interval is negligible [31].
A further scaling step was required to evaluate a representative value of the stand transpiration over a plot, T (mm), from that of individual plants.
In particular, T p was upscaled by considering, as a proximal variable, the ratio between the average plot LAI and the LAI of the single plants, LAI p (m 2 •m −2 ), in which sap flow was monitored, and where overbars denote average values (Equation ( 6)).

Eddy Covariance and Sap-Flow Diurnal Behavior and Time Lag between EC and SF Measurements
Comparing sap flow and eddy covariance measurements is necessary for analyzing the daily behavior of transpiration T and ET.Although EC data are sampled at high frequency (20 Hz), measurements are processed on a half-hour scale.Thus, the SF time series were averaged accordingly.Data normalization was applied to take into account upscaling inaccuracies of SF measurements and scale mismatching among the tree canopies sap flow and the eddy covariance footprint.By dividing each measurement by the total (e.g., the sum of measurements of the whole day), this normalization standardizes and scales the measurements.This allows for the comparison between the diurnal behavior of ET and T during the studied period.The average daily behavior (1-day cycle) of ET vs. T was derived by averaging 30 min measurements.This was completed for a 5-week period from 29 June, 12:00 p.m., to 2 August, 11:00 p.m. UTC (all times are reported in UTC).Finally, the determination coefficient, r 2 (-), between the daytime measures of ET and T, was maximized by shifting forward in time the latter time series, allowing for the quantification of the time lag between T at the trunk and ET at the canopy.

The Surface Energy Balance Model Validation
The micro-meteorological measurements of λET were used to validate remote sensing surface energy balance estimations over the flux tower footprint.To figure out whether soil evaporation was negligible, sap flow measures (i.e., T) were compared with FT measures (i.e., ET).
The absolute difference, AD (-), was calculated as evaluation metric (7): where AD (7a) indicates λET absolute differences compared to eddy covariance λET closed with the Bowen ratio approach; AD Λ and AD Rs indicate absolute daily ET D difference based on Equations ( 3) and (4), respectively.For the absolute difference calculations, eddy covariance λET and ET D were adjusted, which forced the energy closure by preserving the value of the Bowen ratio, β R .

Spatial Resolution Analysis
Both remote-sensing-based λET and ET D retrievals vary depending on the spatial resolution, which impacts the validation of surface energy balance models with micro-meteorological or sap flow data.To quantify these effects, original multi-spectral and thermal images were downscaled from the raw resolutions to coarser resolutions characterizing the canopy sizes of the orchard trees.Thus, the average of the footprint area was compared with FT fluxes and with upscaled T obtained using SF measures.
To verify whether λET values can represent hourly values, a coefficient of variation, CV (unitless), was evaluated over half-hourly EC λET measures.The CV is computed over time by applying a 1 h moving window (3 values); CV is then assigned to the median acquisition time.

The Experimental Farm
The "Tenuta Rocchetta-Azienda Agricola Angela Consiglio" experimental farm (Figure 1) is located along the southwest coast of Sicily, Italy.The farm lies several kilometers southeast of the Castelvetrano town by the lower Belice river valley.The district is mostly agricultural and is characterized by a typical Mediterranean climate with moderate rainfall in the autumn and winter seasons and where the summer remains hot and dry.The farm has a 2 ha olive orchard (Nocellara del Belice cultivar).Olive is the main crop in the district.
Trees are planted on a regular grid of 8 × 5 m (250 plants•ha −1 and a plant pertinence area equal to 40 m 2 ).The olive canopies are characterized by an average height of ≈3.3 m, an average fraction of vegetation cover of ≈0.35 m, and a canopy size ranging between ≈2.1 and ≈5.6 m (on average ≈4.5 m).A flat landscape and homogeneous soil characterize the study area.The soil is called silty clay loam (according to the United States Department of Agriculture-USDA classification), with 60% of the soil in the form of sand and 24% in the form of clay.Due to significant evaporation in the summer season and the near absence of rainfall during the crop season, agricultural fields are occasionally irrigated (aid irrigation).Water is supplied by a drip irrigation system (four emitters per plant, each flowing ~8 lit•h −1 of water).An eddy covariance system and sap flow sensors are installed in the olive orchard (Figure 1) to monitor the evapotranspiration and transpiration processes, respectively.
Remote Sens. 2018, 10, x FOR PEER REVIEW 6 of 18 of water).An eddy covariance system and sap flow sensors are installed in the olive orchard (Figure 1) to monitor the evapotranspiration and transpiration processes, respectively.Eco-physiological and micrometeorological data were pre-processed in order to dispose of a data set distributed during the mid-season from June to the first week of September, which corresponds from the initial fruit development stage to the berry growth stage.

In Situ Measurements: Sensors and Data
Proximity sensing platform-The acquisition system encompassed an MCA-II multispectral camera (TETRACAM Inc.: C ， hatsworth CA, USA) and an A320G thermal camera (FLIR system Inc.: Wilsonville, Oregon, USA).The multispectral camera is an 8 bit 6 bands CCD sensor with customized optical filters set up at 450, 650, 720, and 800 nm.The field of view was 43° × 35°.The A320G micro-bolometric thermal camera acquired thermal infrared images in the spectral broadband range of 7.5-13 μm when the field of view was 45° × 34°.
A remote laptop, which was connected to an on-board mini personal computer, controlled both cameras.The acquisition instruments were installed on an 8 m 3 helium balloon (Figure 2) flying at 120 m above ground level (a.g.l.).The balloon was secured to the ground through three Kevlar ® ropes.Afterward, the balloon was slowly let go.An ASD Field Spec Hand-Held spectroradiometer (Analytical Spectral Device, Inc., Boulder, CO, USA) enabled the study of ground reflectance between 325 and 1075 nm, which were used for in-reflectance calibration at the ground level.Eco-physiological and micrometeorological data were pre-processed in order to dispose of a data set distributed during the mid-season from June to the first week of September, which corresponds from the initial fruit development stage to the berry growth stage.

In Situ Measurements: Sensors and Data
Proximity sensing platform-The acquisition system encompassed an MCA-II multispectral camera (TETRACAM Inc.: Chatsworth, CA, USA) and an A320G thermal camera (FLIR system Inc.: Wilsonville, Oregon, USA).The multispectral camera is an 8 bit 6 bands CCD sensor with customized optical filters set up at 450, 650, 720, and 800 nm.The field of view was 43 • × 35 • .The A320G micro-bolometric thermal camera acquired thermal infrared images in the spectral broadband range of 7.5-13 µm when the field of view was 45 • × 34 • .
A remote laptop, which was connected to an on-board mini personal computer, controlled both cameras.The acquisition instruments were installed on an 8 m 3 helium balloon (Figure 2) flying at 120 m above ground level (a.g.l.).The balloon was secured to the ground through three Kevlar ® ropes.Afterward, the balloon was slowly let go.An ASD Field Spec Hand-Held spectroradiometer (Analytical Spectral Device, Inc., Boulder, CO, USA) enabled the study of ground reflectance between 325 and 1075 nm, which were used for in-reflectance calibration at the ground level.
Images were acquired on 1 August 2009 at spatial resolutions of 0.12 m (multispectral) and 0.50 m (thermal infrared).Reflectance was re-sampled using pixel combinations at 2, 5, and 10 m to analyze the "suitable" structural resolutions for evapotranspiration retrievals at plant and plot scales.Radiant exitance emitted by the soil-vegetation system (assumed as a grey body) and recorded by the thermal sensor was re-sampled at lower spatial resolutions and converted into kinetic temperature.Images were acquired on 1 August 2009 at spatial resolutions of 0.12 m (multispectral) and 0.50 m (thermal infrared).Reflectance was re-sampled using pixel combinations at 2, 5, and 10 m to analyze the "suitable" structural resolutions for evapotranspiration retrievals at plant and plot scales.Radiant exitance emitted by the soil-vegetation system (assumed as a grey body) and recorded by the thermal sensor was re-sampled at lower spatial resolutions and converted into kinetic temperature.
Flux tower-The flux tower installed in the olive orchard (37.64294°Lat., 12.84710° Lon.) allowed for the monitoring of H and λET via the eddy covariance technique.The flux tower was equipped with several probes.A three-dimensional sonic anemometer (CSAT3-3D, Campbell Scientific Inc. Logan, UT, USA) was used to measure high-frequency winds at 5.6 m a.g.l.An infrared open-path gas analyzer (IRGA LI7500, Li-COR Biosciences Inc. Lincoln, NE, USA) measured high-frequency densities of water vapor and carbon dioxide at the same height of the sonic anemometer.
A four-component net radiometer (NR-Lite-L, Kipp and Zonen-Delft, Netherlands) positioned at an elevation of 8.5 m a.g.l. was used to measure direct shortwave and longwave radiations along with reflected shortwave and longwave radiations.Two soil heat flux plates (HFP01SC self-calibrating heat flux sensor, Hukseflux Thermal sensors-Delft, Netherlands) installed at a depth of 0.1 m below ground level (b.g.l.) were used to measure the ground heat flux.Two pyranometers (IRTS-P, IR Precision Infrared Thermocouple Sensor, Apogee Instruments, Inc.-Logan, UT, USA) installed at 4.80 m a.g.l. were used to measure both soil and vegetation radiometric temperature.
Sap flow-Sap flow measurements provided direct estimates of transpiration every 15 min on three olive trees.In particular, the heat dissipation method measured half an hour of sap velocity.Three olive trees (namely P1, P2, and P3) located close to the flux tower (Figure 1) were monitored.Measurements were carried out using two standard thermal dissipation probes.The main probe used was the Sapflow Sensor (SFS2 Type M-M; UP GmbH), which was installed into the trees' trunks at a height of ≈0.4 m a.g.l.According to prior footprint analysis, trees were selected at an average distance from the FT approximately equal to the daily maximum Xmax.Trees were selected according to their trunk diameter [32] to be representative of the experimental plot.The trunk diameters of P1, P2, and P3 were 0.18, 0.23, and 0.30 m, respectively.
The probes were installed on the north side of the olive trees trunk and then insulated (Figure 3b) to avoid direct sun exposure.Sap flow-Sap flow measurements provided direct estimates of transpiration every 15 min on three olive trees.In particular, the heat dissipation method measured half an hour of sap velocity.Three olive trees (namely P 1 , P 2 , and P 3 ) located close to the flux tower (Figure 1) were monitored.Measurements were carried out using two standard thermal dissipation probes.The main probe used was the Sapflow Sensor (SFS2 Type M-M; UP GmbH), which was installed into the trees' trunks at a height of ≈0.4 m a.g.l.According to prior footprint analysis, trees were selected at an average distance from the FT approximately equal to the daily maximum X max .Trees were selected according to their trunk diameter [32] to be representative of the experimental plot.The trunk diameters of P 1 , P 2 , and P 3 were 0.18, 0.23, and 0.30 m, respectively.
The probes were installed on the north side of the olive trees trunk and then insulated (Figure 3b) to avoid direct sun exposure.
To convert sap fluxes into tree transpiration, a relationship between the cross-sectional area of conducting sapwood and trunk diameter was found by coloring the sapwood area.Indeed, at the end of the experiment, the sapwood area was determined using a color-based method on a total of six wood cores extracted with a Pressler gimlet on the monitored trees.To enhance the difference between sapwood and heartwood, the conductive section of the wood adding methyl-orange to the wood core was identified.Each image of colored wood core was then analyzed with software Image-Pro Plus 6.0 (Media Cybernetics Inc.: Rockville, Maryland, USA) to define the sapwood depth.To convert sap fluxes into tree transpiration, a relationship between the cross-sectional area of conducting sapwood and trunk diameter was found by coloring the sapwood area.Indeed, at the end of the experiment, the sapwood area was determined using a color-based method on a total of six wood cores extracted with a Pressler gimlet on the monitored trees.To enhance the difference between sapwood and heartwood, the conductive section of the wood adding methyl-orange to the wood core was identified.Each image of colored wood core was then analyzed with software Image-Pro Plus

Fluxes Temporal Behavior
The 30 min SF values were averaged over a week and normalized over the daily total.These values show a circadian pattern characterized by a saturation-type response.T and ET diurnal changes (Figure 4, black and empty dots, respectively) are nearly similar during the 24 h cycle with a considerable time shift.At noon, SF fluxes were nearly steady, while FT ones followed the daily trend of atmospheric evapotranspiration demand.
In particular, the daily behavior of T is smoother than that observed for ET, which rapidly varies due to meteorological forces.The observed time lag confirmed the outcomes of other studies [33][34][35].If there is a large capacitive exchange between the transpiration stream and stem water storage, there is a significant time lag between SF measured at the base of the stem and T (at the canopy leafs).Moreover, the differences between sap flow transpiration and flux tower evapotranspiration highlight an imbalance between canopy transpiration and tree water uptake.
Figure 4 shows that the diurnal cycles may be represented by two phases.The first phase, before 2:30 p.m., where ET is higher than T, can be attributed to limited soil water content or tree capacitance when the tree is using water stored in the trunk.The second phase is after 2:30 p.m., where ET decreases and T preserves a steady phase and high values.These differences in water use dynamics could be interpreted as an effect of tree capacitance.

Fluxes Temporal Behavior
The 30 min SF values were averaged over a week and normalized over the daily total.These values show a circadian pattern characterized by a saturation-type response.T and ET diurnal changes (Figure 4, black and empty dots, respectively) are nearly similar during the 24 h cycle with a considerable time shift.At noon, SF fluxes were nearly steady, while FT ones followed the daily trend of atmospheric evapotranspiration demand.The amount of water stored in the stem introduces uncertainty when estimating diurnal whole-tree canopy transpiration and SF transpiration measured at the base of the stem.Motisi et al. [31] showed that this imbalance between canopy transpiration and tree water uptake is recovered during the night and is likely attributable to tree capacitance.Other studies have shown that, when the stem capacitance is large, an abundance of nighttime water recharges into the stem because whole-tree transpiration is higher than the stem q around mid-day [28,34,[36][37][38].
Since T exhibits almost a flat behavior close to its maximum value (Figure 4), surface energy balance outcomes can be directly measured using SF values if remote sensing images are acquired close to midday (without taking into account any time lag).It would be different if imaging and In particular, the daily behavior of T is smoother than that observed for ET, which rapidly varies due to meteorological forces.The observed time lag confirmed the outcomes of other studies [33][34][35].If there is a large capacitive exchange between the transpiration stream and stem water storage, there is a significant time lag between SF measured at the base of the stem and T (at the canopy leafs).Moreover, the differences between sap flow transpiration and flux tower evapotranspiration highlight an imbalance between canopy transpiration and tree water uptake.
Figure 4 shows that the diurnal cycles may be represented by two phases.The first phase, before 2:30 p.m., where ET is higher than T, can be attributed to limited soil water content or tree capacitance when the tree is using water stored in the trunk.The second phase is after 2:30 p.m., where ET decreases and T preserves a steady phase and high values.These differences in water use dynamics could be interpreted as an effect of tree capacitance.
The amount of water stored in the stem introduces uncertainty when estimating diurnal whole-tree canopy transpiration and SF transpiration measured at the base of the stem.Motisi et al. [31] showed that this imbalance between canopy transpiration and tree water uptake is recovered during the night and is likely attributable to tree capacitance.Other studies have shown that, when the stem capacitance is large, an abundance of nighttime water recharges into the stem because whole-tree transpiration is higher than the stem q around mid-day [28,34,[36][37][38].
Since T exhibits almost a flat behavior close to its maximum value (Figure 4), surface energy balance outcomes can be directly measured using SF values if remote sensing images are acquired close to midday (without taking into account any time lag).It would be different if imaging and subsequent measurement were performed at the beginning or the end of the diurnal period (e.g., the 8:25 a.m.image acquisition of this experiment).Finally, uncertainties about using T measurements to validate the surface energy balance estimates from remote sensing arise from the individual shape of the two curves (T p and ET).
The daily cumulative upscaled sap flow T was compared to the flux tower λET D for the day of remote sensing acquisition.The diurnal changes of T (Figure 5, black dots) and ET derived by FT (Figure 5, empty dots) show, again, the time lag between those different measurements.ET and T had a considerable time shift; moreover, ET is slightly greater than T at the end of the day when the daily cumulative values are 2.7 mm and 2.55 mm for ET and T, respectively, where T ≈ 95.5% of ET.The amount of water stored in the stem introduces uncertainty when estimating diurnal whole-tree canopy transpiration and SF transpiration measured at the base of the stem.Motisi et al. [31] showed that this imbalance between canopy transpiration and tree water uptake is recovered during the night and is likely attributable to tree capacitance.Other studies have shown that, when the stem capacitance is large, an abundance of nighttime water recharges into the stem because whole-tree transpiration is higher than the stem q around mid-day [28,34,[36][37][38].
Since T exhibits almost a flat behavior close to its maximum value (Figure 4), surface energy balance outcomes can be directly measured using SF values if remote sensing images are acquired close to midday (without taking into account any time lag).It would be different if imaging and subsequent measurement were performed at the beginning or the end of the diurnal period (e.g., the 8:25 a.m.image acquisition of this experiment).Finally, uncertainties about using T measurements to validate the surface energy balance estimates from remote sensing arise from the individual shape of the two curves (Tp and ET).
The daily cumulative upscaled sap flow T was compared to the flux tower λETD for the day of remote sensing acquisition.The diurnal changes of T (Figure 5, black dots) and ET derived by FT (Figure 5, empty dots) show, again, the time lag between those different measurements.ET and T had a considerable time shift; moreover, ET is slightly greater than T at the end of the day when the daily cumulative values are 2.7 mm and 2.55 mm for ET and T, respectively, where T ≈ 95.5% of ET.Dichio et al. [39] pointed out an adaptive advantage for plants growing in semi-arid climates, such as olive trees.This adaptation allows them to maintain constant evapotranspiration rates even with low soil water content.Indeed, the high values of T/ET is consistent with the olive trees large capacity for osmotic adjustment under low soil water content [40], letting the plant keep high transpiration rates.
The forward time-shifting of T led to measuring the time lag between transpiration flux passing through the trunk section and latent heat flux through the canopy.In particular, the maximum determination coefficient took place for a time shift of ≈2.5 h (Figure 6), which can be assessed as a daily average time lag.This finding highlights that the measurement of remote sensing estimated through SF transpiration values should consider the time lag required by the sap to reach the canopy and finally evaporate.Thus, it should be considered canopy-specific.
The forward time-shifting of T led to measuring the time lag between transpiration flux passing through the trunk section and latent heat flux through the canopy.In particular, the maximum determination coefficient took place for a time shift of ≈2.5 h (Figure 6), which can be assessed as a daily average time lag.This finding highlights that the measurement of remote sensing estimated through SF transpiration values should consider the time lag required by the sap to reach the canopy and finally evaporate.Thus, it should be considered canopy-specific.Surely, T measurement can be used for validation if the soil evaporation is negligible, or if the validation of ET is on a daily scale base.Additionally, T seems less "noisy" and more stable compared to the FT measurement.As such, it has to be highlighted that, in this case, SF measurement needs to be accurately upscaled.
Since the average footprint that determine 70% of the fluxes, FP0.70, ranges between 60 and 110 m, and is much larger than the images' spatial resolution (i.e., original spatial resolutions were of 0.12 m and 0.50 m for multi-spectral and thermal images, respectively), spatially averaged λET values at the plant scale were compared to hourly transpiration T and its daily value (measured through SF).
The comparison was limited to 2 and 5 m spatial resolutions since fully vegetated pixels cannot be selected using 10 m spatial resolution images.During the experiment, the day cycle of Λ exhibited a typical concave-up shape (Figure 7) with an average daily value of 0.58 (if calculated between 5:00 a.m. and 5:00 p.m.).Indeed, Λ behaves according to the relative amplitude and phase of the H and λET harmonics [41].In particular, when no phase difference occurs between the turbulent fluxes, diurnal Λ exhibits a typical symmetric concave-up shape.However, when the phase difference occurs, the behavior becomes asymmetric.The Λ values are nearly constant during mid-day with two peaks near sunrise and sunset, which is when Λ deviates and increases sharply.Excluding the sunrise and sunset peaks, the average value was about 0.50.At 11:00 a.m., Λ has a minimum of 0.15, which might be attributed to several factors, including excessive solar radiation, deviation from the optimal air temperature, and/or very low soil water content.Surely, T measurement can be used for validation if the soil evaporation is negligible, or if the validation of ET is on a daily scale base.Additionally, T seems less "noisy" and more stable compared to the FT measurement.As such, it has to be highlighted that, in this case, SF measurement needs to be accurately upscaled.
Since the average footprint that determine 70% of the fluxes, FP 0.70 , ranges between 60 and 110 m, and is much larger than the images' spatial resolution (i.e., original spatial resolutions were of 0.12 m and 0.50 m for multi-spectral and thermal images, respectively), spatially averaged λET values at the plant scale were compared to hourly transpiration T and its daily value (measured through SF).
The comparison was limited to 2 and 5 m spatial resolutions since fully vegetated pixels cannot be selected using 10 m spatial resolution images.During the experiment, the day cycle of Λ exhibited a typical concave-up shape (Figure 7) with an average daily value of 0.58 (if calculated between 5:00 a.m. and 5:00 p.m.).Indeed, Λ behaves according to the relative amplitude and phase of the H and λET harmonics [41].In particular, when no phase difference occurs between the turbulent fluxes, diurnal Λ exhibits a typical symmetric concave-up shape.However, when the phase difference occurs, the behavior becomes asymmetric.The Λ values are nearly constant during mid-day with two peaks near sunrise and sunset, which is when Λ deviates and increases sharply.Excluding the sunrise and sunset peaks, the average value was about 0.50.At 11:00 a.m., Λ has a minimum of 0.15, which might be attributed to several factors, including excessive solar radiation, deviation from the optimal air temperature, and/or very low soil water content.According to Cammalleri et al. [36], comparing the stomata closure hourly behavior of different days, characterized by a regular Λ shape, revealed that the minimum value is not attributable to any stoma closure.
Thus, this minimum could be due to environmental factors that showed high variability during the morning.In particular, wind speed and direction changed abruptly at approximately 10:00-10:30 a.m., with the former increasing from 1.1 ms −1 on the average between 8:00 and 10:00 a.m. to 2.9 ms −1 on the average between 10:30 a.m. and 12:00 p.m.A direct consequence of the wind change is a According to Cammalleri et al. [36], comparing the stomata closure hourly behavior of different days, characterized by a regular Λ shape, revealed that the minimum value is not attributable to any stoma closure.
Thus, this minimum could be due to environmental factors that showed high variability during the morning.In particular, wind speed and direction changed abruptly at approximately 10:00-10:30 a.m., with the former increasing from 1.1 ms −1 on the average between 8:00 and 10:00 a.m. to 2.9 ms −1 on the average between 10:30 a.m. and 12:00 p.m.A direct consequence of the wind change is a difference in the footprint's size and position where turbulent fluxes originate.FP 0.70 increased from 32 to 98 m between the first and last acquisitions due to an increase in wind speed (Table 1).During the morning, FP 0.70 was 40 m between approximately 8.00 and 9:30 a.m., and it increased to 81 m between 10:30 a.m. and midday.Moreover, the average direction changed from the southwest (230 • clockwise, CW, from the north) to the west (279 • CW).The wind change affected the distance making the maximum contribution to the measured turbulent fluxes.This was 7 m in the first part of the morning and then increased up to 15 m during the second part.
The time required by the turbulent fluxes to reach the open path sensor from the source area range between 5 (fluxes originating at X max at 12.00 p.m., Table 1) and 40 s (fluxes originating at FP 0.70 at 8.00 a.m.).Thus, the quantified time shift should not be considered as an EC measurement artifact.

Time
Wind Speed Wind Direction X max FP 0.70

Diachronic Analysis
The 30 min weekly averaged diurnal behavior of λET, based on the R s approach (Figure 8, empty dots) and on the Λ approach (black dots), follows the same pattern of changes during the daytime (weekly average of the five selected weeks), as the hourly λET based on the Λ approach overestimated λET before 7:00 a.m. and after 4:00 p.m.However, the hourly λET based on the R s approach underestimated λET during the same times.By considering an acceptable variability (dashed lines) of ±15% of the EC λET βR (Figure 8, black line), the resulting best remote sensing acquisition time ranges between 10:00 a.m. and 4:00 p.m.
empty dots) and on the Λ approach (black dots), follows the same pattern of changes during the daytime (weekly average of the five selected weeks), as the hourly λET based on the Λ approach overestimated λET before 7:00 a.m. and after 4:00 p.m.However, the hourly λET based on the Rs approach underestimated λET during the same times.By considering an acceptable variability (dashed lines) of ±15% of the EC λETβR (Figure 8, black line), the resulting best remote sensing acquisition time ranges between 10:00 a.m. and 4:00 p.m. Coherently, the analysis of the temporal variability of CV (Figure 9) reveals that, between 8:30 a.m. and 2:00 p.m., half-hourly λET is representative of hourly λET (CV < 0.15).CV is lower (0.02) between 10:00 and 11:30 a.m., which aligns with findings from previous studies.

Spatial Resolution Analysis
Once the surface energy balance was applied, λET values were compared to those measured with the flux tower and added at the image acquisition times (194.8, 224.3, 189.7, and 253.3 Wm −2 at 8:23, 9:14; 10:00, and 10:18 a.m., respectively).ETD was evaluated through both the Λ self-preservation and Rs-based approaches at re-sampled spatial resolutions characterizing the canopy sizes and orchard structure (2, 5, and 10 m).Results show (Figure 10) that the second approach generally gives better results.approach underestimated λET during the same times.By considering an acceptable variability (dashed lines) of ±15% of the EC λETβR (Figure 8, black line), the resulting best remote sensing acquisition time ranges between 10:00 a.m. and 4:00 p.m. Coherently, the analysis of the temporal variability of CV (Figure 9) reveals that, between 8:30 a.m. and 2:00 p.m., half-hourly λET is representative of hourly λET (CV < 0.15).CV is lower (0.02) between 10:00 and 11:30 a.m., which aligns with findings from previous studies.

Spatial Resolution Analysis
Once the surface energy balance was applied, λET values were compared to those measured with the flux tower and added at the image acquisition times (194.8, 224.3, 189.7, and 253.3 Wm −2 at 8:23, 9:14; 10:00, and 10:18 a.m., respectively).ETD was evaluated through both the Λ self-preservation and Rs-based approaches at re-sampled spatial resolutions characterizing the canopy sizes and orchard structure (2, 5, and 10 m).Results show (Figure 10) that the second approach generally gives better results.

Spatial Resolution Analysis
Once the surface energy balance was applied, λET values were compared to those measured with the flux tower and added at the image acquisition times (194.8, 224.3, 189.7, and 253.3 W•m −2 at 8:23, 9:14; 10:00, and 10:18 a.m., respectively).ET D was evaluated through both the Λ self-preservation and R s -based approaches at re-sampled spatial resolutions characterizing the canopy sizes and orchard structure (2, 5, and 10 m).Results show (Figure 10) that the second approach generally gives better results.
When Λ does not exhibit regular behavior, the self-preservation approach provided overestimated retrievals.In particular, compared to the ET D, βR flux tower measurements (2.7 mm, surface energy balance closure), good alignments were found at 10:00 a.m.(2.83 mm, R s approach) and 10:18 a.m.(2.76 mm, Λ approach) at the 5 m spatial resolution.However, at the 2 and 10 m spatial resolutions, no good alignments were found, likely due to the application of a single source model.
The data do not fit the in situ SF measurements at lower spatial resolution (2 and 5 m), with the values closest to the daily flow (2.55 mm) at 9:14 a.m.(2.37 mm, R s approach) and 10:18 a.m.(2.76 mm, Λ approach) at 5 m of spatial resolution.When Λ does not exhibit regular behavior, the self-preservation approach provided overestimated retrievals.In particular, compared to the ETD,βR flux tower measurements (2.7 mm, surface energy balance closure), good alignments were found at 10:00 a.m.(2.83 mm, Rs approach) and 10:18 a.m.(2.76 mm, Λ approach) at the 5 m spatial resolution.However, at the 2 and 10 m spatial resolutions, no good alignments were found, likely due to the application of a single source model.
The data do not fit the in situ SF measurements at lower spatial resolution (2 and 5 m), with the values closest to the daily flow (2.55 mm) at 9:14 a.m.(2.37 mm, Rs approach) and 10:18 a.m.(2.76 mm, Λ approach) at 5 m of spatial resolution.

Absolute Difference Analysis (FT)
Flux tower observations were used to measure the absolute errors of λET and ETD.Errors affecting daily evapotranspiration were computed over both the self-preservation and the Rs techniques.Absolute errors within the range of ±0.15 (Figure 10, bands bounded by the dashed grey lines) were compatible to flux tower closure errors acceptability [28].
Acquisition at 8:23 a.m.led to worse results (Figure 11a), whereas acquisitions at 9:14 and 10:00 a.m.gave better estimates (at 2 and 5 m spatial resolutions).Integrating daily values based on the hypothesis of the constant evaporative fraction confirms the best results at 9:14 and 10:00 a.m.However, more accurate ADΛ values (Figure 11b) were found at 5 and 10 m spatial resolutions.Inaccuracies are probably due to the diurnal variability of evaporative fraction during the image acquisition day at 8:23 and 9:14 a.m.
More accurate daily evapotranspiration estimates are obtained through the RS approach (at all the acquisition times), with best estimates at 10:00 a.m. at 5 m spatial resolution (ADRs = 0.05).Even though error analysis cannot be considered conclusive, it seems to indicate the best acquisition time at 10:00 a.m. with a spatial resolution close to the average plant crown size.

Absolute Difference Analysis (FT)
Flux tower observations were used to measure the absolute errors of λET and ET D .Errors affecting daily evapotranspiration were computed over both the self-preservation and the R s techniques.Absolute errors within the range of ±0.15 (Figure 10, bands bounded by the dashed grey lines) were compatible to flux tower closure errors acceptability [28].
Acquisition at 8:23 a.m.led to worse results (Figure 11a), whereas acquisitions at 9:14 and 10:00 a.m.gave better estimates (at 2 and 5 m spatial resolutions).Integrating daily values based on the hypothesis of the constant evaporative fraction confirms the best results at 9:14 and 10:00 a.m.However, more accurate AD Λ values (Figure 11b) were found at 5 and 10 m spatial resolutions.Inaccuracies are probably due to the diurnal variability of evaporative fraction during the image acquisition day at 8:23 and 9:14 a.m.
More accurate daily evapotranspiration estimates are obtained through the R S approach (at all the acquisition times), with best estimates at 10:00 a.m. at 5 m spatial resolution (AD Rs = 0.05).Even though error analysis cannot be considered conclusive, it seems to indicate the best acquisition time at 10:00 a.m. with a spatial resolution close to the average plant crown size.When Λ does not exhibit regular behavior, the self-preservation approach provided overestimated retrievals.In particular, compared to the ETD,βR flux tower measurements (2.7 mm, surface energy balance closure), good alignments were found at 10:00 a.m.(2.83 mm, Rs approach) and 10:18 a.m.(2.76 mm, Λ approach) at the 5 m spatial resolution.However, at the 2 and 10 m spatial resolutions, no good alignments were found, likely due to the application of a single source model.
The data do not fit the in situ SF measurements at lower spatial resolution (2 and 5 m), with the values closest to the daily flow (2.55 mm) at 9:14 a.m.(2.37 mm, Rs approach) and 10:18 a.m.(2.76 mm, Λ approach) at 5 m of spatial resolution.

Absolute Difference Analysis (FT)
Flux tower observations were used to measure the absolute errors of λET and ETD.Errors affecting daily evapotranspiration were computed over both the self-preservation and the Rs techniques.Absolute errors within the range of ±0.15 (Figure 10, bands bounded by the dashed grey lines) were compatible to flux tower closure errors acceptability [28].
Acquisition at 8:23 a.m.led to worse results (Figure 11a), whereas acquisitions at 9:14 and 10:00 a.m.gave better estimates (at 2 and 5 m spatial resolutions).Integrating daily values based on the hypothesis of the constant evaporative fraction confirms the best results at 9:14 and 10:00 a.m.However, more accurate ADΛ values (Figure 11b) were found at 5 and 10 m spatial resolutions.Inaccuracies are probably due to the diurnal variability of evaporative fraction during the image acquisition day at 8:23 and 9:14 a.m.
More accurate daily evapotranspiration estimates are obtained through the RS approach (at all the acquisition times), with best estimates at 10:00 a.m. at 5 m spatial resolution (ADRs = 0.05).Even though error analysis cannot be considered conclusive, it seems to indicate the best acquisition time at 10:00 a.m. with a spatial resolution close to the average plant crown size.

Conclusions
This research focused on the use of the eddy covariance latent heat flux and sap flow transpiration to validate a surface energy balance model.Two time upscaling approaches to estimate daily evapotranspiration were tested at different spatial resolutions and acquisition times.Daily values were retrieved by assuming the self-preservation of the evaporative fraction in a diurnal cycle or by considering the proportionality with shortwave solar radiation at ground level.
Measures of transpiration via sap flow appear stable if compared with latent heat flux measured via a flux tower, even though sap flow measurements need to be accurately upscaled to manage transpiration in order to be used for validation.
The analysis confirmed that eddy covariance latent heat flux is suitable for the validation at both the instantaneous and the daily rate.A time lag was found between latent heat flux at the canopy and sap flow at the trunk; it lasted approximately 2.5 h within this study.
The time lag plays a role during validation since sap flow transpiration can be compared in a straightforward manner to surface energy balance estimates, if surface energy balance is applied on images acquired close to solar noon or if daily evapotranspiration is the aim of the validation.
Diachronic and spatial resolution analyses show that only a limited number of combinations can provide for accurate estimates of hourly and/or daily actual evapotranspiration.The upscaling approach based on diurnal evaporative fraction produces inaccuracies when estimating daily evapotranspiration, which is likely due to the irregular behavior of the evaporative fraction.However, the solar radiation technique yields better results at all acquisition times.On the whole, the analysis suggests a best acquisition time at ≈10:00 a.m. with a spatial resolution comparable to the average tree crown size (5 m).As such, the CV of EC λET reveals that half-hourly λET is representative of hourly values between 10:00 and 11:30 a.m.where CV <0.02.
Further investigations are required to determine how time lag varies with conductance reduction factors, including the plant's available soil water content.β R Bowen ratio (-)

Figure 1 .
Figure 1.(a) Geographical location of the experimental fields; the continuous black boundary has average coordinates of lat.: 37.6424°; long: 12.8464°; UTM WGS 84, which is close to the town of Castelvetrano, in the island of (b) Sicily, Italy.The blue dot and green square indicate the positions of FT and the average position of SF stations, respectively.

Figure 1 .
Figure 1.(a) Geographical location of the experimental fields; the continuous black boundary has average coordinates of lat.: 37.6424 • ; long: 12.8464 • ; UTM WGS 84, which is close to the town of Castelvetrano, in the island of (b) Sicily, Italy.The blue dot and green square indicate the positions of FT and the average position of SF stations, respectively.

Figure 2 .
Figure 2. The tethered helium balloon during the launching phase.In the inset (lower right corner), there is an image of the proximity-sensing platform, which includes an MCA-II multispectral camera, an A320G thermal camera, a mini PC, and a wireless Wi-Fi omnidirectional antenna.

Figure 2 .
Figure 2. The tethered helium balloon during the launching phase.In the inset (lower right corner), there is an image of the proximity-sensing platform, which includes an MCA-II multispectral camera, an A320G thermal camera, a mini PC, and a wireless Wi-Fi omnidirectional antenna.
6.0 (Media Cybernetics Inc.: Rockville, Maryland, USA) to define the sapwood depth.The values of LAI and LAIp were estimated from in situ observations collected by means of a Licor LAI-2000 Plant Canopy Analyzer (LI-COR ® Biosciences): LAI was ≈0.89 m 2 •m −2 , whereas LAIp was ≈1.27 m 2 •m −2 , on average.

Figure 4 .
Figure 4.The 30 min normalized upscaled Tp measured with SF (black dots) and 30 min normalized diurnal ET measured with the EC tower (empty dots) averaged over a week.Data points represent the three-tree average transpiration averaged over a week.Nocturnal values of normalized EC λET are not shown because of the inaccurate measurements of small fluxes.

Figure 4 .
Figure 4.The 30 min normalized upscaled Tp measured with SF (black dots) and 30 min normalized diurnal ET measured with the EC tower (empty dots) averaged over a week.Data points represent the three-tree average transpiration averaged over a week.Nocturnal values of normalized EC λET are not shown because of the inaccurate measurements of small fluxes.

Figure 4 .
Figure 4.The 30 min normalized upscaled Tp measured with SF (black dots) and 30 min normalized diurnal ET measured with the EC tower (empty dots) averaged over a week.Data points represent the three-tree average transpiration averaged over a week.Nocturnal values of normalized EC λET are not shown because of the inaccurate measurements of small fluxes.

Figure 5 .
Figure 5.The daily cumulative upscaled SF Tp and EC ET during the image acquisition day.Each data point of transpiration represents the upscaled average of three trees.

Figure 5 .
Figure 5.The daily cumulative upscaled SF Tp and EC ET during the image acquisition day.Each data point of transpiration represents the upscaled average of three trees.

Figure 6 .
Figure 6.The coefficient of determination, r 2 , between EC evapotranspiration and time-shifted upscaled SF transpiration during the image acquisition day.

Figure 6 .
Figure 6.The coefficient of determination, r 2 , between EC evapotranspiration and time-shifted upscaled SF transpiration during the image acquisition day.

Figure 7 .
Figure 7.Diurnal Λ behavior during the image acquisition day.

ETD (mm•d − 1 )Figure 8 .
Figure 8.The weekly 30 min averaged diurnal changes of (black dots) and (empty dots) with the ETD measured with EC tower and the range of variability ±15%.

Figure 9 .
Figure 9.Half hourly temporal variability of the coefficient of variation, CV (unitless), evaluated over the time from 29 June to 2 August.The over-imposed dash line indicates a 15% limit of acceptability.

Figure 8 .
Figure 8.The weekly 30 min averaged diurnal changes of (black dots) and (empty dots) with the ET D measured with EC tower and the range of variability ±15%.

ETD (mm•d − 1 )Figure 8 .
Figure 8.The weekly 30 min averaged diurnal changes of (black dots) and (empty dots) with the ETD measured with EC tower and the range of variability ±15%.

Figure 9 .
Figure 9. Half hourly temporal variability of the coefficient of variation, CV (unitless), evaluated over the time from 29 June to 2 August.The over-imposed dash line indicates a 15% limit of acceptability.

Figure 9 .
Figure 9. Half hourly temporal variability of the coefficient of variation, CV (unitless), evaluated over the time from 29 June to 2 August.The over-imposed dash line indicates a 15% limit of acceptability.

Figure 10 .
Figure 10.Spatially averaged ETD on the footprint area, calculated using the integration factors based on Λ (grey bars) and Rs (black bars), during the image acquisition day for different input spatial resolutions: 2 m (a); 5 m (b); 10 m (c).The dashed lines indicate the ETD measured via the FT (2.7 mm, black line) and the T24 measured via the SF method (2.55 mm, grey line).

Figure 11 .
Figure 11.Absolute difference, AD (unitless), of λET (a), and ETD based on Λ scaling factor (b), and Rs scaling factor (c) during the image acquisition day.The two dashed grey lines indicate an acceptable range of variability of ±15%.

Figure 10 .
Figure 10.Spatially averaged ET D on the footprint area, calculated using the integration factors based on Λ (grey bars) and R s (black bars), during the image acquisition day for different input spatial resolutions: 2 m (a); 5 m (b); 10 m (c).The dashed lines indicate the ET D measured via the FT (2.7 mm, black line) and the T24 measured via the SF method (2.55 mm, grey line).

Figure 10 .
Figure 10.Spatially averaged ETD on the footprint area, calculated using the integration factors based on Λ (grey bars) and Rs (black bars), during the image acquisition day for different input spatial resolutions: 2 m (a); 5 m (b); 10 m (c).The dashed lines indicate the ETD measured via the FT (2.7 mm, black line) and the T24 measured via the SF method (2.55 mm, grey line).

Figure 11 .
Figure 11.Absolute difference, AD (unitless), of λET (a), and ETD based on Λ scaling factor (b), and Rs scaling factor (c) during the image acquisition day.The two dashed grey lines indicate an acceptable range of variability of ±15%.

Figure 11 .
Figure 11.Absolute difference, AD (unitless), of λET (a), and ET D based on Λ scaling factor (b), and R s scaling factor (c) during the image acquisition day.The two dashed grey lines indicate an acceptable range of variability of ±15%.
ET D actual daily evapotranspiration (mm•d −1 ) ET D, βR ET D adjusted by preserving β R (mm•d −1 ) ET D,Rs ET D based on R s as integration factor (mm•d −1 ) ET D,Λ ET D based on the self-preservation hypothesis of Λ (mm•d −1 ) FP 0.70 average footprint determining 70% of the fluxes (m) G 0 soil heat flux at ground level at the acquisition time (W•m −2 ) H sensible heat flux (W•m −2 ) LAI leaf Area Index (m 2 •m −2 ) LAI p LAI of a single plant (m 2 •m −2 ) R n net radiation at the acquisition time (W•m −2 ) R n,D daily net radiation (W•m −2 ) R s hourly incoming shortwave radiation (W•m −2 ) S the cross-sectional area of conducting sapwood (m 2 ) T hourly transpiration of an entire field (mm•d −1 ) T p hourly transpiration of a single plant (mm•d −1 ) X max distance determining the maximum flux contribution (m) α multiplicative coefficient of the sap flow Equation (5) (m•s −1 )