Search USGSSearch

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 109 records · Page 6Linked to original sources

Modeling tritium transport through a deep unsaturated zone in an arid environment

Understanding transport of tritium ( 3 H) in unsaturated zones is critical to evaluating options for waste isolation. Tritium typically is a large component of low-level radioactive waste (LLRW). Studies at the U.S. Geological Survey's Amargosa Desert Research Site (ADRS) in Nevada investigate 3 H transport from a closed LLRW facility. Two boreholes are 100 and 160 m from the nearest waste trench and extend to the water table at 110 m. Soil-water vapor samples from the deep boreholes show elevated levels of 3 H at all depths. The objectives of this study were to (i) test source thermal and gas-advection mechanisms driving 3 H transport and (ii) evaluate model sensitivity to these mechanisms and to selected physical and hydraulic properties including porosity, tortuosity, and anisotropy. A two-dimensional numerical model incorporated a non-isothermal, heterogeneous domain of the unsaturated zone and instantaneous isotopic equilibrium. The TOUGH2 code was used; however, it required modification to account for temperature dependence of both the Henry's law equilibrium constant and isotopic fractionation with respect to tritiated water. Increases in source temperature, pressure, and porosity enhanced 3 H migration, but failed to match measured 3 H distributions. All anisotropic simulations with a source pressure component resembled, in shape, the upper portion of the 3 H distribution of the nearest borehole. Isotopic equilibrium limited migration of 3 H, while effects of radioactive decay were negligible. A 500 Pa pressure increase above ambient pressure in conjunction with a high degree of anisotropy (1:100) was necessary for simulated 3 H transport to reach the nearest borehole.

Vadose Zone Journal

The processes of preferential flow in the unsaturated zone

Preferential flow, a major influence in unsaturated soil and rock almost everywhere, occurs by multiple phenomenologically distinct hydraulic processes. For the mode known as funneled flow, concentrated in particularly conductive portions of the medium, the surface-tension/viscous-flow processes of traditional unsaturated flow theory predominate. Fingered flow, through conductive paths of higher water content than surrounding material, requires amendments to traditional theory concerning instabilities and dynamic flow-regime boundaries. Macropore flow, the most recognized preferential flow mode, poses unanswered questions and major difficulties in practice. Accumulated evidence shows that water flows preferentially mostly through macropores that are (a) only partially filled with water, and (b) surrounded by matrix material that is drier, sometimes much drier, than saturation. With partial filling, geometric characteristics such as aperture have much less influence than was previously thought, and the intra-macropore configuration of the flowing water phase, about which little is conclusively known, is then a dominant controlling influence. With unsaturated surroundings, macropore/matrix exchange interactions control, for given input and medium, the initiating circumstances, conveyed flux, and duration of macropore flow. The multiple processes in play during such interactions have different sensitivities to the matrix water state and different directions of influence. The net influence of matrix water content on macropore flow is thus highly complex and a major research need. Additional high-priority topics are: flowpath connectivity, for watersheds as well as small scales; intra-macropore processes, to discern their importance and possible means of quantification; and the identification of measurable soil and rock properties that can be utilized predictively.

Soil Science Society of America Journal

Maximum magnitude earthquakes induced by fluid injection

Analysis of numerous case histories of earthquake sequences induced by fluid injection at depth reveals that the maximum magnitude appears to be limited according to the total volume of fluid injected. Similarly, the maximum seismic moment seems to have an upper bound proportional to the total volume of injected fluid. Activities involving fluid injection include (1) hydraulic fracturing of shale formations or coal seams to extract gas and oil, (2) disposal of wastewater from these gas and oil activities by injection into deep aquifers, and (3) the development of enhanced geothermal systems by injecting water into hot, low-permeability rock. Of these three operations, wastewater disposal is observed to be associated with the largest earthquakes, with maximum magnitudes sometimes exceeding 5. To estimate the maximum earthquake that could be induced by a given fluid injection project, the rock mass is assumed to be fully saturated, brittle, to respond to injection with a sequence of earthquakes localized to the region weakened by the pore pressure increase of the injection operation and to have a Gutenberg-Richter magnitude distribution with a b value of 1. If these assumptions correctly describe the circumstances of the largest earthquake, then the maximum seismic moment is limited to the volume of injected liquid times the modulus of rigidity. Observations from the available case histories of earthquakes induced by fluid injection are consistent with this bound on seismic moment. In view of the uncertainties in this analysis, however, this should not be regarded as an absolute physical limit.

