USGS Science⌕ Search

SEARCH · USGS Science

Results for “Integrated Environmental Assessment and Management”

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

Strong shaking from past Cascadia Subduction Zone earthquakes encoded in coastal landforms

Strong earthquakes along subduction zones are often devastating events, but sparse records along some tectonic margins limit our understanding of seismic hazards. Constraining shaking intensities is critical, especially in subduction zones with infrequent but large-magnitude earthquakes like the Cascadia Subduction Zone (CSZ), where the lack of recorded ground motions has led to uncertainty in the severity and potential impacts of future earthquakes. Here we fill this observational gap with a novel inventory of quantitative estimates of past shaking intensities from geotechnical modeling of coastal landforms. One hundred fifty-four deep-seated landslides and 65 fragile geologic features constrain minimum and maximum peak ground accelerations, respectively. These estimates are broadly consistent with model predictions of M9 ruptures, suggesting strong shaking of 0.4–0.8 g during past CSZ earthquakes. Local discrepancies between our geologic shaking constraints and earthquake simulations may inform past rupture behavior, leading to better predictions of shaking intensity for future earthquakes.

California, Oregon, Washington↗
USGS2024

Predicting hydrothermal reservoir depth from chemical geothermometers using a three-dimensional temperature model in the Great Basin, USA

Recent work in the Great Basin region of the western United States has made it possible to predict the depth of hydrothermal reservoirs (i.e., the depth at which heat is accumulated prior to ascent via hydrothermal upflow) identified through geochemistry and to contextualize the spatial patterns of these reservoir depths. Chemical geothermometers use the chemical and mineral constituents of hydrothermal fluids to predict the temperature at which fluids equilibrated with the host rocks at depth. Assuming that most of the Great Basin is dominated by conductive conditions until a vertically connected hydrothermal flow path is created (e.g., by faulting), geothermometers reflect the chemical and thermal conditions at the depth interval that the fluid has conductively equilibrated over a long period before a vertical conduit allows convective upflow. By pairing geothermometer temperature estimates with our recent three-dimensional temperature model of conductive heat flow in the Great Basin, we estimate the corresponding reservoir depths and construct a map of circulation depths. The predicted depths from geothermometers have spatial patterns across the Great Basin that relate to patterns seen in other geologic and geophysical data. Deeper springs generally occur disproportionately in areas with higher strain rates and in basins. We posit that current elevated strain rates reflect patterns of historic deformation where ongoing tectonic activity maintains permeable pathways to deeper reservoirs, some of which are estimated to exceed 6 km depth. Basins, as expected, contain a disproportionate number of these deep systems, because the underlying aquifers are closer to the surface in basins, thus requiring less water pressure to reach the surface than in mountain ranges. Most springs estimated to have their source in a deep reservoir occur at places known to host a hydrothermal system; these refined depth estimates of the source reservoir can help to better constrain the source depth for many known hydrothermal systems across the Great Basin.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2025

Long-term dynamics of earthquake swarms in the Yellowstone caldera

The factors controlling the spatial distribution and temporal evolution of earthquake swarms in volcanic systems remain unclear. We leverage leading-edge deep learning algorithms and a detailed three-dimensional velocity model to construct a 15-year high-resolution earthquake catalog of the Yellowstone caldera region. More than half of the region’s earthquakes are clustered into swarm-like families characterized by episodes of hypocenter expansion and migration. Adjacent earthquake swarms, separated by long quiescent periods, are found to be a dominant feature. We suggest that these swarms are controlled by the interplay between slowly diffusing aqueous fluids and rapid episodic fluid injections, which may result from the breaking of permeability seals. Our analyses also indicate that clustered seismicity beneath the caldera occurs on relatively immature, rougher fault structures, compared to more planar faults outside. Our results provide additional context for understanding seismicity in hydrothermal systems, highlighting the key role played by long-term fluid diffusion processes in driving the occurrence of earthquake swarms.

