AN EVALUATION OF FLOOD FREQUENCY ESTIMATES BASED ON RAINFALL/RUNOFF MODELING1

ABSTRACT: An evaluation of flood frequency estimates simulated from a rainfall/runoff model is based on (1) computation of the equivalent years of record for regional estimating equations based on 50 small stream sites in Oklahoma and (2) computation of the bias for synthetic flood estimates as compared to observed estimates at 97 small stream sites with at least 20 years of record in eight eastern states. Because of the high intercorrelation of synthetic flood estimates between watersheds, little or no regional (spatial) information may be added to the network as a result of the modeling activity. The equivalent years of record for the regional estimating equations based totally on synthetic flood discharges is shown to be considerably less than the length of rainfall record used to simulate the runoff. Furthermore, the flood estimates from the rainfall/runoff model consistently underestimate the flood discharges based on observed record, particularly for the larger floods. Depending on the way bias is computed, the synthetic estimate of the 100-year flood discharge varies from 11 to 29 percent less than the value based on observed record. In addition, the correlation between observed and synthetic flood frequency estimates at the same site is also investigated. The degree of correlation between these estimates appears to vary with recurrence interval. Unless the correlation between these two estimates is known, it is not possible to compute a weighted estimate with minimum variance.     

CLIMATE FACTOR FOR SMALL-BASIN FLOOD FREQUENCY1

ABSTRACT: A climate factor, C_T, (T = 2–, 25-, and 100-year recurrence intervals) that delineates regional trends in small-basin flood frequency was derived using data from 71 long-term rainfall record sites. Values of C_T at these sites were developed by a regression analysis that related rainfall-runoff model estimates of T-year floods to a sample set of 50 model calibrations. C_T was regionalized via kriging to develop maps depicting its geographic variation for a large part of the United States east of the 105th meridian. Kriged estimates of C_T and basin-runoff characteristics were used to compute regionalized T-year floods for 200 small drainage basins. Observed T-year flood estimates also were developed for these sites. Regionalized floods are shown to account for a large percentage of the variability in observed flood estimates with coefficients of determination ranging from 0.89 for 2-year floods to 0.82 for 100-year floods. The relative importance of the factors comprising regionalized flood estimates is evaluated in terms of scale (size of drainage area), basin-runoff characteristics (rainfall. runoff model parameters), and climate (C_T).     

UNCERTAINTY OF EXTREME FLOOD ESTIMATES INCORPRATING HISTORICAL DATA1

ABSTRACT: The ability to predict extreme floods is an important part of the planning process for any water project for which failure will be very costly. The length of a gage record available for use in estimating extreme flows is generally much shorter than the recurrence interval of the desired flows, resulting in estimates having a high degree of uncertainty. Maximum likelihood estimators of the parameters of the three parameter lognormal (3PLN) distribution, which make use of historical data, are presented. A Monte Carlo study of extreme flows estimated from samples drawn from three hypothetical 3PLN populations showed that inclusion of historical flows with the gage record reduced the bias and variance of extreme flow estimates. Asymptotic theory approximations of parameter variances and covariances calculated using the second and mixed partial derivatives of the log likelihood function agreed well with Monte Carlo results. First order approximations of the standard deviations of the extreme flow estimates did not agree with the Monte Carlo results. An alternative method for calculating those standard deviations, the “asymptotic simulation” method, is described. The standard deviations calculated by asymptotic simulation agree well with the Monte Carlo results.     

CONTINENTAL SCALE SIMULATION OF THE HYDROLOGIC BALANCE1

This paper describes the application of a continuous daily water balance model called SWAT (Soil and Water Assessment Tool) for the conterminous U.S. The local water balance is represented by four control volumes; (1) snow, (2) soil profile, (3) shallow aquifer, and (4) deep aquifer. The components of the water balance are simulated using “storage” models and readily available input parameters. All the required databases (soils, landuse, and topography) were assembled for the conterminous U.S. at 1:250,000 scale. A GIS interface was utilized to automate the assembly of the model input files from map layers and relational databases. The hydrologic balance for each soil association polygon (78,863 nationwide) was simulated without calibration for 20 years using dominant soil and land use properties. The model was validated by comparing simulated average annual runoff with long term average annual runoff from USGS stream gage records. Results indicate over 45 percent of the modeled U.S. are within 50 mm of measured, and 18 percent are within 10 mm without calibration. The model tended to under predict runoff in mountain areas due to lack of climate stations at high elevations. Given the limitations of the study, (i.e., spatial resolution of the data bases and model simplicity), the results show that the large scale hydrologic balance can be realistically simulated using a continuous water balance model.     

