Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Hydraulic Research”

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 19 recordsLinked to original sources

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↗

A unique approach to estimating lateral anisotropy in complex geohydrologic environments

Aquifers in fractured rock or karstic settings are likely to have anisotropic transmissivity distributions. Aquifer tests that are performed in these settings also are frequently affected by leakage from adjacent confining units. Finite-difference models such as MODFLOW are convenient tools for estimating the hydraulic characteristics of the stressed aquifer and adjacent confining units but are poor tools for the estimation of lateral anisotropy. This limitation of finite-difference methods can be overcome by application of the spin method, a technique whereby the positions of the observation wells are rotated about the production well to estimate anisotropy and orientation. Formal parameter estimation is necessary to analyze aquifer tests because of the number of parameters that are estimated. As a test, transmissivity, anisotropy, and orientation were successfully estimated for a simple hypothetical problem with known properties. The technique also was applied to estimate hydraulic properties of the Santee Limestone/Black Mingo (SL/BM) aquifer and a leaky confining unit beneath Charleston, South Carolina. A 9-day aquifer test with an average discharge of 6441/min was analyzed numerically. Drawdowns in the SL/BM aquifer and confining unit were simulated with a 12-layer MODFLOW model that was discretized into 81 rows of 81 columns. Simulated drawdowns at seven observation wells that ranged from 23 to 2700 m from the production well were matched to measured drawdowns. Transmissivity estimated along the minor axis ranged from 10 to 15 m 2 /day and along the major axis ranged from 80 to 10Om 2 /day. The major axis of transmissivity was oriented along compass heading 116° (degrees clockwise from north), which agrees with geologic interpretations. Vertical hydraulic conductivity and specific storage estimates for the overlying confining unit were 4 x 10 -5 m/day and 2 x 10 -4 1/m, respectively.

South Carolina↗

Modeling downstream fining in sand-bed rivers. I: Formulation

In this paper a numerical modeling formulation is presented for simulation of the development of the longitudinal profile and bed sediment distribution in sand-bed rivers. The objective of the model application, which is presented in the companion paper (Wright and Parker, 2005), is to study the development of two characteristics of large, low-slope, sand-bed rivers: (1) a downstream decrease in bed slope (i.e. concave upward longitudinal profile) and (2) a downstream decrease in characteristic bed sediment diameter (e.g. the median bed surface size D 50 ). Three mechanisms that lead to an upward concave profile and downstream fining are included in the modeling formulation: (1) a delta prograding into standing water at the downstream boundary, (2) sea-level rise, and (3) tectonic subsidence. In the companion paper (Wright and Parker, 2005) the model is applied to simulate the development of the longitudinal profile and downstream fining in sand-bed rivers flowing into the ocean during the past 5000 years of relatively slow sea-level rise.

Journal of Hydraulic Research↗

Modeling downstream fining in sand-bed rivers. II: Application

In this paper the model presented in the companion paper, Wright and Parker (2005) is applied to a generic river reach typical of a large, sand-bed river flowing into the ocean in order to investigate the mechanisms controlling longitudinal profile development and downstream fining. Three mechanisms which drive downstream fining are studied: a delta prograding into standing water, sea-level rise, and tectonic subsidence. Various rates of sea-level rise (typical of the late Holocene) and tectonic subsidence are modeled in order to quantify their effects on the degree of profile concavity and downstream fining. Also, several other physical mechanisms which may affect fining are studied, including the relative importance of the suspended versus bed load, the effect of the loss of sediment overbank, and the influence of the delta bottom slope. Finally, sensitivity analysis is used to show that the grain-size distribution at the interface between the active layer and substrate has a significant effect on downstream fining.

Journal of Hydraulic Research↗

Flow over bedforms in a large sand-bed river: A field investigation

An experimental field study of flows over bedforms was conducted on the Missouri River near St. Charles, Missouri. Detailed velocity data were collected under two different flow conditions along bedforms in this sand-bed river. The large river-scale data reflect flow characteristics similar to those of laboratory-scale flows, with flow separation occurring downstream of the bedform crest and flow reattachment on the stoss side of the next downstream bedform. Wave-like responses of the flow to the bedforms were detected, with the velocity decreasing throughout the flow depth over bedform troughs, and the velocity increasing over bedform crests. Local and spatially averaged velocity distributions were logarithmic for both datasets. The reach-wise spatially averaged vertical-velocity profile from the standard velocity-defect model was evaluated. The vertically averaged mean flow velocities for the velocity-defect model were within 5% of the measured values and estimated spatially averaged point velocities were within 10% for the upper 90% of the flow depth. The velocity-defect model, neglecting the wake function, was evaluated and found to estimate thevertically averaged mean velocity within 1% of the measured values.

