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.     

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.     

SOME SUSPENDED SEDIMENT YIELDS FROM SOUTH ISLAND CATCHMENTS,NEW ZEALAND1

: Estimates of specific annual suspended sediment yields, some of which rank among the highest reported in the world, are presented for 33 basins of South Island, New Zealand. Yield from each basin was determined by combining a suspended sediment concentration rating with the complete flow record of each catchment stream gaging station. A multiple regression analysis between sediment yields and climatic, hydrologic and physiographic parameters of each basin demonstrates that most of the variance in yields is explained by catchment mean rainfall. Geology apparently has little influence on sediment yield as suspended sediment concentration ratings, from rivers draining catchments of differing lithology, and regolith, are indistinguishable. Specific suspended sediment yield prediction equations are given for four defined regions covering in area almost all South Island; and except for one area, feature rainfall as the principle independent variable. Differences between regions may be due to variations in intensity, frequency, and duration patterns of storms. It is proposed that a simple power law relationship between yield and rainfall provides useful suspended sediment yield estimates in mountainous regions of temperate maritime climate, provided catchments have not been modified extensively by man.     

ALTERNATIVE APPROACH TO THE FORMULATION OF BASIN HYDRAULIC GEOMETRY EQUATIONS1

ABSTRACT: Along a drainage network, there is a systematic variation of average flow parameters (width, depth, and velocity) at flows having the same flow duration. Hydraulic geometry equations mathematically express this interdependent relationship of stream-flow characteristics for a basin for annual flow durations varying from 10 to 90 percent. However, the equations proposed so far have had rather poor predictive performance for low flows. An independent investigation of the variation of discharge with drainage area and annual flow duration demonstrates a consistent relationship between these parameters. The relationship for the high to median-flow range differs, however, from that for the median— to low-flow range. The proposed equations provide a better predictive performance for low flows than previous formulations and a versatile means of estimating flow parameters for streams throughout a basin. The improved basin hydraulic geometry equations have a wide range of applications in areas such as stream habitat assessment, water quality modeling, channel design, and stream restoration projects.     

MONTHLY RUNOFF DISTRIBUTION RELATED TO BASIN GEOGRAPHY1

ABSTRACT: Mean monthly runoff from ungaged drainage basins that have significant snowpacks each year can be estimated quite well by assuming that the time duration between snowfall and snowmelt is the predominant factor in temporal runoff distribution. That time span is related to basin temperatures which are, in turn, functions of basin elevation and latitude. Regional hydrologic analyses of gaged basin data create regression equations for estimating runoff distribution by month. These equations then can be applied to ungaged basins. Basin latitude and mean elevation are two independent variables that can be used in estimating monthly runoff distributions.     

PREDICTION OF PEAK FLOWS ON SMALL WATERSHEDS IN OREGON FOR USE IN CULVERT DESIGN1

ABSTRACT: Using data from 80 Oregon watersheds that ranged in size from 0.54 km² to 27.45 km², equations were developed to predict peak flows for use in culvert design on forest roads. Oregon was divided into six physiographic regions based on previous studies of flood frequency. In each region, data on annual peak flow from gaging stations with more than 20 years of record were analyzed using four flood frequency distributions: type 1 extremal, two parameter-log normal, three parameter-log normal, and log-Pearson type III. The log-Pearson type III distribution was found to be suitable for use in all regions of the State, based on the chi-square goodness-of-fit-test. Flood magnitudes having recurrence intervals of 10, 25, 50, and 100 years were related to physical and climatic characteristics of drainage basins by multiple regression. Drainage basin size was the most important variable in explaining the variation of flood peaks in all regions. Mean basin elevation and mean annual precipitation were also significantly related to flood peaks in two regions of western Oregon. The standard error of the estimate for the regression relationships ranged from 26 to 84 percent.     

PREDICTION OF MEAN ANNUAL FLOWS IN THE FRASER RIVER CATCHMENT,BRITISH COLUMBIA1

ABSTRACT: A method of predicting annual flows is presented and is applied to the Fraser River catchment. Statistical tests show the annual flow records to be stationary and aerially independent and can be adequately approximated by Gaussian distributions. Estimates are made of the Gaussian parameters for each subbasin. The spatial variations of these parameters are described by third order trend surfaces. The fitted surfaces can then be used to predict parameters of ungaged basins using the latitude and longitude of the basin centroids. The predicted parametric values are substituted into the Gaussian distribution to generate flows of various return periods.     

