HYDRAULIC GEOMETRY AND STREAM CHANNEL BEHAVIOR: A UNCERTAIN LINK1

ABSTRACT: The at-a-station hydraulic geometry of stream channels can serve as a predictor of alluvial stream channel behavior. This geometry is the empirical relations describing changes in water surface width, mean depth, and mean velocity with changing discharge. The exponent values are correlated with channel morphology and behavior such as scour and fill, flow resistance, bank resistance, and competence. Channel behavior and morphology are apparently related, but some causes for effects are uncertain. Several studies, using empirical and theoretical bases, are reviewed here to illustrate the relation between hydraulic geometry and channel behavior, but the relations are not always consistent. Hydraulic geometry variables are easy to measure and readily available, but they do not always reflect what may be more important ones such as turbulence, the velocity distribution profile, and distribution and cohesion of sediment particles. This paper illustrates some of these problems, provides some solutions, and addresses need for more work to better predict stream channel behavior from hydraulic geometry     

INFLUENCES OF A FOREST ON THE HYDRAULIC GEOMETRY OF TWO MOUNTAIN STREAMS1

ABSTRACT The influence of a forest on the formation of steps in two small streams of the Colorado Rocky Mountains was studied. Steps provided by logs fallen across the channel added to flow energy reduction. The streams required additional gravel bars to adjust to slope. Average step length between logs and gravel bars was strongly related to channel gradient and median bed material size. Based on the average number of log steps per 50 feet of channel, an average of 116 percent of gravel bars were added at Fool Creek and 60 percent at Deadhorse Creek. The latter had 52 percent more logs in the channel and therefore required less bed material movement than the former. Although these are “rushing mountain streams,” most flow velocities ranged between 0.5 and 2.5 f.p.s. Exponents of a function relating rate of change of depth or velocity to discharge indicated that dynamic stream equilibrium was attained. Implications for forest management are that sanitation cuts (removal of dead and dying trees) would not be permissible where a stream is in dynamic equilibrium and bed material movement should be minimized.     

USING CURVE NUMBERS TO DETERMINE BASELINE VALUES OF GREEN-AMPT EFFECTIVE HYDRAULIC CONDUCTIVITIES1

ABSTRACT: Since the trend in infiltration modeling is currently toward process-based approaches such as the Green-Ampt equation, more emphasis is being placed on methods of determining appropriate parameters for this approach. The SCS curve number method is an accepted and commonly used empirical approach for estimating surface runoff, and is based on numerous data from a variety of sources. The time and expense of calibrating process-based infiltration parameters to measured data are often prohibitive. This study uses curve number predictions of runoff to develop equations to estimate the “baseline” hydraulic conductivities (K_b) for use in the Green-Ampt equation. Curve number predictions of runoff were made for 43 soils. K_b values in the Water Erosion Prediction Project (WEPP) model were then calibrated so that the annual runoff predicted by WEPP was equal to the curve number predictions. These calibrated values were used to derive an equation that estimated K_b based on the percent sand, percent clay, and cation exchange capacity of the soil. Estimated values of K_b from this equation compared favorably with measured values and values calibrated to measured natural runoff plot data. WEPP predictions of runoff using both optimized and estimated values of K_b were compared to curve number predictions of runoff and the measured values. The WEPP predictions using the optimized values of K_b were the best in terms of both average error and model efficiency. WEPP predictions using estimated values of K_b were shown to be superior to predictions obtained from the curve number method. The runoff predictions all tended to be biased high for small events and low for larger events when compared to the measured data. Confidence intervals for runoff predictions on both an annual and event basis were also developed for the WEPP model.     

INFLUENCES OF INCREASED SAND DELIVERY ON THE MORPHOLOGY OF SAND AND GRAVEL CHANNELS1

ABSTRACT A flume study was conducted to examine (1)changes in the particle-size distribution of sediments in riffles due to the proportion of sand in transport and the total rate of bedload transport at the time the riffle is deposited and (2) the effect of high sand transport rates on the stability of gravel riffles. The median particle size of sediment deposited in the riffle was larger than that of the sediment in transport. Small but significant (a = 0.05) decreases in the median particle size of riffle sediments resulted as the sand-to-gravel ratio. Increased concentrations of sand in transport caused previously stable gravel riffles to undergo scour. These results, in combination with information from other studies, suggest that an alluvial channel with pool-riffle sequences and with sand and gravel beds may respond to an increased delivery of sand by reducing form roughness. Form roughness can be reduced by degrading riffles and filling pools. Subsequent responses may be increases in width-to-depth ratio and slope.     