Combined effects of constant versus variable intensity simulated rainfall and reduced tillage management on cotton preemergence herbicide runoff

Pesticide runoff research relies heavily on rainfall simulation experiments. Most are conducted at a constant intensity, i.e., at a fixed rainfall rate; however, large differences in natural rainfall intensity is common. To assess implications we quantified runoff of two herbicides, fluometuron and pendimethalin, and applied preemergence after planting cotton on Tifton loamy sand. Rainfall at constant and variable intensity patterns representative of late spring thunderstorms in the Atlantic Coastal Plain region of Georgia (USA) were simulated on 6-m2 plots under strip- (ST) and conventional-tillage (CT) management. The variable pattern produced significantly higher runoff rates of both compounds from CT but not ST plots. However, on an event-basis, runoff totals (% applied) were not significantly different, with one exception: fluometuron runoff from CT plots. There was about 25% more fluometuron runoff with the variable versus the constant intensity pattern (P = 0.10). Study results suggest that conduct of simulations using variable intensity storm patterns may provide more representative rainfall simulation-based estimates of pesticide runoff and that the greatest impacts will be observed with CT. The study also found significantly more fluometuron in runoff from ST than CT plots. Further work is needed to determine whether this behavior may be generalized to other active ingredients with similar properties [low K(oc) (organic carbon partition coefficient) approximately 100 mL g(-1); high water solubility approximately 100 mg L(-1)]. If so, it should be considered when making tillage-specific herbicide recommendations to reduce runoff potential.     

ACCURACY AND PRECISION OF NRCS MODELS FOR SMALL WATERSHEDS1

ABSTRACT: In the last 30 years, the National Resource Conservation Service's TR‐55 and TR‐20 models have seen a dramatic increase in use for stormwater management purposes. This paper reviews some of the data that were originally used to develop these models and tests how well the models estimate annual series peak runoff rates for the same watersheds using longer historical data record lengths. The paper also explores differences between TR‐55 and TR‐20 peak runoff rate estimates and time of concentration methods. It was found that of the 37 watersheds tested, 25 were either over‐ or under‐predicting the actual historical watershed runoff rates by more than 30 percent. The results of this study indicate that these NRCS models should not be used to model small wooded watersheds less than 20 acres. This would be especially true if the watershed consisted of an area without a clearly defined outlet channel. This study also supports the need for regulators to allow educated hydrologists to alter pre‐packaged model parameters or results more easily than is currently permitted.     

MATHEMATICAL SIMULATIONS OF THE TOCCOA FALLS,GEORGIA, DAM-BREAK FLOOD1

ABSTRACT: Four dam-break models were selected for testing with an observed data set from the November 6, 1977, disaster at Toccoa Falls, Georgia. The Kelly Barnes Dam failure occurred with a 35-ft head of water and produced a peak discharge of 23,000 ft³/s. The selected models included: (1) Modified Puls (MP), (2) U. S. Army Corps of Engineers Gradually Varied Unsteady Flow Profiles (USTFLO), (3) National Weather Service's Dam-Break Flood Forecast (DBFF), and (4) U. S. Geological Survey's method of characteristics (MOC) coupled with a general purpose streamflow simulation (J879DB). Achieving a successful simulation was easiest with the MP model. The DBFF model required a moderate effort while the MOC-J879DB models required some data alterations and considerable effort. The USTFLO model failed to simulate this test case. In the stream segment near the dam, the computed peak stages were generally within 5 feet of the observed high water marks. Elsewhere, the peak stage results were much better, generally within 2 feet. The peak discharges computed by the models were generally within 20 percent of discharges estimated by slope area and contracted opening measurements, except near the dam where the MOC-J879DB model's results was 80 percent too high.     

Estimation of Catchment Nutrient Loads in New Zealand Using Monthly Water Quality Monitoring Data

