Heterogeneity of Hydrologic Response in Urban Watersheds1

Meierdiercks, Katherine L., James A. Smith, Mary Lynn Baeck, and Andrew J. Miller, 2010. Heterogeneity of Hydrologic Response in Urban Watersheds. Journal of the American Water Resources Association (JAWRA) 46(6):1221–1237. DOI: 10.1111/j.1752-1688.2010.00487.x Abstract: The changing patterns of streamflow associated with urbanization are examined through analyses of discharge and rainfall records for the study watersheds of the Baltimore Ecosystem Study (BES). Analyses utilize a decade (1999-2008) of observations from a dense U.S. Geological Survey stream gaging network and Hydro-NEXRAD radar rainfall fields. The principal study watershed of the BES is Gwynns Falls (171 km²). Focus is given to two Gwynns Falls basins with contrasting patterns and histories of development, Dead Run and Upper Gwynns Falls. The sharp contrasts in streamflow properties between the basins reflect the differences in urban development prior to implementation of stormwater management regulations (much of Dead Run) and development for which stormwater management is an integral part of the hydrologic system (Upper Gwynns Falls). The mean annual runoff in Dead Run (558 mm) is 35% larger than that of Upper Gwynns Falls; larger contrasts in runoff properties typify the “warm season” and are linked to storm event hydrologic response. Spatial heterogeneities in storm event response are reflected in seasonal and diurnal properties of streamflow. Analyses of storm event response are presented for June 2006, during which monthly rainfall over the BES region ranged from less than 150 to more than 500 mm. Baisman Run, the BES forest reference watershed, and Moores Run, a highly urbanized watershed in Baltimore City, provide “end-member” representations of urban impacts on streamflow.     

PENN STATE URBAN RUNOFF MODEL TO PINPOINT FLOOD PEAK SOURCE LOCATION1

ABSTRACT: The Penn State Urban Runoff Model, developed in 1976, is described in this paper. Aside from locating infiltration and detention basin operation in an unconventional manner, the model includes a peak flow presentation table which identifies watershed subareas chiefly responsible for the occurrence of flooding conditions at certain points in the watershed. The results of a case study on an urban drainage basin in the Philadelphia area is discussed, and preferred sites for retention ponds are suggested. The simplicity of the Penn State model is pointed out and computer run costs between 10 and 20 percent of the corresponding cards for HEC-I and SWMM are cited.     

Modeling Hydrologic Benefits of Low Impact Development: A Distributed Hydrologic Model of The Woodlands,Texas

Low Impact Development (LID) is alternative design approach to land development that conserves and utilizes natural resources to minimize the potential negative environmental impacts of development, such as flooding. The Woodlands near Houston, Texas is one of the premier master‐planned communities in the United States. Unlike in a typical urban development where riparian corridors are often replaced with concrete channels, pervious surfaces, vegetation, and natural drainage pathways were preserved as much as possible during development. In addition, a number of detention ponds were strategically located to manage runoff on site. This article uses a unique distributed hydrologic model, Vflo?, combined with historical (1974) and recent (2008 and 2009) rainfall events to evaluate the long‐term effectiveness of The Woodlands natural drainage design as a stormwater management technique. This study analyzed the influence of LID in The Woodlands by comparing the hydrologic response of the watershed under undeveloped, developed, and highly urbanized conditions. The results show that The Woodlands drainage design successfully reflects predeveloped hydrologic conditions and produces peak flows two to three times lower than highly urbanized development. Furthermore, results indicate that the LID practices employed in The Woodlands successfully attenuate the peak flow from a 100‐year design event, resulting in flows comparable to undeveloped hydrologic conditions.     

Identifying Urban Features from LiDAR for a High‐Resolution Urban Hydrologic Model

