USGS Science⌕ Search

SEARCH · USGS Science

Results for “Science of the Total Envionrment”

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,687 records · Page 20Linked to original sources
USGS2026

Assessing environmental drivers of denitrification in restored riverine floodplains

Restoration of impaired floodplains is an increasingly prevalent strategy for alleviating water quality concerns and reducing downstream flooding at watershed scales. Floodplains temporarily store water and slow flow velocity to promote sedimentation during overbank flooding and remove inorganic nitrogen from floodwater and groundwater via denitrification. Evaluating the impacts of different restoration strategies on denitrification can inform more strategic investments into floodplain modifications that improve water quality outcomes. Our research investigates how denitrification rates in floodplains respond to environmental factors that are actionable from an engineering perspective through design and water resources management. We seasonally measured soil denitrification enzyme activity and various environmental characteristics in 4 floodplains with different restoration design and management approaches at the confluence of the Wabash and Tippecanoe Rivers in Indiana, United States. Our results showed that denitrification rates in an agricultural floodplain were significantly lower than in restored floodplains with native vegetation. Certain soil conditions characteristic of floodplain wetlands were associated with higher denitrification, particularly elevated total nitrogen, moisture, silt, and organic matter contents. Vegetation species composition was correlated with denitrification rates. This link may reflect the direct effects of vegetation on soil conditions, such as supplying labile organic carbon, or indirect effects, such as vegetation acting as an indicator of hydrologic regime and land use. Denitrification seasonally varied, peaking in winter when nitrate supply from rivers draining agricultural watersheds in the region is also high. Substrate limitation of soil denitrification enzyme activity was most significant during the summer when overbank flooding, which replenishes soil nitrogen stocks, rarely occurs. Our findings indicate that denitrification capacity will likely be maximized in riverine floodplains that are restored as wetlands with diverse native vegetation and enhanced hydrologic connectivity. Such restoration activities promote higher denitrification rates via elevated moisture, fine sediment deposition, and soil organic matter.

Tippecanoe River, Wabash River↗
USGS2024

Post-glacial stratigraphy and late Holocene record of great Cascadia earthquakes in Ozette Lake, Washington, USA

Ozette Lake is an ~100-m-deep coastal lake located along the outer coast of the Olympic Peninsula (Washington, USA); it is situated above the locked portion of the northern Cascadia megathrust but also relatively isolated from active crustal faults and intraslab earthquakes. Here we present a suite of geophysical and geological evidence for earthquake-triggered mass transport deposits (MTDs) and related turbidite deposition in Ozette Lake since ca. 14 ka. Comprehensive high-resolution bathymetry data, seismic reflection profiles, and sediment cores are used to characterize the post-glacial stratigraphic framework and examine paleoseismic evidence in the lacustrine sediments. Stacked sequences of MTDs along the steep eastern flanks of the lake appear to grade basin-ward from thick, chaotic, blocky masses to thin, parallel-bedded turbidite beds. The discrete turbidite event layers are separated by fine-grained (silt and clay) lake sedimentation. The event layers are observed throughout the lake, but the physical characteristics of the deposits vary considerably depending on proximity to primary depocenters, steep slopes, and subaqueous deltas. A total of 30–34 event deposits are observed in the post-glacial record. Radiometric dating was used to reconstruct a detailed sedimentation history over the last ~5.5 k.y., develop an age model, and estimate the recurrence (365–405 yr) for the most recent 12 event layers. Based on sedimentological characteristics, temporal overlap with other regional paleoseismic chronologies, and recurrence estimates, at least 10 of the dated event layers appear to be sourced from slope failures triggered by intense shaking during megathrust ruptures; the recurrence interval for these 10 events is 440–560 yr. Thus, Ozette Lake contains one of the longest and most robust geological records of repeated shaking along the northern Cascadia subduction zone.

Washington↗
USGS2024

Slip rate for the Rose Canyon fault through San Diego, California, based on analysis of GPS data: Evidence for a potential Rose Canyon–San Miguel-Vallecitos fault connection?

