ESTIMATING THE Q100 IN BRITISH COLUMBIA: A PRACTICAL PROBLEM IN FOREST HYDROLOGY1

ABSTRACT. Estimates of peak flows, with specified return periods, are needed in practice for the design of works that affect streams in forested areas. In the province of British Columbia (B.C.), Canada, the new Forest Practices Code specifies the 100-year instantaneous peak flow (Q100) for the design of bridges and culverts for stream crossings under forest roads; and many practitioners are engaged in making such estimates. The state of the art is still quite primitive, very similar to the state of urban hydrology 30 years ago, when popular estimating techniques were used with little consideration given to their applicability. Urban hydrology then evolved on a much more scientific basis, such that within about a 10-year period, standard approaches to design were developed. Forest hydrology should follow the same pattern, at least as far as estimating design flows is concerned. Popular present day design procedures include the rational method and other empirical approaches based on rainfall data, as use of the standard flood frequency approach is limited by the paucity of relevant flow data. Estimating procedures based on peak streamflow measurements and statistics are likely to evolve, and these will include distinctions for rain, snowmelt, and rain on snow floods. Guidelines will also be developed for selecting and applying appropriate procedures for particular areas.     

Examining Modeling Approaches for the Rainfall‐Runoff Process in Wildfire‐Affected Watersheds: Using San Dimas Experimental Forest

Wildfire can significantly change watershed hydrological processes resulting in increased risks for flooding, erosion, and debris flow. The goal of this study was to evaluate the predictive capability of hydrological models in estimating post‐fire runoff using data from the San Dimas Experimental Forest (SDEF), San Dimas, California. Four methods were chosen representing different types of post‐fire runoff prediction methods, including a Rule of Thumb, Modified Rational Method (MODRAT), HEC‐HMS Curve Number, and KINematic Runoff and EROSion Model 2 (KINEROS2). Results showed that simple, empirical peak flow models performed acceptably if calibrated correctly. However, these models do not reflect hydrological mechanisms and may not be applicable for predictions outside the area where they were calibrated. For pre‐fire conditions, the Curve Number approach implemented in HEC‐HMS provided more accurate results than KINEROS2, whereas for post‐fire conditions, the opposite was observed. Such a trend may imply fundamental changes from pre‐ to post‐fire hydrology. Analysis suggests that the runoff generation mechanism in the watershed may have temporarily changed due to fire effects from saturation‐excess runoff or subsurface storm dominated complex mechanisms to an infiltration‐excess dominated mechanism. Infiltration modeling using the Hydrus‐1D model supports this inference. Results of this study indicate that physically‐based approaches may better reflect this trend and have the potential to provide consistent and satisfactory prediction.     

MODIFICATION OF NEILL'S EQUATION FOR STORM‐TIDE DESIGN FLOW ANALYSIS1

