USGS Science⌕ Search

SEARCH · USGS Science

Results for “Remote Sensing Applications: Society and Environment”

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

Forecasting volcanic activity in Germany—A multi-criteria approach

Igneous activity, including shallow intrusions and volcanism, has the potential to disrupt underground critical infrastructure. Notably, future underground infrastructure projects like high-level radioactive waste repositories must be sited in areas of extremely low disruption probability by igneous activity. In Germany, according to the Repository Site Selection Act of 2017 (Standortauswahlgesetz, or StandAG), areas in which Quaternary volcanism is either present or future volcanic activity is expected within the next 1 million years (m.y.) must be excluded from the site selection process. Although the locations of regions with Quaternary volcanism are reasonably well known in Germany, forecasting potential igneous activity at intraplate volcanic fields is challenging, as many processes and their interactions control the spatial distribution of volcanic centers. Here, a semi-quantitative, multi-criteria approach is proposed for a regional evaluation of the relative potential of future igneous activity in Germany. A variety of geoscientific indicators are used, including seismic anomalies in Earth’s mantle, gravity data, tectonic activity, sutures, ground motion, earthquakes, mantle degassing centers, and geochronological data of volcanic rocks. The indicators describe the sequence of processes from potential melt generation in Earth’s mantle, through ascent and accumulation of melt within the lithosphere, to eruption at Earth’s surface. In total, 15 out of 30 proposed geoscientific indicators are selected and quantified using 20 total assigned parameters. Defined threshold values are used to spatially delimit relevant parameter properties to focus on areas with higher potential of future magmatic activity. To consider uncertainties of parameters and their underlying processes, which are usually more spatially extensive below ground, buffer zones are defined in which values of relevance decrease with increasing distance from the initial lateral shape of a parameter. Normalized parameters are combined into an index, whose spatial value distribution is used to differentiate the relative potential of future igneous activity (within the next 1 m.y.). The sensitivity of the results is shown by varying the weighting factors for the relevant parameters in country-wide index maps. Thereby, profiles illustrate the distribution of the resulting index values and respective index fractions of various parameters. Different index maps for the relative potential of future igneous activity are presented and can be used for hazard assessments.

Professional Paper↗
USGS2025

Groundwater surveillance of swine pathogens from private wells supplying swine farms in Iowa

Biosecurity practices are essential for maintaining pig health and productivity. Despite these measures, pathogen spread still occurs. Water is one of the largest daily inputs on swine farms by volume and is not routinely tested or disinfected before it is consumed by the animals [1-3], making it a poorly understood biosecurity risk. Groundwater from privately-owned wells is a common water source for swine farms. Pathogens in the landscape, such as bacteria, viruses, and protozoa, can reach groundwater more rapidly through soil macropores, maintaining viability and facilitating transmission of pathogens into aquifers [3-13].

Iowa↗
USGS2024

Rainfall intensification amplifies exposure of American Southwest to conditions that trigger postfire debris flows

Short-duration, high-intensity rainfall can initiate deadly and destructive debris flows after wildfire. Methods to estimate the conditions that can trigger debris flows exist and guidance to determine how often those thresholds will be exceeded under the present climate are available. However, the limited spatiotemporal resolution of climate models has hampered efforts to characterize how rainfall intensification driven by global warming may affect debris-flow hazards. We use novel, dynamically downscaled (3.75-km), convection-permitting simulations of short-duration (15-min) rainfall to evaluate threshold exceedance for late 21st-century climate scenarios in the American Southwest. We observe significant increases in the frequency and magnitude of exceedances for regions dominated by cool- and warm-season rainfall. We also observe an increased frequency of exceedance in regions where postfire debris flows have not been documented, and communities are unaccustomed to the hazard. Our findings can inform planning efforts to increase resiliency to debris flows under a changing climate.

Arizona, California, Nevada, New Mexico, Utah↗
USGS2025

Rebuilding a volcano one lava flow at a time—Visualizing the lava dome-building eruption in the crater of Mount St. Helens, 1982–1986

Between 1980 and 1986, the U.S. Geological Survey made a series of 1:2,000-scale topographic contour maps from aerial photographic surveys to monitor the eruption. These maps were made for operational purposes and were not intended for publication. Since then, advances in technology made it possible to digitize the original, highly detailed hardcopy maps and derive new digital data elevation models of the surface of the lava dome. These digital elevation models allow for the visualization of the progression of the eruption and reveal the rubbly, chaotic surface of the lava flows and dome. Additionally, these new data help fill gaps in the long-term record of topographic changes that have occurred at the volcano since the May 18, 1980, eruption.

Washington↗
USGS2024

Examining 22 years of ambient seismic wavefield at Mount St. Helens

