USGS Science⌕ Search

SEARCH · USGS Science

Results for “International Journal of Environmental Science and Technology”

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

Bathymetric and acoustic-backscatter mapping of Lake Sammamish, Washington, during USGS field activity 2021-656-FA

In 2021, scientists from the U.S. Geological Survey (USGS) acquired bathymetric and acoustic-backscatter data within Lake Sammamish, Washington. Mapping was completed to find evidence of past earthquakes such as underwater landslides and is part of a larger USGS project to understand the overall geologic hazards history of the Cascadia Margin region. The survey was conducted using the USGS research vessel Parke Snavely outfitted with an interferometric sidescan-sonar system for swath mapping. Three estimates of bathymetric uncertainty were calculated from the post-processed bathymetric data. This report provides a summary of the mapping mission.

Washington↗
USGS2024

Shallow faulting and folding beneath south‐central Seattle, Washington State, from land‐based high‐resolution seismic‐reflection imaging

The geologic framework of the Seattle fault zone (SFZ) has been extensively studied, but the structure and fault strand locations in the central portion of the fault zone through the city of Seattle have remained controversial. Much of what is known about the SFZ has come from light detection and ranging (lidar)‐topographic surveys and paleoseismic investigations of fault scarps primarily west of Puget Sound, regional gravity and aeromagnetic modeling, and multiscale marine seismic imaging in waters both west and east of Seattle. We analyze ∼24 km of land‐based P ‐wave seismic‐reflection data that fill in a critical gap in our understanding of the SFZ beneath the urban areas of West Seattle, south‐central Seattle, and Mercer Island. These data image deformed strata in the upper 1 km, including upwarped Tertiary rock and younger sediments. Collectively, these data provide evidence for multiple Quaternary‐active thrust faults, back thrusts, and sub‐basins within the SFZ beneath the city of Seattle. The results indicate that multiple and potentially active back thrusts in the upper ∼500 m extend across the length of the SFZ and the entire urban corridor that may be analogous to those on Bainbridge Island west of Puget Sound.

Washington↗
USGS2024

Status of water-level altitudes and long-term and short-term water-level changes in the Chicot and Evangeline (undifferentiated) and Jasper aquifers, greater Houston area, Texas, 2023

Since the early 1900s, groundwater withdrawn from the primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers—has been an important source of water in the greater Houston area, Texas. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of water-level altitudes and water-level changes in these aquifers in the greater Houston area. In this report, the Chicot and Evangeline aquifers are treated as a single aquifer for the purposes of providing annual assessments of regional-scale water-level altitudes and water-level changes over time. In 2023, shaded depictions of water-level altitudes for the Chicot and Evangeline aquifers (undifferentiated) ranged from about 286 feet (ft) below the North American Vertical Datum of 1988 (NAVD 88) to about 169 ft above NAVD 88. The largest decline in water-level altitudes indicated by the 1977–2023 long-term water-level-change map was in south-central Montgomery County southeast of The Woodlands. In comparison, the 1990–2023 long-term water-level-change map depicts the largest declines in water-level altitudes in localized areas at or near certain wells in parts of northwestern Harris County and south-central Montgomery County. The largest rise in water-level altitudes for 1977–2023 is depicted in a relatively large area in southeastern Harris County, whereas the largest rise in water-level altitudes for 1990–2023 is depicted in a relatively large area in central Harris County. The 5-year short-term water-level-change map depicts the largest declines at three wells in northern Fort Bend County, one well in western Harris County, and three wells in south-central Montgomery County and the largest rise at one well in central Harris County. The 1-year short-term water-level-change map depicts the largest declines at one well in northern Fort Bend County and two wells in southwestern Harris County and the largest rises at one well in northern Brazoria County and one well in south-central Montgomery County. In 2023, shaded depictions of water-level altitudes for the Jasper aquifer ranged from about 242 ft below NAVD 88 to about 218 ft above NAVD 88. The 2000–23 long-term water-level-change map depicts water-level declines throughout the study area where water-level-measurement data from the aquifer were collected, with the largest declines in north-central Harris County and south-central Montgomery County south of The Woodlands. The 5-year short-term water-level-change map depicts the largest declines at two wells in central Montgomery County near Conroe and two wells in south-central Montgomery County southeast of The Woodlands and the largest rise at one well in western Montgomery County. The 1-year short-term water-level-change map depicts the largest declines at four wells in south-central Montgomery County southeast of The Woodlands and one well in central Montgomery County near Conroe and the largest rises at two wells in western Montgomery County.

