SIZING OF SURFACE WATER RUNOFF DETENTION PONDS FOR WATER QUALITY IMPROVEMENT1

ABSTRACT: Historically, storm water management programs and criteria have focused on quantity issues related to flooding and drainage system design. Traditional designs were based on large rainfall‐runoff events such as those having two‐year to 100‐year return periods. While these are key criteria for management and control of peak flows, detention basin designs based on these criteria may not provide optimal quality treatment of storm runoff. As evidenced by studies performed by numerous public and private organizations, the water quality impacts of storm water runoff are primarily a function of more frequent rainfall‐runoff events rather than the less frequent events that cause peak flooding. Prior to this study there had been no detailed investigations to characterize the variability of the more frequent rainfall events on Guam. Also, there was a need to develop some criteria that could be applied by designers, developers, and agency officials in order to reduce the impact of storm water runoff on the receiving bodies. The objectives of this paper were three‐fold: (1) characterize the hourly rainfall events with respect to volume, frequency, duration, and the time between storm events; (2) evaluate the rainfall‐runoff characteristics with respect to capture volume for water quality treatment; and (3) prepare criteria for sizing and designing of storm water quality management facilities. The rainfall characterization studies have provided insight into the characteristics of rainstorms that are likely to produce non‐point source pollution in storm water runoff. By far the most significant fmdings are the development of a series of design curves that can be used in the actual sizing of storm water detention and treatment facilities. If applied correctly, these design curves could lead to a reduction of non‐point source pollution to Guam's streams, estuaries, and coastal environments.     

EFFECT OF SURFACE MINING ON STORM FLOW AND PEAK FLOW FROM SIX SMALL BASINS IN EASTERN KENTUCKY1

ABSTRACT: Hydrologic records from six small Eastern Kentucky watersheds were analyzed to determine the effect of surface mining on storm flows and peak flows. Average storm flow volumes were not changed by surface mining, whereas average peak flows were increased 36 percent. Peak flow increases were only in the summer. Smaller peak flows are doubled; moderate ones are increased by about a third; peak flows around 100 csm seem to be largely unaffected; and the larger peak flows may have been reduced by surface mining. The maximum annual storm flows, usually in winter or spring, appeared slightly reduced. No time trend in either storm flows or peak flows could be detected in five years of postmining record. Surface mining is not a serious floodwater discharge problem.     

Preferential flow estimates to an agricultural tile drain with implications for glyphosate transport

Agricultural subsurface drains, commonly referred to as tile drains, are potentially significant pathways for the movement of fertilizers and pesticides to streams and ditches in much of the Midwest. Preferential flow in the unsaturated zone provides a route for water and solutes to bypass the soil matrix and reach tile drains faster than predicted by traditional displacement theory. This paper uses chloride concentrations to estimate preferential flow contributions to a tile drain during two storms in May 2004. Chloride, a conservative anion, was selected as the tracer because of differences in chloride concentrations between the two sources of water to the tile drain, preferential and matrix flow. A strong correlation between specific conductance and chloride concentration provided a mechanism to estimate chloride concentrations in the tile drain throughout the storm hydrographs. A simple mixing analysis was used to identify the preferential flow component of the storm hydrograph. During two storms, preferential flow contributed 11 and 51% of total storm tile drain flow; the peak contributions, 40 and 81%, coincided with the peak tile drain flow. Positive relations between glyphosate [N-(phosphonomethyl)glycine] concentrations and preferential flow for the two storms suggest that preferential flow is an important transport pathway to the tile drain.     

MEASURED AND PREDICTED VELOCITY AND LONGITUDINAL DISPERSION AT STEAI)Y AND UNSTEADY FLOW,COLORADO RIVER,GLEN CANYON DAM TO LAKE MEAD1

