RELATING TRENDS OF PRINCIPAL COMPONENTS TO TRENDS OF WATER-QUALITY CONSTITUENTS1

ABSTRACT: The seasonal Kendall test is used for detecting water-quality trend or lack of trend for monthly data of 15 water-quality constituents at 15 sampling stations in the Arkansas River, the Neosho River, and the Verdigris River basins. Trends of individual constituents and the trends of the first four principal components for the correlation matrix of water-quality data at each station are determined, and the relationships between the trends of constituents and the trends of principal components are established. Using the principal components not only reduces the high dimensionality of the original data to a few principal components, but also presents an overall picture of water-quality trend of these river basins.     

REGIONAL FREQUENCY ANALYSIS OF HYDROLOGIC MULTIYEAR DROUGHTS1

This paper is concerned with regional frequency analysis of hydrologic multiyear droughts. A drought event is defined by three parameters: severity, duration, and magnitude. A method is proposed here to standardize drought severities with a duration adjustment to enable comparison among drought events. For purposes of a regional study, the index drought method is selected and applied to standardized droughts to give a regional frequency curve. However, the recurrence intervals of the drought events obtained from index drought method are limited to the historic period of record. Therefore, by taking advantage of random variations of droughts in both time and space, a multivariate simulation model is used to estimate exceedence probabilities associated with regional drought maxima. This method, named the regional extreme drought method, is capable of generating a series of drought events which, although they have not occurred historically, are more severe than historic events. By combining the results of the index drought method and regional extreme drought analysis, a regional drought probability graph is constructed which ranges from severe droughts to more frequent droughts. This procedure is applied to the mean annual flow records of streams located in the San Joaquin Valley of California, and drought-severity-frequency plots are prepared for 1-year, 2-year, and 3-year durations.     

APPLICATION OF PARAMETER OPTIMIZATION METHOD FOR CALIBRATING TANK MODEL1

ABSTRACT: Many automatic calibration processes have been proposed to efficiently calibrate the 16 parameters involved in the four‐layered tank model. The Multistart Powell and Stuffed Complex Evolution (SCE) methods are considered the best two procedures. Two rainfall events were designed to compare the performance and efficiency of these two methods. The first rainfall event is short term and the second designed for long term rainfall data collection. Both rainfall events include a lengthy no‐rainfall period. Two sets of upper and lower values for the search range were selected for the numerical tests. The results show that the Multistart Powell and SCE methods are able to obtain the true values for the 16 parameters with a sufficiently long no‐rainfall period after a rainfall event. In addition, by using two selected objective functions, one based on root mean square error and one based on root mean square relative error criteria, it is found that the no‐rainfall period lengths necessary to obtain the converged true values for the 16 parameters are roughly the same. The SCE method provides a more efficient search based on an appropriate preliminary search range. The Multistart Powell method, on the other hand, leads to more accurate search results when there is no suitable search range selected based on the parameter calibration experience.     

A COMPARISON OF STATISTICAL TESTS ON THE EXTREME-VALUE TYPE-1 DISTRIBUTION1

ABSTRACT: The effects of the selections of plotting position formulae and class division schemes on goodness of fit tests are investigated for the extreme-value type-1 distribution using annual flood and annual maximum daily rainfall data. It was found that the plotting position formulae have a minor influence on the tests which involve their use, while the class division schemes may have a pronounced effect on the X²-test. The study also recommended the maximum likelihood method for fitting purposes and a new test which is independent of the afore-mentioned selections for judging the goodness of fit.     

DETECTION OF FRESH WATER AQUIFERS IN THE GLACIAL DEPOSITS OF NORTHWESTERN MISSOURI BY GEOELECTRICAL METHODS1

Geoelectrical investigations in Grundy County of northwestern Missouri, where the groundwater resources of the glacial deposits have already been examined through an extensive drilling program by the Missouri Geological Survey and Water Resources, indicate that water-bearing gravel deposits can be distinguished from glacial deposits containing appreciable amounts of clay and limited amounts of water. The Schlumberger method used for the geoelectric depth soundings in the vicinity of the Survey's drillholes demonstrates the exploratory usefulness of the method in that it can partly replace the more expensive procedure of drilling. The method also provides improved interpretation between drillholes. Results of the investigation show that, in the area, clay has a resistivity below 20in, that the fresh water-bearing gravel at the bottom of the buried glacial stream channels has a resisitivity of 40 to 50fim, and that the near surface glacial gravel deposits have a resistivity above lOOfim. Interpretation of the depth soundings and the conductivity of water obtained from a local well implies that its water is drawn from the saline water of the bedrock. A recommendation is made for the quality improvement of this particular well.     

APPLICATION OF CRITICAL PATH METHOD TO WATER RESOURCES PLANNING1

