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 88Linked to original sources
USGS1976

Diamond drilling at the Ma'milah Mine, Kingdom of Saudi Arabia

The Ma'millah gold mine is about 90 km southeast of At Ta'if at lat 21°03'N., long 41°18'E., in southwestern Saudi Arabia. The deposit, which was worked extensively by ancient miners, consists of several veins and massive lenses of quartz along a regional fault. The extensive ancient mine workings, large size of some quartz outcrops, and significant gold content of several samples indicated that subsurface exploration was warranted. Accordingly, two holes were drilled into the deposit during the interval December 1972 to February 1973. The drill findings were disappointing; the grade of vein material was found to be too low to be of economic interest and the quartz bodies were found to pinch out both along strike and downdip. The drilling program therefore was terminated. The deposit is not worth further exploration

Ma'milah Mine↗
USGS2024

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications↗
USGS2025

Submarine avalanche deposits hold clues to past earthquakes

Earthquakes and other natural events sometimes shake the seafloor near coastlines severely enough to cause underwater avalanches that rush down steep slopes, scouring the seabed and carrying sediment to greater depths. These fast-moving sediment-laden flows, called turbidity currents , have at times damaged underwater infrastructure like pipelines and communications cables, as they did, for example, in snapping transatlantic cables off the coast of Newfoundland after the 1929 Grand Banks earthquake.

EOS Transactions↗
USGS2025

An unexplained tsunami: Was there megathrust slip during the 2020 Mw7.6 Sand Point, Alaska, earthquake?

On October 19, 2020, the M w 7.6 Sand Point earthquake struck south of the Shumagin Islands in Alaska. Moment tensors indicate the earthquake was primarily strike-slip, yet the event produced an enigmatic tsunami that was larger and more widespread than expected for an earthquake of that magnitude and mechanism. Using a suite of hydrodynamic, seismic, and geodetic modeling techniques, we explore plausible causes of the tsunami. We find that strike-slip models consistent with the moment tensor orientation cannot produce the observed tsunami. Hydrodynamic inversion of sea surface deformation from deep ocean and tide gauge data suggest seafloor deformation more closely matches a megathrust, rather than a strike-slip, source. Static slip inversions, using sea level and Global Navigation Satellite System data, allow for a portion of co-seismic megathrust slip that can explain tsunamigenesis. Combining all available geophysical datasets to model the kinematic rupture, we show that considerable, relatively slow, megathrust slip is allowable in the Shumagin segment, concurrent with strike-slip faulting. We hypothesize that the slow megathrust rupture does not contribute much seismic radiation allowing it to previously go unnoticed with traditional seismic monitoring.

Alaska↗
USGS2025

To heal or not to heal?: 2. The moment-recurrence time behavior of repeating earthquakes in the 2011 Prague, Oklahoma aftershock sequence is consistent with laboratory healing rates