ABSTRACT: The effect of unsteadiness of dam releases on velocity and longitudinal dispersion of flow was evaluated by injecting a fluorescent dye into the Colorado River below Glen Canyon Dam and sampling for dye concentration at selected sites downstream. Measurements of a 26-kilometer reach of Glen Canyon, just below Glen Canyon Dam, were made at nearly steady dam releases of 139, 425, and 651 cubic meters per second. Measurements of a 380-kilometer reach of Grand Canyon were made at steady releases of 425 cubic meters per second and at unsteady releases with a daily mean of about 425 cubic meters per second. In Glen Canyon, average flow velocity through the study reach increased directly with discharge, but dispersion was greatest at the lowest of the three flows measured. In Grand Canyon, average flow velocity varied slightly from subreach to subreach at both steady and unsteady flow but was not significantly different at steady and unsteady flow over the entire study reach. Also, longitudinal dispersion was not significantly different during steady and unsteady flow. Long tails on the time-concentration curves at a site, characteristic of most rivers but not predicted by the one-dimensional theory, were not found in this study. Absence of tails on the curves shows that, at the measured flows, the eddies that are characteristic of the Grand Canyon reach do not trap water for a significant length of time. Data from the measurements were used to calibrate a one-dimensional flow model and a solute-transport model. The combined set of calibrated flow and solute-transport models was then used to predict velocity and dispersion at potential dam-release patterns.     

Comparison of Unsteady and Quasi‐Unsteady Flow Models in Simulating Sediment Transport in an Ephemeral Arizona Stream1

Hummel, Ryan, Jennifer G. Duan, and Shiyan Zhang, 2012. Comparison of Unsteady and Quasi‐Unsteady Flow Models in Simulating Sediment Transport in an Ephemeral Arizona Stream. Journal of the American Water Resources Association (JAWRA) 48(5): 987‐998. DOI: 10.1111/j.1752‐1688.2012.00663.x Abstract: Hydrodynamic and sediment transport models are useful engineering tools for predicting unsteady flood flow and sediment transport. Many models such as HEC‐RAS, HEC‐6, and IALLUVIAL apply quasi‐unsteady flow model, whereas others apply the unsteady flow model. It remains unknown if a quasi‐unsteady flow model is sufficiently accurate for simulating sediment transport in rapidly varied unsteady flood events, especially in ephemeral rivers in arid and semiarid regions. This study compared the quasi‐unsteady HEC‐RAS 4.1 model with one‐dimensional (1D) Finite Volume Method (FVM) based model in simulating flood flow and sediment transport in the Pantano Wash, a dryland river in the state of Arizona. The objective is to determine which sediment transport method is appropriate in predicting bed elevation changes in an ephemeral stream, Pantano Wash, and if an unsteady model is more accurate than a quasi‐unsteady flow model in predicting sediment transport. Results showed that the quasi‐unsteady HEC‐RAS model and the 1D FVM yielded similar results of bed degradation and aggradation for this dryland stream, although the FVM model predicted better flood hydrographs. Among the seven sediment transport formulas embedded in HEC‐RAS, Yang’s and Engelund‐Hansen’s equations gave the best matches with the field measurements for this particular case study.     

HYDROLOGIC IMPACTS OF ALTERNATIVE APPROACHES TO STORM WATER MANAGEMENT AND LAND DEVELOPMENT1

Low impact development (LID) and other land development methods have been presented as alternatives to conventional storm water management and site design. Low impact development encourages land preservation and use of distributed, infiltration‐based storm water management systems to minimize impacts on hydrology. Such systems can include shallow retention areas, akin to natural depression storage. Other approaches to land development may emphasize land preservation only. Herein, an analysis of four development alternatives is presented. The first was Traditional development with conventional pipe/pond storm water management and half‐acre lots. The second alternative was Cluster development, in which implementation of the local cluster development ordnance was assumed, resulting in quarter‐acre lots with a pipe/pond storm water management system and open space preservation. The “Partial” LID option used the same lot layout as the Traditional option, with a storm water management system emphasizing shallow depression storage. The “Full” LID used the Cluster site plan and the depression storage‐based storm water management system. The alternatives were compared to the hydrologic response of existing site conditions. The analysis used two design storms and a continuous rainfall record. The combination of land preservation and infiltration‐based storm water management yielded the hydrologic response closest to existing conditions, although ponds were required to control peak flows for the design storms.     

DEBRIS FLOWS THROUGH DIFFERENT FOREST AGE CLASSES IN THE CENTRAL OREGON COAST RANGE1