Light Detection and Ranging (LiDAR), is relatively inexpensive, provides high spatial resolution sampling at great accuracy, and can be used to generate surface terrain and land cover datasets for urban areas. These datasets are used to develop high‐resolution hydrologic models necessary to resolve complex drainage networks in urban areas. This work develops a five‐step algorithm to generate indicator fields for tree canopies, buildings, and artificial structures using Geographic Resources Analysis Support System (GRASS‐GIS), and a common computing language, Matrix Laboratory. The 54 km² study area in Parker, Colorado consists of twenty‐four 1,500 × 1,500 m LiDAR subsets at 1 m resolution with varying degrees of urbanization. The algorithm correctly identifies 96% of the artificial structures within the study area; however, application success is dependent upon urban extent. Urban land use fractions below 0.2 experienced an increase in falsely identified building locations. ParFlow, a three‐dimensional, grid‐based hydrological model, uses these building and artificial structure indicator fields and digital elevation model for a hydrologic simulation. The simulation successfully develops the complex drainage network and simulates overland flow on the impervious surfaces (i.e., along the gutters and off rooftops) made possible through this spatial analysis process.     

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.     

Evaluation of the National Land Cover Database for Hydrologic Applications in Urban and Suburban Baltimore,Maryland1

Smith, Monica Lipscomb, Weiqi Zhou, Mary Cadenasso, Morgan Grove, and Lawrence E. Band, 2010. Evaluation of the National Land Cover Database for Hydrologic Applications in Urban and Suburban Baltimore, Maryland. Journal of the American Water Resources Association (JAWRA) 46(2):429-442. DOI: 10.1111/j.1752-1688.2009.00412.x Abstract: We compared the National Land Cover Database (NLCD) 2001 land cover, impervious, and canopy data products to land cover data derived from 0.6-m resolution three-band digital imagery and ancillary data. We conducted this comparison at the 1 km², 9 km², and gauged watershed scales within the Baltimore Ecosystem Study to determine the usefulness and limitations of the NLCD in heterogeneous urban to exurban environments for the determination of land-cover information for hydrological applications. Although the NLCD canopy and impervious data are significantly correlated with the high-resolution land-cover dataset, both layers exhibit bias at <10 and >70% cover. The ratio of total impervious area and connected impervious area differs along the range of percent imperviousness – at low percent imperviousness, the NLCD is a better predictor of pavement alone, whereas at higher percent imperviousness, buildings and pavement together more resemble NLCD impervious estimates. The land-cover composition and range for each NLCD urban land category (developed open space, low-intensity, medium-intensity, and high-intensity developed) is more variable in areas of low-intensity development. Fine-vegetation land-cover/lawn area is incorporated in a large number of land use categories with no ability to extract this land cover from the NLCD. These findings reveal that the NLCD may yield important biases in urban, suburban, and exurban hydrologic analyses where land cover is characterized by fine-scale spatial heterogeneity.     

Assessing the Hydrologic and Water Quality Benefits of a Network of Stormwater Control Measures in a SE U.S. Piedmont Watershed

The hydrologic and water quality benefits of an existing engineered stormwater control measures (SCMs) network, along with the alternative stormwater control simulations, were assessed in the rapidly urbanizing Beaverdam Creek watershed located in SE U.S. Piedmont region through the use of distributed Model of Urban Stormwater Improvement Conceptualization stormwater model. When compared with predevelopment conditions, the postdevelopment watershed simulation without SCMs indicated a 2 times increase in total runoff volume, 3 times average increase in peak flow for 1.5‐3.2 cm 6‐h storm events, and 30 times, 12 times, and 3 times higher total suspended solids (TSS), total phosphorous (TP), and total nitrogen (TN) loadings, respectively. The existing SCMs network, in comparison with the postdeveloped watershed without SCMs, reduced the average peak flow rates for 1.5‐3.2 cm 6‐h storm events by 70%, lowered the annual runoff volume by 3%, and lowered TSS, TP, TN annual loads by 57, 51, and 10%, respectively. A backyard rain garden simulation resulted in minimal additional reduction in TSS (1.6%), TP (0.4%), and TN (4%). Model simulations indicate that mandatory 85% TSS and 70% TP annual load reductions in comparison with the predevelopment levels would require the diversion of runoff from at least 70% of the contributing drainage areas runoff into additional offline bioretention basins.     