The timing and failure conditions of an earthquake are governed by the interplay between fault reloading and restrengthening. The moment-recurrence time behavior of repeating earthquakes can give observational estimates of fault healing rates; however, it is difficult to link these observed healing rates to laboratory studies of frictional healing in part because of uncertainty in lithology. Here, we study the 2011 Prague earthquake sequence, which includes repeating earthquakes in the Arbuckle group and the granitic basement, and compare them to laboratory experiments on samples of the Arbuckle and Troy granite (representative of the basement rock) (Okamoto et al., 2025, https://doi.org/10.1029/2024JB030573 ). We find three spatially distinct groups of repeating earthquakes with different moment-recurrence behavior: (a) constant moment-recurrence time in the Arbuckle group, (b) scattered moment-recurrence time at the intersection of the foreshock-mainshock fault in the granitic basement, and (c) moment-predictable behavior outside of the foreshock-mainshock fault intersection also in the granitic basement. Our observation of stagnant healing for repeating sequences in the Arbuckle group is consistent with laboratory observations of low healing rates for moderately high pore fluid pressures in Arbuckle samples. For the moment-predictable group, the source radius that is required in order to match healing rates is consistent with source radius estimations when taking into account reasonable attenuation of the P- pulse width. Overall, we observe diverse healing behaviors in the seismic families that are consistent with laboratory healing rates, providing seismic evidence that contact-scale frictional mechanisms are relevant to large-scale earthquake dynamics.

Oklahoma↗
USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗
USGS2024

Development of USGS NSHMs: Do small changes in hazard imply small changes in risk?

One of the flagship products from the U.S. Geological Survey (USGS) is the National Seismic Hazard Model (NSHM). Since 1976, the NSHM has been periodically updated to reflect newly published earthquake science and provide probabilistic estimates of seismic hazard for the United States. During each update cycle, alternative models are deliberated, analyzed, and documented through logic trees and their corresponding logic tree branch weights. For example, the decision to modify a logic tree branch weight may be influenced by sensitivity analyses of the logic tree branches in their effects on the mean hazard. However, do small changes in traditional measures of hazard imply small changes in risk? In this study, we make use of two update cycles of the USGS NSHMs and a National Bridge Inventory (NBI) from the Federal Highway Administration (FHWA) to explore the preceding question. Specifically, we first identify geographic locations in the conterminous United States in which the change in hazard from one cycle to another is relatively small. Next, we model the seismic risk to highway bridges for these locations and for each update cycle, while simultaneously distinguishing low hazard environments from high hazard environments. These data enable quantitative analysis of how much changes in risk can be expected given small changes in hazard, investigating the importance of risk models in decision-making during development of the USGS NSHMs.

Conference Paper↗
USGS2026

Beaver dam analogs as nature-based solutions to mitigate snowpack loss in northern New Mexico

Reductions in snow have left many streams in northern New Mexico dry or with very low flows during the summer. Base flow, defined as the contribution of groundwater to streamflow, can sustain streamflow during the summer and during drought conditions. Beaver dam analogs can be used to “slow the flow,” or increase the infiltration of rain-based runoff to replenish shallow groundwater reservoirs and increase base flow and potentially increase summertime streamflow.

New Mexico↗
USGS2024

Predicted exposure of communities in southeastern United States to climate-related coastal hazards

A rigorous analysis of 21st Century multi-hazard exposure for U.S. Southeast Atlantic coastal communities indicates that up to 70% of residents will be exposed daily to shallow and emerging groundwater by ~2100, 15 times higher than from surficial flooding alone. This threat further exacerbates other coastal stressors, such as flooding, subsidence, and beach erosion, that impact these vulnerable coastal communities and ecosystems under expected future climate change scenarios.

Florida, Georgia, North Carolina, South Carolina, ↗
USGS2025

U.S. Geological Survey Groundwater Climate Response Network, 2024

As of October 2024, the U.S. Geological Survey (USGS) operated 588 sites across the United States and its territories as part of the Groundwater Climate Response Network (CRN). The CRN is comprised of wells selected to monitor the effects of climate variability, such as droughts, on groundwater levels nationwide. The CRN includes nearly 500 locations with real-time data and more than 100 sites with non-real-time data available to the public on the CRN web mapper and the USGS National Water Dashboard.

General Information Product↗
USGS2025

Migration of seismicity from the mantle to the upper crust beneath Harrat Lunayyir volcanic field, Saudi Arabia

Harrat Lunayyir is a volcanic field in Saudi Arabia that experienced a Mw~5.4 earthquake driven by an upper-crustal dike intrusion in May 2009. This volcanic field has exhibited numerous forms of volcanic seismicity both prior to and since the 2009 dike intrusion. Significantly, earthquakes within the lithospheric mantle and, rarely, the lower crust are present in the two-decade long seismicity catalog of Harrat Lunayyir. Here we analyze 24 years of volcanic seismicity at Harrat Lunayyir from 1998 to 2022. We find that: 1) precursory seismicity began at least eight years prior to the 2009 event, with a particularly notable seismic episode one year prior; 2) lithospheric mantle seismicity is highly localized in space and in time, largely occurring in discrete sequences lasting on the order of a few hours to a few days; 3) one seismic sequence clearly migrates upward from the lithospheric mantle to the upper crust, including seismicity within the nominally ductile lower crust; 4) crustal seismicity has been slowly declining over time; and 5) lithospheric-mantle seismicity does not show any apparent decline with time. From these observations we infer that the seismicity is driven by magmatic fluids or volatiles, and seismic monitoring of this volcanic field should continue into the future.

