USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Hydrology”

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.

1,659 records · Page 33Linked to original sources
USGS2025

Factors affecting the distribution of water-bearing fractures in the bedrock aquifers of West Virginia

Bedrock aquifers cover 23,601 square miles within the State of West Virginia and comprise 97.4 percent of the surficial area within the State; the remaining 2.6 percent (621 square miles) consists of alluvial sand-and-gravel and glacial outwash aquifers bordering the State’s major rivers. While West Virginia’s alluvial aquifers have been studied extensively, bedrock aquifers have only been characterized for studies completed in a few areas in Jefferson, McDowell, and Monroe Counties. Bedrock aquifers are water supplies for public supply, agriculture, industry, and residential homeowner use. In this study, the U.S. Geological Survey, in cooperation with the West Virginia Department of Environmental Protection Division of Water and Waste Management, provides a statewide assessment of the occurrence and distribution of fractures within bedrock aquifers of the State and the various topographic, physiographic, and lithologic influences controlling the occurrence and distribution of bedrock fractures. The results of this study provide an increased understanding of the distribution of fractures in bedrock aquifers in West Virginia and help to verify trends that have been suspected for many years but were never well documented or verified by data. The results confirmed that the density of fractures and those that were determined to be water bearing decrease significantly with depth. A statistically significant difference in the density of fractures was observed at a depth of 215 feet for wells in the Appalachian Plateaus Physiographic Province’s and in the Valley and Ridge Physiographic Province’s aquifers; a higher density of fractures and water-bearing fractures were above a depth of 215 feet than below that depth. This is an important consideration when drilling wells for residential, commercial, industrial, or agricultural water supply. Abandoned underground coal mines are commonly believed to form large pools of water in the interconnected mine entries in abandoned room and pillar coal mines. Such pools of water can and do exist in abandoned underground coal mines, but many mines lack open entries and are held up by overburden strata and pillars that can collapse and form aquifers comprised of vast interconnected rubble zones (gob), especially in older mines. Data assessed for this study showed that shale-corrected values of effective porosity for limestone aquifers in West Virginia had a median value of 2 percent and an average value of 4 percent and generally are mineralized with low effective porosity. Argillaceous or sandy limestone has a median shale-corrected porosity of 4 percent and an average shale-corrected porosity of 5 percent. The median and average shale-corrected porosity of sandstone aquifers was estimated to be 14 percent, but the median shale-corrected porosity for argillaceous or calcareous sandstone was 5 percent and the average shale-corrected porosity for argillaceous or calcareous sandstone was 6 percent. Even though shale has a relatively high total sonic porosity compared to other lithologies, shale and siltstone had relatively low shale-corrected porosity, ranging from 0 to 2 percent. Well yields were previously documented to be highest in valley settings, lowest on hilltops, and intermediate on hillsides. Transmissivity data provided by this study confirm this general pattern within the Appalachian Plateaus Province; however, the Valley and Ridge Province does not follow this pattern. While still lowest on hilltop settings, the highest well yields were in hillside settings. The trend for the Valley and Ridge Province was likely skewed because of 9 high-yield wells specifically targeting deeper thin limestone units, such as the Tonoloway and Helderberg Limestones, at depths with transmissivity in excess of 2,000 feet squared per day in Mineral County, West Virginia, or targeting karst aquifers in Berkeley, Jefferson, or Greenbrier Counties, West Virginia. Finally, water-bearing fractures have been hypothesized to comprise a small number of all fractures within a typical bedrock well in West Virginia. Data collected for this study support this theory. A total of 3,403 fractures were identified during this study; 3,151 (92.6 percent) of those fractures are low-transmissive, and only 252 (7.4 percent) fractures are water-bearing. Even though a well may contain many fractures, less than 8 percent are considered water-bearing fractures.

West Virginia↗
USGS1976

Detailed investigation of a coastal plain-Piedmont fault contact in northeastern Virginia

