A HYDROLOGIC RESPONSE MAP FOR THE STATE OF GEORGIA1

Hydrologic response, defined as the ratio annual direct runoff divided by the annual precipitation, was calculated from 178 years of record on 55 watersheds less than 200 mi² in the State of Georgia. Direct runoff was determined from Geological Survey records by a universal method of hydrograph separation. Regression analysis showed that the effect of area and the deviation of actual from normal annual precipitation can be removed from the response ratio, revealing the average capacity of watershed source areas for releasing or detaining potential flood waters. A map of Georgia reveals a four-fold range in the response ratio between major provinces, and another map of a small mountain watershed shows an eight-fold range over a distance of three miles. The response ratio is proposed as a new mapping unit and its use in watershed planning and education is discussed.     

Weighting Nitrogen and Phosphorus Pixel Pollutant Loads to Represent Runoff and Buffering Likelihoods

Watershed models often estimate annual nitrogen (N) or phosphorus (P) pollutant loads in rural areas with export coefficient (EC) (kg/ha/yr) values based on land cover, and in urban areas as the product of spatially uniform event mean concentration (EMC) (mg/L) values and runoff volume. Actual N and P nonpoint source (NPS) pollutant loading has more spatial complexity due to watershed variation in runoff likelihood and buffering likelihood along surface and subsurface pathways, which can be represented in a contributing area dispersal area (CADA) NPS model. This research develops a CADA NPS model to simulate how watershed properties of elevation, land cover, and soils upslope and downslope of each watershed pixel influence nutrient loading. The model uses both surface and subsurface runoff indices (RI), and surface and subsurface buffer indices (BI), to quantify the runoff and buffering likelihood for each watershed pixel, and generate maps of weighted EC and EMC values that identify NPS pollutant loading hotspots. The research illustrates how CADA NPS model maps and pixel loading values are sensitive to the spatial resolution and accuracy of elevation and land cover data, and model predictions can represent the lower and upper bounds of NPS loading. The model provides managers with a tool to rapidly visualize, rank, and investigate likely areas of high nutrient export.     

THE USE OF FLOOD POTENTIAL INDICES FOR FLOOD PEAK ESTIMATION ON UNGAGED WATERSHEDS1

: The construction of a flood peak index map was attempted for use by hydrologists in the simple format of rainfall maps. Since flood peaks are highly dependent on watershed area, the effect of area was removed. By regression analysis flood peaks of 2.33 and 100-year return periods were found to be proportional to watershed area to the 0.8 and 0.7 powers, respectively. Therefore, indices C_{2 33}= Q_{2 33}/_A0.7 were completed at each gage and plotted on a Pennsylvania map. It was attempted to further remove some of the scatter by regression of C with several other watershed parameters like slope, percent forest cover, and watershed shape, but no significant correlation could be found. The index maps, drawn without attenuation of the scatter, can be used by hydrologists to compute flood peaks as Q = CAⁿ (with n = 0.8 and 0.7 for the 2.33 and 100-year flood peaks, respectively). Flood peak safety factors can be based on visual observation of the index variation in the vicinity of the location for which the flood peak estimate is needed.     

Classification and spatial mapping of riparian habitat with applications toward management of streams impacted by nonpoint source pollution

