Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Hydraulic Engineering”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5Linked to original sources

Interrelations among the soil-water retention, hydraulic conductivity, and suction-stress characteristic curves

The three fundamental constitutive relations that describe fluid flow, strength, and deformation behavior of variably saturated soils are the soil-water retention curve (SWRC), hydraulic conductivity function (HCF), and suction-stress characteristic curve (SSCC). Until recently, the interrelations among the SWRC, HCF, and SSCC have not been well established. This work sought experimental confirmation of interrelations among these three constitutive functions. Results taken from the literature for six soils and those obtained for 11 different soils were used. Using newly established analytical relations among the SWRC, HCF, and SSCC and these test results, the authors showed that these three constitutive relations can be defined by a common set of hydromechanical parameters. The coefficient of determination for air-entry pressures determined independently using hydraulic and mechanical methods is >0.99, >0.98 for the pore size parameter, and 0.94 for the residual degree of saturation. One practical implication is that one of any of the four experiments (axis-translation, hydraulic, shear-strength, or deformation) is sufficient to quantify all three constitutive relations.

Journal of Geotechnical and Geoenvironmental Engin

Long-term behavior of water content and density in an earthen liner

An extensively instrumented compacted earthen liner was constructed at the Illinois State Geological Survey facility in Champaign, III. in 1987. A pond of water 0.31 m deep was maintained on top of the 7.3 m ?? 14.6 m ?? 0.9 m thick liner for 14 years. One of the goals of the project was to evaluate the long-term performance of a compacted earthen liner by monitoring the long-term changes in water content and density. The water content of the earthen liner showed no trend with depth or time. The liner density remained essentially constant from construction through excavation in 2002. The liner did not become fully saturated. Upon excavation of the liner, the degree of saturation was 80.0??6.3% after 14 years of ponding under a hydraulic head of 0.31 m. The results imply that properly designed and constructed earthen liners may reduce the possibility of pollutants leaching from municipal solid waste containment facilities by remaining partially saturated for years and maintaining the placement density. Journal of Geotechnical and Geoenvironmental Engineering ?? ASCE.

Journal of Geotechnical and Geoenvironmental Engin

Effect of hydraulic hysteresis on the stability of infinite slopes under steady infiltration

Hydraulic hysteresis, including capillary soil water retention (SWR), air entrapment SWR, and hydraulic conductivity, is a common phenomenon in unsaturated soils. However, the influence of hydraulic hysteresis on suction stress, and subsequently slope stability, is generally ignored. This paper examines the influence of each of these three types of hysteresis on slope stability using an infinite slope stability analysis under steady infiltration conditions. First, hypothetical slopes for representative silty and sandy soils are examined. Then a monitored hillslope in the San Francisco Bay Area, California is assessed, using observed rainfall conditions and measured hydraulic and geotechnical properties of the colluvial soil. Results show that profiles of suction stress and the corresponding factor of safety are generally strongly affected by hydraulic hysteresis. Results suggest that each of the three types of hydraulic hysteresis may play a major role in the occurrence of slope failure, indicating that ignoring hydraulic hysteresis will likely lead to underestimates of failure potential and hence to inaccurate slope stability analysis.

Journal of Geotechnical and Geoenvironmental Engin

Performance of bedload sediment transport formulas applied to the Lower Minnesota River

Despite limitations in reproducing complex bedload sediment transport processes in rivers, formulas have been preferred over collection and analysis of field data due to the high cost and time-consuming nature of bedload discharge measurements. However, the performance of such formulas depends on the hydraulic and sedimentological conditions they attempt to describe. The availability of field measurements provides a unique opportunity to test bedload transport formulas to better guide formula selection. Hydraulic parameters and bedload discharge data from the Lower Minnesota River and two of its tributaries were used to evaluate nine bedload transport formulas using three different indices. The bedload data for the different sites were collected by the United States Geological Survey (USGS) from 2011 through 2014, with bed material varying from very coarse to medium sand. The formulas calculated higher bedload rates than were measured due to a combination of site-specific physical characteristics, including the presence of bed forms (dunes), and sampling uncertainties. Because of the lack of reproducibility of the tested formulas, five power functions, based on the relation between the specific unit power (independent hydraulic variable) and the USGS measured data (dependent variable), were derived as provisional equations to estimate the bedload discharge on the Lower Minnesota River and tributaries.