A trench across the Dumfries fault of the Stafford fault zone near Stafford, Virginia (see Mixon and Newell, 1976) exposed a high-angle reverse fault juxtaposing Paleozoic Quantico Slate and Cretaceous Potomac Group rocks. A system of subsidiary reverse faults and small scale normal faults disrupt the Coastal Plain strata adjacent to the main fault. Cross-cutting relations and stratigraphic details indicate that the subsidiary faults in the Cretaceous rocks document a history of episodic deformation, in part active during accumulation of the Coastal Plain strata. Slickensides and stereo net analysis of theorientation of subsidiary reverse faults, normal faults, and bedding indicate that, through time, faulting consistently included right lateral separation.

Virginia↗
USGS2025

Global maps of critical mineral production in 2023

Introduction The global production of many mineral commodities, especially critical minerals, is concentrated in a few countries that have mineral resources and the infrastructure necessary to mine and process those resources. For this reason, the type and amount of mineral production differ by country. For example, many countries produce such metallic ores as gold and silver, whereas only a few countries produce magnesium, niobium, platinum-group metals, and rare earths. The concentration of mining and processing in certain countries necessitates the existence of a global supply chain. A mineral supply chain is the sequence of mining and processing of minerals and manufacturing of products. Mineral supply chains are global in scale, complex, and dynamic. Supply chain data can be used to understand how a country’s mineral resources and various economic, technical, and environmental factors affect the complexity of global supply chains. This fact sheet summarizes the world’s leading countries (those accounting for 5 percent or more of a commodity’s global production in 2023) for production of select mineral commodities (mainly critical minerals) in the mining and processing stages. These countries and the minerals they produce are synthesized on global maps to communicate the status of, and potential risk to, mineral commodity supply chains from geographic production concentration. Trade data from United Nations Statistics Division (2025) is used to support assessments of the observed production data.

Fact Sheet↗
USGS2025

Physicochemical properties and bioreactivity of sub-10 µm geogenic particles: Comparison of volcanic ash and desert dust

Exposure to ambient particulate matter (PM) with an aerodynamic diameter of <10 μm (PM 10 ) is a well-established health hazard. There is increasing evidence that geogenic (Earth-derived) particles can induce adverse biological effects upon inhalation, though there is high variability in particle bioreactivity that is associated with particle source and physicochemical properties. In this study, we investigated physicochemical properties and biological reactivity of volcanic ash from the April 2021 eruption of La Soufrière volcano, St. Vincent, and two desert dust samples: a standardized test dust from Arizona and an aeolian Gobi Desert dust sampled in China. We determined particle size, morphology, mineralogy, surface texture and chemistry in sub-10 μm material to investigate associations between particle physicochemical properties and observed bioreactivity. We assessed cellular responses (cytotoxic and pro-inflammatory effects) to acute particle exposures (24 hr) in monocultures at the air-liquid interface using two types of cells of the human airways: BEAS-2B bronchial epithelial cells and A549 alveolar type II epithelial cells. In acellular assays, we also assessed particle oxidative potential and the presence of microorganisms. The results showed that volcanic ash and desert dust exhibit intrinsically different particle morphology, surface textures and chemistry, and variable mineralogical content. We found that Gobi Desert dust is more bioreactive than freshly erupted volcanic ash and Arizona test dust, which is possibly linked to the presence of microorganisms (bacteria) and/or nanoscale elongated silicate minerals (potentially clay such as illite or vermiculite) on particle surfaces.

GeoHealth↗
USGS1976

A comparison of electronic particle counting and pipette techniques in routine mud analysis

Electronic particle counting and pipette size-analysis techniques are compared for their relative effectiveness in large-scale routine mud analysis. Optimum laboratory procedures for conducting electronic analyses with a Model TA Coulter Counter are presented. Comparative electronic and pipette analyses conducted on the same samples resulted in equivalent precision for both techniques. Results of electronic analyses consistently showed size distributions significantly coarser than those obtained by pipette analyses. The laboratory efficiency of the relatively rapid electronic technique is far superior to that of the pipette technique, although more extensive precautionary measures are required to maintain a high level of analytical accuracy and precision. The quality of electronic analyses largely depends upon the effectiveness of standardized operational procedures established within the laboratory.

Open-File Report↗
USGS2025

Pre-Acadian tectonics of the eastern Orange-Milford Belt, south-central Connecticut

