USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Toxicology and Chemistry”

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

ShakeAlert Earthquake Early Warning System performance during the Mw 7.0 offshore Cape Mendocino earthquake

The 5 December 2024 M w 7.0 Offshore Cape Mendocino earthquake was a challenging test of the U.S. West Coast ShakeAlert earthquake early warning system due to its offshore epicenter and limited near‐source station coverage. We analyzed real‐time performance of all components of the ShakeAlert system, including the seismic algorithms (earthquake point‐source integrated code [EPIC] and Finite‐fault rupture Detector [FinDer]), the geodetic algorithm (Geodetic First Approximation of Size and Time–peak ground displacement [GFAST‐PGD]), and network telemetry during the event. EPIC created the first solution for this earthquake 15 s after origin time with an initial magnitude estimate of M 5.6 and location error of 10 km from the Advanced National Seismic System epicenter. An early spurious trigger from station CE.89101 fortuitously maintained location accuracy and, correspondingly, magnitude accuracy. FinDer contributed its first solution at 18 s with a location estimate closer to the seismic network and produced two distinct rupture geometries, leading to minor fluctuations in estimated intensity contours. GFAST‐PGD did not meet alerting thresholds but otherwise performed as expected. Network latencies were <2 s for most stations, supporting the rapid detection of this earthquake by the system. Roughly five million alerts were delivered to cell phone devices in California and Oregon during this event. This was also the first instance of a school district‐wide ShakeAlert‐powered system being activated. Comparisons to recorded seismograms demonstrate that the maximum warning times before potentially damaging shaking (intensity 6+) were in the range of 5–55 s. Although the ShakeAlert system provided accurate solutions and useful alert delivery, this earthquake raised awareness of potential issues within the system, including the need for improved offshore location estimates, a combination of solutions from ShakeAlert servers, and handling of spurious triggers.

California↗
USGS2024

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington↗
USGS2024

Time varying crustal anisotropy at Whakaari/White Island volcano

Whakaari/White Island has been the most active New Zealand volcano in the 21st century, producing small phreatic and phreatomagmatic eruptions, which are hard to predict. The most recent eruption occurred in 2019, tragically claiming the lives of 22 individuals and causing numerous injuries. We employed shear-wave splitting analyses to investigate variations in anisotropy during quiescence, unrest, and an eruption between 2018 and 2020. We examined spatial and temporal variations in 7972 shear-wave splitting measurements and derived 1402 V p /IV s ratio measurements. Comparing shear-wave splitting parameters from similar earthquake paths across different times indicates that the observed temporal changes are unlikely to result from variations in earthquake paths through media with spatial variability. Instead, these changes may stem from variations in anisotropy over time, likely caused by changes in crack alignment due to stress or varying fluid content.

Whakaari/White Island volcano↗
USGS2025

Framework for assessing intranational risk of critical mineral sectors: Case study of India

Recent geopolitical actions regarding the supply of critical minerals have highlighted both the threats posed by countries dominating a critical mineral supply chain, and the potential mitigations served by alternative supplier countries. However, no systematic quantification exists for the intranational factors that determine such threats or mitigations. This study outlines a methodological framework for quantifying the intranational risk profile of critical mineral sectors and compares them with the sector's international potential for leverage. The four intranational factors of concentration, synchronization, private sector involvement, and regulatory obstacles were analyzed across multiple scales in order to determine an intranational risk score that quantifies the resilience or fragility of the critical mineral sector in question. Network analysis was also separately conducted to identify the most important entities operating across all critical mineral sectors. Finally, this study assessed the critical mineral sector's global share of reserves and production to determine an international risk score that quantifies the country's critical mineral sector's leverage potential. A case study of India showed that two critical mineral sectors with high leverage potential, chromite and barite, were respectively intranationally resilient and fragile. Of the remaining sectors with minimal leverage, five sectors such as aluminum were intranationally resilient and may serve as viable mitigation options for future critical mineral supply. Three sectors such as titanium were intranationally fragile, but recent regulatory developments may both decrease intranational risk and increase international risk, potentially making them major critical mineral sectors in the future.

Resources Policy↗
USGS2026

Evidence for a biological origin of uranium-rich carbon masses within the Ediacaran Salt Range Formation of Pakistan

