USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Research of the U.S. Geological Survey”

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

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics↗
USGS2025

Extremophile hotspots linked to containerized industrial waste dumping in a deep-sea basin

Decaying barrels on the seafloor linked to DDT contamination have raised concerns about the public health implications of decades old industrial waste dumped off the coast of Los Angeles. To explore their contents, we collected sediment cores perpendicular to five deep-sea barrels. The concentration of DDT and its breakdown products were highly elevated relative to control sites yet did not vary with distance from the barrels, suggesting that they were not associated with the contamination. Sediment cores collected through white halos surrounding three barrels were enriched in calcite and had elevated pH. The associated microbial communities were low diversity and dominated by alkalophilic bacteria with metagenome-assembled genomes adapted to high pH. A solid concretion sampled between a white halo and barrel was composed of brucite, a magnesium hydroxide mineral that forms at high pH. Based on these findings, we postulate that leakage of containerized alkaline waste triggered the formation of mineral concretions that are slowly dissolving and raising the pH of the surrounding sediment pore water. This selects for taxa adapted to extreme alkalinity and drives the precipitation of “anthropogenic” carbonates forming white halos, which serve as a visual identifier of barrels that contained alkaline waste. Remarkably, containerized alkaline waste discarded >50 years ago represents a persistent pollutant creating localized mineral formations and microbial communities that resemble those observed at some hydrothermal systems. These formations were observed at one-third of the visually identified barrels in the San Pedro Basin and have unforeseen, long-term consequences for benthic communities in the region.

California↗
USGS2025

Chemical and isotopic fractionation during melt inclusion formation

Melt inclusions are used to study the origin and evolution of magmas. The extent to which they represent equilibrium melt compositions, however, critically hinges on the ratio of crystal growth rate to diffusion rate in melts. If the rate of crystal growth is limited by the supply of nutrients and the dissipation of unwanted components, the trapped melt will be depleted in compatible elements and enriched in incompatible elements. Despite widespread recognition of the potential for melt inclusions to trap disequilibrium compositions, their identification remains a challenge. Here, we demonstrate how stable isotopes provide a solution to this problem. Melt inclusions in basaltic tephra from Kīlauea volcano, Hawai‘i, display Mg isotope fractionation up to 0.49‰ ( 26 Mg/ 24 Mg) relative to bulk rock. These observations indicate that kinetic processes, such as chemical diffusion, are at play, as equilibrium mineral-melt isotope fractionation at high temperatures is an order of magnitude smaller. We propose that the heavy Mg isotopic compositions of the melt inclusions were generated by the faster diffusion and thus preferential incorporation of 24 Mg by the growing olivine, leaving the slower diffusant ( 26 Mg) behind. The incompatible elements (e.g., Ca, Ti, and K) are predicted to display similarly large fractionations when boundary layer effects are significant. These findings show that stable isotopes can be a useful tool in identifying disequilibrium compositions in melt inclusions.

Geology↗
USGS2026

Seasonal migrations and other movements

In the past 25 years new information has been obtained on the migrations and movements of mountain sheep (bighorn [Ovis canadensis], thinhorn [Ovis dalli]). This chapter provides a comprehensive overview of mountain sheep migration and other movements across their broad distribution in western North America. Across the range of mountain sheep, migrations and other seasonal movements define a complex movement portfolio that relates to and supports all aspects of mountain sheep ecology. Bighorn and thinhorn species have examples of migratory behaviors that span the continuum between annual residents and long-distance migrants. Migratory behaviors can be characterized as low- or high-elevation residents; elevational migrations and other variants; and geographic migrants. Native populations that have been extant on the landscape without notable human intervention have greater migratory propensity and more diverse migratory portfolios due to the maintenance of migration through cultural learning and social transmission. Restored and augmented populations, where the population-level knowledge of migration has been lost or greatly reduced, are largely nonmigratory, although translocations show some ability to restore short-distance elevational migrations. Seasonal spring and fall migrations are less common in desert bighorn sheep (O. canadensis spp.) or bighorn sheep living in canyon or prairie breaks landscapes. Mineral lick visitation is important and common across the range of mountain sheep. Managing for migratory diversity can help to sustain migratory behavior in the face of climate change and other anthropogenic pressures, which can limit landscape connectivity between seasonal ranges or alter the spatiotemporal dynamics of regional phenology with cascading effects to other biotic and abiotic interactions such as the need to balance forage and predation risk or the spatial refugia with increased temperatures.

