Next Article in Journal
Real-Time Flood Control by Tree-Based Model Predictive Control Including Forecast Uncertainty: A Case Study Reservoir in Turkey
Next Article in Special Issue
A Simplified Nitrogen Assessment in Tagus River Basin: A Management Focused Review
Previous Article in Journal
Groundwater Chemistry Regulated by Hydrochemical Processes and Geological Structures: A Case Study in Tongchuan, China
Previous Article in Special Issue
Basin Irrigation Design with Multi-Criteria Analysis Focusing on Water Saving and Economic Returns: Application to Wheat in Hetao, Yellow River Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Implementation of a Two-Source Model for Estimating the Spatial Variability of Olive Evapotranspiration Using Satellite Images and Ground-Based Climate Data

1
Centro de Investigación y Transferencia en Riego y Agroclimatología (CITRA) and Research Program on Adaptation of Agriculture to Climate Change (A2C2), Universidad de Talca, Casilla 747, Talca 3460000, Chile
2
Faculty of Natural Resources, Universidad Católica de Temuco, Temuco 4780000, Chile
*
Author to whom correspondence should be addressed.
Water 2018, 10(3), 339; https://doi.org/10.3390/w10030339
Submission received: 6 November 2017 / Revised: 14 March 2018 / Accepted: 16 March 2018 / Published: 19 March 2018
(This article belongs to the Special Issue Innovation Issues in Water, Agriculture and Food)

Abstract

:
A study was carried out to evaluate the potential use of the two-source Shuttleworth and Wallace (SW) model to compute the intra-orchard spatial variability of actual evapotranspiration (ET) of olive trees using satellite images and ground-based climate data. The study was conducted in a drip-irrigated olive orchard using satellite images (Landsat 7 ETM+), which were acquired on clear sky days during the main phenological stages (2009/10 growing season). The performance of the SW model was evaluated using instantaneous latent heat flux (LE) measurements that were obtained from an eddy correlation system. At the time of satellite overpass, the estimated values of net radiation ( Rn i ) and soil heat flux ( G i ) were compared with ground measurements from a four-way net radiometer and soil heat flux plates, respectively. The results indicated that the SW model subestimated instantaneous LE (W m−2) and daily ET (mm d−1), with errors of 12% and 10% of observed values, respectively. The root mean square error (RMSE) and mean absolute error (MAE) values for instantaneous LE were 26 and 20 W m−2, while those for daily values of ET were 0.31 and 0.28 mm d−1, respectively. Finally, the submodels computed Rn i and G i with errors of between 4.0% and 8.0% of measured values and with RMSE and MAE between 25 and 39 W m−2.

1. Introduction

Better water management of irrigation is required to optimize the water productivity of olive oil production due to water scarcity. Under this scenario, sophisticated irrigation water management will be required to maintain sufficient levels of productivity and quality [1,2]. For these objectives, it is necessary to have an accurate estimation of daily actual evapotranspiration (ET) that is generally computed as a function of reference evapotranspiration ( ET o ) and crop coefficient (Kc) [3]. Several researchers also suggested using the dual crop coefficient approach to describe the ratio of ET to ET o by separating Kc into basal crop coefficient (Kcb) and soil evaporation coefficient (Ke) [3,4,5]. However, the values of Kc and Kcb reported in the literature for heteronomous canopies require local adjustment because they depend on canopy architecture and non-linear interaction of soil, cultivar, and climate [1,6,7]. For a hedge-pruned olive orchard, Martínez-Cob and Faci [8] indicated that Kc values depend on the geometric characteristics of the canopy (canopy shape, distance between trees, etc.) and fractional cover (fc). For a super intensive olive orchard, Paço et al. [7] suggested that the Kc value is affected by several factors, including the canopy architecture, fractional cover, crop management practices, and rainfall variability. For a drip-irrigated olive orchard, López-Olivari et al. [9] indicated that the ratios of transpiration and evaporation to ET o could be significantly affected by the irrigation systems, which determine the percentage of wetted area (Aw) of the soil surface. Allen and Pereira [10] and Paço et al. [7] recommended adjusting the ET/ ET o ratios using fc and tree height, according to a density coefficient (Kd) that describes the increase in Kc with an increase in canopy size. Pôças et al. [11] for olive orchard improved the estimation of crop coefficients using remotely sensed vegetation indices and the SIMdualKc soil water balance model. Cammalleri et al. [6] also evaluated a remote sensing-based approach to estimate olive ET, combining a modified version of the standard FAO-56 dual crop coefficient procedure and Penman-Monteith (PM) equation with actual canopy characteristics (i.e., leaf area index, albedo, and canopy height) that were derived from optical remote sensing data. Finally, Er-Raki et al. [12] evaluated the potential of assimilating ET derived from satellite thermal infrared observations to improve the ET simulation performances from the FAO-56 single crop coefficient approach.
Several studies indicated that the two-layer model of Shuttleworth and Wallace (SW) could be used to estimate ET directly without using crop coefficients [13,14,15,16,17]. The SW model can estimate soil evaporation and transpiration separately and has been widely used to estimate ET for homogeneous (sorghum, corn, and wheat) and heterogeneous (vineyards and orchards) canopies [14,18,19,20,21]. For heterogeneous canopies, such as vineyards and orchards, the SW approach has been able to estimate ET with errors ranging between 6% and 25%. In a furrow-irrigated vineyard (fc = 0.35), Zhang et al. [20] indicated that the SW model overestimated latent heat flux (LE) with an index of agreement ( I a ) = 0.75 and a mean absolute error (MAE) = 39 W m−2, while Zhang et al. [22] noted that the SW model overestimated ET by 25.2% in a furrow-irrigated vineyard. For drip-irrigated vineyards (fc = 0.30), Ortega-Farías et al. [14] observed that the SW model was able to predict ET with a root-mean-square error (RMSE) = 0.51 mm d−1 and an MAE = 0.41 mm d−1. In a furrow-irrigated vineyard (fc = 0.30), Zhao et al. [23] observed that the SW model estimated ET with a RMSE and MAE of 0.68 and 0.52 mm d−1, respectively. In a drip-irrigated olive orchard, Ortega-Farías and Lopez-Olivari [24] observed that the SW model overestimated LE by approximately 2% of measured values, with a RMSE = 28 W m−2 and I a = 0.98. This study also indicated that the SW model was very sensitive to errors in the values of stomatal resistance ( r st ) and leaf area index ( LAI ). On this matter, Brenner and Incoll [25] indicated that the SW model overestimated crop transpiration of sparsely vegetated shrublands due to the over-estimation of the radiation that is absorbed by the canopy.
The traditional approach of SW model does not consider the effect of intra-orchard spatial variability of soil and canopy vigor on the estimation of ET. Canopy covers of commercial olive orchards are generally incomplete as a result of the canopy geometry that depends on canopy size, leaf area index, and plant density [26]. The canopy training system and associated canopy geometry may significantly affect the partitioning of Rn into LE, H, or G and thus orchard water requirements. To include the effect of intra-field spatial variability into the SW model, the instantaneous values of net radiation ( Rn i ) and soil heat flux ( G i ) over orchards can be estimated using remote sensing images. In this regard, several studies have suggested that Rn i and G i can be adequately estimated in heterogeneous canopies using satellite images. In a drip-irrigated olive orchard, Ortega-Farías et al. [27] indicated that Rn i was estimated with an error of 3% of measured values when using thermal and multispectral images from the Landsat platform. In a drip-irrigated vineyard, Carrasco-Benavides et al. [28] using Landsat 7 images observed that the METRIC model estimated Rn i and G i with an error of 11% and 5% of measured values, respectively. However, the estimation of ET from the METRIC model requires the selection of the two “anchor” pixels, which is subjective and depends on the ability of the operator in search and isolate of the most appropriate hot and cold pixels [29,30,31]. For drip-irrigated olive orchards (fc < 0.3), the selection of the hot and cold pixels is critical because sensible heat flux (H), which is produced at the soil surface, is the main component of energy balance and plays a key role in the tree transpiration and stomatal closure [24].
Dhungel et al. [32] successfully used remote sensing tools to facilitate a comparison of different parameterizations to estimate LE over sparse canopies. This comparison was conducted using gridded weather and Landsat satellite data, such as instantaneous values of ET, surface roughness (Zom), emissivity ( ε 0 ), albedo ( α ), and leaf area index ( LAI ), to parameterize the Penman-Monteith (PM) equation. According to our knowledge, there is no available information about the estimation of the intra-orchard spatial variability of olive ET using ground-based weather measurements and remote sensing images as inputs to the SW model. Thus, the main objective of this study is to evaluate the potential use of the two-source SW model to compute the intra-orchard spatial variability of ET of a drip-irrigated olive orchard using satellite images and ground-based climate data. Moreover, the sub-models that were used to estimate the instantaneous values of latent heat flux ( LE i ), available energy ( A i ), net radiation ( Rn i ) , and soil heat flux ( G i ) were also evaluated for the main phenological growth stages of olive trees.

