A TRIANGULAR MODEL OF DIMENSIONLESS RUNOFF PRODUCING RAINFALL HYETOGRAPHS IN TEXAS1

ABSTRACT: A synthetic triangular hyetograph for a large data base of Texas rainfall and runoff is needed. A hyetograph represents the temporal distribution of rainfall intensity at a point or over a watershed during a storm. Synthetic hyetographs are estimates of the expected time distribution for a design storm and principally are used in small watershed hydraulic structure design. A data base of more than 1,600 observed cumulative hyetographs that produced runoff from 91 small watersheds (generally less than about 50 km²) was used to provide statistical parameters for a simple triangular shaped hyetograph model. The model provides an estimate of the average hyetograph in dimensionless form for storm durations of 0 to 24 hours and 24 to 72 hours. As a result of this study, the authors concluded that the expected dimensionless cumulative hyetographs of 0 to 12 hour and 12 to 24 hour durations were sufficiently similar to be combined with minimal information loss. The analysis also suggests that dimensionless cumulative hyetographs are independent of the frequency level or return period of total storm depth and thus are readily used for many design applications. The two triangular hyetographs presented are intended to enhance small watershed design practice in applicable parts of Texas.     

DEVELOPMENT OF DESIGN STORM HYETOGRAPHS FOR CINCINNATI,OHIO1

ABSTRACT A synthetic storm rainfall hyetograph for a one-year design frequency is derived from the one-year intensity-duration curve developed for Cincinnati, Ohio. Detailed rainfall data for a three-year period were collected from three raingages triangulating the Bloody Run Sewer Watershed, an urban drainage areas of 2380 acres'in Cincinnati, Ohio. The advancement of the synthetic storm pattern is obtained from an analysis of the antecedent precipitation immediately preceding the maximum period of three selected durations. Rains which produced excessive runoff at least for some duration were considered only. The same approach can be used for other design frequencies. The purpose of this study is to provide synthetic storm hyetographs to be used as input in deterministic mathematical models simulating urban storm water runoff for the design, analysis and possible surcharge prediction of sewer systems.     

Assessment of the Extreme Rainfall Event at Nashville,TN and the Surrounding Region on May 1–3, 2010

This paper analyzes the May 1–3, 2010 rainfall event that affected the south‐central United States, including parts of Mississippi, Tennessee, and Kentucky. The storm is evaluated in terms of its synoptic setting, along with the temporal distributions, and spatial patterns of the rainfall. In addition, the recurrence interval of the storm is assessed and the implications for hydrologic structure designs are discussed. The event was associated with an upper‐level trough and stationary frontal boundary to the west of the rainfall region, which remained quasi‐stationary for a period of 48 h. Heavy rainfall was produced by two slow‐moving mesoscale convective complexes, combined with abundant atmospheric moisture. Storm totals exceeding 330 mm occurred within a large elongated area extending from Memphis to Nashville. Isolated rainfall totals over 480 mm were reported in some areas, with NEXRAD weather radar rainfall estimates up to 501 mm. An extreme value analysis was performed for one‐ and two‐day rainfall totals at Nashville and Brownsville, Tennessee, as well as for gridded rainfall estimates for the entire region using the Storm Precipitation Analysis System. Results suggest maximum rainfall totals for some durations during the May 1–3, 2010 event exceeded the 1,000‐year rainfall values from National Oceanic and Atmospheric Administration Atlas 14 for a large portion of the region and reached up to 80% of the probable maximum precipitation values for some area sizes and durations.     

Differences in Dissolved Organic Carbon and Nitrogen Responses to Storm‐Event and Ground‐Water Conditions in a Forested,Glaciated Watershed in Western New York1