Impediments and Solutions to Sustainable,Watershed-Scale Urban Stormwater Management: Lessons from Australia and the United States

In urban and suburban areas, stormwater runoff is a primary stressor on surface waters. Conventional urban stormwater drainage systems often route runoff directly to streams and rivers, thus exacerbating pollutant inputs and hydrologic disturbance, and resulting in the degradation of ecosystem structure and function. Decentralized stormwater management tools, such as low impact development (LID) or water sensitive urban design (WSUD), may offer a more sustainable solution to stormwater management if implemented at a watershed scale. These tools are designed to pond, infiltrate, and harvest water at the source, encouraging evaporation, evapotranspiration, groundwater recharge, and re-use of stormwater. While there are numerous demonstrations of WSUD practices, there are few examples of widespread implementation at a watershed scale with the explicit objective of protecting or restoring a receiving stream. This article identifies seven major impediments to sustainable urban stormwater management: (1) uncertainties in performance and cost, (2) insufficient engineering standards and guidelines, (3) fragmented responsibilities, (4) lack of institutional capacity, (5) lack of legislative mandate, (6) lack of funding and effective market incentives, and (7) resistance to change. By comparing experiences from Australia and the United States, two developed countries with existing conventional stormwater infrastructure and escalating stream ecosystem degradation, we highlight challenges facing sustainable urban stormwater management and offer several examples of successful, regional WSUD implementation. We conclude by identifying solutions to each of the seven impediments that, when employed separately or in combination, should encourage widespread implementation of WSUD with watershed-based goals to protect human health and safety, and stream ecosystems.     

Modeling Nitrate-Nitrogen Load Reduction Strategies for the Des Moines River,Iowa Using SWAT

The Des Moines River that drains a watershed of 16,175 km² in portions of Iowa and Minnesota is impaired for nitrate-nitrogen (nitrate) due to concentrations that exceed regulatory limits for public water supplies. The Soil Water Assessment Tool (SWAT) model was used to model streamflow and nitrate loads and evaluate a suite of basin-wide changes and targeting configurations to potentially reduce nitrate loads in the river. The SWAT model comprised 173 subbasins and 2,516 hydrologic response units and included point and nonpoint nitrogen sources. The model was calibrated for an 11-year period and three basin-wide and four targeting strategies were evaluated. Results indicated that nonpoint sources accounted for 95% of the total nitrate export. Reduction in fertilizer applications from 170 to 50 kg/ha achieved the 38% reduction in nitrate loads, exceeding the 34% reduction required. In terms of targeting, the most efficient load reductions occurred when fertilizer applications were reduced in subbasins nearest the watershed outlet. The greatest load reduction for the area of land treated was associated with reducing loads from 55 subbasins with the highest nitrate loads, achieving a 14% reduction in nitrate loads achieved by reducing applications on 30% of the land area. SWAT model results provide much needed guidance on how to begin implementing load reduction strategies most efficiently in the Des Moines River watershed.     

NATURAL RESTABILIZATION OF STREAM CHANNELS IN URBAN WATERSHEDS1

ABSTRACT: Stream channels are known to change their form as a result of watershed urbanization, but do they restabilize under subsequent conditions of constant urban land use? Streams in seven developed and developing watersheds (drainage areas 5–35 km²) in the Puget Sound lowlands were evaluated for their channel stability and degree of urbanization, using field and historical data. Protocols for determining channel stability by visual assessment, calculated bed mobility at bankfull flows, and resurveyed cross‐sections were compared and yielded nearly identical results. We found that channel restabilization generally does occur within one or two decades of constant watershed land use, but it is not universal. When (or if) an individual stream will restabilize depends on specific hydrologic and geomorphic characteristics of the channel and its watershed; observed stability is not well predicted by simply the magnitude of urban development or the rate of ongoing land‐use change. The tendency for channel restabilization suggests that management efforts focused primarily on maintaining stability, particularly in a still‐urbanizing watershed, may not always be necessary. Yet physical stability alone is not a sufficient condition for a biologically healthy stream, and additional rehabilitation measures will almost certainly be required to restore biological conditions in urban systems.     

