A remote sensing approach for estimating the location and rate of urban irrigation in semi-arid climates

Urban irrigation is an important component of the hydrologic cycle in many areas of the arid and semiarid western United States. This paper describes a new approach that uses readily available datasets to estimate the location and rate of urban irrigation. The approach provides a repeatable methodology at 1/3 km² resolution across a large urbanized area (500 km²). For this study, Landsat Thematic Mapper satellite imagery, air photos, climatic records, and a land-use map were used to: (1) identify the fraction of irrigated landscaping in urban areas, and (2) estimate the monthly rate of irrigation being applied to those areas. The area chosen for this study was the San Fernando Valley in Southern California.

Identifying irrigated areas involved the use of 29 satellite images, air photos, and a land-use map. The fraction of a pixel that consists of irrigated landscaping (F_irr) was estimated using a linear-mixture model of two land-cover endmembers (selected pixels within a satellite image that represent a targeted land-cover). The two endmembers were impervious and fully-irrigated landscaping. In the San Fernando Valley, we used airport buildings, runways, and pavement to represent the impervious endmember; golf courses and parks were used to represent the fully irrigated endmember. The average F^irr using all 29 satellite scenes was 44%. F_irr calculated from hand-digitizing using air photos for 13 randomly selected single-family-residential neighborhoods showed similar results (42%).

Estimating the rate of irrigation required identification of a third endmember: areas that consisted of urban vegetation but were not irrigated. This "nonirrigated" endmember was used to compute a Normalized Difference Vegetation Index (NDVI) surplus, defined as the difference between the NDVI signals of the irrigated and nonirrigated endmembers. The NDVI signals from irrigated areas remains relatively constant throughout the year, whereas the signal from nonirrigated areas rises and falls seasonally due to precipitation. The areas between airport runways were chosen to represent the nonirrigated endmember. Water-delivery records from 65 spatially-distributed single-family neighborhoods, consisting of nearly 1800 homes, were correlated with the NDVI surplus. The results show a strong exponential correlation (r² = 0.94).

In the absence of water-delivery records, which can be difficult to obtain, a surrogate was identified: the landscape evapotranspiration rate (ET_L). ET_L was used to scale NDVI surplus (which is dimensionless) to irrigation rates using an exponential scaling function. The monthly irrigation rates calculated from satellite and climatic data compared well with irrigation rates calculated from actual water-delivery data using a paired Wilcoxan signed-rank test (p = 0.0063).

Identification of F_irr at the pixel scale, along with identification of the irrigation rate for a fully-irrigated pixel, allows for mapping of urban irrigation over large areas. Maps showing the location and rate of monthly irrigation for the San Fernando study area were computed for January and August 1997.