USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Science and Technology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,659 records · Page 83Linked to original sources
USGS2026

Rupture into slow-slip fault regime during the 2018 Mw 6.9 Island of Hawaiʻi earthquake is followed by modest postseismic slip

On 4 May 2018, a M w 6.9 earthquake occurred on the south flank of Kīlauea, in the midst of an historic event that included a voluminous eruption from Kīlauea’s lower East Rift zone and caldera collapse at its summit. The earthquake was a consequence of both short‐ and long‐term stress buildup due to magmatic activity associated with the eruption and steady flank motion, respectively, and it revealed features of Kīlauea’s décollement fault that can inform understanding of future earthquake activity. We used geodetic data to determine the distributions of slip during the coseismic and postseismic periods and compared these with areas of known fault slip during past earthquakes and slow‐slip events (SSEs). The 2018 earthquake ruptured into an area of the décollement fault that was active during quasi‐regular SSEs that occurred in the two decades prior to 2018 but that have not been observed since. The coseismic slip model indicates that the amount of motion on the décollement fault was several times greater than what typically occurred during SSEs, suggesting that it may take decades for the fault to rebuild stress to the point at which SSEs will occur again. Postseismic afterslip also occurred in an area of the fault known to experience slow slip; however, unlike at other creeping faults, postseismic afterslip was rapid, being largely over within 2–3 days. The rapid nature and small magnitude of the postseismic afterslip may be due to the lack of a viscoelastic relaxation component, which is possibly a result of the shallow dip of the décollement fault not transferring stress efficiently into the lower crust.

Hawaii↗
USGS2025

Igneous differentiation in arcs

This chapter assesses crystallization-differentiation of convergent-margin magmas both from the standpoints of natural rock suites and from laboratory phase-equilibrium experiments. Crystallization-differentiation of dry subalkaline and alkaline basalts are summarized at the beginning for reference and context. Experiments crystallizing arc basalts through andesites at middle to upper crustal pressures (100–400 MPa), saturated either with H 2 O or mixed H 2 O-CO 2 , produce liquid lines of descent similar to arc igneous suites, with higher f O 2 s giving results closer to calc-alkaline magma series. Early enrichment of melts in SiO 2 with little or no enrichment in FeO* (total Fe as FeO) comes about by the combined effects of H 2 O suppressing crystallization of plagioclase relative to olivine, shifting plagioclase to more calcic – therefore lower SiO 2 – compositions, and suppressing silicate minerals to greater extents than Fe-rich spinel. Also contributing to SiO 2 enrichment is crystallization near the liquidus of higher proportions of olivine relative to calcic clinopyroxene, compared with MORB, but it is unclear if this is brought about by H 2 O or elevated pressure. Exhumed arc terranes, however, generally fail to expose evidence for primitive liquids having reached and differentiated at middle to upper crustal depths, at least in any abundance. Instead, most differentiation in arcs appears to take place in the lower crust or near the crust-mantle transition, and although hydrous experiments at such pressures produce liquids broadly similar to arc suites in their SiO 2 enrichment, the experimental liquids differ in becoming peraluminous at early stages of differentiation while melts have basaltic, basaltic andesitic, or andesitic SiO 2 concentrations. Crystallization of abundant calcic clinopyroxene close to the liquidus causes the high-pressure experimental liquids to become peraluminous with only small increases in SiO 2 concentration. In contrast, natural arc volcanic and plutonic suites become gradually less metaluminous with increasing SiO 2 , crossing to peraluminous at dacitic to rhyolitic SiO 2 concentrations. Hybridism, encompassing magma-mixing, assimilation of old crust, and incorporation of near-solidus melts from only shortly older deep, antecedent arc intrusions may accompany natural crystallization-differentiation, preventing melts from attaining peraluminous compositions at low-SiO 2 concentrations. If so, then nearly all arc suites are to substantial extents composites, rather than preserving simple, progressive liquid lines of descent. Some comments and assessments are also made on the participation of garnet in arc differentiation and on the production of arc batholiths.

Book chapter↗
USGS2026

Cruise summary—Samoa Basin abyssal mapping—Box coring leg, Pago Pago, Territory of American Samoa to Pago Pago, American Samoa, April 11– May 1, 2026

Expedition Summary U.S. Geological Survey scientists led a box coring effort to the Samoa Basin to characterize minerals and the surrounding abyssal sediments and fauna. Thirty-eight box cores were deployed between April 13, 2026, and April 28, 2026. Thirty-six box cores recovered sufficient material to determine nodule density, and 35 recovered sufficient material for subcores to be collected. The purpose of this Data Report is to provide a summary of samples collected, initial results that were obtained shipboard, and briefly mention planned future analyses from this expedition.