Idaho, Montana, Wyoming↗
USGS2024

Preliminary observations of the April 5th, 2024, Mw4.8 New Jersey earthquake

On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical Extreme Events Reconnaissance Association and the National Institute of Standards and Technology documented structural and nonstructural damage, including substantial damage to a historic masonry building in Lebanon, New Jersey. The USGS National Earthquake Information Center reported a focal depth of about 5 km, consistent with a lack of signal in Interferometric Synthetic Aperture Radar data. The focal mechanism solution is strike slip with a substantial thrust component. Neither mechanism’s nodal plane is parallel to the primary northeast trend of geologic discontinuities and mapped faults in the region, including the Ramapo fault. However, many of the relocated aftershocks, for which locations were augmented by temporary seismic deployments, form a cluster that parallels the general northeast trend of the faults. The aftershocks lie near the Tewksbury fault, north of the Ramapo fault.

New Jersey↗
USGS2024

The advantages of electric-grade heat storage in long closed-loop wells

Though long (e.g., tens of kilometers) closed-loop geothermal wells are under consideration for the extraction of subsurface heat, these wells might also serve as an efficient energy storage mechanism for electricity generation. Using a semi-analytic model solution, the potential for electric-grade heat storage as a function of ambient temperature (e.g., corresponds to depth of loop), well length, well diameter, and flow rate is evaluated. The following simplified cases for comparison of standard closed loop operation with energy storage operation were considered: [1] constant flow rate and constant temperature (90 °C) injection to represent the standard closed-loop base-case; and [2] constant flow rate and annual cycle sinusoidal temperature (90-150 °C) to represent seasonal (summer) charging while solar resources are peak. Initial temperatures for all scenarios considered herein are a uniform 175 °C. Electric-grade heat is assumed to be delivered whenever temperatures at the extraction point exceed 90 °C. For calculation purposes only, if temperature falls below 100 °C, it is assumed that heat delivered is sub-economic, so no electricity would be produced. For all scenarios, temperatures at the extraction well asymptotically approach the flow-weighted average injection temperature, but energy storage scenarios exhibited a damped time-varying signal that diminishes in magnitude with length of the loop. The asymptotic approach depends on initial temperatures in the rock and the heat extraction rate (a function of well diameter and flowrate). This analysis demonstrates that shorter closed loops can produce more electricity over time than longer closed-loops previously proposed for electricity production over typical engineering design lifetimes (e.g., 30 years). Although only a high-temperature scenario is considered herein, rock that is initially below boiling temperature would not host a standard closed-loop resource, but injection of hot water seasonally would asymptotically heat this low-temperature system to temperatures capable of electricity production. In other words, regardless of initial temperatures, closed loops could be used to store electricity with no critical minerals in the geothermal battery.

Conference Paper↗
USGS2025

What 25+ years of "Did You Feel It" intensities tell us about shaking in California

“When will the Big One happen?” is a question that people often have for earthquake scientists. But while waiting for the “Big One” to occur, people will usually experience frightening or damaging shaking from multiple relatively smaller‐magnitude earthquakes. Given this context, it raises the question: “Where does most of the damage come from?” Could smaller, yet more frequent, earthquakes account for the majority of reported impactful shaking? To explore this question, we consider reports of earthquake damage and felt shaking experiences from a catalog of community‐collected intensity values from the U.S. Geological Survey’s “Did You Feel It?” system. Comparing these intensities to expectations from a ground‐motion model, we find that earthquakes of magnitudes smaller than expected are responsible for most reported intensities of community decimal intensities (CDI) 4.5 and above (moderate and higher shaking intensity levels). (Here “expected value” is meant in its mathematical sense of the mean or equivalently median shaking intensity.) We also present a regional analysis of observed earthquake shaking for specific areas, detailing the maximum intensity experienced within a predetermined area. We identify several instances of M < 4.5 events that generated maximum intensities of CDI > 5 in regions in California surrounding the cities of Eureka, San Francisco, Los Angeles, and San Diego. Our results motivate the need to include smaller‐magnitude earthquakes in communications about earthquake hazard and risk reduction.

