USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Hydrology, Regional Studies”

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,662 records · Page 70Linked to original sources
USGS2025

Sensitivity of benthic biota and toxicity of cadmium, cobalt, copper, nickel, lead, and zinc mixtures in Washington, United States, and British Columbia, Canada

Relative sensitivities and responses of juvenile white sturgeon ( Acipenser transmontanus ), Hyalella azteca , two families of mayfly (Ephemerellidae, Heptageniidae), one family of caddisfly (Brachycentridae), and a natural community of benthic macroinvertebrates (BMI) to multiple metals are predicted using previously collected laboratory and field samples and a metal mixture model. Biological responses in single metal exposures are used to parameterize toxicity functions, which include accumulations of hydrogen and selected metals on biological receptors, intrinsic potencies of hydrogen and metals, sensitivities of organisms, and times of exposure. The model then is used to predict responses in multiple metal laboratory exposures and field-collected porewater. The following sensitivity sequence in porewater was determined based on endpoints of survival or total abundance: juvenile white sturgeon greater than (>) Ephemerellidae family > Hyalella azteca > Heptageniidae family about equal to (~) benthic macroinvertebrate community > Brachycentridae family. The fraction of porewater samples that are predicted to have adverse impacts on benthic biota (20-percent or greater negative response) depends on organism sensitivities and metal toxicities, and ranges from 44 to 48 percent for juvenile white sturgeon, 23 to 26 percent for the Ephemerellidae family, 16 to 22 percent for Hyalella azteca , 5 to 8 percent for the Heptageniidae family and BMI community, and 0 percent for the caddisfly family. The most toxic porewater in the upper Columbia River Basin (UCR) is at the backwater bar site at Deadmans Eddy and China Bend. The model also indicates that the element responsible for the most toxic conditions in UCR porewater is copper for all organisms, except Hyalella azteca and the metal-insensitive Brachycentridae family. Copper and lead result in the most toxic conditions for Hyalella azteca . This approach and results can aid in assessing metal toxicity and its potential risk to aquatic biota in ecosystems impacted by historical mining activities.

British Columbia, Washington↗
USGS2025

Geologic map of the Buckner 7.5-Minute quadrangle, Louisa County, Virginia

The Buckner 7.5-minute quadrangle straddles three terrane boundaries in the Piedmont Physiographic Province in central Virginia: the Chopawamsic terrane, the Elk Hill Complex, and the Goochland terrane. In much of the map area, the Elk Hill Complex separates the Chopawamsic and Goochland terranes. Rocks of the Chopawamsic terrane include Ordovician metavolcanic and metasedimentary rocks of the Chopawamsic Formation, Ordovician to Silurian granodiorite sheet intrusions, and Paleozoic mafic intrusions. Silurian to Devonian rocks of the Quantico Formation, mostly garnet-mica schist, crop out in the northwesternmost part of the map area, and are in unconformable contact with rocks of the Chopawamsic Formation on the southeastern limb of the Quantico synclinorium. The main map unit in the Elk Hill Complex is Neoproterozoic mica gneiss, which is in pre-metamorphic fault contact with rocks of the Chopawamsic Formation to the west. The main map unit of the Goochland terrane is the Maidens Gneiss. Except for Jurassic diabase dikes, all rocks on the Buckner 7.5-minute quadrangle were metamorphosed to amphibolite facies during the Alleghanian orogeny and preserve multiple compositional and phyllosilicate penetrative foliations. Evidence of amphibolite-facies metamorphism during the Taconic orogeny is preserved in rocks of the Elk Hill Complex. The entire width of the Maidens Gneiss on the Buckner 7.5-minute quadrangle is within the Spotsylvania high-strain zone and amphibolite-facies mylonitic textures are pervasive. Quartz veins and Jurassic diabase dikes crosscut all older rocks of the quadrangle. Multiple levels of terrace deposits are present along and near the major streams of the quadrangle. The lower terrace deposits are likely remnants of former positions of the Little River on the landscape, whereas higher deposits may be remnants of former deposits of the Atlantic Coastal Plain that covered this portion of the Piedmont Province. A linear cluster of aftershocks from the magnitude 5.8 earthquake that occurred near Mineral, Virginia, in 2011 defines the Fredericks Hall fault, which is at depth on the Buckner quadrangle. Most of the aftershocks occurred in the core of the Elk Creek antiform and have no relation to faults mapped at the surface. Several abandoned crushed stone and building stone quarries, as well as a mica prospect, exist in the quadrangle.