Harrat Lunayyir volcanic field↗
USGS2024

On algorithmically determined versus traditional macroseismic intensity assignments

The utility of macroseismic data, defined as the effects of earthquakes on humans and the built environment, has been increasingly recognized following the advent of online systems that now produce unprecedented volumes of macroseismic intensity information. Contributed reports from the U.S. Geological Survey “Did You Feel It?” (DYFI) system ( Wald et al ., 1999 ) are used to generate intensity values with an algorithm based on seminal work by Dengler and Dewey (1998) . The algorithm was developed initially to reproduce intensity values assigned by expert opinion using questionnaire results collected by telephone survey. In this article, I discuss reasons why intensity values from (self‐selected) DYFI responses can differ from values that would be assigned by expert opinion given more complete data from randomly selected participants. For example, with the data used by Dengler and Dewey (1998) , intensities near 4 could be determined from the percentage of people who felt shaking in each town. With less spatially rich data from self‐selected participants, this percentage often cannot be determined reliably. Audible noises are key additional diagnostic criteria for modified Mercalli intensity (MMI) 4, but, although the DYFI system includes a question about noise, following Dengler and Dewey (1998) , the DYFI algorithm does not include a noise indicator. At the upper end of the scale, as defined the DYFI algorithm yields a maximum intensity value of 9.05, nominally corresponding to peak ground acceleration of 75% g . These and other factors can result in DYFI values that are low compared to traditional MMI values assigned using expert opinion, even absent factors that can bias traditional MMI assignments. Modern ground‐motion intensity conversion equations determined using DYFI intensities are expected to be appropriate for DYFI intensities, but the results of this study suggest that biases may be introduced if DYFI and traditional intensities are assumed to be interchangeable.

Seismological Research Letters↗
USGS2025

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California↗
USGS2026

Remote compositional analyses of space-weathered lunar maria

Visible-to-shortwave infrared (VSWIR) reflectance spectroscopy has revolutionized our understanding of planetary surface compositions. However, space-weathering processes on airless bodies complicate quantitative compositional analyses. Here, we present a framework to isolate the signatures of space weathering in VSWIR spectra of lunar maria by leveraging radiative transfer modeling under the assumptions that (i) a space-weathered target can be expressed as a mixture of fresh and fully space-weathered components and (ii) remaining signatures can be modeled by including agglutinates as an end-member component. We first validate this approach against laboratory spectra of space-weathered Apollo mare soils of known mineral compositions using a probabilistic Markov Chain Monte Carlo implementation of the Hapke radiative transfer model. Second, we illustrate how this approach can be applied to orbital Moon Mineralogy Mapper data. The proposed space-weathering correction workflow for lunar maria could be expanded to other lunar lithologies and applied to existing and future data sets.

Planetary Science Journal↗
USGS2025

Simulating human behavior under earthquake early warning

Earthquakes are a rapid-onset hazard where advance planning and learning plays a key role in mitigating injuries and death to individuals. Recent advances in earthquake detection have resulted in the development of earthquake early warning (EEW) systems. These systems can provide advance warning to predetermined geographic regions that an earthquake is in progress, which may result in individuals receiving warning seconds before significant shaking is felt at their location. This additional time could allow individuals to take more effective protective actions during the immediate disaster. To demonstrate this, we created an agent-based simulation of a basic apartment that allowed us to randomly and repeatedly simulate an individual receiving and responding to an EEW message. The results of our preliminary simulation show that, in our study environment, earthquake early warning alerts have the potential to allow for sufficient time for individuals to take protective actions.