California↗
USGS2026

A Lake Charr pangenome reveals highly conserved Ohnologs as drivers of phenotypic diversity

Whole-genome duplication (WGD) is hypothesized to spur evolutionary diversification by producing genome-wide duplicate gene sets (Ohnologs) that are initially functionally redundant but can diverge markedly as the effects of relaxed selection accumulate over time. However, the underlying mechanisms remain unclear, in part because genomic studies often reconstruct Ohnolog evolution over millions of years, during which subsequent mutations can obscure deep-time signals. Investigating the relationship between Ohnolog evolution and diversification on a contemporary timescale offers clearer insights. We explore this relationship in Lake Charr ( Salvelinus namaycush ), where ∼10% of genes are retained highly conserved polyploid duplicates following the Salmonid-Specific Fourth Round WGD. Using 31 chromosome-level assemblies of Lake Charr from morphologically and ecologically diverse populations, joined into a pangenome graph, we characterized 189,555 structural variants (SVs) that were significantly less likely to affect genes retained as sequence-conserved Ohnolog pairs, nuancing the hypothesis that gene redundancy, relaxed selection, and functional diversification are intertwined. However, we found that SVs affecting such conserved Ohnologs may be potent drivers of adaptive evolution. Notably, we identified a putative 938-Kb interchromosomal translocation containing 25 genes with highly conserved Ohnologs in a paralogous (but untranslocated) genomic block. This putative translocation appears to have facilitated Ohnolog divergence in ankrd11 and hp , genes putatively linked to craniofacial and lipid metabolic diversity in sympatric Lake Superior morphs. This research reveals that conserved Ohnologs previously presumed to be redundant remain a reservoir for adaptive change.

BioRxiv↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS2025

The Mammoth magnetic anomaly, Pinal County, Arizona

A high-resolution Earth Mapping Resources Initiative airborne geophysical survey was flown in the southwest North American porphyry copper province to improve bedrock geologic maps and to identify areas that have unrecognized critical mineral resource potential. During the review of the aeromagnetic data, a distinctly monopolar-shaped, negative magnetic anomaly was observed at a flight elevation of 200 m above the ground with a maximum amplitude of –9500 nT. We have named this the Mammoth magnetic anomaly (MMA) because it is centered 12 km northeast of the town of Mammoth, Arizona, USA. The total field anomaly (TFA) contour of –500 nT enclosing the MMA defines an elongate shape measuring 2.5 km long by 1 km wide that trends northwest–southeast. Given the striking nature of this negative, monopolar-shaped magnetic anomaly, we conducted a ground campaign in May 2025 to determine its authenticity and potential relationship to critical mineral endowment. The MMA was confirmed on the ground with a TFA approaching –46,000 nT. Total magnetic intensity (TMI) observations routinely fell below the 18,000 nT operating floor of an industry-standard cesium-vapor total field magnetometer, and extremely low TMI measurements were corroborated along coincident traverse lines using two high dynamic range, but lower sensitivity, smartphone vector magnetometers. The lowest TMI values recorded by both smartphone magnetometers were 1000 nT and confirmed with multiple adjacent and crossing lines. Field observations suggest that this magnetic feature is caused by strong remanent magnetization within fine-grained magnetite hosted within locally altered Pinal Schist.

Arizona↗
USGS2025

From critical minerals to food security, the benefits of data collaboration

The volume of data in the public geoscience sphere is rapidly and continually expanding. At Geoscience Australia (GA) we saw an over 500% increase in data points within our relational databases between 2018 and 2024, over the life of the Exploring for the Future (EFTF) program. With the Resourcing Australia’s Prosperity initiative, a continued increase in data quantity will be seen for the next 10 to 35 years. At the same time, a broadening audience for geoscience data is increasing the desire to enhance the diversity of delivery streams. This ranges from data-dense highly technical outputs for geoscience specialists to curated interpretive products for people who are non-geoscientists. Development of these curated outputs has contributed to our awareness of the need for data to be collected and compiled in a way that ensures its reuse, with a focus on quality metadata and data provenance.

