USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Surface Processes and Landforms”

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,665 records · Page 90Linked to original sources
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↗
USGS2025

Searching for seismic precursors - The Barry Landslide hazard clean up

The Barry Landslide, located in Barry Arm of Prince William Sound, Alaska, poses a major hazard due to its steep, unstable slopes and the potential for a massive landslide-generated tsunami. With an estimated volume of 500–700 million cubic meters, the Barry Landslide could trigger highly destructive waves. In this study, we focus on seismic signals from the Barry Landslide, which are critical for providing timely tsunami warnings. Since the summer of 2020, the region has been instrumented to monitor the landslide, but the seismic record is complicated by the presence of nearby glaciers and frequent regional earthquakes. Among these signals, we analyze a specific class of short-duration, high-frequency seismic events that exhibit strong seasonal variability, increasing in rate from late summer to mid-winter before ceasing abruptly in late winter or early spring. Our analysis suggests that the source of these signals is likely near or beneath Cascade Glacier, adjacent to the landslide, rather than within the landslide mass itself. We apply detection algorithms to construct a time history for this signal type, which we then compare with environmental factors like precipitation, temperature, and slope displacement data from ground-based radar and remote sensing. Correlations indicate that these seismic events may be driven by seasonal hydrological changes, particularly the freeze-up of subglacial water pathways. While these events are not directly linked to landslide motion, they serve as indirect markers of subsurface hydrological conditions that influence slope stability. Our findings highlight the complex interplay between glaciers, groundwater, and landslide dynamics, emphasizing the need for multi-parameter monitoring to assess evolving geohazards in the region.

Alaska↗
USGS2026

Methodology for quantitative spatial sensitivity analysis of volcanic geodetic networks

Introduction This report introduces a methodology for assessing the state of the U.S. Geological Survey Volcano Observatories’ geodetic monitoring networks that measure how volcanoes deform or change shape. This new method uses a model-based approach that considers the uniqueness of the instrument environments at each volcano. This report focuses on simplified volcanic sources, is independent of the shape or size of the volcano, or the network geometry, and thus highlights the strengths and potential vulnerabilities of each volcano’s geodetic network in an actionable visual format. This analysis can help observatories to make informed decisions about whether volcanoes have an adequate level of geodetic monitoring and indicate where improvements are needed.

Lassen Peak, Mount Shasta↗
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↗
USGS2026

Roadway runoff induced acute mortality in juvenile coho salmon during spring storm events

Extensive mortalities of adult coho salmon ( Oncorhynchus kisutch ), often called “Urban Runoff Mortality Syndrome” (URMS), have been documented during the fall in creeks where water quality has been degraded by roadway runoff. The primary cause of mortality is 6PPD-quinone (6PPDQ; N -(1,3-dimethylbutyl)- N ′-phenyl- p -phenylenediamine quinone)–an ozone transformation product that forms on all vehicle tires. Laboratory studies have shown that juvenile coho salmon are highly sensitive to 6PPDQ exposure. Unlike adults, juveniles reside in impacted watersheds year-round, including during the spring when 6PPDQ concentrations can frequently exceed lethal thresholds during storms. To assess the potential incidence of URMS in springtime rearing habitats for juvenile coho salmon, we conducted a paired water quality and toxicology study at Miller Creek, a runoff-impacted watershed in Normandy Park, WA, USA. Using a small field facility, three naïve groups of juvenile coho salmon ( N = 720) were exposed to either creek water or groundwater ( N = 120 per treatment per storm), across three spring storms while comparing water quality and mortality end points. In creek water during exposures, peak 6PPDQ concentrations reached 73–110 ng/L, exceeding reported median lethal concentrations (LC 50 ) for coho salmon. Over each 24–73 h storm exposure period, ∼80% of Miller Creek-exposed juvenile salmon died. No mortality occurred among control fish exposed to groundwater. These results indicate previously unidentified mortality risks for juvenile life stages of coho salmon during spring storms, suggesting substantial and year-round water quality impediments to coho salmon health and recovery across roadway runoff-impacted spawning, rearing, and migratory habitats.

Washington↗
USGS2025

Preventing overfitting when using tree-based methods for mapping hydrothermal favorability

Ensemble tree-based algorithms are robust tools for estimating sparsely distributed resources with non-linear dependencies (e.g., hydrothermal systems). These algorithms naturally accommodate the threshold conditions necessary to enable and support hydrothermal systems (e.g., having sufficient heat and permeability) and are simpler than many other non-linear machine learning strategies (e.g., artificial neural networks), which is an advantage when working with few labeled examples from which to learn. In previous work, we used eXtreme Gradient Boosting (XGBoost) to produce regional prediction and uncertainty maps of hydrothermal favorability; however, recent studies suggest that, even when properly applied, XGBoost has some risk of overfitting when there are few labeled examples from which to learn. To evaluate overfitting when constructing hydrothermal favorability maps with tree-based methods, we compare XGBoost with Extremely Randomized Trees (ExtraTrees), another ensemble tree-based algorithm that has the potential to underfit when using few labeled examples. We hold all other modeling parameters constant, resulting in two contrasting favorability maps of conventional geothermal resources for the Great Basin. Our results indicate that ExtraTrees demonstrably reduces overfitting compared with XGBoost. After considering overall performance, we conclude that ExtraTrees provides a more suitable modeling approach than XGBoost for the purposes of conventional hydrothermal resource assessments.