VERIFICATION AND APPLICATION OF REGIONAL EQUATIONS FOR THE STORAGE FUNCTION RUNOFF MODEL1

The storage function model is a nonlinear rainfall-runoff model that has been developed for and applied to flood runoff analysis in Japan. This paper extends the model applicability by developing practical equations for estimating model parameters which are appropriate on a regional basis, i.e., so-called regional equations. Previously, the parameters were computed from historical data for a specific basin or from relationships that do not account for land use and topography. To develop the regionalized equations, model parameters were identified for 91 flood events from 22 watersheds in Japan by applying a mathematical optimization technique. Results from 39 of these events were statistically compared and regional relationships were determined as a function of land use, basin area and rainfall intensity. The utility of the estimated equations were tested by computing runoff hydrographs for lumped basins. The estimated parameters were also applied in a distributed watershed model formulation. Both applications showed acceptable results that validate the use of the regionalized relationships.     

A MODIFIED RATIONAL FORMULA FOR FLOOD DESIGN IN SMALL BASINS1

ABSTRACT: New formulas and procedures under the framework of the Rational Formula are presented that are applicable to flood design problems for a small basin if the geometry of the basin can be approximated as an ellipse or a rhombus. Instead of making the assumption in the traditional rational formula that the rainfall is uniformly distributed in the whole duration (D_w) of a design storm, the new method modifies that assumption as: the rainfall is uniformly distributed only in each time interval CD) of the design storm hyetograph, thus extending the rational formula applicable to the case that the rainfall duration is less than the basin concentration time (T_c). The new method can be applied to estimate the flood design peak discharge, and to generate the flood hydrograph simultaneously. The derivation of the formulas is provided in detail in this paper, and an example is also included to illustrate how to apply the new formulas to the flood design problems in small basins.     

AN OBJECTWE TEST FOR HYDROLOGIC SCALE1

ABSTRACT: Improving the reliability of parametric hydrologic models (sometimes called cenceptual rainfall-runoff models) in the continuous simulation of runoff from ungaged catchments has been frustrated by difficulties in estimating model parameters from catchment characteristics. An underlying problem is that these models use parameters to represent catchments as a whole, whereas data on catchment characteristics are collected at multiple field locations and are difficult to transform into one measure of collective impact. Subdividing the catchment and calibrating a stochastic parametric model to estimate distributions for the parameters that covered the range of observed streamflow values was found to improve the simulations. This paper presents an optimization of the amount of subdivision to use in simulation with a version of the Stanford Watershed Model using available climatological data. The calibration process assumes that catchment heterogeneity introduces errors that can be reduced by calibrating parameters as spatial distributions rather than single values. Calibrations for three diverse small gaged catchments located in California and in Virginia found the optimal number of subdivisions to range from 4 to 25 and the optimal scale to range from 0.3 to 2.1 mi².     

WATER LOSS EQUATIONS AND COEFFICIENTS TO ESTIMATE RUNOFF FROM RAINFALL1

Infiltration models are based on physical characteristics of the soil and initial soil moisture. For a given soil it is based on the initial soil moisture distribution. A computer simulation model for flood runoff systems (FH-Model) was used to analyze 39 sets of rainfall-runoff data on four small watersheds ranging in size from 17 to 342 square kilometers located in the Yamaska River basin in Quebec. From these analyses, parameters and coefficients have been determined for a water loss (infiltration) equation. A method for determining the loss parameters, using a nonlinear least square curve fitting technique, is presented. Expressions were made to relate the loss parameters to antecedent precipitation. The equations were tested on 11 storm rainfall and runoff events on a watershed located in the same region and close agreements were found.     

ROLE OF CARBONATE ROCKS IN MODIFYING FLOOD FLOW BEHAVIOR1

