THE DEPENDENCE OF THE RESIDUAL GRAVITY ON HYDRAULIC CONSTANTS IN GLACIAL DEPOSITS1

ABSTRACT: Buried glacial stream channels contain large and easily accessible groundwater resources. Gravity surveys have been frequently applied for their location. A gravity survey in the geohydrologically explored Wood River Valley Area of southern Rhode Island shows extreme lows of -2 mgals over channel depths of maximal 300 feet. Three gravity profiles were observed in east-west direction across a north-south striking stream channel. The bedrock depth increases rapidly towards the south from 130 to 300 feet. The gravity lows observed across each profile are not related to the bedrock depth but rather to the saturated thickness of the main quifer and its hydraulic transmissivity. Well logs indicate that the large change of bedrock depth is solely due to an increase of till of low permeability. The volume of the glacial outwash, which is the major groundwater resource, changes little underneath the three profiles. The gravity lows apear to be directly related to the density contrast between glacial outwash and till. The response to the hydraulically more pertinent units renews the interest in the gravity method as it may have a potential to estimate yields of hydrologically complex aquifers     

Prediction of the Fate and Transport Processes of Atrazine in a Reservoir

The fate and transport processes of a toxic chemical such as atrazine, an herbicide, in a reservoir are significantly influenced by hydrodynamic regimes of the reservoir. The two-dimensional (2D) laterally-integrated hydrodynamics and mass transport model, CE-QUAL-W2, was enhanced by incorporating a submodel for toxic contaminants and applied to Saylorville Reservoir, Iowa. The submodel describes the physical, chemical, and biological processes and predicts unsteady vertical and longitudinal distributions of a toxic chemical. The simulation results from the enhanced 2D reservoir model were validated by measured temperatures and atrazine concentrations in the reservoir. Although a strong thermal stratification was not identified from both observed and predicted water temperatures, the spatial variation of atrazine concentrations was largely affected by seasonal flow circulation patterns in the reservoir. In particular, the results showed the effect of flow circulation on spatial distribution of atrazine during summer months as the river flow formed an underflow within the reservoir and resulted in greater concentrations near the surface of the reservoir. Atrazine concentrations in the reservoir peaked around the end of May and early June. A good agreement between predicted and observed times and magnitudes of peak concentrations was obtained. The use of time-variable decay rates of atrazine led to more accurate prediction of atrazine concentrations, while the use of a constant half-life (60 days) over the entire period resulted in a 40% overestimation of peak concentrations. The results provide a better understanding of the fate and transport of atrazine in the reservoir and information useful in the development of reservoir operation strategies with respect to timing, amount, and depth of withdrawal.     

Thermal and Electrical Analysis in a Solid Oxide Fuel Cell Stack Using a Curved Profile for the Inlet Flow Rate Along the Stacking Direction

This study analyzes the cell temperature and current density fields of a solid oxide fuel cell stack, in which the inlet flow rate along the stacking direction of the fuel is either a curved profile derived through hydraulics or a progressively decreasing profile. The results show that the temperature field was neglibibly affected by the non-uniform inlet profile of the fuel, but the current density distribution was apparently affected. The difference in current density distribution between the curved profile and the progressively decreasing profile was small in the case of an solid oxide fuel cell (SOFC) with 10 stacks, but became larger in an SOFC with 20 stacks. In the 10-stack SOFC, the difference of current density between the uniform inlet flow and the curved inlet flow was less than 3%, whereas the deviation in the curved profile was less than 12%. However, in a 20-stack SOFC, the difference of current density was over 20%, whereas the deviation in the curved profile was more than 50%. Comparing the results between the curved and the progressively decreasing inlet profile indicates that the current density difference was over 20% in the 20-stack SOFC, but below 1% in the 10-stack SOFC. Therefore, a progressively decreasing inlet profile can replace a curved profile for analyzing an SOFC stack when the number of stacks is less than 10.     

