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.     

Effects of Impervious Cover at Multiple Spatial Scales on Coastal Watershed Streams1

Abstract: The spatial scale and location of land whose development has the strongest influence on aquatic ecosystems must be known to support land use decisions that protect water resources in urbanizing watersheds. We explored impacts of urbanization on streams in the West River watershed, New Haven, Connecticut, to identify the spatial scale of watershed imperviousness that was most strongly related to water chemistry, macroinvertebrates, and physical habitat. A multiparameter water quality index was used to characterize regional urban nonpoint source pollution levels. We identified a critical level of 5% impervious cover, above which stream health declined. Conditions declined with increasing imperviousness and leveled off in a constant state of impairment at 10%. Instream variables were most correlated (0.77 ≤ |r| ≤ 0.92, p < 0.0125) to total impervious area (TIA) in the 100‐m buffer of local contributing areas (～5‐km² drainage area immediately upstream of each study site). Water and habitat quality had a relatively consistent strong relationship with TIA across each of the spatial scales of investigation, whereas macroinvertebrate metrics produced noticeably weaker relationships at the larger scales. Our findings illustrate the need for multiscale watershed management of aquatic ecosystems in small streams flowing through the spatial hierarchies that comprise watersheds with forest‐urban land use gradients.     

Bison and cattle grazing management,bare ground coverage,and links to suspended sediment concentrations in grassland streams

This study quantified the impact of bison and cattle grazing management practices on bare ground coverage at the watershed, riparian, and forested riparian scales within the Flint Hills ecoregion in Kansas. We tested for correlations between bare ground coverage and fluvial suspended sediment concentrations during base‐flow and storm‐flow events. We used remotely sensed imagery combined with field surveys to classify ground cover and quantify the presence of bare ground. Base‐flow water samples were collected bi‐monthly during rain‐free periods and 24 h following precipitation events. Storm‐flow water samples were collected on the rising limb of the hydrograph, using single‐stage automatic samplers. Ungrazed treatments contained the lowest coverage of bare ground at the watershed, riparian, and forested riparian scales. Bison treatments contained the highest coverage of bare ground at the watershed scale, while high‐density cattle treatments contained the highest coverage of bare ground at the riparian and forested riparian scales. In bison and cattle‐grazed treatments, a majority of bare ground was located near fence lines, watershed boundaries, and third‐ and fourth‐order stream segments. Inorganic sediment concentrations at base flow were best predicted by riparian bare ground coverage, while storm‐flow sediment concentrations were best predicted by watershed scale bare ground coverage.     

Estimating the Variability and Confidence of Land Use and Imperviousness Relationships at a Regional Scale1

Abstract: Impervious cover is a commonly used metric to help explain or predict anthropogenic impacts on aquatic resources; often it is used as a surrogate for intensity of human impacts when evaluating effects on aquatic resources. The most common way to estimate imperviousness is based on relationships with land use. Few studies have evaluated how the relationship between impervious surface and land use varies among geographies with different levels of development and between types of imagery used to assign land use type. In this study, we assess variability in estimates of imperviousness based on two locally available land use datasets: one based on aerial imagery (2‐m resolution) and another based on satellite imagery (30‐m resolution). The ranges and variability in imperviousness within land use categories were assessed at several spatial scales, including within counties, between counties, and between watersheds. Results indicate that there was considerable variability for all developed land use types. Estimated impervious cover often varied over a range of 20‐40% points within a land use category. Furthermore, there were clear spatial patterns both between and within counties, with impervious cover for a given land use type being higher near the urban centers and lower at the margins of development. Estimates of imperviousness for 12 study watersheds indicated that variability increased with increasing watershed development, making it difficult to confidently set management or regulatory targets based on impervious cover. This study suggests that locally derived, high resolution satellite or aerial imagery should be used to estimate imperviousness when a high level of accuracy and precision is required for regulatory or management decisions. Furthermore, the error associated with impervious land use relationships should be accounted for when using impervious cover in runoff or water quality models, or when making management decisions regarding stream health.     

Discussion