Data Report↗
USGS2025

Footprints of past mining in Alaska (USA) derived from high-resolution satellite imagery

Mapping the land area used for mining in the past is essential for guiding the remediation of affected landscapes and assessing the resource potential of related waste products. Despite significant recent progress delineating footprints of active and inactive mining globally, the known inventory of such mine lands remains incomplete. Here, I describe a new map dataset of footprints of land surface disturbance and waste at sites of past mining in Alaska (USA) based on visual interpretation of satellite imagery. This dataset maps 6–14 times the area of previous regional and global mine footprint maps in Alaska and is the first in the region to explicitly delineate mine waste landforms (e.g., tailings piles). The data are publicly available from the U.S. Geological Survey under a “no rights reserved” Creative Commons (CC0) license agreement.

Alaska↗
USGS2026

USGS Geochron Database

Introduction Geochronology helps us understand Earth’s history by determining when important events, like volcanic eruptions, the rise of mountains, the formation of mineral resources, and changes in the landscape, happened. Geochronological data directly support geologic mapping and can inform decisions about geologic hazard mitigation, natural resource management, and infrastructure resilience. The U.S. Geological Survey (USGS) Geochron database provides access to more than 300,000 published, publicly available age measurements from more than 40,000 geological samples. This database is the result of a collaborative effort with State geological surveys and geoscientists from across the globe. The USGS Geochron database is the most comprehensive collection of geochronological data available for the United States. Users can view data through an interactive map explorer, download datasets, and integrate data into geospatial software or other analysis tools.

Fact Sheet↗
USGS2024

Separating signals in elevation data improves supervised machine learning predictions for hydrothermal favorability

A recent study identified topography (land surface elevation above sea level) as an important input dataset (feature) for predicting the location of hydrothermal systems in the Great Basin in Nevada. Yet, topography is generally a result of more than one geological process and may consequently contain multiple distinct signals. For example, the geologic evolution of the Great Basin has produced both crustal thickening (i.e., regional-scale trends in elevation) and thinning via Basin and Range extensional faulting (i.e., valley-scale topographic relief). We postulate that these geologic processes may affect the occurrence of hydrothermal systems differently. Therefore, we separate the regional trend from the valley-scale signal in the Great Basin, and then use them separately to evaluate the importance of each as predictors for hydrothermal favorability. Our prior work applying supervised machine learning (ML) using the data from the Nevada Machine Learning Project demonstrated that employing a training strategy that randomly selects negative training sites produces better performing models for predicting hydrothermal favorability than a training strategy that uses expert-selected negatives. The models created using both training strategies exhibited a west-east geographic trend in the predictions for the favorability of hydrothermal resources. These models generally predicted higher favorability in western Nevada and lower favorability in eastern Nevada. This west-east trend in predicted favorability correlates with elevation across the Great Basin, which trends higher from west to east. By separating the original elevation feature into distinct features for elevation trend (i.e., regional-scale topography) and detrended elevation (i.e., valley-scale or local relative topography), we find that models using the separated topographic signals consistently outperform competing models that use the original elevation feature. Although western Nevada still exhibits higher favorability than eastern Nevada, using separated signals for regional elevation and local structure reduces the west-east prediction trend in the region and emphasizes structures associated with hydrothermal upflow. This work emphasizes how carefully engineering features to represent geological conditions relevant to hydrothermal systems allows ML algorithms to detect important patterns for predicting hydrothermal resource favorability and leads to better model performance.

Conference Paper↗
USGS2025

Long‐period ground motions from dynamic rupture simulations of large earthquakes on the creeping Hayward–Calaveras–Rodgers Creek fault system

he Hayward, Calaveras, and Rodgers Creek faults in the San Francisco Bay region of California have a high probability of producing a large earthquake in the next decades. Although these faults creep, the creep is insufficient to keep up with their relatively rapid slip rates on their deepest sections, so they have been storing tectonic strain since their last large earthquakes, with the Hayward’s and Rodgers Creek’s more than 150 yr ago. We do not know what the next large Hayward–Calaveras–Rodgers Creek earthquakes will look like or how strongly they will shake the San Francisco Bay region. Harris et al. (2021) used the 3D dynamic (spontaneous) rupture method to simulate large earthquakes on these creeping faults. In this article, we examine the resulting simulated long‐period ( T > 1 s) ground shaking from 0 to 50 km distance, for earthquakes nucleating on the Hayward fault and earthquakes nucleating on the Rodgers Creek fault. We compare these simulated long‐period ground motions with the Boore et al. (2014) well‐established empirically based ground‐motion model suitable for the slowest material velocity in our 3D velocity structure. We find that the simulated long‐period ground motions from the creeping‐fault earthquake scenarios produce a reasonable agreement with the empirical expectations if frictional cohesion is included only where it is appropriate.