Testing of the Storm Water Management Model Low Impact Development Modules

Stormwater infrastructure designers and operators rely heavily on the United States Environmental Protection Agency’s Storm Water Management Model (SWMM) to simulate stormwater and wastewater infrastructure performance. Since its inception in the late 1970s, improvements and extensions have been tested and evaluated rigorously to verify the accuracy of the model. As a continuation of this progress, the main objective of this study was to quantify how accurately SWMM simulates the hydrologic activity of low impact development (LID) storm control measures. Model performance was evaluated by quantitatively comparing empirical data to model results using a multievent, multiobjective calibration method. The calibration methodology utilized the PEST software, a Parameter ESTimation tool, to determine unmeasured hydrologic parameters for SWMM’s LID modules. The calibrated LID modules’ Nash–Sutcliffe efficiencies averaged 0.81; average percent bias (PBIAS) ?9%; average ratio of root mean square error to standard deviation of measured values 0.485; average index of agreement 0.94; and the average volume error, simulated vs. observed, was +9%. SWMM accurately predicted the timing of peak flows, but usually underestimated their magnitudes by 10%. The average volume reduction, measured outflow volume divided by inflow volume, was 48%. We had more difficulty in calibrating one study, an infiltration trench, which identified a significant limitation of the current version of the SWMM LID module; it cannot simulate lateral exfiltration of water out of the storage layers of a LID storm control measure. This limitation is especially severe for a deep LIDs, such as infiltration trenches. Nevertheless, SWMM satisfactorily simulated the hydrologic performance of eight of the nine LID practices.     

NETWORK AND SUBWATERSHED PARAMETERS EXTRACTED FROM DIGITAL ELEVATION MODELS: THE BILLS CREEK EXPERIENCE1

ABSTRACT: An automated extraction of channel network and sub-watershed characteristics from digital elevation models (DEM) is performed by model DEDNM. This model can process DEM data of limited vertical resolution representing low relief terrain. Such representations often include ill-defined drainage boundaries and indeterminate flow paths. The application watershed is an 84 km² low relief watershed in southwestern Oklahoma. The standard for validation is the network and subwatershed parameters defined by the blue line method on USGS 7.5–minute topographic maps. Evaluation of the generated and validation networks by visual comparisons shows a high degree of correlation. Comparison of selected network parameters (channel length, slope, drainage density, etc.) and of drainage network composition (bifurcation, length, slope, and area ratios) shows that, on the average, the generated parameters are within 5 percent of those derived from the validation network. The largest discrepancies were found for the channel slope values. The results of this application demonstrate that DEDNM effectively addresses network definition problems often encountered in low relief terrain and that it can generate accurate network and subwatershed parameters under those conditions.     

CHANNEL NETWORK EXTENSION BY LOGGING ROADS IN TWO BASINS,WESTERN CASCADES,OREGON1

ABSTRACT: Based on field surveys and analysis of road networks using a geographic information system (GIS), we assessed the hydrologic integration of an extensive logging-road network with the stream network in two adjacent 62 and 119 km² basins in the western Cascades of Oregon. Detailed surveys of road drainage for 20 percent of the 350 km road network revealed two hydrologic flow paths that link roads to stream channels: roadside ditches draining to streams (35 percent of the 436 culverts examined), and roadside ditches draining to culverts with gullies incised below their outlets (23 percent of culverts). Gully incision is significantly more likely below culverts on steep (< 40 percent) slopes with longer than average contributing ditch length. Fifty-seven percent of the surveyed road length is connected to the stream network by these surface flowpaths, increasing drainage density by 21 to 50 percent, depending on which road segments are assumed to be connected to streams. We propose a conceptual model to describe the hydrologic function of roads based on two effects: (1) a volumetric effect, increasing the volume of water available for quickflow and (2) a timing effect, altering flow-routing efficiency through extensions to the drainage network. This study examines the second of these two effects. Future work must quantify discharge along road segments connected to the stream network in order to more fully explain road impacts on basin hydrology.     