This excursion presents a reinterpretation of mapping and new analytical data from the eastern Orange-Milford belt (OMB) in south-central Connecticut. The OMB is a fault-bound terrane of argillites and mafic rocks of anomalously low metamorphic grade—and of poorly constrained ages and tectonic affinity—wedged between kyanite/sillimanite-grade peri- Laurentian rocks to the west and anatectic peri-Gondwanan rocks to the east (Fig. 1A). Our data demonstrate that Ordovician(?) igneous and sedimentary rocks of the OMB were variably metamorphosed in the Ordovician and Silurian but escaped regionally pervasive, high-grade Devonian and later metamorphism. Previous interpretations (Fritts 1963a, 1965a, 1965b; Burger, 1967; Burger and others, 1968; Rodgers, 1985) described these rocks as a conformable sequence of low-grade, Ordovician to Devonian metasediments and metavolcanics. Our results reveal that the “metavolcanics” are not extrusive rocks but rather slivers of lower oceanic crust with complicated high- and low-grade metamorphic fabrics, intruded by a swarm of Silurian sheeted basalt dikes, and in fault contact with the surrounding metasediments. The purpose of this trip is to show evidence of early Paleozoic (Taconic) deformation and metamorphism preserved in rocks of the eastern OMB. These rocks remained shallow, cool, and sufficiently dry during the regionally dominant Acadian and Alleghanian orogenies to have avoided significant overprinting. As such, these rocks serve as windows into a geologic history otherwise unavailable between anatectic rocks of the peri-Gondwanan Bronson Hill, Avalon, and Gander terranes east of the Hartford basin and sillimanite-grade rocks of the peri-Laurentian Hartland and gneiss dome belts west of the OMB. We present major and trace element geochemistry including rare-earth element patterns for all mafic units in the eastern OMB as well as 40 Ar/ 39 Ar age spectra of amphibole, muscovite, and K-feldspar from rocks of the Maltby Lakes complex (of Deasy and others, 2017), Savin Schist, and Wepawaug Schist. Our evidence demonstrates that the units of the OMB have been assembled by faulting or intrusion, and that no stratigraphic relationships exist between the argillaceous schists and the metaigneous rocks.

Connecticut↗
USGS1976

Diamond drilling at the Ma'milah Mine, Kingdom of Saudi Arabia

The Ma'millah gold mine is about 90 km southeast of At Ta'if at lat 21°03'N., long 41°18'E., in southwestern Saudi Arabia. The deposit, which was worked extensively by ancient miners, consists of several veins and massive lenses of quartz along a regional fault. The extensive ancient mine workings, large size of some quartz outcrops, and significant gold content of several samples indicated that subsurface exploration was warranted. Accordingly, two holes were drilled into the deposit during the interval December 1972 to February 1973. The drill findings were disappointing; the grade of vein material was found to be too low to be of economic interest and the quartz bodies were found to pinch out both along strike and downdip. The drilling program therefore was terminated. The deposit is not worth further exploration

Ma'milah Mine↗
USGS2025

Capturing the uncertainty of seismicity observations in earthquake rate estimates: Implications for probabilistic seismic hazard analysis and the USGS National Seismic Hazard Model

The rate of earthquakes in a region is a fundamental input to Probabilistic Seismic Hazard Analysis. We present a Monte Carlo method for computing that rate from seismicity catalogs while including a range of data and analysis uncertainties. This method is applied to regions for which the b value is assumed to be spatially invariant. Each region is broken down into epochs for which each epoch is estimated to have a uniform magnitude of completeness ( ⁠⁠ M c ). The distribution of earthquake rates for M ≥ M c is determined for each epoch by considering the Poisson likelihood of rates given the number of observed earthquakes with M ≥ M c ⁠ . We use a Monte Carlo process to include the uncertainty in b , ⁠ M c , and individual event magnitudes. The result for each epoch is the joint distribution of the Poisson rate of earthquakes with magnitudes larger than the minimum value used to calculate hazard ( ⁠⁠ M 1 ) and the Gutenberg–Richter b values, which control the extrapolation to other magnitudes. The rate for each region is either the duration‐weighted average over the epochs or, to better capture temporal variations, we also consider mixture models. The mixture models also provide an avenue to allow temporal variations in b values. To implement this joint distribution in a logic tree, we use the mean and 95% confidence branches, each of which is parameterized with an M ≥ M 1 rate and b value. We explore different ways of defining those branches, as well as non‐Gutenberg–Richter branches, and their impact on hazard estimates. The mean hazard, but not the fractiles, is robust with respect to these choices. To illustrate these new methods, we use synthetic data and catalogs from recent U.S. Geological Survey National Seismic Hazard Models for the Central and Eastern United States and for Puerto Rico and the U.S. Virgin Islands.