Minnesota

Coupled semivariogram uncertainty of hydrogeological and geophysical data on capture zone uncertainty analysis

This study investigates capture zone uncertainty that relates to the coupled semivariogram uncertainty of hydrogeological and geophysical data. Semivariogram uncertainty is represented by the uncertainty in structural parameters (range, sill, and nugget). We used the beta distribution function to derive the prior distributions of structural parameters. The probability distributions of structural parameters were further updated through the Bayesian approach with the Gaussian likelihood functions. Cokriging of noncollocated pumping test data and electrical resistivity data was conducted to better estimate hydraulic conductivity through autosemivariograms and pseudo-cross-semivariogram. Sensitivities of capture zone variability with respect to the spatial variability of hydraulic conductivity, porosity and aquifer thickness were analyzed using ANOVA. The proposed methodology was applied to the analysis of capture zone uncertainty at the Chicot aquifer in Southwestern Louisiana, where a regional groundwater flow model was developed. MODFLOW-MODPATH was adopted to delineate the capture zone. The ANOVA results showed that both capture zone area and compactness were sensitive to hydraulic conductivity variation. We concluded that the capture zone uncertainty due to the semivariogram uncertainty is much higher than that due to the kriging uncertainty for given semivariograms. In other words, the sole use of conditional variances of kriging may greatly underestimate the flow response uncertainty. Semivariogram uncertainty should also be taken into account in the uncertainty analysis. ?? 2008 ASCE.

Journal of Hydrologic Engineering

Conservation-form equations of unsteady open-channel flow

The unsteady open-channel flow equations are typically expressed in a variety of forms due to the imposition of differing assumptions, use of varied dependent variables, and inclusion of different source/sink terms. Questions often arise as to whether a particular equation set is expressed in a form consistent with the conservation-law definition. The concept of conservation form is developed to clarify the meaning mathematically. Six sets of unsteady-flow equations typically used in engineering practice are presented and their conservation properties are identified and discussed. Results of the theoretical development and analysis of the equations are substantiated in a set of numerical experiments conducted using alternate equation forms. Findings of these analytical and numerical efforts demonstrate that the choice of dependent variable is the fundamental factor determining the nature of the conservation properties of any particular equation form.

Journal of Hydraulic Research

Time-averaged near-bed suspended sediment concentrations under waves and currents: Comparison of measured and model estimates

Profiles of suspended sediment concentration and velocity were measured over a 15-day period at a near-shore site off Santa Cruz, CA in Monterey Bay. The concentration and velocity data were collected with an Acoustic Backscattering System (ABS) and Acoustic Current Profiler (ACP) that were mounted on a bottom tripod. High-resolution bottom scanning sonar was also attached to the tripod to provide images of bed features during the experiment. Hourly time-averaged near-bed concentrations of suspended sediment were calculated from three models and compared with the measurements. Surface waves and currents that were generated by a storm of moderate intensity caused bed stresses that exceeded threshold stress for D 50 =0.02 cm, the median size of the moderately well-sorted bottom sediment, over a period of about 7 days. Estimates of the concentration at 1 cm above the bottom, C a1 , were obtained using the ABS measurements. These observations have been compared with predictions for the concentration at 1 cm above the bottom, C 1 . Nielsen's models for reference concentration C o [Nielsen, P., 1986. Suspended sediment concentrations under waves. Coastal Engineering 10, 32-31; Nielsen, P., 1992. Coastal Bottom Boundary Layers and Sediment Transport, Advanced Series on Ocean Engineering. World Scientific, Hackensack, NJ.] are purely wave-based and do not include effects of bottom currents on bed stress and bedform scales. C 1 calculated from this model compared well with measured C a1 when currents were weak and small oscillatory ripples were observed in the sonar images. However, during the 3-day period of highest bottom stresses modeled C 1 did not compare well to C a1 . The other two models for C 1 , Glenn and Grant [Glenn, S.M., Grant, W.D., 1987. A suspended sediment stratification correction for combined wave and current flows. Journal of Geophysical Research 92(C8), 8244-8264.] and van Rijn and Walstra [Van Rijn, L.C., Walstra, D.J.R., 2004. Description of TRANSPOR2004 and implementation in Delft3D-ONLINE. Interim Report prepared for DG Rijkswaterstaat, Rijksinstituut voor Kust en Zee. Delft Hydraulics Institute, The Netherlands.], accounted for combined wave-current stresses and included different formulations for predicting bedform scales. C 1 calculated from these models compared well with measurements throughout the early and middle period of the storm, but underpredicted measured values during the latter part of the experiment. An alternative method based on dimensional analysis provided a straightforward way to relate C 1 to bed stress and sediment parameters. This new relationship based on dimensional analysis explained about 77% of the variability in C a1 derived from the ABS data.