Spatial Relationships of Levees and Wetland Systems within Floodplains of the Wabash Basin,USA

Given the unique biogeochemical, physical, and hydrologic services provided by floodplain wetlands, proper management of river systems should include an understanding of how floodplain modifications influence wetland ecosystems. The construction of levees can reduce river–floodplain connectivity, yet it is unclear how levees affect wetlands within floodplains, let alone the cumulative impacts within an entire watershed. This paper explores spatial relationships between levee and floodplain wetland systems in the Wabash Basin, United States. We used a hydrogeomorphic floodplain delineation technique to map floodplain extents and identify wetlands that may be hydrologically connected to river networks. We then spatially examined the relationship between levee presence, wetland area, and other river network attributes within discrete subbasins. Our results show that cumulative wetland area is relatively constant in subbasins that contain levees, regardless of maximum stream order within the subbasin. In subbasins that do not contain levees, cumulative wetland area increases with maximum stream order. However, we found that wetland distributions around levees can be complex, and further studies on the influence of levees on wetland habitat may need to consider finer resolution spatial scales.     

CHANGING RAINFALL-RUNOFF RELATIONSHIPS IN THE URBANIZING PEACHTREE CREEK WATERSHED,ATLANTA, GEORGIA1

ABSTRACT: Peachtree Creek is a gaged watershed that has experienced a substantial increase in urbanization. The relationships of runoff to rainfall were studied for total annual flows, low flows, and peak flows. For each type of flow the relationship in the later, more urbanized period was compared to that in the earlier, less urbanized period. An increase in total runoff in wet years was observed as urbanization increased, but a decrease occurred during dry years. For low flows a similar decrease of runoff in dry years was found. An increase in peak runoff was observed over most of the range of precipitation. Increasing peak flows and declining low flows can be adequately explained by urban hydrologic theoryshed. which focuses on the effects of urban impervious surfaces upon direct runoff and infiltration. However, a decline of total runoff in dry years can be explained only by taking into account evapotranspiration as well. The concept of advectively assisted urban evapotranspiration, previously discovered by climatologists, is needed to explain such a loss of total runoff. Urban hydrologic theory must take into account vegetation and evapotranspiration, as well as impervious surfaces and their direct runoff, to explain the magnitude of total annual flows and low flows. Urban stormwater management should address the restoration of low flows, as well as the control of floods.     

DUAL URBAN AND RURAL HYDROGRAPH SIGNALS IN THREE SMALL WATERSHEDS1

ABSTRACT: Many studies can be found in the literature pertaining to the effects of urbanization on surface runoff in small watersheds and the hydrologic response of undeveloped watersheds. However, an extensive literature review yielded few published studies that illustrate differing hydrologic responses from multiple source areas within a watershed. The concepts discussed here are not new, but the methods used provide a unique, basic procedure for investigating stormwater hydrology in topographically diverse basins. Six storm hydrographs from three small central Pennsylvania watersheds were analyzed for this paper; five are presented. Two important conclusions are deduced from this investigation. First, in all cases we found two distinct peaks in stream discharge, each representing different contributing areas to direct discharge with greatly differing curve numbers and lags representative of urban and rural source regions. Second, the direct discharge represents only a small fraction of the total drainage area with the urban peak becoming increasingly important with respect to the rural peak with the amount of urbanization and as the magnitude of the rain event decreases.     