Epps et al. (2013) derived Curve Number (CN) values for two forested headwater watersheds in the Lower Coastal Plain of South Carolina during the 2008‐2011 period from rainfall‐runoff data, resulting in 23 events for the Upper Debidue Creek (UDC) watershed and in 20 events for Watershed 80 (WS80). D'Asaro and Grillone analyzed the P, CN data of the UDC watershed finding an evident “complacent” behavior, characterized by a declining CN with increasing P but without approaching a stable value at large storms. In this case, the traditional runoff CN equation does not fit well with the rainfall‐runoff data that indicate a partial source area watershed behavior and are more aptly modeled by the equation introduced by D'Asaro and Grillone (2012), who introduced a C parameter in the well‐known runoff CN formula. The C value, that represents the source area (fraction of drainage area) of the basin that produces runoff with a fixed CN < 100, is here considered to vary with the water table elevation (WTE), to which it is well correlated.     

A NEW FLASHINESS INDEX: CHARACTERISTICS AND APPLICATIONS TO MIDWESTERN RIVERS AND STREAMS1

ABSTRACT: The term flashiness reflects the frequency and rapidity of short term changes in streamflow, especially during runoff events. Flashiness is an important component of a stream's hydrologic regime. A variety of land use and land management changes may lead to increased or decreased flashiness, often to the detriment of aquatic life. This paper presents a newly developed flashiness index, which is based on mean daily flows. The index is calculated by dividing the pathlength of flow oscillations for a time interval (i.e., the sum of the absolute values of day‐to‐day changes in mean daily flow) by total discharge during that time interval. This index has low interannual variability, relative to most flow regime indicators, and thus greater power to detect trends. Index values were calculated for 515 Midwestern streams for the 27‐year period from 1975 through 2001. Statistically significant increases were present in 22 percent of the streams, primarily in the eastern portion of the study area, while decreases were present in 9 percent, primarily in the western portion. Index values tend to decrease with increasing watershed area and with increasing unit area ground water inputs. Area compensated index values often shift at ecoregion boundaries. Potential index applications include evaluation of programs to restore more natural flow regimes.     

EVALUATION OF INFLOW TO MIRROR LAKE,NEW HAMPSHIRE1

ABSTRACT: Measured stream discharge plus calculated ground water discharge (total measured runoff) were compared with runoff calculated by the unit-runoff method for the two largest watersheds of Mirror Lake for 1981–1983. Runoff calculated by the unit-runoff method, using Hubbard Brook watershed 3 as the index watershed, was greater than the total measured runoff into Mirror Lake during periods of high flow and slightly less than the total measured runoff into Mirror Lake during periods of low flow. Annual calculated unit runoff was 17 to 37 percent greater than total measured runoff. Differences in monthly runoff are far greater, ranging from 0 to greater than 100 percent. For high flows the calculated unit runoff is about 2 times greater than total measured runoff. For low flows the northwest basin of Mirror Lake has the greatest ground water contribution compared to the other two basins. In contrast, Hubbard Brook watershed 3 has the least ground water contribution.     

Thermal Characteristics of Stormwater Runoff from Asphalt and Sod Surfaces1

Abstract: Urban impervious surfaces absorb and store thermal energy, particularly during warm summer months. During a rainfall/runoff event, thermal energy is transferred from the impervious surface to the runoff, causing it to become warmer. As this higher temperature runoff enters receiving waters, it can be harmful to coldwater habitat. In an urban watershed, impervious asphalt surfaces (roads, parking lots, and driveways) and pervious residential lawns comprise a significant portion of the watershed area. A paired asphalt‐turfgrass sod plot was constructed to compare the thermal runoff characteristics between asphalt and turfgrass sod surfaces, to identify meteorological variables that influence these thermal characteristics, and to evaluate evaporative heat loss for runoff from asphalt surfaces. Rainfall simulations were conducted during the summers of 2004 and 2005 under a range of climatic conditions. Asphalt surface temperatures immediately prior to rainfall simulations averaged 43.6°C and decreased an average of 12.3°C over 60 min as rain cooled the surface. In contrast, presimulation sod surface temperatures averaged only 23.3°C and increased an average of 1.3°C throughout the rainfall events. Heat transferred from the asphalt to the runoff resulted in initial asphalt runoff temperatures averaging 35.0°C that decreased by an average of 4.1°C at the end of the event. Sod runoff temperatures averaged only 25.5°C and remained fairly constant throughout the simulations. Multivariable regression equations were developed to predict (1) average asphalt surface temperature (R² = 0.90) and average asphalt runoff temperature (R² = 0.92) as a function of solar radiation, rain temperature, and wind speed, and (2) average sod surface temperature (R² = 0.85) and average sod runoff temperature (R² = 0.94) as a function of solar radiation, rain temperature, rain intensity, and wind speed. Based on a heat balance analysis, existing evaporation equations developed from studies on lakes were not adequate to predict evaporation from runoff on a heated impervious surface. The combined heat from the asphalt and sod plots was an average of 38% less than the total heat had the total area consisted solely of asphalt.     