Management of riparian habitats has been recognized for its importance in reducing instream effects of agricultural nonpoint source pollution. By serving as a buffer, well structured riparian habitats can reduce nonpoint source impacts by filtering surface runoff from field to stream. A system has been developed where key characteristics of riparian habitat, vegetation type, height, width, riparian and shoreline bank slope, and land use are classified as discrete categorical units. This classification system recognizes seven riparian vegetation types, which are determined by dominant plant type. Riparian and shoreline bank slope, in addition to riparian width and height, each consist of five categories. Classification by discrete units allows for ready digitizing of information for production of spatial maps using a geographic information system (GIS). The classification system was tested for field efficiency on Tom Beall Creek watershed, an agriculturally impacted third-order stream in the Clearwater River drainage, Nez Perce County, Idaho, USA. The classification system was simple to use during field applications and provided a good inventory of riparian habitat. After successful field tests, spatial maps were produced for each component using the Professional Map Analysis Package (pMAP), a GIS program. With pMAP, a map describing general riparian habitat condition was produced by combining the maps of components of riparian habitat, and the condition map was integrated with a map of soil erosion potential in order to determine areas along the stream that are susceptible to nonpoint source pollution inputs. Integration of spatial maps of riparian classification and watershed characteristics has great potential as a tool for aiding in making management decisions for mitigating off-site impacts of agricultural nonpoint source pollution.     

UNCERTAINTY ANALYSIS OF RUNOFF ESTIMATES FROM A RUNOFF CONTOUR MAP1

ABSTRACT: The U.S. Environmental Protection Agency (EPA) in cooperation with the U.S. Geological Survey (USGS) conducted an analysis to quantify the uncertainty associated with interpolating runoff to specific sites using a runoff contour map. We interpolated runoff to 93 gaged watersheds from a runoff contour map using (1) hand interpolation to the watershed outlet, (2) a computer interpolation to the watershed outlet, and (3) hand interpolation to the watershed centroid. We compared the interpolated values to the actual gaged values and found that there was a bias in the average interpolated value for runoff estimated at basin outlets, with interpolated values being less than the actual. We found no significant difference between the hand interpolation method and the computer interpolation method except that the computer method tended to have higher variability due to factors inherent to the software used. There were no strong spatial correlations or regional patterns in the runoff interpolations, which indicates that there are no regional biases introduced in the development of the contour map. We determined that we could estimate runoff, on the average, within approximately 8.9 cm (3.5 in; 15 percent) of the measured value using the three methods. The results of this work indicate that runoff contour maps can he used in regional studies to estimate runoff to ungaged systems with quantifiable uncertainty.     

DEFINING REGIONAL POPULATIONS OF LAKES FOR THE ASSESSMENT OF SURFACE WATER QUALITY1

ABSTRACT: Topographic maps are commonly used to define populations of lakes in regional surveys of surface water quality. To illustrate the effect of different maps on that process, we compared the lakes represented on the 1:250,000-scale maps used for the Northeast Region of the Eastern Lake Survey—Phase I (ELS-I) to the lakes on a sample of large-scale maps (1:24,000 or 1:62,500). Lake areas at or near the lower limit of representation delimited “smallest-lake” values for the compared 1:250,000-scale maps. The regional median for these values was 4.5 hectares (ha) and ranged from 0.6 to 24.8 ha. Lake representation is influenced by cartographic limitations such as map scale, age, and complexity as well as the inherent variability of waterbodies (e.g., water level fluctuations or the creation of reservoirs, beaver impoundments, and oxbows). The total number of lakes on large-scale maps increased markedly as lake area decreased. Approximately 15,700 of the estimated 29,000 lakes in the EPA's Northeast Region were 1 to 4 ha in area. Because maps affect the size distribution of lakes included in a regional survey and because lake areas are thought to modify lake chemistry, maps ultimately affect the estimates of regional surface water quality.     

Estimating Watershed Level Nonagricultural Pesticide Use From Golf Courses Using Geospatial Methods1