Journal of Geophysical Research B: Solid Earth

The importance of diverse data types to calibrate a watershed model of the Trout Lake Basin, Northern Wisconsin, USA

As part of the USGS Water, Energy, and Biogeochemical Budgets project and the NSF Long-Term Ecological Research work, a parameter estimation code was used to calibrate a deterministic groundwater flow model of the Trout Lake Basin in northern Wisconsin. Observations included traditional calibration targets (head, lake stage, and baseflow observations) as well as unconventional targets such as groundwater flows to and from lakes, depth of a lake water plume, and time of travel. The unconventional data types were important for parameter estimation convergence and allowed the development of a more detailed parameterization capable of resolving model objectives with well-constrained parameter values. Independent estimates of groundwater inflow to lakes were most important for constraining lakebed leakance and the depth of the lake water plume was important for determining hydraulic conductivity and conceptual aquifer layering. The most important target overall, however, was a conventional regional baseflow target that led to correct distribution of flow between sub-basins and the regional system during model calibration. The use of an automated parameter estimation code: (1) facilitated the calibration process by providing a quantitative assessment of the model's ability to match disparate observed data types; and (2) allowed assessment of the influence of observed targets on the calibration process. The model calibration required the use of a 'universal' parameter estimation code in order to include all types of observations in the objective function. The methods described in this paper help address issues of watershed complexity and non-uniqueness common to deterministic watershed models. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Hydrology

The chemistry and isotopic composition of waters in the low-enthalpy geothermal system of Cimino-Vico Volcanic District, Italy

Geothermal energy exploration is based in part on interpretation of the chemistry, temperature, and discharge rate of thermal springs. Here we present the major element chemistry and the δD, δ 18 O, 87 Sr/ 86 Sr and δ 11 B isotopic ratio of groundwater from the low-enthalpy geothermal system near the city of Viterbo in the Cimino-Vico volcanic district of west-Central Italy. The geothermal system hosts many thermal springs and gas vents, but the resource is still unexploited. Water chemistry is controlled by mixing between low salinity,HCO 3 -rich fresh waters (<24.2°C) flowing in shallow volcanic rocks and SO 4 -rich thermal waters (25.3°C to 62.2°C) ascending from deep, high permeability Mesozoic limestones. The (equivalent) SO 4 /Cl (0.01–0.02), Na/Cl (2.82–5.83) and B/Cl ratios (0.02–0.38) of thermal waters differs from the ratios in other geothermal systems from Central Italy, probably implying a lack of hydraulic continuity across the region. The δ 18 O (−6.6‰ to −5.9‰) and δD (−40.60‰ to −36.30‰) isotopic composition of spring water suggest that the recharge area for the geothermal system is the summit region of Mount Cimino. The strontium isotope ratios ( 87 Sr/ 86 Sr) of thermal waters (0.70797–0.70805) are consistent with dissolution of the Mesozoic evaporite-carbonate units that constitute the reservoir, and the ratios of cold fresh waters mainly reflect shallow circulation through the volcanic cover and some minor admixture (<10%) of thermal water as well. The boron isotopic composition (δ 11 B) of fresh waters (−5.00 and 6.12‰) is similar to that of the volcanic cover, but the δ 11 B of thermal waters (−8.37‰ to −4.12‰) is a mismatch for the Mesozoic reservoir rocks and instead reflects dissolution of secondary boron minerals during fluid ascent through flysch units that overlie the reservoir. A slow and tortuous ascent enhances extraction of boron but also promotes conductive cooling, partially masking the heat present in the reservoir. Overall data from this study is consistent with previous studies that concluded that the geothermal system has a large energy potential.

Journal of Volcanology and Geothermal Research

Using borehole geophysics and cross-borehole flow testing to define hydraulic connections between fracture zones in bedrock aquifers