Nutrient and Sediment Production,Watershed Characteristics,and Land Use in the Tahoe Basin,California‐Nevada1

Abstract: In efforts to control the degradation of water quality in Lake Tahoe, public agencies have monitored surface water discharge and concentrations of nitrogen, phosphorus, and suspended sediment in two separate sampling programs. The first program focuses on 20 watersheds varying in size from 162 to 14,000 ha, with continuous stream gaging and periodic sampling; the second focuses on small urbanized catchments, with automated sampling during runoff events. Using data from both programs, we addressed the questions (1) what are the fluxes and concentrations of nitrogen and phosphorus entering the lake from surface runoff; (2) how do the fluxes and concentrations vary in space and time; and (3) how are they related to land use and watershed characteristics? To answer these questions, we calculated discharge‐weighted average concentrations and annual fluxes and used multiple regression to relate those variable to a suite of GIS‐derived explanatory variables. The final selected regression models explain 47‐62% of the variance in constituent concentrations in the stormwater monitoring catchments, and 45‐72% of the variance in mean annual yields in the larger watersheds. The results emphasize the importance of impervious surface and residential density as factors in water quality degradation, and well‐developed soil as a factor in water quality maintenance.     

Modeling Watershed‐Wide Bioretention Stormwater Retrofits to Achieve Thermal Pollution Mitigation Goals

Stream ecosystems are increasingly at risk for thermal impairment as urbanization intensifies, resulting in more heated runoff from impervious cover that is less likely to be cooled naturally. While several best management practices, including bioretention filters, have been able to reduce thermal pollution, success has been limited. The extent of thermal mitigation required to prevent ecological damage remains unknown. A calibrated runoff temperature model of a case study watershed in Blacksburg, VA was developed to determine the cumulative treatment volume of bioretention filters required to reduce thermal impacts caused by runoff from development in the watershed to regulated biologically acceptable levels. A future build out scenario of the study watershed was also analyzed. Results from this study established that runoff thermal pollution cannot be fully reduced to goal thresholds during all storms using bioretention filter retrofits. While retrofitting significantly decreased temperatures and heat exports relative to the controls, increasing treatment volumes did not really enhance mitigation. Alternate thermal mitigation methods that actively remove runoff volume should be considered where more thermal mitigation is required.     

A COMPARISON OF PHOSPHORUS AND WATER CONTRIBUTIONS BY SNOWFALL AND RAIN IN NORTHEAST OHIO1

ABSTRACT A record snowfall of 55.8 centimeters occurred on December 1 and 2, 1974 in Portage County, Ohio. An early winter thaw melted the greater part of the snow by December 22, 1974, and a two-day rain fell from December 23 to December 25. These weather events provided an opportunity to compare snowmelt and rainfall contribution to runoff and phosphorus loading to the Twin Lakes Watershed. Phosphorus concentrations of the snow and rain were determined. Six lake inflows and two lake outflows were measured daily for volume and phosphorus concentration. The snow added 217,000 cubic meters of water and 2.2 kilograms of total phosphorus to the watershed. The rain added 74,000 cubic meters of water and 1.6 kilograms of total phosphorus. Total water discharge from the watershed during December was 244,537 cubic meters and total phosphorus output was 20.3 kilograms. The snow provided 49.9% of the discharge and 8% of the phosphorus whereas the rainfall contributed 28% of the discharge and 6% of the phosphorus. These results indicate that while snow is a significant source of water, it is not a large source of phosphorus. The greatest contribution of phosphorus comes from fine sediment carried by storm runoff.     