Virginia↗
USGS2026

Importance of fish in the diet of Adélie penguins across multiple life stages

Over recent decades, the Adélie penguin ( Pygoscelis adeliae ) population has grown across most of its range, the exception being the northern coast of the western Antarctic Peninsula (WAP). In the Ross Sea, the very large Cape Crozier colony grew to reach ~9% of the global population. Demographic factors driving the Cape Crozier trend have yet to be identified. However, low chick fledging mass - a negative influence in the northern WAP where the diet is mostly less-energy-dense krill - is not showing an effect in the Cape Crozier colony. Previously, we hypothesized that Cape Crozier fledglings, and post-breeding adults, must be finding sufficient higher-quality prey once free of the prey-depleted colony foraging area, with subsequent higher survival. Here we explore penguin diet within and outside the colony foraging area using stable isotope analysis (SIA) of feathers: those grown by post-breeding adults outside the foraging area, prior to moult, compared to those of near-fledged chicks raised on meals obtained within the foraging area. Second, we explore whether a major difference occurs between breeding and post-breeding diet of adults from a small, nearby colony (Cape Royds) that exploits a smaller, little-depleted preyscape. Finally, we compare the results to SIA analyses published by others to explore whether the pre-moult diet of Ross Island adults, mostly foraging within the north-eastern Ross Sea, compares with that of adults nesting in the southern WAP (Bellingshausen Sea coast). At the latter, populations are not decreasing, with diet containing an appreciable contribution of fish. The results confirmed that fish contributed greatly to chick diets at both Ross Island colonies (δ 15 N 10.6–11.4). Once adults were released from central-place foraging, SIA levels (δ 15 N 9.7–12.9) were similar to those of the southern WAP. Foraging on energy-dense fish in coastal Antarctic waters could be one reason Adélie penguins are among the most abundant and therefore most ecologically successful penguin species.

Antarctic Science↗
USGS1993

Hydraulic properties of the Madison aquifer system in the western Rapid City area, South Dakota

Available information on hydrogeology, data from borehole geophysical logs, and aquifer tests were used to determine the hydraulic properties of the Madison aquifer. From aquifer-test analysis, transmissivity and storage coefficient were determined for the Minnelusa and Madison aquifers, and vertical hydraulic conductivity (K v ') along with specific storage (Ss') for the Minnelusa confining bed. Borehole geophysical well logs were used to determine the thickness and location of the Minnelusa aquifer, the lower Minnelusa confining bed, and the Madison aquifer within the Madison Limestone. Porosity values determined from quantitative analysis of borehole geophysical well logs were used in analyzing the aquifer-test data. The average porosity at the two aquifer-test sites is about 10 percent in the Minnelusa aquifer, 5 percent in the lower Minnelusa confining bed, and 35 percent in the Madison aquifer. The first aquifer test, which was conducted at Rapid City production well #6, produced measured drawdown in the Minnelusa and Madison aquifers. Neuman and Witherspoon's method of determining the hydraulic properties of leaky two-aquifer systems was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. Analysis of the aquifer test for the Minnelusa aquifer yielded a transmissivity value of 12,000 feet squared per day and a storage coefficient of 3 x 10 -3 . The specific storage of the Minnelusa confining bed was 2 x 10 -7 per foot, and its vertical hydraulic conductivity was 0.3 foot per day. The transmissivity of the Madison aquifer at this site was 17,000 feet squared per day, and the storage coefficient was 2 x 10 -3 . The second aquifer test, which was conducted at Rapid City production well #5 (RC-5) produced measured drawdown only in the Madison aquifer. Hantush and Jacob's method of determining the hydraulic properties of leaky confined aquifers with no storage in the confining bed was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. The analysis of data from the RC-5 aquifer test showed that transmissivity was not equal in all directions. Hantush's method was used to determine the direction of radial anisotropy and magnitude of the major and minor axes of transmissivity. The major axis of transmissivity is at an angle of 42° east of north, and the transmissivity along this axis is about 56,000 feet squared per day. The minor axis of transmissivity is at an angle of 48° west of north, and the transmissivity along this axis is about 1,300 feet squared per day. The major axis of transmissivity intersects Cleghorn Springs, a large resurgent spring on the west edge of Rapid City. The shape of the potentiometric contours of the Madison aquifer near RC-5 agree with the orientation of the transmissivity ellipse. The average value of the storage coefficient from the isotropic analysis of the aquifer-test data was 3.5 x 10 -4 , and the average vertical hydraulic conductivity of the lower Minnelusa confining bed was 9.6 x 10 -3 foot per day.