California

Sequestration and transformation in chemically enhanced treatment wetlands: DOC, DBPPs and Nutrients

We examined the effectiveness of chemically enhanced treatment wetlands (CETWs), wetlands that received water treated with coagulants, to remove dissolved organic carbon (DOC), disinfection byproduct precursors (DBPPs), nutrients and metals from agricultural drain water. Wetlands consisted of controls with no coagulant addition, ferric sulfate dosed and polyaluminum chloride dosed treatments. CETWs were more effective in removal of DOC, DBPPs, phosphate, dissolved organic nitrogen and metals than control wetlands. Coagulation treated wetlands removed 245 – 349 g/m2yr DOC, whereas control wetlands produced 51 g/m2yr. Wetland passage released DOC in the controls and treatments; this release was directly correlated to temperature and suggested thermally dependent mechanism(s) were partly responsible. A first-order plug flow reactor kinetic model that considered hydraulic retention time (HRT), temperature and concentration was tested for wetland DOC processing. Models indicate that operating CETWs at high coagulant dosing and low temperature can result in lowest DOC release with additional release suppression. Operating at the lowest HRT to meet discharge targets help overcome wetland processes that increase DOC release and provide the smallest footprint needed for treatment.

Journal of Environmental Engineering

Slope instability caused by small variations in hydraulic conductivity

Variations in hydraulic conductivity can greatly modify hillslope ground-water flow fields, effective-stress fields, and slope stability. In materials with uniform texture, hydraulic conductivities can vary over one to two orders of magnitude, yet small variations can be difficult to determine. The destabilizing effects caused by small (one order of magnitude or less) hydraulic conductivity variations using ground-water flow modeling, finite-element deformation analysis, and limit-equilibrium analysis are examined here. Low hydraulic conductivity materials that impede downslope ground-water flow can create unstable areas with locally elevated pore-water pressures. The destabilizing effects of small hydraulic heterogeneities can be as great as those induced by typical variations in the frictional strength (approximately 4°–8°) of texturally similar materials. Common “worst-case” assumptions about ground-water flow, such as a completely saturated “hydrostatic” pore-pressure distribution, do not account for locally elevated pore-water pressures and may not provide a conservative slope stability analysis. In site characterization, special attention should be paid to any materials that might impede downslope ground-water flow and create unstable regions.

Journal of Geotechnical and Geoenvironmental Engin

Evaluation of electrical and electromagnetic geophysical techniques to inspect earthen dam and levee structures in Arkansas

Within the state of Arkansas there is an increasing number of aging dams and levees that have little to no documentation concerning their construction or composition. Surface geophysical surveys offer a non-intrusive method for investigating these structures: To describe their lithologic makeup, to evaluate the materials that they were constructed upon, and to identify potential flow paths through them. Techniques such as electrical resistivity tomography, seismic refraction, and electromagnetic induction have all been used to image dams and levees and require additional information from geologic outcrops, geotechnical borings, or drill cores in order to make informed geologic interpretations of the geophysical models. These geologic models then allow the owners of these structures to make more informed decisions about their operation and maintenance. Between 2011 and 2018, the U.S. Geological Survey conducted geophysical and geotechnical investigations of three earthen structures within the state of Arkansas. Electrical and electromagnetic geophysical data were used to develop lithologic models of these structures and the underlying geology. Self-potential surveys were utilized to detect the movement of water through these structures indicating possible seepage pathways. Geotechnical methods such as electric and hydraulic direct-push well logs and cores acted as both a control on the geophysical interpretations and a confirmation of anomalies. This integrated approach detected the lack of an impermeable core within a levee, imaged a change in lithology of the bedrock forming the seal beneath a gravity dam, and identified a potential seepage feature within the core of an earthen dam. These results further support that this method of extending known lithologic features via surface and borehole geophysics is a useful approach for characterizing earthen water control structures.