ABSTRACT: Debris flows in the Pacific Northwest can play a major role in routing sediment and wood stored on hillslopes and in first‐through third‐order channels and delivering it to higher‐order channels. Field surveys following a large regional storm event investigated 53 debris flows in the central Oregon Coast Range to determine relationships among debris flow characteristics and the age class of the surrounding forest. The volume of sediment and wood delivered by debris flows was strongly correlated with runout length. Debris flows that initiated at roads were significantly longer than nonroad related failures, and road related landslides were an order of magnitude larger than nonroad related landslides. Clearcuts and roads tended to have more numerous contributing landslides relative to second growth and mature forests. No statistically significant difference in the average debris flow runout length was detected among the forest age classes, although debris flows initiating in clearcuts and mixed forest and at roads occasionally supported extremely long runout lengths that were outside the range of variability observed in completely forested basins. The size of wood in deposits was not correlated with the size of trees on the adjacent slopes, suggesting that the majority of wood in debris flow deposits was from remobilization of wood previously stored in low order channels.     

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.     

TMDL Balance: A Model for Coastal Water Pollutant Loadings

Bacterial contamination accounts for more than 60% of the impairments included on the 2008 Texas 303(d) List. Many of these bacterial impairments are along the Texas Gulf Coast because coastal waters often are regulated for oyster harvesting, which have strict water quality standards. Under the Clean Water Act, each one of these impaired waterbodies requires a total maximum daily load (TMDL) study to be performed. A recent, statewide study recommended the development and application of simple modeling approaches to address the majority of Texas's bacteria TMDLs, including “… simple load duration curve, GIS [geographic information systems], and/or mass balance models.” We developed the TMDL Balance model in response to this recommendation. TMDL Balance is a steady state, mass balance, GIS‐based model for simulating pollutant loads and concentrations in coastal systems. The model uses plug‐flow reactor and continuously‐stirred tank reactor equations to route spatially distributed point and nonpoint source loads through a watershed via overland flow, non‐tidal flow, and tidal flow, decaying the loads via first‐order kinetics. In this paper, we explain the development of the watershed loading portion of the TMDL Balance model, demonstrating the methodology through a case study: computing bacterial loads in the Copano Bay watershed of southeast Texas. The application highlights an example of distributing bacterial sources spatially based on land use data.     

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.     

LAKE MICHIGAN DIVERSION - STREAM QUALITY PLANNING1

ABSTRACT: By United States Supreme Court action, the diversion of water from Lake Michigan and the Lake Michigan Drainage Basin in the Metropolitan Chicago Area is regulated at an annual maximum rate of 3,200 cfs. Approximately 1,700 cfs of this diversion is used for water supply, and the remaining 1,500 cfs consists primarily of stormwater runoff with lesser amounts of direct lake diversion, such as lockage and leakage, navigational makeup water flows, and discretionary diversion needed to maintain water quality standards in the Metropolitan Sanitary District of Greater Chicago's basic waterways. In order to assess the schedule of its discretionary diversion needs, the District, using a computer model of its basic waterway system, has calculated the minimum discretionary diversion requirements for projected water quality conditions as successive elements of the District's water pollution control program are completed. The results of these analyses can be used as a basis for developing plans for future allocations of the limited supply of Lake Michigan water to other uses such as domestic water supply, when and if such supplies become available.     

Evaluating the Applicability of a Two‐dimensional Flow Model of a Highly Heterogeneous Domain to Flow and Environmental Management

Two‐dimensional simulation of highly heterogeneous domains, especially those with disparate length scales, roughness conditions, and geometries, often leads to challenges such as long computation times and numerical instability. Simulation of challenging domains is often needed to guide flood management and environmental regulation agencies in operation and potential domain modifications. This work evaluates the ability of a two‐dimensional unsteady hydrodynamic model to represent long‐duration transient flows over a domain with highly heterogeneous roughness, geometric characteristics, and length scales through bed roughness representation. The domain includes 13 km of Cache Creek and the 14.5 km² Cache Creek Settling Basin, which traps both sediment and mercury. Calibration under different bed roughness methods, validation, and modeling results of bathymetric modification scenarios are presented. The modeling approach's performance supports its application as a tool for management of similar domains, such as settling basins, leveed floodplains, and reservoirs. Accurate representation of flow dynamics can also inform environmental management that involves transport of sediments, nutrients, and heavy metals. This study found that a two‐dimensional unsteady flow model can accurately represent long‐duration transient flow in a large settling basin with highly heterogeneous characteristics without parsing of the domain or flow events simulated.     

INFLUENCE OF DIFFERENT TEMPORAL SAMPLING STRATEGIES ON ESTIMATING TOTAL PHOSPHORUS AND SUSPENDED SEDIMENT CONCENTRATION AND TRANSPORT IN SMALL STREAMS1