South Dakota↗
USGS1993

Geochemistry of and radioactivity in ground water of the Highland Rim and Central Basin aquifer systems, Hickman and Maury counties, Tennessee

A reconnaissance of the geochemistry of and radioactivity in ground water from the Highland Rim and Central Basin aquifer systems in Hickman and Maury Counties, Tennessee, was conducted in 1989. Water in both aquifer systems typically is of the calcium or calcium magnesium bicarbonate type, but concentrations of calcium, magnesium, sodium, potassium, chloride, and sulfate are greater in water of the Central Basin system; differences in the concentrations are statistically significant. Dissolution of calcite, magnesium-calcite, dolomite, and gypsum are the primary geochemical processes controlling ground-water chemistry in both aquifer systems. Saturation-state calculations using the computer code WATEQF indicated that ground water from the Central Basin system is more saturated with respect to calcite, dolomite, and gypsum than water from the Highland Rim system. Geochemical environments within each aquifer system are somewhat different with respect to dissolution of magnesium-bearing minerals. Water samples from the Highland Rim system had a fairly constant calcium to magnesium molar ratio, implying congruent dissolution of magnesium-bearing minerals, whereas water samples from the Central Basin system had highly variable ratios, implying either incongruent dissolution or heterogeneity in soluble constituents of the aquifer matrix. Concentrations of radionuclides in water were low and not greatly different between aquifer systems. Median gross alpha activities were 0.54 picocuries per liter in water from each system; median gross beta activities were 1.1 and 2.3 picocuries per liter in water from the Highland Rim and Central Basin systems, respectively. Radon-222 concentrations were 559 and 422 picocuries per liter, respectively. Concentrations of gross alpha and radium in all samples were substantially less than Tennessee’s maximum permissible levels for community water-supply systems. The data indicated no relations between concentrations of dissolved radionuclides (uranium, radium-226, radium-228, radon-222, gross alpha, and gross beta) and any key indicators of water chemistry, except in water from the Highland Rim system, in which radon-222 was moderately related to pH and weakly related to dissolved magnesium. The only relation among radiochemical constituents indicated by the data was between radium-226 and gross alpha activity; this relation was indicated for water from both aquifer systems.

Tennessee↗
USGS2025

An exploration of the relative influence of physical models for Omori’s law

Omori’s law states that the rate of aftershocks decays as a function of inverse time. There are multiple physical explanations that we reduce into a nonlinear mixed effects relation of three terms: (1) a Rate/State expression that can account for static/dynamic and viscoelastic triggering caused directly by the mainshock, (2) a fluid diffusion triggering term, and (3) a randomized secondary triggering (cascade) term. We fit free physical-model parameters to an observed aftershock sequence through two nonlinear regression methods to find the relative contributions of physics-based models in an observed aftershock sequence. Results from both methods show that Rate/State models overpredict aftershock rates by ∼0–30%. Secondary aftershocks cause a net negative contribution (seismicity rate reduction that corrects overprediction by other terms) ranging between ∼0 and 30%. All regression solutions yield negative secondary triggering contributions without being guided to do so. A physical explanation for this is that aftershock occurrence relieves stress from the crust, ultimately causing the sequence to extinguish itself. Fluid diffusion triggering contributions range from ∼0 to 20%. Diffusion processes are observed to be shorter in time than the full duration of an aftershock sequence and they are also spatially limited, diminishing their influence. Our results apply to an aftershock decay curve from the 2016 Central Apennines earthquake sequence, meaning that our specific results may not be general. Our primary conclusion is that any one physical model cannot alone fit the observed sequence as well as the combination of three we investigated.