NUMERICAL MODELING OF THERMAL STRATIFICATION IN A RESERVOIR WITH LARGE DISCHARGE-TO-VOLUME RATIO1

ABSTRACT: A numerical model study of thermal stratification in a high discharge-to-volume reservoir is described. Predicted temperature profiles are compared with field data for two different years. The model accurately predicts the date of fall turnover, and predicts degree of stratification and depth of the thermocline within about 20% for both years simulated. A parametric study of stratification mechanics for a high flow reservoir indicated that diffusion was the predominant heat transport mechanism in the hypolimnion, while surface effects dominated the epilimnion. Flow effects for the particular case studied, in which all inflows and outflows occur in the epilimnion, did not significantly affect stratification behavior.     

MIXING INDUCED BY AN INTERNAL HYDRAULIC JUMP1

ABSTRACT The efficiency of an inverted internal hydraulic jump as a mixing and dispersion mechanism in an aquatic environment was examined. The flow considered was a two-dimensional buoyant flow from a shallow channel over a sloping bottom into a deep reservoir. It could be seen that a rapidly varied flow associated with violent turbulent mixing occurred near the point of discharge if specific discharge conditions and downstream controls were met. Downstream from the mixing zone the flow was stably stratified. The main object of the study was to find the conditions under which a turbulent mixing zone occurred and the rate of turbulent entrainment. Energy loss and length of the mixing zone were also investigated. The independent variables were the outlet densimetric Froude number, the density differential between outfall water and receiving water, the relative depths of the upper and lower layers in the stratified flow portion, the total depth, and the slope of the transition. Theoretical calculations had to be confined to a step increase in depth, but experiments in a laboratory flume showed that results obtained with slopes of 23° and 90° were quite similar. Both theory and experiments showed, for example, that dilution (entrainment rates) up to 2:1 (2 parts heated water to one part cold water) can be achieved at very little energy expense and with downstream depths approximately 8 times the outlet depth. Theoretically, any amount of dilution can be obtained, but large depths may be required. Experimental results also indicate that the length of the mixing zone was frequently ten times the value of the outlet densimetric Froude number.     

A PROGRAM TO CALCULATE CHANNEL SCOUR AND FILL1

ABSTRACT: A computer program written in BASIC calculates net changes in stream channel cross-sections. Calculations are based on dividing the channel cross-section into discrete regions of scour and fill. Internal boundaries of these regions (along the x-axis of the cross-section) are determined by the location of vertical depth measurements taken at two distinct times. The right and left boundaries of the cross-section can be specified so that scour or fill can be calculated for any portion of the profile desired.     

THERMAL AND DISSOLVED OXYGEN CHARACTERISTICS OF A SOUTH CAROLINA COOLING RESERVOIR1

ABSTRACT: Temperature and dissolved oxygen concentrations were measured monthly from January 1971 to December 1982 at 1-m depth intervals at 13 stations in Keowee Reservoir in order to characterize spatial and temporal changes associated with operation of the Oconee Nuclear Station. The reservoir water column was i to 4°C warmer in operational than in non-operational years. The thermo-dine was at depths of 5 to 15 m before the operation of Oconee Nuclear Station, but was always below the upper level of the intake (20 m) after the station was in full operation; this suggests that pumping by the Oconee Nuclear Station had depleted all available cool hypolimnetic water to this depth. As a result summer water temperatures at depths greater than 10 m were usually 10°C higher after plant operation began than before. By fall the reservoir was nearly homothemious to a depth of 27 m, where a thermocine developed. Seasonal temperature profiles varied with distance from the plant; a cool water plume was evident in spring and a warm water plume was present in the summer, fall, and winter. A cold water plume also developed in the northern section of the reservoir due to the operation of Jocassee Pumped Storage Station. Increases in the mean water temperature of the reservoir during operational periods were correlated with the generating output of the power plant. The annual heat load to the reservoir increased by one-third after plant operations began. The alteration of the thermal stratification of the receiving water during the summer also caused the dissolved oxygen to mix to greater depths.     

