Assessment of Landsat-Based Evapotranspiration Using Weighing Lysimeters in the Texas High Plains

: Evapotranspiration (ET) is one of the largest data gaps in water management due to the limited availability of measured evapotranspiration data, and because ET spatial variability is di ﬃ cult to characterize at various scales. Satellite-based ET estimation has been shown to have great potential for water resource planning and for estimating agricultural water use at ﬁeld, watershed, and regional scales. Satellites with low spatial resolution, such as NASA’s MODIS (Moderate Resolution Imaging Spectroradiometer), and those with higher spatial resolution, such as Landsat (Land Satellite), can potentially be used for irrigation water management purposes and other agricultural applications. The objective of this study is to assess satellite based-ET estimation accuracy using measured ET from large weighing lysimeters. Daily, seven-day running average, monthly, and seasonal satellite-based ET data were compared with corresponding lysimeter ET data. This study was performed at the USDA-ARS Conservation and Production Research Laboratory (CPRL) in Bushland, Texas, USA. The daily time series Landsat ET estimates were characterized as poor for irrigated ﬁelds, with a Nash Sutcli ﬀ e ﬃ ciency (NSE) of 0.37, and good for monthly ET, with an NSE of 0.57. For the dryland managed ﬁelds, the daily and monthly ET estimates were unacceptable with an NSE of − 1.38 and − 0.19, respectively. There are various reasons for these results, including uncertainties with remotely sensed data due to errors in aerodynamic resistance surface roughness length estimation, surface temperature deviations between irrigated and dryland conditions, poor leaf area estimation in the METRIC model under dryland conditions, extended gap periods between satellite data, and using the linear interpolation method to extrapolate daily ET values between two consecutive scenes (images).