A CRITICAL APPROACH TO THE CALIBRATION OF A WATERSHED MODEL1

ABSTRACT: A complex watershed-scale water quality simulation model, the Hydrological Simulation Program-FORTRAN (HSPF) model, was calibrated for a 16 km² catchment. The simulation step size was 0.33 hours with predicted and recorded hydrologic flows compared on an annual and monthly basis during a total calibration period of four years. Unguided numerical optimization when applied alone did not yield a model parameter set with acceptable predictive capability; instead, it was necessary to apply a critical process that included sensitivity analysis, numerical optimization, and testing of derived model parameter sets to evaluate their performance for periods other than those for which they were determined. Using this critical calibration process, the model was proven to have significant predictive capability. Numerical optimization is an aid for model calibration, but it must not be used blindly.     

INFILTRATION OF ATRAZINE AND METABOLITES FROM A STREAM TO AN ALLUVIAL AQUIFER1

ABSTRACT: The infiltration of atrazine, deethylatrazine, and deisopropylatrazine from Walnut Creek, a tributary stream, to the alluvial valley aquifer along the South Skunk River in central Iowa occurred where the stream transects the river's flood plain. A preliminary estimate indicated that the infiltration was significant and warrants further investigation. Infiltration was estimated by measuring the loss of stream discharge in Walnut Creek and the concentrations of atrazine and its metabolites deethylatrazine and deisopropylatrazine, in ground water 1 m beneath the streambed. Infiltration was estimated before application of agrichemicals to the fields during a dry period on April 7, 1994, and after application of agrichemicals during a period of small runoff on June 8, 1994. On April 7, the flux of atrazine, deethylatrazine, and deisopropylatrazine from Walnut Creek into the alluvial valley aquifer ranged from less than 10 to 60 (μg/d)/m², whereas on June 8 an increased flux ranged from 270 to 3060 (μg/d)/m². By way of comparison, the calculated fluxes of atrazine beneath Walnut Creek, for these two dates, were two to five orders of magnitude greater than an estimated flux of atrazine to ground water caused by leaching from a field on a per-unit-area basis. Furthermore, the unit-area flux rates of water from Walnut Creek to the alluvial valley aquifer were about three orders of magnitude greater than estimated recharge to the alluvial aquifer from precipitation. The large flux of chemicals from Walnut Creek to the alluvial valley aquifer was due in part to the conductive streambed and rather fast ground water velocities; average vertical hydraulic conductivity through the streambed was calculated as 35 and 90 m/d for the two sampling dates, and estimated ground water velocities ranged from 1 to 5 m/d.     

AN APPROACH TO MARGINAL ECONOMIC ANALYSIS OF HYDROMETRIC DATA COLLECTION1

ABSTRACT: A procedure for computing the benefit/cost ratio of a hydrometric network is outlined. It consists of two steps: firstly, establishing a relationship between hydrometric station density and error in hydrologic parameters and, secondly, relating hydrologic eror to changes in project cost. The procedure was applied to both the whole Canadian hydrometric network and the provincial networks.     

STREAM FLOW IN RELATION TO OHIA FOREST DECLINE ON THE ISLAND OF HAWAII1

ABSTRACT: Dieback of the ohia forest over a large part of the Hilo watershed and adjacent areas has resulted in a severe loss of the overstory crown component of the vegetation. The decline could cause serious damage to the watershed. To evaluate possible changes in conditions in the Hilo area since the decline began, stream flow, water quality, and precipitation data from 1929 to 1980 were collected and analyzed. The limited data available do not indicate that the ohia decline has resulted in a significant change in either annual stream flow or peak stream flow of a stream discharging from the ohia forest in the Hilo area. Suspended sediment production of two streams remains well within accepted levels as does the chemical makeup of components dissolved in the streams.     

INTERPRETATION OF THE DUAL PROGRAM FOR OPTIMAL GROUNDWATER HYDRAULIC CONTROL1

ABSTRACT: Optimization formulations for hydraulic control that take the form of linear programs possess a corresponding dual linear program. The economic and physical interpretations of the dual linear program are examined for formulations in which hydraulic head in groundwater systems is constiained. In each case it is shown that the dual linear program has a physically meaningful interpretation. For a hydraulic gradient control formulation used for remedial analysis it is shown that the dual variable can be interpreted as the remedial benefit due to each gradient control constraint. The dual linear program maximizes the remedial benefit. The value of the dual variable can be used to compute such useful properties as the total remedial benefit of pumping at a specific location. For a formulation that optimizes aquifer yield while constraining drawdown the dual variable can be used to measure the total cost of drawdown capacity consumption per unit of pumping at a specific location. The dual program minimizes the cost of drawdown capacity consumption. By examining the meaning of the dual linear program an alternate statement of the problem under study is revealed. Quantities arising from the dual program add to the value of the optimization approach. Significant new information can be derived from existing linear optimization formulations with minimal additional computational effort.     