ABSTRACT. The interrelationships between the runoff characteristics of watersheds (expressed as the mean annual flood), standard basin parameters (area, drainage properties, and relief), and the parameters which describe the solutional modification of the basins (carbonate rock fractions, sinkhole development, and measures of internal drainage) were used to group 62 carbonate watersheds. Simple binary correlations were obtained by direct plotting of the data. This was followed by multivariate analyses: factor and cluster analyses. Following the cluster analysis, which separated the basins into three groups, the variance within each group was examined again by binary correlations and by factor analysis. Prediction equations for those basins underlain by dolomite rock [QBAR = 12.4 TOT^1.01] and for those basins underlain by carbonate rock with very little surface expression [QBAR = 43.5 TOT^0.87] were proposed. Basins underlain by karstic limestone had a large amount of variance within the data set; therefore no prediction equation could be obtained. (QBAR = mean annual flood, cfs; TOT = total length of all blue lines shown on topographic maps, miles.)     

SPATIAL AND TEMPORAL VARIABILITY IN SUSPENDED SEDIMENT YIELDS OF NORTH ISLAND BASINS,NEW ZEALAND1

ABSTRACT: Specific annual suspended sediment yields and their standard deviations are presented for 47 basins of North Island, New Zealand. Most of the variance in yields is explained by catchment mean rainfall. Rivers with similar flow range have similar suspended sediment concentration ratings, independent of differing watershed lithology and regolith, except for six basins having an abundance of soft fine sediments. Prediction equations for yield and its standard deviation are derived for four essentially arbitrary regions. AU feature rainfall as the independent variable. Differences between regions may owe to variations in intensity, frequency, and duration patterns of storms and, in one area, to bed material size as well. The temporal distribution of annual yields from a basin m be modeled by a two-parameter lognormal function: the prediction equations above may be used to evaluate this function at a site for which suspended sediment data are unavailable.     

Examining Flood Frequency Distributions in the Midwest U.S.1

Villarini, Gabriele, James A. Smith, Mary Lynn Baeck, and Witold F. Krajewski, 2011. Examining Flood Frequency Distributions in the Midwest U.S. Journal of the American Water Resources Association (JAWRA) 47(3):447‐463. DOI: 10.1111/j.1752‐1688.2011.00540.x Abstract: Annual maximum peak discharge time series from 196 stream gage stations with a record of at least 75 years from the Midwest United States is examined to study flood peak distributions from a regional point of view. The focus of this study is to evaluate: (1) “mixtures” of flood peak distributions, (2) upper tail and scaling properties of the flood peak distributions, and (3) presence of temporal nonstationarities in the flood peak records. Warm season convective systems are responsible for some of the largest floods in the area, in particular in Nebraska, Kansas, and Iowa. Spring events associated with snowmelt and rain‐on‐snow are common in the northern part of the study domain. Nonparametric tests are used to investigate the presence of abrupt and slowly varying changes. Change‐points rather than monotonic trends are responsible for most violations of the stationarity assumption. The abrupt changes in flood peaks can be associated with anthropogenic changes, such as changes in land use/land cover, agricultural practice, and construction of dams. The trend analyses do not suggest an increase in the flood peak distribution due to anthropogenic climate change. Examination of the upper tail and scaling properties of the flood peak distributions are examined by means of the location, scale, and shape parameters of the Generalized Extreme Value distribution.     

SENSITWITY ANALYSIS OF SIMULATED CONTAMINATED SEDIMENT TRANSPORT1

ABSTRACT: A simulation analysis of contaminated sediment transport involves model selection, data collection, model calibration and verification, and evaluation of uncertainty in the results. Sensitivity analyses provide information to address these issues at several stages of the investigation. A sensitivity analysis of simulated contaminated sediment transport is used to identify the most sensitive output variables and the parameters most responsible for the output variable sensitivity. The output variables included are streamflow and the flux of sediment and Cs¹³⁷. The sensitivities of these variables are measured at the field and intermediate scales, for flood and normal flow conditions, using the HSPF computer model. A sensitivity index was used to summarize and compare the results of a large number of output variables and parameters. An extensive database was developed to calibrate the model and conduct the sensitivity analysis on a 6.2 mi² catchment in eastern Tennessee. The fluxes of sediment and Cs¹³⁷ were more sensitive than streamflow to changes in parameters for both flood and normal flow conditions. The relative significance of specific parameters on output variable sensitivity varied according to the type of flow condition and the location in the catchment. An implications section illustrates how sensitivity analysis results can help with model selection, planning data collection, calibration, and uncertainty analysis.     

Mechanistic Simulation of Tree Effects in an Urban Water Balance Model1