ABSTRACT. The task of resource management, in this case water resources, is rapidly becoming more complex, particularly because decision making is often contingent upon various prior activities and sets of data. Comprehensive planning is required in order to prevent misallocation of resources or mismanagement in resource development. Such planning involves five general phases which are applicable to any problem faced by society: (1) problem identification; (2) formulation of alternatives; (3) evaluation of alternatives; (4) implementation; (5) review. There have been many attempts to simplify the planning process and effectively carry out these five phases. The experience of the Alberta Water Resources Division has been that the Critical Path Method is one of the most useful tools available today for planning. It involves two basic steps: (1) preparation of a network diagram which (a) identifies all the activities necessary for the completion of a project, (b) correctly sequences these activities, (c) allocates resources; and (2) mathematical computations for scheduling the activities. In other words, this approach breaks a task down into smaller units or activities for easier organization, scheduling, and performance for eventual completion of the project. This paper will illustrate the effectiveness of the Critical Path Method by discussing its application to actual water resources projects.     

DIGITAL SIMULATION OF THE EFFECT OF THERMAL DISCHARGE ON STREAM WATER QUALITY1

A modified transient version of the Streeter-Phelps model along with the energy balance equation is employed to analyze the effects of waste heat discharge from power plants on stream water quality. Analysis is also made to examine the effects of the upstream water quality and stream velocity on the downstream DO concentration level. The resulting coupled nonlinear hyperbolic partial differential equations representing the energy, BOD and DO concentrations are solved by the method of characteristics and simulated on a digital computer. Final numerical results indicate that the allowable quantity of thermal discharge does heavily depend on the upstream quality.     

SUITABILITY ANALYSIS FOR SELECTION OF PIPELINE ROUTES1

ABSTRACT: A brief review of environmental assessment methods is presented as an introduction to a study of environmental management problems in a coastal area of southern Louisiana. The assessment method used in the study is called suitability analysis and involves the use of matrices. The matrix approach is simply a convenient method for organizing all factors representing environmental concerns in a region. The matrices proved to be very useful in defining the present state of an environmental regime and in maintaining an inventory of regional features and potential impacts. The method is described in detail and illustrated by evaluating the potential impacts of a major pipeline construction job in the area. Of major concern is the potential impact on the water resources of the region. Disruption to the natural drainage could be detremental, and affect the suitability of available water for a given use. Limitations of the method are discussed and potential areas for improvement and research are noted. Coupled with an economic assessment, the suitability analysis method should prove beneficial for selecting development sites for industries, residential areas, and other construction activities.     

PARAMETER ESTIMATION OF THE NONLINEAR MUSKINGUM MODEL USING HARMONY SEARCH1

ABSTRACT: A newly developed heuristic algorithm, Harmony Search, is applied to the parameter estimation problem of the nonlinear Muskingum model. Harmony Search found better values of parameters in the nonlinear Muskingum model than five other methods including another heuristic method, genetic algorithm, in terms of SSQ (the sum of the square of the deviations between the observed and routed outflows), SAD (the sum of the absolute value of the deviations between the observed and routed outflows), DPO (deviations of peak of routed and actual flows), and DPOT (deviations of peak time of routed and actual outflow). Harmony Search also has the advantage that it does not require the process of assuming the initial values of design parameters. The sensitivity analysis of Harmony Memory Considering Rate showed that relatively large values of Harmony Memory Considering Rate makes the Harmony Search converge to a better solution.     

EFFECTS OF TIME STEP SIZE IN IMPLICIT DYNAMIC ROUTING1

ABSTRACT The effects of the size of the Δt time step used in the integration of the implicit difference equations of unsteady open-channel flow are determined for numerous typical hydrographs with durations in the order of days or even weeks. Truncation errors related to the size of the Δt time step cause a numerical distortion (dispersion and attenuation) of the computed transient. The magnitude of the distortion is related directly to the size of the time step, the length of channel reach, and the channel resistance and inversely to the time of rise of the hydrograph. The type of finite difference expression which replaces spatial derivatives and non-derivative terms in the partial differential equations of unsteady flow has an important influence on the magnitude of the numerical distortion, as well as the numerical stability of the implicit difference equations. Time step sizes in the range of 3 to 6 hrs generally tend to minimize the combination of required computation time and numerical distortion of transients having a time of rise of the order of several days.     

ESTIMATION OF POTENTIAL REDUCTIONS IN RECREATIONAL BENEFITS DUE TO SEDIMENTATION1