Abstract: Differences in the storm‐event responses of dissolved organic carbon (DOC) and nitrogen (DON) in streamflow and ground water were evaluated for a glaciated forested watershed in western New York. Eight to ten storm events with varying rainfall amounts, intensities, and antecedent moisture conditions were studied for three catchments (1.6, 3.4, and 696 ha) over a three‐year period (2003‐2005). Concentrations of DOC in streamflow exiting the catchments were significantly higher for storm events following a dry period, whereas no similar response was observed for DON. Highest DON concentrations in streamflow were typically associated with storm events following wet antecedent moisture conditions. In addition to antecedent moisture conditions, DOC concentrations were also positively correlated with precipitation amounts, while DON did not reveal a consistent pattern. Streamwater and ground‐water concentrations of DOC during storm events were also strongly correlated with riparian ground‐water depths but a similar relationship was not observed for DON. Ground‐water DON concentrations were also more variable than DOC. We hypothesized that the differences in DOC and DON responses stemmed from the differences in catchment sources of these solutes. This study suggests that while DOC and DON are intrinsically linked as components of dissolved organic matter, their dynamics and exports from watersheds may be regulated by a different set of mechanisms and factors. Identifying these differences is critical for developing more reliable and robust models for transport of dissolved organic matter.     

Estimating Lag to Peak between Rainfall and Peak Streamflow with a Mixed‐Effects Model

We test the use of a mixed‐effects model for estimating lag to peak for small basins in Maine (drainage areas from 0.8 to 78 km²). Lag to peak is defined as the time between the center of volume of the excess rainfall during a storm event and the resulting peak streamflow. A mixed‐effects model allows for multiple observations at sites without violating model assumptions inherent in traditional ordinary least squares models, which assume each observation is independent. The mixed model includes basin drainage area and maximum 15‐min rainfall depth for individual storms as explanatory features. Based on a remove‐one‐site cross‐validation analysis, the prediction errors of this model ranged from ?42% to +73%. The mixed model substantially outperformed three published models for lag to peak and one published model for centroid lag for estimating lag to peak for small basins in Maine. Lag to peak estimates are a key input to rainfall–runoff models used to design hydraulic infrastructure. The improved accuracy and consistency with model assumptions indicates that mixed models may provide increased data utilization that could enhance models and estimates of lag to peak in other regions.     

PROBLEMS IN WEIGHTING OF HYETOGRAPHS1

ABSTRACT: It was found that the conventional weighting factor application to hyetograph ordinates results in artificially attenuated storm patterns. A modified weighting procedure is suggested which allows adjustments in the storm timing, peak intensity, and volume but conserves the storm pattern observed at the raingage nearest to the watershed point of interest. The systematic underestimation of peak flood flows, which result from conventional hyetograph weighting, can be avoided by conserving the hyetograph shape from the raingage nearest to any subarea of a modeled watershed and merely applying weighting factors to the rainfall volumes and temporal center of gravity of several hyetographs.     

Evaluation of Gauge‐Radar Merging Methods Using a Semi‐Distributed Hydrological Model in the Upper Thames River Basin,Canada

Gauge‐radar merging methods combine rainfall estimates from rain gauges and radar to capitalize on the strengths of the individual instruments. The performance of four well‐known gauge‐radar merging methods, including mean field bias correction, Brandes spatial adjustment, local bias correction using kriging, and conditional merging, are examined using Environment Canada radar and the Upper Thames River Basin in southwestern Ontario, Canada, as a case study. The analysis assesses the effect of gauge‐radar merging methods on: (1) the accuracy of predicted rainfall accumulations; and (2) the accuracy of predicted streamflows using a semi‐distributed hydrological model. In addition, several influencing factors (i.e., gauge density, storm type, basin type, proximity to the radar tower, and time‐step of adjustment) are analyzed to determine their effect on the performance of the rainfall estimation techniques. Confirming results of previous studies, the merging methods provide an increase in the accuracy of both rainfall accumulation estimations and predicted streamflows. The results also indicate specific factors such as gauge density, rainfall intensity, and time‐step of adjustment can reduce the accuracy of merging methods and play a key role in the examination of its use for operational purposes. Results provide guidance for hydrologists and engineers assessing how best to apply corrected radar products to improve rainfall estimation and hydrological modeling accuracy.     