The Influence of Forest Regrowth on the Stream Discharge in the North Carolina Piedmont Watersheds

This study focuses on the relationships of watershed runoff with historical land use/land cover (LULC) and climate trends. Over the 20th Century, LULC in the Southeast United States, particularly the North Carolina Piedmont, has evolved from an agriculture dominated to an extensively forested landscape with more recent localized urbanization. The regrowth of forest has an important influence on the hydrology of the region as it enhances ecosystem interaction with recent climate change. During 1920‐2009, the amount of precipitation in some parts of the North Carolina Piedmont forest regrowth area showed increasing trends without corresponding increments in runoff. We employed the Soil and Water Assessment Tool (SWAT) to backcast long‐term hydrologic behavior of watersheds in North Carolina with different LULC conditions: (1) LULC conversion from agricultural to forested area and (2) long‐term stable forested area. Comparing U.S. Geological Survey‐measured stream discharge with SWAT‐simulated stream discharge under the assumption of constant 2006 LULC, we found significant stream discharge underprediction by SWAT in two LULC conversion watersheds during the early simulation period (1920s) with differences gradually decreasing by the mid‐1970s. This model bias suggests that forest regrowth on abandoned agricultural land was a key factor contributing to mitigate the impact of increased precipitation on runoff due to increasing water consumption driven by changes in vegetation.     

Physically Based,Hydrologic Model Results Based on Three Precipitation Products1

Abstract: The main objective of the study is to examine the accuracy of and differences among simulated streamflows driven by rainfall estimates from a network of 22 rain gauges spread over a 2,170 km² watershed, NEXRAD Stage III radar data, and Tropical Rainfall Measuring Mission (TRMM) 3B42 satellite data. The Gridded Surface Subsurface Hydrologic Analysis (GSSHA), a physically based, distributed parameter, grid‐structured, hydrologic model, was used to simulate the June‐2002 flooding event in the Upper Guadalupe River watershed in south central Texas. There were significant differences between the rainfall fields estimated by the three types of measurement technologies. These differences resulted in even larger differences in the simulated hydrologic response of the watershed. In general, simulations driven by radar rainfall yielded better results than those driven by satellite or rain‐gauge estimates. This study also presents an overview of effects of land cover changes on runoff and stream discharge. The results demonstrate that, for major rainfall events similar to the 2002 event, the effect of urbanization on the watershed in the past two decades would not have made any significant effect on the hydrologic response. The effect of urbanization on the hydrologic response increases as the size of the rainfall event decreases.     

Temporal Trends in the Spatial Distribution of Impervious Cover Relative to Stream Location

Use of impervious cover is transitioning from an indicator of surface water condition to one that also guides and informs watershed planning and management, including Clean Water Act (33 U.S.C. §1251 et seq.) reporting. Whether it is for understanding surface water condition or planning and management, impervious cover is most commonly expressed as summary measurement (e.g., percentage watershed in impervious cover). We use the National Land Cover Database to estimate impervious cover in the vicinity of surface waters for three time periods (2001, 2006, 2011). We also compare impervious cover in the vicinity of surface waters to watershed summary estimates of impervious cover for classifying the spatial pattern of impervious cover. Between 2001 and 2011, surface water shorelines (streams and water bodies) in the vicinity of impervious cover increased nearly 10,000 km. Across all time periods, approximately 27% of the watersheds in the continental United States had proximally distributed impervious cover, i.e., the percentage of impervious cover in the vicinity of surface waters was higher than its watershed summary expression. We discuss how impervious cover spatial pattern can be used to inform watershed planning and management, including reporting under the Clean Water Act.     