Book chapter↗
USGS2024

Seismic techniques and suggested instrumentation to monitor volcanoes

Introduction Changes in the pressure or location of magma can stress or break surrounding rocks and trigger flow of nearby waters and gases, causing seismic signals, such as discrete earthquakes and tremor. These phenomena are types of seismic unrest that commonly precede eruption and can be used to forecast volcanic activity. Mass movements at the surface, including avalanches, debris flows, and lahars, may also generate seismic signals that are specifically addressed in chapter H , this volume (Thelen and others, 2024). Our focus in this chapter is to determine the levels of instrumentation recommended to produce high-quality, well-constrained seismic observations important for early warning of impending eruptions, detecting changes in ongoing eruptions, and characterizing other hazardous volcanic events. There are emerging techniques and new types of instrumentation, such as distributed acoustic sensing or rotational seismometers, that we do not consider here. These types of instrumentation show promise for monitoring but still require maturation before being considered more generally in volcano monitoring. Most of the capabilities mentioned below are universal for all types of volcanic systems, although some are best applied to stratovolcanoes with an apical single vent. In some settings, such as calderas or shield volcanoes, we must broaden coverage to include multiple possible storage regions or vent locations. As an example, Thelen (2014) discretized the long rift zones of shield volcanoes in Hawaiʻi as a set of evenly spaced “vents.” In this construct, each vent comes with recommendations, and several thousand network configurations were simulated to assess the effect on network quality levels and to determine the most efficient network design. The same process could be applied in a caldera setting or a volcanic field, where an evenly spaced grid of potential vents is considered. Localized recommendations for each unique system are beyond the scope of this report and left up to local experts to assess based on the conditions, restrictions, and requirements of each volcano.

Scientific Investigations Report↗
USGS2025

Constraining large magnitude event source and path effects using ground motion simulations

The purpose of this study is to use ground motion simulations to investigate ways in which source and path effects for large magnitude events can be represented in non-ergodic GMMs. While we initially developed computation techniques using CyberShake simulations, the range of magnitudes and source-site combinations is not adequate to replicate what is observed empirically. We therefore designed a new ground motion simulation study, which includes earthquakes with a large range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of simulations (M4-M7), we then develop a non-ergodic GMM with the simulation data. We find that the within-site residuals are dominated by the radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within-site residuals among multiple source realizations. Finally, after removing the source effects from the within-site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site with the large events, while the second approach considers all small events on the fault plane. The results indicate that the path effects of large events cannot be satisfactorily approximated with that of small events using either approach.

Conference Paper↗
USGS1998

Peak-discharge frequency and potential extreme peak discharge for natural streams in the Brazos River basin, Texas

The 2-, 5-, 10-, 25-, 50-, and 100-year peak discharges were estimated for 186 streamflow-gaging stations with at least 8 years of data for natural streams in and near the Brazos River Basin, Texas. Multiple regression equations were developed to estimate peak-discharge frequency for the 2-, 5-, 10-, 25-, 50-, and 100-year recurrence intervals for each of three hydrologic regions that compose the Brazos River Basin. The equations for each region are a function of significant basin characteristics (explanatory variables). The significant explanatory variables among six that were tested are the contributing drainage area and stream slope for regions 1 and 2 and the contributing drainage area for region 3. For the three sets of equations, the coefficient of determination ranges from 0.59 to 0.93, and the standard error ranges from 0.184 to 0.391 log units. A larger coefficient of determination and a lower standard error generally are associated with the equations for hydrologic regions 2 and 3. Statistics from the regression analysis allow computation of the prediction interval associated with a given significance level for a peak-discharge frequency estimate. The regression equations can be used to estimate peak discharges for sites at, near, or away from sites with streamflow-gaging stations. The potential extreme peak-discharge curves as related to contributing drainage area were estimated for each of the three hydrologic regions from measured extreme peaks of record at 186 sites with streamflow-gaging stations and from measured extreme peaks at 37 sites without streamflow-gaging stations in and near the Brazos River Basin. The potential extreme peak-discharge curves generally are similar for hydrologic regions 1 and 2, and the curve for region 3 consistently is below the curves for regions 1 and 2, which indicates smaller peak discharges.

