CAUSES OF PEAK FLOWS IN NORTHWESTERN MONTANA AND NORTHEASTERN IDAHO1

ABSTRACT: Both catchment experiments and a review of hydrologic processes suggest a varying effect of forest harvest on the magnitude of peak flows according to the cause of those peak flows. In northwestern Montana and Northeastern Idaho, annual maximum flows can result from spring snowmelt, rain, mid-winter rain-on-snow, or rain-on-spring-snowmelt. Meteorologic and physical data were used to determine the cause of annual maximum flows in six basins which had the necessary data and were smaller than 150 mi². Rain-on-spring-snowmelt was the most frequent cause of annual maximum flows in all six basins, although there was a strong gradient in the magnitude and cause of peak flows from southwest to northeast. Less frequent mid-winter rain-on-snow events caused the largest flows on record in four basins. Mid-winter rain-on-snow should be distinguished from rain-on-spring-snowmelt because of differences in seasonal timing, the relative contributions of rain vs. snowmelt, and the projected effects of forest harvest. The effects of mixed flood populations on the flood-frequency curve varied from basin to basin. Annual maximum daily flows could not be reliably predicted from rainfall and snowmelt data.     

GEOMORPHIC PARAMETERS PREDICT HYDROGRAPH CHARACTERISTICS IN THE SOUTHWEST1

ABSTRACT: Critical design characteristics of ephermal runoff such as hydrograph rise time, duration, mean peak discharge, volume, peak-volume ratio, and maximum flood were related to physical basin parameters such as area, shape, slope, drainage density, basin relief, stream length, and combinations of these in intermontane watersheds representative of the Mexican Highland section of the Basin and Range Province. Parameters used were restricted to those easily obtainable from maps or aerial photographs. A parameter expressing basin shape and size was developed which proved to be as accurate a predictor as others used in existing prediction equations tested and was simpler and faster to derive. Simple prediction equations derived for hydrograph characteristics were all significant except for volume at the 5% level; three were significant at the 1% level. Relationships determined are applicable in semi-arid basins of the Southwest up to 60 square miles (155 km²) in area.     

FLOOD FREQUENCY CHARACTERISTICS OF SOME ARIZONA WATERSHEDS1

The flood frequency characteristics of 18 watersheds in southeastern Arizona were studied using the log-Boughton and the log-Pearson Type 3 distribution. From the flood frequency study, a generalized envelope for Q₁₀₀ for watersheds 0.01 to 4000 mi² in area has been produced for southeastern Arizona. The generalized envelope allows comparisons to be made among the relative flood characteristics of the watersheds used in the study and provides a conservative estimate of Q₁₀₀ for ungaged watersheds in the region.     

SAMPLING STRATEGIES FOR ESTIMATING ACUTE AND CHRONIC EXPOSURES OF PESTICIDES IN STREAMS1

ABSTRACT: The Food Quality Protection Act of 1996 requires that human exposure to pesticides through drinking water be considered when establishing pesticide tolerances in food. Several systematic and seasonally weighted systematic sampling strategies for estimating pesticide concentrations in surface water were evaluated through Monte Carlo simulation, using intensive datasets from four sites in northwestern Ohio. The number of samples for the strategies ranged from 4 to 120 per year. Sampling strategies with a minimal sampling frequency outside the growing season can be used for estimating time weighted mean and percentile concentrations of pesticides with little loss of accuracy and precision, compared to strategies with the same sampling frequency year round. Less frequent sampling strategies can be used at large sites. A sampling frequency of 10 times monthly during the pesticide runoff period at a 90 km² basin and four times monthly at a 16,400 km² basin provided estimates of the time weighted mean, 90^th, 95^th, and 99^th percentile concentrations that fell within 50 percent of the true value virtually all of the time. By taking into account basin size and the periodic nature of pesticide runoff, costs of obtaining estimates of time weighted mean and percentile pesticide concentrations can be minimized.     