Frontiers in Earth Science↗
USGS2026

Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults

Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.

Callifornia↗
USGS2026

Methodology for quantitative spatial sensitivity analysis of volcanic geodetic networks

Introduction This report introduces a methodology for assessing the state of the U.S. Geological Survey Volcano Observatories’ geodetic monitoring networks that measure how volcanoes deform or change shape. This new method uses a model-based approach that considers the uniqueness of the instrument environments at each volcano. This report focuses on simplified volcanic sources, is independent of the shape or size of the volcano, or the network geometry, and thus highlights the strengths and potential vulnerabilities of each volcano’s geodetic network in an actionable visual format. This analysis can help observatories to make informed decisions about whether volcanoes have an adequate level of geodetic monitoring and indicate where improvements are needed.

Lassen Peak, Mount Shasta↗
USGS2026

Geologic map of the Sierra Nevada, California and western Nevada

THE GEOLOGIC MAP OF THE SIERRA NEVADA is a core component of the Sierra Nevada Earth Science Atlas, which also includes geophysical, neotectonic, economic, and geochronologic data. The map illustrates the distribution of geologic units across the Sierra Nevada and related adjacent areas. Geologic units are grouped by type and age into three categories: Principally Paleozoic and Mesozoic metasedimentary and metavolcanic wall rocks, most of which are grouped into terranes; Paleozoic and Mesozoic plutons and intrusive suites, which intrude the wall rocks and form the core batholith of the range; and Late Cretaceous and Cenozoic sedimentary and volcanic rocks and surficial deposits that unconformably overlap the older units. Related rock units were combined and simplified for presentation at a scale of 1:400,000, as shown in the list of map units and associated correlation of map units (Plate 1B) and description of map units (Appendix A). Tectonic faulting, folding and uplift have overprinted the rocks and have strongly influenced the spatial distribution of units and the distinct morphology of the Sierra Nevada as we see it today. The Atlas is the result of collaborative work by scientists and mapmakers from the California Geological Survey and the U.S. Geological Survey. The Atlas was originally envisioned by the late geologist Warren Nokleberg (1939-2021), who contributed much to the initial geologic map compilation

California, Nevada↗
USGS2024

Anaerobic biodegradation of perfluorooctane sulfonate (PFOS) and microbial community composition in soil amended with a dechlorinating culture and chlorinated solvents

Perfluorooctane sulfonate (PFOS), one of the most frequently detected per - and polyfluoroalkyl substances (PFAS) occurring in soil, surface water, and groundwater near sites contaminated with aqueous film-forming foam (AFFF), has proven to be recalcitrant to many destructive remedies, including chemical oxidation. We investigated the potential to utilize microbially mediated reduction (bioreduction) to degrade PFOS and other PFAS through addition of a known dehalogenating culture, WBC-2, to soil obtained from an AFFF-contaminated site. A substantial decrease in total mass of PFOS (soil and water) was observed in microcosms amended with WBC-2 and chlorinated volatile organic compound (cVOC) co-contaminants — 46.4 ± 11.0 % removal of PFOS over the 45-day experiment. In contrast, perfluorooctanoate (PFOA) and 6:2 fluorotelomer sulfonate (6:2 FTS) concentrations did not decrease in the same microcosms. The low or non-detectable concentrations of potential metabolites in full PFAS analyses, including after application of the total oxidizable precursor assay, indicated that defluorination occurred to non-fluorinated compounds or ultrashort-chain PFAS. Nevertheless, additional research on the metabolites and degradation pathways is needed. Population abundances of known dehalorespirers did not change with PFOS removal during the experiment, making their association with PFOS removal unclear. An increased abundance of sulfate reducers in the genus Desulfosporosinus (Firmicutes) and Sulfurospirillum (Campilobacterota) was observed with PFOS removal, most likely linked to initiation of biodegradation by desulfonation. These results have important implications for development of in situ bioremediation methods for PFAS and advancing knowledge of natural attenuation processes.

Science of the Total Environment↗
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↗
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↗
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↗