Abstract: Limited information exists on pesticide use for nonagricultural purposes, making it difficult to estimate pesticide loadings from nonagricultural sources to surface water and to conduct environmental risk assessments. A method was developed to estimate the amount of pesticide use on recreational turf grasses, specifically golf course turf grasses, for watersheds located throughout the conterminous United States (U.S.). The approach estimates pesticide use: (1) based on the area of recreational turf grasses (used as a surrogate for turf associated with golf courses) within the watershed, which was derived from maps of land cover, and (2) from data on the location and average treatable area of golf courses. The area of golf course turf grasses determined from these two methods was used to calculate the percentage of each watershed planted in golf course turf grass (percent crop area, or PCA). Turf‐grass PCAs derived from the two methods were used with recommended application rates provided on pesticide labels to estimate total pesticide use on recreational turf within 1,606 watersheds associated with surface‐water sources of drinking water. These pesticide use estimates made from label rates and PCAs were compared to use estimates from industry sales data on the amount of each pesticide sold for use within the watershed. The PCAs derived from the land‐cover data had an average value of 0.4% of a watershed with minimum of 0.01% and a maximum of 9.8%, whereas the PCA values that are based on the number of golf courses in a watershed had an average of 0.3% of a watershed with a minimum of <0.01% and a maximum of 14.2%. Both the land‐cover method and the number of golf courses method produced similar PCA distributions, suggesting that either technique may be used to provide a PCA estimate for recreational turf. The average and maximum PCAs generally correlated to watershed size, with the highest PCAs estimated for small watersheds. Using watershed specific PCAs, combined with label rates, resulted in greater than two orders of magnitude over‐estimation of the pesticide use compared to estimates from sales data.     

Urbanization and the Loss of Resource Lands in the Chesapeake Bay Watershed

We made use of land cover maps, and land use change associated with urbanization, to provide estimates of the loss of natural resource lands (forest, agriculture, and wetland areas) across the 168,000 km² Chesapeake Bay watershed. We conducted extensive accuracy assessments of the satellite-derived maps, most of which were produced by us using widely available multitemporal Landsat imagery. The change in urbanization was derived from impervious surface area maps (the built environment) for 1990 and 2000, from which we estimated the loss of resource lands that occurred during this decade. Within the watershed, we observed a 61% increase in developed land (from 5,177 to 8,363 km²). Most of this new development (64%) occurred on agricultural and grasslands, whereas 33% occurred on forested land. Some smaller municipalities lost as much as 17% of their forest lands and 36% of their agricultural lands to development, although in the outlying counties losses ranged from 0% to 1.4% for forests and 0% to 2.6% for agriculture. Fast-growing urban areas surrounded by forested land experienced the most loss of forest to impervious surfaces. These estimates could be used for the monitoring of the impacts of development across the Chesapeake Bay watershed, and the approach has utility for other regions nationwide. In turn, the results and the approach can help jurisdictions set goals for resource land protection and acquisition that are consistent with regional restoration goals.     

Tools for Groundwater Protection Planning: An Example from McHenry County, Illinois, USA

/ This paper presents an approach for producing aquifer sensitivity maps from three-dimensional geologic maps, called stack-unit maps. Stack-unit maps depict the succession of geologic materials to a given depth, and aquifer sensitivity maps interpret the successions according totheir ability to transmit potential contaminants. Using McHenry County, Illinois, as a case study, stack-unit maps and an aquifer sensitivity assessment were made to help land-use planners, public health officials, consultants, developers, and the public make informed decisions regarding land use. A map of aquifer sensitivity is important for planning because the county is one of the fastest growing counties in the nation, and highly vulnerable sand and gravel aquifers occur within 6 m of ground surface over 75% of its area. The aquifer sensitivity map can provide guidance to regulators seeking optimal protection of groundwater resources where these resources are particularly vulnerable. In addition, the map can be used to help officials direct waste-disposal and industrial facilities and other sensitive land-use practices to areas where the least damage is likely to occur, thereby reducing potential future liabilities. KEY WORDS: Geologic mapping; Groundwater; Aquifers; Aquifer sensitivity; Land-use planning     

EXTRACTION AND UTILIZATION OF SPACE ACQUIRED PHYSIOGRAPHIC DATA FOR WATER RESOURCES DEVELOPMENT1