Missouri↗

Field methods, quality-assurance, and data management plan for water-quality activities and water-level measurements, Idaho National Laboratory, Idaho

Introduction Water-quality activities and water-level measurements conducted by the U.S. Geological Survey (USGS) Idaho National Laboratory (INL) Project Office coincide with the USGS mission of evaluating the quantity and quality of the Nation’s water resources. The activities are conducted in cooperation with the U.S. Department of Energy’s (DOE) Idaho Operations Office. Results of water-quality and hydraulic head research efforts are presented in various USGS and scientific journal publications (refer to Fisher, 2022). These data are stored internally in the Aquarius Time Series and Aquarius Samples databases and are publicly accessible through National Water Quality Monitoring Council (2025) and U.S. Geological Survey (2025). Data collected from our studies are used by researchers, Federal and State agencies, water management and regulatory organizations, as well as the public. This quality assurance plan (QAP) describes the methods and processes for field methods, data collection, data management, data auditing, and equipment management for both the water-quality and water-level programs at the USGS INL Project Office (hereto referred to as INL Project Office). A comprehensive quality assurance (QA) plan ensures that the processes defined in this document will guide the program staff to collect and publish reliable, useful, and defensible data products for stakeholders. This QAP supersedes previous versions of this document and is intended to complement the Quality Assurance and Data Management (QADM) Plan for the Idaho Water Science Center (IDWSC; Christopher Mebane and Lauren Zinsser, written commun., 2024).

Idaho↗

Effects of dredged channels on trace-metal migration in an estuary

Determination of trace-metal levels in the sediments of the Matagorda Bay system revealed anomalously high mercury values. The distribution of the mercury-rich sediment deposits is the result of the sedimentological regime of the bay system produced by the tidal currents in the dredged channel. According to this model, the oxygenated open gulf water pushed into the bay by tidal currents activates the mercury and reintroduces it into the sediment regime of the bay in an area where the turbidity maximum is most prevalent. Within this region, the absorption sites are at a maximum, tying up the mercury. The mercury-enriched suspended material is then transported and deposited according to the hydraulic regime within the bay system.

Texas↗

Non-double-couple mechanisms of microearthquakes induced by hydraulic fracturing

We have inverted polarity and amplitude information of representative microearthquakes to investigate source mechanisms of seismicity induced by hydraulic fracturing in the Carthage Cotton Valley, east Texas, gas field. With vertical arrays of four and eight three-component geophones in two monitoring wells, respectively, we were able to reliably determine source mechanisms of the strongest events with the best signal-to-noise ratio. Our analysis indicates predominantly non-double-couple source mechanisms with positive volumetric component consistent with opening cracks oriented close to expected hydraulic fracture orientation. Our observations suggest the induced events are directly the result of opening cracks by fluid injection, in contrast to many previous studies where the seismicity is interpreted to be primarily shearing caused by pore pressure diffusion into the surrounding rock or associated with shear stresses created at the hydraulic fracture tip.

Journal of Geophysical Research B: Solid Earth↗

Field methods and quality-assurance plan for water-quality activities and water-level measurements, U.S. Geological Survey, Idaho National Laboratory, Idaho

Water-quality activities and water-level measurements by the personnel of the U.S. Geological Survey (USGS) Idaho National Laboratory (INL) Project Office coincide with the USGS mission of appraising the quantity and quality of the Nation’s water resources. The activities are carried out in cooperation with the U.S. Department of Energy (DOE) Idaho Operations Office. Results of the water-quality and hydraulic head investigations are presented in various USGS publications or in refereed scientific journals and the data are stored in the National Water Information System (NWIS) database. The results of the studies are used by researchers, regulatory and managerial agencies, and interested civic groups. In the broadest sense, quality assurance refers to doing the job right the first time. It includes the functions of planning for products, review and acceptance of the products, and an audit designed to evaluate the system that produces the products. Quality control and quality assurance differ in that quality control ensures that things are done correctly given the “state-of-the-art” technology, and quality assurance ensures that quality control is maintained within specified limits.