Development and Operational Testing of a Super‐Ensemble Artificial Intelligence Flood‐Forecast Model for a Pacific Northwest River

Coastal catchments in British Columbia, Canada, experience a complex mixture of rainfall‐ and snowmelt‐driven contributions to flood events. Few operational flood‐forecast models are available in the region. Here, we integrated a number of proven technologies in a novel way to produce a super‐ensemble forecast system for the Englishman River, a flood‐prone stream on Vancouver Island. This three‐day‐ahead modeling system utilizes up to 42 numerical weather prediction model outputs from the North American Ensemble Forecast System, combined with six artificial neural network‐based streamflow models representing various slightly different system conceptualizations, all of which were trained exclusively on historical high‐flow data. As such, the system combines relatively low model development times and costs with the generation of fully probabilistic forecasts reflecting uncertainty in the simulation of both atmospheric and terrestrial hydrologic dynamics. Results from operational testing by British Columbia's flood forecasting agency during the 2013‐2014 storm season suggest that the prediction system is operationally useful and robust.     

Runoff Curve Numbers at the Agricultural Field‐Scale and Implications for Continuous Simulation Modeling

This study used monitoring in the waterways of agricultural fields to understand the use of the runoff curve number (CN) in continuous simulation models. The CN has a long history as a design tool for estimating runoff volumes for large, single storms on small watersheds, but its use in continuous simulation models to describe runoff from smaller storms and relatively small areas is more recent and controversial. We examined 788 nonwinter rainfall events on four agricultural fields over five years (2004‐2008) during which runoff was generated in 87 events. The largest 20 runoff events on each field generated approximately 90% of the total runoff volume. The runoff event CNs showed an inverse correlation with storm depth that could not consistently be explained by previous precipitation. We review how small areas of higher runoff generation within larger areas will systematically increase the apparent CN of the larger area as the storm size decreases. If this variation is not incorporated into a model explicitly, continuous simulation modelers must understand that when source areas are aggregated or when runoff generation is spatially variable, the overall CN is not unique when smaller storms are included in the calibration set.     

Long‐Term High‐Resolution Radar Rainfall Fields for Urban Hydrology

Accurate records of high‐resolution rainfall fields are essential in urban hydrology, and are lacking in many areas. We develop a high‐resolution (15 min, 1 km²) radar rainfall data set for Charlotte, North Carolina during the 2001‐2010 period using the Hydro‐NEXRAD system with radar reflectivity from the National Weather Service Weather Surveillance Radar 1988 Doppler weather radar located in Greer, South Carolina. A dense network of 71 rain gages is used for estimating and correcting radar rainfall biases. Radar rainfall estimates with daily mean field bias (MFB) correction accurately capture the spatial and temporal structure of extreme rainfall, but bias correction at finer timescales can improve cold‐season and tropical cyclone rainfall estimates. Approximately 25 rain gages are sufficient to estimate daily MFB over an area of at least 2,500 km², suggesting that robust bias correction is feasible in many urban areas. Conditional (rain‐rate dependent) bias can be removed, but at the expense of other performance criteria such as mean square error. Hydro‐NEXRAD radar rainfall estimates are also compared with the coarser resolution (hourly, 16 km²) Stage IV operational rainfall product. Stage IV is adequate for flood water balance studies but is insufficient for applications such as urban flood modeling, in which the temporal and spatial scales of relevant hydrologic processes are short. We recommend the increased use of high‐resolution radar rainfall fields in urban hydrology.     

RESPONSE OF GREEN‐AMPT MEIN‐LARSEN SIMULATED RUNOFF VOLUMES TO TEMPORALLY AGGREGATED PRECIPITATION1