Causes of variation between loads estimated using alternative calculation methods and their repeatability were investigated using 20 years of daily flow and monthly concentration samples for 77 rivers in New Zealand. Loads of dissolved and total nitrogen and phosphorus were calculated using the Ratio, L5, and L7 methods. Estimates of loads and their precision associated with short‐term records of 5, 10, and 15 years were simulated by subsampling. The representativeness of the short‐term loads was quantified as the standard deviation of the 20 realizations. The L7 method generally produced more realistic loads with the highest precision and representativeness. Differences between load estimates were shown to be associated with poor agreement between the data and the underlying model. The best method was shown to depend on the match between the model and functional and distributional characteristics of the data, rather than on the contaminant. Short‐term load estimates poorly represented the long‐term load estimate, and deviations frequently exceeded estimated imprecision. The results highlight there is no single preferred load calculation method, the inadvisability of “unsupervised” load estimation and the importance of inspecting concentration‐flow, unit load‐flow plots and regression residuals. Regulatory authorities should be aware that the precision of loads estimated from monthly data are likely to be “optimistic” with respect to the actual repeatability of load estimates.     

FOREST GROWTH IMPACT ON PEAK SNOW PACK MEASUREMENTS AT LONG TERM SNOW COURSES1

ABSTRACT: Snow course surveys in late winter provide stream‐flow forecasters with their best information for making water supply and flood forecasts for the subsequent spring and summer runoff period in mountainous regions of western North America. Snow survey data analyses are generally based on a 30‐year “normal” period. It is well documented that forest cover changes over time will affect snow accumulation on the ground within forests. This paper seeks to determine if forest cover changes over decades at long term snow courses decrease measured peak snow water equivalent (SWE) enough to affect runoff prediction. Annual peak SWE records were analyzed at four snow courses in two different forest types having at least 25 years of snowpack data to detect any decreases in SWE due to forest growth. No statistically significant decreases in annual peak SWE over time were found at any of these four snow courses. The wide range of annual winter precipitation and correspondingly highly variable peak snowpack accumulation, as well as many other weather and site variables, masked any minor trends in the data.     

A MONTE CARLO TEST OF LOAD CALCULATION METHODS,LAKE TAHOE BASIN,CALIFORNIA‐NEVADA1

ABSTRACT: The sampling of streams and estimation of total loads of nitrogen, phosphorus, and suspended sediment play an important role in efforts to control the eutrophication of Lake Tahoe. We used a Monte Carlo procedure to test the precision and bias of four methods of calculating total constituent loads for nitrate‐nitrogen, soluble reactive phosphorus, particulate phosphorus, total phosphorus, and suspended sediment in one major tributary of the lake. The methods tested were two forms of the Beale's Ratio Estimator, the Period Weighted Sample, and the Rating Curve. Intensive sampling in 1985 (a dry year) and 1986 (a wet year) provided a basis for estimating loads by the “worked record” method for comparison with estimates based on resampling actual data at the lower intensity that characterizes the present monitoring program. The results show that: (1) the Period Weighted Sample method was superior to the other methods for all constituents for 1985; and (2) for total phosphorus, particulate phosphorus, and suspended sediment, the Rating Curve gave the best results in 1986. Modification of the present sampling program and load calculation methods may be necessary to improve the precision and reduce the bias of estimates of total phosphorus loads in basin streams.     

Impact of 20 Years of Land-Cover Change on the Hydrology of Streams in the Southeastern United States1

Cruise, James F., Charles A. Laymon, and Osama Z. Al-Hamdan, 2010. Impact of 20 Years of Land-Cover Change on the Hydrology of Streams in the Southeastern United States. Journal of the American Water Resources Association (JAWRA) 46(6):1159–1170. DOI: 10.1111/j.1752-1688.2010.00483.x Abstract: Land-cover changes for portions of Alabama, Georgia, and Tennessee were estimated for the years 1980, 1990, and 2000 using classified Landsat images, and associated with hydrologic indices for 12 watersheds in the region. Rainfall-adjusted mean annual streamflow, an ET proxy (precipitation minus runoff), frequency of inundation above thresholds, and duration of inundation were used to characterize the hydrologic response of the test basins over the two-decade study period. Results indicate that several of the watersheds had undergone significant (>20%) reductions in agricultural land cover with a coincident increase in forested land. Attempts to correlate the hydrologic results with the land-cover changes were only partially successful. Watersheds with the largest land conversion from agriculture to forest (20% or more) did show significant trends in hydrologic indices indicating decreasing streamflow; however, other basins evidenced ambiguous results. The net conclusions of the study are that land-cover effects on hydrologic variables may be nuanced and can sometimes be only indirectly evident, and that a rigorous and detailed land-cover classification effort along with a battery of statistical tests with the same objective may be necessary to uncover these effects.     