ABSTRACT: ERTS-1 satellite imagery has been evaluated as a means of providing useful watershed physiography information. From these data physiographic parameters such as drainage basin area and shape, drainage density, stream length and sinuosity, and the percentage of a watershed occupied by major land use types were obtained in three study areas. The study areas were: (1) Southwestern Wisconsin; (2) Eastern Colorado; and (3) portions of the Middle Atlantic States Using ERTS-1 imagery at 1:250,000 and 1:100,000 scales it was found that drainage basin area and shape and stream sinuosity were comparable (within 10%) in all study areas to physiographic measurements derived from conventional topographic maps at the same scales Land use information can be usefully extracted for watersheds as small as 30 mi²(78 km²) in area. Improved drainage network and density information is obtained from ERTS-1 imagery in dissected areas such as Southwestern Wisconsin, but in heavily vegetated areas (Middle Atlantic States) or areas with little physical relief (Eastern Colorado) low order streams are difficult to detect and the derived drainage densities are significantly smaller than those obtained from standard maps. It is concluded that ERTS-1 imagery can be employed to advantage in mean annual runoff prediction techniques and in providing or maintaining land use information used in the calibration and operation of watershed models.     

Usefulness of natural regions for lake management: Analysis of variation among lakes in northwestern Wisconsin,USA

A map of summer total phosphorus in lakes was compiled recently for a three-state area of the upper Midwest for purposes of identifying regional patterns of total phosphorus in lakes and attainable lake trophic state. Spatial patterns in total phosphorus from approximately 3000 lakes were studied in conjunction with maps of geographic characteristics that tend to affect phosphorus balance in lakes to identify regions of similarity in phosphorus concentrations in lakes or similarity in the mosaic of values as compared to adjacent areas. While degrees of relative homogeneity are apparent at many scales, the map was designed at a scale that would yield regions with sufficient homogeneity to be useful for lake management throughout the area. In this study, data from 210 lakes in a 1560-mi² area in northwestern Wisconsin, sampled by the Wisconsin Department of Natural Resources in the spring of 1988 (subsequent to the compilation of the phosphorus map), were examined to: (1) substantiate the existence of the regions depicted on the map in northwest Wisconsin, (2) determine the nature and relative precision of the regional boundaries, (3) determine the relative importance of natural and anthropogenic watershed characteristics, lake types, lake area, and lake depth in explaining within-region differences in lake phosphorus, and (4) demonstrate how the regions might be used by local lake managers.     

WATERSHED WEIGHTING OF EXPORT COEFFICIENTS TO MAP CRITICAL PHOSPHOROUS LOADING AREAS1

ABSTRACT: The Export Coefficient model (ECM) is capable of generating reasonable estimates of annual phosphorous loading simply from a watershed's land cover data and export coefficient values (ECVs). In its current form, the ECM assumes that ECVs are homogeneous within each land cover type, yet basic nutrient runoff and hydrological theory suggests that runoff rates have spatial patterns controlled by loading and filtering along the flow paths from the upslope contributing area and downslope dispersal area. Using a geographic information system (GIS) raster, or pixel, modeling format, these contributing area and dispersal area (CADA) controls were derived from the perspective of each individual watershed pixel to weight the otherwise homogeneous ECVs for phosphorous. Although the CADA‐ECM predicts export coefficient spatial variation for a single land use type, the lumped basin load is unaffected by weighting. After CADA weighting, a map of the new ECVs addressed the three fundamental criteria for targeting critical pollutant loading areas: (1) the presence of the pollutant, (2) the likelihood for runoff to carry the pollutant offsite, and (3) the likelihood that buffers will trap nutrients prior to their runoff into the receiving water body. These spatially distributed maps of the most important pollutant management areas were used within New York's West Branch Delaware River watershed to demonstrate how the CADA‐ECM could be applied in targeting phosphorous critical loading areas.     

RELATIONSHIP BETWEEN ANNUAL RUNOFF AND WATERSHED AREA FOR THE EASTERN UNITED STATES1

