USGS Science⌕ Search

SEARCH · USGS Science

Results for “Applied Computing and Geosciences”

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

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

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

Communications Earth & Environment↗
USGS2026

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

Re-examining the growing degree day minimum for grass carp spawning initiation in Lake Erie: Have we missed the mark?

Understanding the timing of life history processes of invasive fish species, such as adult aggregations during spawning, allows natural resource managers to optimize targeted population control efforts. Temperature-based estimates of the timing of spawning initiation for grass carp ( Ctenopharyngodon idella ) exist from its native range, but have not previously been developed for grass carp reproduction in its invaded range in the Laurentian Great Lakes of North America. We accounted for differences among methods across studies and calculated cumulative growing degree days (GDD) for the dates of grass carp egg capture in Lake Erie tributaries (Sandusky, Maumee, and Huron rivers) during 2015–2024 and compared these values to the existing GDD base 15°C requirement estimate from the native range (633 GDD–15°C). The earliest egg detection in Lake Erie (Sandusky River) was 19 GDD–15°C or 361 GDD–5°C (latest egg detection was 2285 GDD–5°C from the Sandusky River), indicating that spawning was initiated up to 25 days earlier than would be predicted by the native range GDD. Moreover, spawning in Lake Erie tributaries occurred earlier than predicted for 58% of spawning events examined based on the native range GDD. The GDD values for spawning initiation in Lake Erie tributaries can be used to refine the timing of control efforts targeting pre-spawn and spawning adult fish, inform the timing of egg sampling, and provide linkages to existing early life stage models to infer the length of the growing season available for age-0 fish to reassess tributary suitability.

Lake Erie, Maumee, River, Sandusky river↗
USGS2026

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

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

Washington↗
USGS2025

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

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

California↗
USGS2025

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

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

Montana↗
USGS2026

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

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

Physical Review Letters↗
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↗
USGS2025

Rare earth elements on the Moon

Rare earth elements (REEs) are a scarce but vital resource for our modern economies and lifestyles. Since the late 1990s, China has supplied the vast majority of the world’s refined REEs. Increasing global demand has broadened the search for REE deposits to unconventional places, including the Moon. Although most lunar rocks have very low REE concentrations, Apollo samples showed that one type of lunar rock containing potassium (K), REEs, and phosphorus (P)—known by the acronym KREEP—has high concentrations of REEs. Data from orbiting satellites have identified locations where substantial deposits of KREEP are likely. The viability of mining these deposits depends on the evolution of REE economics, the development of the Earth-Moon infrastructure, and the findings from future lunar mineral exploration missions.

Fact Sheet↗
USGS2024

Geoelectric evidence for a wide spatial footprint of active extension in central Colorado

Three-dimensional magnetotelluric (MT) imaging in central Colorado reveals a set of north-striking high-conductivity tracks at lower-crustal (50–20 km) depths, with conductive finger-like structures rising off these tracks into the middle crust (20–5 km depth). We interpret these features to represent saline aqueous fluids and partial melt that are products of active extensional tectonomagmatism. These conductors are distributed over a wider region than the narrow corridor along which Rio Grande rift structures are traditionally mapped at the surface, and they consequently demarcate regions of the lower crust where accommodation of bulk extensional strain has concentrated conductive phases. Our observations reveal limitations in existing models of Rio Grande rift activity and may reflect unrecognized spatiotemporal variations in rift system evolution globally.

Colorado↗
USGS2025

State shifts in the deep Critical Zone drive landscape evolution in volcanic terrains

Understanding the near-surface environment where atmospheric and solid earth processes interact, often termed the “Critical Zone,” is important for assessing resources and building resilient societies. Here, we examine a volcanic landscape in the Oregon Cascade Range, an understudied Critical Zone setting that is host to major regional water resources, pervasive silicate weathering, and significant geohazards. We leverage a bedrock age chronosequence to show that the volcanic Critical Zone undergoes a structural shift, from depth extents of >1 km to meters, over timescales of ~1 My. We map an active groundwater volume comparable to major continental lakes, stored at the Cascade Range crest. This state shift makes volcanic landscape evolution a unique probe of deep coupling between Earth systems.

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

A ground-motion model derived using a generalized mean rupture distance for large slab interface earthquakes

Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup ⁠ . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.

Seismological Research Letters↗
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↗
USGS2025

Anomalous shear stress variation in wet granular medium: Implications for landslide lateral faults

Landslide assessments typically focus on the mechanical properties of the basal shear zone, but lateral faults are frequently overlooked, possibly due to their lower normal stresses and variably saturated conditions. Using double-cylinder shear experiments on wet granular systems as analogs for landslide lateral faults, we observe anomalous shear stress variations with fluid volume fractions, defying an expected unimodal relationship associated with capillary cohesion. At low fluid volume fractions, shear strength weakens as the wet grain assembly experiences reduced lateral pressure and increased boundary slip. This boundary slip subsequently vanishes, with an abrupt strengthening due to the dilation of the grain assembly against fluid surface tension as saturation approaches. Strike-slip motion and confinement in this system explain the strength anomaly, highlighting a critical role of lateral faults in landslide stability, particularly in cases where dynamics cannot be adequately explained by monitored pore-water pressure or basal friction.

Geophysical Research Letters↗
USGS1977

Acid-rain weathering of a metasedimentary rock basin, Herkimer County, New York

A geochemical balance of major chemical constituents of precipitation inflow with those of surface-water outflow was used to evaluate the effect of acid rainfall (average pH = 4.3) on weathering rates of metasedimentary rocks underlying a surface-water basin in New York State. Precipitation samples and surface-water samples show that approximately 45 percent of the net cation outflow is due to weathering by precipitation acidity. The remaining 55 percent of the cation load in the outflow is attributed to interactions between atmosphere, rock, and soil water in the unsaturated zone. Comparison of the acidity of acid precipitation and air-equilibrated distilled water suggests that since preindustrial time, precipitation acidity has probably increased weathering rates for this basin approximately 75 percent. The increased acidity of the precipitation may also increase the rate of dissolution of heavy metals from water-piping systems.

New York↗