ABSTRACT: A travel cost model is developed to estimate the potential reductions in recreational benefits from sedimentation in Reelfoot Lake in northwestern Tennessee. In addition to the consumer surplus estimates generated by the model, three other aspects of the study were significant. First, the study applied a relatively untested methodology for deriving the opportunity cost of travel time. The study resulted in a value that is less than one-half of the Water Resource Council's “one-third of the wage rate” rule-of-thumb. Second, water quality perceptions were unsuccessfully incorporated into the model as a demand shifter. This raised questions as to the appropriate manner in which perceptions could be included in a travel cost model. Finally, a simple methodology was outlined by which estimates of the recreational value of Reelfoot Lake could be used to suggest how much cost could be justified for soil erosion control on agricultural land surrounding the lake.     

THE CHUBB/BAUMAN METHOD FOR ESTIMATING THE RECREATION POTENTIAL OF RIVERS: A CRITICAL REVIEW1

ABSTRACT: The Chubb/Bauman (Ch/B) method for making quantitative estimates of recreation potential for rivers is based on the 1968/ 69 Leopold method for quantitative assessment of the scenic beauty of rivers. Both use classifications of environmental variables as the database. Unlike the Leopold method, the classifications used in the Ch/B method consistently reflect human preferences. The Ch/B method collects information on 67 variables, and uses a computer program to produce estimates of potential for 16 common recreation activities. This critique evaluates selected concepts and procedures of the Ch/B method partly by comparison with other available methods of recreation resource inventory. It considers the validity and utility of numerical weighting of variables, the use of numbers derived from place in a classification, and the transformation process. The quantitative techniques of the method exhibit serious flaws. Much of the data produced by the method appears to be quantitative but in fact is not, and it does not produce truly quantitative estimates of recreation potential. Classifications of generalized geographic or environmental variables are shown to have serious defects as a basis for evaluation of recreational potential.     

RISK OF EXTREME EVENTS IN A MULTIOBJECTIVE FRAMEWORK1

ABSTRACT: This paper summarizes advances made in risk-based decisionmaking in water resources through use of the partitioned multiobjective risk method (PMRM). (Risk is a measure of the probability and severity of adverse effects.) In the PMRM, the risk of extreme events is differentiated from risk involving less extreme damage severity and is evaluated within a multiobjective framework. Study of the extreme-event risk function f₄(*) has addressed the following issues: methods for calculating f₄(*); the sensitivity of f₄(*) to various parameters, particularly to the partitioning point of the extreme-event range and the selection of probability density functions; insight provided by the statistics of extremes; and the impact of f₄(*) on risk management, for example, in the application of the PMRM to water resources problems. In particular, this paper shares with the reader recent research results on the PMRM, the relationship between the statistics of extremes and the conditional expected value, derived formulas for f₄(*), distribution-free estimates of f₄(*), documented case studies in dam safety, and future research directions.     

EVALUATION OF CURVE NUMBER PROCEDURES TO PREDICT RUNOFF IN GLEAMS1

ABSTRACT: Proper selection of curve number values will improve the capability of the SCS-Curve Number procedure in predicting runoff. Both CREAMS and GLEAMS models use the Smith and Williams (1980) approach of converting CNII (curve number value for average antecedent moisture conditions) into CNI (curve number value for dry antecedent moisture conditions) in calculating the soil retention parameter (S). CREAMS and GLEAMS have been found to under predict runoff because of the internal conversion of CNII to CNI. This study shows modifications of the GLEAMS model using CMI without converting it to CM and it also shows the seasonal curve number approaches with and without converting CNII to CNI. Results indicate that using CNII without internal conversion to CNI provides better runoff and erosion predictions than the original version of GLEAMS and versions with seasonal curve numbers when tested with four years of field data in the Coastal Plain physiographic region of Maryland.     

USE OF THE DELPHI METHOD IN RESOLVING COMPLEX WATER RESOURCES ISSUES1

ABSTRACT: The tri‐state river basins, shared by Georgia, Alabama, and Florida, are being modeled by the U.S. Fish and Wildlife Service and the U.S. Army Corps of Engineers to help facilitate agreement in an acrimonious water dispute among these different state governments. Modeling of such basin reservoir operations requires parallel understanding of several river system components: hydropower production, flood control, municipal and industrial water use, navigation, and reservoir fisheries requirements. The Delphi method, using repetitive surveying of experts, was applied to determine fisheries' water and lake‐level requirements on 25 reservoirs in these interstate basins. The Delphi technique allowed the needs and requirements of fish populations to be brought into the modeling effort on equal footing with other water supply and demand components. When the subject matter is concisely defined and limited, this technique can rapidly assess expert opinion on any natural resource issue, and even move expert opinion toward greater agreement.     

EVALUATION OF METHODS FOR ESTIMATING STREAM WATER QUALITY PARAMETERS IN A TRANSIENT MODEL FROM STOCHASTIC DATA1