Abstract: A semidistributed, physical‐based Urban Forest Effects – Hydrology (UFORE‐Hydro) model was created to simulate and study tree effects on urban hydrology and guide management of urban runoff at the catchment scale. The model simulates hydrological processes of precipitation, interception, evaporation, infiltration, and runoff using data inputs of weather, elevation, and land cover along with nine channel, soil, and vegetation parameters. Weather data are pre‐processed by UFORE using Penman‐Monteith equations to provide potential evaporation terms for open water and vegetation. Canopy interception algorithms modified established routines to better account for variable density urban trees, short vegetation, and seasonal growth phenology. Actual evaporation algorithms allocate potential energy between leaf surface storage and transpiration from soil storage. Infiltration algorithms use a variable rain rate Green‐Ampt formulation and handle both infiltration excess and saturation excess ponding and runoff. Stream discharge is the sum of surface runoff and TOPMODEL‐based subsurface flow equations. Automated calibration routines that use observed discharge has been coupled to the model. Once calibrated, the model can examine how alternative tree management schemes impact urban runoff. UFORE‐Hydro model testing in the urban Dead Run catchment of Baltimore, Maryland, illustrated how trees significantly reduce runoff for low intensity and short duration precipitation events.     

A HYDROCLIMATIC APPLICATION STRATEGY FOR THE POISSON PARTIAL DURATION MODEL1

ABSTRACT: A strategy for formulating and testing the Poisson partial duration extreme value model is presented. The procedure is demonstrated using recorded Streamflow series from a humid subtropical region of the southern United States. The observed data series are partitioned by climatic causes and tested for both the Poisson assumption and the validity of the exponential as marginal distributions. Several statistical tests are utilized in making these determinations. Some important aspects of the model as applied to humid climates are demonstrated. It was found that a majority of Streamflow series could be represented by the model and that significant differences do exist between the arrival structures of floods resulting from different climatic mechanisms. However, these differences generally do not exist in the distribution of the flood magnitudes. In addition, it is possible that model validity is restricted by drainage basin size.     

Rehabilitation and Flood Management Planning in a Steep,Boulder-Bedded Stream

This study demonstrates the integration of rehabilitation and flood management planning in a steep, boulder-bedded stream in a coastal urban catchment on the South Island of New Zealand. The Water of Leith, the primary stream flowing through the city of Dunedin, is used as a case study. The catchment is steep, with a short time of concentration and rapid hydrologic response, and the lower stream reaches are highly channelized with floodplain encroachment, a high potential for debris flows, significant flood risks, and severely degraded aquatic habitat. Because the objectives for rehabilitation and flood management in urban catchments are often conflicting, a number of types of analyses at both the catchment and the reach scales and careful planning with stakeholder consultation were needed for successful rehabilitation efforts. This included modeling and analysis of catchment hydrology, fluvial geomorphologic assessment, analysis of water quality and aquatic ecology, hydraulic modeling and flood risk evaluation, detailed feasibility studies, and preliminary design to optimize multiple rehabilitation and flood management objectives. The study showed that all of these analyses were needed for integrated rehabilitation and flood management and that some incremental improvements in stream ecological health, aesthetics, and public recreational opportunities could be achieved in this challenging environment. These methods should be considered in a range of types of stream rehabilitation projects.     

EFFECT OF RAINFALL EXCESS CALCULATIONS ON MODELED HYDROGRAPH ACCURACY AND UNIT-HYDROGRAPH PARAMETERS1

ABSTRACT: Two methods of computing rainfall excess in the U.S. Army Corps of Engineers’flood hydrograph package (HEC-1), the Initial and Uniform method and the Exponential method, are compared to evaluate the effects on modeled hydrograph accuracy. Two computed unit-hydrograph parameters, time of concentration and storage coefficient, were also compared. Rainfall and runoff data from 209 storms in 32 gaged basins in Illinois were used to calibrate the HEC-1 model. Three hydrograph characteristics - sum of incremental flows, peak discharge, and time of peak discharge - were used to evaluate modeled hydrograph accuracy. Mean percent error for each basin and hydrograph characteristic was computed. An evaluation of the mean errors indicates that, although some bias in modeled hydrograph accuracy is evident, rainfall excess computed using either method results in a computed hydrograph accuracy that is within generally accepted limits. Application of a linear-regression model shows no significant differences in computed values of unit-hydrograph parameters.