THE USE OF FLOOD POTENTIAL INDICES FOR FLOOD PEAK ESTIMATION ON UNGAGED WATERSHEDS1

: The construction of a flood peak index map was attempted for use by hydrologists in the simple format of rainfall maps. Since flood peaks are highly dependent on watershed area, the effect of area was removed. By regression analysis flood peaks of 2.33 and 100-year return periods were found to be proportional to watershed area to the 0.8 and 0.7 powers, respectively. Therefore, indices C_{2 33}= Q_{2 33}/_A0.7 were completed at each gage and plotted on a Pennsylvania map. It was attempted to further remove some of the scatter by regression of C with several other watershed parameters like slope, percent forest cover, and watershed shape, but no significant correlation could be found. The index maps, drawn without attenuation of the scatter, can be used by hydrologists to compute flood peaks as Q = CAⁿ (with n = 0.8 and 0.7 for the 2.33 and 100-year flood peaks, respectively). Flood peak safety factors can be based on visual observation of the index variation in the vicinity of the location for which the flood peak estimate is needed.     

REGIONALIZATION AND PREDICTION OF FLOODS IN THE FRASER RIVER CATCHMENT,B.C.1

Probability distributions that model the return periods of flood characteristics derived from partial duration series are proposed and tested in the Fraser River catchment of British Columbia. Theoretical distributions describing the magnitude, duration, frequency and timing of floods are found to provide a goof fit to the observed data. The five estimated parameters summarizing the flood characteristics of each basin are entered into a discriminant analysis procedure to establish flood regions. Three regions were identified, each displaying flood behavior closely related to the physical conditions of the catchment. Within each region, regression equations are obtained between parameter values and basin climatic and physiographic variables. These equations provide a satisfactory prediction of flood parameters and this allows the estimation of a comprehensive set of flood characteristics for areas with sparse hydrologic information.     

REGIONALIZATION AND PREDICTION OF ANNUAL FLOODS IN THE FRASER RIVER CATCHMENT,BRITISH COLUMBIA1

A method of predicting probability distributions of annual floods is presented and is applied to the Fraser River catchment of British Columbia. The Gumbel distribution is found to adequately describe the observed flood frequency data. Using the estimated Gumbel parameters, discriminant analysis is performed to separate basins into flood regions. Within each region, regression analysis is used to relate physiographic and climatic variables to the means and standard deviations of the annual flood series. The regression equations are applied to four test basins and the results indicate that the method is suitable for an estimation of annual floods.     

A NEW LOOK AT FLOOD RISK DETERMINATION1

ABSTRACT: The probability distributions of annual peak flows used in flood risk analysis quantify the risk that a design flood will be exceeded. But the parameters of these distributions are themselves to a degree uncertain and this uncertainty increases the risk that the flood protection provided will in fact prove to be inadequate. The increase in flood risk due to parameter uncertainty is small when a fairly long record of data is available and the annual flood peaks are serially independent, which is the standard assumption in flood frequency analysis. But standard tests for serial independence are insensitive to the type of grouping of high and low values in a time series, which is measured by the Hurst coefficient. This grouping increases the parameter uncertainty considerably. A study of 49 annual peak flow series for Canadian rivers shows that many have a high Hurst coefficient. The corresponding increase in flood risk due to parameter uncertainty is shown to be substantial even for rivers with a long record, and therefore should not be neglected. The paper presents a method of rationally combining parameter uncertainty due to serial correlation, and the stochastic variability of peak flows in a single risk assessment. In addition, a relatively simple time series model that is capable of reproducing the observed serial correlation of flood peaks is presented.     

CORRELATION OF ANNUAL PEAK FLOWS FOR PENNSYLVANIA STREAMS1