Texas↗
USGS2025

Conceptual and numerical groundwater flow model of the Iowa River alluvial aquifer near Tama County, Iowa, 1980 through 2022

The Iowa River alluvial aquifer is an important source of water on the Meskwaki Settlement in Tama County, Iowa, which is land owned by the Sac & Fox Tribe of the Mississippi in Iowa (commonly known as the Meskwaki Nation). The U.S. Geological Survey constructed a groundwater flow model, including a conceptual and numerical model, of the Iowa River alluvial aquifer and underlying hydrogeologic units near the Meskwaki Settlement in Tama County, Iowa, for the period of January 1980–August 2022 to estimate the fraction of water pumped from the Iowa River alluvial aquifer by Meskwaki Settlement wells that is derived from streamflow depletion in the Iowa River and its tributaries. Streamflow depletion is a reduction in streamflow caused by groundwater pumping and includes the interception by groundwater production wells of water that otherwise would have been discharged to streams (called “captured groundwater discharge”) and induced infiltration of streamflow to the production wells. Calibrated model runs were performed with no simulated pumping and simulated pumping only at Meskwaki Settlement wells, and the change in simulated flow rates between the groundwater system and streams for the two model runs represents the amount of streamflow depletion in the Iowa River and tributary streams resulting from pumping at the Meskwaki Settlement wells. Streamflow depletion in the Iowa River and its tributaries as a percentage of simulated pumping at the Meskwaki Settlement wells was calculated by dividing this difference by the total simulated pumping rate for the Meskwaki Settlement wells. The model results demonstrate that the mean monthly streamflow depletion, including induced infiltration and captured discharge, in the Iowa River and its tributary streams as a percentage of mean monthly pumping at the Meskwaki Settlement wells was 97.4 percent and ranged from 65.4 to 112 percent. Of the total streamflow depletion, mean monthly induced recharge was 20.9 percent and ranged from 4.9 to 37.2 percent. Mean monthly captured discharge was 76.5 percent and ranged from 57.1 to 97.1 percent. These results indicate that most of the water pumped from the Meskwaki Settlement wells is the result of streamflow depletion, in the form of both induced infiltration and captured discharge.

Iowa↗
USGS2025

Dynamic baseflow storage estimates and the role of topography, geology and evapotranspiration on streamflow recession characteristics in the Neversink Reservoir Watershed, New York

Estimates of dynamic groundwater volumes supplying baseflow to streams are important for water availability projections during extended periods of drought. The primary goals of this study were to provide dynamic storage volume estimates, inferred from streamflow recession analysis, for baseflow regimes within seven gaged catchments within the Neversink Reservoir Watershed (NRW), a critical municipal water source for New York City. Additionally, geomorphological properties, surficial geology and hydro-meteorological processes were quantified and described in relation to time and spatially variable recession behaviour and storage estimates across the NRW. To explore these relationships, we (1) evaluated seasonal trends in streamflow recession behaviour in relation to modelled potential evapotranspiration (PET) and catchment runoff rates, (2) derived empirical streamflow models for cool-season runoff using both linear and nonlinear reservoir assumptions for baseflow and (3) calculated metrics related to the geology and geomorphology of each catchment and compared these metrics to area normalised baseflow dynamic storage estimates. Results show that baseflow recession behaves as a nonlinear reservoir, and applying linear groundwater reservoir assumptions may underestimate the total dynamic storage volumes compared to what would be predicted for a nonlinear reservoir. Increases in PET caused decreases in storage conditions that resulted in increased recession rates and nonlinearity in streamflow recession during the growing season. Additionally, we found that while no single physical catchment characteristic solely predicted catchment storage dynamics, sediment volume and stream gradients were stronger predictors of normalised storage volumes than catchment surface area or surface topography alone. Within the NRW, catchments with the highest sediment volume exhibited the lowest recession rates and higher dynamic storage volumes, while the smallest catchment, mostly devoid of sediment, had the fastest recession rate and lowest dynamic storage volume.