EFFECTS OF MULTI‐SCALE ENVIRONMENTAL CHARACTERISTICS ON AGRICULTURAL STREAM BIOTA IN EASTERN WISCONSIN1

ABSTRACT: The U.S. Geological Survey examined 25 agricultural streams in eastern Wisconsin the determine relations between fish, invertebrate, and algal metrics and multiple spatial scales of land cover, geologic setting, hydrologic, aquatic habitat, and water chemistry data. Spearman correlation and redundancy analyses were used to examine relations among biotic metrics and environmental characteristics. Riparian vegetation, geologic, and hydrologic conditions affected the response of biotic metrics to watershed agricultural land cover but the relations were aquatic assemblage dependent. It was difficult to separate the interrelated effects of geologic setting, watershed and buffer land cover, and base flow. Watershed and buffer land cover, geologic setting, reach riparian vegetation width, and stream size affected the fish IBI, invertebrate diversity, diatom IBI, and number of algal taxa; however, the invertebrate FBI, percentage of EPT, and the diatom pollution index were more influenced by nutrient concentrations and flow variability. Fish IBI scores seemed most sensitive to land cover in the entire stream network buffer, more so than watershed‐scale land cover and segment or reach riparian vegetation width. All but one stream with more than approximately 10 percent buffer agriculture had fish IBI scores of fair or poor. In general, the invertebrate and algal metrics used in this study were not as sensitive to land cover effects as fish metrics. Some of the reach‐scale characteristics, such as width/depth ratios, velocity, and bank stability, could be related to watershed influences of both land cover and geologic setting. The Wisconsin habitat index was related to watershed geologic setting, watershed and buffer land cover, riparian vegetation width, and base flow, and appeared to be a good indicator of stream quality Results from this study emphasize the value of using more than one or two biotic metrics to assess water quality and the importance of environmental characteristics at multiple scales.     

A Modeling System to Assess Land Cover Land Use Change Effects on SAV Habitat in the Mobile Bay Estuary

Estuarine ecosystems are largely influenced by watersheds directly connected to them. In the Mobile Bay, Alabama watersheds we examined the effect of land cover and land use (LCLU) changes on discharge rate, water properties, and submerged aquatic vegetation, including freshwater macrophytes and seagrasses, throughout the estuary. LCLU scenarios from 1948, 1992, 2001, and 2030 were used to influence watershed and hydrodynamic models and evaluate the impact of LCLU change on shallow aquatic ecosystems. Overall, our modeling results found that LCLU changes increased freshwater flows into Mobile Bay altering temperature, salinity, and total suspended sediments (TSS). Increased urban land uses coupled with decreased agricultural/pasture lands reduced TSS in the water column. However, increased urbanization or agricultural/pasture land coupled with decreased forest land resulted in higher TSS concentrations. Higher sediment loads were usually strongly correlated with higher TSS levels, except in areas where a large extent of wetlands retained sediment discharged during rainfall events. The modeling results indicated improved water clarity in the shallow aquatic regions of Mississippi Sound and degraded water clarity in the Wolf Bay estuary. This integrated modeling approach will provide new knowledge and tools for coastal resource managers to manage shallow aquatic habitats that provide critical ecosystem services.     

Modeling the Effect of Summertime Heating on Urban Runoff Temperature1

Abstract: The summertime heating of runoff in urban areas is recognized as a common and consistent urban climatological phenomenon. In this study, a simple thermal urban runoff model (TURM) is presented for the net energy flux at the impervious surfaces of urban areas to account for the heat transferred to runoff. The first step in developing TURM consists of calculating the various factors that control how urban impervious areas absorb heat and transfer it to moving water on the surface. The runoff temperature is determined based on the interactions of the physical characteristics of the impervious areas, the weather, and the heat transfer between the moving film of runoff and the impervious surface common in urban areas. Key surface and weather factors that affect runoff temperature predictions are type of impervious surface, air temperature, humidity, solar radiation before and during rain, rainfall intensity, and rainfall temperature. Runoff from pervious areas is considered separately and estimated using the Green‐Ampt Mein‐Larson rainfall excess method. Pervious runoff temperature is estimated as the rainfall temperature. Field measurements indicate that wet bulb temperature can be used as a surrogate for rainfall temperature and that runoff temperatures from sod average just 2°C higher than rainfall temperatures. Differences between measured and predicted impervious runoff temperature average approximately 2°C, indicating that TURM is a useful tool for determining runoff temperatures for typical urban areas.     