EFFECTS OF THE URBAN ENVIRONMENT ON HEAVY RAINFALL DISTRIBUTION1

ABSTRACT: A network of 225 recording raingages was operated over an area of 5200 km² in the St. Louis region during 1971-1975, in conjunction with an extensive investigation of urban effects on precipitation. Study of urban-induced effects on the frequency of heavy rainstorms has revealed a pronounced increase in the occurrence of storms producing 25 mm (1 inch) or more of rain. The increase is greatest in an area that is frequently in the path of storms passing across two urban-industrial regions. Analyses of raincells (rain intensity centers) within heavy convective storms shows a pronounced increase in water yield from cells exposed to potential urban effects, compared with those exposed only to the surrounding rural environment. Naturally-occurring heavy cells tend to undergo the greatest enhancement from urban exposure. Other analyses indicate an above-average frequency of excessive rain rates for periods of five minutes to two hours downwind of the urban-industrial complex. It is concluded that urban-induced intensification of short-duration rainstorms is sufficient to merit inclusion in the design and operation of urban-area hydrologic systems that control the flow of surplus storm water.     

ASSESSMENT OF POSSIBLE IMPACTS OF CLIMATE CHANGE IN AN URBAN CATCHMENT1

ABSTRACT: Stationarity of rainfall statistical parameters is a fundamental assumption in hydraulic infrastructure design that may not be valid in an era of changing climate. This study develops a framework for examining the potential impacts of future increases in short duration rainfall intensity on urban infrastructure and natural ecosystems of small watersheds and demonstrates this approach for the Mission/Wagg Creek watershed in British Columbia, Canada. Nonstationarities in rainfall records are first analyzed with linear regression analysis, and the detected trends are extrapolated to build potential future rainfall scenarios. The Storm Water Management Model (SWMM) is used to analyze the effects of increased rainfall intensity on design peak flows and to assess future drainage infrastructure capacity according to the derived scenarios. While the framework provided herein may be modified for cases in which more complex distributions for rainfall intensity are needed and more sophisticated stormwater management models are available, linear regressions and SWMM are commonly used in practice and are applicable for the Mission/Wagg Creek watershed. Potential future impacts on stream health are assessed using methods based on equivalent total impervious area. In terms of impacts on the drainage infrastructure, the results of this study indicate that increases in short duration rainfall intensity may be expected in the future but that they would not create severe impacts in the Mission/Wagg Creek system. The equivalent levels of imperviousness, however, suggest that the impacts on stream health could be far more damaging.     

Coupling PCSWMM and WASP to Evaluate Green Stormwater Infrastructure Impacts to Storm Sediment Loads in an Urban Watershed

We coupled rainfall–runoff and instream water quality models to evaluate total suspended solids (TSS) in Wissahickon Creek, a mid‐sized urban stream near Philadelphia, Pennsylvania. Using stormwater runoff and instream field data, we calibrated the model at a subdaily scale and focused on storm responses. We demonstrate that treating event mean concentrations as a calibration parameter rather than a fixed input can substantially improve model performance. Urban stormwater TSS concentrations vary widely in time and space and are difficult to represent simply. Suspended and deposited sediment pose independent stressors to stream biota and model results suggest that both currently impair stream health in Wissahickon Creek. Retrofitting existing detention basins to prioritize infiltration reduced instream TSS loads by 20%, suggesting that infiltration mitigates sediment more effectively than detention. Infiltrating stormwater from 30% of the watershed reduced instream TSS loads by 47% and cut the frequency of TSS exceeding 100 mg/L by half. Settled loads and the frequency of high TSS values were reduced by a smaller fraction than suspended loads and duration at high TSS values. A widely distributed network of infiltration‐focused projects is an effective stormwater management strategy to mitigate sediment stress. Coupling rainfall–runoff and water quality models is an important way to integrate watershed‐wide impacts and evaluate how management directly affects urban stream health.