OVERPUMPED ARTESIAN WELLS AMONG A WELL GROUP1

ABSTRACT. Overpumped artesian wells are defined and explained. Unified well equations which can be used for artesian, gravity and overpumped artesian wells are derived under both the steady and transient conditions of flow. The location of the boundary between artesian and gravity steady and transient flow systems in artesian aquifers, due to multiple wells among which overpumped artesian wells may exist, is explained.     

GROUND WATER REPLENISHMENT FROM RIVERFLOW1

ABSTRACT: Ground water contours and gradients were determined within and adjacent to a gravelly riverbed with about 100 m active channel width and about 800 m flood plain width during low flow, making use of about 60 shallow dug wells and one deeper cased well Ground water slopes away from the riverbed to the south and slopes to the riverbed from the north. Replenishment of ground water to the south is past a “shell,” a zone of low permeability at the margin of and below the riverbed, (1) into a shallow aquifer from “sideway-spill-over” and (2) into an underlying aquifer that is seemingly confined by the “shell.” Seepage from river channels to underflow takes place all along channel lengths except below riffles where the flow direction is reversed. Half the seepage to underflow was lost to ground water replenishment Underflow seems to constitute merely a few percent of total riverflow but makes a disproportionate contribution to ground water replenishment.     

MANAGEMENT OF WATER QUALITY IN RELEASES FROM SOUTHWESTERN IMPOUNDMENTS1

ABSTRACT. The objective of this investigation was to determine the selectivity of withdrawal which is possible in southwestern reservoirs. Two stratified flow solutions were examined to test their applicability under field conditions. Although both appeared capable of accurate prediction of the outflow velocity profile, the Bohan-Grace solution, which required less input data, was utilized to predict the chemical constituents of single and simultaneous releases from several southwestern impoundments. Prediction of outflow water quality was within fifteen percent for southwestern reservoirs as shallow as fifty-five feet. The withdrawal layer thickness for the subject Texas impoundments included the entire hypolimnion or epilmnion depending on outlet location. The sensitivity of the velocity profile to seasonal changes, reservoir discharge rate and withdrawal port dimensions also is illustrated.     

SEMI-INFINITE SOLID MODEL FOR PREDICTION OF TEMPERATURE IN DEEP RESERVOIRS AND LAKES1

ABSTRACT: A lake or reservoir, heated mainly through the air-water interface, is visualized as a solid which is bounded by the interface plane (x = 0) and extends to infinity in the positive x-direction. Analytical expressions give the surface temperature, the temperature profile (with depth) and the position of the thermocline as a function of time (Julian Day). Preliminary model results indicate that the general shape of the temperature profile and the position of the thermocline are influenced primarily by the air-water heat exchange and the time lapse since an isothermal condition existed in the water body.     

COLIFORMS AND WATER QUALITY: USE OF DATA IN PROJECT DESIGN AND OPERATION1

ABSTRACT: Knowledge of coliform transport and disappearance may provide information for project design and operation that minimizes potential water quality problems such as the violation of body contact recreation standards. Storm events were sampled in the Caddo River above DeGray Reservoir, Arkansas, and then tracked through the reservoir using the increased turbidity associated with the storm flows. Fecal coliforms were sampled both in the river and throughout the water column in the reservoir. In general, increased fecal coliform concentrations were closely associated with the increased turbidity resulting from the storm flows. This association existed for all three types of turbidity plume movement - overflow, interflow, and underflow. As the turbidity plume moved down the reservoir, fecal coliform concentrations decreased due to die-off, settling, and dilution. With several assumptions, it is possible to use this information to assist in locating recreational sites in a reservoir or to anticipate possible body contact standard violations at existing recreation sites.     

GIS‐Based Predictive Models of Hillslope Runoff Generation Processes1