ABSTRACT: Data collection frequency in automated systems is user determined and can range from seconds to hours or days. Currently, there is no standard or recommended frequency interval for collecting precipitation data from automated systems for input to event‐based models such as Green‐Ampt Mein‐Larsen (GAML). Data from 47 storm events at seven locations were used to simulate the response of GAML excess rainfall to temporally aggregated precipitation data. No difference in model efficiency was recognized when comparing one‐minute interval data (R²= 1.00) to five‐minute data (R²= 1.00). Very little model efficiency was lost at a 10‐minute (R²= 0.96) interval. After 10‐minutes, decline in efficiency became more rapid with R²= 0.16 at one hour. The combined effect of time interval with respect to drainage area, hydraulic conductivity, maximum 30‐minute intensity, and total precipitation also revealed similar results.     

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.     

Effects of Rainfall and Ground‐Water Pumping on Streamflow in Mākaha,O’ahu,Hawai’i1

Abstract: Land‐use/land‐cover changes in Mākaha valley have included the development of agriculture, residential dwellings, golf courses, potable water supply facilities, and the introduction of alien species. The impact of these changes on surface water and ground water resources in the valley is of concern. In this study, streamflow, rainfall, and ground‐water pumping data for the upper part of the Mākaha valley watershed were evaluated to identify corresponding trends and relationships. The results of this study indicate that streamflow declined during the ground‐water pumping period. Mean and median annual streamflow have declined by 42% (135 mm) and 56% (175 mm), respectively, and the mean number of dry stream days per year has increased from 8 to 125. Rainfall across the study area appears to have also declined though it is not clear whether the reduction in rainfall is responsible for all or part of the observed streamflow decline. Mean annual rainfall at one location in the study area declined by 14% (179 mm) and increased by 2% (48 mm) at a second location. Further study is needed to assess the effect of ground‐water pumping and to characterize the hydrologic cycle with respect to rainfall, infiltration, ground‐water recharge and flow in the study area, and stream base flow and storm flow.     

RAINFALL POINT INTENSITIES IN AN AIR MASS THUNDERSTORM ENVIRONMENT: WALNUT GULCH,ARIZONA1

ABSTRACT: Point rainfall intensities for a given return period are often used to formulate design storms for rainfall/runoff models to simulate design floods. These design floods are in turn used to design bridges, culverts, and a variety of drainage and flood control structures. The projected rapid growth in the southwestern United States will require very substantial monetary investments in drainage infrastructure. Accurate estimates of point rainfall intensities are critical to ensure both safe designs while not wasting dollars in overdesign. Rainfall point intensities (accumulated rainfall depth over a specified duration) for 5‐, 15‐, 30‐, and 60‐minute durations for the 2‐, 5‐, 10‐, 25‐, 50‐, and 100‐year return periods were determined for southeast Arizona. Thirty‐five years of rainfall record (1961 to 1995) were used in this study. The records came from 20 stations that were grouped into five sets of four independent stations to extend the rainfall records. The stations are in the USDA‐ARS Walnut Gulch Experimental Watershed (WGEW), which is representative of large portions of the Southwest whose runoff generation is dominated by air‐mass thunderstorms. The 5‐, 15‐, 30‐, and 60‐minute maximum intensities per year followed log‐normal distributions. The mean point rainfall intensities of the five sets of gages are very close (between 0 and 11 percent) to the NOAA values of the 5‐, 15‐, 30‐, and 60‐minute durations for all return periods. Much larger differences between the mean point rainfall intensities for all durations were found when these results were compared to those of a previous study done with a shorter rainfall record (between 14 and 33 percent for the 25‐, 50‐, and 100‐year return‐periods). The difference between the largest and the smallest values of point rainfall intensities recorded by each group, for all durations, usually increases as the return period increases.     

REGIONAL RAINFALL INTENSITY-DURATION-FREQUENCY CURVES FOR PENNSYLVANIA1