ABSTRACT: Various temporal sampling strategies are used to monitor water quality in small streams. To determine how various strategies influence the estimated water quality, frequently collected water quality data from eight small streams (14 to 110 km²) in Wisconsin were systematically subsampled to simulate typically used strategies. These subsets of data were then used to estimate mean, median, and maximum concentrations, and with continuous daily flows used to estimate annual loads (using the regression method) and volumetrically weighted mean concentrations. For each strategy, accuracy and precision in each summary statistic were evaluated by comparison with concentrations and loads of total phosphorus and suspended sediment estimated from all available data. The most effective sampling strategy depends on the statistic of interest and study duration. For mean and median concentrations, the most frequent fixed period sampling economically feasible is best. For maximum concentrations, any strategy with samples at or prior to peak flow is best. The best sampling strategy to estimate loads depends on the study duration. For one‐year studies, fixed period monthly sampling supplemented with storm chasing was best, even though loads were overestimated by 25 to 50 percent. For two to three‐year load studies and estimating volumetrically weighted mean concentrations, fixed period semimonthly sampling was best.     

Hydrologic Impact Assessment of Land Cover Change and Stormwater Management Using the Hydrologic Footprint Residence

Urbanization impacts the stormwater regime through increased runoff volumes and velocities. Detention ponds and low impact development (LID) strategies may be implemented to control stormwater runoff. Typically, mitigation strategies are designed to maintain postdevelopment peak flows at predevelopment levels for a set of design storms. Peak flow does not capture the extent of changes to the hydrologic flow regime, and the hydrologic footprint residence (HFR) was developed to calculate the area and duration of inundated land during a storm. This study couples a cellular automata land cover change model with a hydrologic and hydraulic framework to generate spatial projections of future development on the fringe of a rapidly urbanizing metropolitan area. The hydrologic flow regime is characterized for existing and projected land cover patterns under detention pond and LID‐based control, using the HFR and peak flow values. Results demonstrate that for less intense and frequent rainfall events, LID solutions are better with respect to HFR; for larger storms, detention pond strategies perform better with respect to HFR and peak flow.     

Water quality in headwater catchments with deer wallows

The pastoral grazing of farmed red deer (Cervus elaphus) is common in New Zealand. However, red deer have a natural instinct to seek out water and wallow in it. Often, in headwater catchments, they will create a wallow in a wet area connected to a waterway. Aesthetically, wallowing areas can be unpleasant and give the impression they are significant sources of contaminants entering waterways. This paper aimed to quantify their contribution to loads of contaminants lost from three headwater catchments (4.1 to 32.1 ha). Monthly water samples were taken of base flow and of all storm flow events and analyzed for nitrogen (N) and phosphorus (P) species, suspended sediment (SS), and the fecal indicator bacteria-E. coli. Median concentrations were generally in excess of recommended guidelines for lowland water quality and contact recreation in New Zealand (guidelines=9 microg dissolved reactive P L(-1), 30 microg total P L(-1), 444 microg nitrate and nitrite N L(-1), 0.9 mg NH4+-N L(-1) at pH 7, 4 mg SS L(-1), and 260 cfu 100 mL(-1)). Loads of P (up to c. 3 kg P ha(-1)) and SS (up to 4.5 Mg ha(-1)) exceeded the highest loads measured (1.7 kg P and 2 Mg SS ha(-1)) for a range of pastoral catchments in New Zealand, including deer-farmed catchments without many wallows connected to waterways. More losses occurred during storm flow than base flow but, more importantly, the majority of losses only occurred when deer were in the paddock and wallowing. Hence, to mitigate most contaminant losses, management should focus on discouraging wallowing and/or breaking the connectivity between wallows and waterways.     

Prediction of Total Sediment Load in Sand‐Bed Rivers in Korea Using Lateral Distribution Method

A new method for numerically predicting the total sediment load in a river is proposed. The method can be used to predict the total sediment load with information on channel geometry and slope, flow, and bed materials. The conventional method uses a 1D approach that assumes the channel has a wide rectangular shape. However, the proposed method computes depth‐averaged velocity over the width and predicts the total sediment load based on the flow computations. The new method, therefore, is expected to predict better if the flow changes significantly in the lateral direction. The proposed method was applied to three large sand‐bed rivers in Korea, where information is available regarding suspended sediment. Five formulas were tested of use in making total sediment load computations, namely Engelund‐Hansen's, Ackers‐White's, Yang's, Brownlie's, and Karim's formulas. The predicted total sediment loads are compared not only with measured data but also with results calculated using the 1D approach. Discrepancy ratios between the predicted and measured total sediment loads are given and the results are discussed.     