USING SOIL TEXTURE TO ESTIMATE SATURATED HYDRAULIC CONDUCTIVITY AND THE IMPACT ON RAINFALL-RUNOFF SIMULATIONS1

ABSTRACT: In this paper a new set of soil texture data is used to estimate the spatial distribution of saturated hydraulic conductivity values for a small rangeland catchment. The estimates of conductivity are used to re-excite and re-evaluate a quasi-physically based rainfall-runoff model. The performance of the model is significantly reduced with conductivity estimates gleaned from soil texture data rather than the infiltration data used in our previous efforts.     

THE ROLE OF MINIMUM VARIANCE THEORY IN DEFINING THE REGIME CHARACTERISTICS OF THE LOWER NAMOI-GWYDIR DRAINAGE BASIN1

ABSTRACT: While minimum variance theory appears to offer an explanation for the hydraulic behaviour and regularity among channel systems it is not entirely successful in predicting the regime of a channel system. In the case of the Namoi-Gwydir river system the hydraulic variables velocity, depth, width, slope, friction, and shear appear to govern the behaviour and, hence regime of the channels. The significance of sediment load in determining regime could not be assessed.     

KALMAN FILTER ESTIMATION AND PREDICTION OF DAILY STREAM FLOWS: II. APPLICATION TO THE POTOMAC RIVER1

Results are reported from an application of the state space formulation and the Kalman filter to real-time forecasting of daily river flows. It is shown that the application of filtering techniques improves the overall forecasting performance of the model. As is true for most hydrologic systems, the model is not completely known. Therefore, the procedures pertaining to on-line parameter and noise statistics estimation, as presented in the first paper, are implemented. The example in this paper shows that these techniques also perform satisfactorily when applied to a real-world situation.     

MODELING AND CONTROL OF THE POLLUTION OF WATER RESOURCES SYSTEMS VIA MULTILEVEL APPROACH1

ABSTRACT. The problem of modeling and control of water pollution is considered. A general mathematical model, where the pollution effluent is discharged directly into the river, into the lake, or into a bypass pipe leading to an advanced Waste Water Treatment (AWT) plant, is developed. The Water Resource System (WRS) under consideration is decomposed into N subsystems. The pollution effluent input vector to each subsystem includes the water quantity and different water characteristics such as BOD, DO, pH, conductivity, temperature, algae, phosphates, nitrates, etc. Treatment cost functions and quality transition functions as well as system model constraints are introduced, where all functions can be nonlinear. A system Lagrangian is formed to incorporate the system constraints and coupling. The Lagrangian is decomposed into N independent subsystems, and a two level optimization methodology is introduced. Each subsystem is independently and separately minimized at the first level assuming known Lagrange multipliers. At the second level, the total Lagrangian is maximized with respect to the Lagrange multipliers using optimal values for effluent inputs from all subsystems obtained from the first level. Economic interpretation on the Lagrange multipliers reveals that they are merely prices imposed by the central authority (second level) for the pollution caused by the subsystems. Advantages of the multilevel approach are discussed.     

A SYSTEMS RESOURCES APPROACH TO CITIZEN PARTICIPATION: THE CASE OF THE CORPS OF ENGINEERS1

ABSTRACT: The thesis of this paper is that the citizen participation process provides necessary, but not sufficient conditions to affect substantive change in federal water resource management agencies' planning and decisionmaking. That is, in its present form, the citizen participation process has been observed to occur outside of the normal decision arenas of federal resource management agencies. The paper reviews concepts of citizen participation and defines some theoretical problems inherent with them. Then, a strategy for the structuring of a citizen participation process is proposed. This strategy is based upon the notion that citizen support for federally sponsored programs are essential if such projects are to be implemented. Therefore, an approach which integrates citizen valves with those held by other institutions and the agency is suggested. In conclusion, the relevancy of actively developing and including citizen input to the water resource planning process is illustrated by a discussion of three cases of the Corps of Engineers and Urban Studies planning process, in different metropolitan regions.     

PHOSPHORUS ASSIMILATION IN A STREAM SYSTEM OF THE LAKE OKEECHOBEE BASIN1