Abstract: Successful nonpoint source pollution control using best management practice placement is a complex process that requires in‐depth knowledge of the locations of runoff source areas in a watershed. Currently, very few simulation tools are capable of identifying critical runoff source areas on hillslopes and those available are not directly applicable under all runoff conditions. In this paper, a comparison of two geographic information system (GIS)‐based approaches: a topographic index model and a likelihood indicator model is presented, in predicting likely locations of saturation excess and infiltration excess runoff source areas in a hillslope of the Savoy Experimental Watershed located in northwest Arkansas. Based on intensive data collected from a two‐year field study, the spatial distributions of hydrologic variables were processed using GIS software to develop the models. The likelihood indicator model was used to produce probability surfaces that indicated the likelihood of location of both saturation and infiltration excess runoff mechanisms on the hillslope. Overall accuracies of the likelihood indicator model predictions varied between 81 and 87% for the infiltration excess and saturation excess runoff locations respectively. On the basis of accuracy of prediction, the likelihood indicator models were found to be superior (accuracy 81‐87%) to the predications made by the topographic index model (accuracy 69.5%). By combining statistics with GIS, runoff source areas on a hillslope can be identified by incorporating easily determined hydrologic measurements (such as bulk density, porosity, slope, depth to bed rock, depth to water table) and could serve as a watershed management tool for identifying critical runoff source areas in locations where the topographic index or other similar methods do not provide reliable results.     

DEEP RESERVOIR THERMAL STRATIFICATION MODEL1

ABSTRACT. A method of predicting reservoir temperature profiles is presented. The model requires monthly averaged meteorological and hydrological data along with reservoir physical characteristics. The one-dimensional profiles are predicted by solving the time-dependent energy equation with an assumed epilimnetic depth and a diffusion coefficient applicable below the epilimnion. The model is verified using an existing reservoir.     

Constraints on the in situ density of CO2 within the Utsira formation from time-lapse seafloor gravity measurements

At Sleipner, CO₂ is being separated from natural gas and injected into an underground saline aquifer for environmental purposes. Uncertainty in the aquifer temperature leads to uncertainty in the in situ density of CO₂. In this study, gravity measurements were made over the injection site in 2002 and 2005 on top of 30 concrete benchmarks on the seafloor in order to constrain the in situ CO₂ density. The gravity measurements have a repeatability of 4.3 μGal for 2003 and 3.5 μGal for 2005. The resulting time-lapse uncertainty is 5.3 μGal. Unexpected benchmark motions due to local sediment scouring contribute to the uncertainty. Forward gravity models are calculated based on both 3D seismic data and reservoir simulation models. The time-lapse gravity observations best fit a high temperature forward model based on the time-lapse 3D seismics, suggesting that the average in situ CO₂ density is about to 530 kg/m³. Uncertainty in determining the average density is estimated to be ±65 kg/m³ (95% confidence), however, this does not include uncertainties in the modeling. Additional seismic surveys and future gravity measurements will put better constraints on the CO₂ density and continue to map out the CO₂ flow.     

CRITICAL STATE OF SALT-WATER UPCONING BENEATH ARTESIAN DISCHARGE WELLS1

ABSTRACT: An approach, based on the realization of the vertical components of flow is presented to determine the free surface of gravity wells and the shape of salt-water upconing in artesian aquifers. The transitional stages from pumping fresh water to pumping salt water at the critical condition are discussed. Vertical hydraulic resistances, being determined and included in the flow domain, allow the use of the finite-difference approximation in a simple form. The suggested approach implies the use of either simple modeling or computing techniques.     

BANKFULL HYDRAULIC GEOMETRY RELATIONSHIPS AND RECURRENCE INTERVALS FOR NORTH CAROLINA'S COASTAL PLAIN1