Arkansas

New insights on scale-dependent surface-groundwater exchange from a floating self-potential Dipole

In south-central Texas the lower Guadalupe River has incised into the outcrop of the Carrizo-Wilcox aquifer. The river and the aquifer are hydraulically connected across the outcrop, although the connectivity is obscured at the surface by alluvium and surface-water and groundwater exchange dynamics are currently poorly understood. To investigate surface-water and groundwater exchange dynamics between the lower Guadalupe River and the Carrizo-Wilcox aquifer, a geophysical study was completed along a 14.86 km reach of the river by using water-borne gradient self-potential (SP) profiling and two-dimensional direct-current electric resistivity tomography. This paper explores the applicability of these water-borne geoelectric methods in delineating gaining and losing channel reaches, and demonstrates that geoelectric signals in the form of total electric field strength can be logged with an electric dipole and decomposed into component SP signals depicting regional and local groundwater flow patterns attributable to regional and localized hydraulic gradients. Localized SP anomalies of several tens of millivolts, indicative of hyporheic exchange flows, are observed and superimposed upon a 124 mV regional SP anomaly indicative of ambient groundwater exchange flows between the river and the aquifer. The observed SP signals are interpreted through two-dimensional finite-element modeling of streaming potentials attributable to ambient groundwater exchange and hyporheic exchange flow patterns. Variables of the channel environment such as temperature and concentration gradients, depth, and velocity are considered and subsequently eliminated as alternative sources of the SP signals that are presented.

Texas

Investigation of scale-dependent groundwater/surface-water exchange in rivers by gradient self-potential logging: Numerical modeling and field experiments

Exchanges of groundwater and surface-water are fundamental to a wide range of water-supply and water-quality management issues but challenging to map beyond the reach scale. Waterborne gradient self-potential (SP) measurements are directly sensitive to water flow through riverbed sediments and can be used to infer exchange locations, direction (gain versus loss), scale, and relative changes, but to date applications to river corridor hydrology are limited. Numerical modeling and field experiments were therefore performed herein, each emphasizing waterborne gradient SP logging for identifying and locating focused vertical groundwater discharge (surface-water gain) and recharge (surface-water loss) in a river. Two and three-dimensional numerical models were constructed to simulate the polarities, appearances, and peak amplitudes of streaming-potential and electric-field anomalies on a riverbed and in the surface-water that were attributable to steady-state vertical fluxes of groundwater through high-permeability conduits in the riverbed. Effects of varied hydraulic length-scale of exchange and surface-water depth were tested through numerical modeling. Modeling results aided in data acquisition and interpretation for three repeated field experiments performed along a 1.5–2.0 km reach of the Quashnet River in Cape Cod, Massachusetts, where focused, meter-scale groundwater discharges occur at discrete locations within otherwise ubiquitous and more diffuse groundwater upwelling conditions. Strong gradient SP anomalies were repeatedly measured in the Quashnet River at previously confirmed locations of focused groundwater discharge, showing the efficacy of waterborne gradient SP logging in identifying and characterizing groundwater/surface water exchange dynamics at multiple river network scales.

Massachusetts

Well predictive performance of play-wide and Subarea Random Forest models for Bakken productivity

In recent years, geologists and petroleum engineers have struggled to clearly identify the mechanisms that drive productivity in horizontal, hydraulically-fractured oil wells producing from the middle member of the Bakken formation. This paper fills a gap in the literature by showing how this play’s heterogeneity affects factors that drive well productivity. It is important because understanding the relative strength of productivity drivers and how predictors vary spatially facilitates best-practices for well site selection and well completion design. The paper describes an application of the Random Forest (RF) machine learning technique to identify these mechanisms and to evaluate their importance across 9 subareas of the North Dakota portion of the Bakken play. The study examined productivity of 7311 wells initiating production from 2010 through 2017. Well productivity varied considerably across the 9 subareas within the play, so it was not surprising that the dominant predictors, the initial 180-day water cut and the 30-day initial gas production, vary spatially to mirror local conditions that strongly affect well productivity. The relative importance of well completion predictor variables, that is, the numbers of fractures stages per well, volume of injected proppant per stage, volume of injected fluids per stage, and lateral length, varied considerably across the subareas. Statistical permutation tests are presented that generally confirm the importance rankings. Subarea Random Forest models explained from 50 percent to 82 percent of the variation in productivity test samples while the play-wide model explained 73 percent of the test sample well productivity. Weakness in the predictive ability of the Random Forest models are traced to the limited variability in the training data. Implications of the empirical findings regarding the Bakken play for operators and for research and government institutions are discussed in the concluding section.