New York↗
USGS2000

Assessment of potential for natural attenuation of chlorinated ethenes and ethanes in ground water at a petrochemical reclamation site, Harris County, Texas

Redox conditions in the Numerous Sand Channels Zone beneath a petrochemical reclamation site in Harris County, Texas, range from sulfate reducing to methanogenic as indicated by the presence of methane in ground water and the range of molecular hydrogen concentrations. Assessment of the potential for reductive dechlorination using BIOCHLOR as a screening tool indicated conditions favoring anaerobic degradation of chlorinated organic compounds in the Numerous Sand Channels Zone. Evidence supporting reductive dechlorination includes apparently biogenic cis -1,2-dichloroethene; an increased ratio of 1,2-dichloroethane to 1,1,2-trichloroethane downgradient from the assumed contaminant source area; ethene and methane concentrations greater than background concentrations within the area of the contaminant plume; and a positive correlation of the ratio of ethene to vinyl chloride as a function of methane concentrations. The body of evidence presented in this report argues for hydrogenolysis of trichloroethene to cis -1,2-dichloroethene; of 1,1,2-trichloroethane to 1,2-dichloroethane; and of vinyl chloride to ethene within the Numerous Sand Channels Zone. Simulations using BIOCHLOR yielded apparent first-order decay constants for reductive dechlorination in the sequence tetrachloroethene trichloroethene cis -1,2-dichloroethene vinyl chloride ethene within the range of literature values reported for each compound and apparent first-order decay constants for reductive dechlorination in the sequence 1,1,2-trichloroethane 1,2-dichloroethane slightly greater than literature values reported for each compound along the upgradient segment of a simulated ground-water flowpath. Except for vinyl chloride, apparent rates of reductive dechlorination for all simulated species show a marked decrease along the downgradient segment of the simulated ground-water flowpath. Evidence for reductive dechlorination of chlorinated ethenes within the Numerous Sand Channels Zone indicates potential for natural attenuation of chlorinated ethenes. Reductive dechlorination of chlorinated ethanes apparently occurs to a lesser extent, indicating relatively less potential for natural attenuation of chlorinated ethanes. Additional data are needed on the concentrations and distribution of chlorinated ethenes and ethanes in individual fine sand intervals of the Numerous Sand Channels Zone. This information, combined with lower minimum reporting levels for future chloroethane analyses, might enable a more complete and quantitative assessment of the potential for natural attenuation at the site.

Texas↗
USGS2000

Quality assurance and analysis of water levels in wells on Pahute Mesa and vicinity, Nevada Test Site, Nye County, Nevada

Periodic and continual water-level data from 1963 to 1998 were compiled and quality assured for 65 observation wells on Pahute Mesa and vicinity, Nye County, Nevada. As part of the quality assurance of all water levels, ancillary data pertinent to computing hydraulic heads in wells were compiled and analyzed. Quality-assured water levels that were not necessarily in error but which did not represent static heads in the regional aquifer system, or required some other qualification, were flagged. Water levels flagged include those recovering from recent pumping or well construction, water levels affected by nuclear tests, and measurements affected by borehole deviations. A cursory examination of about 30 wells with available water-level and down-hole temperature data indicate that water levels in most wells on Pahute Mesa would not be significantly affected by temperature if corrected to 95 degrees Fahrenheit. Wells with large corrections (greater than 10 feet) are those with long water columns (greater than 1,500 feet of water above the assumed point of inflow) in combination with mean water-column temperatures exceeding 105 degrees Fahrenheit. Water-level fluctuations in wells on Pahute Mesa are caused by several factors including infiltration of precipitation, barometric pressure, Earth tides, ground-water pumpage, and seismic events caused by tectonic activity and underground nuclear testing. No observed water-level fluctuations were attributed to a naturally occurring earthquake. The magnitude and duration of changes in water levels caused by nuclear tests are affected by the test size and the distance from a well to the test. Identifying water levels that might be affected by past nuclear tests is difficult because pre-testing water-level data are sparse. Hydrologically significant trends were found in 13 of 25 wells with multiple years of water-level record. The largest change in water levels (1,029 feet in 25 years) occurred in well U-19v PS 1D as a result of the Almendro nuclear test. Likely explanations for trends in most of the wells are either changes in precipitation patterns that affect recharge rates to the ground-water system, pumping effects from water-supply well U-20 WW, or a combination of these two factors.