Texas↗
USGS2026

Morphodynamics of snow barchans at Concordia Station, Antarctica

Barchans are crescent-shaped dunes that form and migrate under unimodal wind regimes in areas of low sediment availability. While sand barchan morphodynamics has been extensively studied, snow barchans remain poorly documented. Here, using hourly time-series images taken near Concordia Station, Antarctica, we measure the sizes (lengths 0.3 m to 13 m, widths 0.4 m to 7.5 m) and migration rates (0.13 m h −1 to 3.5 m h −1 ) of 44 populations of snow barchans during 2017–2022. We compare this dataset with established height–velocity and length-width relationships for sand barchans in both air and water. Our results show that snow barchans follow the same scaling as sand barchans, allowing us to estimate a characteristic sintering (bonding of snow particles) time of about one day. This sintering process limits the size of snow barchans when active, reduces their aspect ratios, and, once wind falls below the transport threshold, explains their preservation over a few years. During wind events, migrating barchans are too small to maintain their shape and orientation, making them highly responsive to fluctuations in wind direction. Thus, studying snow barchans improves understanding of an elementary bedform encountered in multiple environments and also of the surface properties of snow on different timescales, offering insights into snow accumulation and wind-blown snow in ice cores.

Geomorphica↗
USGS1997

Hydrologic and water-quality conditions in the Horse Creek basin, west-central Florida, October 1992-February 1995

A baseline study of the 241-square-mile Horse Creek basin was undertaken from October 1992 to February 1995 to assess the hydrologic and water-quality conditions of one of the last remaining undeveloped basins in west-central Florida. During the period of the study, much of the basin remained in a natural state, except for limited areas of cattle and citrus production and phosphate mining. Rainfall in 1993 and 1994 in the Horse Creek basin was 8 and 31 percent, respectively, above the 30-year long-term average. The lowest and highest maximum instantaneous peak discharge of the six daily discharge stations occurred at the Buzzard Roost Branch and the Horse Creek near Arcadia stations with 185 to 4,180 cubic feet per second, respectively. The Horse Creek near Arcadia station had the lowest number of no-flow days with zero days and the Brushy Creek station had the highest number with 113 days. During the study, the West Fork Horse Creek subbasin had the highest daily mean discharge per square mile with 30.6 cubic feet per second per square mile, and the largest runoff coefficient of 43.7 percent. The Buzzard Roost Branch subbasin had the lowest daily mean discharge per square mile with 5.05 cubic feet per second per square mile, and Brushy Creek and Brandy Branch shared the lowest runoff coefficient of 0.6 percent. Brandy Branch had the highest monthly mean runoff in both 1993 and 1994 with 11.48 and 19.28 inches, respectively. During the high-baseflow seepage run, seepage gains were 8.87 cubic feet per second along the 43-mile Horse Creek channel. However, during the low-baseflow seepage run, seepage losses were 0.88 cubic foot per second. Three methods were used to estimate average annual ground-water recharge in the Horse Creek basin: (1) well hydrograph, (2) chloride mass balance, and (3) streamflow hydrograph. Estimated average annual recharge using these three methods ranged from 3.6 to 8.7 inches. The high percentage of carbonate plus bicarbonate analyzed at the Carlton surficial aquifer well could indicate an upward ground-water flow from the underlying intermediate aquifer system. Based on constituent concentrations in water samples from the six daily discharge stations, concentrations generally are lower in the upper three subbasins, West Fork Horse Creek, Upper Horse Creek, and Brushy Creek than in the lower three subbasins. Typically, concentrations were highest for major ions at Buzzard Roost Branch and nutrients at Brushy Creek.