Montana, North Dakota, South Dakota

Integrating surface and borehole geophysics in ground water studies: An example using electromagnetic soundings in south Florida

Time domain surface electromagnetic soundings, borehole induction logs, and other borehole logging techniques are used to construct a realistic model for the shallow subsurface hydraulic properties of unconsolidated sediments in south Florida. Induction logs are used to calibrate surface induction soundings in units of pore water salinity by correlating water sample specific electrical conductivity with the electrical conductivity of the formation over the sampled interval for a two‐layered aquifer model. Geophysical logs are also used to show that a constant conductivity layer model is appropriate for the south Florida study. Several physically independent log measurements are used to quantify the dependence of formation electrical conductivity on such parameters as salinity, permeability, and clay mineral fraction. The combined interpretation of electromagnetic soundings and induction logs was verified by logging three validation boreholes, confirming quantitative estimates of formation conductivity and thickness in the upper model layer, and qualitative estimates of conductivity in the lower model layer.

Journal of Environmental & Engineering Geophysics

Using the KINEROS2 modeling framework to evaluate the increase in storm runoff from residential development in a semi-arid environment

The increase in runoff from urbanization is well known; one extreme example comes from a 13 hectare residential neighborhood in southeast Arizona where runoff was 27 times greater than an adjacent grassland watershed over a forty‐month period from 2005 to 2008. Rainfall‐runoff modeling using the newly‐described KINEROS2 urban element and tension infiltrometer measurements indicate that 17±14 percent of this increase in runoff is due to a 53 percent decrease in the saturated hydraulic conductivity of constructed pervious areas, as compared to the undeveloped grassland. Directly connected impervious areas, primarily streets and driveways, cause 56 percent of the increase in runoff, and indirectly connected impervious areas, primarily rooftops and sidewalks, and a decrease in canopy interception account for the remaining 27 percent increase. Tension infiltrometer measurements show that saturated hydraulic conductivity (K s ) is about double in the grassland watershed than in the urban watershed, 6.2 ± 3.5mm/hr and 2.9 ± 1.6mm/hr, respectively. K s in the urban watershed identified from calibrating the rainfall‐runoff model to measured runoff is 9.5 ± 2.8 mm/hr—higher than what was measured but much lower than the 26 mm/hr value indicated by a soil‐texture based KINEROS2 parameter look‐up table. A new component of the KINEROS2 modeling framework, the urban element, forms the basis for the model by simulating a contiguous row of houses and the adjoining street as a series of pervious and impervious overland flow planes. Tests using different levels of discretization found that watershed geometry can be represented in a simplified manner, although more detailed discretization led to better model performance.

Arizona

Use of Doppler velocity radars to monitor and predict debris and flood wave velocities and travel times in post-wildfire basins