Nevada↗
USGS2025

The δ13C signature of dissolved organic and inorganic carbon reveals complex carbon transformations within a salt marsh

Coastal wetlands have high rates of atmospheric CO 2 uptake, which is subsequently respired back to the atmosphere, stored as organic matter within flooded, anoxic soils, or exported to the coastal ocean. Transformation of fixed carbon occurs through a variety of subsurface aerobic and anaerobic microbial processes, and results in a large inventory of dissolved carbon. Carbon source and the roles of aerobic respiration, sulfate reduction, and methane cycling were evaluated within salt marsh peat and the underlying sandy subterranean estuary. There is a large increase in dissolved inorganic carbon (DIC, 7,350 ± 3,900 μmol L −1 ), dissolved organic carbon (DOC, 1,040 ± 1,480 μmol L −1 ) and CH 4 (14.5 ± 33.3 μmol L −1 ) within the marsh porewaters compared to creek waters. Alkalinity production (5,730 ± 2,170 μeq L −1 ) and sulfate removal (1,810 ± 1,970 μmol L −1 ) indicate anaerobic respiration, however, relative contributions from the various decomposition pathways cannot be identified due to overlapping geochemical signatures. The δ 13 C of the DOC (−29.0 ± 3.7‰) and DIC (−11.2 ± 1.1‰) produced within the marsh differed from the bulk soil organic matter δ 13 C (−14.5 ± 0.2‰). We explore a variety of mechanisms that could result in co-occurring depleted δ 13 C-DOC and enriched δ 13 C-DIC compared to the bulk soil organic carbon pool and salt marsh vegetation, including selective mineralization, production of δ 13 C-depleted bacterial biomass, and methane-derived DOC. While important questions remain about carbon cycling pathways, we found evidence of a cryptic methane cycle. Alteration of the δ 13 C of carbon species complicates source attribution in solid and dissolved phases and careful consideration should be used when carbon is partitioned between in situ salt marsh production and external marine and terrestrial sources.

Massachusetts↗
USGS2026

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter↗
USGS2024

Seismic response comparison of a historical masonry church subject to real and simulated ground motions

In recent years, advanced numerical models and high-performance computing have facilitated the utilization of ground motion time series in the assessment of the non-linear dynamic behavior of historic masonry structures. Since recorded accelerograms can be sparse for specific analysis conditions, stochastic ground motion simulations have become a viable alternative to overcome this limitation. This study simulates the recorded acceleration time series of the Central Italy 2016 earthquake event at the closest station to the town of Macerata using a site-based stochastic approach. The simulated motions are seismologically evaluated using a goodness-of-fit method in terms of various intensity measures. The simulated records, in conjunction with real records, are used to study the non-linear dynamic behavior of San Filippo Neri church located in Macerata. The church of San Filippo represents an important example of Baroque religious architecture in central Italy, which was damaged and closed off to the public after the 2016 earthquake events. The construction was investigated with a vast diagnostic campaign which included on-site testing and dynamic identification tests. The collected data is used to calibrate the dynamic response of a three-dimensional finite element model of the church. The model is finally used to compare the non-linear seismic responses under real and simulated ground motions with the site recorded damage. The results of structural responses demonstrate a strong agreement between the real and simulated records, providing evidence to support the validation of the site-based stochastic simulation.

Macerata↗
USGS2008

Late Jurassic igneous rocks in south-central Arizona and north-central Sonora: Magmatic accompaniment of crustal extension