The Rose Canyon fault is the southern extension of the larger Newport–Inglewood–Rose Canyon fault system, which represents a major structural boundary in the Inner Continental Borderland (ICB) offshore of southern California. Ten to fifteen percent of total plate boundary motion in southern California is thought to be accommodated by the faults of the ICB, but the exact distribution of slip is uncertain. With an onshore segment, the Rose Canyon fault offers an opportunity to measure the slip rate using traditional geodetic methods. In this study, we use Global Positioning System (GPS) surface velocities from a combined campaign and continuous GPS network to constrain elastic models of the Rose Canyon fault. We then compare the observed surface velocities with proposed conceptual models of regional fault connections that facilitate the transfer of slip into the Rose Canyon fault to assess how well the observations are explained by the models. The results of elastic half‐space models suggest that the Rose Canyon fault may be slipping toward the higher end of geologic estimates, with the preferred model indicating a slip rate of 2.4 ± 0.5 mm/yr. Although limited in terms of near‐fault benchmarks, we find an improved model fit using an asymmetrical elastic half‐space model and a higher slip rate, suggesting a potential rheological contrast across the Rose Canyon fault, similar to observations from the northern Newport–Inglewood fault segments. Observed GPS surface velocities, background seismicity, and gravity anomalies south of San Diego Bay point toward a more easterly trace for the Rose Canyon fault, suggesting a possible connection with the San Miguel–Vallecitos fault system. Such a connection could increase the potential rupture lengths of future earthquakes and have important consequences for regional seismic hazards.

California↗
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

Characterizing Meteor Crater impact melts through geochemistry and textural analysis

The U.S. Geological Survey Astrogeology Science Center houses the Meteor Crater sample collection, an assemblage of over 2,500 meters of cuttings from 161 drill holes into Meteor Crater’s rim, flanks, and ejecta blanket. We have utilized this unique collection to study the composition and spatial distribution of impact-generated materials from within the ejecta blanket. Meteor Crater has historically been known to have generated only a relatively small amount of impact melt compared to other terrestrial craters of similar size. A detailed compositional and textural dataset of impact-derived melts from this impact can therefore be a useful asset in improving our understanding of crater formation, and in particular impact melt formation. We have characterized 42 impact-melt particles from Meteor Crater using a scanning electron microscope and an electron microprobe for textural and compositional analysis. We analyzed samples from six drill holes in the ejecta blanket, situated to the northwest, southeast, south, and southwest of the crater (ejecta northeast of the crater is devoid of impact melts). Impact melts were collected from drill cuttings at various depths within the ejecta blanket, ranging from a few centimeters below the surface down to ~6.5 meters. Backscattered electron (BSE) images were acquired for each analyzed impact-melt particle. To characterize the various textures and phases present in each impact melt, we also took many detailed BSE images. Our geochemical analyses include full spectral profiles using energy dispersive X-ray spectrometry and well-calibrated wavelength dispersive spectrometry for a number of phases, including minerals (olivine, pyroxene, and so on), pristine glass, and metallic inclusions. The full dataset is available in ScienceBase as a data release (Gullikson and others, 2024), accessible at https://doi.org/10.5066/P9OGAJ8P . Our goal for this Open-File Report is to provide a summary of this immense dataset, details on data collection, descriptions of the different phases observed within impact-melt particles (both geochemically and texturally), and observable trends.

Arizona↗
USGS2026

Scenarios to assess the future water availability in the Mississippi River Valley Alluvial Aquifer for the Cache River and Grand Prairie Regions of Arkansas