2. Materials and Methods

2.1. Study Site Description

Measurements were collected in a drip-irrigated commercial olive orchard (Olea europaea L. cv Arbequina) that was located in Pencahue, Maule Region, Chile (35°23′ S, 71°44′ W, WGS84, 96 m above sea level) during the 2009–2010 season. The study area has a Mediterranean climate with a medium average temperature of 14.8 °C and an accumulated ET o of 1013 mm. The average annual rainfall in the region is approximately 602 mm, which is falling mainly during the winter (May to September). The summer period (December to March) is generally hot and dry with a high water vapor pressure deficit. The soil at the experimental site was classified as Quepo series with a clay loam soil. The trees were established using a hedgerow system with 5.0 m between rows and 1.5 m between trees (1333 trees ha−1). The orchard was irrigated by two drippers per tree (2.0 L h−1) spaced at intervals of 0.75 m along the rows. After irrigation, the percentage of wetted area ( A W ) by the drippers located under the tree canopy was 4.5% of the total area.
To evaluate the irrigation management, the midday stem water potential (Ψmd) was measured using a pressure chamber (PMS Instruments Co., Model 1000, Corvallis, OR, USA). A chosen shoot (two per tree, one tree per replicate) containing five to six pairs of leaves was encased in a plastic bag and then wrapped in aluminum foil for (at least) 2 h before being cut [33,34]. In addition, a LI-COR gas analyzer (Li-6400, LI-COR Inc., Lincoln, NE, USA) was used to measure stomatal resistance ( r st ) on two leaves that were directly exposed to the sun and were located on the mid-section of the tree [35]. The Ψmd and r st measurements were collected from exposed leaves on 10 trees that were located on both sides of the olive rows. Finally, the LAI was measured weekly using a plant canopy analyzer (LAI-2000, LI-COR, Lincoln, NE, USA), which was calibrated according to Lopez-Olivari et al. [9]. Also, it is important to indicate that Ψmd and r st were measured at the time of satellite overpass.

2.2. Measurements of Climate and Energy Balance Data

Two towers (4.8 m height) were installed on a flat and homogeneous plot (6.45 ha) in an olive orchard to measure the surface energy balance components and micrometeorological variables at intervals of 30 min. Hobo sensors (Onset Computer, Inc., Bourne, MA, USA) were used to measure relative humidity ( RH ) and air temperature ( Ta ). An anemometer (03101-5, R. M. Young Co., Traverse City, MI, USA) and a pluviometer (A730RAIN, Adcon Telemetry, Klosterneuburg, Austria) were used to measure wind speed ( u ) and precipitation (Pp), respectively. The net radiation (Rn) was measured with a four-way net radiometer (CNR1, Kipp & Zonen, Delft, The Netherlands). Finally, the sensors to measure u, Pp, Ta, RH, and Rn were located at 1.9 m above the tree canopy.
Soil heat flux was estimated through the use of eight flux plates installed on either side of the rows (four plates in the inter row and four plates below the row). This arrangement considers the row shade effect throughout the day [36]. The flux plates with a constant thermal conductivity (HFT3, Campbell Sci., Logan, UT, USA) were installed at a depth of 0.08 m. Additionally, two averaging thermocouple probes (TCAV, Campbell Sci., Logan, UT, USA) that were used to measure soil temperature were installed above each flux plate at depths of 0.02 and 0.06 m. All of the thermocouple probe signals were recorded on an electronic datalogger (CR3000, Campbell Sci, Logan, UT, USA) with a thirty-minute interval. Finally, soil heat flux was calculated at each position by adding the measured flux at 0.08 m to the heat that was stored in the layer above the heat flux plates [9,37].
The latent (LE) and sensible (H) heat fluxes were measured with an eddy covariance (EC) system, which is composed of a three-dimensional sonic anemometer (CSAT3, Campbell Scientific Inc., Logan, UT, USA) and an infrared gas analyzer (LI-7500, LI-COR Inc., Lincoln, NE, USA). The fluxes were registered at 1.6 m above the canopy at intervals of 10 Hz and were stored on a data logger (CR5000, Campbell Scientific Inc., Logan, UT, USA). Finally, a post-processing step was conducted to correct the sonic temperature due to crosswind influences [38] and water vapor density due to the influences of the fluctuations in temperature and humidity [39].
To reduce the uncertainty that is associated with the errors in the LE and H measurements, entire days were excluded from the study when the ratios of ( H + LE ) to ( Rn G ) were outside the range between 0.8 and 1.2 [9]. Assuming that the measurements of Rn and G were representative of the available energy over the olive orchard, the fluxes of H and LE were forced to close the energy balance using the Bowen ratio approach ( B = H / LE ) [8,24,36,40]:
LE B = ( Rn G ) ( 1 + B )
H B = ( Rn G ) ( 1 + B 1 )
On a daily basis, the olive ET was calculated, as follows:
ET EC = n = 1 24 LE B   n λ ρ w 1.8
where ET EC is the actual evapotranspiration of the olive orchard measured by the EC system (mm d−1), 1.8 is a conversion factor, λ is the latent heat of vaporization (1013 MJ kg−1), ρ w is the water density (1000 kg m−3), and n is the number of measurements during a 24-h period. The subscript B indicates that turbulent fluxes were recalculated using the Bowen ratio approach.

2.3. Shuttleworth and Wallace Model Description

At the time of satellite overpass, the SW model can estimate latent heat flux over olive orchards as the sum of the Penman-Monteith equation for evaporation and transpiration weighted by a set of coefficients that account for the combination of soil and canopy resistances. This can be done as follows [41]:
LE i = T i + E i
T i = C c Δ A i + ( ρ a C p D i Δ r a c A si r a a + r a c ) Δ + γ ( 1 + r s c ( r a a + r a c ) )
E i = C s Δ A i + ( ρ a C p D i Δ r a s ( A i A si ) r a a + r a s ) Δ + γ ( 1 + r s s ( r a a + r a s ) )
where LE i is the instantaneous latent heat flux computed from the SW model (W m−2), T i is the latent heat flux from the tree transpiration (W m−2), E i is the latent heat flux from soil evaporation (W m−2), C C is the canopy resistance coefficient (dimensionless), C S is the soil surface resistance coefficient (dimensionless), Δ is the slope of the saturation vapor pressure curve at the mean temperature (kPa °C−1), A i is the available energy leaving the complete canopy (W m−2), A Si is the available energy at the soil surface (W m−2), C P is the specific heat of the air at a constant pressure (1013 J kg−1 K−1), D i is the water vapor pressure deficit at the reference height (kPa), r a c is the bulk boundary layer resistance of the vegetative elements in the canopy (s m−1), r a a is the aerodynamic resistance between the canopy source height and reference level (s m−1), γ is the psychrometric constant (kPa °K−1), r s c is the canopy resistance (s m−1), r a s is the aerodynamic resistance between the soil and canopy source height (s m−1), and r s s is the soil surface resistance (s m−1). Subscripts “i” denote the values that are computed at the time of satellite overpass.