Open-File Report↗

Field methods, quality-assurance, and data management plan for water-quality activities and water-level measurements, Idaho National Laboratory, Idaho

Introduction Water-quality activities and water-level measurements conducted by the U.S. Geological Survey (USGS) Idaho National Laboratory (INL) Project Office coincide with the USGS mission of appraising the quantity and quality of the Nation’s water resources. The activities are conducted in cooperation with the U.S. Department of Energy’s (DOE) Idaho Operations Office. Results of water-quality and hydraulic head investigations are presented in various USGS publications or in refereed scientific journals, and the data are stored in the National Water Information System (NWIS) database. The results of the studies are used by researchers, regulatory and managerial agencies, and civic groups. In its broadest sense, “quality assurance” refers to doing the job right the first time. It includes the functions of planning for products, review and acceptance of the products, and an audit designed to evaluate the system that produces the products. Quality control and quality assurance differ in that quality control ensures that things are done correctly given the “state-of-the-art” technology, and quality assurance ensures that quality control is maintained within specified limits.

Idaho↗

Reconstructing depositional processes and history from reservoir stratigraphy: Englebright Lake, Yuba River, northern California

Reservoirs provide the opportunity to link watershed history with its stratigraphic record. We analyze sediment cores from a northern California reservoir in the context of hydrologic history, watershed management, and depositional processes. Observations of recent depositional patterns, sediment‐transport calculations, and 137 Cs geochronology support a conceptual model in which the reservoir delta progrades during floods of short duration (days) and is modified during prolonged (weeks to months) drawdowns that rework topset beds and transport sand from topsets to foresets. Sediment coarser than 0.25–0.5 mm deposits in foresets and topsets, and finer material falls out of suspension as bottomset beds. Simple hydraulic calculations indicate that fine sand (0.063–0.5 mm) is transported into the distal bottomset area only during floods. The overall stratigraphy suggests that two phases of delta building occurred in the reservoir. The first, from dam construction in 1940 to 1970, was heavily influenced by annual, prolonged >20 m drawdowns of the water level. The second, built on top of the first, reflects sedimentation from 1970 to 2002 when the influence of drawdowns was less. Sedimentation rates in the central part of the reservoir have declined ∼25% since 1970, likely reflecting a combination of fewer large floods, changes in watershed management, and winnowing of stored hydraulic mining sediment.

Journal of Geophysical Research F: Earth Surface↗

Hydrogeology of the Hawaii Scientific Drilling Project borehole KP-1 1. Hydraulic conditions adjacent to the well bore

Temperature and formation resistivity logs obtained in borehole KP-1 of the Hawaii Scientific Drilling Project indicate that the adjacent formation is characterized by several zones of distinctly different average temperature and water salinity. A series of hydraulic analyses and water sampling programs were conducted to rule out the possibility of local hydraulic effects associated with the presence of the borehole in the generation of these apparent groundwater zones. Hydraulic tests and sampling with the borehole cased to a depth of 710 m and open below that depth indicate that the deep aquifer contains seawater at a temperature nearly identical to that of the open ocean at the same depth. Various analyses give estimates of aquifer transmissivity of about 10-3 m2/s in the vicinity of the borehole. Isolation of this deeper aquifer from the overlying groundwater zones was investigated by perforating the casing at six locations and then measuring the changes in water level in the borehole, in the salinity of the fluid column, in the temperature profile of the fluid column, and in the rate of flow in the fluid column induced by the perforations. These results positively confirm that the zones of distinctly different formation properties indicated on the temperature and resistivity logs are not caused by flow in or around casing. Flow and fluid column salinity induced by the perforations also confirm significant differences between the hydraulic heads and geochemistry of the different groundwater zones inferred from the well logs.

Journal of Geophysical Research B: Solid Earth↗

Constraining the permeability and outer-rise hydration at the Central America margin