EFFECTS OF URBANIZATION ON CHANNEL INSTABILITY1

ABSTRACT: Channel instability and aquatic ecosystem degradation have been linked to watershed imperviousness in humid regions of the U.S. In an effort to provide a more process‐based linkage between observed thresholds of aquatic ecosystem degradation and urbanization, standard single event approaches (U.S. Geological Survey Flood Regression Equations and rational) and continuous hydrologic models (HSPF and CASC2D) were used to examine potential changes in flow regime associated with varying levels of watershed imperviousness. The predicted changes in flow parameters were then interpreted in concert with risk‐based models of channel form and instability. Although low levels of imperviousness (10 to 20 percent) clearly have the potential to destabilize streams, changes in discharge, and thus stream power, associated with increased impervious area are highly variable and dependent upon watershed‐specific conditions. In addition to the storage characteristics of the pre‐development watershed, the magnitude of change is sensitive to the connectivity and conveyance of impervious areas as well as the specific characteristics of the receiving channels. Different stream types are likely to exhibit varying degrees and types of instability, depending on entrenchment, relative erodibility of bed and banks, riparian condition, mode of sediment transport (bedload versus suspended load), and proximity to geomorphic thresholds. Nonetheless, simple risk‐based analyses of the potential impacts of land use change on aquatic ecosystems have the potential to redirect and improve the effectiveness of watershed management strategies by facilitating the identification of channels that may be most sensitive to changes in stream power.     

Curve Number Derivation for Watersheds Draining Two Headwater Streams in Lower Coastal Plain South Carolina,USA

The objective of this study was to assess curve number (CN) values derived for two forested headwater catchments in the Lower Coastal Plain (LCP) of South Carolina using a three‐year period of storm event rainfall and runoff data in comparison with results obtained from CN method calculations. Derived CNs from rainfall/runoff pairs ranged from 46 to 90 for the Upper Debidue Creek (UDC) watershed and from 42 to 89 for the Watershed 80 (WS80). However, runoff generation from storm events was strongly related to water table elevation, where seasonally variable evapotranspirative wet and dry moisture conditions persist. Seasonal water table fluctuation is independent of, but can be compounded by, wet conditions that occur as a result of prior storm events, further complicating flow prediction. Runoff predictions for LCP first‐order watersheds do not compare closely to measured flow under the average moisture condition normally associated with the CN method. In this study, however, results show improvement in flow predictions using CNs adjusted for antecedent runoff conditions and based on water table position. These results indicate that adaptations of CN model parameters are required for reliable flow predictions for these LCP catchments with shallow water tables. Low gradient topography and shallow water table characteristics of LCP watersheds allow for unique hydrologic conditions that must be assessed and managed differently than higher gradient watersheds.     

BASE FLOW TRENDS IN URBANIZING WATERSHEDS OF THE DELAWARE RIVER BASIN1

Rapid land development is raising concern regarding the ability of urbanizing watersheds to sustain adequate base flow during periods of drought. Long term streamflow records from unregulated watersheds of the lower to middle Delaware River basin are examined to evaluate the impact of urbanization and imperviousness on base flow. Trends in annual base flow volumes, seven‐day low flows, and runoff ratios are determined for six urbanizing watersheds and four reference watersheds across three distinct physiographic regions. Hydrograph separation is used to determine annual base flow and stormflow volumes, and nonparametric trend tests are conducted on the resulting time series. Of the watersheds examined, the expected effects of declining base flow volumes and seven‐day low flows and increasing stormflows are seen in only one watershed that is approximately 20 percent impervious and has been subject to a net water export over the past 15 years. Both interbasin transfers and hydrologic mechanisms are invoked to explain these results. The results show that increases in impervious area may not result in measurable reductions in base flow at the watershed scale.