ABSTRACT: A statistical analysis of all available continuous hourly and 15-minute duration rainfall records for Pennsylvania was performed to develop an updated procedure to estimate design storms. As a resuit of this study, Pennsylvania was divided into five homogeneous rainfall regions and a set of rainfall intensity-duration curves developed for each region, for return periods of 1 to 100 years and durations ranging from 5 minutes to 24 hours. The PDT-IDF curves were judged to be a better representation of Pennsylvania rainfall than the nationwide TP-40 maps, particularly for storm events of 10-years and lower return periods. The average time distribution of 24-hour storms in Pennsylvania was found to be well represented by the SCS Type II distribution. The Corps of Engineers SPS 24-hour distribution was found to differ appreciably from both the SCS Type H and the Pennsylvania 24-hour storm distribution. For storm durations between 15 and 90 minutes the standard Yarnell intensity-duration curves closely resemble Pennsylvania storm distributions.     

AN INTEGRATED ONE‐DIMENSIONAL AND TWO‐DIMENSIONAL URBAN STORMWATER FLOOD SIMULATION MODEL1

ABSTRACT: Flash flooding is the rapid flooding of low lying areas caused by the stormwater of intense rainfall associated with thunderstorms. Flash flooding occurs in many urban areas with relatively flat terrain and can result in severe property damage as well as the loss of lives. In this paper, an integrated one‐dimensional (1‐D) and two‐dimensional (2‐D) hydraulic simulation model has been established to simulate stormwater flooding processes in urban areas. With rainfall input, the model simulates 2‐D overland flow and 1‐D flow in underground stormwater pipes and drainage channels. Drainage channels are treated as special flow paths and arranged along one or more sides of a 2‐D computational grid. By using irregular computation grids, the model simulates unsteady flooding and drying processes over urban areas with complex drainage systems. The model results can provide spatial flood risk information (e.g., water depth, inundation time and flow velocity during flooding). The model was applied to the City of Beaumont, Texas, and validated with the recorded rainfall and runoff data from Tropical Storm Allison with good agreement.     

Characterizing Storm Hydrograph Rise and Fall Dynamics With Stream Stage Data1

Abstract: Storm‐flow transients (i.e., hydrograph rise and fall dynamics) may represent an important aspect of understanding streamflow dynamics. However, little is known about how temporal resolution of transient data and climate variability may color these potential indicators of hydrologic pattern or condition. Warm‐season stream stage and rainfall were monitored continuously (5 min) during the 2002 water year in eight tributaries of the Little Miami River (Ohio), which drain 17‐58 km² catchments. Rise rates generated using 5‐min data were different than those generated with mean daily data [calculated with the Indicators of Hydrologic Alteration (IHA) software], though fall rates were similar for fine and coarse temporal data. This result suggests that data with low temporal resolution may not be adequate to fully represent the dynamics of storm rise rates. Conversely, fall rates based on daily stage data (via IHA) were similar to those based on the 5‐min data, and so daily mean data may be appropriate for characterizing fall rates. We next analyzed the possible correlations between rainfall variability and storm‐flow stage dynamics. We derived rise and recession rates from storm stage hydrographs by assuming exponential rise and decay of a runoff peak. We found that raw rise rates (R_raw) were correlated with both the maximum rainfall rate and the time to the centroid of a rain event. We subsequently removed the trend based on these rainfall characteristics, which yielded new representations of rise rates abbreviated as R_rate and R_tcent, respectively, and that had lower variability than the uncorrected (raw) data. Fall rates were found to be independent of rainfall characteristics. Due to the predominant influence of stream hydrology upon aquatic biota and nutrient fluxes, our work suggests that these stage data analysis protocols can refine or otherwise reduce variability in these indices by accounting for relevant factors such as rainfall forcing. These protocols for derivation of transient indices should be tested for their potential to improve correlations between stream hydrology and temporally aligned biotic data and dissolved nutrient fluxes in streams.     

Chloride persistence in a deiced access road drainage system