Subduction zones are important drivers of Earth's geochemical evolution; the presence of water within the subduction system has substantial implications for geohazards, including the style of seismicity and the volume of magma generated at arc volcanoes. We utilize estimates of pore fluid distributions in the outer-rise region of the Central American margin imaged using controlled-source electro-magnetic surveys and novel two-dimensional forward models simulating poro-elastic deformation of slab bending to constrain the background permeability structure of the Cocos Plate and within outer-rise fault damage zones. The model results suggest that a relatively low background permeability structure is required to match the observed porosity at depths greater than 3 km (less than ~10 −21 m 2 . Similarly, increasing the permeability within the damage zones of outer-rise faults by more than a factor of 10 also over-predicts observed porosity distributions at depths greater than 3 km. However, the simulations consistently underpredict pore fluid distributions between depths of 1–3 km, which may reflect nonlinear elastic or anelastic deformation processes not included in the simulations or inherent uncertainty in the conversion of resistivity measurements to porosity distributions. The best-fitting depth-permeability relationship is relatively low, which has important implications for the timing of pore fluid release down-dip of the trench, as it restricts shallow fluid release, likely leading to higher pore fluid pressures, more widespread hydraulic fracturing, and deeper release of pore fluid along the plate interface.

Journal of Geophysical Research, Solid Earth↗

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↗

Linking morphodynamic response with sediment mass balance on the Colorado River in Marble Canyon: issues of scale, geomorphic setting, and sampling design

Measurements of morphologic change are often used to infer sediment mass balance. Such measurements may, however, result in gross errors when morphologic changes over short reaches are extrapolated to predict changes in sediment mass balance for long river segments. This issue is investigated by examination of morphologic change and sediment influx and efflux for a 100 km segment of the Colorado River in Grand Canyon, Arizona. For each of four monitoring intervals within a 7 year study period, the direction of sand-storage response within short morphologic monitoring reaches was consistent with the flux-based sand mass balance. Both budgeting methods indicate that sand storage was stable or increased during the 7 year period. Extrapolation of the morphologic measurements outside the monitoring reaches does not, however, provide a reasonable estimate of the magnitude of sand-storage change for the 100 km study area. Extrapolation results in large errors, because there is large local variation in site behavior driven by interactions between the flow and local bed topography. During the same flow regime and reach-average sediment supply, some locations accumulate sand while others evacuate sand. The interaction of local hydraulics with local channel geometry exerts more control on local morphodynamic response than sand supply over an encompassing river segment. Changes in the upstream supply of sand modify bed responses but typically do not completely offset the effect of local hydraulics. Thus, accurate sediment budgets for long river segments inferred from reach-scale morphologic measurements must incorporate the effect of local hydraulics in a sampling design or avoid extrapolation altogether.

Arizona↗

Measurement of flow under ice covers in North America

A substantial proportion of natural streams in the United States and Canada are affected by ice cover during the winter. To substantiate the currently used procedures for measuring streamflow during the winter, the U.S. Geological Survey (USGS) and Water Survey of Canada (WSC) began independent, coordinated programs for research and development related to the measurement of streamflow under an ice cover. Detailed measurements of vertical velocity profiles under ice covers in field settings were collected by each agency in accordance with standardized guidelines. The data were then compiled into a joint database. This paper presents a description of the two measurement programs, describes the structure and format of the joint database, and provides preliminary summaries of the data. Ongoing research efforts by the USGS and WSC are described briefly to give examples of the use of the joint database.

Journal of Hydraulic Engineering↗

A fault constitutive relation accounting for thermal pressurization of pore fluid

The heat generated in a slip zone during an earthquake can raise fluid pressure and thereby reduce frictional resistance to slip. The amount of fluid pressure rise depends on the associated fluid flow. The heat generated at a given time produces fluid pressure that decreases inversely with the square root of hydraulic diffusivity times the elapsed time. If the slip velocity function is crack-like, there is a prompt fluid pressure rise at the onset of slip, followed by a slower increase. The stress drop associated with the prompt fluid pressure rise increases with rupture propagation distance. The threshold propagation distance at which thermally induced stress drop starts to dominate over frictionally induced stress drop is proportional to hydraulic diffusivity. If hydraulic diffusivity is 0.02 m 2 /s, estimated from borehole samples of fault zone material, the threshold propagation distance is 300 m. The stress wave in an earthquake will induce an unknown amount of dilatancy and will increase hydraulic diffusivity, both of which will lessen the fluid pressure effect. Nevertheless, if hydraulic diffusivity is no more than two orders of magnitude larger than the laboratory value, then stress drop is complete in large earthquakes.

Journal of Geophysical Research B: Solid Earth↗