The U.S. Geological Survey, as part of the Arkansas Groundwater Initiative, developed forecast scenarios using previously calibrated MODFLOW 6 groundwater models that focused on the Cache and Grand Prairie Critical Groundwater Areas to assess the impact of future climate and water management strategies on the Mississippi River Valley alluvial aquifer. A Soil Water Balance model was used to forecast recharge and irrigation water use. The forecast scenario period was from January 1, 2019, through December 31, 2055, with monthly stress periods. Twenty scenarios were simulated and included seven alternate climate forecasts, five 13 general groundwater pumping reduction scenarios (round 1), and groundwater pumping reduction scenarios by crop type and for the Bayou Meto Water Management Project and Grand Prairie Area Demonstration Project (round 2). Declines in saturated thickness within the Cache Critical Groundwater Area were larger for 18 of the 20 scenarios as compared to outside of the Critical Groundwater Area. The largest average increase in saturated thickness inside the Critical Groundwater Area was 6.4 m which occurred for the round 1, 50 percent reduction scenario. Automatic reductions in groundwater pumping by MODFLOW 6 in the Cache simulation ranged from 0.02 to 13.1 percent of total groundwater pumping. For the Grand Prairie model domain, the average change in saturated thickness of the Mississippi River Valley alluvial aquifer inside the Critical Groundwater Area for the forecast period ranged between -6.6 to 1.7 m. The average saturated thickness of the Mississippi River Valley alluvial aquifer inside the Grand Prairie Critical Groundwater Area declined for 16 of the 20 scenarios. The average reduction in requested groundwater pumping for all scenarios inside the Grand Prairie Critical Groundwater Area was 25.1 percent, and the largest reduction was 46.5 percent.

ESS Open Archive↗
USGS2026

Effects of human activity on Khumbu Glacier: Towards a sustainable Everest Base Camp

The Everest Base Camp (EBC), a critical staging area for mountaineers attempting to summit Mount Everest, has seen a significant increase in human activity over the past few decades. The increasing number of climbers, trekkers, and support staff, particularly at the campsite area, has intensified environmental impacts on the Khumbu Glacier. This study investigates the effects of human activity at EBC, focusing on three contributors to glacier melt: global climate change, local fossil fuel consumption, and human urinary discharge. Through field data collected during the spring 2023 climbing season, this study examines the energy released through the use of liquefied petroleum gas (LPG), kerosene, and petrol for cooking, heating, and electricity, as well as the heat generated by human urinary discharge. The total heat energy released from these activities at EBC in the spring of 2023 was 849,174 ± 179,774 MJ, which would be sufficient to melt 2492 ± 528 tons of glacier ice and snow. Trends in land surface temperature at EBC and adjacent environments during the 1991–2023 period were calculated from Landsat 5, 7, and 9 satellite data. The 32-year Landsat record reveals that EBC surface temperature increased by 0.28 °C year −1 , which was roughly twice the rate of warming at the surface of Khumbu Glacier adjacent to EBC and about 15% higher than the temperature increase of the debris-covered terrain immediately to the north of EBC. The findings suggest that addressing these anthropogenic influences could help to preserve the Khumbu Glacier and support the sustainability of mountaineering in the region. This paper also takes a transformative approach and explores stable and safer locations in case there is a necessity to relocate the current EBC. Two sites have been identified southwest of the existing EBC, which, unlike the current supraglacial site, are situated on stable ground.

Everest Base Camp, Khumbu Glacier↗
USGS2025

Flood-inundation maps for 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, 2023

Digital flood-inundation map libraries for two reaches that constitute 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, were created by the U.S. Geological Survey (USGS) in cooperation with the Papio-Missouri River Natural Resource District. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://www.usgs.gov/mission-areas/water-resources/science/flood-inundation-mapping-fim-program , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at Little Papillion Creek at Irvington, Nebr. (USGS station 06610750), Little Papillion Creek at Ak-Sar-Ben at Omaha, Nebr. (USGS station 06610765), and Big Papillion Creek at Q Street at Omaha, Nebr. (USGS station 06610770) streamgages. Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System database at https://doi.org/10.5066/F7P55KJN or from the National Weather Service Advanced Hydrologic Prediction Service at https://water.weather.gov/ahps/ . Flood profiles were computed for two different reaches that constitute 14.8 miles of stream length in the study area by using hydraulic models. The models were calibrated by adjusting roughness coefficients to best represent the current (2022) stage-streamflow relation at the streamgages within the study reach. The hydraulic models were then used to compute water-surface profiles at 1-foot stage intervals for selected stage ranges to represent various flooding scenarios at the streamgages in each reach. The simulated water-surface profiles then were combined with a digital elevation model using a geographic information system, which had a 10-foot grid spacing to delineate the flooding extents and water depths for each stage. The availability of these flood-inundation maps, along with information regarding current stage from the USGS streamgages, can provide emergency management personnel and residents with information that is critical for flood response activities and post flood recovery efforts.