In Middle Jurassic time, the region of south-central Arizona and north-central Sonora was part of a continental magmatic arc. In early Late Jurassic time, this magmatic arc gradually gave way to a extensional regime. Late Jurassic rifting was accompanied by magmatism, typically small in volume relative to the preceding arc phase, and with mildly alkaline compositions. These igneous rocks, the 158- to 146-Ma Ko Vaya Suite, make up much or most of the Comobabi, Artesa, Quijotoa, and Brownell Mountains and Sierra del Cobre. The bimodal Ko Vaya Suite constitutes volcanic to shallow plutonic complexes comprising chiefly trachyandesite and volcanic wacke; monzodiorite; rhyolite porphyry; distinctive, compositionally and texturally heterogeneous, quartzpoor granite (the Ko Vaya Granite); A-type perthite granite and quartz syenite; aphyric intrusive rhyolite; and locally abundant hematite veins. Alteration, particularly potassic alteration, is widespread. Minimally altered samples of granite and quartz syenite have the composition of trachyte or trachydacite and alkali rhyolite. The Ko Vaya Suite is nonconformably overlain by the latest Jurassic Sand Wells Formation, approximately correlative with the Glance Conglomerate, the basal unit of the Bisbee Group. The Ko Vaya Suite is in some ways analogous to the coeval Independence dike swarm of eastern California, but less diffuse. We infer that the Ko Vaya Suite represents especially intense or focused Late Jurassic extensional magmatism.

Arizona, Sonora↗
USGS2025

Status of water-quality conditions in the United States, 2010–20

Degradation of water quality can make water harmful or unusable for humans and ecosystems. Although many studies have assessed the effect of individual constituents or narrow suites of constituents on freshwater systems, no consistent, comprehensive assessment exists over the wide range of water-quality effects on water availability. Using published studies, data, and models completed at regional or national scales in the United States during 2010–20, this chapter moves towards a comprehensive assessment by summarizing how selected anthropogenic and geogenic water-quality constituents affect national-scale water availability for human and ecosystem needs. Several types of human health, agricultural, ecological, and beneficial-use standards or thresholds were used to provide context for categorizing surface-water and groundwater quality. Water availability for human and ecological use is limited by elevated concentrations of geogenic and anthropogenic constituents in surface and groundwater. Elevated concentrations of five geogenic constituents (arsenic, manganese, strontium, radium, and adjusted gross alpha) are common in groundwater and collectively affect the drinking water supply to over 30 million people. Surface water sourced drinking water supplies are impaired in about a third of assessed stream miles, most commonly because of non-mercury metals and salinity. Health-based violations at community water systems may disproportionately affect socially vulnerable communities. Ecological water uses are predominantly limited by nutrients, sediment, temperature, pathogens, salinity, and pesticides. Water availability for human and ecological use is adversely affected by human activities including human contaminant sources (for example, wastewater, agriculture), processes (for example, dredging, groundwater pumping), or permanent landscape modifications (for example, dams, urbanization). Primary contaminant sources vary spatially and include fertilizer and manure, atmospheric deposition, wastewater treatment plants, urban land, and a range of natural sources. Contaminants of emerging concern, contaminants without regulatory thresholds, and mixtures of geogenic and anthropogenic water contaminants also contribute to ecological degradation and human exposure.

conterminous United States↗
USGS2026

Tectonically driven integration of the 4.8 Ma Colorado River USA tracked with detrital sanidine and fish genetics

The development of the continental-scale Colorado River system, western USA, from 8 to 4.8 Ma, is revealed using 60-40 Ma detrital sanidine tracer grains and fish phylogeny. Here we show that precursor paleoriver segments became integrated north to south as traced by 60-40 Ma sand grains that were derived from the north and sequentially appeared in the 25-8 Ma Browns Park Formation of Utah, 7-6 Ma upper Bidahochi Formation of Arizona, and 4.8 Ma Bouse Formation of the lower Colorado River and proto Gulf of California. This timing is mimicked by molecular clock estimates of divergence times among fish lineages. River integration was a response to headwater uplifts in the Yellowstone hotspot track and Rocky Mountains. 40 Ar/ 39 Ar ages refine the timing for mantle-drips that caused subsidence, then uplift, of depositional basins that influenced the integration pathway and tempo. The ~ 3-million-year timescale suggests that multiscale mantle-driven uplift, rather than lake spillover, was the primary driver for integration of the proto-Colorado River through Grand Canyon.

Colorado River↗
USGS2024

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado↗
USGS2025

Characterization of stream water quality and groundwater levels in the Central Pine Barrens region, Suffolk County, New York, 2017–23