A common assumption in flood frequency analysis is that annual peak flows are independent events. This study was undertaken to investigate the validity of this assumption with regard to Pennsylvania streams by statistically analyzing the dependence between annual peak flows and to determine if basin carryover effects relate to the degree of dependence. Five tests of dependence, the autocorrelation test, the median crossing test, the turning points test, the rank difference test, and the Spearman rank order serial correlation coefficient test were applied to the series of annual peak flows for 57 streams. Of the 57 streams analyzed, only two exhibited signs of dependence by at least two of the tests performed, and the baseflow component of annual peak flows was found to be unrelated to the degree of dependence exhibited between annual peak flows. It was concluded that the assumption of independence of annual peak flows is valid in flood frequency analysis for Pennsylvania streams.     

COMPARISON OF METHODS FOR ESTIMATING FLOOD MAGNITUDES ON SMALL STREAMS IN GEORGIA1

ABSTRACT: The U.S. Geological Survey has collected flood data for small, natural streams at many sites throughout Georgia during the past 20 years. Flood-frequency relations were developed for these data using four methods: (1) observed (log-Pearson Type HI analysis) data, (2) rainfall-runoff model, (3) regional regression equations, and (4) map-model combination. The results of the latter three methods were compared to the analyses of the observed data in order to quantify the differences in the methods and determine if the differences are statistically significant. Comparison of regression-estimates with observed discharges for sites having 20 years (1966 to 1985) and 10 years (1976 to 1985) of record at different sites of annual peak record indicate that the regression-estimates are not significantly different from the observed data. Comparison of rainfall-runoff-model simulated estimates with observed discharges for sites having 10 years and 20 years of annual peak record indicated that the model-simulated estimates are significantly and not significantly different, respectively. The biasedness probably is due to a “loss of variance” in the averaging procedures used within the model and the short length of record as indicated in the 10 and 20 years of record. The comparison of map-model simulated estimates with observed discharges for sites having 20 years of annual-peak runoff indicate that the simulated estimates are not significantly different. Comparison of “improved” map-model simulated estimates with observed discharges for sites having 20 years of annual-peak runoff data indicate that the simulated estimates are different. The average adjustment factor suggested by Lichty and Liscum to calculate the “improved” map-model overestimates in Georgia by an average of 20 percent for three recurrence intervals analyzed.     

FLOOD FREQUENCY ANALYSIS FOR SMALL WATERSHEDS IN SOUTHERN ILLINOIS1

ABSTRACT: Twenty-two gaging stations were selected for developing a regional flood frequency curve for small (area less than 2 square miles) watersheds in southern Illinois. Five probability functions were compared, and the extreme value type I function was selected to develop the regional flood curve. The curve was generated with the index flood method and also another empirical method that related the function parameters to the watershed area. Estimated peak discharges with various return periods were compared with the results obtained from multiple regression analysis.     

SUSPENDED SEDIMENT PRODUCTION FROM FORESTED WATERSHEDS ON OAHU,HAWAII1

Suspended sediment from forested and agricultural watersheds was sampled over a five-year period on the island of Oahu. A variety of storm conditions were sampled, giving a measure of the extreme variability in suspended sediment production. Total annual suspended sediment from all watersheds sampled ranged from 8400 kg/km² to 617,000 kg/km². Normally, about 90 percent of the total suspended sediment was produced during less than 2 percent of the time. Suspended sediment concentrations rapidly increased during rising stream flow resulting from rain storms. Time to peak of less than two hours is common, with a similarly rapid return to prestorm conditions. The data presented indicate the great variability of suspended sediment yields, making establishment of effective standards difficult.     

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.     

REGIONAL AND TEMPORAL PATTERNS OF TOTAL SOLUTES IN THE SASKATCHEWAN RWER BASIN1