central and eastern United States, Puerto Rico, U.↗
USGS1995

Water-quality assessment of the Ozark Plateaus study unit, Arkansas, Kansas, Missouri, and Oklahoma: Analysis of information on nutrients, suspended sediment, and suspended solids, 1970-92

Water-quality data collected during water years 1970-90 (October 1 to September 30) for 83 surface-water sites and during 1970-92 for 395 ground-water sites in the 48,000 square mile Ozark Plateaus study unit of the National Water Quality Assessment Program were analyzed using selected descriptive and statistical methods. The water- quality data include nutrient (nitrogen and phosphorus), suspended sediment, and suspended- solids data, and ancillary information such as fertilizer use, animal waste, sewage treatment plant, and land use. Statistically significant differences exist in surface-water quality that can be attributed to physiography, land use, and other effects. The sites that were considered to be substantially affected by sewage treatment plants had the largest concentrations of nutrients. Nutrient concentrations generally were larger at sites associated with agricultural basins than at sites associated with forested basins. Statistically significant differences existed in the quality of ground water that can be attributed to hydrogeologic and land-use effects. Nutrient concentrations generally were largest where the water source is indicated to be shallow in origin and where parts of the hydrogeologic units are in agricultural land-use areas. Water quality has changed at several surface-water sites since 1970. Nutrient concentrations appear to have increased at some sites and decreased at other sites. Causes of these apparent trends are not known, but many of the sites with apparent trends are in agricultural areas. Surface-water loads of nutrients and suspended sediment were affected by several factors including streamflow, climate, drainage area, reservoir operation, and inputs from point and nonpoint sources. Annual loads were largest in large basins, with large inputs of nutrients or sediment during periods of high streamflows at locations where reservoir operation effects are not substantial.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS2025

A 700-year rupture sequence of great eastern Aleutian earthquakes from tsunami modeling of stratigraphic records

Great Aleutian underthrusting earthquakes produced destructive tsunamis impacting Hawaiʻi in 1946 and 1957. Prior modeling of the 1957 tsunami deposit and runup records on eastern Aleutian and Hawaiian Islands jointly with tide-gauge observations across the Pacific Ocean constrained a rupture model with shallow slip up to 26 m along 600 km of the plate boundary. Here we implement this modeling approach to older deposits and show alternating deep and shallow megathrust slip up to 26, 32, and 22 m for great earthquakes along the same segment in the 18 th , 15 th , and 14 th centuries. All three modeled prehistoric Aleutian earthquakes produce tsunami inundation in Hawaiʻi with the most severe, 14 th century event having impacts exceeding the 1957 event. The along-dip variability of these four ruptures spanning seven centuries provides insights on earthquake cycles for engineering design and hazard assessment. The 15 th century and 1957 rupture models provide evidence for recurrence of tsunami earthquakes, which can produce disproportionately large tsunamis for a given moment magnitude due to reduced rigidity in the shallow megathrust. The 14 th and 18 th century events likely ruptured deeper regions that did not slip in 1957, suggesting potential for corresponding deeper failure in the next great eastern Aleutian earthquake.

Alaska, Hawaii↗
USGS1976

Guidelines for collection and field analysis of ground-water samples for selected unstable constituents

The unstable nature of many chemical and physical constituents in ground water requires special collection procedures and field analysis immediately after collection. This report describes the techniques and equipment commonly used in the collection and field analysis of samples for pH, temperature, carbonate, bicarbonate, specific conductance, Eh, and dissolved oxygen.

Open-File Report↗
USGS2024

Paleoseismology of the Sawtooth fault and implications for fault behavior in the epicentral region of the 2020 Mw 6.5 Stanley, Idaho, earthquake