ABSTRACT: The ability to predict how streams and wetlands retain phosphorus (P) is critical to the management of watersheds that contribute nutrients to adjacent aquatic systems such as lakes. Field and laboratory experiments were conducted to determine the P assimilatory capacity of a stream (Otter Creek) in the Taylor Creek/Nubbin Slough Basin located north of Lake Okeechobee, Florida. Dominant soils in this basin are sandy Spodosols; landuse is primarily dairy farms and beef cattle pastures. Estimates of P assimilation show that sediments assimilate approximately 5 percent of the P load. Phosphorus assimilation rates in the stream were estimated using first-order relationships based on the total P concentration of the water column as a function of distance from the primary source. This method assumes minimal lateral inputs. Stream lengths required for one turnover in P assimilation were estimated to be in the range of 3–16 km. Laboratory studies using intact sediment cores indicated a P assimilation rate of 0.025 m day^?1, and equilibrium P concentration of 0.16 ± 0.03 mg L^?1 in the water column. Dissolved P concentration gradients in the sediments showed upward flux of P at water column P concentration of <0.16 mg L^?1. Approximately 56–77 percent of the P assimilated in the above-ground vegetation during active growth was released or translocated within six months of senesence, suggesting short-term storage in above-ground vegetation. Bottom sediments and recalcitrant detrital plant tissue provide for long-term P assimilation in the creek. Although stream sediments have the potential to adsorb P, high flow rate and low contact period between water and sediment limits this process.     

A WATER QUALITY ASSESSMENT OF DEVELOPMENT IN THE SENEGAL RIVER BASIN1

A study was made to determine the impact on water quality due to water resource development in a large river basin in a semi-arid region of West Africa. Mathematical modeling and the examination of case histories were used to project impacts. The impacts associated with changes in water quality were shown to be slight assuming that modern basin and agricultural management practices are adopted. Analytical techniques normally implemented in studies of more highly developed basins are useful for analysis of water quality impacts in relatively undeveloped basins.     

CONTROL OF NATURAL BRINE SPRINGS IN BRAZOS RIVER BASIN PART II: BRINE DISPOSAL1

ABSTRACT: Water quality in the Brazos River of Texas is seriously degraded by natural salt pollution. Two thousand tons/day of total dissolved solids emanate from brine springs and seeps in the Upper Brazos River drainage. Approximately 45 percent of the total salt load comes from a relatively small flow in the Dove Creek area. The companion paper demonstrates that a system of wells pumping brine at a constant rate of about 2 cfs from the near surface aquifer should eliminate the brine springs in this area. In this paper, injection into deep brine aquifers is shown to be a feasible brine disposal alternative. Four brine aquifers were determined from the literature to be possible injection zones. Accurate net aquifer thickness maps were generated in a 23 by 14 mile area centered on the Dove Creek area for three of the aquifers from an interpretation of 41 well logs. Constant injection for a project life of 100 years was simulated using the SWIFT/486 software. Modeling suggests that one well would be sufficient to inject the entire disposal volume into either the Strawn or Ellenburger Formation.     

OPTIMAL WATER ALLOCATION IN THE LAKES BASIN OF NICARAGUA1

ABSTRACT: The water resources of Nicaragua are largely undeveloped and distributed unequally throughout the nation. In addition, there is a significant geographical imbalance between the abundant water supply in the Atlantic Basin and the less abundant supply in the Pacific Basin which accounts for most of the water demand. The Lakes Basin, comprising Lakes Managua and Nicaragua, could be manipulated to solve the imbalance. Specifically, this paper presents a scheme to transfer water from Lake Nicaragua to Lake Managua and, subsequently, direct water from each lake for irrigation and hydroelectric generation. The scheme has been designed to maximize economic benefits and requires environmental impact analysis.     

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.     

ESTIMATION OF SURFACE WATER LAG TIME FROM THE KINEMATIC WAVE EQUATIONS1

ABSTRACT Numerous concepts of surface water lag time have been developed and applied in the past. In this report, hydraulic solutions of a lag time derived by Overton [1970] are presented for several idealistic surfaces using the kinematic wave equations. These surfaces are: (1) a uniform plane; (2) hillslope as a cascade of planes; (3) V-shaped watershed; (4) V-shaped watershed with hill-slopes; (5) converging surface; (6) concave surface. The lag times are shown to be related to roughness, length and catchment slope, and the input rate. These relations may be used immediately in predicting lag time as the parameter in a unit response function. A lag relation has been developed for a nonuniform catchment in terms of the lag of a uniform plane and a convergence factor. A numerical procedure is shown whereby the convergence factor can be evaluated for any nonuniform catchment from observed input and output data.