Florida↗
USGS2025

Long-term trends in microseismicity during operational shut-ins at the Coso Geothermal Field, California

Pausing injection and production can lead to induced seismicity in a variety of settings, with some of the largest events occurring during these so-called shut-ins. In geothermal fields, shut-ins are periodically conducted for maintenance on wells and surface infrastructure, thereby offering recurring means of estimating stress changes in the subsurface that lead to increased seismicity rates. Here, we examine daily production and microseismicity data from the Coso Geothermal Field (CGF) in California between 1996 and 2010 to investigate the repetitive trends of operational shut-in microseismicity. Using 24 local seismic stations, we first analyze spatial and temporal trends of over 60,000 earthquakes with magnitudes between -0.4 to 3.8. We find that the northern region exhibits no significant seismicity changes during shut-ins, whereas the rest of the field experiences induced seismicity during almost every shut-in with an increasing intensity towards the southern and eastern portions of the field, highlighting local differences in stress within the CGF. Additionally, we cluster the seismicity using waveform cross-correlation, revealing several earthquake clusters primarily occurring during shut-in periods. These observations suggest that certain fracture and fault sections respond quicker to changes in pore pressure and poroelastic stresses within the geothermal system, possibly highlighting main fluid pathways.

California↗
USGS2025

Rainfall thresholds for postfire debris-flow initiation vary with short-duration rainfall climatology

The size, frequency, and geographic scope of severe wildfires are expanding across the globe, including in the Western United States. Recently burned steeplands have an increased likelihood of debris flows, which pose hazards to downstream communities. The conditions for postfire debris-flow initiation are commonly expressed as rainfall intensity-duration thresholds, which can be estimated given sufficient observational history. However, the spread of wildfire across diverse climates poses a challenge for accurate threshold prediction in areas with limited observations. Studies of mass-movement processes in unburned areas indicate that thresholds vary with local climate, such that higher rainfall rates are required for initiation in climates characterized by frequent intense rainfall. Here, we use three independent methods to test whether initiation of postfire runoff-generated debris flows across the Western United States varies similarly with climate. Through the compilation of observed thresholds at various fires, analysis of the spatial density of observed debris flows, and quantification of feature importance at different spatial scales, we show that postfire debris-flow initiation thresholds vary systematically with short-duration rainfall-intensity climatology. The predictive power of climatological data sets that are readily available before a fire occurs offers a much-needed tool for hazard management in regions that are facing increased wildfire activity, have sparse observational history, and/or have limited resources for field-based hazard assessment. Furthermore, if the observed variation in thresholds reflects long-term adjustment of the landscape to local climate, rapid shifts in rainfall intensity related to climate change will likely induce spatially variable shifts in postfire debris-flow likelihood.

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

Wet meadow regeneration through restoration of biophysical feedbacks