Washington↗
USGS2024

Chloride concentrations in groundwater from the western part of the Southern Hills regional aquifer system, Louisiana, 2021–22

Groundwater is heavily used for public supply and industrial uses in the Baton Rouge, Louisiana, area. Lowered water levels resulting from groundwater withdrawals have induced the movement of saltwater towards wells in East Baton Rouge and West Baton Rouge Parishes. Saltwater intrusion has the potential to affect water supply infrastructure, reduce water availability for some uses, and increase treatment costs. To document current conditions, samples were collected from 161 wells screened in 10 aquifers of the Southern Hills regional aquifer system during November 2021 through February 2022. The results were compared with historical data to identify where chloride concentrations are increasing, which could indicate that saltwater intrusion is occurring. Saltwater intrusion, to varying degrees and areal extents, was observed in most of the 10 aquifers. The limited availability of monitoring wells near or within some of the known saltwater plume areas restricts tracking of the movement or delineation of the plumes’ current extents.

Louisiana↗
USGS2025

Controls on water quality below a reclaimed surface coal mine, southeastern Montana

Coal mining and reclamation can have a profound influence on hydrogeologic systems, with clear consequences for groundwater quality, yet their long-term influence on downgradient water quality over time following reclamation is less well documented. Geochemical trends were evaluated in water quality downgradient of a fully reclaimed landscape at the former Big Sky Mine in the Rosebud Creek watershed (southeastern Montana, USA), over a 3-year period (2020–2022), including bond release in 2022. Within 6 km downgradient from the reclaimed area, sulfate concentrations decreased from approximately 3500 to 1800 mg l −1 within the Miller Coulee alluvial aquifer. Major ions, δ 34 S SO4 values, and residence time tracers suggest that the observed decreases in sulfate concentration result from a combination of dilution by mixed-age inflows and incomplete transit of the high salinity plume from the mine boundary. Both bedrock and alluvial aquifers of the Rosebud Creek corridor contained contributions of millennia-old regional groundwater, which may serve to mitigate mine-derived high salinity waters. Rosebud Creek, which traverses the outflow zone of Miller Coulee in the study area, exhibited high sulfate concentrations during low flows and consistent downgradient increases in sulfate concentration. The possibility of plume dynamics in Miller Coulee suggests that the greatest water quality impacts may not yet have reached Rosebud Creek.

Montana↗
USGS2024

Characterizing directivity in small (M 2.4-5) aftershocks of the Ridgecrest sequence

Directivity, or the focusing of energy along the direction of an earthquake rupture, is a common property of earthquakes of all sizes and can cause increased hazard due to azimuthally dependent ground‐motion amplification. For small earthquakes, the effects of directivity are generally less pronounced due to reduced rupture size, yet the directivity in small events can bias source property estimates and provide important insights into general regional faulting patterns. However, due to observational limitations, directivity is usually only measured and modeled for large events. As such, many studies of small earthquakes either ignore directivity altogether or assume a constant rupture direction for all events in a cluster. In our study, we apply a refined directivity fitting method constrained with two separate methods of source deconvolution to the dataset of aftershocks of the 2019 Ridgecrest earthquakes, which contain a large number of well‐recorded small‐to‐mid sized earthquakes occurring in close proximity to each other. The revealed directivity of 100+ small (M 2.4–5) earthquakes is highly heterogeneous and primarily oblique to and away from the main fault strike, suggesting a complex postseismic stress redistribution. In addition, the energy focusing effect of directivity appears to bias the selection of high‐quality data from stations in the direction of rupture, leading to average stress‐drop increases of 50% if directivity is not accounted for.

California↗