Thucholites are unique organic structures found in igneous and sedimentary rocks composed of a U-C-rich interior enclosed by an organic outer shell. Their formation and occurrence have perplexed scientists for over 100 years. Typically, thucholites are sparse in sedimentary rocks but where found in abundance, they may be the result of rapid paleoecological disruptions, e.g., volcanic ashfall. Here, we evaluated thucholites from the Ediacaran Salt Range Formation of the Indus Basin in Pakistan using field emission scanning electron microscopy (FESEM) and nanoscale secondary ion mass spectroscopy (NanoSIMS) to propose a two-stage mechanism for thucholite development. NanoSIMS results suggest organic matter in thucholite cores formed by biological fractionation (represented by the presence of lighter 13 C/ 12 C and 34 S/ 32 S) while the outer organic mantle formed via radiolysis-induced polymerization. FESEM elemental analysis confirms compositional differences between the two thucholite components (core and mantle), further implying their contrasting origins. Dimensional comparison of thucholite cores and their U-bearing mineral morphologies to ancient and modern U-biomineralized microbes suggests that the thucholite cores formed from metal-biological interactions. The presence of volcanogenic biotite, alkali feldspar, and clay spherules (interpreted as devitrified volcanic glass) with thucholite suggests that the thucholite cores in the Salt Range Formation are biological responses to rapid paleoenvironmental change from volcanic eruption which preserved the fossilized morphologies of Proterozoic microorganisms.

Salt Range↗
USGS2025

The 1912 Ms 7.2 earthquake in the Denali region of central Alaska

The 2002 M w 7.9 earthquake in central Alaska ruptured 340 km on three faults—Susitna Glacier thrust fault, Denali fault, Totschunda fault—crossing both the Richardson Highway and the Alaska Pipeline. Its occurrence prompted renewed interest in historical large earthquakes that possibly originated on the Denali fault. One of these earthquakes was a M s 7.2 event on July 7, 1912, which we revisit with two approaches: (1) probabilistic relocation of the epicenter using globally recorded arrival times, and (2) compilation and reassessment of shaking intensity reports to estimate a macroseismic epicenter. Our preferred instrumental epicenter is west of the Parks Highway and in agreement with the maximum‐reported shaking, which was from the Parker–Browne expedition of Denali. We also relocated a M s 6.4 aftershock, whose epicenter is 11 km from the mainshock. Candidate faults for the 1912 earthquake include the Denali fault, the McLeod Creek thrust fault, and the Kantishna Hills thrust fault. Future analysis of active faults, paleoseismic results, 1912 instrumental data, and 1912 felt reports may help in interpreting the fault and mechanism of the 1912 earthquake.

Alaska↗
USGS1983

Technique for estimating depth of floods in Tennessee

Estimates of flood depths are needed for design of roadways across flood plains and for other types of construction along streams. Equations for estimating flood depths in Tennessee were derived using data for 150 gaging stations. The equations are based on drainage basin size and can be used to estimate depths of the 10-year and 100-year floods for four hydrologic areas. A method also was developed for estimating depth of floods having recurrence intervals between 10 and 100 years. Standard errors range from 22 to 30 percent for the 10-year depth equations and from 23 to 30 percent for the 100-year depth equations.

Tennessee↗
USGS2024

The new self-anchored suspension (SAS) San Francisco Bay Bridge- Its response to a small earthquake

This paper presents a summary of previously published work (Celebi 2023) related to the new Self-Anchored Suspension (SAS) bridge that went into service within the last decade as a replacement for the older truss bridge spanning between Yerba Buena Island and Oakland, California, within the San Francisco Bay Area. During the October 19, 1989 M6.9 Loma Prieta earthquake, which occurred ~100 km south of the Bay Bridge, a section of the upper deck of the truss bridge fell onto the lower deck – thus closing this important lifeline between San Francisco and Oakland. The SAS is unique, self-anchored, and suspended by a single tower that is pivotal in trafficking the cable and hanger system to support the decks. The SAS bridge is extensively instrumented by the California Geological Survey’s Strong Motion Instrumentation Program (CSMIP). There are approximately 85 channels of accelerometers in the seismic monitoring system that recorded the October 14, 2019 Mw4.6 Pleasant Hill earthquake. The data allow a complex but identifiable coupled response of the deck, tower, and cable system. Both acceleration and displacement time-history data are used to extract significant frequencies using system identification methods, including spectral analyses. Results are compared to those from finite-element-model (FEM) analyses carried out during the design and analysis process of the bridge in 2002 (Nader et al. 2002). There are differences between FEM analyses results and those from the low amplitude shaking caused by a seismic event. An apparent frequency (period) of the SAS bridge is assessed (approximately 5.2 seconds). In a plot of deck length versus period, there is an almost linear relationship with periods of other regular suspension bridges, such as the Golden Gate Bridge and the Carquinez Bridge, both in the San Francisco Bay.