An increase in seismic activity precedes most volcanic eruptions. Whereas event-based forecasting approaches have been successful, some eruptions remain unanticipated, resulting in casualties and damage. Our study leverages the recent advancements in ambient field seismology. We explore features extracted from continuous ambient fields using traditional methods, for example, peak ground velocity, peak ground acceleration, root mean square, root median square, real-time seismic amplitude measurement, and novel methods (displacement seismic amplitude ratio and spectral width). In addition, we explore unsupervised learning of higher order wavelet features using scattering networks. We find that combining all the methods was necessary to disentangle the effects of seismic sources from structural changes at Mount St. Helens. Although the ambient wavefield-based approach does not yield additional or more significant precursory signals than event-based methods at Mount St. Helens, our study demonstrates that the ambient wavefield provides supplementary information, mainly about structural changes and complements traditional methods. The ambient seismic wavefield offers additional insights into long-lasting processes. We find enhanced wave attenuation correlating with geochemical measurements. We interpret this as ongoing structural changes, such as dome growth or the evolution of the volcanic conduit system. On annual and decadal timescales, we interpret seasonal seismic attenuation in the shallow subsurface as groundwater fluctuations, corroborated by observations at the nearby Spirit Lake level. This multimethod approach at Mount St. Helens sheds light on a volcanic system’s underlying dynamics and structure.

Washington↗
USGS2025

Refinements to the attenuated propagation of local earthquake shaking (APPLES) ground-motion-based earthquake early warning algorithm

We refined the Attenuated ProPagation of Local Earthquake Shaking (APPLES) ground-motion-based earthquake early warning (EEW) approach, and directly compare APPLES performance with that of the source-characterization-based U.S. ShakeAlert EEW system for a suite of historical earthquakes in the U.S. West Coast and Japan. APPLES is an extension of the Propagation of Local Undamped Motion (PLUM) algorithm in which observed shaking intensity at seismic stations is used to forward-predict intensity distributions to surrounding areas using an attenuation model derived from an intensity prediction equation. We test new configuration options within APPLES, such as using the second highest estimated ground motion rather than the maximum, to better match median ground-motion observations and reduce alerts for small magnitude earthquakes, both of which are key alerting priorities within ShakeAlert. We evaluate these configurations alongside ShakeAlert by comparing the ground-motion estimation accuracy and available warning times relative to station observations and ShakeMap distributions. Our preferred APPLES configuration produces accurate ground-motion estimates and corresponds better with median observations compared to ShakeAlert’s estimates. This preferred configuration substantially reduces alert issuance for M < 5.0 earthquakes compared to the previous APPLES configuration, and alert-release criteria can further restrict alerts to primarily M ≥ 5.5 earthquakes without requiring magnitude estimation. Prioritizing matching median-observed ground motions may reduce APPLES warning times compared to configurations that were tuned to avoid missed alerts (such as those that use the maximum estimated ground motions), which can lead to shorter warning times compared to ShakeAlert for the same alert threshold. However, station-based warning time assessments demonstrate that APPLES can outperform ShakeAlert for high target thresholds. APPLES is a simple, independent EEW approach that may improve the robustness of EEW for the West Coast of the U.S.

California, Oregon, Washington↗
USGS2026

Late Miocene Colorado River arrival in the Bidahochi basin supports spillover origin of Grand Canyon

The timing and mechanism of the integration of the Colorado River and incision of the Grand Canyon remain among geology’s enduring controversies. A key question is the configuration of the upper Colorado River watershed between 11 and 6 million years ago. In this study, we present new evidence from zircon uranium-lead geochronology for the arrival of distinctive Colorado–Green River sediment in the Bidahochi basin by 6.6 million years ago derived from the Browns Park Formation. This is coeval with an order-of-magnitude increase in depositional rate, an increase in carbonate strontium isotope ( 87 Sr/ 86 Sr) ratios, the appearance of large fish species characteristic of fast-flowing waters, and other sedimentological changes. This evidence is consistent with the Colorado River supplying water and sediment to the Bidahochi basin before spillover integration of the river through the Grand Canyon.

Arizona, Nevada, Utah↗
USGS2026

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth↗
USGS2025

Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part 2: Xenoliths