California↗
USGS2025

Site-specific, extended ShakeMaps for earthquake engineering applications

The U.S. Geological Survey (USGS) routinely produces ShakeMaps of shaking intensity across the globe. Due to practical constraints, the number of response spectral periods was limited to three standard periods (0.3, 1.0, and 3.0 sec). We have recently developed the tools that are necessary to expand this functionality to include 22 periods (matching the current U.S. National Seismic Hazard Model periods) as well as the orientation-independent components (e.g., “RotD50”). We refer to ShakeMap products that include these extensions as “extended ShakeMaps.” The added level of complexity motivated us to also develop a user-friendly tool called the “ShakeMap Sampling Tool” (SST) that gives all the estimated shaking metrics for a specific location (or list of locations). Additionally, we develop a web application where users can input locations of interest and view/download the SST results. We further familiarize users with the concept of “Composite ShakeMaps.” For earthquakes sequences such as a mainshock and larger foreshocks and aftershocks, this provides a map of the maximum value of each shaking metric, which is useful for overall loss estimates, the full extent of ground failure triggering potential, and a better portrayal of the repeated shaking levels at a given point for a series of earthquakes. Such a site-specific shaking history facilitates earthquake forensics at building or infrastructure sites for which damage may be of concern, as described in the Disproportionate Damage Earthquake trigger specified in the IEBC (2018, Section 405.2.2) and in developing ATC-145 guidelines (Guidelines for Post-Earthquake Assessment, Repair, and Retrofit of Buildings). The composite ShakeMap can be combined with the SST for a variety of earthquake-hazard applications, such as systematically inferring triggering shaking estimates at specific sites of geotechnical interest for landsliding, liquefaction, and lateral-spreading hazards.

Kahramanmaraş↗
USGS2026

Measuring stress In high pressure deformation experiments with high speed fiber-optics

High‐pressure, high‐temperature rock deformation experiments are essential for understanding deep Earth processes, but accurately measuring stress is challenging due to the inaccuracy introduced by seal friction within the apparatus and large inertia, which affects earthquake process measurements made far from the experimental fault. To overcome these limitations, we developed and implemented a simple, compact fiber‐optic sensor based on an External Cavity Fabry–Perot Interferometer for in situ load measurement on a piston inside the pressure vessel of a Griggs‐type high‐pressure apparatus. The sensor can be used at sample temperatures up to 800°C for both slow creep and fast rupture testing with bandwidth capability from DC to 6 MHz. Two important experimental results are described: (1) the first direct measurements of seal friction in this apparatus type; these measurements record seal friction approximately one‐tenth of the confining pressure under both low (150 MPa) and high (1 GPa) pressure conditions; and (2) high‐temporal resolution capture of dynamic stress drops during stick‐slip faulting events at high‐confining pressure, revealing near‐fault wave propagation details, high‐frequency oscillations, and implying high‐transient slip rates (6–12 m/s) previously inaccessible with external sensors.

The Seismic Record↗
USGS2024

Cross-fade sampling: Extremely efficient Bayesian inversion for a variety of geophysical problems

This paper introduces cross-fade sampling, a computationally efficient Markov Chain Monte Carlo simulation method that uses a semi-analytical approach to quickly solve Bayesian inverse problems that do not themselves have an analytical solution. Cross-fading is efficient in two ways. First, it requires fewer samples to obtain the same quality simulation of the target probability density function (PDF). Secondly, it is much faster to evaluate the posterior probability of each sample than conventional sampling methods for simulating Bayesian posterior PDFs. Conventional methods require evaluating the prior probability (which describes your a priori constraints) and data likelihood (which describes the fit between the observations and the predictions of the model) for each sample model. However, cross-fading does not require evaluating the data likelihood, meaning that ‘big data’ can be fit with zero additional computational cost. Further, the cross-fading approach can be used to calculate the marginal likelihood associated with a model design, facilitating model comparison and Bayesian model averaging. Topics covered in this paper include derivation of the cross-fade approach and how it can be used to simulate Bayesian posterior PDFs and compute the marginal likelihood, discussion of the class of problems to which cross-fading can be applied (with examples from earthquake statistics, earthquake ground motion modelling, volcanic eruption forecasting, and finite fault slip modelling), demonstration of efficiency relative to existing sampling methods and discussion of how cross-fading can be used to account for prediction errors (i.e. epistemic errors) as part of the geophysical inverse problem.