2.4. Estimation of Resistances

Values of r s c and r a c are calculated as:
r s c = r st LAI
r a c = r b LAI
where r st is the mean stomatal resistance (s m−1), r b is the mean boundary layer resistance (s m−1) and LAI is the leaf area index (m2 m–2). For this study, r b was assumed to be equal to 25 s m−1, and a constant value of r s s = 2000 sm−1 was used because of the dryness of the soil surface (between rows) during the study periods [41]. Additionally, a mean value of measured r st = 235 (±60.84) s m−1 was used for the simulation period. Finally, general descriptions of r a a , r a s , C s , and C c can be found in Shuttleworth and Wallace [41] and Ortega-Farías et al., 2007 [21].

2.5. Estimation of Available Energy at the Time of the Satellite Overpass

The instantaneous values of available energy over the canopy ( A i ) and soil surface ( A Si ) can be calculated as:
A i = Rn i G i
A si = Rn Si G i
where Rn i is the net radiation (W m−2), G i is the soil heat flux (W m−2) and Rn Si is the net radiation on the surface of the ground (W m−2), which can be calculated using Beer’s law, through the next expression:
Rn Si = Rn i × exp ( C × LAI )
where C is the extinction coefficient of the net radiation of the crop (dimensionless). In this study, a mean value of measured LAI = 1.29 (±0.07) was used as a constant while C was assumed to be 0.66 [24]. Values of Rn i were obtained, as follows [42]:
Rn i = ( 1 α i ) × R s + R L R L ( 1 ε 0 ) × R L
where α i is the surface albedo, R s is the short-wave incident radiation (W m−2), R L and R L are the incoming and outgoing long-wave radiation, respectively, and ε 0 corresponds to the superficial emissivity. The values of R L and R L were obtained as:
R L = ε 0 × σ × T s 4
R L = ε a × σ × T a 4
where σ is the Stefan-Boltzmann constant (5.67 × 10−8 W m−2 K−4), T s is the superficial temperature (°K), ε 0 is the “broad-band” surface emissivity (dimensionless), T a is the near-surface air temperature (°K), and ε a is the effective atmospheric emissivity (dimensionless). The methodologies to estimate ε 0 ,   ε a ,   and   T s are indicated by Allen et al. [42].
Gi was calculated using the following empirical relation [26]:
G i = 0.3236 × Rn i 51.52

2.6. Instantaneous to Daily Extrapolation of Actual Evapotranspiration

Using multispectral and thermal images from the Landsat satellite and ground-based weather data, the daily ET was computed for each image pixel using the following equations [28,42]:
ET i = 3600 LE i ρ w λ
E F = ET i ET oi
ET sw = E F × ET o
where ET i is the instantaneous ET computed from the SW model (mm h−1), ρ w is the density of water (1000 kg m−3), λ is the latent heat of vaporization (J kg−1), ET oi is the instantaneous reference evapotranspiration (mm h−1), ET sw is the daily actual evapotranspiration estimated using the SW model (mm d−1), ET o is the daily reference evapotranspiration (mm d−1), and E F is the reference evapotranspiration fraction.

2.7. Image Processing