Three years of field data, classical linear reservoir theory, and a new dissolution model confirm the hypothesis that residual chloride from highway deicing applications dissolves into precipitation throughout the year. The measured input includes 52 storm hyetographs and logs of salt and premix applications on an access road with a closed drainage system subject to runoff, interflow, and baseflow. The output data feature discharge and conductivity in an outlet weir measured continuously from February 1998 to May 2000. Individual storm hydrographs and pollutographs yield calibrated first flush dissolved chloride concentrations and residual solid chloride loads that persist at appreciable levels over the entire period of record. The storm calibrations imply a source strength w of 2.01 x 10(-6) s(-1) that accurately models chloride dissolution kinetics through three salt seasons on the access road. This w rests on physically plausible values for the depression storage depth zeta (3 mm) and porosity n (0.40) that store the residual chloride.     

Effects of Forest Composition and Spatial Patterns on Storm Flows of a Small Watershed1

Abstract: The PRMS_Storm model was built as a storm event, distributed hydrological model for studying the hydrological effects of forest composition and spatial distribution on storm‐flow volume and peakflow rates in the Xiangshuixi Watershed in the Three Gorges Reservoir Area, in the Yangtze River Basin in southwestern China. We developed three simulation scenarios based on forest composition and their spatial arrangements across the watershed, including all mixed conifer‐evergreen broadleaf forests (Scenario 1), all mixed evergreen broadleaf forests (Scenario 2), and mixed conifer + evergreen broadleaf + shrub forests (Scenario 3). We examined 11 storm events observed during 2002‐2005. Compared with the existing forest covers, modeling results suggested that the amount of overland flow was reduced by 21, 23, and 22%, and the interflow increased by 16, 88, and 30%, for Scenarios 1, 2, and 3, respectively. During the same time, peakflow rates were reduced by 20.8, 9.6, and 18.9%, respectively. The reduction of peakflow rates was most significant when rainfall intensity exceeded 0.8 mm/min and events with a short duration and effect was minor when rainfall intensity was below 0.5 mm/min. In general, we found that Scenarios 1 and 3 were preferred for reducing storm‐flow volume and peakflow rates due to their higher interception rates, large soil water holding capacity, and higher soil infiltration capacity. The modeled results suggested soil properties are important in affecting the flow processes and thus forest composition and forest spatial distributions will affect storm‐flow volume and peakflow rates at the watershed scale. To maximize flood reduction functions of a watershed, high priority should be given to those forest types (Scenarios 1 and 3) in reforestation practices in the study region. This study suggests both forest composition and spatial pattern are important reforestation designs for flood reduction in the Three Gorges Reservoir Area.     

A STRIP MODEL APPROACH TO PARAMETERIZE A COUPLED GREEN‐AMPT KINEMATIC WAVE MODEL1

ABSTRACT: Infiltration processes at the plot scale are often described and modeled using a single effective hydraulic conductivity (Kg) value. This can lead to errors in runoff and erosion prediction. An integrated field measurement and modeling study was conducted to evaluate: (1) the relationship among rainfall intensity, spatially variable soil and vegetation characteristics, and infiltration processes; and (2) how this relationship could be modeled using Green and Ampt and a spatially distributed hydrologic model. Experiments were conducted using a newly developed variable intensity rainfall simulator on 2 m by 6 m plots in a rangeland watershed in southeastern Arizona. Rainfall application rates varied between 50 and 200 mm/hr. Results of the rainfall simulator experiments showed that the observed hydrologic response changed with changes in rainfall intensity and that the response varied with antecedent moisture condition. A distributed process based hydrologic simulation model was used to model the plots at different levels of hydrologic complexity. The measurement and simulation model results show that the rainfall runoff relationship cannot be accurately described or modeled using a single Kg value at the plot scale. Multi‐plane model configurations with infiltration parameters based on soil and plot characteristics resulted in a significant improvement over single‐plane configurations.