Abstract

Accurate estimation of crop evapotranspiration (ET) is a key factor in agricultural water management including irrigated agriculture. The objective of this study was to compare ET estimated from the satellite-based remote sensing METRIC model to in situ atmometer readings. Atmometer readings were recorded from three sites in eastern South Dakota every morning between 8:15 and 8:30 AM for the duration of the 2016 growing season. Seven corresponding clear sky images from Landsat 7 and Landsat 8 (Path 29, Row 29) were processed and used for comparison. Three corn fields in three sites were used to compare actual evapotranspiration (). The results showed a good relationship between estimated by the METRIC model (-METRIC) and estimated with atmometer (-atm) ( = 0.87, index of agreement of 0.84, and RMSE = 0.65 mm day−1). However, -atm values were consistently lower than -METRIC values. The differences in daily between the two methods increase with high wind speed values (>4 m s−1). Results from this study are useful for improving irrigation water management at local and field scales.

1. Introduction

With increasing demands placed on freshwater resources worldwide, it is necessary to accurately estimate crop water consumption efficiently. Crop water use information is needed for a range of applications, including improving agricultural water management, irrigation and crop selection, water resource planning, water rights management, and water regulations [13]. Irrigated agriculture produces 40% of global food and fiber supply from 20% of the world’s croplands [4]. In arid areas, up to 90% of all water withdrawals may be for irrigation purposes [5]. With increasing population and water uses a scarce water supply is put under additional pressure and other water users relying on the same water supply may experience insufficient water allocations. At the same time, a reduction of irrigation water supply may result in loss of production and, ultimately, threatened food security. There is an opportunity, however, to optimize the management of water in agricultural production systems, and the accurate estimation of evapotranspiration (ET) is critical in that regard.

ET is the loss of water from the land surface to the atmosphere through two processes, namely, evaporation (E) from soil and water surfaces and transpiration (T) from vegetative surfaces [6, 7]. ET rates are affected by weather conditions such as solar radiation, air temperature, wind speed and air vapor pressure deficit, and plant and soil characteristics conditions [8, 9]. Different methods, direct and indirect, exist to estimate ET. Direct methods include weighing lysimeters and soil water balance estimations, while indirect methods include pan evaporation, atmometer, Bowen Ratio Energy Balance System (BREBS), Eddy Covariance (EC), scintillometer, sap flow, and remote sensing [1]. An attractive property of satellite-based remote sensing ET estimates using Landsat imagery is its coverage on a field by field basis at a regional scale [1, 10, 11].

Several models have been developed to estimate ET using remote sensing. One of them is the Mapping EvapoTranspiration at High Resolution using Internalized Calibration (METRIC) model. METRIC utilizes the innovative Surface Energy Balance Algorithm for Land (SEBAL) method for estimating sensible heat flux. METRIC uses the near surface to air temperature gradient () for each pixel within an image based on a regression relationship between the and radiometric surface temperature () of two anchor pixels. The anchor pixels represent the conditions of an agricultural field with full vegetation cover and maximum crop ET (cold condition) and a bare agricultural field with no vegetation cover (hot condition) [10, 12].

One of the advantages of the METRIC model compared to previous surface energy balance-based models for use in arid areas is that it utilizes reference evapotranspiration () for estimating actual evapotranspiration () at the cold pixel condition [10]. Because is based on ground-based meteorological measurements and because is calibrated to account for atmospheric conditions common in arid and semiarid conditions, such as horizontal advection, the METRIC model is particularly useful for estimations under arid or semiarid conditions [10].

Previous studies have compared METRIC model outputs to other methods of estimation such as weighing lysimeter [10], soil water balance [13], Bowen Ratio Energy Balance System (BREBS) [1417], Eddy Covariance (EC), for example, [1821], Large Aperture Scintillometer (LAS) [22], and the METRIC-MODIS method [23]. These studies showed from moderate to strong relationships between observed and METRIC-estimated , indicating that the METRIC model is a useful tool for estimating accurate at local and field scales. In addition, the METRIC model has been compared with other models such as water balance model [24], trapezoid interpolation model (TIM) [25, 26], two-source energy balance model (TSEB) [27, 28], SIMDualKc model [29], and the Landsat-MODIS fusion model [30]. However, comparison of METRIC model outputs to an atmometer for estimation has not yet been attempted.