ABSTRACT: In this study the estimation of parameters in water quality models represented by linear first order partial differential equations is investigated. Two sets of simulated input-output data, one with input noise and the other with output measurement error, were used. The parameters were estimated by a gradient technique (Bard's method) and a pattern search technique. The results indicate that the output measurement error significantly affects the values of parameter estimates as compared to the noise added to the input. Bard's method consistently gave results with a smaller sum of square value.     

EVALUATION OF THE MAXEY-EAKIN METHOD FOR ESTIMATING RECHARGE TO GROUND-WATER BASINS IN NEVADA1

AESTRACT An evaluation of the Maxey-Eakin method for calculating recharge to ground-water basins in Nevada was performed. The evaluation consisted of comparing Maxey-Eakin estimates with independent estimates of recharge, and analyzing the nature of the differences between the groups of estimates. In the comparison with the Maxey-Eakin estimates, two different groups of independent estimates were used: (1) 40 recharge estimates that were identified from water budgets contained in reports by the Nevada Department of Conservation and Natural Resources and (2) 27 recharge estimates that were identified from previous studies that used models. The results of the comparisons indicate generally good agreement between the Maxey-Eakin estimates and both groups of independent estimates. To quantify this agreement, an analysis was conducted to estimate the uncertainty in the Maxey-Eakin method. The analysis produced an upper bound on the standard deviation of the Maxey-Eakin estimate for a given basin. For the group of 40 water-budget estimates, the upper bound on the standard deviation for an individual basin is 4,800 acre-ft/yr, and the corresponding coefficient of variation of the Maxey-Eakin estimate is no greater than 44 percent. For the group of 27 model estimates, the upper bound on the standard deviation is 4,100 acre-ft/yr, and the corresponding coefficient of variation is no greater than 24 percent.     

OPTIMUM DESIGN OF BORDER IRRIGATION SYSTEM1

ABSTRACT: Border irrigation systems like most of the other surface irrigation systems, do not need too much energy and special equipment. Thus, many farmers have used this system for a long time. On the other hand, design of surface irrigation systems including border irrigation requires many input parameters, and need intensive engineering calculations. The burden of these requirements probably led the users to experimental design of the systems with low efficiencies. However, accurate design and high quality optimization of the border irrigation system that can result in a highly efficient system is possible. In this study, an optimization model for border irrigation system is presented. The Hook Jeev's pattern search optimization method in conjunction with a general mathematical model of border irrigation is used to maximize the irrigation application efficiency. The border irrigation storage and distribution efficiencies, border slope and length, inflow rate, cutoff time, and the Manning's roughness coefficient are selected as constraints. The model is applied to field-measured data. The results show that it is possible to select a suitable combination of the border system's parameters (border's length, inflow rate, and cutoff time) to obtain a maximum application efficiency.     

DEVELOPMENT OF LARGE SCALE GRIDDED RIVER NETWORKS FROM VECTOR STREAM DATA1

ABSTRACT: The Network Tracing Method (NTM) has been developed to determine gridded coarse river networks for modeling large hydrologic systems. For a coarse resolution grid, the NTM determines the downstream cell of each cell and the distance along the actual meandering flow paths between them. Unlike previously developed methods, the NTM uses fine resolution vector river networks as the source of information of the flow patterns rather than digital elevation models. The main advantage of using vector river networks as input is that they capture the hydrologic terrain features better than topographic data do, particularly in areas of low topographic relief. The NTM was applied to South America with a grid resolution of 1 degree by 1 degree and to the globe with a resolution of 2.815 degrees by 2.8125 degrees. Overall, the method captured the flow patterns well. Generated digital river networks and drainage divides showed minor disagreement with those obtained from existing maps, and most of them were consistent with the resolution of the coarse river network. The majority of estimated basin areas were also close to documented values. River lengths calculated with the NTM, however, were consistently underpredicted.     

OPTIMAL CAPACITIES OF WATER SUPPLY RESERVOIRS IN SERIES AND PARALLEL1

ABSTRACT: The planning of water supply reservoirs has traditionally been based on the Rippl or sequent peak analysis which applies to the design of a single reservoir. This paper incorporates the sequent peak method as the central feature in establishing a procedure for determining the sizes of several potential reservoirs located in a system of one or more rivers. Separate algorithms are developed for sites on parallel streams and for sites on the same stream. In both cases the approach is to find the combination of reservoirs which can satisfy a given constant monthly demand at a minimum total construction cost. It is shown that both problems can be cast in the form of a dynamic programming problem. A more complex system is then a combination of reservoirs in parallel and in series. An extension is given if the monthly demand is not constant but each reservoir satisfies a constant fraction of the monthly demand.