Geophysical Journal International↗
USGS2025

Application of Hydrologic Simulation Program—FORTRAN (HSPF) as part of an integrated hydrologic model for the Salinas Valley, California

The U.S. Geological Survey (USGS), in cooperation with the Monterey County Water Resources Agency, conducted studies to help evaluate the surface-water and groundwater resources of the Salinas Valley study area, consisting of the entire Salinas River watershed and several smaller, adjacent coastal watersheds draining into Monterey Bay. The Salinas Valley study area is a highly productive agricultural region that depends on the coordinated use of surface water and groundwater to meet demand for irrigation and public water supply. To continue to meet these demands, a better understanding of the historical water balance and the effects of water-resource development on the long-term sustainability of water resources in the Salinas Valley study area is needed.

California↗
USGS2024

Limited preservation of strike-slip surface displacement in the geomorphic record

Offset geomorphic markers are commonly used to interpret slip history of strike-slip faults and have played an important role in forming earthquake recurrence models. These data sets are typically analyzed using cumulative probability methods to interpret average amounts of slip in past earthquakes. However, interpretation of the geomorphic record to infer surface slip history is complicated by slip variability, measurement uncertainty, and modification of offset features in the landscape. To investigate how well geomorphic data record surface slip, we use offset measurements from recent strike-slip surface ruptures ( n = 39), faults with geomorphic evidence of multiple strike-slip earthquakes ( n = 29), and synthetic slip distributions with added noise ( n > 10,000) to examine the constraints of the geomorphic record and the underlying assumptions of the cumulative offset probability distribution analysis method. We find that the geomorphic record is unlikely to resolve more than two paleo-slip distributions, except in specific cases with low slip variability, high slip-per-event, and semiarid climate. In cases where site-specific conditions allow for interpretation of more than two earthquakes, lateral extrapolation along a fault is not straightforward because on-fault displacement and distributed deformation may be spatially variable in each earthquake. We also find that average slip in modern earthquakes is adequately recovered by probability methods, but the reported prevalence of strike-slip faults with characteristic slip history is not supported by geomorphic data. We also propose updated methods to interpret slip history and construct uncertainty bounds for paleo-slip distributions.

Journal of Geophysical Research: Solid Earth↗
USGS2025

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗
USGS1990

Preconditioned Conjugate-Gradient 2 (PCG2), a computer program for solving ground-water flow equations

This report documents PCG2: a numerical code to be used with the U.S. Geological Survey modular three-dimensional, finite-difference, ground-water flow model. PCG2 uses the preconditioned conjugate-gradient method to solve the equations produced by the model for hydraulic head. Linear or nonlinear flow conditions may be simulated. PCG2 includes two reconditioning options: modified incomplete Cholesky preconditioning, which is efficient on scalar computers; and polynomial preconditioning, which requires less computer storage and, with modifications that depend on the computer used, is most efficient on vector computers. Convergence of the solver is determined using both head-change and residual criteria. Nonlinear problems are solved using Picard iterations. This documentation provides a description of the preconditioned conjugate gradient method and the two preconditioners, detailed instructions for linking PCG2 to the modular model, sample data inputs, a brief description of PCG2, and a FORTRAN listing.

Water-Resources Investigations Report↗
USGS2026

Rapid characterization of the 2025 Mw 8.8 Kamchatka, Russia earthquake

The 29 July 2025 M w 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near‐source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real‐time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite‐fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The M w 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

Kuril-Kamchatka subduction zone↗
USGS2026

Incorporating location uncertainty improves inference with stop-level North American Breeding Bird Survey data