An atmometer is a device that measures the amount of water evaporated from wet porous surface to the atmosphere [31]. Atmometers are simple and inexpensive devices, which consist of a ceramic evaporation plate (Bellani plate) covered by a green canvas, mounted on top of a cylindrical water reservoir, to provide a visual interpretation of atmospheric demand for pulling water out of the vegetation and the soil [32, 33]. The standard model with number 54 green canvas is recommended for measuring alfalfa similar to the alfalfa-based Penman-Monteith , while number 30 green canvas is designed to simulate grass similar to the grass-based Penman-Monteith [32, 34].

Research demonstrated that estimated with atmometers was moderately correlated ( = ±0.70) with weighing lysimeters values [35, 36], strongly correlated ( = 0.90) with pan evaporation values [37, 38], and strongly correlated ( = 0.92) with agrometeorological data values, for example, [3943].

While METRIC is a well-recognized and widely used method to estimate ET, as noted by you, it is a complex and time intensive method that requires an experienced user to be run. For practical applications of crop water use estimations such as for irrigation scheduling, other in situ methods may be more appropriate. This is where atmometers may serve a purpose since they are simple and intuitive to use. They may be less accurate compared to other methods, but for applications where slightly higher uncertainty is acceptable, atmometers may be a good option. The objective of this study was to compare estimated from satellite-based remote sensing METRIC model to estimated with atmometers in corn fields in eastern South Dakota.

2. Material and Methods

2.1. Study Area

The study was carried out at three sites in eastern South Dakota at Brookings (44°19′N, 96°46′W), Volga (44°18′N, 96°55′W), and Oak Lake (44°30′N, 96°31′W) at elevations 500, 497, and 574 m above sea level, respectively. Three corn fields near to each atmometer (within 5 km) were selected and considered to estimate (Figure 1). The population density was approximately 78,000 plants ha−1 in all fields. Corn fields in Brookings and Volga had 0–2% slope, while the Oak Lake had 2–6% slope [44]. All fields used in this study are in corn-soybean crop rotation system. All corn fields emerged in early May and harvested in late October. The average annual precipitation is 533 mm, of which 3/4 typically falls during the growing season April through October. The mean annual maximum temperature is 12.3°C, minimum 0.3°C, and mean 6.3°C. The climate of the study area is classified as moist subhumid according to the Thornthwaite climate classification system [45].

2.2. Landsat Images