Nearly a decade of intensive geophysical logging at fractured rock hydrology research sites indicates that geophysical logs can be used to identify and characterize fractures intersecting boreholes. However, borehole-to-borehole flow tests indicate that only a few of the apparently open fractures found to intersect boreholes conduct flow under test conditions. This paper presents a systematic approach to fracture characterization designed to define the distribution of fractures along boreholes, relate the measured fracture distribution to structure and lithology of the rock mass, and define the nature of fracture flow paths across borehole arrays. Conventional electrical resistivity, gamma, and caliper logs are used to define lithology and large-scale structure. Borehole wall image logs obtained with the borehole televiewer are used to give the depth, orientation, and relative size of fractures in situ. High-resolution flowmeter measurements are used to identify fractures conducting flow in the rock mass adjacent to the boreholes. Changes in the flow field over time are used to characterize the hydraulic properties of fracture intersections between boreholes. Application of this approach to an array of 13 boreholes at the Mirror Lake, New Hamsphire site demonstrates that the transient flow analysis can be used to distinguish between fractures communicating with each other between observation boreholes, and those that are hydraulically isolated from each other in the surrounding rock mass. The Mirror Lake results also demonstrate that the method is sensitive to the effects of boreholes on the hydraulic properties of the fractured-rock aquifer. Experiments conducted before and after the drilling of additional boreholes in the array and before and after installation of packers in existing boreholes demonstrate that the presence of new boreholes or the inflation of packers in existing boreholes has a large effect on the measured hydraulic properties of the rock mass surrounding the borehole array.

Journal of Applied Geophysics

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

Three Mw ≥ 4.7 earthquakes within the Changning (China) shale gas field ruptured shallow faults intersecting with hydraulic fracturing wells

From 2017 to 2019, three destructive earthquakes (27 January 2017 Mw 4.7, 16 December 2018 Mw 5.2, and 3 January 2019 Mw 4.8) occurred in the Changning shale gas field in the southwest Sichuan Basin, China. Previous seismological studies attributed these events to hydraulic fracturing (HF), but were unable to identify the causative seismogenic faults and their slip behaviors. Here, we use Sentinel-1 synthetic aperture radar data to measure surface deformation triggered by the three events and conduct geodetic inversions to characterize their rupture models. The resulting coseismic interferograms show prominent surface deformation with the maximum line-of-sight displacements of up to 4 cm. The inversion results show that all three earthquakes mainly ruptured sedimentary formations above the shale gas bed, in the upper 3 km of the crust, with slip magnitudes ranging from 8.5 to 15 cm, and stress drops ranging from ∼1.8 to ∼3.3 MPa. Their source faults intersect with horizontal HF wells, but do not root in the crystalline basement. Combined with the reported difficulty of increasing HF operation pressures prior to the three events, we argue that they were most likely induced by direct injection of pressurized fluids into the fault zones. Crustal deformation patterns inferred from regional topography and GPS velocities highlight that the Changning field is located within a triple junction region near the southeastern margin of the Tibetan Plateau with large deformation gradients; such conditions are not only favorable to the development of critically stressed faults, but also facilitate the occurrence of at least moderate magnitude earthquakes.

Changning shale gas field, Sichuan Basin

Evaluation of an experimental LiDAR for surveying a shallow, braided, sand-bedded river

Reaches of a shallow (<1.0m), braided, sand-bedded river were surveyed in 2002 and 2005 with the National Aeronautics and Space Administration's Experimental Advanced Airborne Research LiDAR (EAARL) and concurrently with conventional survey-grade, real-time kinematic, global positioning system technology. The laser pulses transmitted by the EAARL instrument and the return backscatter waveforms from exposed sand and submerged sand targets in the river were completely digitized and stored for postflight processing. The vertical mapping accuracy of the EAARL was evaluated by comparing the ellipsoidal heights computed from ranging measurements made using an EAARL terrestrial algorithm to nearby (<0.5m apart) ground-truth ellipsoidal heights. After correcting for apparent systematic bias in the surveys, the root mean square error of these heights with the terrestrial algorithm in the 2002 survey was 0.11m for the 26 measurements taken on exposed sand and 0.18m for the 59 measurements taken on submerged sand. In the 2005 survey, the root mean square error was 0.18m for 92 measurements taken on exposed sand and 0.24m for 434 measurements on submerged sand. In submerged areas the waveforms were complicated by reflections from the surface, water column entrained turbidity, and potentially the riverbed. When applied to these waveforms, especially in depths greater than 0.4m, the terrestrial algorithm calculated the range above the riverbed. A bathymetric algorithm has been developed to approximate the position of the riverbed in these convolved waveforms and preliminary results are encouraging. ?? 2007 ASCE.

Journal of Hydraulic Engineering

Predicted pH at the domestic and public supply drinking water depths, Central Valley, California