Conference Paper↗
USGS2024

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California↗
USGS2026

pySATSI: A Python package for computing focal mechanism stress inversions

We introduce pySATSI, a Python package for computing earthquake focal mechanism stress inversions. This algorithm can handle a wide variety of types of stress inversion problems with a single script and can duplicate many capabilities of preceding methodologies. We also add new capabilities that include spatiotemporally variable inversion grids, damped stress estimates for clusters with few or no focal mechanisms, and variable fault‐plane ambiguities that the user can assign to individual events. In addition, we added the ability to use damped stress inversions with fault‐plane ambiguity probabilities that are weighted by fault instabilities. Our algorithm is computationally efficient with faster runtimes than previous algorithms, scales well for large datasets, and can be easily parallelized.

Seismological Research Letters↗
USGS2025

Groundwater age estimates for the Mississippi River Valley alluvial aquifer based on tracer data collected during 2018–20

This study characterized groundwater age across the Mississippi River Valley alluvial aquifer (MRVA). Groundwater samples from 69 MRVA wells and 19 wells in Tertiary units of the Mississippi embayment aquifer system (MEAS) were analyzed for sulfur hexafluoride (SF 6 ), tritium ( 3 H), helium (He), and (or) carbon-14 of dissolved inorganic carbon ( 14 C). The age distributions of 89 samples were estimated by fitting lumped parameter models to processed tracer concentrations with the U.S. Geological Survey software TracerLPM. Mean ages of MRVA groundwater samples ranged from 12 to 22,000 years, with a median of 140 years. Mean ages of MEAS groundwater samples ranged from 230 to 52,000 or more years, with a median of 13,500 years. The spatial distribution of MRVA groundwater ages was found to be influenced by depth, inflow of groundwater from deeper units, and soil saturated hydraulic conductivity. In parts of the MRVA, the spatial distribution of MRVA groundwater ages was found to be influenced by annual recharge and (or) annual groundwater pumpage.

Alabama, Arkansas, Illinois, Kentucky, Louisiana, ↗
USGS2024

Intense alteration on early Mars revealed by high-aluminum rocks at Jezero Crater

The NASA Perseverance rover discovered light-toned float rocks scattered across the surface of Jezero crater that are particularly rich in alumina ( ~ 35 wt% Al 2 O 3 ) and depleted in other major elements (except silica). These unique float rocks have heterogeneous mineralogy ranging from kaolinite/halloysite-bearing in hydrated samples, to spinel-bearing in dehydrated samples also containing a dehydrated Al-rich phase. Here we describe SuperCam and Mastcam-Z observations of the float rocks, including the first in situ identification of kaolinite or halloysite on another planet, and dehydrated phases including spinel and apparent partially dehydroxylated kaolinite. The presence of spinel in these samples is likely detrital in origin, surviving kaolinitization, pointing to an ultramafic origin. However, the association of low hydration with increased Al 2 O 3 abundances suggests heating-induced dehydration which could have occurred during the lithification or impact excavation of these rocks. Given the orbital context of kaolinite-bearing megabreccia in the Jezero crater rim, we propose an origin for these rocks involving intense aqueous alteration of the parent material, followed by dehydration/lithification potentially through impact processes, and dispersion into Jezero crater through flood or impact-related processes.

Communications Earth & Environment↗
USGS2026

Deformed submarine terraces in Puget Sound, Pacific Northwest, indicate only one M >~7.5 earthquake on the Seattle fault zone in the past 11,000 yr