Wet meadows are globally significant ecosystems that provide critical hydrological, ecological, and biogeochemical functions, yet their extent has declined dramatically due to land use changes and hydrologic alteration. These sedge-dominated wetlands exist at the drier end of the wetland gradient, maintained by shallow groundwater and periodic inundation. This paper is a global synthesis of the ecological, geomorphic, and hydrological dynamics of wet meadows, with an emphasis on alluvial systems, to inform effective restoration strategies. We compare wet meadows to other wetlands, classify them into palustrine, lacustrine, and alluvial types, then focus on alluvial wet meadows and discuss how their formation and persistence depend on ground and surface water interactions, sediment deposition and flow obstructions, all mediated by biological processes. In particular, we highlight the role of hydric graminoids in resisting erosion and maintaining soil cohesion, how beaver promote meadow persistence, and the significance of wet meadows as carbon sinks. We also present stratigraphic evidence demonstrating that incision, often triggered by anthropogenic activity or changing climate, is the primary mechanism of alluvial wet meadow degradation, resulting in water table decline and shifts in vegetation composition. Restoration requires reversing these incisional processes through techniques that elevate water tables, disperse flow and retain sediment—methods traditionally associated with either soil conservation or stream restoration. These include nature-based solutions that create obstructions such as beaver dams and their analogues, rock and wood-based obstructions and incision trench or gully filling and grading. Given their multifunctional value—including but not limited to flood attenuation, biodiversity support, and carbon sequestration—wet meadows warrant a focused restoration framework. This review advocates for a valley-floor scale restoration paradigm that integrates hydrological reconnection, sediment retention, and biological reinforcement to ensure long-term resilience of these systems in the face of changing climate and land use pressures.

Frontiers in Environmental Science↗
USGS2025

Fiber-imaged supershear dynamics in the 2024 Mw 7 Mendocino Fault earthquake

Fault structure and rupture physics are deeply intertwined, and observations of this coupling are critical for understanding earthquake behavior. Rupture propagation is observable at fine scales using dense seismic networks. Fiber-optic sensing allows for long-term deployments of ultradense arrays that enable high-resolution measurements of infrequent, large earthquakes. We recorded the 2024 moment magnitude ( M w ) 7 Mendocino Fault earthquake with a nearby fiber-optic array and imaged its behavior with seismic beamforming. The rupture propagated to the east at subshear velocity; stagnated near the Mendocino Triple Junction, a zone of structural complexity; and subsequently transitioned to supershear velocity. The correlation between source physics and structure shows how lithospheric heterogeneity affects first-order characteristics of earthquake ruptures. Our results also demonstrate the potential for fiber-optic sensing to improve real-time estimation of key parameters for early warning.

California↗
USGS2025

Triple-oxygen isotopic evidence of prolonged direct bioleaching of pyrite with O2

Sulfate is often touted as containing atmospheric oxygen whose isotopic signature can constrain redox, environmental conditions, and biological activity. Yet, the amount and isotopic fractionation associated with air-O 2 incorporation during sulfate formation is still debated, making its verification difficult. In this study, we identify a distinct, microbially dominated environment with the potential to preserve maximum signals of air-O 2 in sulfate. We report triple-oxygen isotope data for sulfate produced from pyrite oxidation in microbial and abiotic experiments, and from natural dissolved sulfate from the Rio Tinto, Spain, an acid mine drainage site. The oxygen isotope systematics of sulfate in these environments define a unique kinetic isotope effect associated with initial stage pyrite oxidation by Acidithiobacillus ferrooxidans that preserves >80 % oxygen from air-O 2 in sulfate. Unlike experiments, which evolve toward water-oxygen dominated sulfate on short time scales, Rio Tinto, Spain hosts a microbe rich environment with distinct geochemistry that maintains high O 2 -oxygen in sulfate. Therefore, in addition to containing isotopic records from water and air, sulfates can also contain a biosignature that is promising for understanding conditions on Mars and early Earth.

Earth and Planetary Science Letters↗
USGS2025

Young explosive eruptions from the Clear Lake volcanic field

The Clear Lake volcanic field is the northernmost and youngest field in a chain of volcanic fields in and near the California Coast Range mountains. For 2 million years, numerous eruptions have happened around (and through) Clear Lake. The most recent period of activity in the Clear Lake volcanic field probably started around 40,000 years ago and was mainly explosive eruptions concentrated on and near faults in and around the lake. The combination of hot magma and groundwater created explosive releases of hot steam. The craters, called maars, threw out volcanic ash, pumice, and lava fragments, which fell back to the ground and draped the surrounding landscape. The deposits are visible in many places around the southeast end of the lake.

California↗
USGS1989

Techniques for simulating flood hydrographs and estimating flood volumes for ungaged basins in east and west Tennessee