Mauna Loa volcano erupts crystal-poor material at its summit and more crystal-rich material on its rift zones. Some of the more olivine-rich lava flows contain xenoliths with diverse mineralogy, including cumulate harzburgites with high-Mg# orthopyroxenes and high-Fo olivines (both > 84). Previous experimental work and thermodynamic modelling has proposed that high-Mg# orthopyroxenes only crystallize from Mauna Loa melts at high pressures (> 6 kbar, > 20 km), leading to suggestions that there is a region of sub-Moho magma storage at Mauna Loa in addition to the geophysically imaged magma reservoir at 2–5 km depth below the summit. We use melt and fluid inclusion barometry combined with thermodynamic models to further investigate this suggestion. Fluid inclusion data from harzburgites and dunitic xenoliths yield storage depths remarkably similar to those found in non-xenolithic crystals from lavas and tephras, with a clear peak at ~ 2–3 km (below the summit). Depths from melt inclusions in these xenoliths overlap with fluid inclusion pressures, ruling out the possibility of fluid inclusion re-equilibration during a period of stalling in a shallower reservoir. We examine five different thermodynamic models and find that the minimum pressure of olivine-orthopyroxene co-saturation varies by ~ 4 kbar (~ 12 km). These models also fail to predict that orthopyroxene is stable in ~ 15–80% of compositionally relevant experimental charges which grew orthopyroxene. Overall, this shows that phase stability modelling is an unreliable method of determining magma storage depth at Mauna Loa. We suggest that model discrepancies reflect a lack of experimental constraints on orthopyroxene stability at > 1200 ℃ and 0.01–5 kbar. Based on the presence of large oikocrystic orthopyroxenes completely enclosing rounded olivine chadacrysts, we suggest that these harzburgitic xenoliths formed through the reaction of intruding melts with olivine mush piles within the Mauna Loa edifice at ~ 3 km depth below the summit, with no need for a deeper storage reservoir. The predominance of pre-eruptive shallow storage means that there is more chance of detecting reservoir destabilization with geophysical monitoring techniques compared to a scenario where melts are supplied from sub-Moho reservoirs.

Hawaii↗
USGS2025

Carbonatite-hosted residual REE deposits

Rare earth elements (REEs) occur in magmatic rocks but are especially enriched in carbonatite and alkaline silicates. If these rocks are chemically weathered, then the REEs may become further enriched within the regolith developed from these rocks. Primary magmatic REE minerals, as well as the various carbonate minerals and apatite, provide the REEs which, under pervasive chemical weathering, are incorporated within low-temperature REE minerals forming within the regolith. Many of these minerals, as well as their textures, are characteristic of this mode of formation. Lateritic conditions of weathering are instrumental in producing a thick, weathered, or regolith, profile, and the roles of sulfide oxidation, fluctuating groundwater tables, and downward mass wasting due to carbonate dissolution are identified as the most important controls on REE enrichment in the regolith.

Book chapter↗
USGS2023

Cursed? Why one does not simply add new data sets to supervised geothermal machine learning models

Recent advances in machine learning (ML) identifying areas favorable to hydrothermal systems indicate that the resolution of feature data remains a subject of necessary improvement before ML can reliably produce better models. Herein, we consider the value of adding new features or replacing other, low-value features with new input features in existing ML pipelines. Our previous work identified stress and seismicity as having less value than the other feature types (i.e., heat flow, distance to faults, and distance to magmatic activity) for the 2008 USGS hydrothermal energy assessment; hence, a fundamental question regards if the addition of new but partially correlated features will improve resulting models for hydrothermal favorability. Therefore, we add new maps for shear strain rate and dilation strain rate to fit logistic regression and XGBoost models, resulting in new 7-feature models that are compared to the old 5-feature models. Because these new features share a degree of correlation with the original relatively uninformative stress and seismicity features, we also consider replacement of the two lower-value features with the two new features, creating new 5-feature models. Adding the new features improves the predictive skill of the new 7-feature model over that of the old 5-feature model; albeit, that improvement is not statistically significant because the new features are correlated with the old features and, consequently, the new features do not present considerable new information. However, the new 5-feature XGBoost model has a statistically significant increase in predictive skill for known positives over the old 5-feature model at p = 0.06. This improved performance is due to the lower-dimensional feature space of the former than that of the latter. In higher-dimensional feature space, relationships between features and the presence or absence of hydrothermal systems are harder to discern (i.e., the 7-feature model likely suffers from the “curse of dimensionality”).

Geothermal Resources Council Transactions↗
USGS2025

Modeling interconnected minerals markets with multicommodity supply curves: Examining the copper-cobalt-nickel system

Demand for many of the metals used in the energy transition is expected to grow rapidly. Many of these are by-products, often considered critical because their production responds weakly to prices and is instead tied to the economics of the host mineral. We present a model of prices and production for jointly produced commodities that accounts for interconnectivity between host and by-product markets at the mine level. We demonstrate this method using the copper–cobalt–nickel system, in which approximately 99% of cobalt is a by-product of copper or nickel mining. Our results show that the model more accurately captures the economic benefits of diversified mine outputs than previous approaches. Furthermore, changes in demand drivers for any two commodities produce non-linear effects on production and price. We challenge the prior best-practice assumption that cobalt cannot impact the copper or nickel markets. Recognizing the importance of both copper and cobalt for future electrification, we emphasize that incentivizing the copper industry to reduce cobalt supply risks could inadvertently undermine copper supply.

Nature Communications↗
USGS2025

Stratigraphic notes—Volume 2, 2025

This is the second 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 .