Suspended sediment, nitrogen and phosphorus concentrations and exports during storm-events to the Tuross estuary, Australia

This paper presents a process for estimating pollutant loads from water quality data, to improve catchment-scale modelling in the region for resource management purposes. It describes a program to estimate suspended sediment, total and dissolved nitrogen and phosphorus loads to the Tuross estuary from the Tuross River catchment (1810 km(2)) of coastal southeast Australia. Event-based water quality sampling results obtained during storm events in 2005 are presented. Event 1, during July 2005 was the largest storm event in terms of peak flow for 3.5 years. Other events monitored were also in July, November and December 2005. The early July 2005 event had a flow-weighted mean suspended sediment (SS) concentration during the first 4 days of 63 mg L(-1). Of the events monitored, this was unusual as it was preceded by drought and had the largest SS concentrations (peaking at 180 mg L(-1)) during the rising-stage. In contrast, the November event had a much lower flow-weighted SS mean (28 mg L(-1)), even though peak flow magnitudes were similar. The July and November 2005 events had peak flows of 12,360 and 11,330 ML day(-1). Low-cost rising-stage siphon samplers were used to collect samples during the rapidly rising phase of these events. The use of such samplers and consideration of time-lead/lag flow adjustments, quantified using cross-correlation analysis to account for hysteresis effects, were incorporated into the load estimation techniques. The technique is a potentially useful approach for understanding relationships between water quality concentrations and flow for modelling catchment source strengths and transport processes.     

Simulated interior floods and the flood experience in urban communities

The Hydrologic Engineering Center (HEC-1) model was used to construct synthetic hydrographs for isolated interior urban floods. Flood peak and lag time were very well preserved in simulated flows. Total volume was not adequately expressed. Lag time varied inversely with both urban development and storm intensity. Peak discharge varied with storm intensity, but this variability was well defined only at very high urbanization levels. An 175% increase in storm intensity produced a change of about 15% in peak discharge. Claims for flood damage correlated well with estimates of peak flow and lag time combined. Other measures of flood experience also correlated with the two features. Within the range of storms utilized, urban development factors consistently outranked storm intensity as a determining factor in flood damage.     

USING INFILTRATION ENHANCEMENT AND SOIL WATER MANAGEMENT TO REDUCE DIAZINON IN RUNOFF1

ABSTRACT: Pesticide runoff from dormant sprayed orchards is a major water quality problem in California's Central Valley. During the past several years, diazinon levels in the Sacramento and San Joaquin Rivers have exceeded water quality criteria for aquatic organisms. Orchard water management, via post‐application irrigation, and infiltration enhancement, through the use of a vegetative ground cover, are management practices that are believed to reduce pesticide loading to surface waters. Field experiments were conducted in Davis, California, to measure the effectiveness of these management practices in reducing the toxicity of storm water runoff. Treatments using a vegetative ground cover significantly reduced peak concentrations and cumulative pesticide mass in runoff for first flush experiments compared with bare soil treatments. Post‐application irrigation was found to be an effective means of reducing peak concentrations and cumulative mass in runoff from bare soil treatments, but showed no significant effect on vegetated treatments.     

CHALLENGES IN MODELING HYDROLOGIC AND WATER QUALITY PROCESSES IN RIPARIAN ZONES1

This paper presents key challenges in modeling water quality processes of riparian ecosystems: How can the spatial and temporal extent of water and solute mixing in the riparian zone be modeled? What level of model complexity is justified? How can processes at the riparian scale be quantified? How can the impact of riparian ecosystems be determined at the watershed scale? Flexible models need to be introduced that can simulate varying levels of hillslope‐riparian mixing dictated by topography, upland and riparian depths, and moisture conditions. Model simulations need to account for storm event peak flow conditions when upland solute loadings may either bypass or overwhelm the riparian zone. Model complexity should be dictated by the level of detail in measured data. Model algorithms need to be developed using new macro‐scale and meso‐scale experiments that capture process dynamics at the hillslope or landscape scales. Monte Carlo simulations should be an integral part of model simulations and rigorous tests that go beyond simple time series, and point‐output comparisons need to be introduced. The impact of riparian zones on watershed‐scale water quality can be assessed by performing simulations for representative hillsloperiparian scenarios.