California↗
USGS2025

Towards mobile wind measurements using joust configured ultrasonic anemometer for applications in gas flux quantification

Small uncrewed aerial systems (sUASs) can be used to quantify emissions of greenhouse and other gases, providing flexibility in quantifying these emissions from a multitude of sources, including oil and gas infrastructure, volcano plumes, wildfire emissions, and natural sources. However, sUAS-based emission estimates are sensitive to the accuracy of wind speed and direction measurements. In this study, we examined how filtering and correcting sUAS-based wind measurements affects data accuracy by comparing data from a miniature ultrasonic anemometer mounted on a sUAS in a joust configuration to highly accurate wind data taken from a nearby eddy covariance flux tower (aka the Tower). These corrections had a small effect on wind speed error, but reduced wind direction errors from 50° to >120° to 20–30°. A concurrent experiment examining the amount of error due to the sUAS and the Tower not being co-located showed that the impact of this separation was 0.16–0.21 ms &#x2212; 1 "> ms − 1 , a small influence on wind speed errors. Lower wind speed errors were correlated with lower turbulence intensity and higher relative wind speeds. There were also some loose trends in diminished wind direction errors at higher relative wind speeds. Therefore, to improve the quality of sUAS-based wind measurements, our study suggested that flight planning consider optimizing conditions that can lower turbulence intensity and maximize relative wind speeds as well as include post-flight corrections.

Alaska↗
USGS2025

Approximate inland extent of saltwater intrusion at the base of the Biscayne aquifer, Miami-Dade County, Florida, 2022

Miami-Dade County is part of a densely populated urban corridor in southeastern Florida. The Biscayne aquifer serves as Miami-Dade County’s primary drinking water source and is characterized by highly permeable karstic limestone and carbonate sand. The aquifer’s coastal location and permeable nature make it susceptible to saltwater intrusion. Monitoring the current inland extent and the rate of movement of the saltwater front in the aquifer can inform management strategies for conserving the long-term sustainability of the county’s water supply. In the 1950s, the U.S. Geological Survey published a map of the inland extent of saltwater intrusion in the Biscayne aquifer and has continued to update this map to monitor changes over time, with the most recent update published in 2018. An updated map has been created showing the approximate inland extent of saltwater intrusion in the Biscayne aquifer in eastern Miami-Dade County in 2022, with the 2018 extent shown for comparison. The inland extent of saltwater intrusion was mapped through the interpretation of borehole electromagnetic induction logs and measurements of chloride and specific conductance in groundwater samples. The location of the saltwater interface at the base of the Biscayne aquifer was represented by the 1,000-milligram-per-liter isochlor. This report describes changes in the location of the saltwater interface from 2018 to 2022. By 2022, the saltwater interface had moved farther inland in both the northern and southern parts of the county, advancing by as much as 0.3 kilometer in the north and up to 0.8 kilometer in the Model Land Area to the south. However, it remained relatively unchanged from its 2018 position in the east-central part of the county.

Florida↗
USGS2026

A geographic dataset of rocky reefs habitat areas of particular concern for the United States West Coast

The United States National Marine Fisheries Service determines “essential fish habitat (EFH)” for federally managed species in coordination with regional fishery management councils, considers adverse effects to those habitats, and provides information to further habitat conservation and enhancement. Identifying discrete subsets of EFH as “habitat areas of particular concern (HAPC)” can help focus conservation, management, and research efforts. In 2006, the Pacific Fishery Management Council designated rocky reefs along the United States (U.S.) West Coast as HAPCs for groundfishes because of their ecological significance, sensitivity to human impacts, and relative rarity. To better understand where rocky reefs occur, we (1) located rocky reef areas that were not included in the 2006 rocky reef dataset, and (2) incorporated best available data into a refined geographic dataset that enables visualization. Our update shows that rocky reefs are distributed throughout the U.S. West Coast continental margin, are patchier than previously known, and comprise 8% of the extent of all data inputs. This updated dataset will inform resource management decisions in coastal and marine environments.

California, Oregon, Washington↗
USGS2025

Origin of the high Pd/Pt ratio of the J-M Reef, Stillwater Complex Montana USA