Introduction
Evapotranspiration (ET) is an important component of the water cycle and the surface energy balance, where water changes phase from liquid to vapor, and a change in energy takes place [1]. The development of reliable long-term estimates of ET is needed to improve agricultural water use efficiency (WUE). WUE is defined as the assimilated carbon amount as biomass/grain produced per consumed unit of water by a crop. Paulson [2] reported that climate elements and wind conditions affect regional and seasonal ET estimates. For the irrigation scheduling decision-making process, knowledge of ET variability is important for proper water resources management. This is especially true in arid regions, where crop water requirements exceed rainfall, and irrigation is essential for crop production.
To sustain agricultural production in regions that are dependent on irrigation, blue water (surface water) is extracted from streams and aquifers (groundwater) to help sustain the crop and to meet water demand. Over time, farmers are required to dig deeper wells and extract water from deeper aquifers to ensure adequate water supply for their irrigation needs [3]. Unfortunately, the recharge rate of some aquifers is not sufficient to meet this demand, and historical declines in underground water table have occurred due to extensive water withdrawal [4]. In some cases, this has resulted in land subsidence and localized infrastructure damage [3]. ET is the main driver in irrigation water management and planning. The more accurate the ET estimates, the more likely associated management strategies are to achieving potential crop yields. ET is also an important component in estimating soil water to facilitate improved water use efficiency.
Remote sensing-based ET models are effective for crop water requirement estimations at the field and regional scales [5]. Several ET algorithms have been developed to utilize airborne and satellite data for irrigation scheduling and management purposes. ET can be measured over a surface using the Bowen ratio (BR), the eddy covariance (EC), and lysimeter systems at the field-scale. In all of these cases, spatial variability does not apply, since each method represents a very small scale measurement (typically less than 150 m in an agricultural setting). In addition, these methods only provide a single, averaged value that may not adequately capture the variability across a region. Satellite-based ET models produce regional-scale crop water use [6]. Many remote sensing-based ET algorithms have been developed, assessed, and widely used for estimating regional ET [6][7][8][9].
Park et al. [10], Jackson [11], and Choudhury et al. [12] reported that regional and watershed-scale spatially distributed ET can be better represented using remote sensing compared to traditional ET estimation methods. There are typically two approaches that are used to provide remote sensing ET estimates: (1) the land surface energy balance (EB) approach, and (2) the reflectance-based crop coefficient (Kc) and reference ET approach. The first approach is based on ET being a change of the state of water, using available energy in the environment for vaporization [13].
Multiple satellite platforms are available to use with energy balance ET models, such as Land Satellite (Landsat), The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Geostationary Operational Environmental Satellite (GEOS), moderate resolution imaging spectroradiometer (MODIS), and others, based on surface energy balance [14]. However, the spatiotemporal resolution is complex and cannot be combined with only one satellite. For daily estimates of ET with a high spatial resolution (field-scale), results depend on data assimilation, model applicability, and accuracy [15]. ET calibration based on remote sensing ET estimates has been conducted for a single source method such as the Surface Energy Balance Algorithm for Land (SEBAL) and the Mapping Evapotranspiration with Internalized Calibration model (METRIC) [16,17], and two source methods such as the Two-Source Energy Balance model (TSEB) [18,19]. Energy balance models based on thermal remote sensing include TSEB [18,20], SEBAL [8], METRIC [6,17], and several others [7,21].
A detailed assessment of existing EB models [7,22] stated that daily ET estimates varied 3 to 35% in comparison to Bowen ratio and EC ET measurements. Error sources included (a) modeling uncertainties, and (b) measurement errors and discrepancies in model-measurement scales. Other studies that assessed multiple airborne high-resolution sensor platforms indicated good agreement with measured ET [23].
The BEAREX08 experiment [15] was a robust remote sensing experiment that involved mass balance measurements of ET using four weighing lysimeters [24,25]. A wide array of instrumentation was installed for supportive data, including neutron probe (NP) access tubes for soil water measurements, multi-level canopy temperature measurements, and above and below canopy irradiance measurements. The purpose of these measurements was to provide as many measured data as possible to reduce estimations and provide more points of comparison. These data were used to address the lack of studies where EB model estimates of ET are compared with mass balance measurements, reducing uncertainty sources in EB models. Gowda et al. [7] reported that few remotely sensed ET estimates are used in irrigation scheduling and in-field management due to the absence of daily data with field-scale resolution. Many methods have been proposed to overcome this issue, such as use of infrequent, high spatial resolution data, including Landsat with 60 m spatial resolution, 16-day temporal resolution, in combination with lower spatial resolution and more frequent data, such as the Moderate Resolution Imaging Spectroradiometer (MODIS) with 1000 m spatial resolution and 1-day temporal resolution. The Disaggregated Atmosphere-Land Exchange Inverse (DisALEXI) model demonstrated this concept [14,26,27], but no evaluation of this approach for management has been successful thus far.
Spatial and temporal resolution problems could be resolved using aircraft. However, cost, data processing, and lack of experienced users have prevented their widespread use in providing such platforms as potential imagery sources for water management. The EB method typically provides an instantaneous value of ET, and interpolated to daily ET using the evaporative fraction method, [28], or reference ET approach [6,29]. Each ET estimation approach has its uncertainties; however, measurement data quality assessment and quality control reduce such modeling uncertainties [30,31]. Multi-time scale remotely-sensed ET estimation accuracy is essential for water management and irrigation planning. Extensive studies have been conducted to assess the METRIC model performance on irrigated fields [32][33][34][35]; however, none of these studies evaluated the model performance under dryland conditions for extended periods with various crops and temporal resolutions [36]. This study assessed daily interpolated ET estimation accuracy using the linear interpolation method on dryland and irrigated lysimeters for a ten year period with various crops at multiple time scales. Such assessment is crucial for water management policies, researchers, hydrologists, and irrigators when using such technologies for real time irrigation decisions.
The objective of this research was to assess the daily, seven-day running average, monthly, and seasonal satellite-based ET estimation accuracy versus large weighing lysimeter ET measurements in the Texas High Plains under irrigated and dryland conditions during the growing and non-growing seasons.

Study Site
Data used in this study from 2001-2010 were obtained at the USDA-ARS Conservation and Production Research Laboratory (CPRL) in Bushland, Texas, USA (35.19 • N, 102.10 • W). Four square fields were selected for this study; each field was~4.7 ha. Four large precise weighing lysimeters were installed towards the center of each field [24,25]. Two lysimeters were managed as irrigated (NE and SE), and the other two lysimeters were managed as dryland (NW and SW). The irrigated lysimeter field was equipped with a linear-move irrigation system with Nelson sprinklers (Nelson Irrigation Corporation, Walla Walla, WA, USA) [37]. Irrigation scheduling was performed using neutron probe data, and in 2001 50% of crop water requirements (CWR) was replenished for cotton (deficit irrigation), and 100% of CWR was added for the remaining study period. Crop management data were collected and summarized by the CPRL [38,39].
Leaf area index (LAI) data were collected from the study site between 2001 and 2010. Plant samples were collected from locations within fields near the lysimeters biweekly during the growing season through destructive plant sampling. Leaves were separated from stems, and the average LAI values were obtained from at least three samples. Plant samples were not collected from the lysimeters, as the sampling was destructive and would impact lysimeter ET measurements. A digital scanning bed leaf area meter (model LI-3100) was used to measure the leaf area index of the leaf samples. The LAI values were calculated as the ratio of the upper side leaf area (m 2 ) to the ground area (m 2 ) [22,40].
For analysis simplicity, one lysimeter was selected for each management condition (irrigated and dryland). Statistical assessment was achieved for the estimated ET values for the dryland lysimeter (NW) and the irrigated lysimeter (NE). The soil characteristics for the study field are deep, well-drained Pullman silty clay loam (fine, mixed, superactive, thermic torrertic paleustoll) [40]. The local climate is classified as semi-arid, with large daily air temperature variations. Cotton, soybean, grain, and silage sorghum, sunflower, and cotton were the predominant crops for the research fields during the ten year study period [38,41].

Bushland Evapotranspiration and Agricultural Remote Sensing (BEARS)
Bushland Evapotranspiration and Agricultural Remote Sensing (BEARS) is an image processing and geographic information system (GIS) software developed by researchers at the USDA ARS CPRL in Bushland, TX, used for deriving hourly, daily, and seasonal ET maps, and other energy exchanges between land and atmosphere using Landsat 5,7, and 8 [42]. It is an open-source Java software, Version 1.0.1 available for download at: https://data.nal.usda.gov/dataset/bushlandevapotranspiration-and-agricultural-remote-sensing-system-bears-software. The software allows for custom models and equations but provides the option to select one of five default energy balance-based ET methods: Mapping Evapotranspiration at High Resolution with Internalized Calibration (METRIC), Surface Energy Balance Algorithm for Land (SEBAL), Surface Energy Balance System (SEBS), Two-Source Model (TSM), and Simplified Surface Energy Balance (SSEB) [42].

Image Analysis
In this study, the METRIC model was used to analyze Landsat satellite imagery, producing ET time-series datasets for the study period. The METRIC model description and required inputs based on Landsat satellite datasets are summarized in the literature [6,41]. Several hourly and daily outputs were obtained at the completion of each image analysis. Four daily outputs were obtained including daily ET (mm day −1 ), evaporation fraction (EF) (unitless), leaf area index (LAI) (m 2 m −2 ), and normalized difference vegetation index (NDVI) (unitless). Hashem et al. [41] assessed hourly estimation accuracy for ET (mm hr −1 ), net radiation (Rn) (W m −2 ), soil heat flux (Go) (W m −2 ), and surface air temperature (Ts) ( • C). A detailed description of the hourly assessment procedure was summarized [38,41].

Landsat Satellite Dataset and Processing
Landsat 5 TM cloud-free images were selected for analysis, which were obtained through Earth Explorer  Table 1. Processed image date and day of the year (DOY) are summarized in Table 1.

Landsat ET Gap Filling
Linear interpolation methods were used to estimate the satellite-based ET during the gap period using the satellite-based daily evaporation fraction between consecutive images. The gap between consecutive images was split equally into two periods of eight days; the first eight days used the first image evaporation fraction value, and the second eight-day period used the next image evaporation fraction value.
To estimate the actual ET time series, the evaporation fraction (k c ) for each day was multiplied by the reference ET [43]. The daily reference ET was estimated using the ASCE−2005 [44] standardized reference evapotranspiration equation. Hourly ET was calculated for the study period, and the daily ET was obtained by accumulating the hourly ET in a given day. The average frequency and the maximum gap periods (days) were calculated using Equations (1) and (2) Table 2. Areas of interest (AOIs) of 3 by 3 grids (30 m spatial resolution), with a total of 9 pixels and a surface area of 8100 m 2 each, were created, and each lysimeter was located at the center of each AOI. ArcGIS 10.2.2 model builder tool software was utilized to extract the average ET for each AOI.

Dry and Wet Pixel Determination
Obtaining the most accurate wet (cold) and dry (hot) agricultural pixels is a significant step that is essential to produce accurate ET maps [34,45]. Wet and dry pixels represent the extreme conditions per scene. The manual selection method was used, along with the surface temperature and NDVI histogram distribution, to obtain these pixels for each image. NDVI and surface temperature threshold values are summarized in Table 3 to obtain the hot and dry pixel characteristics. The hot pixel is defined as bare agricultural soil with high temperature and low ET value, and the cold pixel is defined as cultivated agricultural soil with low temperature and high ET value [6,46]. The histogram distribution for surface temperature and NDVI help to determine the range of temperature and NDVI values per image.

Statistical Analysis
Statistical comparisons were performed including multiple statistical coefficients, including the coefficient of determination (R 2 ), regression line slope and intercept, root mean square error (RMSE), percent of root mean square error (% RMSE error), mean bias error (MBE), and Nash Sutcliff efficiency (NSE), between the measured and the calculated satellite-based ET. A detailed description of each coefficient and its range and interpretation of values were provided [38,41,47]. The mean bias error (MBE) and % MBE provided a means to determine the deviation between the measured and satellite-based estimates, with MBE = 0 indicating no bias in the estimation within the satellite data. The MBE and % MBE were calculated using Equations (3) and (4), respectively.
where y t = the ith measured value,ŷ t = the ith simulated value, and n = the total number of observations.
where x obs is mean of measured values.

Dryland Lysimeter ET Estimation
The northwest dryland lysimeter daily measured ET, Landsat ET, and seven-day running average values were plotted from 2001 to 2005, as shown in Figure 1a, and from 2006 to 2010 as shown in Figure 1b. The seven-day running average was used for visual comparison. The 1:1 graph between the daily measured and calculated Landsat ET can be seen in Figure 2, and the daily and monthly summary statistics are summarized in Table 4. The monthly average ET values were plotted, as shown in Figure 3. Detailed growing and non-growing season summary statistics are reported in Table 5.      The measured mean ET was 1.3 mm d −1 , and the mean Landsat ET was 1.7 mm d −1 . Due to the daily ET value deviations, the summary statistics provided a poor match during the comparison period with an R 2 value of 0.01, NSE value of −1.38, RMSE of 1.8, and an RMSE error of~144.3%, which is considered a high error value. These statistical parameters indicated that there was no correlation between measured and calculated values [48].
The growing season dryland measured LAI was plotted versus calculated Landsat estimates for the days where Landsat images were available, as shown in Figure 4. It can be clearly seen that the Landsat LAI underestimated the LAI during the study period for all cultivated crops. Dryland conditions in Bushland, TX typically have much less vegetation than irrigated areas. The lower amounts of plant biomass allow the soil background to be more prominent at the Landsat spatial resolution. The increased soil background can skew the reflectance and LAI, causing lower values. Another potential source of errors is due to uncertainties with the remotely sensed datasets under dryland conditions, due to surface roughness length and ET extrapolation methods that have been incorporated in the METRIC model [32,38,49]. Chavez et al. [7] reported that the METRIC model estimation error was (0.7 ± 0.9 mm d −1 ), explaining that the variations were due to errors associated with the surface roughness length and aerodynamic resistance. Another potential source of uncertainties is that the METRIC model uses a SURFACE NDVI vs. LAI relationship, where it is generated by fitting a generalized equation to six LAIs compared to NDVI functions [35], defined in the MODIS LAI backup [45], indicating a higher source of LAI estimation under dryland conditions.  Figure 5a,b presents the northeast irrigated lysimeter daily measured ET, Landsat ET, and seven-day running average from 2001 to 2010. The 1:1 graph between the daily measured and calculated Landsat ET is shown in Figure 6, and the daily and monthly summary statistics are summarized in Table 6. The growing and non-growing season summary statistics are reported in Table 7. The daily 1:1 graph is shown in Figure 6, and the monthly ET is shown in Figure 7.      The daily mean ET was 2.4 mm d −1 and 2.4 mm d −1 for the Landsat and measured, respectively. The summary statistics improved with the irrigated field and provided a weak correlation with an R 2 value of 0.38, NSE of 0.37, RMSE of 2.1, and RMSE~86.4%.

Irrigated Lysimeter ET Estimation
The growing season irrigated measured LAI was plotted versus Landsat estimates for the days where Landsat images were available during the year, as shown in Figure 8. Landsat better estimated LAI under irrigated conditions compared to the dryland conditions. Consequently, higher NDVI values were obtained [38], producing higher LAI values for irrigated fields, and resulted in better ET estimates.

Dryland Daily ET Comparison
The relationship between measured and Landsat ET for the dryland lysimeter showed significant deviation with periods of both over and underestimation of ET throughout the year for the entire study period. The satellite-based LAI was assessed versus the measured LAI (Figure 4), and the LAI assessment summary is summarized in Hashem [38]. The daily time series ET deviations were related to errors in LAI estimation [38,41,47], where Landsat LAI estimates were significantly lower than measured LAI during the growing season for the dryland lysimeter. The higher the NDVI values, the more the LAI values increase, resulting in greater ET values.
In 2002, 2005, and 2009, the lysimeter field was fallow, and Landsat overestimated ET in each of the three years (Figure 1), and these results agree with Allen et al. [46]. Cotton was cultivated in 2001 and 2008, and Landsat estimates of ET closely matched the measured ET at the beginning of each year. However, towards the end of 2001, Landsat significantly underestimated ET due to low NDVI values and, consequently, underpredicting LAI [38]. ET data in 2004 and 2007 were omitted from the analysis, as the Landsat data overestimated the ET compared to the measured ET due to the large gap period in Landsat data and linear interpolation method used to fill the gap. In 2003, when sorghum was cultivated, the satellite-based-ET overestimated measure ET in both the beginning and towards the end of the year, and underestimated towards the middle of the growing season.
A detailed statistical analysis was performed for the growing and non-growing seasons (Table 5), where the growing season was defined as the days between planting and harvest. Monthly statistics showed better statistical performance, with monthly RMSE and NSE of −0.19 and 1.2 compared to values of −1.38 and 1.8 mm for the daily assessment [38,46,50]. The RMSE during the growing season was greater compared to the non-growing season (Table 5), with values almost double for the growing season compared to the non-growing season due to low measured ET values during the non-growing season. Hence, there was less variation between the measured and satellite-based ET values. However, the %RMSE error was higher during the non-growing season than the growing season, and these results agree with Allen et al. [46].
The satellite was able to distinguish between bare soil and vegetation in the field, providing useful information on when the field was fallow versus when a crop was growing. However, the overall LAI estimation from Landsat was lower than the measured LAI for all cultivated crops during this study under dryland conditions. Potential reasons for the LAI undercalculations are the water stress during the growing season producing low NDVI values under dryland conditions, uncertainties with aerodynamic resistance surface roughness length [36], long gap periods, and using the linear interpolation method to generate daily ET time series [38].

Irrigated Daily ET Comparison
The relationship between measured and Landsat ET for the irrigated lysimeter provided overall better agreement compared to the dryland field [38,41,51,52]. The Landsat ET estimates were closely matched most of the year, except the middle of the growing seasons, during the peak crop water requirements. The satellite-based approach underestimated ET toward the middle of the growing season for cotton and soybeans, and overestimated the ET early and late during the growing season.
A detailed statistical analysis was performed for the daily and monthly ET ( Table 6). The irrigated daily ET estimates were considered poor with an NSE of 0.37, RMSE of 2.1 mm d −1 , and %RMSE of 86.4%. However, there was a statistical improvement with the monthly ET values with an NSE of 0.57, RMSE of 1.5 mm d −1 , and % RMSE of 56.7%. Similar to the dryland lysimeter, the RMSE during the growing season was greater compared to the non-growing season (Table 7), with values almost double for the growing season compared to the non-growing season due to low ET measured values during the non-growing season. Hence, there was less variation between the measured and satellite-based ET values. However, the %RMSE error was higher during the non-growing season than the growing season, and these results agree with Allen et al. [46]. Allen et al. [6] illustrated that the use of reference ET considers the advective effects on a METRIC model performance, which can make the METRIC ET overestimate the ET from irrigated fields, exceeding daily net radiation in arid and semi-arid conditions. Allen et al. [46] reported that daily ET had the largest differences due to ET fluctuating the most during the growing season, and the monthly and season ET lumped most of the daily variations [38,46], and this is in agreement with the current study results.
Similarly, for the dryland lysimeter, the deviation between Landsat and measured ET was related to higher LAI estimation [32,36,38,47], advective condition effects under irrigated conditions [6,36], and extended gap periods [38]. In addition, most of the studies conducted evaluated the ET on the current scene (image) days with minimal EB closure errors [36,49,53], and no studies evaluated the extrapolated daily ET assessment for dryland conditions with clumped crops [36]. However, the irrigated field difference magnitude was far less than for the dryland field.
Landsat LAI estimates were better for the irrigated lysimeter, as the METRIC model performance was affected with the wet and cold pixel determination, and the METRIC model performed better with full canopy (full irrigated), compared to dryland (partial canopy) [32,36,38,49]. As irrigated fields produced more vegetation vigor, higher NDVI values were obtained [36,38], and consequently higher estimates of LAI were obtained and resulted in better estimates of Landsat ET for areas managed under irrigated conditions (Figure 8) [32,36,38,41].
The overall Landsat LAI estimation somewhat matched the measured LAI for most cultivated crops from 2001 to 2010 (2004 and 2007 omitted due to large gaps in clear Landsat data during the growing season). Three indicators that the satellite imagery was able to differentiate between irrigated (full canopy) and dryland fields (partial canopy) as well as identify the growing season were as follows: 1.
The estimated LAI values for the irrigated field were much larger than that of the dryland field. 2.
All LAI values were zero in the beginning and end of each year, and this reflects that the field was bare soil. However, there were LAI values recorded during the growing season for the same field, providing useful information on when the field was fallow versus when a crop was growing.