Conference Paper↗
USGS2026

Biome-scale spatial patterns of avian abundance reveal proactive conservation opportunities in North American grasslands

North American grassland birds have experienced steeper population declines than any other avian guild, yet conservation efforts remain largely reactive and fragmented. We used nearly four decades of North American Breeding Bird Survey data to identify biome-scale spatial patterns (clustering) of grassland bird abundance for the Great Plains. Our results reveal an ecological core in the north-central Plains where community-level abundance is either increasing by >100% or remains high and stable, providing a strategic roadmap for a “Defend the Core” conservation approach. This approach flips the script from reactive triage centered on isolated population fragments to a proactive strategy of maintaining large-scale ecosystem integrity. Conversely, we found that population losses are more spatially clustered than wins, reflecting the relentless, one-way movement of woody encroachment and agricultural conversion. This asymmetry supports prioritizing intact landscapes, as current restoration rates are often outpaced by the scale of habitat loss. Notably, we found that community-level spatial clustering is a more robust indicator of biome condition than trends of individual flagship species, suggesting that managing for ecosystem integrity provides a more effective multi-species umbrella. Given our results, there is an opportunity for operationalizing a Great Plains Conservation Design that is ecosystem-centric and rooted in the sustainability of the private-land cattle production that maintains these open spaces. By leveraging avian abundance as a biological sensor, managers and producers can deploy a shared vision that matches the spatial scale of the threats, moving from reactive triage to proactive defense of core working grasslands in North America.

Great Plains biome↗
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

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↗
USGS2026

Testing characteristic magnitude distributions in modern PSHA models

The characteristic magnitude distribution hypothesis predicts a higher rate of large earthquakes than a Gutenberg–Richter extrapolation of the small‐earthquake rate would imply. Characteristic magnitude distributions have been commonly applied to faults in probabilistic seismic hazard analysis (PSHA), and in modern models they can emerge from the way short‐term seismicity constraints are combined with long‐term geologic and geodetic constraints. We test the characteristic magnitude distribution hypothesis by comparing the fault‐based magnitude distributions from the 2023 update to the National Seismic Hazard Model (NSHM23) in the Western United States with observed seismicity over the past 93 yr. We find that observed magnitude distributions fall outside the model‐predicted confidence bounds in regions where NSHM23 produces characteristic magnitude distributions: in these regions, the model predicts higher rates of large earthquakes than are observed. An analysis of the earlier California model (Uniform California Earthquake Rupture Forecast, version 3) also reveals discrepancies between the modeled and observed magnitude distributions. In addition, we find that observed magnitude distributions near modeled faults are not significantly different from those in background regions. These results challenge the prevalence of characteristic magnitude distributions in fault‐based seismic hazard models and call for a reassessment of how disparate data sets are integrated in PSHA.

western United States↗
USGS2026

Preliminary geology of the North Meadow Creek Area, Tobacco Root Mountains, southwest Montana: The North Meadow Creek Fault

This report documents a previously unmapped fault, informally referred to herein as the North Meadow Creek fault, on the east flank of the Tobacco Mountains of southwestern Montana. This fault has an apparently long and complex history, including Quaternary offset of an older alluvial terrace, offset of a debris flow of presumed Pleistocene age, offset of rhyolite lava flows of presumed Tertiary age from their presumed source several kilometers east, on the opposite side of the fault, and offset and possible rotation of blocks of Archean gneiss on northeast side of the fault. The western part of the fault is covered by glacial moraines of two apparent ages, both Pleistocene.

Montana↗
USGS1996

Techniques to estimate generalized skew coefficients of annual peak streamflow for natural basins in Texas

This report presents two techniques to estimate generalized skew coefficients used for log-Pearson Type III peak-streamflow frequency analysis of natural basins in Texas. A natural basin has less than 10 percent impervious cover, and less than 10 percent of its drainage area is controlled by reservoirs. The estimation of generalized skew coefficients is based on annual peak and historical peak streamflow for all U.S. Geological Survey streamflow-gaging stations having at least 20 years of annual peak-streamflow record from natural basins in Texas. Station skew coefficients calculated for each of 255 Texas stations were used to estimate generalized skew coefficients for Texas. One technique to estimate generalized skew coefficients involved the use of regression equations developed for each of eight regions in Texas, and the other involved development of a statewide map of generalized skew coefficients. The weighted mean of the weighted mean standard errors of the regression equations for the eight regions is 0.36 log 10 skew units, and the weighted mean standard error of the map is 0.35 log 10 skew units. The technique based on the map is preferred for estimating generalized skew coefficients because of its smooth transition from one region of the State to another.

Texas↗
USGS2026