The J-M Reef of the Stillwater Complex exhibits a high and consistent Pd/Pt ratio (~3.8). This ratio results from the equilibration of an immiscible sulfide liquid with a relatively high Pd/Pt silicate melt rather than an unusually Pd- and Pt-enriched parental melt. Numerical modeling suggests that the original silicate melt contained typical mantle-derived concentrations of Pd and Pt (~10–20 ppb Pd/Pt ~1). The partitioning of Pt and Pd between sulfide liquid and silicate melt alone cannot explain the consistently high Pd/Pt ratio across variable melt-to-sulfide mass ratios (R factors). Instead Pt-depletion caused by the early fractionation of Pt-alloy from S-undersaturated silicate magma likely established the high Pd/Pt signature. High Pd/Pt ratios can form through batch equilibration of sulfide liquid with silicate melt if partition coefficients are extremely high (>10⁶). Alternatively Pd enrichment may result from sulfide upgrading within the resident footwall mush under smaller partition coefficients (10⁴–10⁶) in this model the instantaneous R factors remain low (R ≈ 100–700). This limits the impact of Pt and Pd partitioning on sulfide composition and helps explain the Pd-enriched character of the J-M Reef.

Montana↗
USGS2025

Cooling perspectives on the risk of pathogenic viruses from thawing permafrost

Climate change is inducing wide-scale permafrost thaw in the Arctic and subarctic, triggering concerns that long-dormant pathogens could reemerge from the thawing ground and initiate epidemics or pandemics. Viruses, as opposed to bacterial pathogens, garner particular interest because outbreaks cannot be controlled with antibiotics, though the effects can be mitigated by vaccines and newer antiviral drugs. To evaluate the potential hazards posed by viral pathogens emerging from thawing permafrost, we review information from a diverse range of disciplines. This includes efforts to recover infectious virus from human remains, studies on disease occurrence in polar animal populations, investigations into viral persistence and infectivity in permafrost, and assessments of human exposure to the enormous viral diversity present in the environment. Based on currently available knowledge, we conclude that the risk posed by viruses from thawing permafrost is no greater than viruses in other environments such as temperate soils and aquatic systems.

mSystems↗
USGS1998

Ground-water flow in the surficial aquifer system and potential movement of contaminants from selected waste-disposal sites at Cecil Field Naval Air Station, Jacksonville, Florida

As part of the Installation Restoration Program, Cecil Field Naval Air Station, Jacksonville, Florida, is considering remedialaction alternatives to control the possible movement of contaminants from sites that may discharge to the surface. This requires a quantifiable understanding of ground-water flow through the surficial aquifer system and how the system will respond to any future stresses. The geologic units of interest in the study area consist of sediments of Holocene to Miocene age that extend from land surface to the base of the Hawthorn Group. The hydrogeology within the study area was determined from gamma-ray and geologists' logs. Ground-water flow through the surficial aquifer system was simulated with a seven-layer, finite-difference model that extended vertically from the water table to the top of the Upper Floridan aquifer. Results from the calibrated model were based on a long-term recharge rate of 6 inches per year, which fell in the range of 4 to 10 inches per year, estimated using stream hydrograph separation methods. More than 80 percent of ground-water flow circulates within the surficial-sand aquifer, which indicates that most contaminant movement also can be expected to move through the surficial-sand aquifer alone. The surficial-sand aquifer is the uppermost unit of the surficial aquifer system. Particle-tracking results showed that the distances of most flow paths were 1,500 feet or less from a given site to its discharge point. For an assumed effective porosity of 20 percent, typical traveltimes are 40 years or less. At all of the sites investigated, particles released 10 feet below the water table had shorter traveltimes than those released 40 feet below the water table. Traveltimes from contaminated sites to their point of discharge ranged from 2 to 300 years. The contributing areas of the domestic supply wells are not very extensive. The shortest traveltimes for particles to reach the domestic supply wells from their respective contributing areas ranged from 70 to 200 years.

Florida↗
USGS1983

Use of digital land-cover data from the Landsat satellite in estimating streamflow characteristics in the Cumberland Plateau of Tennessee

Characteristics of Streamflow at ungaged sites in coal-mining areas of the Cumberland Plateau may be estimated with regression equations that relate streamflow characteristics to the physiographic and climatic characteristics of the corresponding drainage basins. An experiment was performed to compare the accuracy of equations using basin characteristics derived from maps and climatological records (control group equations) with the accuracy of equations using basin characteristics derived from digital processing of Landsat spectral data as well as maps and climatological records (experimental group equations). Results of this experiment show that (with the exception of low flows and four of six annual peak logarithms) drainage area can explain more than 90 percent of the variance in all streamflow characteristics in both groups of equations. Seventeen of 39 experimental group equations that have two basin characteristics each are different from the corresponding control group equations. Five of the 17 differing experimental group equations have no counterpart in the control group because of lack of significance in control group basin characteristics, seven are slightly more accurate, and five are not measurably different in accuracy. When the equations in both groups are arranged into six flow categories, there is no substantial difference in accuracy between equations using basin characteristics derived from maps and climatological records (control group) and equations using basin characteristics derived from Landsat tapes as well as maps and climatological records (experimental group) for this particular study area, the Cumberland Plateau of Tennessee.