The 2020 moment magnitude ( M w ⁠⁠ ) 6.5 Stanley, Idaho, earthquake raised questions about the history and extent of complex faulting in the northwestern Centennial Tectonic Belt (CTB) and its relation to the Sawtooth normal fault and Eocene Trans‐Challis fault system (TCFS). To explore faulting in this area, we excavated a paleoseismic trench across the Sawtooth fault along the western margin of the CTB, and compared an early Holocene (9.1 ± 2.1 ka, 1δ ⁠ ) rupture at the site with lacustrine paleoseismic data and fault mapping in the 2020 epicentral region. We find: (1) a history of partial to full rupture of the Sawtooth fault ( M w ⁠ 6.8–7.4), (2) that shorter ruptures ( M w ≤ ⁠⁠6.9 ) are likely along distributed and discontinuous faults in the epicentral region, (3) that this complex system that hosted the 2020 earthquake is not directly linked to the Sawtooth fault, (4) that the northeast‐trending TCFS likely plays a role in controlling fault length and rupture continuity for adjacent faults, and (5) that parts of the TCFS may facilitate displacement transfer between normal faults that accommodate crustal extension and rotation. Our results help unravel complex faulting in the CTB and imply that relict structures can help inform regional seismic hazard assessments.

Idaho↗
USGS1976

Transport and dispersion of fluorescent tracer particles for the dune-bed condition, Atrisco Feeder Canal near Bernalillo, New Mexico

A tracer technique in which mineral particles were coated with fluorescent dyes was used to study the rates of transport and dispersion of sediment particles`of various diameters and specific gravities for a dune-bed condition in an alluvial channel. The experiment was conducted in the Atrisco Feeder Canal near Bernalillo, N. Mex., between May 1 and July 14, 1967. A continuous point source of tracers, approximated by injections at 5- or 10-minute intervals, was maintained for 7 days so that the steady-dilution procedure could be used to calculate the transport rate of bed material. After termination of the injection process, the spatial-integration procedure was used to follow the movement of the tracers downstream and to calculate the transport rate. Samples of the bed material in transport and the accompanying tracers moving along the surface of the dune bed were obtained periodically throughout the study with the "dustpan" sampler especially designed for fluorescent tracer studies. In addition, the spatial distributions of the tracers in the dune bed were determined three times during the study by core sampling.

New Mexico↗
USGS2026

High-resolution magnetic survey using an unoccupied aerial vehicle to constrain buried lava flow geometry, volume, and eruptive history of Little Cones, Crater Flat, Nevada

Magnetic surveys are an important tool used to augment geologic mapping in distributed volcanic fields. Using magnetic anomalies, it is possible to model the geometry of shallowly buried volcanic features, such as conduits, sills, and lava flows. This subsurface mapping is important for understanding eruption dynamics and emplacement of lava flows, and it sometimes reveals buried volcanoes no longer visible at the surface. These data are critical to better interpret the numbers, styles, and magnitudes of eruptions in distributed volcanic fields and their associated volcanic hazards. New advances in unoccupied aerial vehicles (UAVs) offer an attractive middle range of resolution and aerial coverage between ground-based magnetic surveys and aeromagnetic surveys. Here, we present the results of a UAV fluxgate magnetic survey of the Little Cones, Nevada, scoria cones, which have been the target of previous ground and aeromagnetic surveys. The magnetic anomalies at Little Cones are of interest because the surrounding alluvium conceals lava flows that erupted from Little Cones, making it very difficult to understand the volume and morphology of lava flows from geologic mapping alone. Nonlinear inversion of UAV-collected magnetic data were used to model the thickness and morphology of buried Little Cones’ lava flows with higher precision than achieved previously. The sequence of events and calculated flow characteristics are then interpreted. The total volume of Little Cones, including concealed lava flows, is approximately 0.016 cubic kilometer, and the initial sheet flow erupted in less than 24 hours. The findings presented herein demonstrate that UAV-based magnetic surveys are a reliable method of data collection and an efficient alternative to other survey methods, facilitating development of a three-dimensional perspective of distributed volcanic fields.

Nevada↗