The magnitude and timing of extreme events such as debris and floodflows (collectively referred to as floodflows) in post-wildfire basins are difficult to measure and are even more difficult to predict. To address this challenge, a sensor ensemble consisting of noncontact, ground-based (near-field), Doppler velocity (velocity) and pulsed (stage or gage height) radars, rain gages, and a redundant radio communication network was leveraged to monitor flood wave velocities, to validate travel times, and to compliment observations from NEXRAD weather radar. The sensor ensemble (DEbris and Floodflow Early warNing System, DEFENS) was deployed in Waldo Canyon, Pike National Forest, Colorado, USA, which was burned entirely (100 percent burned) by the Waldo Canyon fire during the summer of 2012 ( MTBS, 2020 ). Surface velocity, stage, and precipitation time series collected during the DEFENS deployment on 10 August 2015 were used to monitor and predict flood wave velocities and travel times as a function of stream discharge (discharge; streamflow). The 10 August 2015 event exhibited spatial and temporal variations in rainfall intensity and duration that resulted in a discharge equal to 5.01 cubic meters per second (m 3 /s). Discharge was estimated post-event using a slope-conveyance indirect discharge method and was verified using velocity radars and the probability concept algorithm. Mean flood wave velocities – represented by the kinematic celerity ( c k = 2.619 m e t e r s p e r s e c o n d , m / s ± 0.556 p e r c e n t "> 𝑐 𝑘 = 2.619 𝑚𝑒𝑡𝑒𝑟𝑠 𝑝𝑒𝑟 𝑠𝑒𝑐𝑜𝑛𝑑, m/s ± 0.556 𝑝𝑒𝑟𝑐𝑒𝑛𝑡) and dynamic celerity ( c d = 3.533 m / s ± 0.181 p e r c e n t a n d t h e i r u n c e r t a i n t i e s "> 𝑐 𝑑 = 3.533 m/s ± 0.181 𝑝𝑒𝑟𝑐𝑒𝑛𝑡) and their uncertainties were computed. L-moments were computed to establish probability density functions (PDFs) and associated statistics for each of the at-a-section hydraulic parameters to serve as a workflow for implementing alert networks in hydrologically similar basins that lack data. Measured flood wave velocities and travel times agreed well with predicted values. Absolute percent differences between predicted and measured flood wave velocities ranged from 1.6 percent to 49 percent and varied with water slope, hydraulic radius, and depth. The kinematic celerity was a better predictor for steep slopes and wide flood plains associated with the Upper Waldo and Middle Waldo radar streamgages; whereas, the dynamic celerity was a better surrogate for shallow slopes and incised channels such as the Lower Waldo radar streamgage. The method demonstrates the potential extensibility of a post-wildfire warning system by (1) leveraging multiple systems (i.e., weather radar, near-field velocity and stage radars, and rain gages) for accurate and timely warnings of debris and floodflows, (2) establishing an order of operations to site, install, and operate near-field radars and conventional rain gages to record floodflows, forecast travel times, and document geomorphic change in this basin and hydrologically similar basins that lack data, and (3) communicating data operationally with the Colorado Department of Transportation engineering staff, National Weather Service forecasters, and emergency managers.

Colorado

Volatilization of organic compounds from streams

Mass-transfer coefficients for the volatilization of ethylene and propane were correlated with the hydraulic and geometric properties of seven streams, and predictive equations were developed. The equations were evaluated using a normalized root-mean-square error as the criterion of comparison. The two best equations were a two-variable equation containing the energy dissipated per unit mass per unit time and the average depth of flow and a three-variable equation containing the average velocity, the average depth of flow, and the slope of the stream. Procedures for adjusting the ethylene and propane coefficients for other organic compounds were evaluated. These procedures are based on molecular diffusivity, molecular diameter, or molecular weight. Because of limited data, none of these procedures have been extensively verified. Therefore, until additional data become available, it is suggested that the mass-transfer coefficient be assumed to be inversely proportional to the square root of the molecular weight.

Journal of Environmental Engineering

An engineering economic analysis of a program for artificial groundwater recharge

This study describes and demonstrates two alternate methods for evaluating the relative costs and benefits of artificial groundwater recharge using percolation ponds. The first analysis considers the benefits to be the reduction of pumping lifts and land subsidence; the second considers benefits as the alternative costs of a comparable surface delivery system. Example computations are carried out for an existing artificial recharge program in Santa Clara Valley in California. A computer groundwater model is used to estimate both the average long term and the drought period effects of artificial recharge in the study area. For the example problem, the benefits of reduced average annual pumping lifts and reduced incremental subsidence are greater than the total costs of continuing the existing artificial recharge program. Benefits for reduced subsidence are strongly dependent on initial aquifer conditions. The second analysis compares the costs of continuing the artificial recharge program with the costs of a surface system which would achieve the same hydraulic effects. Results indicate that the costs of artificial recharge are considerably smaller than the alternative costs of an equivalent surface system. In evaluating a particular program, consideration should also be given to uncertainties in future supplies and demands for water as well as to the probability of extreme events such as droughts

California