Estimating aftershock risk for entry into earthquake-damaged buildings

We present a simple method to estimate the risk of experiencing strong shaking from aftershocks during entry into earthquake-damaged buildings. We compute wait times until the probability of strong ground shaking from aftershocks reaches a predefined risk threshold; for example, a 0.4 percent probability of experiencing Modified Mercalli Intensity 7 or greater shaking during the planned building entry. We also develop a relation between aftershock probability and the probability of strong shaking, so that users can reference the U.S. Geological Survey aftershock forecast during an ongoing aftershock sequence to determine if the risk threshold has been met. We apply our method to active continental regions (for example, the Western United States), stable continental regions (for example, the Central and Eastern United States), and subduction zones (for example, Cascadia or Alaska).

Open-File Report↗
USGS1976

Uranium analysis of single drops of natural waters using the fission-track technique

The fission-track technique has been applied to uranium analysis of natural waters using a single drop of water. Samples may be collected in the field by evaporating a drop of water on a detector, and samples may then be analyzed by returning the detector to the laboratory. A sample collecting kit has been assembled, which facilitates the designed sampling method. Some unique advantages of this technique are picogram sensitivity (10 -12 g) for uranium on single drops of water, maintenance of sample integrity during transport and storage, and the ability to discern whether the uranium is in solution or in the suspended material. Counting the tracks is the most time-consuming part of the analysis, but it can be reduced by using automatic counting systems. The analytical error is generally ± 10 percent over the entire range of concentrations.

Open-File Report↗
USGS2026

Microtextural characteristics of adularia in banded quartz veins from the Midas low-sulfidation epithermal deposit, Nevada

High-grade ores at the Miocene Midas low-sulfidation epithermal deposit in northern Nevada are confined to crustiform quartz veins containing abundant adularia. Micro-X-ray fluorescence elemental mapping reveals that adularia is a common gangue mineral occurring in colloform bands, bands showing bladed textures, and bands with dendritic terminations. The adularia aggregates have delicate shapes and are comprised of stacked, submillimeter crystals hosted by fine-grained quartz. The textural evidence suggests that the adularia aggregates originally formed within a gel-like, noncrystalline silica matrix, which subsequently transformed into quartz. This indicates that the adularia did not precipitate in open space along the vein walls. Correlative microscopy, involving scanning electron microscopy-based automated mineralogy and optical petrography, demonstrates that bands containing abundant adularia are not the primary host to ore minerals. The ore minerals occur in different bands within the crustiform veins, implying that adularia and ore mineral precipitation did not always occur simultaneously. It is hypothesized here that fluid flow at Midas involved intermittent short-lived events of fluid flashing, causing rapid solute supersaturation in the liquid. During each flashing event, different amounts of vapor were produced along a given vein. Compositional differences between adjacent bands in the crustiform quartz veins may, therefore, be linked to variations in the amount of vapor formed during each flash event.

Nevada↗
USGS2025

Analysis of the potential effects of Uzbekistan’s mineral endowment on the critical mineral supply of tungsten

Tungsten appears on the 2018 and 2022 U.S. Geological Survey critical mineral lists in part because of a very high global production concentration in China, which produces almost 83 percent of the world’s mined tungsten. Using known parameters and values from other tungsten mining operations, we created hypothetical scenarios in which three tungsten deposits in Uzbekistan are considered for development. Our results show that all three deposits are likely to be economically viable to develop under 2024 market conditions. If the three studied tungsten deposits were put into production, Uzbekistan could become the third-leading tungsten-producing country in the world and increase world output of tungsten by 2.7 percent. Putting these tungsten deposits in Uzbekistan into production could slightly reduce the tungsten global market concentration, therefore reducing the supply disruption potential for tungsten.

Open-File Report↗
USGS2024

Simulated ground-motion records for the seismic assessment of monumental masonry structures

Earthquakes are natural disasters that can cause widespread devastation and loss of life. Simulated ground-motion records can be useful in regions with limited seismic stations or a history of damaging but infrequent earthquakes. This is especially true in areas with a high concentration of heritage masonry structures, which are especially susceptible to damage, as simulated records can be crucial in predicting their seismic response and protecting these buildings from seismic damage. Despite the importance of simulated earthquakes, few studies have investigated how effective they are compared to real earthquakes when assessing the structural response of heritage buildings. To address this knowledge gap, we employ two different simulation methods of the Mw 6.2 Faial earthquake, which occurred on July 9th, 1998, in the Azores, to replicate the recorded time-series at four available stations within an epicentral distance of 150 km. The study has two objectives: first, to validate the simulated records of the 1998 Faial earthquake using alternative stochastic ground-motion simulation approaches, and second, to determine how these approaches impact the seismic assessment of historic masonry structures. To accomplish these objectives, this study uses real and simulated ground-motion datasets to conduct non-linear response history analyses of the São Francisco Church, a monumental structure in Horta that sustained damage during the Faial earthquake. The results show that both simulation approaches yield structural responses similar to the observed records.

Conference Paper↗