This scientific investigations map is a product of the U.S. Geological Survey (USGS) National Water-Quality Assessment (NAWQA) project modeling and mapping team. The prediction grids depicted in this map are of continuous pH and are intended to provide an understanding of groundwater-quality conditions at the domestic and public supply drinking water zones in the groundwater of the Central Valley of California. The chemical quality of groundwater and the fate of many contaminants is often influenced by pH in all aquifers. These grids are of interest to water-resource managers, water-quality researchers, and groundwater modelers concerned with the occurrence of natural and anthropogenic contaminants related to pH. In this work, the median well depth categorized as domestic supply was 30 meters below land surface, and the median well depth categorized as public supply is 100 meters below land surface. Prediction grids were created using prediction modeling methods, specifically boosted regression trees (BRT) with a Gaussian error distribution within a statistical learning framework within the computing framework of R ( http://www.r-project.org/ ). The statistical learning framework seeks to maximize the predictive performance of machine learning methods through model tuning by cross validation. The response variable was measured pH from 1,337 wells and was compiled from two sources: USGS National Water Information System (NWIS) database (all data are publicly available from the USGS: http://waterdata.usgs.gov/ca/nwis/nwis ) and the California State Water Resources Control Board Division of Drinking Water (SWRCB-DDW) database (water quality data are publicly available from the SWRCB: http://www.waterboards.ca.gov/gama/geotracker_gama.shtml ). Only wells with measured pH and well depth data were selected, and for wells with multiple records, only the most recent sample in the period 1993–2014 was used. A total of 1,003 wells (training dataset) were used to train the BRT model, and 334 wells (hold-out dataset) were used to validate the prediction model. The training r-squared was 0.70, and the root-mean-square error (RMSE) in standard pH units was 0.26. The hold-out r-squared was 0.43, and RMSE in standard pH units was 0.37. Predictor variables consisting of more than 60 variables from 7 sources were assembled to develop a model that incorporates regional-scale soil properties, soil chemistry, land use, aquifer textures, and aquifer hydrology. Previously developed Central Valley model outputs of textures (Central Valley Textural Model, CVTM; Faunt and others, 2010) and MODFLOW-simulated vertical water fluxes and predicted depth to water table (Central Valley Hydrologic Model, CVHM; Faunt, 2009) were used to represent aquifer textures and groundwater hydraulics, respectively. In this work, wells were attributed to predictor variable values in ArcGIS using a 500-meter buffer. Faunt, C.C., ed., 2009, Groundwater availability in the Central Valley aquifer, California: U.S. Geological Survey Professional Paper 1776, 225 p., accessed at https://pubs.usgs.gov/pp/1766/ . Faunt, C.C., Belitz, K., and Hanson, R.T., 2010, Development of a three-dimensional model of sedimentary texture in valley-fill deposits of Central Valley, California, USA: Hydrogeology Journal, v. 18, no. 3, p. 625–649, https://doi.org/10.1007/s10040-009-0539-7 .

California

Inverse modeling of BTEX dissolution and biodegradation at the Bemidji, MN crude-oil spill site

The U.S. Geological Survey (USGS) solute transport and biodegradation code BIOMOC was used in conjunction with the USGS universal inverse modeling code UCODE to quantify field-scale hydrocarbon dissolution and biodegradation at the USGS Toxic Substances Hydrology Program crude-oil spill research site located near Bemidji, MN. This inverse modeling effort used the extensive historical data compiled at the Bemidji site from 1986 to 1997 and incorporated a multicomponent transport and biodegradation model. Inverse modeling was successful when coupled transport and degradation processes were incorporated into the model and a single dissolution rate coefficient was used for all BTEX components. Assuming a stationary oil body, we simulated benzene, toluene, ethylbenzene, m , p -xylene, and o -xylene (BTEX) concentrations in the oil and ground water, respectively, as well as dissolved oxygen. Dissolution from the oil phase and aerobic and anaerobic degradation processes were represented. The parameters estimated were the recharge rate, hydraulic conductivity, dissolution rate coefficient, individual first-order BTEX anaerobic degradation rates, and transverse dispersivity. Results were similar for simulations obtained using several alternative conceptual models of the hydrologic system and biodegradation processes. The dissolved BTEX concentration data were not sufficient to discriminate between these conceptual models. The calibrated simulations reproduced the general large-scale evolution of the plume, but did not reproduce the observed small-scale spatial and temporal variability in concentrations. The estimated anaerobic biodegradation rates for toluene and o -xylene were greater than the dissolution rate coefficient. However, the estimated anaerobic biodegradation rates for benzene, ethylbenzene, and m , p -xylene were less than the dissolution rate coefficient. The calibrated model was used to determine the BTEX mass balance in the oil body and groundwater plume. Dissolution from the oil body was greatest for compounds with large effective solubilities (benzene) and with large degradation rates (toluene and o -xylene). Anaerobic degradation removed 77% of the BTEX that dissolved into the water phase and aerobic degradation removed 17%. Although goodness-of-fit measures for the alternative conceptual models were not significantly different, predictions made with the models were quite variable.

Minnesota

Ice sheet load cycling and fluid underpressures in the Eastern Michigan Basin, Ontario, Canada

Strong fluid underpressures have been detected in Paleozoic strata in the eastern Michigan Basin, with hydraulic heads reaching ~400&thinsp;m below land surface (~4&thinsp;MPa underpressure) and ~200&thinsp;m below sea level in strata where unusually low permeabilities (~10 &minus;20 &ndash;10 &minus;23 &thinsp;m 2 ) were measured in situ. Multiple glaciations, including three with as much as 3&thinsp;km of ice cover at the site in the last 120&thinsp;ka, suggest a causal link with the underpressures. We examined this possibility using a one-dimensional groundwater flow model incorporating mechanical loading from both ice weight and lithospheric flexure. Because hydrologic and mechanical changes during glaciation are not well characterized and subsurface properties are imperfectly known, the model was used inversely to estimate flexural loads and loosely constrained permeabilities by matching observed pressures. Acceptable matches were obtained for a surprisingly wide range of scenarios with permeabilities close to measured values and plausible flexural loads. Matches were not obtained when too many parameters were preselected, or when permeabilities were constrained to be significantly larger than measured values. In successful model runs groundwater expulsion under glacial-mechanical loads caused the underpressuring, and flexural loads were important if aquifer and sub-glacial pressures were significantly elevated during glaciation. Simulated fluid pressures in the low-permeability strata fluctuated by 30&ndash;40&thinsp;MPa during glacial cycles but resulted in advective transport of only tens of meters or less. Although other mechanisms cannot be ruled out, we conclude that glacial-mechanical forcing of a water-saturated system can explain the observed underpressures.

Ontario

Impacts of spontaneous waterfall development on bedrock river longitudinal profile morphology

River profiles are shaped by climatic and tectonic history, lithology, and internal feedbacks between flow hydraulics, sediment transport and erosion. In steep channels, waterfalls may self-form without changes in external forcing (i.e., autogenic formation) and erode at rates faster or slower than an equivalent channel without waterfalls. We use a 1-D numerical model to investigate how self-formed waterfalls alter the morphology of bedrock river longitudinal profiles. We modify the standard stream power model to include a slope threshold above which waterfalls spontaneously form and a rate constant allowing waterfalls to erode faster or slower than other fluvial processes. Using this model, we explore how waterfall formation alters both steady state and transient longitudinal profile forms. Our model predicts that fast waterfalls create km-scale reaches in a dynamic equilibrium with channel slope held approximately constant at the threshold slope for waterfall formation, while slow waterfalls can create local channel slope maxima at the location of slow waterfall development. Furthermore, slow waterfall profiles integrate past base level histories, leading to multiple possible profile forms, even at steady-state. Consistency between our model predictions and field observations of waterfall-rich rivers in the Kings and Kaweah drainages in the southern Sierra Nevada, California, supports the hypothesis that waterfall formation can modulate river profiles in nature. Our findings may help identify how bedrock channels are influenced by waterfall erosion and aid in distinguishing between signatures of external and internal perturbations, thereby strengthening our ability to interpret past climate and tectonic changes from river longitudinal profiles.

Journal of Geophysical Research - Earth Surface

On extracting sediment transport information from measurements of luminescence in river sediment

Accurately quantifying sediment transport rates in rivers remains an important goal for geomorphologists, hydraulic engineers, and environmental scientists. However, current techniques for measuring long-time scale (102–106 years) transport rates are laborious, and formulae to predict transport are notoriously inaccurate. Here we attempt to estimate sediment transport rates by using luminescence, a property of common sedimentary minerals that is used by the geoscience community for geochronology. This method is advantageous because of the ease of measurement on ubiquitous quartz and feldspar sand. We develop a model from first principles by using conservation of energy and sediment mass to explain the downstream pattern of luminescence in river channel sediment. We show that the model can accurately reproduce the luminescence observed in previously published field measurements from two rivers with very different sediment transport styles. The model demonstrates that the downstream pattern of river sand luminescence should show exponential-like decay in the headwaters which asymptotes to a constant value with further downstream distance. The parameters from the model can then be used to estimate the time-averaged virtual velocity, characteristic transport lengthscale, storage time scale, and floodplain exchange rate of fine sand-sized sediment in a fluvial system. The sediment transport values predicted from the luminescence method show a broader range than those reported in the literature, but the results are nonetheless encouraging and suggest that luminescence demonstrates potential as a sediment transport indicator. However, caution is warranted when applying the model as the complex nature of sediment transport can sometimes invalidate underlying simplifications.

Journal of Geophysical Research F: Earth Surface

In situ stress, natural fracture distribution, and borehole elongation in the Auburn Geothermal Well, Auburn, New York

Hydraulic fracturing stress measurements and a borehole televiewer survey were conducted in a 1.6‐km‐deep well at Auburn, New York. This well, which was drilled at the outer margin of the Appalachian Fold and Thrust Belt in the Appalachian Plateau, penetrates approximately 1540 m of lower Paleozoic sedimentary rocks and terminates 60 m into the Precambrian marble basement. Analysis of the hydraulic fracturing tests indicates that the minimum horizontal principal stress increases in a nearly linear fashion from 9.9±0.2 MPa at 593 m to 30.6±0.4 MPa at 1482 m. The magnitude of the maximum horizontal principal stress increases in a less regular fashion from 13.8±1.2 MPa to 49.0±2.0 MPa over the same depth range. The magnitudes of the horizontal principal stresses relative to the calculated overburden stress are somewhat lower than is the norm for this region and are indicative of a strike‐slip faulting regime that, at some depths, is transitional to normal faulting. As expected from the relative aseismicity of central New York State, however, analysis of the magnitudes of the horizontal principal stresses indicates, at least to a depth of 1.5 km, that frictional failure on favorably oriented preexisting fault planes is unlikely. Orientations of the hydraulic fractures at 593 and 919 m indicate that the azimuth of the maximum horizontal principal stress at Auburn is N83°E±15°, in agreement with other stress field indicators for this region. The borehole televiewer log revealed a considerable number of planar features in the Auburn well, the great majority of which are subhorizontal (dips < 5°) and are thought to be bedding plane washouts or drill bit scour marks. In addition, a smaller number of distinct natural fractures were observed on the borehole televiewer log. Of these, the distinct steeply dipping natural fractures in the lower half of the sedimentary section at Auburn tend to strike approximately east‐west, while those in the upper part of the well and in the Precambrian basement exhibit no strong preferred orientation. The origin of this east‐west striking fracture set is uncertain, as it is parallel both to the contemporary direction of maximum horizontal compression and to a late Paleozoic fracture set that has been mapped to the south of Auburn. In addition to these planar features the borehole televiewer log indicates paired dark bands on diametrically opposite sides of the borehole throughout the Auburn well. Processing of the borehole televiewer data in the time domain revealed these features to be irregular depressions in the borehole wall. As these depressions were consistently oriented in a direction at right angles to the direction of maximum horizontal compression, we interpret them to be the result of stress‐induced spalling of the borehole wall (breakouts).

New York

Sediment supply versus local hydraulic controls on sediment transport and storage in a river with large sediment loads

The Rio Grande in the Big Bend region of Texas, USA, and Chihuahua and Coahuila, Mexico, undergoes rapid geomorphic changes as a result of its large sediment supply and variable hydrology; thus, it is a useful natural laboratory to investigate the relative importance of flow strength and sediment supply in controlling alluvial channel change. We analyzed a suite of sediment transport and geomorphic data to determine the cumulative influence of different flood types on changing channel form. In this study, physically based analyses suggest that channel change in the Rio Grande is controlled by both changes in flow strength and sediment supply over different spatial and temporal scales. Channel narrowing is primarily caused by substantial deposition of sediment supplied to the Rio Grande during tributary-sourced flash floods. Tributary floods have large suspended-sediment concentrations, occur for short durations, and attenuate rapidly downstream in the Rio Grande, depositing much of their sediment in downstream reaches. Long-duration floods on the mainstem have the capacity to enlarge the Rio Grande, and these floods, released from upstream dams, can either erode or deposit sediment in the Rio Grande depending upon the antecedent in-channel sediment supply and the magnitude and duration of the flood. Geomorphic and sediment transport analyses show that the locations and rates of sand erosion and deposition during long-duration floods are most strongly controlled by spatial changes in flow strength, largely through changes in channel slope. However, spatial differences in the in-channel sediment supply regulate sediment evacuation or accumulation over time in long reaches (greater than a kilometer).

Big Bend National Park

Non-double-couple microearthquakes at Long Valley caldera, California, provide evidence for hydraulic fracturing

Most of 26 small (0.4≲ M ≲3.1) microearthquakes at Long Valley caldera in mid-1997, analyzed using data from a dense temporary network of 69 digital three-component seismometers, have significantly non-double-couple focal mechanisms, inconsistent with simple shear faulting. We determined their mechanisms by inverting P - and S -wave polarities and amplitude ratios using linear-programming methods, and tracing rays through a three-dimensional Earth model derived using tomography. More than 80% of the mechanisms have positive (volume increase) isotropic components and most have compensated linear-vector dipole components with outward-directed major dipoles. The simplest interpretation of these mechanisms is combined shear and extensional faulting with a volume-compensating process, such as rapid flow of water, steam, or CO 2 into opening tensile cracks. Source orientations of earthquakes in the south moat suggest extensional faulting on ESE-striking subvertical planes, an orientation consistent with planes defined by earthquake hypocenters. The focal mechanisms show that clearly defined hypocentral planes in different locations result from different source processes. One such plane in the eastern south moat is consistent with extensional faulting, while one near Casa Diablo Hot Springs reflects en echelon right-lateral shear faulting. Source orientations at Mammoth Mountain vary systematically with location, indicating that the volcano influences the local stress field. Events in a ‘spasmodic burst’ at Mammoth Mountain have practically identical mechanisms that indicate nearly pure compensated tensile failure and high fluid mobility. Five earthquakes had mechanisms involving small volume decreases, but these may not be significant. No mechanisms have volumetric moment fractions larger than that of a force dipole, but the reason for this fact is unknown.

California

Passage and behaviour of cultured Lake Sturgeon in a prototype side-baffle fish ladder: I. Ladder hydraulics and fish ascent

Research and development of a fish ladder for sturgeons requires understanding ladder hydraulics and sturgeon behaviour in the ladder to insure the ladder is safe and provides effective passage. After years of research and development, we designed and constructed a full‐scale prototype side‐baffle ladder inside a spiral flume (38.3 m long × 1 m wide × 1 m high) on a 6% (1 : 16.5) slope with a 1.92‐m rise in elevation (bottom to top) to test use by sturgeons. Twenty‐eight triangular side baffles, each extending part way across the flume, alternated from inside wall to outside wall down the ladder creating two major flow habitats: a continuous, sinusoidal flow down the ladder through the vertical openings of side‐baffles and an eddy below each side baffle. Ascent and behaviour was observed on 22 cultured Lake Sturgeon = LS ( Acipenser fulvescens ) repeatedly tested in groups as juveniles (as small as 105.1 cm TL, mean) or as adults (mean TL, 118 cm) during four periods (fall 2002 and 2003; spring 2003 and 2007). Percent of juveniles entering the ladder that ascended to the top was greater in spring (72.7%) than in fall (40.9–45.5%) and 90.9% of 11 adults, which ascended as juveniles, ascended to the top. Six LS (27.3%) never swam to the top and seven (31.8%) swam to the top in all tests, indicating great variability among individuals for ascent drive. Some LS swam directly to the top in <1 min, but most rested in an eddy during ascent. Juveniles swimming through outside wall baffle slots (mean velocity, 1.2 m s −1 ) swam at 1.8–2.2 body lengths s −1 and 3.2–3.3 tail beats s −1 , either at or approaching prolonged swimming speed. The side‐baffle ladder was stream‐like and provided key factors for a sturgeon ladder: a continuous flow and no full cross‐channel walls, abundant eddies for resting, an acceptable water depth, and a water velocity fish could ascend swimming 2 bl s −1 . A side‐baffle ladder passes LS and other moderate‐swimming fishes.

Journal of Applied Ichthyology