Tennessee↗
USGS2026

Ice sheet dynamics drive pronounced changes in the subsurface freshwater-saltwater interface

Saltwater is migrating into freshwater aquifers globally with water quality and biogeochemical implications, yet saltwater intrusion in glaciated regions is sparsely investigated. Field observations suggest that groundwater head in glaciated systems is influenced by ice sheet forcings and provides evidence that seawater infiltrated into offshore aquifers during past deglaciation events. To understand links between ice sheet dynamics, groundwater head, and saltwater intrusion, we use numerical models to explore the effects of deglaciation on nearshore head and salinity distributions. We find that ice sheet thinning diminishes groundwater head, and the resulting shift in subsurface pressure gradients drives rapid landward movement of the subsurface freshwater-saltwater interface up to 4.0 km or 1.3 m per m ice sheet loss. Results highlight an overlooked saltwater intrusion mechanism that aligns with field observations and affects glaciated coastlines undergoing ice sheet retreat, underscoring the need to consider this mechanism in studies of contemporary coastal water quality.

Geophysical Research Letters↗
USGS2024

On connecting hydro-social parameters to vegetation greenness differences in an evolving groundwater-dependent ecosystem

Understanding groundwater-dependent ecosystems (i.e., areas with a relatively shallow water table that plays a major role in supporting vegetation health) is key to sustaining water resources in the western United States. Groundwater-dependent ecosystems (GDEs) in Colorado have non-pristine temporal and spatial patterns, compared to agro-ecosystems, which make it difficult to quantify how these ecosystems are impacted by changes in water availability. The goal of this study is to examine how key hydrosocial parameters perturb GDE water use in time and in space. The temporal approach tests for the additive impacts of precipitation, surface water discharge, surface water mass balance as a surrogate for surface–groundwater exchange, and groundwater depth on the monthly Landsat normalized difference vegetation index (NDVI). The spatial approach tests for the additive impacts of river confluences, canal augmentation, development, perennial tributary confluences, and farmland modification on temporally integrated NDVI. Model results show a temporal trend (monthly, 1984–2019) is identifiable along segments of the Arkansas River at resolutions finer than 10 km. The temporal impacts of river discharge correlate with riparian water use sooner in time compared to precipitation, but this result is spatially variable and dependent on the covariates tested. Spatially, areal segments of the Arkansas River that have confluences with perennial streams have increased cumulative vegetation density. Quantifying temporal and spatial dependencies between the sources and effects of GDEs could aid in preventing the loss of a vulnerable ecosystem to increased water demand, changing climate, and evolving irrigation methodologies.

Colorado↗
USGS2024

Characteristics of the fault damage zone From high-resolution seismic imaging along the Palos Verdes Fault, California

The distribution and intensity of fault damage zones provides insight into fault activity and its relationship to fluid flow in the crust. Presently, measures of the in-situ distribution of fault damage remain limited and along-strike studies are rare. This study focuses on an offshore section Palos Verdes Fault damage zone that spans 28 km, near Los Angeles, California. To investigate the previously unresolved shallow (∼400 m below the seafloor) fault damage zone we use densely spaced (∼500 m line separation) newly collected sparker multichannel seismic lines and sub-bottom profiles. The combination of high-resolution acquisition methods and specialized seismic processing workflows provide improved imaging of shallow faulting. We apply a multi-trace similarity technique to identify discontinuities in the seismic data that may be attributed to faults and fractures. This fault detection approach reveals diverse fault damage patterns on adjacent seismic profiles. However, a discernible damage zone pattern emerges by stacking multiple damage detection profiles along strike. We find that peak damage identified in this way corresponds to the active main fault strand, confirmed in this study, and thus the technique may be useful for identifying active fault strands elsewhere. Additionally, we observe that the variable width of the damage zone along strike is controlled by fault obliquity. Furthermore, our observations reveal a correlation between fault damage and seafloor fluid seeps visible in the water column, suggesting that damage plays a role in controlling fluid flow around the fault.

California↗