ABSTRACT: Design of bridges spanning tidal estuaries or bays requires an estimate of peak tidal flow. One common approach to estimating these flows (Neill's method) uses a first‐order approximation of uniform water surface rise in the water body. For larger water bodies, the assumptions of this method are decreasingly valid. This study develops a simple modification that accounts for the spatial variability in the response of tidal waterways to storm surge flows. The peak tidal flow predicted by Neill's equation is compared to the peak flow determined by numerical simulation of estuaries with simple geometries, ranging from 1 to 25 km in length, using the U.S. Army Corps of Engineers one‐dimensional unsteady flow model, UNET. Results indicate that, under certain conditions, it may be appropriate to apply a correction factor to the peak discharge and peak velocity predicted by Neill's method. An algorithm, developed by nonlinear regression, is presented for computing correction factors based on estuary length, shape, mean depth, and storm‐tide characteristics. The results should permit the design of more reliable, cost‐effective structures by providing more realistic estimates of the potential for bridge scour in tidal waterways, especially when a full solution of the unsteady flow equations is impractical.     

UNCERTAINTIES IN ESTIMATING THE WATER BALANCE OF LAKES1

ABSTRACT: Evaluation of hydrologic methodology used in a number of water balance studies of lakes in the United States shows that most of these studies calculate one or more terms of the budget as the residual. A literature review was made of studies in which the primary purpose was error analysis of hydrologic measurement and interpretation. Estimates of precipitation can have a wide range of error, depending on the gage placement, gage spacing, and areal averaging technique. Errors in measurement of individual storms can be as high as 75 percent. Errors in short term averages are commonly in the 15-30 percent range, but decrease to about 5 percent or less for annual estimates. Errors in estimates of evaporation can also vary widely depending on instrumentation and methodology. The energy budget is the most accurate method of calculating evaporation; errors are in the 10–15 percent range. If pans are used that are located a distance from the lake of interest, errors can be considerable. Annual pan-to-lake coefficients should not be used for monthly estimates of evaporation because they differ from the commonly used coefficient of 0.7 by more than 100 percent. Errors in estimates of stream discharge are often considered to be within 5 percent. If the measuring section, type of flow profile, and other considerations, such as stage discharge relationship, are less than ideal errors in estimates of stream discharge can be considerably greater than 5 percent. Errors in estimating overland (nonchannelized) flow have not been evaluated, and in most lake studies this component is not mentioned. Comparison of several lake water balances in which the risdual consists solely of errors in measurement, shows that such a residual, if interpreted as ground water, can differ from an independent estimate of ground water by more than 100 percent.     

ADAPTATION OF FLOOD PEAKS AND DESIGN HYDROGRAPHS FROM GAGED TO NEARBY UNGAGED WATERSHEDS1

ABSTRACT: Equations were developed to transform peak flows and to adapt design hydrographs and unit hydrographs from gaged watersheds to ungaged watersheds with similar hydrologic characteristics. Dimensional analysis was used to develop adjustment equations for peak flow and time base, and these equations were reinforced with results from regional flood frequency research. The authors believe that the use of these transformation equations should yield more reliable flood peak values and hydrogrphs than the common use of empirical flood estimating curves or equations.     

The Importance of the Chosen Technique to Estimate Diffuse Solar Radiation by Means of Regression

The Ordinary Least Squares (OLS) is one of the widely used methods which is used for estimating the diffuse solar radiation. However, in order to use the OLS method in the estimation, the dataset must provide certain assumptions. In this study, alternative robust methods have been described and they were compared with the OLS method, which is used for estimating diffuse radiation frequently in an application. At the end of the analysis, the R² value obtained by the OLS method is less than the values obtained by M regression models. In other words, the explanation of the dependent value is weak when the OLS method is used. Finally, it can be said that the most appropriate method is Andrews for estimating the diffuse solar radiation.     

PHYSICAL HYDROLOGY AND THE EFFECTS OF FOREST HARVESTING IN THE PACIFIC NORTHWEST: A REVIEW1

The Pacific Northwest encompasses a range of hydrologic regimes that can be broadly characterized as either coastal (where rain and rain on snow are dominant) or interior (where snowmelt is dominant). Forest harvesting generally increases the fraction of precipitation that is available to become streamflow, increases rates of snowmelt, and modifies the runoff pathways by which water flows to the stream channel. Harvesting may potentially decrease the magnitude of hyporheic exchange flow through increases in fine sediment and clogging of bed materials and through changes in channel morphology, although the ecological consequences of these changes are unclear. In small headwater catchments, forest harvesting generally increases annual runoff and peak flows and reduces the severity of low flows, but exceptions have been observed for each effect. Low flows appear to be more sensitive to transpiration from vegetation in the riparian zone than in the rest of the catchment. Although it appears that harvesting increased only the more frequent, geomorphically benign peak flows in several studies, in others the treatment effect increased with return period. Recovery to pre‐harvest conditions appeared to occur within about 10 to 20 years in some coastal catchments but may take many decades in mountainous, snow dominated catchments.     

SOME TECHNIQUES FOR USING FREQUENCY ANALYSIS AND REALTIME DATA TO INTERPRET FLOOD POTENTIAL DATA1

ABSTRACT: Flood potential data can be effectively interpreted if simple frequency analysis concepts are used to explain the significance of flood potential. Instead of simply presenting data as a quantitative amount or as a percentage of the average condition, predictions can be discussed in terms of their probabilities of exceedance, or return periods. Criteria are presented for evaluating the significance of various return periods. Frequency interpretations are applied to snow course data, peak flow forecasts, and streamflow volume forecasts in northern Utah to illustrate these concepts. In addition, access to realtime data allows tracking of snowmelt progression and identification of any deviations from the forecast flood potential situation. Several data elements, including snowpack, streamfiow volume and peak, and realtime data are jointly evaluated to assess potential hazard and probable risk.     

THE DETERMINATION OF PRIMARY PRODUCTION IN A STREAM USING AN EXACT SOLUTION TO THE OXYGEN BALANCE EQUATION1

ABSTRACT. A promising technique for recognition of early stages of cultural eutrophication relies on determining production and respiration in streams. The most successful and most widely used method of estimating production of a segment of a stream is the upstream-downstream, diurnal curve method introduced by Odum [1956]. This technique is equivalent to obtaining an approximate solution to the oxygen balance equation. We report here an exact solution of the balance equation as a method for calculating primary production. Data presented by Owens [1966] are analyzed; effects of depth and oxygen saturation are studied. A major advantage of the method described here is that continuous temporal variation of net production may be rigorously handled. The method is shown to be well suited to our ultimate goal of studying energy budgets of streams, and thereby the eutrophication process.     

APPLICATION OF ADVANCES IN FLOOD FREQUENCY ANALYSIS1

ABSTRACT: Flood frequency analyses are frequently being made using widely available computer programs. Serious errors can result from blind acceptance of such results. Visual interpretation of observed flood series can be used for evaluation on frequency paper with compatible scales. Such frequency papers are presented in the paper. In ephemeral streams, more infrequent floods may constitute a separate set from the more frequent floods because (a) runoff producing storms cover only a portion of the contributing area, (b) transmission losses in the normally dry streambed may reduce the peak flow, and (c) some runoff may be stored in stock water ponds which therefore leads to partial area runoff. The Cunnane plotting position used in this paper is superior to the more widely used Weibull equation, having a mathematically sound basis for locating observed floods on an assumed probability.     

A METHOD FOR ESTIMATING PEAK WATER DEMAND OF MULTIFAMILY RESIDENCES1

ABSTRACT: The importance of estimating peak water demands for determining the dimensions of pipe size and meters which provide household water to multifamily residences is examined. Several methods utilized in North America and Europe are examined. The analysis makes clear the necessity of studying the peak water demand through statistics based on local data concerning multifamily residences. For different periods of return, the peak demand of a given apartment building is related to its size (the number of apartments) in order to compare the results obtained with existing formula. By use of Ridge-regression technique, the relationship between peak water demand, and building size (number of apartments) and available pressure is established. It can be concluded that peak demand can be estimated with the square root of number of apartment units in the building and the cube root of water pressure.     

INSTANTANEOUS PEAK FLOW ESTIMATION PROCEDURES FOR NEWFOUNDLAND STREAMS1

ABSTRACT: A procedure of estimating instantaneous flood flows for various return periods on the Island of Newfoundland is presented. The procedure is based on annual maximum instantaneous flows rather than annual maximum daily-mean flows, as the latter requires the conversion of estimated daily-mean flows into instantaneous flows. Regression equations were developed for each of three homogeneous regions for the desired return periods. The flood flow estimation capability of the presented procedure is demonstrated to be better than any other currently available procedure on the Island.     

CALIBRATION OF SNYDER COEFFICIENTS FOR PENNSYLVANIA1

ABSTRACT Unit hydrograph theory is one of the most widely used techniques to predict surface runoff. The present study is concerned with the Snyder unit hydrograph and the calibration of the Snyder coefficients for Pennsylvania. Twenty-seven study basins were selected, located randomly across the state. With the rainfall and runoff recorded for several events for each basin (more than 500 events were analyzed) unit hydrographs were calculated and the Snyder coefficients determined. A map of the coefficients was drawn to illustrate the variability in the coefficients and two equations using multiple regression theory were developed. The unexplained variability of the coefficients suggests that upper and lower bounds on the peak flow might be placed on storm hydrographs developed for ungaged watersheds.     

Peak Runoff Model Comparison on Central Illinois Watersheds1

ABSTRACT: Four peak runoff rate models were tested with 183 gage years of record to determine the model most applicable to small watersheds of mild topography in east-central Illinois. The Cypress Creek, Rational, Chow, and SCS peak runoff models were evaluated for their performance. Statistical analyses indicated the Soil Conservation Service model was most appropriate for the watersheds tested.     

A MULTICOMPONENT MODEL OF PHOSPHORUS DYNAMICS IN RESERVOIRS1

The models available for simulating phosphorus dynamics and trophic state in impoundments vary widely. The simpler empirically derived phosphorus models tend to be appropriate for long-term, steady or near steady state analyses. The more complex ecosystem models, because of computational expense and the importance of input parameter uncertainty, are impractical for very long-term simulation and most applicable for time-variable water quality simulations generally of short to intermediate time frames. An improved model for time variable, long-term simulation of trophic state in reservoirs with fluctuating inflow and outflow rates and volume is needed. Such a model is developed in this paper representing the phosphorus cycle in two-layer (i.e., epilimnion and hypolimnion) reservoirs. The model is designed to simulate seasonally varying reservoir water quality and eutrophication potential by using the phosphorus state variable as the water quality indicator. Long-term simulations with fluctuating volumes and variable influent and effluent flow rates are feasible and practical. The model utility is demonstrated through application to a pumped storage reservoir characteristic of these conditions.     

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.     

MODELING DEVELOPED COASTAL WATERSHEDS WITH THE AGRICULTURAL NON-POINT SOURCE MODEL1

ABSTRACT: Many coastal states are facing increasing urban growth along their coast lines. The growth has caused urban non-point source nitrogen runoff to be a major contributor to coastal and estuarine enrichment. Water resource managers are responsible for evaluating the impacts from point and non-point sources in developed watersheds and developing strategies to manage future growth. Non-point source models provide an effective approach to these management challenges. The Agricultural Non-Point Source Model (AGNPS) permits the incorporation of important spatial information (soils, landuse, topography, hydrology) in simulating surface hydrology and nitrogen non-point source runoff. The AGNPS model was adapted for developed coastal watersheds by deriving urban coefficients that reflect urban landuse classes and the amount of impervious surface area. Popperdam Creek watershed was used for model parameter development and model calibration. Four additional watersheds were simulated to validate the model. The model predictions of the peak flow and total nitrogen concentrations were close to the field measurements for the five sub-basins simulated. Measured peak flow varied by 30 fold among the sub-basins. The average simulated peak flow was within 14 percent of the average measured peak flow. Measured total nitrogen loads varied over an order of magnitude among the sub-basins yet error between the measured and simulated loads for a given sub-basin averaged 5 percent. The AGNPS model provided better estimates of nitrogen loads than widely used regression methods. The spatial distribution of important watershed characteristics influenced the impacts of urban landuse and projecting future residential expansion on runoff, sediment and nitrogen yields. The AGNPS model provides a useful tool to incorporate these characteristics, evaluate their importance, and evaluate fieldscale to watershed-scale urban impacts.     

Comparisons and critical assessment of global and diffuse solar irradiation estimation methodologies

Defining better methodologies of accurate predictions of the amount of monthly mean daily global and diffuse solar irradiation exposed is of utmost importance in order to determine the potential for utilizing the solar energy. This study compares and discusses the main methodologies, databases, and software that are used in estimating the solar irradiation to be used for the short- and long-term performances and feasibilities of solar energy systems, especially photovoltaic power plants in Turkey, and addresses the best one to be used to make the most accurate estimations. The comparisons are carried out between the recent methodologies developed by the authors, some models taken from the related literature that are concluded to be better, a widely used database, namely Meteonorm, and a widely used software, namely EU PVGIS. The reference data to develop the methodologies and to make comparison are provided from the State Meteorological Service of Turkey, which is the responsible body in Turkey to make measurements of solar irradiation. The comparisons are based on monthly mean daily values of global and diffuse solar irradiation and are carried out by statistical errors: mean bias error and root mean square error. The results showed that the methodology developed by the authors has shown better performances in estimating the monthly mean daily global and diffuse solar irradiation amount for Turkey.     

Enhancing Hyrdological Simulation Program – FORTRAN Model Channel Hydraulic Representation1

Abstract: The hydrological simulation program – FORTRAN (HSPF) is a comprehensive watershed model that employs depth‐area‐volume‐flow relationships known as the hydraulic function table (FTABLE) to represent the hydraulic characteristics of stream channel cross‐sections and reservoirs. An accurate FTABLE determination for a stream cross‐section site requires an accurate determination of mean flow depth, mean flow width, roughness coefficient, longitudinal bed slope, and length of stream reach. A method that uses regional regression equations to estimate mean flow depth, mean flow width, and roughness coefficient is presented herein. FTABLES generated by the proposed method (Alternative Method) and FTABLES generated by Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) were compared. As a result, the Alternative Method was judged to be an enhancement over the BASINS method. First, the Alternative Method employs a spatially variable roughness coefficient, whereas BASINS employs an arbitrarily selected spatially uniform roughness coefficient. Second, the Alternative Method uses mean flow width and mean flow depth estimated from regional regression equations whereas BASINS uses mean flow width and depth extracted from the National Hydrography Dataset (NHD). Third, the Alternative Method offers an option to use separate roughness coefficients for the in‐channel and floodplain sections of compound channels. Fourth, the Alternative Method has higher resolution in the sense that area, volume, and flow data are calculated at smaller depth intervals than the BASINS method. To test whether the Alternative Method enhances channel hydraulic representation over the BASINS method, comparisons of observed and simulated streamflow, flow velocity, and suspended sediment were made for four test watersheds. These comparisons revealed that the method used to estimate the FTABLE has little influence on hydrologic calibration, but greatly influences hydraulic and suspended sediment calibration. The hydrologic calibration results showed that observed versus simulated daily streamflow comparisons had Nash‐Sutcliffe efficiencies ranging from 0.50 to 0.61 and monthly comparisons had efficiencies ranging from 0.61 to 0.84. Comparisons of observed and simulated suspended sediments concentrations had model efficiencies ranging from 0.48 to 0.56 for the daily, and 0.28 to 0.70 for the monthly comparisons. The overall results of the hydrological, hydraulic, and suspended sediment concentration comparisons show that the Alternative Method yielded a relatively more accurate FTABLE than the BASINS method. This study concludes that hydraulic calibration enhances suspended sediment simulation performance, but even greater improvement in suspended sediment calibration can be achieved when hydrological simulation performance is improved. Any improvements in hydrological simulation performance are subject to improvements in the temporal and spatial representation of the precipitation data.