ABSTRACT: As part of the U.S. Environmental Protection Agency's effort to determine the long-term effects of acidic deposition on surface water chemistry, annual runoff was estimated for about 1000 ungaged sites in the eastern U.S. using runoff contour maps. One concern in using contour maps was that a bias may be introduced in the runoff estimates due to the size of the 1000 ungaged sites relative to the size of the watersheds used in developing the maps. To determine if a bias was present the relationship between the annual runoff (expressed as depth) and the watershed area for the Northeast (NE) and Southern Blue Ridge Province (SBRP) was tested using five regional data bases. One short-term data base (1984 Water Year, n = 531) and two long-term data bases (1940–57, n = 134 and 1951–80, n = 342) were used in the NE. In the SBRP one short-term database (1984 Water Year, n = 531) and one long-term data base (1951–80, n = 60) were used. For the NE and the SBRP, runoff was not directly correlated with watershed area using the five regional databases. Also, runoff normalized by precipitation was not related to watershed area.     

ESTIMATED RUNOFF FROM MAN-MADE SNOW1

ABSTRACT: The consumptive loss from man-made snowmaking at six Colorado ski areas is calculated. The focus of the procedures in this investigation is on the consumptive loss that occurs to man-made snow particles during the period they reside on or in the snowpack until spring snowmelt (termed the watershed loss). Calculated watershed losses under a variety of precipitation and temperature conditions at six ski areas varied from 7 to 33 percent. These calculations were made using the calibrated Subalpine Water Balance Simulation Model (Leaf and Brink, 1973a, 1973b). The watershed loss of 7 to 33 percent indicates the range of likely watershed losses that can be expected at Colorado ski areas. A previous paper by the authors (Eisel et al., 1988) provided estimates of the mean consumptive loss during the snowmaking process (termed initial loss) for conditions existing at Colorado ski areas to be 6 percent of water applied. Therefore, based on the mean initial loss, the total consumptive loss from man-made snowmaking under conditions found at Colorado ski areas could be expected to range from 13 to 37 percent. These results demonstrate the range of total consumptive losses that could be expected in various years and for various watershed conditions. These total percentage losses cannot be extrapolated directly to other specific sites because the total consumptive loss is dependent on temperature during actual snowmaking, temperature and precipitation throughout the winter at the specific ski area, and watershed conditions at the ski area. Consumptive losses to man-made snow for a specific ski area should be estimated using the handbook procedures developed especially for this purpose (Colorado Ski Country USA, 1986b).     

EFFECTS OF IMPROVED PRECIPITATION ESTIMATES ON AUTOMATED RUNOFF MAPPING: EASTERN UNITED STATES1

ABSTRACT: We evaluated maps of runoff created by means of two automated procedures. We implemented each procedure using precipitation estimates of both 5-km and 10-km resolution from PRISM (Parameter-elevation Regressions on Independent Slopes Model). Our goal was to determine if using the 5-km PRISM estimates would improve map accuracy. Visual inspection showed good general agreement among our runoff maps, as well as between our maps and one produced using a manual method. A quantitative uncertainty analysis comparing runoff interpolated from our maps with gage data that had been withheld showed slightly smaller actual and percentage interpolation errors for the 5-km PRISM-based maps. Our analyses suggest a modest region-wide improvement in runoff map accuracy with the use of PRISM-based precipitation estimates of 5-km (compared to 10-km) resolution.     

ANALYSIS AND PREDICTION OF SURFACE RUNOFF IN AN URBANIZING WATERSHED USING SATELLITE IMAGERY1