Nebraska↗
USGS2026

A ground-motion model derived using a generalized mean rupture distance for large slab interface earthquakes

Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup ⁠ . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.

Seismological Research Letters↗
USGS2025

An evaluation of the effects of different deicing salt application rates on three watersheds in Essex County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Transportation, evaluated the effects of different deicing salt application rates on surface water, groundwater, and highway runoff quality near State highways in northern New York. Three reaches of State highways were tested with different deicing treatments between October 2019 and November 2022: a salt-sand mixture (Treatment A), a salt mixture applied at a lower rate (Treatment B), and a control mixture consistent with typical deicing salt amounts and application rates. Data on pavement conditions and the quality of surface water, highway runoff, and groundwater were collected. Surface electromagnetic data were also collected. Surface-water and groundwater quality downgradient from the State highways were compared with water quality at upgradient locations. The percentage of snow or ice coverage was used to evaluate the effectiveness of the salt applications. This report provides an overview of the transport of deicing salt. The Treatment B watershed had deicing mixture applied more frequently than other highway reaches, which caused it to have the highest annual total chloride application. Despite differences in chloride application, flow-weighted mean chloride concentrations in highway runoff were comparable across treatments. Chloride concentrations were elevated in surface water and groundwater downgradient from highways relative to chloride concentrations upgradient from highways. A chloride mass balance, calculated for one treatment watershed, indicated that groundwater affected by legacy deicing practices may be contributing additional chloride to surface water. Spatial patterns from electromagnetic surveys show a shallow saline plume alongside the highway in that area. Differences in winter severity and pavement-surface conditions drove deicing salt applications in the treatment areas. This study found that several factors affect chloride loads in the watersheds, including variable winter conditions, adaptive snow and ice management, legacy management practices, and area-specific aquifer and groundwater conditions.

New York↗
USGS2024

Mapping karst groundwater flow paths and delineating recharge areas for springs in the Little Sequatchie and Pryor Cove watersheds, Tennessee

The Little Sequatchie River and Pryor Cove Branch, in southern Tennessee, drain the eastern escarpment of the Cumberland Plateau to the Sequatchie River near the southern end of the Sequatchie Valley. The Little Sequatchie River is the largest tributary to the Sequatchie River by drainage area, covering over 120 square miles. The hydrology of the two drainage areas has been largely altered by karst processes, which has caused the majority of the streams to sink at the contact between the Mississippian Pennington Formation and the underlying Mississippian Bangor Limestone. The U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service and Tennessee Department of Environment and Conservation, initiated a study in 2021 to map the karst groundwater pathways in both watersheds in order to delineate recharge areas for several springs. One of these springs, Sequatchie Cave, represents a significant habitat for two Species of Greatest Conservation Need, the Glyphopsyche sequatchie (Sequatchie caddisfly) and the federally endangered Marstonia ogmorhaphe (royal marstonia). Springs and springflow-dominated streams in the Little Sequatchie River valley and Pryor Cove also provide water for agricultural practices and serve as a drinking water source for nearby communities. During the study, a total of 25 dye injections were conducted over eight rounds from January 2022 through March 2023. Dye traces from these injections helped to delineate recharge areas for six major springs, ranging from 7.3 to 65.2 square miles in area. The majority of the dye traces remained subsurface (from sinkpoint to recovery site) for long distances, with karst groundwater travelling nearly 8 miles before resurfacing. The dye traces also had rapid traveltimes, often travelling hundreds to thousands of feet per hour. The goal of this project was to provide scientific data related to karst groundwater pathways and spring recharge areas to aid State and Federal agencies in making informed decisions to protect and preserve this unique and vulnerable karst system.

Tennessee↗
USGS2024

On-orbit calibration and performance of the EMIT imaging spectrometer