ABSTRACT: In the Saskatchewan River Basin (365,000 km²), which drains the Canadian prairie from the Rocky Mountains east to Manitoba, concentrations of total solutes are usually within the range of 100 to 1000 mg/L. Total solutes levels in tributaries increase markedly from west to east across the basin, as mountain snowmelt and dilute surface runoff are replaced by ion-rich ground water and concentrated prairie runoff as the major influences on solute concentrations. In contrast, total solutes concentrations in main-stem rivers are nearly constant, ranging 200–300 mg/L, with only a small increase across the basin. Dilute mountain runoff dominates solute concentrations in main-stem rivers, despite the influx of increasingly ion-rich water from tributaries. The principal long-term trends in total solute concentrations across the basin, as revealed by linear and sine-curve regressions, were due to the construction of reservoirs, which depress the natural winter maximum in solute concentrations and disrupt the sinusoidal annual pattern, while sharply reducing seasonal variation. These regression methods did not show anticipated anthropogenic increases in salt load in the Red Deer or South Saskatchewan Rivers, but a trend of slowly increasing solutes concentrations (2 mg/L/yr) was detected for autumn flows in the lower Bow River. Municipal wastes from the City of Calgary or irrigation return flows are probably responsible for this increase.     

FORESTS AND OTHER FACTORS ASSOCIATED WITH STREAMFLOWS IN EAST TEXAS1

Effects of proportion of watersheds in forest and watershed physiographic factors on mean annual streamflow (1965-76), median flow, and 12 flood flow characteristics were regionally analyzed for 19 unregulated streams in East Texas. Annual streamflow increased with decreasing proportion of forest area. Differences in annual streamflow between full forest cover and bare watersheds could be as much as 200 mm. Other things being equal, the minimum watershed area required to generate 0.142 cm (5 cfs), a criterion used by the U.S. Corps of Engineering in regulating dredge and fill activity for water pollution abatement in East Texas streams, is 70 km² (27 mi²). Of the 31 physio-climatic parameters analyzed, watershed area, percent forest area, shape index, spring precipitation, and annual temperature were the most significant in affecting streamflow characteristics in East Texas. Using 2-3 of these five variables, all of the 14 streamflow characteristics can be estimated with accuracy ranging from acceptable to excellent levels.     

ANALYZING URBANIZATION IMPACTS ON PENNSYLVANIA FLOOD PEAKS1

ABSTRACT: The impact of man made change on the hydrology of developing watersheds is frequently measured in terms of the ratio: flood peak after development to flood peak before development over a range of return periods. However, the analysis of urbanization effects on flood frequency presents a vexing problem because of a general lack of flood data in urban areas and also because of nonstationarity in the development process. Clearly, the flood peak ratio depends on the impervious fraction and percent of basin sewered and these factors have been taken into account in recent urban flood peak models. In genral, these models are developed either by: (1) split sample analysis of available annual flood data, or (2) by computer simulation using mathematical watershed models capable of representing man made changes. The present paper discusses the results of work in progress to characterize the impact of urbanization on small developing watersheds in Pennsylvania.     

RIVERINE TRANSPORT OF NUTRIENTS AND DETRITUS TO THE APALACHICOLA BAY ESTUARY,FLORIDA1

ABSTRACT: The Applachicola River basin in northwest Florida covers an area of 3,100 square kilometers. Fifteen percent of the area is a dense bottomland hardwood forest which is periodically flooded. The annual leaf-litter fall from the flood-plain trees is a potential source of nutrients and detritus which eventually can flow into Apalachicola Bay. Transport of such material is dependent on the periodic inundation of the flood plain. The U.S. Geological Survey Apalachicola Rim Quality Assessment measured a total organic carbon flux of 2.1 × 10⁵ metric tons during the one-year period from June 3, 1979, to June 2,1980. Fluxes of total nitrogen and phosphorus during the same year were 2.1 × lo⁴ and 1.7 × lo³ metric tons, respectively. Flood characteirstics, such as prior hydrologic conditions, extent, and timing, are important in determining the amount and forms of materials transported. The 1980 spring flood produced a fourfold discharge increase over the annual mean outflow of 800 cubic meters per second. Nutrient concentrations varied little with discharge, but the 86-day spring flood accounted for 53, 60, 48, and 56 percent of the annual flux of total organic carbon, particulate organic carbon, total nitrogen, and total phosphorus, respectively. In 1980, the flood peaks, rather than the rise or recession, accounted for maximum nutrient and detritus transport.     