Ecological models should account for uncertainty to be most effective and useful. Yet, uncertainty from model covariates—unlike that from other sources, such as sampling error or process variability—is seldom explicitly incorporated. This can cause underestimates of uncertainty to cascade through model parameter estimates, predictions, and downstream uses. Burner et al. proposed a method for quantifying uncertainty in covariates and incorporating it into models using informative Bayesian priors. This method was applied to stop-level Breeding Bird Survey (BBS) analyses, where land cover uncertainty at each stop arises from substantial stop location uncertainty. A limited validation of model-estimated land cover, using stops with known locations, indicated the method’s potential effectiveness, but it was not rigorously evaluated. We conduct a robust simulation-based test, generating stop locations, extracting land cover, and simulating bird communities across 210 BBS routes in the upper Midwest. We compare 3 models: a “known” model with true land cover, a “naive” model assuming consistent 800-m stop spacing, and a “full” model using informative priors to estimate land cover. Species parameter estimates and predicted prevalence patterns across gradients in land cover from the full model approached those of the known model and were substantially closer to the true values used in simulations relative to those from the naive model. Naive model parameters were more biased relative to the other models, and credible intervals of predicted species prevalence rarely included the true simulated values. The full model also produced land cover covariate estimates closer to true simulation values relative to the mean informative priors. Our results show that, for the BBS, informative priors enable more accurate stop-level analyses despite location uncertainty. In contrast, naive models that ignore this uncertainty yield poor inferences. More broadly, we demonstrate empirically the utility of informative priors to account for covariate uncertainty in ecological models.

Michigan, Minnesota, Wisconson↗
USGS2023

Stratigraphic notes—Volume 1, 2022

This is the first volume in the U.S. Geological Survey (USGS) series of reports on stratigraphy entitled “Stratigraphic Notes,” which consists of short papers that highlight stratigraphic studies, changes in stratigraphic nomenclature, and explanations of stratigraphic names and concepts used on published geologic maps. “Stratigraphic Notes” is a long-term (multiyear), multivolume publication containing articles that address updates or revisions to stratigraphic nomenclature (and whose content ultimately will be incorporated by National Geologic Map Database personnel into Geolex, https://ngmdb.usgs.gov/Geolex/ ). We welcome papers for the “Stratigraphic Notes” series from geoscientists of the USGS, of State Geological Surveys, and from academicians. Papers can be submitted for publication in “Stratigraphic Notes” by contacting the USGS Geologic Names Committee ( gnc@usgs.gov ). As new “Stratigraphic Notes” volumes are published, links to the volumes will be posted at https://doi.org/10.3133/pp1879 . This first volume ("Stratigraphic notes—Volume 1, 2022") includes articles that provide guidance for those who wish to submit papers to “Stratigraphic Notes,” as well as information on how to make your manuscripts compliant for geologic names reviews and how to organize your paper’s content to facilitate inclusion of new or revised names in Geolex. This volume also includes some specific guidance on conducting geologic names reviews of geologic and hydrogeologic reports.

Professional Paper↗
USGS2025

Groundwater-storage change in the north Phoenix aquifer, Arizona, 2020–23

The city of Phoenix, Arizona, relies primarily on surface water for municipal water supply. The city also maintains wells to withdraw groundwater, particularly in times of drought and reduced surface-water supply, and to recharge groundwater when excess surface water is available. As of 2023, withdrawals from the aquifer in the northeastern part of the city are a small volume of water, less than 3,000 acre-feet in most years. Each year a similar volume of water is recharged through injection wells. The withdrawal wells are permitted to produce in total more than 20,000 acre-feet per year; increased pumping could lead to future groundwater-storage declines. To better understand groundwater-storage change in the north Phoenix aquifer, a repeat microgravity monitoring network was established in 2020. Measurements of changes in Earth’s gravity provide a direct, non-invasive measurement of subsurface mass change. Groundwater-storage changes were small during the 2020–23 study period, consistent with the relatively small volumes of pumping and recharge. Groundwater levels measured in monitoring wells were stable or increased slightly during this period, although the number of monitoring wells within the area of the gravity monitoring network is sparse. In total, about 15,000 acre-feet of water were pumped and 31,000 acre-feet recharged through injection in the north Phoenix aquifer during the 2020–22 period. Within the monitored area, groundwater storage increased by a small amount, about 1,000 acre-feet, in 2020, and decreased by a larger amount, about 6,000 acre-feet, each year in 2021 and 2022. Storage decreased at 89 of 102 stations from 2020 to 2023. Groundwater-storage decreases greater than the volume of net pumping indicate down-gradient subflow from the aquifer is greater than recharge plus incoming subflow, drying of the unsaturated zone resulting from decreased land-surface recharge, or both. At present (2023), the aquifer appears able to store and supply the relatively small amounts of water needed without excessive drawdown or harmful effects, such as land subsidence and surface inundation. If pumping or recharge increases in the future, the established repeat microgravity network is well suited to capture the spatial extent and magnitude of groundwater-storage changes in the aquifer.

Arizona↗