Use of observed hydroclimatic trends to constrain projections of snowmelt season runoff in the Rio Grande headwaters

A method is developed for choosing 21st Century streamflow projections among widely varying results from a large ensemble of climate model-driven simulations. We quantify observed trends in climate–streamflow relationships in the Rio Grande headwaters, which has experienced warming temperature and declining snowpack since the mid-20th Century. Prominent trends in the snowmelt runoff season are used to assess corresponding statistics in downscaled global climate model projections. We define “Observationally Consistent (OC)” simulations as those that reproduce historical changes to linear statistics of diminished snowpack–streamflow coupling in the headwaters and an associated increase in the contribution of spring season (post-peak snowpack) precipitation to streamflow. Only a modest fraction of the ensemble of simulations meets these consistency metrics. The subset of OC simulations projects significant decreases in headwaters flow, whereas the simulations that poorly replicate historical trends exhibit a much wider range of projected changes. These results bolster confidence in model-based projections of declining runoff in the Rio Grande headwaters in the snowmelt runoff season and offer an example of a methodology for evaluating model-based projections in basins with similar hydroclimates that have experienced pronounced climate changes in the recent historical record.     

Impact of Residential Soil Disturbance on Infiltration Rate and Stormwater Runoff1

Woltemade, Christopher J., 2010. Impact of Residential Soil Disturbance on Infiltration Rate and Stormwater Runoff. Journal of the American Water Resources Association (JAWRA) 46(4): 700-711. DOI: 10.1111/j.1752-1688.2010.00442.x Abstract: Soil disturbances such as excavation and compaction in residential developments affect lawn infiltration rates and stormwater runoff. These effects were investigated via measuring saturated infiltration rates at 108 residential sites and 18 agricultural sites near Shippensburg, south-central Pennsylvania, using a double-ring infiltrometer. Residential sites included four neighborhoods distributed across three soil series classified as hydrologic soil group (HSG) B. Additional parcel data included date of house construction, percentage impervious area, lawn condition, and woody vegetation condition. Measured infiltration rates ranged from 0 to >40 cm/hour. Analysis of variance indicated significantly different mean infiltration rates (p < 0.001) for lots constructed pre-2000 (9.0 cm/hour) and those constructed post-2000 (2.8 cm/hour). Test results were used to determine a “field-tested” HSG for each site, representing disturbed soil conditions. Stormwater runoff was estimated from residential lots for a range of 24-hour design storms using the TR-55 model and several alternative methods of determining curve numbers, including five different representations of soil conditions. Curve numbers and stormwater runoff were substantially higher when based on field-tested HSGs for lots constructed post-2000 compared with lots built pre-2000 and when based on the HSG for undisturbed soils, documenting the magnitude of possible error in stormwater runoff models that neglect soil disturbance.     

A HYDROLOGIC RESPONSE MAP FOR THE STATE OF GEORGIA1

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

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.     

FLOOD-RUNOFF FORECASTING WITH HEC1F1

HEC1F is a computer program for making short- to medium-term forecasts of uncontrolled flood runoff. The program employs unit hydrographs and hydrologic routing to simulate runoff from a subdivided basin. Estimates of future rainfall can be accommodated. Runoff parameters for gaged headwater subbasins can be estimated (optimized) in real time. Blending of calculated with observed hydrographs can be performed. HEC1F is a component of an on-line software system that includes capability for data acquisition and processing, precipitation analysis, streamflow forecasting, reservoir system analysis, and graphical display of data and simulation results. The conceptual framework for HEC1F is described, and application of the program is illustrated.     

Phosphorus runoff relationships for Louisiana Coastal Plain soils amended with poultry litter