DRAINAGE GENERATION AND WATER USE IN THE EVERGLADES AGRICULTURAL AREA BASIN1

ABSTRACT: The Everglades Agricultural Area (EAA) covers 2,850 km² in area and is characterized by high water table and organic soil. The area is actively irrigated and drained as a function of weather conditions and crop status. Anthropogenic activities in the basin have resulted in nutrient-enriched drainage water that is discharged to Lake Okeechobee and the Everglades ecosystem. Water quantity and quality issues of the basin have become of increasing interest at local, state, and federal levels, so legislative and regulatory measures have been taken to improve water quality in discharges from the basin. In this study, simulation of hydrologic conditions and soil moisture were conducted using 100 years of daily synthetic rainfall data. From the simulations, the statistical distribution of half-month drainage discharge and supplemental water use in the basin was developed. The mean annual drainage/runoff was 49 cm, the mean supplemental water was 30 cm, and the mean annual a real rainfall was 122 cm. On the average, drainage exceeded supplemental water use in the months of June to September while from December to March drainage and supplemental water use were equivalent. Supplemental water use exceeded drainage in the months of October, November, April, and May. High drainage occurred in June and September; smallest drainage was in February. On the average, the highest supplemental water use occurred in May and November. The 10-year return period of annual drainage during wet and dry cycles were 60 cm and 38 cm per year, respectively. The semi-monthly drainage coefficient of variation (cv) is above 100 percent for the period from the second half of October to end of April. The cv is lower than 100 percent for the remaining season (wet season). The purpose of this paper is to present the magnitude, temporal, and frequency distribution of drainage runoff generation and supplemental water use in the EAA basin. Information on statistics of drainage will contribute to the optimization of the design and operation of drainage water treatment systems.     

AutoRAPID: A Model for Prompt Streamflow Estimation and Flood Inundation Mapping over Regional to Continental Extents

This article couples two existing models to quickly generate flow and flood‐inundation estimates at high resolutions over large spatial extents for use in emergency response situations. Input data are gridded runoff values from a climate model, which are used by the Routing Application for Parallel computatIon of Discharge (RAPID) model to simulate flow rates within a vector river network. Peak flows in each river reach are then supplied to the AutoRoute model, which produces raster flood inundation maps. The coupled tool (AutoRAPID) is tested for the June 2008 floods in the Midwest and the April‐June 2011 floods in the Mississippi Delta. RAPID was implemented from 2005 to 2014 for the entire Mississippi River Basin (1.2 million river reaches) in approximately 45 min. Discretizing a 230,000‐km² area in the Midwest and a 109,500‐km² area in the Mississippi Delta into thirty‐nine 1° by 1° tiles, AutoRoute simulated a high‐resolution (~10 m) flood inundation map in 20 min for each tile. The hydrographs simulated by RAPID are found to perform better in reaches without influences from unrepresented dams and without backwater effects. Flood inundation maps using the RAPID peak flows vary in accuracy with F‐statistic values between 38.1 and 90.9%. Better performance is observed in regions with more accurate peak flows from RAPID and moderate to high topographic relief.     

STREAM TEMPERATURE INCREASES AND LAND USE IN A FORESTED OREGON WATERSHED1

ABSTRACT: The Salmon Creek Watershed drains 325 km² of forested terrain in the Cascade Mountains of western Oregon. Over a 30–year period (from 1955 to 1984) average daily maximum and minimum stream temperatures, calculated from the 10 warmest days of each year, have risen 6°C and 2°C, respectively. In contrast, a small decrease in maximum air temperatures was found over the same period. Regression analysis indicated a highly significant (p < 0.01) relationship between a cumulative index of forest harvesting and maximum stream temperatures. Maximum temperatures also tended to increase for several years following major peak flow events. The interaction between harvest activity (logging and road construction), changing forest and riparian management practices and the occurrence of natural hydrologic events (peak flows and associated mass soil movements) tend to obscure specific cause-and-effect relationships regarding long-term changes in maximum stream temperature.