Seven Landsat (7 ETM+) satellite images were downloaded from the USGS GloVis web platform for free (http://glovis.usgs.gov) using the coordinates WRS-2 Path/Row: 233/85 (Table 1). The SLC instrument suffered a problem in 2003; thus, the images have lines of pixels without information. A method that was proposed by Storey et al. [43] was implemented to fill the pixels without data. This method estimates the value of the pixel based on the neighboring pixels in an image taken on a similar date. The estimations present good results in most of the images, but not in images with abrupt temporal transitions [44]. Also, the image processing for satellite information considered the methodology that was proposed by Allen et al. [42], aiming to generate radiometric and atmospheric corrections.
Multispectral and thermal data from each satellite image were processed pixel by pixel to estimate LE i , Rn i , G i , and ET i . In this regard, the satellite images were obtained from November to February where the main phenological periods were observed in the drip-irrigated olive orchard (Table 1). Finally, it is important to indicate that the number of days with a complete data set (satellite images and EC fluxes) was limited by cloudiness when the satellite overpassed the experimental site and persistent noise of the EC system (instrumental problems and flow distortion through the tower). For these reasons, several researchers have used a limited number of images to evaluate the satellite-based remote sensing (SBRS) models [6,45,46,47,48,49].

2.8. Statistical Analysis

The SW model used to compute the instantaneous latent heat flux ( LE i ) and daily actual evapotranspiration ( ET sw ) was validated using ground-based weather measurements. The validation was carried out using the ratio of the estimated to observed values (b), root-mean-square error (RMSE), mean absolute error (MAE), and index of agreement ( I a ) [50,51]. The t-test was used to determine whether the b value was significantly different from unity at the 95% confidence level. In addition, the sub-models for estimating A i ,   Rn i , and G i were included in the validation. Values of RMSE, MAE, and I a were computed, as follows:
RMSE = [ N 1 i = 1 N ( P i O i ) 2 ] 0.5
MAE = N 1 i = 1 N | P i O i |
I a = 1 [ i = 1 N ( P i O i ) 2 i = 1 N ( | P i O ¯ | + | O i O ¯ | ) 2 ]   0 I a 1
where N is the total number of observations, P i and O i are the estimated and observed values, respectively, and O ¯ is the mean of the observed values. Values of RMSE, MAE, P i , O i , and O ¯ are in W m−2 or mm d−1
Finally, a sensitivity analysis of the inputs parameters (C, d, LAI, n) and resistances ( r b ,   r st ,   r a a ,   r a c ,   r a s ) was conducted to evaluate their effects on the ability of the SW model for estimating ET i . In this case, the percent deviation of the mean ET i was computed when the input values of the parameters and resistances individually varied by +/− 30%.

3. Results

The atmospheric conditions were very dry and hot during the main phenological periods (FS, PH, and FC) of the olive orchard, which were observed from December to February (Table 1). During this period, the values of u, Ta and D were between 0.3 and 4.5 m s−1, 28.3–13.4 °C, and 0.2–2.9 kPa, respectively (Figure 1a,b). During the study, four rainfall events were observed with values of less than 5 mm. The accumulated value of ET o from September to March (growing season) was 907 mm while that from December to March (summer) was 616 mm. During this study, the average values of ET o and ET oi were 6.95 (±0.99) mm d−1 and 0.58 (±0.06) mm h−1, respectively (Table 2). These results indicate that the drip-irrigated olive orchard was under high atmospheric demand for water vapor (Figure 1). Under this atmospheric condition, values of Ψmd ranged between −1.37 and −1.49 MPa, indicating that olive trees were well irrigated during the main phenological periods (F, FS, PH, and FD). Finally, Figure 1 shows that clear days were observed when the satellite passed over the experimental site.
At 30 min time intervals, the ratio of ( Rn G ) to ( H + LE ) was 0.89, suggesting that the orchard energy balance (SEB) was systematically imbalanced by approximately 11% (Figure 2). Several researchers reported a similar situation where the values of (LE + H) were less than those of (Rn − G) above olive orchards when the EC technique was used. In flood-irrigated olive orchards, Williams et al. [52], Ezzahar et al. [53], and Er-Raki et al. [4] reported imbalances ranging between 9% and 26%, while in drip-irrigated olive orchards, Villalobos et al. [54] and Testi et al. [55] indicated imbalances varying between 5% and 17%. The potential errors can be attributed to the uncertainties in the measurements of Rn and G [56,57,58], and energy storage within the olive tree biomass [9,53]. In this case, Twine et al. [40] suggested that LE and H values from de EC system can be corrected using the Bowen ratio (β = H/LE) because the different problems affect the measured values of H and LE in a similar proportion. Using this correction, however, Allen et al. [59] indicated that the lack of closure can be only attributed to errors in the measurements of H and LE without considering the potential bias in Rn-G. According to several researchers, the uncertainties observed in this study are modest and turbulent fluxes were recalculated using the Bowen-ratio approach [8,24,36,40,60].

3.1. Model Validation of Available Energy, Latent Heat Flux and Actual Evapotranspiration

The model validation of Rn i , G i , and A i is indicated in Table 3 when using ground-based weather measurements and satellite images. The results indicated that the sub-model was able to predict values of A i with a RMSE = 25 W m−2 and MAE = 21 W m−2. The RMSE and MAE values for Rn i were between 39 and 32 W m−2, while those for G i were between 33 and 27 W m−2, respectively (Table 3). The t-test showed that the b values were different from unity, indicating that A i ,   Rn i , and G i were underestimated, with errors ranging between 4% and 8% of observed values in the drip-irrigated olive orchard. Finally, the comparison between the measured and modeled values in Figure 3 shows that all points were close to the 1:1 line, suggesting a good performance of the sub-models to estimate Ai, Rn i , and G i . For Rn i and G i , the results that were observed in this study are similar to those that are found in literature for heterogeneous canopies, such as orchards and vineyard. For winegrape (drip irrigated) and table grape (micro sprinkler), Teixeira et al. [61] observed values of coefficient of determination (R2) = 0.94 and RMSE = 17.5 W m−2 for Rn i when using satellite images (Landsat 5 and 7). For a drip-irrigated olive orchard (fc = 0.29), Ortega-Farías et al. [27] found that the METRIC model underestimated Rn i by approximately 3% of the observed values, with RMSE and MAE values of 40 and 33 W m−2, respectively. In a super intensive olive orchard (fc = 0.3), Ortega-Farías et al. [26] indicated that the model was able to estimate Rn i and G i , with errors of 5% and 2%, respectively, when using high-resolution thermal and multispectral data acquired with an unmanned aerial vehicle (UAV). In a drip-irrigated vineyard (fc = 0.3), Carrasco-Benavides et al. [28] noted that the METRIC model was able to estimate Rn i and G i , with errors of 11% and 5%, respectively. In an experiment that was carried out using remote sensing techniques over heterogeneous landscapes, Liu et al. [62] found RMSE and MAE values of 51 and 25 W m−2, respectively, when G i was evaluated.
Table 3 indicates that the SW model underestimated LE i by approximately 12% of observed values, with RMSE = 26 W m−2 and MAE = 20 W m−2 (Table 3). In addition, the comparisons between the observed and estimated values of LE at the time of satellite overpass show that most of the points are close to the 1:1 line (Figure 3). Furthermore, the model validation indicated the SW model was able to simulate the ET with RMSE and MAE values that are equal to 0.31 and 0.28 mm d−1, respectively (Table 3). Additionally, the statistical analysis indicated that the ratio of ET sw / ET EC was significantly different from unity, suggesting that the SW model underestimated ET with an error of 10% of observed values. In this regard, Figure 4 shows that all of the points were distributed below the 1:1 line. In a super intensive olive orchard (fc = 0.3, LAI = 1.3), Ortega-Farías and Lopez-Olivari [24] indicated that the SW model overestimated LE i and ET by approximately 2% and 6% of observed values, respectively, suggesting that the model was very sensitive to errors in the estimation of r st and LAI. Finally, in a furrow-irrigated vineyard (fc = 0.35, IAF = 2), Zhang et al. [20] observed MAE values about of 39 W m−2 when SW model was evaluated.

3.2. Sensitivity Analysis

Table 4 indicates that when Rni varied by ±30%, a variation of ±20% was observed in the estimation of ET i . Also, the results indicated that the sensitivity of the predicted ET i in respect to the uncertainties in C ,   d ,   n ,   and   r b was minimal. In addition, the sensitivity analysis indicated that the estimation of ET i was significantly affected by variations of ±30% in the values of LAI and r s c . In this regard, the values of ET i varied between +15% and −14% and between −13% and +20% when LAI and r s c varied by ±30%, respectively. (Table 4). This sensitivity analysis suggests that the ET i that was estimated using the SW model was very sensitive to Rni, LAI, and r st , but it was not sensitive to errors in the estimations of r a a , r a c , and r a s . Similar results in a drip-irrigated olive orchard were observed by Ortega-Farías and Lopez-Olivari [24], who indicated that the ET predicted by the SW model was sensitive to errors of ±30% in LAI and r st , but was not significantly affected by errors in the estimation of aerodynamic resistances. Therefore, correct estimations of Rn i , LAI and r st become critical to increasing the accuracy of the SW model.

3.3. Spatial Variability of Actual Evapotranspiration

For flowering, pit hardening and fruit colouring, maps of the actual evapotranspiration were generated using the temporal variability of ET sw within the drip-irrigated olive orchard (Figure 5). While differences were observed in the ET sw values over the entire time series (from 1.34 to 3.69 mm), the spatial variability for each of these days was low, with a maximum variation between 2.85 and 3.07 mm d−1 (DOY 22). These low variations of the spatially distributed ET sw maps could be associated with the low intra-orchard spatial variability of A i (Figure 6), which presented a standard deviation between 3 and 6 W m−2 for the entire study period. The highest value of A i was observed on DOY 339 when the A i ranged between 476 and 494 W m−2. On the other hand, when we consider the daily variation, the results reflected a low spatial variability of the available energy. These results are logical when considering that the orchard presents similar management conditions (irrigation practices, training systems, canopy management, fertilization, among others), which are expressed by the low spatial variability at the experimental site.

4. Final Remarks

The results that were obtained in this study were similar to those that were reported in the literature for estimating olive ET when using the traditional approach of SW model [24,26,27]. This approach, however, does not include satellite images to account for the spatial variability of ET. For the main phenological stages, this study suggested that the SW model that uses satellite images and meteorological data could be suitable to estimate olive ET on a pixel-by-pixel basis. For practical application, however, the SW model requires good parameterizations of r st and LAI , which are included in the formulation of r s c . According to Ortega-Farías and Lopez-Olivari [24], the parameterization of r s c depends on the canopy characteristics, which are mainly expressed by the training system and plant water status. In this study, the canopy architecture was maintained almost constant during the study with values of LAI and fc ranging between 1.16–1.38 m2 m−2 and 0.26–0.3, respectively. Also, the drip-irrigated olive orchard was maintained under non-water stress conditions (Ψmd > −1.49 MPa), indicating that the constant value of stomatal resistance (235 ± 60.84 s m−1) used in the SW model was adequate. However, commercial olive orchards that possess different training systems, tree water status, and canopy sizes require the incorporation of spatially distributed values of LAI and stomatal conductance to increase the accuracy of the SW model to simulate ET. To accomplish this task, it is necessary to evaluate sub-models for predicting the spatial variability of Ai, LAI, and r s c of olive orchards with different canopy geometry and tree water status. Finally, the results of this study suggested that the SW model could be used to estimate olive ET as a complement to the satellite-based remote sensing (SBRS) models, like METRIC.

5. Conclusions

Results of this study suggested that the SW model could be a potential tool to compute spatial the variability of olive ET when using satellite images and ground-based climate data. In this case, Rn i and G i were obtained from satellite information and were introduced to the SW model through A i . Simulated values of LEi, Rn i , G i , and A i were generally in good agreement with ground-based measurements within the olive orchard. In this case, the errors and RMSE ranged between 4% and 12% of observed values and 25–39 W m−2, respectively. For daily water consumption, the results indicated that the SW model was able to predict olive ET with errors of 10% and RMSE = 0.3 mm d−1. The observed results were consistent with the literature and are within the acceptable range for applications of remote sensing models in agriculture. Furthermore, the results also encourage the continued testing of this methodology in heterogeneous (vineyards and orchards) canopies with the aim of evaluating the effect of intra-orchard spatial variabilities of soil and canopy vigor on the estimation of ET. However, the practical application of the SW model requires adequate parameterizations of r st and LAI of olive orchards under different training and system.

Acknowledgments

This study was supported by the Chilean government through the projects FONDECYT (No. 1130729) and FONDEF (No. D10I1157) and by the Universidad de Talca through the research program “Adaptation of Agriculture to Climate Change (A2C2)”. The authors would like to thank Manuel Barrera and Alvaro Ried from the “Olivares de Quepu” Company for their technical support.

Author Contributions

Fernando Fuentes-Peñailillo and Samuel Ortega-Farias conceived and designed the experiments, supported the statistical analysis, and wrote the original manuscript draft. Cesar Acevedo-Opazo and David Fonseca-Luengo contributed extensively to results interpretation and discussions and editing of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

List of Symbols

A i Instantaneous available energy leaving the complete canopy (W m−2)
A Si Available energy at the soil surface (W m−2)
A W Percentage of wetted area (%)
bRatio of estimated to observed values (dimensionless)
BBowen ratio (dimensionless)
C C Canopy resistance coefficient (dimensionless)
C S Soil surface resistance coefficient (dimensionless)
CExtinction coefficient of crop for net radiation (dimensionless)
C P Specific heat of the air at constant pressure (1013 J kg−1 °K−1)
CVCoefficient of variation (%)
dZero plane displacement of crop with complete canopy cover (LAI = 4) (m)
D i Water vapor pressure deficit at the reference height (kPa)
EiLatent heat flux from soil evaporation (W m2)
EFReference evapotranspiration fraction (dimensionless)
ETDaily actual evapotranspiration (mm d−1)
ETECDaily actual evapotranspiration obtained from de eddy covariance (EC) method (mm d−1)
ET i Instantaneous actual evapotranspiration computed from the Shuttleworth and Wallace (SW) model (mm h−1)
ET o Reference evapotranspiration (mm d−1)
ET oi Instantaneous reference evapotranspiration (mm h−1)
ET sw Daily actual evapotranspiration estimated using the SW model (mm d−1)
fcFractional crop cover (dimensionless)
G Soil heat flux (W m−2)
G i Instantaneous soil heat flux (W m−2)
HSensible heat flux (W m−2)
H B H corrected using the Bowen ratio approach (W m−2)
I a Index of agreement (dimensionless)
KcCrop coefficient (dimensionless)
Kcbbasal crop coefficient (dimensionless)
Kesoil evaporation coefficient (dimensionless)
LAI Leaf area index (m2 m−2)
LE EC Latent heat flux obtained from the EC method (W m−2)
LE B Instantaneous LE corrected using the Bowen ratio approach (W m−2)
LE i Instantaneous LE estimated using SW (W m−2)
MAEMean absolute error (dimensionless)
N Total number of observations,
n Eddy diffusivity decay constant in crop with complete canopy crop cover (dimensionless)
O i Observed values (W m−2 or mm d−1)
O ¯ Mean of the observed values (W m−2 or mm d−1)
P i Estimated values (W m−2 or mm d−1)
Ppprecipitation (mm)
RHRelative air humidity (%)
Rn Net radiation (W m−2)
Rn avg Average daily values of net radiation (W m2)
Rn i Instantaneous net radiation (W m2)
Rn Si Instantaneous net radiation on the surface of the ground (W m−2)
R s Incident short-wave radiation (W m−2)
R L Incoming long-wave radiation (W m−2)
R L Outgoing long-wave radiation (W m−2)
R2Coefficient of determination (dimensionless)
r b Mean boundary layer resistance per unit area of vegetation (s m−1)
r a c Bulk boundary layer resistance of the vegetative elements in the canopy (s m−1)
r a a Aerodynamic resistance between the canopy source height and reference level (s m−1)
r s c Canopy resistance (s m−1)
r a s Aerodynamic resistance between the soil and canopy source height (s m−1)
r s s Soil surface resistance (s m−1)
r st Mean stomatal resistance (s m−1)
TaAir temperature (°C)
TavgAverage daily values of temperature (°C)
T i Latent heat flux from the olive transpiration (W m−2)
T s Superficial temperature (°K)
T n Near surface air temperature (°K)
uWind speed (m s−1)
uavgAverage daily values of wind speed (m s−1)
ZomSurface roughness (m)
zo Roughness length of crop with complete canopy cover (s m−1)
z ó Roughness length of soil surface (s m−1)
α Surface albedo (dimensionless)
α i Instantaneous surface albedo (dimensionless)
ε a Effective atmospheric emissivity (dimensionless)
ε 0 Superficial emissivity (dimensionless)
λ Latent heat of vaporization (1013 MJ kg−1),
ρ w Water density (1000 kg m−3),
σ Stefan-Boltzmann constant (5.67 × 10−8 Wm−2 °K−4)
γ Psychrometric constant (kPa °K−1)
τ sw Wide-band atmospheric transmissivity (dimensionless)
ΨmdMidday stem water potential (MPa)
Δ Slope of the saturation vapor pressure curve at the mean temperature (kPa °C−1)

References

  1. Ortega-Farías, S.; Irmak, S.; Cuenca, R.H. Special issue on evapotranspiration measurement and modeling. Irrig. Sci. 2009, 28, 1–3. [Google Scholar] [CrossRef]
  2. Fereres, E.; Soriano, M.A. Deficit irrigation for reducing agricultural water use. J. Exp. Bot. 2007, 58, 147–159. [Google Scholar] [CrossRef] [PubMed]
  3. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements—FAO Irrigation and Drainage Paper 56; Food and Agriculture Organization of the United Nations (FAO): Rome, Italy, 1998. [Google Scholar]
  4. Er-Raki, S.; Chehbouni, A.; Boulet, G.; Williams, D.G. Using the dual approach of FAO-56 for partitioning ET into soil and plant components for olive orchards in a semi-arid region. Agric. Water Manag. 2010, 97, 1769–1778. [Google Scholar] [CrossRef] [Green Version]
  5. Er-Raki, S.; Chehbouni, A.; Duchemin, B. Combining Satellite Remote Sensing Data with the FAO-56 Dual Approach for Water Use Mapping in Irrigated Wheat Fields of a Semi-Arid Region. Remote Sens. 2010, 2, 375–387. [Google Scholar] [CrossRef]
  6. Cammalleri, C.; Ciraolo, G.; Minacapilli, M.; Rallo, G. Evapotranspiration from an Olive Orchard using Remote Sensing-Based Dual Crop Coefficient Approach. Water Resour. Manag. 2013, 27, 4877–4895. [Google Scholar] [CrossRef]
  7. Paço, T.A.; Pocas, I.; Cunha, M.; Silvestre, J.C.; Santos, F.L.; Paredes, P.; Pereira, L.S. Evapotranspiration and crop coefficients for a super intensive olive orchard. An application of SIMDualKc and METRIC models using ground and satellite observations. J. Hydrol. 2014, 519, 2067–2080. [Google Scholar] [CrossRef]
  8. Martínez-Cob, A.; Faci, J.M. Evapotranspiration of an hedge-pruned olive orchard in a semiarid area of NE Spain. Agric. Water Manag. 2010, 97, 410–418. [Google Scholar] [CrossRef] [Green Version]
  9. Lopez-Olivari, R.; Ortega-Farías, S.; Poblete-Echeverría, C. Partitioning of net radiation and evapotranspiration over a superintensive drip-irrigated olive orchard. Irrig. Sci. 2016, 34, 17–31. [Google Scholar] [CrossRef]
  10. Allen, R.G.; Pereira, L.S. Estimating crop coefficients from fraction of ground cover and height. Irrig. Sci. 2009, 28, 17–34. [Google Scholar] [CrossRef]
  11. Pôças, I.; Paco, T.A.; Paredes, P.; Cunha, M.; Pereira, L.S. Estimation of Actual Crop Coefficients Using Remotely Sensed Vegetation Indices and Soil Water Balance Modelled Data. Remote Sens. 2015, 7, 2373–2400. [Google Scholar] [CrossRef]
  12. Er-Raki, S.; Chehbouni, A.; Hoedjes, J.; Ezzahar, J.; Duchemin, B.; Jacob, F. Improvement of FAO-56 method for olive orchards through sequential assimilation of thermal infrared-based estimates of ET. Agric. Water Manag. 2008, 95, 309–321. [Google Scholar] [CrossRef]
  13. Anadranistakis, M.; Liakatas, A.; Kerkides, P.; Rizos, S.; Gavanosis, J.; Poulovassilis, A. Crop water requirements model tested for crops grown in Greece. Agric. Water Manag. 2000, 45, 297–316. [Google Scholar] [CrossRef]
  14. Ortega-Farías, S.; Poblete-Echeverría, C.; Brisson, N. Parameterization of a two-layer model for estimating vineyard evapotranspiration using meteorological measurements. Agric. For. Meteorol. 2010, 150, 276–286. [Google Scholar] [CrossRef]
  15. Testi, L.; Villalobos, F.J.; Orgaz, F. Evapotranspiration of a young irrigated olive orchard in southern Spain. Agric. For. Meteorol. 2004, 121, 1–18. [Google Scholar] [CrossRef]
  16. Were, A.; Villagarcía, L.; Domingo, F.; Moro, M.J.; Dolman, A.J. Aggregating spatial heterogeneity in a bush vegetation patch in semi-arid SE Spain: A multi-layer model versus a single-layer model. J. Hydrol. 2008, 349, 156–167. [Google Scholar] [CrossRef]
  17. Zhou, M.C.; Ishidaira, H.; Hapuarachchi, H.P.; Magome, J.; Kiem, A.S.; Takeuchi, K. Estimating potential evapotranspiration using Shuttleworth–Wallace model and NOAA-AVHRR NDVI data to feed a distributed hydrological model over the Mekong River basin. J. Hydrol. 2006, 327, 151–173. [Google Scholar] [CrossRef]
  18. Kato, T.; Kamichika, M. Determination of a crop coefficient for evapotranspiration in a sparse sorghum field. Irrig. Drain. 2006, 55, 165–175. [Google Scholar] [CrossRef]
  19. Gardiol, J.M.; Serio, L.A.; Della Maggiora, A.I. Modelling evapotranspiration of corn (Zea mays) under different plant densities. J. Hydrol. 2003, 271, 188–196. [Google Scholar] [CrossRef]
  20. Zhang, B.; Kang, S.; Li, F.; Zhang, L. Comparison of three evapotranspiration models to Bowen ratio-energy balance method for a vineyard in an arid desert region of northwest China. Agric. For. Meteorol. 2008, 148, 1629–1640. [Google Scholar] [CrossRef]
  21. Ortega-Farías, S.; Carrasco, M.; Olioso, A.; Acevedo, C.; Poblete, C. Latent heat flux over Cabernet Sauvignon vineyard using the Shuttleworth and Wallace model. Irrig. Sci. 2006, 25, 161–170. [Google Scholar] [CrossRef]
  22. Zhang, B.; Kang, S.; Zhang, L.; Tong, L.; Du, T.; Li, F.; Zhang, J. An evapotranspiration model for sparsely vegetated canopies under partial root-zone irrigation. Agric. For. Meteorol. 2009, 149, 2007–2011. [Google Scholar] [CrossRef]
  23. Zhao, P.; Li, S.; Li, F.; Du, T.; Tong, L.; Kang, S. Comparison of dual crop coefficient method and Shuttleworth–Wallace model in evapotranspiration partitioning in a vineyard of northwest China. Agric. Water Manag. 2015, 160, 41–56. [Google Scholar] [CrossRef]
  24. Ortega-Farías, S.; Lopez-Olivari, R. Validation of a Two-Layer Model to Estimate Latent Heat Flux and Evapotranspiration in a Drip-Irrigated Olive Orchard. Trans. ASABE 2012, 55, 1169–1178. [Google Scholar] [CrossRef]
  25. Brenner, A.J.; Incoll, L.D. The effect of clumping and stomatal response on evaporation from sparsely vegetated shrublands. Agric. For. Meteorol. 1997, 84, 187–205. [Google Scholar] [CrossRef]
  26. Ortega-Farías, S.; Ortega-Salazar, S.; Poblete, T.; Kilic, A.; Allen, R.; Poblete-Echeverría, C.; Ahumada-Orellana, L.; Zuniga, M.; Sepulveda, D. Estimation of energy balance components over a drip-irrigated olive orchard using thermal and multispectral cameras placed on a helicopter-based unmanned aerial vehicle (UAV). Remote Sens. 2016, 8, 638. [Google Scholar] [CrossRef]
  27. Ortega-Farías, S.; Ortega-Salazar, S.; Aguilar, R.; De la Fuente, D.; Fuentes, F. Evaluation of a model to estimate net radiation over a drip-irrigated olive orchard using landsat satellite images. Acta Hortic. 2014, 1057, 309–314. [Google Scholar] [CrossRef]
  28. Carrasco-Benavides, M.; Ortega-Farías, S.; Lagos, L.O.; Kleissl, J.; Morales-Salinas, L.; Kilic, A. Parameterization of the satellite-based model (METRIC) for the estimation of instantaneous surface energy balance components over a drip-irrigated vineyard. Remote Sens. 2011, 6, 11342–11371. [Google Scholar] [CrossRef]
  29. Cuenca, R.H.; Ciotti, S.P.; Hagimoto, Y. Application of Landsat to Evaluate Effects of Irrigation Forbearance. Remote Sens. 2011, 5, 3776–3802. [Google Scholar] [CrossRef]
  30. Long, D.; Singh, V.P. A modified surface energy balance algorithm for land (M-SEBAL) based on a trapezoidal framework. Water Resour. Res. 2012, 48. [Google Scholar] [CrossRef]
  31. Long, D.; Singh, V.P. Assessing the impact of end-member selection on the accuracy of satellite-based spatial variability models for actual evapotranspiration estimation. Water Resour. Res. 2013, 49, 2601–2618. [Google Scholar] [CrossRef]
  32. Dhungel, R.; Allen, R.G.; Trezza, R.; Robison, C.W. Comparison of Latent Heat Flux Using Aerodynamic Methods and Using the Penman-Monteith Method with Satellite-Based Surface Energy Balance. Remote Sens. 2014, 6, 8844–8877. [Google Scholar] [CrossRef]
  33. Tognetti, R.; D’Andria, R.; Sacchi, R.; Lavini, A.; Morelli, G.; Alvino, A. Deficit irrigation affects seasonal changes in leaf physiology and oil quality of Olea europaea (cultivars Frantoio and Leccino). Ann. Appl. Biol. 2007, 150, 169–186. [Google Scholar] [CrossRef]
  34. Pérez-López, D.; Gijón, M.C.; Marino, J.; Moriana, A. Water relation response to soil chilling of six olive (Olea europaea L.) cultivars with different frost resistance. Span. J. Agric. Res. 2010, 8, 780–789. [Google Scholar] [CrossRef]
  35. Dichio, B.; Xiloyannis, C.; Sofo, A.; Montanaro, G. Osmotic regulation in leaves and roots of olive trees during a water deficit and rewatering. Tree Physiol. 2006, 26, 179–185. [Google Scholar] [CrossRef] [PubMed]
  36. Poblete-Echeverría, C.; Sepúlveda-Reyes, D.; Ortega-Farías, S. Effect of height and time lag on the estimation of sensible heat flux over a drip-irrigated vineyard using the surface renewal (SR) method across distinct phenological stages. Agric. Water Manag. 2014, 141, 74–83. [Google Scholar] [CrossRef]
  37. Shao, C.; Chen, J.; Li, L.; Xu, W.; Chen, S.; Gwen, T.; Xu, J.; Zhang, W. Spatial variability in soil heat flux at three Inner Mongolia steppe ecosystems. Agric. For. Meteorol. 2008, 148, 1433–1443. [Google Scholar] [CrossRef]
  38. Schotanus, P.; Nieuwstadt, F.T.M.; de Bruin, H.A.R. Temperature measurement with a sonic anemometer and its application to heat and moisture fluxes. Bound.-Layer Meteorol. 1983, 26, 81–93. [Google Scholar] [CrossRef]
  39. Webb, E.K.; Pearman, G.I.; Leuning, R. Correction of flux measurements for density effects due to heat and water vapour transfer. Q. J. R. Meteorol. Soc. 1980, 106, 85–100. [Google Scholar] [CrossRef]
  40. Twine, T.E.; Kustas, W.P.; Norman, J.M.; Cook, D.R. Correcting eddy-covariance flux underestimates over a grassland. Agric. For. Meteorol. 2000, 103, 279–300. [Google Scholar] [CrossRef]
  41. Shuttleworth, W.J.; Wallace, J.S. Evaporation from sparse crops—An energy combination theory. Q. J. R. Meteorol. Soc. 1985, 111, 839–855. [Google Scholar] [CrossRef]
  42. Allen, R.G.; Tasumi, M.; Trezza, R. Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC)—Model. J. Irrig. Drain. Eng. 2007, 133, 380–394. [Google Scholar] [CrossRef]
  43. Storey, J.; Scaramuzza, P.; Schmidt, G.; Barsi, J. Landsat 7 scan line corrector-off gap-filled product development. In Proceedings of the Pecora 16 Global Priorities in Land Remote Sensing, Sioux Falls, SD, USA, 23–27 October 2005. [Google Scholar]
  44. Zeng, C.; Shen, H.; Zhang, L. Recovering missing pixels for Landsat ETM plus SLC-off imagery using multi-temporal regression analysis and a regularization method. Remote Sens. Environ. 2013, 131, 182–194. [Google Scholar] [CrossRef]
  45. Spiliotopoulos, M.; Holden, N.M.; Loukas, A. Mapping evapotranspiration coefficients in a temperate maritime climate using the metric model and landsat TM. Water 2017, 9, 23. [Google Scholar] [CrossRef]
  46. Reyes-González, A.; Kjaersgaard, J.; Trooien, T.; Hay, C.; Ahiablame, L. Comparative Analysis of METRIC Model and Atmometer Methods for Estimating Actual Evapotranspiration. Int. J. Agron. 2017, 2017, 3632501. [Google Scholar] [CrossRef]
  47. Jimenez-Bello, M.A.; Castel, J.R.; Testi, L.; Intrigliolo, D.S. Assessment of a Remote Sensing Energy Balance Methodology (SEBAL) Using Different Interpolation Methods to Determine Evapotranspiration in a Citrus Orchard. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 1465–1477. [Google Scholar] [CrossRef]
  48. Song, L.; Liu, S.; Kustas, W.P.; Zhou, J.; Xu, Z.; Xia, T.; Li, M. Application of remote sensing-based two-source energy balance model for mapping field surface fluxes with composite and component surface temperatures. Agric. For. Meteorol. 2016, 230, 8–19. [Google Scholar] [CrossRef]
  49. Serra, P.; Salvati, L.; Queralt, E.; Pin, C.; Gonzalez, O.; Pons, X. Estimating Water Consumption and Irrigation Requirements in a Long-Established Mediterranean Rural Community by Remote Sensing and Field Data. Irrig. Drain. 2016, 65, 578–588. [Google Scholar] [CrossRef]
  50. Mayer, D.G.; Butler, D.G. Statistical validation. Ecol. Model. 1993, 68, 21–32. [Google Scholar] [CrossRef]
  51. Willmott, C.J.; Ackleson, S.G.; Davis, R.E.; Feddema, J.J.; Klink, K.M.; Legates, D.R.; O’Donnell, J.; Rowe, C.M. Statistics for the Evaluation and Comparison of Models. J. Geophys. Res. Oceans 1985, 90, 8995–9005. [Google Scholar] [CrossRef]
  52. Williams, D.G.; Cable, W.; Hultine, K.; Hoedjes, J.; Yepez, E.A.; Simonneaux, V.; Er-Raki, S.; Boulet, G.; de Bruin, H.; Chehbouni, A.; et al. Evapotranspiration components determined by stable isotope, sap flow and eddy covariance techniques. Agric. For. Meteorol. 2004, 125, 241–258. [Google Scholar] [CrossRef]
  53. Ezzahar, J.; Chehbouni, A.; Hoedjes, J.C.B.; Er-Raki, S.; Chehbouni, A.; Boulet, G.; Bonnefond, J.M.; de Bruin, H.A.R. The use of the scintillation technique for monitoring seasonal water consumption of olive orchards in a semi-arid region. Agric. Water Manag. 2007, 89, 173–184. [Google Scholar] [CrossRef] [Green Version]
  54. Villalobos, F.J.; Testi, L.; Orgaz, F.; García-Tejera, O.; López-Bernal, Á.; Victoria Gonzalez-Dugo, M.; Ballester-Lurbe, C.; Ramon Castel, J.; Jose Alarcon-Cabanero, J.; Nicolas-Nicolas, E.; et al. Modelling canopy conductance and transpiration of fruit trees in Mediterranean areas: A simplified approach. Agric. For. Meteorol. 2013, 171, 93–103. [Google Scholar] [CrossRef]
  55. Testi, L.; Villalobos, F.J.; Orgaz, F.; Fereres, E. Water requirements of olive orchards: I simulation of daily evapotranspiration for scenario analysis. Irrig. Sci. 2006, 24, 69–76. [Google Scholar] [CrossRef]
  56. Lee, X.; Black, T.A. Atmospheric turbulence within and above a douglas-fir stand. Part II: Eddy fluxes of sensible heat and water vapour. Bound.-Layer Meteorol. 1993, 64, 369–389. [Google Scholar] [CrossRef]
  57. Wilson, K.; Goldstein, A.; Falge, E.; Aubinet, M.; Baldocchi, D.; Berbigier, P.; Bernhofer, C.; Ceulemans, R.; Dolman, H.; Field, C.; et al. Energy balance closure at FLUXNET sites. Agric. For. Meteorol. 2002, 113, 223–243. [Google Scholar] [CrossRef]
  58. Leuning, R.; van Gorsel, E.; Massman, W.J.; Isaac, P.R. Reflections on the surface energy imbalance problem. Agric. For. Meteorol. 2012, 156, 65–74. [Google Scholar] [CrossRef]
  59. Allen, R.G.; Pereira, L.S.; Howell, T.A.; Jensen, M.E. Evapotranspiration information reporting: I. Factors governing measurement accuracy. Agric. Water Manag. 2011, 98, 899–920. [Google Scholar] [CrossRef]
  60. Liu, H.; Foken, T. A modified Bowen ratio method to determine sensible and latent heat fluxes. Meteorol. Z. 2001, 10, 71–80. [Google Scholar] [CrossRef]
  61. De Teixeira, A.H.C.; Bastiaanssen, W.G.M.; Ahmad, M.D.; Bos, M.G. Reviewing SEBAL input parameters for assessing evapotranspiration and water productivity for the Low-Middle São Francisco River basin, Brazil. Agric. For. Meteorol. 2009, 149, 462–476. [Google Scholar] [CrossRef] [Green Version]
  62. Liu, S.; Lu, L.; Mao, D.; Jia, L. Evaluating parameterizations of aerodynamic resistance to heat transfer using field measurements. Hydrol. Earth Syst. Sci. 2007, 11, 769–783. [Google Scholar] [CrossRef]