A dimensionless hydrograph developed for a variety of basin conditions in Georgia was tested for its applicability to streams in East and West Tennessee by comparing it to a similar dimensionless hydrograph developed for streams in East and West Tennessee. Hydrographs of observed discharge at 83 streams in East Tennessee and 38 in West Tennessee were used in the study. Statistical analyses were performed by comparing simulated (or computed) hydrographs, derived by application of the Georgia dimensionless hydrograph, and dimensionless hydrographs developed from Tennessee data, with the observed hydrographs at 50 and 75% of their peak-flow widths. Results of the tests indicate that the Georgia dimensionless hydrography is virtually the same as the one developed for streams in East Tennessee, but that it is different from the dimensionless hydrograph developed for streams in West Tennessee. Because of the extensive testing of the Georgia dimensionless hydrograph, it was determined to be applicable for East Tennessee, whereas the dimensionless hydrograph developed from data on streams in West Tennessee was determined to be applicable in West Tennessee. As part of the dimensionless hydrograph development, an average lagtime in hours for each study basin, and the volume in inches of flood runoff for each flood event were computed. By use of multiple-regression analysis, equations were developed that relate basin lagtime to drainage area size, basin length, and percent impervious area. Similarly, flood volumes were related to drainage area size, peak discharge, and basin lagtime. These equations, along with the appropriate dimensionless hydrograph, can be used to estimate a typical (average) flood hydrograph and volume for recurrence-intervals up to 100 years at any ungaged site draining less than 50 sq mi in East and West Tennessee.

Tennessee↗
USGS1991

Effects of septic-tank effluent on ground-water quality in northern Williamson County and southern Davidson County, Tennessee

An investigation of the potential contamination of ground water from septic tank systems blasted in bedrock in Williamson and Davidson Counties, Tennessee, was conducted during 1988-89. Water samples were collected from domestic and observation wells, springs, and surface-water sites in a residential subdivision in the northern part of Williamson County near Nashville. The subdivision has a high density of septic-tank field lines installed into blasted bedrock Water samples also were collected from a well located in an area of Davidson County where field lines were installed in 5 feet of soil. Samples were analyzed for major inorganic constituents, nutrients, total organic carbon, optical brighteners, and bacteria. Although results of analyses of water samples from wells indicate no effect of septic-tank effluent on ground-water quality at these sites, water from two springs located downgradient from the subdivision had slightly larger concentrations of nitrite plus nitrate (2.2 and 2.7 milligrams per liter N), and much larger concentrations of fecal coliform and fecal streptococci bacteria (2,000 to 3,200 and 700 to 900 colonies per 100 milliliters of sample, respectively), than other wells and springs sampled during 1988. Water from one of these springs contained optical brighteners, which indicates that septic-tank effluent is affecting ground-water quality.

Tennessee↗
USGS2026

Temperate wetlands lose climate-cooling capacity under warming

Temperate inland wetlands have long been viewed as climate regulators, with sustained CO₂ sequestration offsetting continuous CH₄ emissions over centennial–millennial timescales. However, the net climatic effect of these opposing fluxes remains uncertain under anthropogenic warming as they differ in both magnitude and atmospheric lifetime. We combine 19,000 CH₄ chamber measurements from 202 temperate inland wetlands with carbon burial rates from 16 dated sediment cores in North America’s Prairie Pothole Region to estimate greenhouse gas dynamics over two centuries. Here we show that CH₄ emissions rise by ~1.39 kg CH₄ ha⁻¹ yr⁻¹ over the historical period and are projected to approach 950 kg CH₄ ha⁻¹ yr⁻¹ under high-forcing scenarios by 2100, while CO₂ sequestration is expected to decline by ~0.3 kg CO₂ ha⁻¹ yr⁻¹, driving a shift toward persistent net CO₂-equivalent source behaviour and a weakening of the long-term cooling effect that wetlands have provided under natural conditions.

Prairie Pothole Region↗