We used seven clear sky images from Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) (Path 29, Row 29), Table 1. The images were downloaded from the United States Geological Survey (USGS) EROS Datacenter (http://glovis.usgs.gov) [46]. The images were selected based on the temporal coverage and cloud-free conditions. Images with cloud located > 10 km from all study sites were considered acceptable. The images were processed using the METRIC model running in ERDAS Imagine Software environment [47]. The time of satellite overpass of both Landsats ranged from 11:11 to 11:14 AM, local time (Table 1). Landsat 7 and Landsat 8 have a pixel resolution of 30 m by 30 m in the shortwave bands and 60 m by 60 m and 100 m by 100 m in the thermal band, respectively.

The wedge-shaped gaps appearing within the Landsat 7 images as a result of the SLC-off issue were removed using the Imagine built-in focal analysis tool. During the process, the gap is filled iteratively based on information from nearby pixels. The gap filling is completed prior to image processing (http://landsat.usgs.gov/gap-filling-landsat-7-slc-single-scenes-using-erdas-imagine-TM) [48]. An example of the process is shown in Figure 2.

2.3. METRIC Model

estimations use the METRIC model approach as described by Allen et al. [1, 6].

METRIC model is a remote sensing image processing model that computes instantaneous ET values as a residual of the surface energy balance equation [1, 6, 49, 50]:where LE is the latent heat flux (W m−2), or ET (mm day−1), is the net radiation (W m−2), is the soil heat flux (W m−2), and is the sensible heat flux (W m−2).

Net radiation () is calculated using surface reflectance and surface temperature () derived by satellite imagery. is the difference between incoming shortwave radiation and outgoing longwave radiation computed aswhere is the incoming shortwave radiation (W m−2) (solar radiation), α is surface albedo (dimensionless), is the incoming longwave radiation (W m−2), is the outgoing longwave radiation (W m−2), and is the surface thermal emissivity (dimensionless). is the fraction of incoming longwave radiation reflected from the surface.

Soil heat flux () is the magnitude of the heat flux stored or released into the soil. was computed using the following equations described by Tasumi [51]:Sensible heat flux () was determined using the aerodynamic based heat transfer equation as follows:where is the air density (kg m−3), is the air specific heat (1004 J kg−1 K−1), is the temperature difference between two heights (0.1 m) and (2 m), and is the aerodynamic resistance to heat transfer (s m−1).

For the estimations, the METRIC model uses the CIMEC (Calibration using Inverse Modeling of Extreme Conditions) procedure described by Bastiaanssen et al. [52] and Allen et al. [6] to calibrate the near surface to air temperature difference for each pixel within an image based on a regression relationship between and of two anchor pixels (hot and cold). The advantage of the CIMEC approach within the METRIC model reduces possible impacts of biases in estimation of aerodynamic stability correction and surface roughness [6].

In this study, hot and cold pixels were selected for each image in agricultural fields near to the weather stations (<15 km). The hot pixel was selected in a bare agricultural field with no vegetation cover, while the cold pixel was selected in an agricultural field with full vegetation cover. An adjustment factor is used to adjust the ET estimation to reflect the actual evaporation amount at the cold and hot pixels. A value of 1.05 is typically used for the cold pixel but may be adjusted down in situations where the no field with vegetation cover is present within the image. For the hot pixel, a value of 0.05 is typically used for situations where the top 10–15 cm of the soil is dry and no rainfall has occurred for several weeks. The bare soil evaporation model by Allen et al. [8] was used to establish the appropriate adjustment factor higher than 0.05 to account for residual soil evaporation from antecedent rainfall.

Based on LE values, the instantaneous values of ET were computed for each pixel aswhere is the hourly instantaneous ET (mm h−1), 3600 is used to convert to hours, LE is the latent heat flux (W m−2) consumed by ET, is the density of water (1000 kg m−3), and is the latent heat of evaporation (j kg−1), which is computed asThe reference ET fraction (F) or crop coefficient () was calculated based on for each pixel and was obtained from local weather data.Daily values of ET (ET24) (mm day−1) for each pixel were calculated as follows:where is the reference ET fraction, is the cumulative alfalfa reference for the day (mm day−1), and ET24 is the actual evapotranspiration for the entire 24-hour period (mm day−1).

Monthly and seasonal are calculated by interpolation daily values of between images and multiplying by for each day and then integrated over the specific month [6]. The interpolation values of are made using a linear interpolation or a curvilinear interpolation function such as a spline function [53]. According to Allen et al. [6] one cloud-free satellite image per month is enough to develop curves for seasonal estimations.

2.4. Meteorological Data

Hourly weather data was used for the internal calibration of the METRIC model. The weather observations were taken from the automatic Brookings, Volga, and Oak Lake stations. All weather stations are located in Brookings County, SD. The weather stations at Brookings and Oak Lake sites are surrounded by grass, whereas the weather station at Volga is surrounded by corn fields.

The values were calculated using weather dataset using the Penman-Monteith equation [8, 54] as follows:where is the alfalfa reference (mm day−1), is the slope pressure versus air temperature curve (kPa°C−1), is the net radiation at the crop surface (MJ m−2 day−1), is the soil heat flux at the soil surface (MJ m−2 day−1), is the mean air temperature at 1.5 to 2.5 m height (°C), is the mean daily wind speed at 2 m height (m s−1), is the saturation vapor pressure of the air (kPa), is the actual vapor pressure of the air (kPa), is the psychrometric constant (0.0671 kPa°C−1), is the vapor pressure deficit (kPa), is the numerator constant (1600 K mm s3 Mg−1 day−1), is the denominator constant (0.38 s m−1) for alfalfa reference, and is the coefficient constant (m2 mm MJ−1).

All weather data were subjected to quality control (QC) prior to being used in any calculations as suggested by Allen et al. [8] and ASCE-EWRI [54]. Hourly QC included the following weather variables such as solar radiation, air temperature (maximum and minimum), wind speed, and air vapor pressure deficit. Therefore, accurate estimations of depend on the quality weather data.

2.5. Atmometers

Three atmometers were used to measure daily . One automated atmometer Model E (ETgage Company, Loveland, Colorado, USA) was placed adjacent to the Oak Lake weather station. The automated atmometer Model E was connected to the automated Oak Lake weather station controlled by a CR1000 datalogger (Campbell Scientific, CSI, Logan, UT, USA), where the evaporated data were recorded every 5 minutes. Two manual atmometers were located adjacent to the Brookings and to the Volga weather stations, respectively. The target corn field was located within 5 km from the weather stations. Manual atmometers were manually recorded every morning between 8:15 and 8:30 AM at the Brookings and Volga sites, respectively. The evaporated water from the green canvas in manual atmometers was measured as the difference between the observed water levels on consecutive days [55].

All atmometers were covered with a number 54 green canvas which simulate evaporation rates of alfalfa reference crop. The atmometers were installed on a vertical wooden post using metal brackets and with the top of the ceramic evaporation surface 1.0 m above the ground surface. The atmometers were surrounded by grass at the Brookings and Oak Lake sites, while at the Volga site the atmometer was surrounded by rainfed corn fields (<5 m) in all directions. Due to rodent damage to the canvas, the wooden post at Volga site was replaced by a metal rod in early August (DOY 217) as shown in Figure 3 which prevented further damage and resulting loss of data.

The atmometer observation period was from May 17, 2016 (DOY 138), to September 18, 2016 (DOY 262); during this period the atmometers were refilled two times with distilled water. Distilled water was used in the atmometer reservoir to prevent accumulation of solutes in and on the top of the plate that can decrease the porosity of the plate and affect the evaporation rates [34].

2.6. Development of Crop Coefficient Curves

curves were developed for each corn field at three sites based on the alfalfa reference crop coefficient using appendix method from Manual 70 [3]. This method uses percent of time from planting to effective cover and days after effective cover to harvest for calculations. In our study, the effective cover for corn fields occurred around 55 days after emergence (DAE). Thus the effective cover was used such as a reference point to calculate local values. Local values for different corn growth stages were calculated and adjusted using DAE. According to Irmak et al. [56] the DAE is more accurate because it ignores the period prior to emergence and is more closely related to the corn growing period, from emergence until physiological maturity.

The curves generated in this study for different corn fields (Figure 6) were multiplied by the obtained from atmometers to estimate (-atm) and compare it with estimated with the METRIC model (-METRIC).

2.7. Statistical Analysis between ETa-METRIC and ETa-Atm

Statistical comparison between -METRIC and -atm was established using a simple linear regression. For the simple regression the model was , where is -atm and is -METRIC. Other statistical evaluations such as mean bias error (MBE) (see (10)), root mean square error (RMSE) (see (11)), coefficient of determination () (see (12)), and Willmott index of agreement “” (see (13)) were used to determine agreement between -atm and -METRIC [57]. where is the observation number, is the estimated value with the METRIC model, is estimated value using atmometer, and the bars above the variables indicate averages.

3. Results and Discussion

3.1. Precipitation and Reference Evapotranspiration (ETr)

The cumulative precipitation values for the growing period were 450 mm, 497 mm, and 380 mm for Brookings, Volga, and Oak Lake, respectively. In 2016 the cumulative precipitation for the three sites was greater than the average rain (±360 mm) that typically falls during the growing season (April-October). The precipitation events had good distribution during the corn growing season due to the major events occurring in development stage, vegetation stage (V5) (June, 17), and tassel stage (VT) (July, 10) (Figure 4).

Daily values of from atmometers (-atm) varied from 0.5 to 10, 0.5 to 9.5, and 0.5 to 7.6 mm day−1 for Brookings, Volga, and Oak Lake, respectively (Figure 4). from Penman-Monteith equation (-PM) varied from 1.3 to 9.1 mm day−1 for Brookings, 1.0 to 10.4 mm day−1 for Volga, and 1.3 to 9.6 mm day−1 for Oak Lake (Figure 4). The highest -atm values recorded in the three sites were in early June (DOY 161) and the lowest values were in early September (DOY 249). The highest -PM values registered in the three sites were in early June (DOY 157) and the lowest values were in middle of September (DOY 259). Even so, moderate correlations between -PM values and -atm values in the three sites were found with of 0.64, 0.59, and 0.67 for Brookings, Volga, and Oak Lake sites, respectively (Figure 5).

3.2. Development of Crop Coefficient Curves

The trends of for each corn field at three different sites during the growing season are shown in Figure 6. The curves showed similar tendencies for all corn fields, where values increased from initial stage (vegetation stage (V3)) to mid-season stage (VT). In this period the values increase as a function of time between 10% of crop cover to 100% of effective cover. In the mid-season the remains constant ( = 1.0), while in the late season the values gradually decreased indicating the crop senescence. At the end of the season (reproductive stage (R6)) the values are low again ( = ±0.6).

The curves depend of vegetation index, soil water content, weather conditions, crop variety, and growing degree days [8, 56, 58, 59]. In this study, curves for corn fields presented little variability because rainfall events, emergence days, and air temperature were almost homogenous in our study area.

The maximum values observed in this study were similar to the (from alfalfa reference) values reported by Djaman and Irmak [58], who reported maximum values from 50 to 70 DAE in corn with rainfed treatment. Also, Wright [53] found maximum values at the 100% of effective full cover for a corn field. However, our values are different from those reported by other researchers [60, 61]. They found the peak values (±1.0) from late July to early August (±70 DAE) for corn fields planted in south central Nebraska.

3.3. ETa Maps and Variation of ETa throughout Growing Season

Spatial and temporal distribution of maps during the growing season were generated by the METRIC model for Brookings, Volga, and Oak Lake (data not showed). However, Figure 7 shows two maps for Brookings site, one map displays high -METRIC value (dark blue color) in a corn field (white rectangle) at mid-season (DOY 194), and another map shows low -METRIC value (light green color) at the end of the season (DOY 258). Generally, high ET rates are related to high crop water demands, which normally occurred in the mid-season period, while low ET rates occurred in early and late growing season when the crop has low vegetation and leaf senescence, respectively. Daily maps helped to explain the variability of crop water requirements throughout the growing season.

The maps developed by the METRIC model in this study were similar to other ET maps generated by the METRIC model and reported by Chávez et al. [13], Santos et al. [24], Folhes et al. [18], Droogers et al. [62], Healey et al. [15], Zhang et al. [20], and Carrillo-Rojas et al. [63], where they reported the spatial and temporal distribution of daily for different crops including corn. In other situations, Chávez et al. [13] reported maximum values (14.1 mm day−1) due to high wind speed values (7.0 m s−1) at the time of satellite overpass in corn field in Texas High Plains, USA.

Figure 8 shows the variation of -METRIC values throughout the growing season at three sites. These results are in agreement with [13, 20, 62]; they reported high ET values at mid-season stage and low ET values at the end of the growing season in agricultural crops using the METRIC model.

3.4. ETa Correlations between the METRIC Model and Atmometer

In the METRIC model was taken from ten randomly selected pixels in three corn fields at three locations for each satellite image, while for the atmometer method was the result of multiplying (from atmometer) by values generated in this study for each corn field.

The linear relationship between -METRIC values and -atmometer (-atm) values in three fields at three sites is shown in Figure 9. Brookings and Volga sites demonstrated good distribution of points around the 1 : 1 line, whereas Oak Lake shows that the points were distributed below the 1 : 1 line; this means that -METRIC values were higher than -atm values except on DOY 154. Even so, strong relationships were observed for Brookings and Oak Lake, and a good relationship was observed for Volga (Table 2). In addition, the sum for all corn fields, the -METRIC values, and the -atm values correlated well (Table 3).

In general, the difference between -METRIC (5.36 mm day−1) and -atm (4.95 mm day−1) at three sites was approximately 8%. The coefficient of determination () and index of agreement (“”) were 0.87 and 0.84, respectively. The corresponding MBE was 0.41 mm day−1 and RMSE was 0.65 mm day−1 (Table 3). According to RMSE value, this can be acceptable assuming an average daily -METRIC value of 5.36 mm day−1 and average daily -atm value of 4.95 mm day−1. The METRIC model has been tested with other methods for estimation. In this study METRIC was compared to atmometers, where the relationship between them was good ( = 0.87) indicating that the METRIC model is a useful tool to estimate at field scale. It should be noted that the comparison between METRIC and the atmometers was done for cloud-free conditions. While daily ETrF values are generated from METRIC enabling the construction of daily ET estimates throughout the growing season, those values were not part of this analysis.

Similar results in (±0.86) were reported by Healey et al. [15], Morton et al. [64], Gordillo et al. [19], French et al. [28], and Liebert et al. [21], who compared estimated with the METRIC model and measured with Bowen Ratio Energy Balance System (BREBS) and Eddy Covariance (EC) methods. Similar results ( = ±0.85) were found in measured with atmometers by other researchers [32, 65], although different results (low = ±0.70) were reported by Chen and Robinson [66] and Lamine et al. [41]. All these authors compared measured using atmometer covered with a No 54 green canvas (alfalfa reference) with estimated using agrometeorological data under different weather conditions.

3.5. ETa Differences between the METRIC Model and Atmometer

The difference between the daily estimated with the METRIC model (-METRIC) and estimated by atmometer (-atm) is presented in Figure 10. Negative values indicated that the -METRIC estimates are lower than -atm, while positive values indicated that the -METRIC estimates exceed -atm.

In Brookings, the daily difference ranged from −0.95 to 1.32 mm day−1, found in field 1 (DOY 154) (V3) and field 2 (DOY 202) (R1), respectively. The more negative values were presented early in the growing season (DOY 154) (V3) due to high value (0.51) used with atmometer method compared to the low value (0.38) used in the METRIC model method. In addition, on DOY 178 (V7) corn field 3 shows negative value (−0.4 mm day−1); this is attributed to hailstorm, which occurred on DOY 169 (V5) nine days before that satellite image overpass. The high positive values (DOY 178, 194, and 202) were related to high wind speed values (>4 m s−1) at the time of satellite overpass (red columns in Figure 12). On DOY 218, 234, and 258 the difference between -METRIC and -atm was small (±0.5 mm day−1).

At the Volga site, the daily difference varied between −1.93 and 1.33 mm day−1. These values were found in field 1 for DOY 194 (VT) and for DOY 178 (V7), respectively. The higher positive values were during the development stage (DOY 178) (V7) when the corn was 1 m tall. The lower negative values were during the mid-season (DOY 194) (VT) when the crop was 2 m tall. This discrepancy was due to not only the wind speed values but also recorded in atmometer on DOY 194 (VT). The recorded in this date was one of the highest values registered during the corn growing season (Figure 4, Volga). After DOY 202 (R1) the difference between -METRIC and -atm was minimal (<0.6 mm day−1), because daily average wind speed values were less than 0.8 m s−1 and -PM were around of 10% lower than the -atm. It is important to mention that at the Volga site the weather station was surrounded by corn fields 2 m tall and the wind speed is reduced by the corn height. So, when low wind speed is used to estimate using the P-M equation, the resulting are too low. For this reason the wind speed did not affect the difference between -METRIC and -atm from mid-season to late-season period.

At the Oak Lake site, the daily difference ranged from −0.62 to 2.61 mm day−1, reported for field 3 (DOY 154) (V3) and for field 1 (DOY 194) (VT), respectively. The negative values were found on DOY 154 (V3) for the three fields; this is attributed to that fact that the -atm was calculated using equal to 0.5, while -METRIC used () equal to 0.35, indicating an overestimation of 30% with -atm method. The higher positive (overestimated) values were observed on DOY 194 (VT) (2.61 mm day−1). At the Oak Lake site, the -METRIC values tend to be higher than the -atm values in almost all corn growing season. This noticeable difference is due to the high wind speed values registered throughout the growing season. These high values of wind speed may be attributed to the elevation of the weather station (574 m above sea level), which is 13% higher than Brookings and Volga elevations. Also, at the Oak Lake site the weather station and automated atmometer were located in smooth hill. On the other hand, we observed different values on DOY 178 (V7) between corn fields. This difference is attributed to hailstorm, which affected the canopy cover in the fields 2 and 3 (Figure 10, Oak Lake).

Similar differences were found by Choi et al. [26], who reported ET difference between −2.2 and 2.5 mm day−1 for different land cover types using METRIC model and trapezoid interpolation model (TIM). They found high discrepancy in ET due to low values of elevation; also they reported that as elevation increases the TIM model slightly overestimates the METRIC ET. On the contrary, low ET difference (±1.0 mm day−1) was reported by Chávez et al. [13], who compare ET estimated by METRIC and ET derived from soil water balance in irrigation corn in Texas. Also, Healey et al. [15] and Hankerson et al. [16] reported that ET difference ranged from −1.0 to 1.0 mm day−1 between the METRIC model and BREBS method in different crops.

In our study, daily differences were attributed to high wind speed values (>4 m s−1) at time of satellite overpass. Generally, as the wind speed increases, the difference increases (Figure 11). Figure 11 shows that at the Volga site (columns with dark gray color) the difference between -METRIC and -atm was minimal after mid-season stage; however at the Oak Lake site (black columns) the difference was major due to high wind speed.

3.6. Hourly Wind Speed at Three Sites

Hourly average wind speed (m s−1) for each overpass date at the Brookings, Volga, and Oak Lake sites is shown in Figure 12. The wind speed recorded at Volga is similar to Brookings and Oak Lake early in the season but then is quite a bit lower later during the growing season. This is because Volga weather station is too close to the corn fields, and the wind speed is reduced by the height of corn crop. On the contrary, at the Oak Lake site higher average wind speed values were found, especially at time of satellite overpass (red column in Figure 12). Also, the Oak Lake weather station recorded the maximum wind speed values throughout the season.

In Maui Island, USA, for example, Anderson et al. [67] and Zhang et al. [20] reported high variation in evapotranspiration due to high wind speeds values in sugarcane fields. Our results are similar to results reported by Westerhoff [68], who found that as the wind speed increases the values slightly increase. In addition, Gleason et al. [65] reported ET underestimation with high wind speed conditions. Mokhtari et al. [69] reported that as the wind speed increases decreases; also they concluded that the METRIC-based ET is highly sensitive to surface temperature but less sensitive to wind speed values.

Based on the results from our study, values from atmometer need to be adjusted during the windy days. The adjustment factors (average ratio of -atm to -METRIC) for Brookings, Volga, and Oak Lake were 0.83, 0.87, and 0.68, respectively. These adjustment factors can be used to adjust the -atm values to get close estimates to the -PM values on windy days (>4 m s−1). These adjustment factors are necessary to correct -atm values to obtain accurate estimations. On the other hand, measured with atmometer can be used by the METRIC model instead of derived from weather data for estimating .

4. Conclusions

The objective was to compare estimated from satellite-based remote sensing METRIC model to estimated with atmometers.

Spatial and temporal values were observed at different crop growth stages using the METRIC model. The values derived from METRIC were higher than values estimated with atmometers. However, the -METRIC and -atm correlated well with coefficient of determination () of 0.87 and index of agreement (“”) of 0.84. The main difference in estimation with both methods was attributed to high wind speed values (>4 m s−1) at the time of satellite image overpass. According to our results, the -atm values need to be adjusted during the windy days to get close estimates to the -PM values to get accurate estimations. In this study adjusted -atm values were not validated.

Based on our results the METRIC model was sensible to wind speed values at the time of satellite overpass. A limitation of satellite derived ET is that is estimated only for clear sky days, while from atmometer is estimated for any weather condition. However, a limitation of from atmometer is that is estimated only at field scale or few miles (±20) from atmometer location, whereas with remote sensing is estimated at field, local, and regional scales.

In conclusion the results of this study can be used by researchers and producers to estimate actual evapotranspiration and improve irrigation water management at regional and field scales using satellite-based remote sensing METRIC model method and at field scale using atmometer method.

Disclosure

This article is a part of a doctoral thesis: http://openprairie.sdstate.edu/cgi/viewcontent.cgi?article=2785&context=etd.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Arturo Reyes-González, Todd Trooien, and Jeppe Kjaersgaard conceived and designed the research. Arturo Reyes-González and Jeppe Kjaersgaard analyzed satellite images and graphs. Arturo Reyes-González, Todd Trooien, and Jeppe Kjaersgaard contributed to the data analysis and discussion of the results. Laurent Ahiablame and Christopher Hay provided a critical revision of the manuscript. All authors contributed to reviewing the final paper.

Acknowledgments

The first author would like to acknowledge the National Council for Science and Technology of México (CONACYT) and the National Institute of Forestry, Agriculture, and Livestock Research (INIFAP) for funding his doctoral scholarship. Additional funding is provided by the South Dakota Agricultural Experiment Station, South Dakota Water Resources Institute, and South Dakota Corn Utilization Council. The authors would like to thank Ryan Vanderleest and Nathan Edwards for providing meteorological datasets.