ABSTRACT: Bank full hydraulic geometry relationships relate stream channel geometry to watershed size for specific physiographic regions. This paper presents bank full hydraulic geometry relationships and recurrence intervals for the Southeastern Plain coercion and the flat woods subtype of the Middle Atlantic Coastal Plain ecoregion found within North Carolina's Coastal Plain physiographic province. Cross‐sectional and longitudinal survey data from gated and unpaged streams were used to compute channel dimension and profile information. Power functions were developed, relating drainage area to bank full discharge, cross‐sectional area, width, and mean depth. Recurrence intervals of bank full events were estimated from gagged streams using both a Log‐Pearson Type III distribution of peak annual discharge and a partial‐duration series of average daily discharge. Results from both methods indicate that average bank full recurrence intervals for the study area are below one year. Determinations of recurrence intervals by the Log‐Pearson Type III distribution were for the most part inconclusive (less than one year), while a partial duration series estimated a 0.19 year average, ranging from 0.11 to 0.31 years.     

WATER QUALITY ALLOCATED COST FOR MULTIPURPOSE RESERVOIRS1

ABSTRACT: Low flow augmentation from multipurpose reservoirs may yield significant water quality benefits. Cost allocation assigns a portion of reservoir expense to water quality consumers, waste water dischargers who benefit from increased receiving flow. Whereas such allocation currently is not authorized for Federal projects, the procedure is increasingly appropriate for efficient multiobjective management. Waste water treatment costs, multipurpose reservoir costs, and water quality are modeled for Oregon's Willamette River. Water quality is expressed as a function of treatment and augmentation levels. Treatment cost necessary to achieve a given water quality without augmentation less treatment cost with augmentation is an alternative cost of water quality maintenance. With a cost allocation procedure, this alternative cost is used to determine water quality's share of reservoir cost. Under current conditions, water quality beneficiaries could be charged approximately seven percent of annualized reservoir expense. This charge is one-fourth the expense of additional treatment facilities required were augmentation not provided.     

WAVE INDUCED FLOW AND TRANSPORT IN SEDIMENT BEDS1

ABSTRACT: The fate of contaminants in large water bodies is highly influenced by the transfer of flow and solutes across the water sediment interface. In this paper, an analytical model is presented where flow in both sediment bed and open channel is coupled at the interface through a boundary layer occupying the upper part of the sediment bed. The presence of this layer allows not only the capture of the inertia effects through a drag term in the generalized Darcy's equation, but also the specification of different soil parameters for the two porous zones. The flow is advective and driven by wave action along the water surface. The resulting system is solved for the pressure and flux in each sediment layer. The generated transport velocity fields are linked to a random walk simulation that is used to examine the trajectories of solute particles. Comparison of these trajectories against experimental tracer tests suggests a pattern very similar to the one attributed to the presence of surface mounds. The results clearly show the significance of the boundary layer and the drag term for soil with high permeability and the impact of both the thickness of the boundary layer and the length of the gravity wave relative to the depth of the water channel on the transport and exchange across the interface. The paper also examines the sensitivity of the mass exchange to the permeability of the two porous zones.     

HYDRODYNAMIC CONTROL OF AN EMERGENT AQUATIC PLANT (SCIRPUS ACUTUS) IN OPEN CHANNELS1

ABSTRACT: Control of emergent aquatic plants such as tule (Scirpus acutus Muhl.; Bigel.) is of direct interest to managers of surface waters in Western North America. Where conditions of water velocity and depth occur that permit this and similar species to colonize and grow, their clonal habit may restrict, or even block, open channels within several seasons after their establishment. Fortunately, sufficient flow depth and velocity naturally prevent these plants from growing into and blocking channels. We investigated physical constraints for tule stem growth with the ultimate intent to apply this knowledge in rehabilitating 60 miles of the diverted Owens River in Eastern California, presently choked with emergent growth. Bending stress resulting from hydrodynamic drag on tule stems was found to induce lodging; permanent deformation and consequent loss of function. The depth-velocity envelope describing this process (at 95 percent confidence) is uD/d= 12.8 where u = average velocity acting upon the stem (m/s), D = local depth of flow (m), and d = tule stem diameter at the point of attachment (m). Maintaining a discharge or reconfiguring a channel so this critical depth-velocity-stem diameter envelope is exceeded (predictable using flow models) through the summer growing period should prevent encroachment into an active channel.