Long-term application of poultry (Gallus gallus domesticus) litter has built high levels of P in certain Coastal Plain soils of north Louisiana. However, soil P/runoff P relationships for soil and environmental conditions of the area have not been examined. This study measured soil P (total, Bray 1, Bray 2, Mehlich 3, resin-exchangeable, and water-extractable) and runoff P (dissolved P, DP; and total P, TP) at four pasture sites previously amended with poultry litter. Sites varied in soil P due to annual litter applications ranging from 1 to more than 20. Three replicated plots at each site were subjected to simulated rainfalls over 2 yr, and concentrations of DP and TP in runoff were measured and related to soil P. This allowed examination of soil P/runoff P relationships and their changes over time. Runoff DP was also related to DP desorbed from surface soil in a miscible displacement experiment. Among measures of soil P, only resin-exchangeable and water-extractable P showed significant decreases over 2 yr. These measures of soil P explained 54 to 64% of the variability in runoff DP data. However, the miscible displacement technique proved the best indicator of runoff DP, explaining 70% of the variability. Runoff varied among sites (decreasing with increasing years of litter application), limiting the predictive capability of the soil extraction methods. Linking runoff characteristics with miscible displacement data may be a useful predictive tool and warrants further examination.     

RELATIONSHIPS BETWEEN FLOOD FREQUENCIES AND RIPARIAN PLANT COMMUNITIES IN THE UPPER KLAMATH BASIN,OREGON1

ABSTRACT: This study was conducted in the Klamath Basin of southwestern Oregon to evaluate the dependency of riparian plant communities upon infrequent flooding. Plant communities were sampled with 1 m² quadrats along established cross‐sections. Data collected for purposes of hydraulic modeling included channel and floodplain elevations (i.e., cross‐sectional profiles) and water surface elevations associated with specific discharges. The elevational distribution of hydrophytic plant communities relative to modeled return periods provided the basis for establishing relationships between these variables for nine sites. Results indicate that, on average, a peak flow frequency of 4.6 years (range of 3.1 to 7.6 years) was needed to sustain riparian plant communities at seven of nine sites. At two sites, results indicated return periods of more than 25 years were needed; these results possibly were influenced by local geomorphic conditions (a narrow steep channel in one case and an incised channel in the other). Overall, these results tend to confirm a strong dependency of riparian plant communities on overbank flows.     

AN EXPANDING ROLE FOR SUBARCTIC WATERSHED RESEARCH1

ABSTRACT: It is increasingly recognized that natural resources research should in many cases be broadened in scope and oriented toward more general “environmental” problems. Locales with a history of “watershed” research can be eminently suited for development of comprehensive, environmental research programs. This is recognized in many research efforts of the International Biological Program (IBP), where watershed research sites have been successfully utilized for intensive investigations of process and function of selected ecosystems or ecosystem components. In the North American Subarctic there is almost no history of “watershed” studies. Basic data on hydrometeorologic parameters such as precipitation amounts and areal and seasonal distribution of runoff are scarce; the data framework within which environmental understanding can be structured is exceedingly sketchy. Opportunity exists in the discontinuous-permafrost settings of central Alaska to begin rectifying this situation. A basic program of multi-agency, multi-discipline research and data acquisition for the most significant hydrologic subregions is being developed, based around several existing environmental research areas (chiefly the Bonanza Creek Experimental Forest, the Caribou-Poker Creeks Research Watershed, the Wickersham Dome Fire Study Area, and a series of outlying sites).     

A Framework to Develop Nationwide Flooding Extents Using Climate Models and Assess Forecast Potential for Flood Resilience

The methods used to simulate flood inundation extents can be significantly improved by high‐resolution spatial data captured over a large area. This paper presents a hydraulic analysis methodology and framework to estimate national‐level floodplain changes likely to be generated by climate change. The hydraulic analysis was performed using existing published Federal Emergency Management Agency 100‐year floodplains and estimated 100‐ and 10‐year return period peak flow discharges. The discharges were estimated using climate variables from global climate models for two future growth scenarios: Representative Concentration Pathways 2.6 and 8.5. River channel dimensions were developed based on existing regional United States Geological Survey publications relating bankfull discharges with channel characteristics. Mathematic relationships for channel bankfull topwidth, depth, and side slope to contributing drainage area measured at model cross sections were developed. The proposed framework can be utilized at a national level to identify critical areas for flood risk assessment. Existing hydraulic models at these “hot spots” could be repurposed for near–real‐time flood forecasting operations. Revitalizing these models for use in simulating flood scenarios in near–real time through the use of meteorological forecasts could provide useful information for first responders of flood emergencies.