Figure 1. (a) Daily variation of net radiation ( Rn avg ) and reference evapotranspiration ( ET o ); (b) average daily values of temperature (Tavg) and average daily values of wind speed (uavg). The vertical bars indicate the time of satellite overpass (SO). Also, phenological stages are indicated by an abbreviation followed by an arrow where F = flowering, FS = fruit Set, PH = pit hardening and FC = fruit colouring.
Figure 1. (a) Daily variation of net radiation ( Rn avg ) and reference evapotranspiration ( ET o ); (b) average daily values of temperature (Tavg) and average daily values of wind speed (uavg). The vertical bars indicate the time of satellite overpass (SO). Also, phenological stages are indicated by an abbreviation followed by an arrow where F = flowering, FS = fruit Set, PH = pit hardening and FC = fruit colouring.
Water 10 00339 g001
Figure 2. Energy balance closure at interval of 30 min for the days when the satellite overpassed the experimental site. H, LE, Rn, and G correspond to sensible heat flux (W m−2), latent heat flux (W m−2), net radiation (W m−2), and soil heat flux (W m−2), respectively. The coefficient of determination (R2) was 0.88 and the ( Rn G )/( H + LE ) ratio was 0.89.
Figure 2. Energy balance closure at interval of 30 min for the days when the satellite overpassed the experimental site. H, LE, Rn, and G correspond to sensible heat flux (W m−2), latent heat flux (W m−2), net radiation (W m−2), and soil heat flux (W m−2), respectively. The coefficient of determination (R2) was 0.88 and the ( Rn G )/( H + LE ) ratio was 0.89.
Water 10 00339 g002
Figure 3. Comparison between estimated and observed values of instantaneous latent heat flux ( LE i ) from the two-source model and instantaneous available energy ( A i ) at the time of satellite overpass. Values of index of agreement (d) for LE i and A i were 0.8 and 0.82, respectively.
Figure 3. Comparison between estimated and observed values of instantaneous latent heat flux ( LE i ) from the two-source model and instantaneous available energy ( A i ) at the time of satellite overpass. Values of index of agreement (d) for LE i and A i were 0.8 and 0.82, respectively.
Water 10 00339 g003
Figure 4. Comparison between estimated and observed values of daily actual evapotranspiration of a drip irrigated olive orchard. ( ET sw ). Index of agreement ( I a ) for ET sw was equal to 0.95. This comparison was done using days when satellite images were available.
Figure 4. Comparison between estimated and observed values of daily actual evapotranspiration of a drip irrigated olive orchard. ( ET sw ). Index of agreement ( I a ) for ET sw was equal to 0.95. This comparison was done using days when satellite images were available.
Water 10 00339 g004
Figure 5. Maps of daily actual evapotranspiration (ET) computed from the SW model using remote sensing and meteorological data.
Figure 5. Maps of daily actual evapotranspiration (ET) computed from the SW model using remote sensing and meteorological data.
Water 10 00339 g005
Figure 6. Maps of instantaneous Available energy (Ai) generated from satellite images and ground-based climate data.
Figure 6. Maps of instantaneous Available energy (Ai) generated from satellite images and ground-based climate data.
Water 10 00339 g006
Table 1. Images selected for processing the Shuttleworth and Wallace model (SW) and main phenological stages of a drip-irrigated olive orchard.
Table 1. Images selected for processing the Shuttleworth and Wallace model (SW) and main phenological stages of a drip-irrigated olive orchard.
DateDay of YearOverpass TimeScene Cloud CoverPhenological Stages
(dd-mm-yy)(DOY)(UTC)(%)
4/02/2009352:24:02 p.m.1FC
20/02/2009512:24:12 p.m.2FC
3/11/20093072:24:42 p.m.7F
5/12/20093392:25:05 p.m.21FS
21/12/20093552:25:21 p.m.1PH
6/01/201062:25:38 p.m.1PH
22/01/2010222:25:52 p.m.0PH
F = flowering; FS = fruit set; PH = pit hardening; FC = fruit colouring.
Table 2. Meteorological conditions during the study period.
Table 2. Meteorological conditions during the study period.
DOYDate E T o E T oi RHTaU
(mm d−1)(mm h−1)(%)(°C)(m s−1)
354/02/20097.880.5651.9422.910.81
5120/02/20096.880.5254.1617.924.11
3073/11/20095.130.5458.8514.122.33
3395/12/20096.470.5756.1216.361.17
35521/12/20096.690.6947.621.292.69
66/01/20107.840.653.5321.332.69
2222/01/20107.750.6152.9922.171.29
Mean6.950.5853.619.442.16
D.E.0.990.063.53.331.15
where ET o = daily reference evapotranspiration, ET oi = instantaneous reference evapotranspiration, RH = Relative air humidity, Ta = Air Temperature and u = wind speed.
Table 3. Model validation of instantaneous values of net radiation ( Rn i ), soil heat flux ( G i ), available energy ( A i ), and latent heat flux ( LE i ) over a drip irrigated olive orchard. Also, daily actual evapotranspiration ( ET SW ) is included in the analysis.
Table 3. Model validation of instantaneous values of net radiation ( Rn i ), soil heat flux ( G i ), available energy ( A i ), and latent heat flux ( LE i ) over a drip irrigated olive orchard. Also, daily actual evapotranspiration ( ET SW ) is included in the analysis.
VariableRMSEMAE I a bt-Test
Rn i (W m−2)39320.790.96F
G i (W m−2)33270.50.92F
A i (W m−2)25210.820.96F
LE i (W m−2)26200.80.88F
ET SW (mm d−1)0.310.280.950.90F
where RMSE = Root mean square error, MAE = Mean absolute error, b = Ratio of estimated to observed values and I a = Index of agreement. T = null hypothesis (b = 1) True; F = alternative hypothesis (b ≠ 1).
Table 4. Sensitivity analysis to instantaneous actual evapotranspiration ( ET i ).
Table 4. Sensitivity analysis to instantaneous actual evapotranspiration ( ET i ).
VariablesSymbolParameters30%−30%
Extinction coefficient of crop for net radiation (dimensionless) C 0.66 *0.090.53
Zero plane displacement of crop with complete canopy cover (LAI = 4) (m) d 2.01 *−0.263.29
Leaf Area Index (m2 m−2) LAI 1.29 *13.93−14.52
Eddy diffusivity decay constant in crop with complete canopy crop cover (dimensionless) n 2.5 *2.16−0.95
Mean boundary layer resistance per unit area of vegetation (s m−1) r b 25 *1.87−1.44
Mean stomatal resistance (s m−1) r st 235 *−13.0619.96
Aerodynamic resistance between canopy source height and reference level (s m−1) r a a **4.51−4.45
Bulk boundary layer resistance of the vegetative elements in the canopy (s m−1) r a c **1.87−1.44
Aerodynamic resistance between the substrate and canopy source height (s m−1) r a s **1.99−1.61
Bulk stomatal resistance of the canopy (s m−1) r s c **−12.8220.3
Surface resistance of the substrate (s m−1) r s s 2000 *−3.195.93
Instantaneous net radiation (W m−2) Rn i **20.60−20.03
Instantaneous Soil heat flux (W m−2) G i **−2.182.75
Instantaneous available energy (W m−2) A i **17.96−17.39
* Constant values; ** Computed values.

Share and Cite

MDPI and ACS Style

Fuentes-Peñailillo, F.; Ortega-Farías, S.; Acevedo-Opazo, C.; Fonseca-Luengo, D. Implementation of a Two-Source Model for Estimating the Spatial Variability of Olive Evapotranspiration Using Satellite Images and Ground-Based Climate Data. Water 2018, 10, 339. https://doi.org/10.3390/w10030339

AMA Style

Fuentes-Peñailillo F, Ortega-Farías S, Acevedo-Opazo C, Fonseca-Luengo D. Implementation of a Two-Source Model for Estimating the Spatial Variability of Olive Evapotranspiration Using Satellite Images and Ground-Based Climate Data. Water. 2018; 10(3):339. https://doi.org/10.3390/w10030339

Chicago/Turabian Style

Fuentes-Peñailillo, Fernando, Samuel Ortega-Farías, César Acevedo-Opazo, and David Fonseca-Luengo. 2018. "Implementation of a Two-Source Model for Estimating the Spatial Variability of Olive Evapotranspiration Using Satellite Images and Ground-Based Climate Data" Water 10, no. 3: 339. https://doi.org/10.3390/w10030339

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop