Remote Sensing and GIS Techniques for Assessing Irrigation Performance: Case Study in Southern California
Publication: Journal of Irrigation and Drainage Engineering
Volume 144, Issue 6
Abstract
This paper presents the potential of remotely sensed data in addressing spatially distributed irrigation equity, adequacy, and sustainability. The surface energy balance algorithm for land (SEBAL) was implemented to map actual evapotranspiration (ET) over an irrigation district in southern California. Potential ET was also mapped based on the Priestley–Taylor method, modified to account for the effect of horizontally transported energy on enhancing/suppressing ET. Remotely sensed products were integrated with ground-based data in a Geographical Information System (GIS) environment to quantify several irrigation and drainage performance indicators. The among- and within-field coefficients of variation of actual ET were comparable to previous studies, suggesting that water consumption was uniform across the irrigation district. The relative ET was high, indicating that irrigation supply was adequate. The extensive network of open drains was also found to be functioning at an optimal level according to the results of two performance indicators based on the magnitude and uniformity of groundwater depth.
Introduction
Historically, the western United States has been known for its arid climate, low precipitation, and frequent droughts, which have made water management a very complex issue in this part of the world. Water governance in the western United States must be performed in such a way as to not only provide for increasing demands from competing sectors, but also to protect ecosystems and critical habitat for flora and fauna. In addition, new concerns on the possible consequences of climate change have added to the complexity of this already challenging task.
Accounting for about 65% of total water withdrawals, irrigation has been the largest use of freshwater in the United States since 1950. Not surprisingly, the majority of agricultural withdrawals (86%) and irrigated area (75%) were in the 17 contiguous western states (Hutson et al. 2004). Therefore, it is of great importance to accurately determine how much water needs to be diverted for irrigation and how it is partitioned into different consumptive and nonconsumptive uses. Quantifying water balance components at irrigation scheme scales has a wide variety of applications, including but not limited to initiating and evaluating water conservation practices, improving irrigation scheduling (Santos et al. 2008), developing irrigation modernization scenarios (Isidoro et al. 2004), assessing biophysical and economical water productivities (Teixeira et al. 2008), and managing soil salinization (Faci et al. 1985; Khan et al. 2006; Marlet et al. 2009).
For irrigation schemes in arid and semiarid regions such as those in the western United States, the most significant water balance component is crop evapotranspiration (ET). Although accurate point measurements of crop ET have been extensively used in managing agricultural water resources, recently developed remote sensing techniques have acceptable levels of accuracies (Gowda et al. 2008) while providing spatially distributed data that enable a higher level of analysis at the scale of individual pixels. Another advantage of air- or space-borne remote sensing data is their objectiveness, an important characteristic that can revolutionize developing and standardizing benchmarks for comparing irrigation schemes from around the world.
When combined with water balance information, remotely sensed ET can also be utilized in evaluating the performance of irrigation and drainage systems (Ahadi et al. 2013). Irrigation performance is traditionally evaluated based on point measurements. A major caveat of this approach is that the results can only provide one value representing the entire area downstream of the measurement location. Considering that water application and partitioning are highly variable from field to field, the results of traditional methods are less useful as the scale of the study increases from field to district and basin. Recent developments in remote sensing and GIS techniques have made it possible to assess scheme-wide performance on a pixel-by-pixel basis (Bastiaanssen and Bos 1999).
Despite the numerous benefits of remote sensing techniques, there is still a gap between research projects and practical applications of these techniques in the real world (Ambast et al. 2002). Bastiaanssen and Bos (1999) and Bastiaanssen et al. (2000) stated that more demonstration projects are needed to increase the level of awareness among water managers about the potential of air- and space-borne imagery and promote the use of these tools for improving irrigation performance. The study presented herein was carried out over an irrigation district in semiarid southern California. A satellite-based energy balance approach was implemented to map ET. The results were integrated in an ArcGIS environment with ground measurements of agrohydrological parameters to estimate several performance indicators addressing irrigation equity and adequacy at different spatial and temporal scales. A new irrigation sustainability indicator was also introduced and another previously defined indicator was evaluated for the first time.
Methods and Materials
Study Area
The Palo Verde Irrigation District (PVID) is located in Imperial and Riverside counties in southeastern California. With about of territory, PVID was privately developed in 1925 to serve local water users. The majority of irrigated fields within PVID are under surface irrigation (laser-graded borders and furrows), receiving water from the Colorado River through an extended network of irrigation canals (400 km) and open drains (230 km). The alluvial soils of PVID were deposited over the years by Colorado River floods. The medium texture of these soils allows them to hold a considerable amount of water and to be easily drained. Alfalfa, cotton, small grains, and vegetables are the most dominant crops in this irrigation district. Fig. 1 shows a false-color satellite image of the PVID and the location of the diversion dam that directs Colorado River water, as well as the network of irrigation canals and open drains.
Input Data
Several types of point and distributed data were collected or estimated to allow for evaluation of irrigation performance across the study area. These included water inflow/outflow, precipitation, groundwater depth, actual ET, and potential ET.
Water Inflow/Outflow
Colorado River water is diverted into PVID’s main canal using a small diversion dam. Since the district is relatively flat and the fields are mostly blocked-end furrows and borders, the surface runoff is not significant, and any generated runoff is directed toward the drains. Therefore, water outflow from PVID consists of two components, drainage and canal spills. The high density of deep open drains (about 0.5 km per every of land), along with the medium texture of PVID’s alluvial soils, significantly enhances the movement of water from the root zone toward the drains. All of the drains merge and discharge into a main outfall drain at the downstream end of PVID. The United States Geological Survey (USGS) measures the flow rates of water diversion, drain discharge, and all canal spills on a daily basis and reports them online. These data were downloaded from the USGS (2010b).
Precipitation
PVID benefits from a network of 32 rain gauges distributed uniformly across the entire district. This means that there is a rain gauge for every of PVID’s cultivated land (). The point measurements of rain gauges were interpolated in ArcGIS to develop precipitation maps for each event.
Actual Groundwater Depth
PVID has a network of 260 piezometers distributed uniformly across the entire district (one piezometer for every of PVID’s cultivated land). Depth to groundwater was measured at each piezometer once every month. These point measurements were collected, analyzed for quality, and then interpolated (Krigin method) in the ArcGIS environment to generate monthly maps of actual groundwater depth ().
Actual Evapotranspiration
The surface energy balance algorithm for land (SEBAL) Idaho implementation was used to map actual evapotranspiration () across the study area (Bastiaanssen et al. 1998b). The performance of SEBAL has been assessed and validated in numerous applications (Bastiaanssen et al. 1998a; Ramos et al. 2009), showing a high accuracy at field and catchment scales (Bastiaanssen et al. 2005). In SEBAL, net radiation () is estimated through quantifying all of the incoming and outgoing short- and long-wave radiation components. Once is determined, soil heat flux (G) is modeled as a ratio of net radiation and a function of surface temperature and fraction of vegetation cover.
The innovation of SEBAL is in mapping sensible heat flux (H). This parameter is estimated for two extreme pixels and then interpolated for all other pixels. One extreme pixel (cold pixel) represents a well-watered crop at full cover, where the available energy is used mainly for changing the phase of water from liquid to gas, leaving a negligible amount of energy to be used in generating a temperature gradient (H about zero). The other extreme pixel (hot pixel) represents a dry, bare soil, where the available energy is used primarily in heating the atmosphere (maximum H). The final step of SEBAL involves estimating the latent heat flux (LE) as the remainder of the surface energy balance equation
(1)
Space- or airborne imagery—as input data to models such as SEBAL—provides only a snapshot of LE at the time of overpass. As a result, instantaneous estimates of ET need to be scaled up to longer periods (daily and seasonal) for most practical purposes. In the earlier versions of SEBAL, instantaneous ET was extrapolated to daily values using the evaporative fraction (EF or ). This concept is based on the assumption that the ratio of instantaneous ET to instantaneous available energy () is constant during the day (Brutsaert and Sugita 1992; Crago 1996), especially under cloud-free conditions (Zhang and Lemeur 1995). Once this ratio is determined, daily ET could be calculated by multiplying the EF ratio and the daily value of available energy.
Although the EF technique has provided reliable results in many studies (Gowda et al. 2008), its accuracy decreases in arid regions because of the occurrence of afternoon advection. Trezza (2002) modified the EF ratio by replacing available energy with alfalfa reference ET (), which encompasses the effect of energy imported from dry neighboring areas. Upscaled daily ET estimates using this new method () have been shown to be significantly improved (Romero 2004; Allen et al. 2007a, b). Alternatively, grass-based reference ET () has been used in extrapolating instantaneous ET values ( method). Colaizzi et al. (2006) compared estimates of five different upscaling techniques with measurements of precision weighing lysimeters at Bushland, Texas, where strong advection of heat is common. For cropped surfaces, the method worked better than and EF methods. Chávez et al. (2008) also indicated that, under advective conditions, outperforms and EF methods. Since PVID is surrounded by dry desert regions, was selected as the up-scaling method in this study.
All cloud-free Landsat TM5 scenes collected between January 2008 and January 2009 were downloaded from the website of the USGS Global Visualization Viewer (GLOVIS) (USGS 2010a). These scenes were all processed using the LPGS processing system, which results in 60-m resolution for the thermal band and 30-m resolution for other bands. PVID is located on the overlap zone of two Landsat paths (38 and 39; Row 37). This allowed the acquisition of six extra scenes from Path 39 in addition to 15 scenes from Path 38 (total of 21 scenes). Required weather parameters (including ) were acquired from a California Irrigation Management Information System (CIMIS) weather station located within PVID (Station 135).
Potential Evapotranspiration
The concept of potential ET () has been used in estimating some performance indicators. However, the exact definition of this parameter and the approach to estimating it have been debated among researchers. In this study, was considered as the level of water consumption that a specific crop could potentially reach. This is similar to the concept of standard (nonstressed) crop ET () in that both are estimated for stress-free conditions, when the amount of crop water consumption is limited by energy availability. Two distributed estimation methods have been used before in evaluating performance over large and heterogeneous irrigation schemes: available energy () method, implemented by Roerink et al. (1997) and Bastiaanssen et al. (1996), and the Priestley–Taylor approach (Priestley and Taylor 1972), utilized by Bastiaanssen et al. (2001), Bandara (2003), and Karatas et al. (2009). A major concern with both methods is that only radiative energy is considered and the advective energy is neglected (Glenn et al. 2007). As explained before, irrigated areas in arid/semiarid regions frequently experience advective conditions. When water is available, the effect of advection on enhancing ET would be large in these regions, resulting in significant errors if neglected.
In this study, was estimated on a pixel-by-pixel basis using the Priestley–Taylor (PT) equation (Priestley and Taylor 1972)where = Priestley–Taylor parameter; = slope of the water vapor saturation curve; = psychrometric constant; and other terms are defined as previously. To account for the effect of advection, the above equation was calibrated for the local condition of PVID through a process similar to what was implemented by Diaz-Espejo et al. (2005). In this process, the original value of 1.26 was modified to force the PT equation to estimate ET rates similar to what was estimated over a reference grass surface, using the Penman-Monteith equation (Allen et al. 1998).
(2)
Performance Indicators
Several performance indicators (PI) were studied at different spatiotemporal resolutions. These PIs can be arranged in three groups based on what aspect of the system they address:
•
Equity: water consumption uniformity (WCU);
•
Adequacy: relative ET (RET) and depleted fraction (DF); and
•
Sustainability: drainage ratio (DR), drainage distribution uniformity (DDU), and relative groundwater depth (RGD).
Water Consumption Uniformity
The WCU is a measure of irrigation equity (Bastiaanssen and Bos 1999), assessed through estimating the coefficient of variation for at two levels, namely, within () and among () irrigation units. In order to extract required statistical parameters for each field in PVID, a crop classification layer was used in ArcGIS as a mask layer. Developed from Landsat TM5 imagery by the U.S. Bureau of Reclamation (USBR), this classification layer defined field boundaries and provided information on the cropping pattern in 2008. Since the thermal band of Landsat TM5, used in mapping , has a coarser resolution than the other six bands, an inner buffer of 60 m was applied to field boundaries to avoid including edge pixels in analysis. The mean and standard deviation of were obtained for every field that had an area larger than 4.0 ha after applying the buffer. This reduced the number of studied units from more than 2,000 fields to 1,485.
Relative Evapotranspiration
RET was estimated as the ratio of to . This indicator was used as a measure of irrigation adequacy to investigate the presence and severity of water shortages. The interpretation of results was accomplished at two levels: including all fields and limiting to fields at full canopy cover. Elimination of the fields that were not close to full cover was achieved by setting a threshold value of 0.65 on the soil adjusted vegetation index [SAVI (Huete 1988)]. All fields with an average SAVI less than this threshold were excluded from the second level of analysis.
Depleted Fraction
DF is a rather new term for an old concept that has been given several names such as irrigation efficiency, water efficiency, water application efficiency, and consumptive use coefficient. All of these terms provide information on the fraction of water that has been depleted from available resources (Jensen 2007). According to Bos et al. (2005), DF is defined aswhere = gross precipitation over the study area; and V = volume of applied () or diverted () water. For the sake of developing benchmarks and comparing irrigation schemes, it is important to differentiate between applied and diverted water, because operational and maintenance constraints may impose a large variation in the amount of canal spills. Estimated DF would be gross () or net (), if or were used in the above equation, respectively.
(3)
Drainage Ratio
The DR is another performance indicator that provides information on what portion of applied water has left the study area in the form of drainage. DR is inversely related to (Bos et al. 2005)
(4)
The leaching fraction (LF) necessary for maintaining a favorable salt balance can be taken as a critical value of DR. If estimated DR is less than required LF, soil salinization may become problematic in the future. A DR value greater than LF means that applied water can be reduced without affecting the current level of agricultural production. Adjusting irrigation management to reduce DR may not have any effect on the quantity of water available to downstream users, but the effect on preserving its quality could be significant.
Drainage Distribution Uniformity
The DDU was introduced in this study for the first time. DDU is assessed by evaluating among and within field coefficient of variation ( and ) of the actual groundwater depth. A low DDU is only an indication of uniform depth to groundwater and does not necessarily mean that water is at a depth that allows an adequate root respiration. Therefore, this indicator should be studied along with other drainage performance indicators, such as relative groundwater depth (explained subsequently). In areas where the crop water requirement is met by controlling the level of groundwater, DDU could serve as a measure of traditional irrigation distribution uniformity.
Relative Groundwater Depth
Bos (1997) proposed the RGD as the ratio of actual groundwater depth () to critical groundwater depth (). Although it was defined over a decade ago, RGD has never actually been estimated over any irrigated area, to the best of the authors’ knowledge. Since for most crops, groundwater should be kept at levels below the root zone to avoid any negative effects caused by waterlogging and/or soil salinization, it is assumed that is equal to the effective root depth for every crop group. Table 22 in Allen et al. (1998) suggests a range of effective depths for each crop. In this study, the upper limit of this range was selected as . These values were assigned to each field in PVID, using the crop classification layer. Then, the minimum value of was obtained for each field from the interpolated maps of depth to groundwater for the months of February and September, when groundwater is farthest from and closest to the surface, respectively. RGD estimates based on minimum and maximum (as opposed to average values) are more conservative because they represent the portion of the field that has the highest potential for waterlogging issues. Values less than unity indicate that the water table is within the crop root zone, a situation that irrigation managers should avoid for most crop types.
Results and Discussion
Input Data
Water Inflow/Outflow
The volume of water diverted from the Colorado River for irrigating PVID lands between January 2008 and January 2009 was or 2,479 mm when divided by the total cultivated area. The monthly average flow rates ranged between 14.5 and in January and June, respectively. With a volume of about , canal operational spills were 11.5% of the total annual water diversion. Monthly average canal spills measured at discharge structures varied from 2.3 to in January and June, respectively. It is very important to take the amount of canal spills into account when evaluating system performance, since this water returns to the river with little to no changes in quality. Monthly average flow rates of drainage discharge measured at the downstream end of the outfall drain ranged between in January and in July, which shows a lag in comparison with the peak flow rate of diversion. With a volume of about , the total annual amount of drainage was 40% of diversion. This portion increases to over 45% if only applied water (diverted minus spills) is considered.
Precipitation
Twenty-five precipitation events were recorded during the study period, all in forms of rainfall. Out of all events, only two had an average depth greater than 10 mm. The total annual precipitation depth was 71 mm, underscoring the aridity of this region.
Actual Groundwater Depth
Depth to groundwater was affected by irrigation intensity. The deepest levels were observed in February 2008, after irrigations had been ceased for some time. Diversion to the main canal was also shut down for maintenance purposes during the first week of January. The average was 3.3 m in February, which was equal to the average depth of open drains. After February, groundwater levels rose until they reached 2.8 m from the surface in September and October 2008, which was toward the end of the main irrigation season. From this point, declined until it reached 3.2 m in January 2009. Fig. 2 demonstrates cumulative frequency curves developed from interpolated maps of for the 3 months of February, June, and September 2008. These months were selected to show three different levels of groundwater from lowest to highest. According to Fig. 2, no PVID groundwater levels were closer than about 2.0 m to the soil surface, an indication of a successfully functioning drainage system.
Actual Evapotranspiration
On average, the annual over PVID’s cultivated land () was 1,286 mm. This was similar to the annual ET estimate of 1,268 mm, reported by Taghvaeian and Neale (2011) over the same period based on a district-wide water balance approach. Annual field-level ranged from less than 70 mm for fallow fields to more than 2,000 mm for alfalfa fields. Water consumption of other dominant crops was less than alfalfa at 1,063 and 1,320 mm for Sudan grass and cotton, respectively. Fig. 3 presents a SEBAL-ET map on September 17, 2008, as an example of the 21 processed images.
Crop coefficients () were determined by dividing SEBAL- by the grass-based reference ET estimated at the CIMIS weather station. The frequency distribution and cumulative frequency of values over PVID are demonstrated in Fig. 4 for two satellite overpass dates of April 8 and July 13, 2008. On both dates, the most frequent occurred at about 1.1, which is about the of alfalfa and cotton at peak water use. Since alfalfa has a year-round growing season, the April 8 distribution was mainly controlled by alfalfa. The frequency distribution changed on July 13 for values smaller than 0.2 and larger than 1.0, but it had minimal differences for values in between these two limits. This evolution can be attributed to the growth of other crops such as cotton. Cotton is planted in mid-March in PVID, so its water consumption is still small in early April ( values less than 0.25). By mid-July, however, most of the cotton fields are at full cover, consuming water at the highest rates. As a result, the low- frequencies move to a range between 1.0 and values larger than 1.2 appeared to belong to recently irrigated fields.
Potential Evapotranspiration
In estimating , the PT parameter, , was adjusted to include the effect of advective energy on enhancing and/or suppressing water consumption. Adjusted values were estimated through dividing the Penman-Monteith by the equilibrium ET. Fig. 5 shows adjusted daily values for the period from March 1 to October 1, a time frame that represented the main agricultural growing season and when advective energy was expected to be significant.
Out of 215 days of analyzed data, only one day had an α value less than 1.26. On this date (May 24) PVID rain gauges recorded an average precipitation depth of 14.4 mm, the most intensive rainfall during the study period. Therefore, it is highly probable that moist and cool air from surrounding deserts converged over PVID, resulting in water consumptions lower than what was predicted by the available energy. Except for 3 days, adjusted values were all above 1.4, suggesting that advection was a major contributor to water consumption of PVID crops. Several previous studies reported values greater than 1.4 (Diaz-Espejo et al. 2005; Pereira and Nova 1992) and 1.5 (Li and Yu 2007), offering that such values indicate the occurrence of enhanced advective conditions.
Performance Indicators
Water Consumption Uniformity
Within field coefficient of variability () of annual ranged from 1.0 to 80%, with an average and median of 7.0 and 3.4%, respectively. Although variabilities as high as 80% were detected, 85% of the fields had a lower than 10%. Using maps of , PVID managers can locate the remaining 15% of the fields with values larger than 10% and focus their attention on these flagged fields rather than the entire district. This is an example of how the distributed nature of remotely sensed indicators can save time and energy in evaluating and improving scheme-wide irrigation performance. The low values of average suggest that, overall, water application at the field level was uniform across PVID. A primary reason could be that almost all PVID fields are precisely leveled, using modern laser grade-control systems. Among the fields, the coefficient of variation () for PVID fields was 38.2%, but this large value does not translate into poor irrigation equity since PVID fields are highly diverse in crop type, growing season, and water requirement. Table 1 compares the CV estimates found in this study with the reported values in the literature.
Publication | Model | RS platform | Number of overpasses | Study area | Study unit | Main crops | (%) | (%) |
---|---|---|---|---|---|---|---|---|
This study | SEBAL | Landsat TM | 21 | California, U.S. | 1,485 fields | Alfalfa, cotton | 38.2 | 7.0 |
Kharrou et al. (2013) | NDVI-ET | Landsat ETM | 8 | Haouz plain, Morocco | 3 zones | Wheat | 17.0 | N/A |
Zwart and Leclert (2010) | SEBAL | Landsat ETM | 12 | Office du Niger, Mali | 5 zones | Rice | 2.4 | 8.9 |
Ahmad et al. (2009) | SEBAL | MODIS | 19 | Rechna Doab, Pakistan | 9 subdivisions | Rice, wheat | 2.4 | 4.0 |
Ahmad et al. (2009) | SEBAL | MODIS | 19 | Rechna Doab, Pakistan | 15 subdivisions | Sugarcane, wheat | 4.9 | 7.5 |
Roerink et al. (1997) | SEBAL | Landsat TM | 1 | Rio Tunuyan, Argentina | 10 secondary units | Orchards, vineyards | 8.6 | N/A |
Roerink et al. (1997) | SEBAL | Landsat TM | 1 | Rio Tunuyan, Argentina | 31 tertiary units | Orchards, vineyards | 6.1 | N/A |
Bastiaanssen et al. (1996) | SEBAL | Landsat TM | 1 | Nile Delta, Egypt | 53 districts | Rice, cotton, maize | 10.0 | 15.0 |
Note: ETM = enhanced thematic mapper; NDVI = normalized difference vegetation index; RS = remote sensing.
Crop-specific water consumption uniformity was also studied over cotton and alfalfa fields. For 22 large cotton fields, and were 3.2 and 8.5%, respectively. The uniformity was also promising for 45 large alfalfa fields in PVID, with and values of 3.1 and 9.4%, respectively. According to Molden and Gates (1990), a of less than 10% is considered good uniformity. Santos et al. (2008) reported a higher variability for 13 cotton fields in the Genil-Cabra Irrigation Scheme (GCIS) in southern Spain ( and ). GCIS fields were under a modern pressurized system with an on-demand delivering regime, while PVID fields were under surface irrigation and a modified-demand water delivery.
Relative Evapotranspiration
Average RET for all PVID fields was 0.97 when an value of 1.26 was used. This is a high RET, suggesting that on average all studied fields were consuming water at rates similar to their potential. When considering only full-cover fields, RET increased further to 1.49, meaning that these fields had water consumption about 50% larger than their potential. Adjusting α for advective conditions reduced the average RET of all fields and full-cover fields to more realistic estimates of 0.69 and 1.06, respectively (Fig. 6). Estimating values that are 6% larger than on average is acceptable since was calibrated with grass-based , while the ET of alfalfa (the most dominant crop in PVID) can be up to 20% larger than grass under the same agroclimatological conditions.
Assuming that was equal to available radiative energy (), Roerink et al. (1997) reported average RET values of about 0.6 for several secondary and tertiary units in the Rio Tunuyan irrigation scheme in Argentina. Roerink et al. suggested that RET values of 0.75 and higher were satisfactory for irrigated agriculture. Using the value of 1.26 (PT original), Karatas et al. (2009) estimated an average RET of 0.7 for the period of May to September over the Lower Gediz Basin (LGB) in Turkey, which is similar to PVID in terms of irrigation method and main crops. It may seem appropriate to compare RET estimates of LGB with the RET estimate of PVID before adjusting α. Such a comparison would favor PVID, with an average RET of 0.97. However, the effect of advection on ET may have been different for these two irrigation schemes. Although the amounts of precipitation during the study periods were similar for both regions (about 30 mm), advective conditions may have been less prevalent under the Mediterranean climate of LGB. If that was the case, comparing the average RET of LGB with the average RET of PVID after adjusting α may be more reasonable (0.7 versus 0.69). The same logic applies to comparing the RET estimate of PVID to that of the Nilo Coelho (0.76), a pressurized irrigation scheme in Brazil (Bastiaanssen et al. 2001). The results of the present study were also similar to those reported for irrigated winter wheat in central Morocco, with an average RET of 0.68 (Kharrou et al. 2013), as well as a flood-irrigated region in Uzbekistan where RET decreased from over 1.0 to less than 0.75 as the distance to the canal intake increased from less than 15 km to over 95 km (Conrad et al. 2013).
Depleted Fraction
Gross and net DF were estimated on a daily basis and averaged for each month in 2008. ranged between 0.33 in December and 0.58 in May, with an annual average of 0.49. As expected, values were larger, varying from 0.41 in December to 0.64 in May, with an annual average of 0.55. Fig. 7 demonstrates the intermonthly variation of and in 2008.
The temporal pattern of monthly was different for irrigated areas in Uzbekistan, with values decreasing from April to June (minimum of about 0.4) and then increasing until September when they exceeded 0.8 (Conrad et al. 2013). For the Nilo Coelho irrigation scheme, Bastiaanssen et al. (2001) reported values ranging between 0.4 in April and 0.85 in November, with an annual average of 0.61. For the condition of this irrigation scheme (perennial orchards under pressurized irrigation), values beyond the operational range of 0.7–1.0 were found to result in at least 10% reduction in the yield. Three irrigation districts in central Morocco had January–February values similar to those estimated in the present study, but significantly larger values (above 0.8) in March and April (Kharrou et al. 2013). For the semiarid LGB in Turkey, Karatas et al. (2009) estimated values larger than 3.0 in May and September. The authors claimed that a depletion that was three times larger than diverted water was achieved through consumption of soil moisture stored in the root zone. values dropped to as low as 0.28 in July, resulting in an average of 0.69 for the whole study period (May to September) over the LGB. Another study found similar average values of 0.60 and 0.72 for cotton and grape fields, respectively, in a subdivision of the LGB (Droogers and Bastiaanssen 2002).
Drainage Ratio
The average annual DR over PVID was 0.45, ranging from 0.36 in May to 0.59 in December. This was smaller than the average seasonal DR reported for three irrigation subdivisions in Uzbekistan, ranging from 0.55 to 0.65 (Conrad et al. 2013). Since no soil salinity study has been carried out in PVID, there is no information on leaching requirements to be used in assessing DR estimates. However, it may be safely assumed that the leaching requirement was not greater than 10–15% in PVID, since its water source had a low salinity (electrical conductivity of the Colorado River water measured by USGS above the diversion dam never exceeded in 2008). If this was the case, the fact that drainage was 45% of applied water may be an indication of overirrigation in this area. Reducing the amount of applied water may have significant effects on downstream ecosystems because of the elevated levels of nutrients in irrigation return flows.
Drainage Distribution Uniformity
The and of the actual groundwater depth were estimated for PVID fields for the 2 months of February and September, representing periods when groundwater was farthest from and closest to the soil surface, respectively. The was 13.4% in February and 15.5% in September. The was significantly smaller, with average values of 0.9 and 1.1% for February and September, respectively. The uniformity of within and across PVID fields suggests that the extent and density of the existing drainage network is appropriate for removing excess water from the top soil. However, the adequacy of drainage depth cannot be assessed based on DDU and requires estimating the next performance indicator.
Relative Groundwater Depth
Similar to DDU, RGD was estimated for the 2 months of February and September. Not a single field had a RGD value less than unity in either month, meaning that the groundwater level was always below the maximum depth that crop roots could reach. However, the distribution of RGD values was different between the 2 months (Fig. 8). For example, only 2.5% of the fields had a RGD between 1.0 and 1.5 in February, but this increased to about 27% in September. No significant correlation was found between RGD and distance to drain, estimated in ArcGIS as the distance between the centroid of each field and the closest drain. This means that expanding the existing drainage network will not have any major effect on the depth of drained soil.
Spatially distributed information on drainage performance indicators generated by GIS techniques can significantly assist irrigation managers with locating the fields that are at higher sustainability risks. Focusing only on these fields would result in a more efficient use of resources to improve drainage performance. Fig. 9 demonstrates two examples of such maps, one based on DDU ( of ) and the other based on RGD, where the red color specifies fields that may have drainage issues. Any mathematical combination of two or more drainage performance indicators can also be generated (e.g., by assigning different weights) based on the importance of each factor in future sustainability of crop production.
Conclusions
All available cloud-free Landsat TM5 imagery acquired from January 2008 through January 2009 were collected and processed using the SEBAL to map actual evapotranspiration () over the Palo Verde Irrigation District in southern California. The results were then combined with ground measurements of precipitation, water inflows/outflows, and depth to groundwater in the ArcGIS environment. Consumptive use of water by PVID crops was 52% of diverted water and 7% of the Colorado River discharge () upstream of the Palo Verde diversion dam during the study period.
The data were used to estimate several irrigation and drainage performance indicators. Within-field water consumption was uniform, most probably because of the precise grading of flood-irrigated fields using laser technology. However, 15% of the fields had variability larger than 10%. Using the distributed information on variability, irrigation managers can locate nonuniform fields and focus their attention on investigating possible reasons and potential solutions. Another indicator, the relative ET, revealed that full-cover fields were consuming water at about 6% greater than their potential rate estimated by the PT method. A relative ET of 1.06 indicates that, on average, PVID fields were provided with adequate water.
Three drainage performance indicators were also estimated over PVID to investigate irrigation sustainability. The drainage ratio was larger than typical leaching requirements of most irrigation schemes, suggesting that many fields were overirrigated. This could increase the load of nutrients in downstream ecosystems receiving the return flows, creating environmental challenges. The depth to water table was both uniform and larger than the maximum root depth of dominant crops at all times. According to these indicators, the PVID drainage system was successfully functioning, removing excess water from the root zone in a uniform fashion. The results showed that remote sensing and GIS techniques can be combined to provide a comprehensive picture of irrigation performance across irrigation schemes.
Acknowledgments
This research was funded in part by the U.S. Bureau of Reclamation under a contract with the Alliance of Universities with a subcontract to Utah State University. Additional funding was provided by Utah Agricultural Experiment Station and the Remote Sensing Services Laboratory, Department of Civil and Environmental Engineering, Utah State University. The authors would like to express their appreciation to Mr. Roger Henning and Ms. Paula Hayden (chief engineer and staff of PVID) for kindly providing us with their databases, and for their valuable comments on the results.
References
Ahadi, R., Samani, Z., and Skaggs, R. (2013). “Evaluating on-farm irrigation efficiency across the watershed: A case study of New Mexico’s Lower Rio Grande Basin.” Agric. Water Manage., 124, 52–57.
Ahmad, M. U. D., Turral, H., and Nazeer, A. (2009). “Diagnosing irrigation performance and water productivity through satellite remote sensing and secondary data in a large irrigation system of Pakistan.” Agric. Water Manage., 96(4), 551–564.
Allen, R. G., et al. (2007a). “Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC)—Applications.” J. Irrig. Drain. Eng., 395–406.
Allen, R. G., Pereira, L. S., Raes, D., and Smith, M. (1998). “Crop evapotranspiration—Guidelines for computing crop water requirements.”, FAO, Rome.
Allen, R. G., Tasumi, M., and Trezza, R. (2007b). “Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC)—Model.” J. Irrig. Drain. Eng., 380–394.
Ambast, S. K., Keshari, A. K., and Gosain, A. K. (2002). “Satellite remote sensing to support management of irrigation systems: Concepts and approaches.” Irrig. Drain., 51(1), 25–39.
ArcGIS version 10 [Computer software]. ESRI, Redlands, CA.
Bandara, K. M. P. S. (2003). “Monitoring irrigation performance in Sri Lanka with high-frequency satellite measurements during the dry season.” Agric. Water Manage., 58(2), 159–170.
Bastiaanssen, W. G. M., et al. (1998a). “A remote sensing surface energy balance algorithm for land (SEBAL). 2: Validation.” J. Hydrol., 212, 213–229.
Bastiaanssen, W. G. M., and Bos, M. G. (1999). “Irrigation performance indicators based on remotely sensed data: A review of literature.” Irrig. Drain. Syst., 13(4), 291–311.
Bastiaanssen, W. G. M., Brito, R. A. L., Bos, M. G., Souza, R. A., Cavalcanti, E. B., and Bakker, M. M. (2001). “Low cost satellite data for monthly irrigation performance monitoring: Benchmarks from Nilo Coelho, Brazil.” Irrig. Drain. Syst., 15(1), 53–79.
Bastiaanssen, W. G. M., Menenti, M., Feddes, R. A., and Holtslag, A. A. M. (1998b). “A remote sensing surface energy balance algorithm for land (SEBAL). 1: Formulation.” J. Hydrol., 212, 198–212.
Bastiaanssen, W. G. M., Molden, D. J., and Makin, I. W. (2000). “Remote sensing for irrigated agriculture: Examples from research and possible applications.” Agric. Water Manage., 46(2), 137–155.
Bastiaanssen, W. G. M., Noordman, E. J. M., Pelgrum, H., Davids, G., Thoreson, B. P., and Allen, R. G. (2005). “SEBAL model with remotely sensed data to improve water-resources management under actual field conditions.” J. Irrig. Drain. Eng., 85–93.
Bastiaanssen, W. G. M., Van der Wal, T., and Visser, T. N. M. (1996). “Diagnosis of regional evaporation by remote sensing to support irrigation performance assessment.” Irrig. Drain. Syst., 10(1), 1–23.
Bos, M. G. (1997). “Performance indicators for irrigation and drainage.” Irrig. Drain. Syst., 11(2), 119–137.
Bos, M. G., Burton, M. A., and Molden, D. J. (2005). Irrigation and drainage performance assessment: Practical guidelines, Cabi Publishing, Wallingford, U.K.
Brutsaert, W., and Sugita, M. (1992). “Application of self-preservation in the diurnal evolution of the surface energy budget to determine daily evaporation.” J. Geophys. Res. Atmos., 97(D17), 18377–18382.
Chávez, J. L., Neale, C. M., Prueger, J. H., and Kustas, W. P. (2008). “Daily evapotranspiration estimates from extrapolating instantaneous airborne remote sensing ET values.” Irrig. Sci., 27(1), 67–81.
Colaizzi, P. D., Evett, S. R., Howell, T. A., and Tolk, J. A. (2006). “Comparison of five models to scale daily evapotranspiration from one-time-of-day measurements.” Trans. ASAE, 49(5), 1409–1417.
Conrad, C., Dech, S. W., Hafeez, M., Lamers, J. P. A., and Tischbein, B. (2013). “Remote sensing and hydrological measurement based irrigation performance assessments in the upper Amu Darya Delta, Central Asia.” Phys. Chem. Earth, Parts A/B/C, 61, 52–62.
Crago, R. D. (1996). “Conservation and variability of the evaporative fraction during the daytime.” J. Hydrol., 180(1–4), 173–194.
Diaz-Espejo, A., Verhoef, A., and Knight, R. (2005). “Illustration of micro-scale advection using grid-pattern mini-lysimeters.” Agric. For. Meteorol., 129(1), 39–52.
Droogers, P., and Bastiaanssen, W. G. M. (2002). “Irrigation performance using hydrological and remote sensing modeling.” J. Irrig. Drain. Eng., 11–18.
Faci, J., Aragüeś, R., Alberto, F., Quilez, D., Machin, J., and Arrue, J. L. (1985). “Water and salt balance in an irrigated area of the Ebro River Basin (Spain).” Irrig. Sci., 6(1), 29–37.
Glenn, E. P., Huete, A. R., Nagler, P. L., Hirschboeck, K. K., and Brown, P. (2007). “Integrating remote sensing and ground methods to estimate evapotranspiration.” Crit. Rev. Plant Sci., 26(3), 139–168.
Gowda, P. H., Chávez, J. L., Colaizzi, P. D., Evett, S. R., Howell, T. A., and Tolk, J. A. (2008). “ET mapping for agricultural water management: Present status and challenges.” Irrig. Sci., 26(3), 223–237.
Huete, A. R. (1988). “A soil-adjusted vegetation index (SAVI).” Remote Sens. Environ., 25(3), 295–309.
Hutson, S. S., Barber, N. L., Kenny, J. F., Linsey, K. S., Lumia, D. S., and Maupin, M. A. (2004). “Estimated use of water in the United States in 2000.” U.S. Geological Survey Circular, Reston, VA.
Isidoro, D., Quılez, D., and Aragüés, R. (2004). “Water balance and irrigation performance analysis: La Violada irrigation district (Spain) as a case study.” Agric. Water Manage., 64(2), 123–142.
Jensen, M. E. (2007). “Beyond irrigation efficiency.” Irrig. Sci., 25(3), 233–245.
Karatas, B. S., Akkuzu, E., Unal, H. B., Asik, S., and Avci, M. (2009). “Using satellite remote sensing to assess irrigation performance in Water User Associations in the Lower Gediz Basin, Turkey.” Agric. Water Manage., 96(6), 982–990.
Khan, S., Tariq, R., Yuanlai, C., and Blackwell, J. (2006). “Can irrigation be sustainable?” Agric. Water Manage., 80(1), 87–99.
Kharrou, M. H., et al. (2013). “Assessment of equity and adequacy of water delivery in irrigation systems using remote sensing-based indicators in semi-arid region, Morocco.” Water Resour. Manage., 27(13), 4697–4714.
Li, L., and Yu, Q. (2007). “Quantifying the effects of advection on canopy energy budgets and water use efficiency in an irrigated wheat field in the North China Plain.” Agric. Water Manage., 89(1), 116–122.
Marlet, S., Bouksila, F., and Bahri, A. (2009). “Water and salt balance at irrigation scheme scale: A comprehensive approach for salinity assessment in a Saharan oasis.” Agric. Water Manage., 96(9), 1311–1322.
Molden, D. J., and Gates, T. K. (1990). “Performance measures for evaluation of irrigation-water-delivery systems.” J. Irrig. Drain. Eng., 804–823.
Pereira, A. R., and Nova, N. A. V. (1992). “Analysis of the Priestley-Taylor parameter.” Agric. For. Meteorol., 61(1–2), 1–9.
Priestley, C. H. B., and Taylor, R. J. (1972). “On the assessment of surface heat flux and evaporation using large-scale parameters.” Mon. Weather Rev., 100(2), 81–92.
Ramos, J. G., Cratchley, C. R., Kay, J. A., Casterad, M. A., Martínez-Cob, A., and Dominguez, R. (2009). “Evaluation of satellite evapotranspiration estimates using ground-meteorological data available for the Flumen district into the Ebro Valley of NE Spain.” Agric. Water Manage., 96(4), 638–652.
Roerink, G. J., Bastiaanssen, W. G. M., Chambouleyron, J., and Menenti, M. (1997). “Relating crop water consumption to irrigation water supply by remote sensing.” Water Resour. Manage., 11(6), 445–465.
Romero, M. G. (2004). “Daily evapotranspiration estimation by means of evaporative fraction and reference evapotranspiration fraction.” Ph.D. dissertation, Utah State Univ., Logan, UT, 194.
Santos, C., Lorite, I. J., Tasumi, M., Allen, R. G., and Fereres, E. (2008). “Integrating satellite-based evapotranspiration with simulation models for irrigation management at the scheme level.” Irrig. Sci., 26(3), 277–288.
Taghvaeian, S., and Neale, C. M. (2011). “Water balance of irrigated areas: A remote sensing approach.” Hydrol. Processes, 25(26), 4132–4141.
Teixeira, A. D. C., et al. (2008). “Energy and water balance measurements for water productivity analysis in irrigated mango trees, northeast Brazil.” Agric. For. Meteorol., 148(10), 1524–1537.
Trezza, R. (2002). “Evapotranspiration using a satellite-based surface energy balance with standardized ground control.” Ph.D. dissertation, Utah State Univ., Logan, UT, 339.
USGS. (2010a). “Global visualization (GloVis) viewer.” U.S. Dept. of the Interior, Reston, VA.
USGS. (2010b). “Lower Colorado river monthly data reports.” U.S. Dept. of the Interior, Reston, VA.
Zhang, L., and Lemeur, R. (1995). “Evaluation of daily evapotranspiration estimates from instantaneous measurements.” Agric. For. Meteorol., 74(1), 139–154.
Zwart, S. J., and Leclert, L. M. (2010). “A remote sensing-based irrigation performance assessment: A case study of the Office du Niger in Mali.” Irrig. Sci., 28(5), 371–385.
Information & Authors
Information
Published In
Copyright
This work is made available under the terms of the Creative Commons Attribution 4.0 International license, http://creativecommons.org/licenses/by/4.0/.
History
Received: Aug 30, 2017
Accepted: Dec 29, 2017
Published online: Apr 10, 2018
Published in print: Jun 1, 2018
Discussion open until: Sep 10, 2018
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.