ABSTRACT: This paper demonstrates how satellite image data [e.g., from Landsat 5 Thematic Mapper (TM)], in conjunction with an urban growth model and simple runoff calculations, can be used to estimate future surface runoff and, by implication, water quality within a watershed. To illustrate the method, predictions of land use change and surface runoff are shown for Spring Creek Watershed, a medium sized urbanizing watershed in Central Pennsylvania. Land cover classifications for this watershed were created from images for summertime 1986 and 1996 and subsequently used as input to the Clarke urban growth model, called SLEUTH, to predict land use changes to the year 2025. Simulations with this model show a progressive growth in the percentage of urban pixels and in impervious surface area in the watershed but also an increase in woodland, primarily in previously clear‐cut areas. Given that woodland area will continue to increase in area, surface runoff into Spring Creek is predicted to remain only slightly above present level. However, should the woodland amount fail to increase, surface runoff is then predicted to increase more significantly during the next 25 years. Finally, the concept of urban sprawl is addressed within the context of predicted increases in urbanization by relating the implied increase in impervious surface area to population density within the watershed.     

Groundwater Vulnerability,Brownfield Redevelopment and Land Use Planning

An understanding of groundwater vulnerability in urban watersheds is important for the prevention of both surface water and groundwater contamination and can therefore be a useful tool in brownfield redevelopment and land use planning. Although industrial activity in southeastern Michigan has historically been restricted to the urbanized sections of metropolitan Detroit, new industrial development is rapidly taking place in rural and undeveloped areas. Although environmentalists and urban planners agree that industrial site recycling in urban centres (a.k.a. brownfield redevelopment) is preferable to developing green areas, many older sites remain undeveloped due to real and perceived risks. Using a PC-based geographic information system, a conceptual model of solute transport in soil was developed to evaluate potential impacts to both groundwater and surface water quality resulting from industrial development. The model was used to create a map of groundwater vulnerability within the Rouge River watershed of southeastern Michigan. The map has been used to pin-point several rural and undeveloped areas where groundwater quality is threatened by proposed development. It has also clearly demonstrated that many older brownfield sites, within the City of Detroit, are located on materials that have a much lower vulnerability to groundwater contamination and may therefore be far less costly to redevelop than greenfield sites in undeveloped areas of the watershed.     

RESERVOIR SEDIMENTATION AND ITS CAUSES IN ILLINOIS1

ABSTRACT: During the years 1930 to 1978 a research program has been carried out in Illinois dealing with reservoir sedimentation. Detailed surveys have been made on 107 lakes. A generalized graph allows annual reservoir capacity loss in percent to be estimated based on size of watershed, lake capacity, and watershed ratio. For the 258 square miles of land on the watershed of Lake Springfield, complete soil maps were measured using a graphic digitizer. A computer readable map was produced. The Universal Soil Loss Equation was solved to provide upland soil losses for each acre of the watershed. Average watershed soil loss was 3.96 tons per acre per year, and 24 percent of this was delivered to the lake.     

A COMPARISON OF HAND- AND SPLINE-DRAWN PRECIPITATION MAPS FOR MOUNTAINOUS MONTANA1

ABSTRACT: Average annual precipitation for the period 1961–1990 was estimated for a mountainous region in Montana with a Laplacian thin-plate spline (ANUSPLIN) and compared to a hand-drawn map. Input data included latitude, longitude, and elevation from a three-arc-second U.S. Geological Survey Digital Elevation Model of the Bozeman and Billings 1 × 2 topographic quadrangles and precipitation data at 96 stations. The two maps are similar in appearance. Digital comparison of the two maps with ARC/INFO's Grid tools shows that mean annual precipitation for the hand-contoured map is 22.9 inches and for the ANUSPLIN map is 23.7 inches. Of the 5,760,000 cells, 53 percent showed no difference between ANUSPLIN and hand-drawn maps; 19 percent showed a two-inch difference, and 28 percent showed more than 2 inches difference. Input data and model output at the same location are not different (standard deviation 1.77, p-value 0.76). Hand-drawn maps show two inches more precipitation during the 1961–90 period than during the 1941–1970 period. Similarly, measured data at 73 sites for the period 1961–1990 are on average 2.4 inches higher than the same stations during the 1941–1970 period. The difference is significant (p-value > 0.0001).