The Earth surface Mineral dust source InvesTigation (EMIT) is a remote visible to shortwave infrared (VSWIR) imaging spectrometer that has been operating onboard the International Space Station since July 2022. This article describes EMIT's on-orbit spectroradiometric calibration and validation. Accurate spectroscopy is vital to achieve consistent mapping results with orbital imaging spectrometers. EMIT takes a unique approach to this challenge, with just six optical elements, no shutter, and no onboard calibration systems. Its simple design focuses on uniformity and stability to enable vicarious spectroradiometric calibration. Our experiments demonstrate that this approach is successful, approaching the fidelity of manual field spectroscopy in some cases, and enabling new and more accurate products across diverse Earth science disciplines. EMIT achieves several notable firsts for an instrument of its class. It demonstrates successful on-orbit adjustments of Focal Plane Array (FPA) alignment with sub-micron precision. It offers spectral uniformity better than 98%. Optical artifacts in the measurement channels are at least three orders of magnitude below the primary solar-reflected surface signals. Its noise performance enables percent-level discrimination in the depths of mineral absorption features. In these aspects, EMIT satisfies the stringent performance needs for the next generation of VSWIR imaging spectrometers to observe the Earth's ecosystems, geology , and water resources.

Remote Sensing of Environment↗
USGS2025

Sources of water and salts for the Zuni Salt Lake in west-central New Mexico

The Zuni Salt Lake is located in a maar in west-central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. There have been several investigations suggesting different sources for the water and salt to the lake. Springs, seeps, and ephemeral streamflow have all been observed to contribute freshwater to the lake, and brackish to hypersaline seeps have been documented along the banks of the lake. This report summarizes the findings of a study that characterizes the lake’s hydrology, its water and salinity sources, and the hydrogeologic conceptual model. Regional groundwater levels indicate that each of the aquifers in the area have the potential to discharge groundwater to the lake. There is also evidence of vertical groundwater flow pathways at the maar that were likely created by the igneous intrusion that fractured the intersecting aquifers. A detailed water budget was constructed from continuous lake stage, precipitation, and evaporation data to estimate the groundwater inflow to the Zuni Salt Lake. It was determined that groundwater inflow to the lake is 441 ±94 acre-feet per year, which composes as much as 77 percent of the total inflows. The high sodium and chloride concentrations measured in two hypersaline samples collected near the lake indicate that the majority of the dissolved solids entering the lake are from a hypersaline groundwater source. The geochemical and isotopic compositions measured in the lake and surrounding features support the interpretation that hypersaline groundwater is the primary source of salts to the lake, which is likely sourced from the older (and deeper) Permian units. The hypersaline groundwater samples collected during this investigation have a unique aqueous chemistry relative to each of the mapped aquifers, and variability in groundwater compositions is interpreted to result from differences in minerology and residence time.

New Mexico↗
USGS2026

Top Elevation of Glacial Till and Thickness of the Big Sioux Aquifer Delineated From Electrical Resistivity Tomography Surveys Near Sioux Falls, South Dakota, 2022 and 2025

The City of Sioux Falls, South Dakota, requested the U.S. Geological Survey perform electrical resistivity surveys on three parcels of land north of the city. Electrical resistivity data were collected along a total of 22 transects during March 14–18, 2022, and November 17–21, 2025. Results from electrical resistivity surveys were used to delineate the top of glacial till deposits for the purpose of characterizing the Big Sioux aquifer near the city. Delineating geologic contacts provides important information on groundwater storage, flow dynamics, well design and placement, contaminant transport, groundwater–surface-water interactions, and regional water modeling. The top elevation of glacial till and the thickness of the Big Sioux aquifer varied among the three survey areas. The interpreted top elevation of glacial till in the North survey area decreases from east to west toward a slough, with elevations ranging from 1,403 to 1,418 feet (ft). The estimated thickness of the Big Sioux aquifer in the North survey area increased from east to west, with thicknesses ranging from 23 to 38 ft. The top elevation of glacial till in the Well 72 survey area generally decreases from northwest to southeast. Top elevations of the glacial till in the Well 72 survey area ranged from 1,400 to 1,409 ft along the southern end of transect W72_2. The estimated thickness of the Big Sioux aquifer in the Well 72 survey area was greatest along a southeast to northwest trending channel, with thicknesses ranging from 28 to 40 ft. The top elevation of glacial till in the Nose survey area generally decreases west toward the Big Sioux River. Top elevations of the glacial till in the Nose survey area ranged from 1,362 to 1,395 ft. The estimated thickness of the Big Sioux aquifer in the Nose survey area ranged from 33 to 70 ft.

