USGS Science⌕ Search

SEARCH · USGS Science

Results for “Communications Earth 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,662 records · Page 89Linked to original sources
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↗
USGS2025

Impact of land subsidence on housing sale values: Evidence from the San Joaquin Valley, California

This study assesses the impact of land subsidence on housing sale values in the San Joaquin Valley, California. The study utilizes home sale transactions and vertical land-surface displacement data from Interferometric Synthetic Aperture Radar techniques. Using fine-scale fixed effects, matching, as well as a repeat-sales approach, our results indicate that land subsidence resulted in a 2.4% to 5.8% reduction in housing sale values, with the largest reductions occurring in areas where substantial subsidence occurred. Such findings may have implications for groundwater management and can potentially help inform policy design to help mitigate the causes and impacts of land subsidence.

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

Groundwater age estimates for the Mississippi River Valley alluvial aquifer based on tracer data collected during 2018–20

This study characterized groundwater age across the Mississippi River Valley alluvial aquifer (MRVA). Groundwater samples from 69 MRVA wells and 19 wells in Tertiary units of the Mississippi embayment aquifer system (MEAS) were analyzed for sulfur hexafluoride (SF 6 ), tritium ( 3 H), helium (He), and (or) carbon-14 of dissolved inorganic carbon ( 14 C). The age distributions of 89 samples were estimated by fitting lumped parameter models to processed tracer concentrations with the U.S. Geological Survey software TracerLPM. Mean ages of MRVA groundwater samples ranged from 12 to 22,000 years, with a median of 140 years. Mean ages of MEAS groundwater samples ranged from 230 to 52,000 or more years, with a median of 13,500 years. The spatial distribution of MRVA groundwater ages was found to be influenced by depth, inflow of groundwater from deeper units, and soil saturated hydraulic conductivity. In parts of the MRVA, the spatial distribution of MRVA groundwater ages was found to be influenced by annual recharge and (or) annual groundwater pumpage.

Alabama, Arkansas, Illinois, Kentucky, Louisiana, ↗
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

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

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

Estimating the probability of elevated nitrate (NO2+NO3-N) concentrations in ground water in the Columbia Basin Ground Water Management Area, Washington

Logistic regression was used to relate anthropogenic (man-made) and natural factors to the occurrence of elevated concentrations of nitrite plus nitrate as nitrogen in ground water in the Columbia Basin Ground Water Management Area, eastern Washington. Variables that were analyzed included well depth, depth of well casing, ground-water recharge rates, presence of canals, fertilizer application amounts, soils, surficial geology, and land-use types. The variables that best explain the occurrence of nitrate concentrations above 3 milligrams per liter in wells were the amount of fertilizer applied annually within a 2-kilometer radius of a well and the depth of the well casing; the variables that best explain the occurrence of nitrate above 10 milligrams per liter included the amount of fertilizer applied annually within a 3-kilometer radius of a well, the depth of the well casing, and the mean soil hydrologic group, which is a measure of soil infiltration rate. Based on the relations between these variables and elevated nitrate concentrations, models were developed using logistic regression that predict the probability that ground water will exceed a nitrate concentration of either 3 milligrams per liter or 10 milligrams per liter. Maps were produced that illustrate the predicted probability that ground-water nitrate concentrations will exceed 3 milligrams per liter or 10 milligrams per liter for wells cased to 78 feet below land surface (median casing depth) and the predicted depth to which wells would need to be cased in order to have an 80-percent probability of drawing water with a nitrate concentration below either 3 milligrams per liter or 10 milligrams per liter. Maps showing the predicted probability for the occurrence of elevated nitrate concentrations indicate that the irrigated agricultural regions are most at risk. The predicted depths to which wells need to be cased in order to have an 80-percent chance of obtaining low nitrate ground water exceed 600 feet in the irrigated agricultural regions, whereas wells in dryland agricultural areas generally need a casing in excess of 400 feet. The predicted depth to which wells need to be cased to have at least an 80-percent chance to draw water with a nitrate concentration less than 10 milligrams per liter generally did not exceed 800 feet, with a 200-foot casing depth typical of the majority of the area.

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