Professional Paper↗
USGS2026

Stratigraphic notes—Volume 3, 2026

This is the third 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 .

Professional Paper↗
USGS2025

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska↗
USGS2026

Earthquake rupture arrest from depth‐dependent frictional stability on the Pütürge segment of the East Anatolian Fault Zone, Türkiye (Turkey)

Determining why earthquake ruptures stop is a central challenge in earthquake science and seismic hazard assessment. The Pütürge segment of the East Anatolian Fault Zone, Türkiye, exhibits shallow creep (<∼3 km depth) yet hosts M > 6.5 earthquake ruptures at greater depth. Here, we evaluate whether variations in frictional stability along this segment aided arrest of the 2020 M 6.7 Elaziğ and 2023 M 7.8 Pazarcık earthquake ruptures. Analysis of Sentinel‐1 Synthetic Aperture Radar imagery indicates the 2023 M 7.8 rupture propagated laterally into a metamorphic massif within the Pütürge segment, where slip rapidly decayed below detection limits. Creepmeters along this segment recorded no significant surface afterslip (<3 mm) in the following year. To investigate this fault‐slip behavior, we conducted triaxial friction experiments on Pütürge fault gouge sampled from an outcrop exposure. The gouge, composed primarily of muscovite, quartz, and calcite, is velocity strengthening at conditions approximating 0–2.5 km depth and velocity weakening at 4–5 km depth. This transition to velocity‐weakening friction is associated with enhanced comminution and shear localization observed microstructurally. Our results suggest that depth‐dependent frictional stability of the Pütürge fault segment facilitates rupture nucleation and propagation at depth while maintaining rupture‐arresting behavior near Earth’s surface.

Pütürge segment of the East Anatolian Fault Zone↗
USGS2024

Evidence for low effective stress within the crust of the subducted Gorda plate from the 2022 December Mw 6.4 Ferndale earthquake sequence

Stress levels on and adjacent to megathrust faults at seismogenic depths remain a key but difficult to constrain parameter for assessing seismic hazard in subduction zones. Although strong ground motions have been observed to be generated from distinct, high-stress regions on the downdip end of the megathrust rupture areas in many great earthquakes, we lack direct constraints on the stress level in the lower seismogenic portion of the Cascadia megathrust. On 2022 December 20, a M w 6.4 strike-slip earthquake occurred near Ferndale, California in southern Cascadia and likely ruptured the Gorda slab crust in the lower seismogenic portion, providing an opportunity to assess the stress level in this region. Here, we relocate the Ferndale mainshock and the first two weeks of aftershocks using a high-resolution 3-D velocity model and estimate rupture dimensions, directivity, and stress drop for several M w 4-5 aftershocks and recent earthquakes. The aftershocks define a strike-slip fault in the slab crust striking ENE, consistent with the mainshock focal mechanism. The orientation of this fault is about 45 ° off the ideally oriented fault plane given the stress state in the slab. The aftershock zone is extensive and broad in the forward direction of the mainshock rupture but still constrained within the volume of high Vp/Vs within the slab crust. Our stress drop estimates are generally lower for M w 4-5 earthquakes located within the slab crust compared to those a few km deeper in the slab mantle. Combined, our results support a relatively low effective stress level in the vicinity of the megathrust in the lower portion of the seismogenic zone in southern Cascadia, likely due to elevated fluid pressures. Consequently, the ground motion in the onshore region above this low-stress seismogenic portion in southern Cascadia may not be as intense as that observed during great earthquakes in other subduction zones.

California↗
USGS2025

Utilizing downhole datasets for modelling the aeromagnetic signature of the Iron Creek Co-Cu deposit in the Idaho Cobalt Belt

The Idaho Cobalt Belt in east-central Idaho is host to some of the largest domestic Co resources, including the informal Iron Creek deposit. The two main ore zones of this deposit, the Iron Creek and the Ruby, are hosted in greenschist-grade interbedded argillite/siltstone and quartz-rich units of the Mesoproterozoic Apple Creek Formation of Lemhi Group. The primary ore mineral is cobalt-bearing pyrite, which occurs with pyrrhotite, chalcopyrite, and magnetite. This study integrates mineralogical and core-scale geophysical measurements with regional-scale aeromagnetic surveys. The high magnetite content within the Ruby zone produces elevated magnetic susceptibility, but the relatively limited spatial footprint of the ore zone would produce a small-scale anomaly that may be overlooked in regional surveys. The low magnetite content in the Iron Creek zone results in low magnetic susceptibility, creating a relatively low amplitude geophysical response. By characterizing the magnetic properties and mineralogy of these ore zones, this study enhances the interpretation of aeromagnetic data, enabling the identification of small or faint anomalies as potential Co targets. These findings improve can improve exploration strategies, both within the Idaho Cobalt Belt as well as for similar deposit types globally.

Idaho↗