3.
The magnitudes of NDVI values for irrigated fields were higher than those for the dryland fields [36,38].
The reason behind this is likely that LAI is better estimated for the irrigated field than the dryland field [36,38] due to more vegetation coverage, resulting in higher NDVI values and consequently ET values. In 2007 and 2009, when forage sorghum and sunflowers were cultivated, respectively, the LAI estimated using Landsat was slightly lower than the measured LAI for the irrigated lysimeter.

Conclusions
Remote sensing-based ET estimation is considered a promising tool for irrigation water management. However, uncertainties associated with satellite-based ET estimation still exist, especially with various remotely sensed platforms due to variations in spatial and temporal resolution. In this study, satellite-based ET was evaluated using Landsat under semi-arid conditions in Texas under irrigated and dryland conditions. Ten years of lysimeter measured ET data were used in this study. The Landsat-based ET overestimated the measured ET early and late in the growing season and underestimated ET during the peak of the growing season. The daily and monthly ET for the dryland lysimeter was unacceptable with negative NSE (−1.38 and −0.19), indicating there was no correlation between the estimates and measured ET; however, the daily and monthly ET for the irrigated lysimeter values showed better statistics with an NSE of 0.37 and 0.57, respectively. Seasonal ET showed more variations during the growing season compared to the non-growing season, because higher ET values were estimated during the growing season.
Under dryland conditions, there was significant LAI underestimation compared to the measured LAI values due to water stress during the growing season. LAI plays a significant role in evapotranspiration; where greater values of NDVI were obtained, consequently greater LAI was obtained under irrigated conditions, resulting in more ET for irrigated conditions. There are several reasons behind uncertainties of LAI and ET estimation, including the following: (1) METRIC model uncertainties with partial canopy estimates, (2) dryland plants' rapid modification of LAI based on available soil water (partial canopy), and (3) uncertainties with aerodynamic resistance surface roughness length as well as surface temperature deviations between irrigated and dryland conditions. Extended gap periods are another significant challenge, and the selection of the filling method can account for ET estimation errors. In this study, gap periods reached up to 184 days in 2004, and the minimum was in 2008 with 40 days. The linear interpolation method was utilized to extrapolate the daily ET estimates between every two consecutive images in this study.
More satellite-based ET assessment under arid and semi-arid conditions is required, where the magnitude and frequency of precipitation are erratic, and irrigation is the only source under arid conditions to replenish crop water needs. With advances in remote sensing, more frequent satellite imagery will be available, with high spatial resolutions. Other extrapolation methods should be considered to generate daily time-series ET datasets. This would likely improve overall ET estimation accuracy by improving the overall spatial and temporal resolution.
Future research opportunities that include the assessment of ET relationship with crop physiology, yield, and yield components (number of flowers, grain quality, etc.) would provide potential information on crop response under dryland and irrigated conditions. Economic analysis of commodity market prices would be another research project due to groundwater decline in the Ogallala aquifer.