The area locally known as the “Central Pine Barrens” region, located in Suffolk County, New York, contains most of Long Island’s preserved and undeveloped land. This region overlays an aquifer system that provides potable groundwater for residents of Suffolk County. Between 2017 and 2023, the U.S. Geological Survey, in cooperation with the Central Pine Barrens Joint Planning & Policy Commission and the Town of Brookhaven, monitored groundwater levels and stream water quality in this region. Groundwater levels were measured monthly at five wells and continuously (15-minute intervals) at a sixth well. Water quality was monitored at five locations in the Carmans River and at two locations in the Peconic River, and samples were analyzed for major ions, trace elements, nutrients, pharmaceuticals, and pesticides. The major ion compositions at the sites were mainly sodium-chloride type waters, and compositions varied the most at the furthest upstream sites in both streams. Concentrations above aquatic-life criteria thresholds also occurred most frequently at the furthest upstream sites. The seasonal patterns of nutrient loads and concentrations varied between the Carmans and Peconic Rivers. Several organic compounds including pharmaceuticals, domestic use products, and pesticides were detected at low concentrations in both streams. Metformin was the most frequently detected pharmaceutical compound, and herbicides were the most frequently detected pesticide class. Water-quality conditions influenced by anthropogenic contributions are a result of current and historical land use, and these contributions include onsite wastewater disposal systems, commercial or domestic fertilizers and pesticides, and urban or industrial contaminants in road runoff. This study characterizes and improves understanding of the current hydrologic conditions in the Central Pine Barrens region and the study findings can help inform the development of plans to manage, protect, and restore water resources.

New York↗
USGS2025

Critical minerals in orogenic (gold) and Coeur d’Alene-type mineral systems of the United States

Orogenic and Coeur d’Alene-type mineral systems are produced by metamorphic devolatilization of thick volcanic or siliciclastic sedimentary rock sequences and the focused flow of hydrothermal fluids upwards along crustal-scale faults. Most orogenic systems are found along the Cordilleran orogen, stretching from California northwards into Alaska, whereas most Coeur d’Alene-type systems occur in the Proterozoic Belt Basin in Idaho and Montana. Although the deposit types in these systems are exploited for precious and base metals, potential exists for the production of critical minerals necessary for current (2025) societal needs in the United States. Publicly available geochemical data compiled for these mineral systems, coupled with mineralogical characteristics, indicate that several critical minerals could potentially be recovered from unmined resources and processed mine waste: arsenic, antimony, tellurium, cobalt, and tungsten from orogenic gold deposits and zinc, antimony, arsenic, and manganese from Coeur d’Alene-type systems. These critical minerals reside primarily in arsenopyrite (arsenic), scheelite (tungsten), siderite (manganese), sphalerite (zinc), tetrahedrite (antimony and arsenic), stibnite (antimony), and telluride (tellurium) minerals.

continental United States↗
USGS2024

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California↗
USGS2024

Utilization of stochastic ground motion simulations for scenario-based performance assessment of geo-structures

Probabilistic seismic performance assessments of engineered structures can be highly sensitive to the seismic input excitation and its variability. In the present study, the scenario-based performance assessment recommended by Federal Emergency Management Agency (FEMA) P-58 guidelines is adopted to estimate seismic fragility of concrete dams for various seismic hazard scenarios. Due to the scarcity of recorded ground motions and thereby their poor representation of uncertainties, stochastic ground motion simulation methods are utilized to obtain the required input excitations. Moreover, to understand the uncertainty in ground motion simulation models, two broadband stochastic simulation models are used to generate input excitations representing six seismic hazard scenarios defined by earthquake magnitude, source-to-site distance, and soil conditions. Optimal intensity measure parameters for each scenario are identified using a systematic procedure that considers criteria such as efficiency, practicality, proficiency, sufficiency, and hazard compatibility. Fragility curves and surfaces are derived using the cloud analysis technique, taking into account various damage measures and limit state functions. The study finds that the derived fragility curves are particularly sensitive to the selection of earthquake scenarios, the choice of records, and the methods used to calculate fragility curves, with less sensitivity observed to different engineering demand parameters. Given this sensitivity, particularly to ground motion selection, the study highlights the necessity of incorporating both model-to-model variability (epistemic uncertainty) and record-to-record variability (aleatory uncertainty), alongside the established material and modeling uncertainties, in the probabilistic seismic assessment.

Reliability Engineering and System Safety (RESS)↗