South Dakota↗
USGS2024

Environmental monitoring of groundwater, surface water, and soil at the Ammonium Perchlorate Rocket Motor Destruction Facility at the Letterkenny Army Depot, Chambersburg, Pennsylvania, 2021

Letterkenny Army Depot in Chambersburg, Pennsylvania, built an Ammonium Perchlorate Rocket Motor Destruction (ARMD) Facility in 2016 to centralize rocket motor destruction and contain all waste during the destruction process. The U.S. Geological Survey has collected environmental samples from groundwater, surface water, and soils at ARMD since 2016. During 2021, samples were collected from four groundwater wells in September, one surface-water site in October, and five soil sites in November near the facility. Samples were analyzed for nutrients, trace metals, major ions, total volatile organic compounds, and perchlorate. Perchlorate was not detected in any 2021 samples. Groundwater results showed no constituents exceeded any U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL). Dissolved arsenic (As) was detected in one well above the reporting detection level (RDL) of 3 micrograms per liter (μg/L) at 5.4 μg/L but below its MCL of 10 μg/L. Dissolved iron (Fe) was the only inorganic constituent measured above an EPA secondary maximum contaminant level (SMCL). All groundwater samples collected in 2021 exceeded the Fe SMCL of 300 μg/L, with concentrations ranging from 390 μg/L to 3,500 μg/L. Surface-water data collected during 2021 showed no measured constituents in the surface-water sample that exceeded any EPA MCL or SMCL. Soil samples collected from 2016 through 2021 showed all concentrations of As exceeded the EPA soil screening levels of 3 milligrams per kilogram (mg/kg) but did not exceed the Pennsylvania medium-specific concentrations for As of 61 mg/kg. Arsenic concentrations in 2021 ranged from 9.1 mg/kg to 12.9 mg/kg. The 2021 results for the ARMD Facility indicate no increases in concentrations of reported compounds compared to data from 2016 to 2020. The contained burn treatment facility for demilitarization of rocket motors during 2021 appears to have operated without elevating concentrations of target compounds compared to previous years.

Pennsylvania↗
USGS2023

Don’t Let Negatives Hold You Back: Accounting for Underlying Physics and Natural Distributions of Hydrothermal Systems When Selecting Negative Training Sites Leads to Better Machine Learning Predictions

Selecting negative training sites is an important challenge to resolve when utilizing machine learning (ML) for predicting hydrothermal resource favorability because ideal models would discriminate between hydrothermal systems (positives) and all types of locations without hydrothermal systems (negatives). The Nevada Machine Learning project (NVML) fit an artificial neural network to identify areas favorable for hydrothermal systems by selecting 62 negative sites where the research team had confidence that no hydrothermal resource exists. Herein, we compare the implications of the expert selection of negatives (i.e., the NVML strategy) with a random sample strategy, where it is assumed that areas outside the favorable structural ellipses defined by NVML are negative. Because hydrothermal systems are sparse, it is highly probable that, in the absence of a favorable geological structure, hydrothermal favorability is low. We compare three training strategies: 1) the positive and negative labeled examples from NVML; 2) the positive examples from NVML with randomly selected negatives in equal frequency as NVML; and 3) the positive examples from NVML with randomly selected negatives reflecting the expected natural distribution of hydrothermal systems relative to the total area. We apply these training strategies to the NVML feature data (input data) using two ML algorithms (XGBoost and logistic regression) to create six favorability maps for hydrothermal resources. When accounting for the expected natural distribution of hydrothermal systems, we find that XGBoost performs better than the NVML neural network and its negatives. Model validation was less reliable using F1 scores, a common performance metric, than comparing probability estimates at known positives, likely because of the extreme natural class imbalance and the lack of negatively labeled sites. This work demonstrates that expert selection of negatives for training in NVML likely imparted modeling bias. Accounting for the sparsity of hydrothermal systems and all the types of locations without hydrothermal systems allows us to create better models for predicting hydrothermal resource favorability.