Submerged marine terraces in Puget Sound, deformed across the Seattle fault zone (SFZ), indicate that only one earthquake as large as M~7.5 has occurred in at least the past 11 kyr. Previous paleoseismic studies document a M~7.5 earthquake between 923–4 CE, which uplifted coastal marine terraces by as much as 8 m. We demonstrate that this earthquake was the only such event since ~11 ka by mapping and quantifying deformation of older marine terraces, now submerged in Puget Sound. The submerged terraces, attributed to a late-glacial sea-level lowstand, record both glacial isostatic rebound and tectonic deformation. Vertical offset of the ~11 ka terraces within the SFZ is comparable to that of the marine terraces uplifted in 923 CE, implying no additional large (M>~7.5) earthquake on the SFZ since ~11 ka. This result implies a longer recurrence interval than current hazard estimates, which assumes recurrence of M>7.1 events every 5 kyr. Our mapping of SFZ deformation since ~11 ka also supports fault segmentation and contiguous block uplift between the Seattle and Tacoma fault zones.

Washington↗
USGS2025

Aftershocks in stress shadows are inconsistent with modeled static Coulomb stress changes

Aftershock triggering is commonly attributed to increases in static Coulomb stress. In some areas, termed "stress shadows", a decrease in Coulomb stress is predicted to suppress earthquake occurrence. However, aftershocks are often observed in the modeled stress shadows. We examine several hypotheses that attempt to reconcile these shadow aftershocks with the static Coulomb stress change model: (1) they appear to be in shadows because of inaccuracy in the stress change calculations, (2) they occur on faults of unusual orientation which actually experienced increased Coulomb stress, (3) they occur on faults with different frictional properties, not modeled well by Coulomb stress, and (4) they are secondary aftershocks triggered by prior aftershocks or afterslip. When tested on the 2016 Mw7.0 Kumamoto, Japan, and 2019 Mw7.1 Ridgecrest, California, aftershock sequences, none of these hypotheses can explain the majority of the shadow aftershocks, and taken together these hypotheses can explain only about half of these aftershocks. This implies that Coulomb stress modeling that lacks small-scale fault zone heterogeneity might be inadequate to fully capture the true static stress changes and/or that other physical triggering models are needed, for example transient processes such as delayed triggering by dynamic stress changes from the passing seismic waves.

California↗
USGS2025

Niobium-rich minerals from the Sheep Creek carbonatite-related deposits, Montana, USA

The Sheep Creek area in southwest Montana contains a number of small, tabular, carbonatite-related deposits that are rich in Nb- and REE-minerals. Two of the “vein-dikes” were mined at a small scale for Nb in the late 1950s. The most abundant Nb-bearing minerals are columbite-(Fe), nioboaeschynite-(Ce), baotite, and two varieties of Ti-Nb-Fe-oxide. An early oxide with Ti/Nb (atomic) near unity later broke down to a mixture of columbite + TiO2 with distinctly lower Nb content. Zircon intergrown with the Nb-REE minerals was dated by U-Pb methods to 116±2 Ma. This age is in agreement with previous dating (Th-Pb on monazite, Re-Os on molybdenite), and indicates that the Sheep Creek deposits are early Cretaceous in age.

Montana↗
USGS2026

Velocity-independent dry friction on mica: A realization of ideal Amontons-Coulomb friction

The Amontons-Coulomb friction law assumes that the frictional force between materials is independent of sliding velocity. However, as Coulomb noted, this is a rough approximation, and a second-order dependence of friction on the logarithm of sliding velocity is incorporated in a commonly used ‘rate- and state-dependent’ friction representation. Here we conduct shear experiments on mica, a layer-structured mineral, at temperatures ranging from 25 to 200ºC and under normal stress of 100 MPa. The friction coefficient clearly depends on the logarithmic sliding velocity at 25ºC, but rate sensitivity decreases with increasing temperature until at 200ºC, the friction coefficient is independent of sliding velocity. Our findings could initiate the development of velocity-independent frictional materials, realizing the ideal Amontons-Coulomb friction.

Physical Review Letters↗