Geothermal Resources Council Transactions↗
USGS2025

Multi-scale geophysical mapping of the brine and bedrock surfaces along the Dolores River, Paradox Valley, Colorado, December 2023

Total dissolved solids derived from salt dome–sourced brine in the underlying alluvial aquifer substantially increase with distance in the reach of the Dolores River that passes through Paradox Valley in southwestern Colorado. The area has been the site of salinity control operations since the 1990s to reduce salt loading to the downstream Colorado River. Previous airborne and ground/water-based electromagnetic (EM) geophysical data have successfully mapped the top of the brine surface, albeit with relatively coarse near-surface resolution and limited spatial coverage. This present December 2023 study used ground-based high-resolution EM and passive seismic (horizontal-to-vertical spectral ratio, HVSR) tools to map in detail the depth and thickness of the brine zone in the alluvial aquifer (top of the brine down to bedrock contact) in areas immediately surrounding the Dolores River where previous airborne EM (AEM) results indicated brine within 10 m of land surface. Results indicate the deepest bedrock is generally associated with the shallowest brine and local depressions in the collapse breccia (caprock to the Paradox Formation salt) may facilitate vertical migration of brine into the alluvial aquifer. Additionally, the ground-based EM mapping corroborated general patterns in depth to brine that were observed in previous AEM results while also revealing additional detail, including suspected focused brine discharge zones to the Dolores River. A river-based survey wherein EM data, channel depth and river water–specific conductance information were collected augmented these findings and indicated specific areas on both the western and eastern sides of the river where focused brine discharge may occur. This study comprises a large-scale, ground- and water-based geophysical mapping effort, including hundreds of HVSR soundings and 100s of kilometres of EM data, which were successfully translated into digital brine and bedrock surfaces that could be incorporated into groundwater modelling efforts, future well siting or other decision-making.

Colorado↗
USGS2023

Sediment thickness map of United States Atlantic and Gulf Coastal Plain Strata, and their influence on earthquake ground motions

With the recent successful accounting of basin depth ground-motion adjustments in seismic hazard analyses for select areas of the western United States, we move toward implementing similar adjustments in the Atlantic and Gulf Coastal Plains by constructing a sediment thickness model and evaluating multiple relevant site amplification models for central and eastern United States seismic hazard analyses. We digitize and combine existing sediment thickness data sets into a composite surface that delineates the base of Cretaceous sediments under the Atlantic Coastal Plain and the base of Mesozoic sediments under the Gulf Coastal Plain. Amplification models dependent on sediment thickness, site natural period, and source-to-site path length are compared with data sets of observed ground motions to evaluate the ability of the new models to improve ground motion estimates. We find that the amplification models can account for observed trends in sediment-thickness and period-dependent residuals, but some tuning is required. For example, the model of Chapman and Guo requires a reference V S 30 , the time-averaged shear-wave velocity within 30 m of the Earth’s surface, for non-Coastal Plain sites, which we estimate to be between about 1 and 2 km/s. Along with our sediment thickness model, we estimate a velocity profile for application to the Harmon et al. site-natural-period-based model in order to best match the Chapman and Guo period dependence for a broad range of sediment thicknesses. The Next Generation of Attenuation models for the eastern United States Gulf Coast path-based adjustment models can also account for seismic attenuation in the Coastal Plain sediments and reduce the standard deviation of total residuals. If enacted in the U.S. Geological Survey National Seismic Hazard Model, these amplification models will reduce predicted short-period (<1 s) and increase predicted long-period (>1 s) ground motions in the Coastal Plains appreciably.

Earthquake Spectra↗