SEARCH · USGS Science

Results for “Hydrology and Earth Systems Science”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,661 records · Page 87Linked to original sources

USGS2024

Lithium resource in the Smackover Formation brines of Southern Arkansas

Lithium-rich brine deposits occur throughout the United States, including in the Smackover Formation. The concentration of lithium in Smackover Formation brines was predicted across southern Arkansas by using a machine-learning model that incorporated lithium concentration data and geologic information. Between 5.1 and 19.0 million metric tons of lithium are calculated to be present in the brines of the Smackover Formation in southern Arkansas. The range in possible total lithium reflects the uncertainty in machine-learning predictions of lithium concentrations and the range of Smackover Formation porosity. This estimate quantifies the in-place lithium resource and does not consider the technological and economic feasibility of extracting the lithium from the brines.

Arkansas↗

USGS1976

Non-radiometric borehole geophysical detection of geochemical halos surrounding sedimentary uranium deposits

Roll-type uranium deposits are formed by the concentration of uranium by ground water in geochemical cells. Non-uranium minerals having different solubilities may be deposited ahead of or behind the uranium minerals, forming halos that surround the ore. In addition, oxidizing and reducing environmental conditions may cause zones of mineral alteration to develop beyond the limits of the uranium deposit. Certain physical property anomalies that are commonly associated with halos can be detected by relatively fast and inexpensive borehole geophysical measurements made either in individual holes, or between two adjacent holes. Borehole measurements that have been found to be useful include electrical resistivity, induced polarization, and magnetic susceptibility. Electrical resistivity is increased by the presence of calcite and other cementing minerals that sometimes create permeability barriers in the neighborhood of uranium deposits. Induced polarization (IP) response is increased by sulfide and clay minerals that are commonly found in anomalous concentrations near roll-type deposits. Magnetic susceptibility is usually decreased by the oxidation of magnetite to hematite or limonite in the zone of chemical alteration that is left as a trail behind roll fronts. Borehole measurements of electrical resistivity, induced polarization and magnetic susceptibility were made in the vicinity of a uranium roll-type deposit in south Texas. Results indicate that mineral halos can be detected by borehole measurements made. in wide-spaced drill holes, and that the total amount of drilling needed to find a deposit can be reduced substantially by this exploration approach.

Open-File Report↗

USGS2025

Site-specific, extended ShakeMaps for earthquake engineering applications

The U.S. Geological Survey (USGS) routinely produces ShakeMaps of shaking intensity across the globe. Due to practical constraints, the number of response spectral periods was limited to three standard periods (0.3, 1.0, and 3.0 sec). We have recently developed the tools that are necessary to expand this functionality to include 22 periods (matching the current U.S. National Seismic Hazard Model periods) as well as the orientation-independent components (e.g., “RotD50”). We refer to ShakeMap products that include these extensions as “extended ShakeMaps.” The added level of complexity motivated us to also develop a user-friendly tool called the “ShakeMap Sampling Tool” (SST) that gives all the estimated shaking metrics for a specific location (or list of locations). Additionally, we develop a web application where users can input locations of interest and view/download the SST results. We further familiarize users with the concept of “Composite ShakeMaps.” For earthquakes sequences such as a mainshock and larger foreshocks and aftershocks, this provides a map of the maximum value of each shaking metric, which is useful for overall loss estimates, the full extent of ground failure triggering potential, and a better portrayal of the repeated shaking levels at a given point for a series of earthquakes. Such a site-specific shaking history facilitates earthquake forensics at building or infrastructure sites for which damage may be of concern, as described in the Disproportionate Damage Earthquake trigger specified in the IEBC (2018, Section 405.2.2) and in developing ATC-145 guidelines (Guidelines for Post-Earthquake Assessment, Repair, and Retrofit of Buildings). The composite ShakeMap can be combined with the SST for a variety of earthquake-hazard applications, such as systematically inferring triggering shaking estimates at specific sites of geotechnical interest for landsliding, liquefaction, and lateral-spreading hazards.

Kahramanmaraş↗

USGS2025

Climate-driven sulfate export in alpine watersheds may stimulate methylmercury production

Climate change is increasing sulfate export and changing wetland extent in mountain regions. These changes may increase microbially mediated production of the neurotoxic substance methylmercury due to enhanced sulfate metabolism in mountain environments. Here, we assess methylmercury concentrations and formation rates across high-elevation wetlands in the Colorado Rocky Mountains. We also investigate sulfate controls on methylmercury production within subalpine peatlands by amending soils with sulfate to mimic increased stream export of sulfate from the alpine zone and measuring methylmercury formation rates for different sulfate treatments. We found that subalpine peatlands have statistically significant higher methylmercury concentrations and formation rates compared to alpine, mineral-soil wetlands. Methylmercury production in subalpine peatlands also increased significantly ( p < 0.05) following sulfate additions; the highest rates occurred in sediments with intermediate extractable sulfate concentrations (∼0.60–1.4 mg sulfate g −1 dry soil). Our study is the first to identify soil sulfate-related thresholds for methylmercury production and sulfate-limitation of methylmercury production in subalpine peatlands. These findings highlight important linkages between climate-driven mineral weathering and mercury cycling in mountain regions globally.

Colorado↗

USGS2024

Cross-fade sampling: Extremely efficient Bayesian inversion for a variety of geophysical problems

This paper introduces cross-fade sampling, a computationally efficient Markov Chain Monte Carlo simulation method that uses a semi-analytical approach to quickly solve Bayesian inverse problems that do not themselves have an analytical solution. Cross-fading is efficient in two ways. First, it requires fewer samples to obtain the same quality simulation of the target probability density function (PDF). Secondly, it is much faster to evaluate the posterior probability of each sample than conventional sampling methods for simulating Bayesian posterior PDFs. Conventional methods require evaluating the prior probability (which describes your a priori constraints) and data likelihood (which describes the fit between the observations and the predictions of the model) for each sample model. However, cross-fading does not require evaluating the data likelihood, meaning that ‘big data’ can be fit with zero additional computational cost. Further, the cross-fading approach can be used to calculate the marginal likelihood associated with a model design, facilitating model comparison and Bayesian model averaging. Topics covered in this paper include derivation of the cross-fade approach and how it can be used to simulate Bayesian posterior PDFs and compute the marginal likelihood, discussion of the class of problems to which cross-fading can be applied (with examples from earthquake statistics, earthquake ground motion modelling, volcanic eruption forecasting, and finite fault slip modelling), demonstration of efficiency relative to existing sampling methods and discussion of how cross-fading can be used to account for prediction errors (i.e. epistemic errors) as part of the geophysical inverse problem.

Geophysical Journal International↗

USGS1990

Preconditioned Conjugate-Gradient 2 (PCG2), a computer program for solving ground-water flow equations

This report documents PCG2: a numerical code to be used with the U.S. Geological Survey modular three-dimensional, finite-difference, ground-water flow model. PCG2 uses the preconditioned conjugate-gradient method to solve the equations produced by the model for hydraulic head. Linear or nonlinear flow conditions may be simulated. PCG2 includes two reconditioning options: modified incomplete Cholesky preconditioning, which is efficient on scalar computers; and polynomial preconditioning, which requires less computer storage and, with modifications that depend on the computer used, is most efficient on vector computers. Convergence of the solver is determined using both head-change and residual criteria. Nonlinear problems are solved using Picard iterations. This documentation provides a description of the preconditioned conjugate gradient method and the two preconditioners, detailed instructions for linking PCG2 to the modular model, sample data inputs, a brief description of PCG2, and a FORTRAN listing.

Water-Resources Investigations Report↗

USGS2026

Rapid characterization of the 2025 Mw 8.8 Kamchatka, Russia earthquake

The 29 July 2025 M w 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near‐source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real‐time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite‐fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The M w 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

Kuril-Kamchatka subduction zone↗

USGS2026

Incorporating location uncertainty improves inference with stop-level North American Breeding Bird Survey data

Ecological models should account for uncertainty to be most effective and useful. Yet, uncertainty from model covariates—unlike that from other sources, such as sampling error or process variability—is seldom explicitly incorporated. This can cause underestimates of uncertainty to cascade through model parameter estimates, predictions, and downstream uses. Burner et al. proposed a method for quantifying uncertainty in covariates and incorporating it into models using informative Bayesian priors. This method was applied to stop-level Breeding Bird Survey (BBS) analyses, where land cover uncertainty at each stop arises from substantial stop location uncertainty. A limited validation of model-estimated land cover, using stops with known locations, indicated the method’s potential effectiveness, but it was not rigorously evaluated. We conduct a robust simulation-based test, generating stop locations, extracting land cover, and simulating bird communities across 210 BBS routes in the upper Midwest. We compare 3 models: a “known” model with true land cover, a “naive” model assuming consistent 800-m stop spacing, and a “full” model using informative priors to estimate land cover. Species parameter estimates and predicted prevalence patterns across gradients in land cover from the full model approached those of the known model and were substantially closer to the true values used in simulations relative to those from the naive model. Naive model parameters were more biased relative to the other models, and credible intervals of predicted species prevalence rarely included the true simulated values. The full model also produced land cover covariate estimates closer to true simulation values relative to the mean informative priors. Our results show that, for the BBS, informative priors enable more accurate stop-level analyses despite location uncertainty. In contrast, naive models that ignore this uncertainty yield poor inferences. More broadly, we demonstrate empirically the utility of informative priors to account for covariate uncertainty in ecological models.

Michigan, Minnesota, Wisconson↗

USGS2023

Stratigraphic notes—Volume 1, 2022

This is the first volume in the U.S. Geological Survey (USGS) series of reports on stratigraphy entitled “Stratigraphic Notes,” which consists of short papers that highlight stratigraphic studies, changes in stratigraphic nomenclature, and explanations of stratigraphic names and concepts used on published geologic maps. “Stratigraphic Notes” is a long-term (multiyear), multivolume publication containing articles that address updates or revisions to stratigraphic nomenclature (and whose content ultimately will be incorporated by National Geologic Map Database personnel into Geolex, https://ngmdb.usgs.gov/Geolex/ ). We welcome papers for the “Stratigraphic Notes” series from geoscientists of the USGS, of State Geological Surveys, and from academicians. Papers can be submitted for publication in “Stratigraphic Notes” by contacting the USGS Geologic Names Committee ( gnc@usgs.gov ). As new “Stratigraphic Notes” volumes are published, links to the volumes will be posted at https://doi.org/10.3133/pp1879 . This first volume ("Stratigraphic notes—Volume 1, 2022") includes articles that provide guidance for those who wish to submit papers to “Stratigraphic Notes,” as well as information on how to make your manuscripts compliant for geologic names reviews and how to organize your paper’s content to facilitate inclusion of new or revised names in Geolex. This volume also includes some specific guidance on conducting geologic names reviews of geologic and hydrogeologic reports.

Professional Paper↗

USGS2025

Fault displacement model for surface principal rupture of strike-slip faults

The probability distribution model for principal displacement accommodated on the surface main trace is a critical input to the fault displacement hazard analysis. This article presents a new model for strike-slip ruptures in the moment magnitude ( M ) range of 6 to 8.3. The new model is the outcome of a multi-year research effort to update the widely used model developed by Petersen and others in 2011. Updates include the adoption of the Fault Displacement Hazard Initiative database and enhancements to rupture and displacement data preparation. Statistical formulation and estimation have also been updated substantially. A three-parameter modified normal distribution that we refer to as the negative Exponentially Modified Gaussian distribution is adopted to model the probability distribution of the natural logarithm of principal displacement. Formulation for the mean parameter of the modified normal includes a random earthquake term, a nonlinear scaling relation with M , and an ellipse function for along-main-trace variation. The aleatory variability of the updated model now depends on M as well as site’s along-main-trace position. These updates not only significantly improve the fit to the distribution of the observed displacements but also yield reasonable 95th percentile predictions for M > 7.5 events. Alternative models representing the estimation uncertainty of the M -scaling relation are also developed. These new models are compared to the previous model in terms of percentile predictions and the calculated hazard curves. The steeper hazard curves from the new models yield a lower exceedance rate than the normal-distribution based model developed previously by Petersen and others.

Earthquake Spectra↗

USGS1986

Potential for leakage among principal aquifers in the Memphis area, Tennessee

The principal aquifers in the Memphis area consist primarily of sand or sand and gravel, and the confining beds consist of clay, silt, sand, and lignite. The Jackson Formation and upper part of the Claiborne Group serve as the confining bed separating the water table aquifers from the Memphis Sand; the Flour Island Formation separates the Memphis Sand from the Fort Pillow Sand. Differences in total hydraulic head among the principal aquifers in the Memphis urban area result in vertical hydraulic gradients which create a potential for inter-aquifer exchange of water. Throughout this area, the gradient is downward from the water table aquifers to the Memphis Sand. In the central part of the Memphis urban area, the vertical hydraulic gradient is upward from the Fort Pillow Sand to the Memphis Sand, and in the eastern and western parts, it is downward from the Memphis Sand to the Fort Pillow Sand. The vertical distribution of carbon-14 data for water from the fluvial deposits, Memphis Sand, and Fort Pillow Sand shows an increase in the relative age of the water with depth. The areal distribution of carbon-14 data for water from the upper part of the Memphis Sand indicates that relatively recent water has been brought into the major cone of depression in the potentiometric surface of the Memphis Sand, either by horizontal movement or from downward vertical leakage. The normal, near-surface geothermal gradient in the Memphis area was determined to be 0.6 C/100 ft. Deviations from the normal geothermal gradient, in areas affected by intense pumping from the Memphis Sand, indicate that downward vertical leakage occurs from the water table aquifers through the Jackson-upper Claiborne confining bed to the Memphis Sand. The velocity of downward vertical leakage of water from the Memphis Sand through the Flour Island confining bed to the Fort Pillow Sand was determined to be 0.0066 ft/day by analysis of borehole temperature data from an observation well in the northeastern part of the Memphis area. From this velocity and the head difference between the Memphis Sand and the Fort Pillow Sand at this locality, the hydraulic conductivity of the Flour Island confining bed was determined to be 0.00114 ft/day. (Lantz-PTT)

Tennessee↗

USGS2024

Prediction of regional broadband strong ground motions using a teleseismic source model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, earthquake

To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( ⁠M w ⁠ ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.

Papanoa↗

USGS2024

Updating the crustal fault model for the 2023 National Seismic Hazard Model for Alaska

We present the crustal fault model for Alaska, based on geologic observations, as a primary input for the 2023 revision of the U.S. Geological Survey National Seismic Hazard Model. We update the 2013 Alaska Quaternary fault and fold database to produce a simplified model of 105 fault sections and four fault zone polygons with basic geologic parameters including slip sense and rate. Significant updates include the following: (1) a slip rate of ∼53 mm/year on the Queen Charlotte Fault indicating it accommodates all of the plate boundary motion; (2) quantified slip rates on megathrust splay faults in the southern Prince William Sound region and near Kodiak Island; (3) improved details of structures in the Chugach-St. Elias orogen; (4) revision of the Castle Mountain Fault from right-lateral slip to a predominantly reverse fault; (5) improved Interior Alaska tectonic models that clarify relationships between the Denali, Totschunda, and thrust faults on both sides of the Alaska Range; (6) identified large earthquake sources in the eastern Brooks Range; and (7) omission of the Chatham Strait section of the Denali Fault. The fault model underscores that the collision of the Yakutat microplate is the dominant driver of active crustal faulting in most of Alaska.

Alaska, Yukon↗

USGS2023

Evaluation of models for estimating hydraulic conductivity in glacial aquifers from NMR logging

Nuclear magnetic resonance (NMR) logging is a promising method for estimating hydraulic conductivity ( K ). During the past ∼60 years, NMR logging has been used for petroleum applications, and different models have been developed for deriving estimates of permeability. These models involve calibration parameters whose values were determined through decades of research on sandstones and carbonates. We assessed the use of five models to derive estimates of K in glacial aquifers from NMR logging data acquired in two wells at each of two field sites in central Wisconsin, USA. Measurements of K , obtained with a direct push permeameter (DPP), K DPP , were used to obtain the calibration parameters in the Schlumberger-Doll Research, Seevers, Timur-Coates, Kozeny-Godefroy, and sum-of-echoes (SOE) models so as to predict K from the NMR data; and were also used to assess the ability of the models to predict K DPP . We obtained four well-scale calibration parameter values for each model using the NMR and DPP measurements in each well; and one study-scale parameter value for each model by using all data. The SOE model achieved an agreement with K DPP that matched or exceeded that of the other models. The Timur-Coates estimates of K were found to be substantially different from K DPP . Although the well-scale parameter values for the Schlumberger-Doll, Seevers, and SOE models were found to vary by less than a factor of 2, more research is needed to confirm their general applicability so that site-specific calibration is not required to obtain accurate estimates of K from NMR logging data.

Wisconsin↗

USGS2026

Miocene evolution of the Humboldt Current

Diatom records from the East Pisco Basin (EPB) of southern Peru and offshore Ocean Drilling Program (ODP) Hole 682 A reveal stepwise increases in the primary productivity of the Humboldt Current during the middle and late Miocene. Although diatoms are present back through the late middle Eocene, successively enhanced diatom production occurs during the Miocene in four steps. The first step between 14.2 and 13.8 Ma marks the onset of diatom deposition in the Pisco-0 sequence. A second step at ∼12.9 Ma coincides with a major drop in global sea level during which diatom deposition ceased in the EPB but continued in ODP 682 A as well as in Chile, Ecuador, and Colombia. Beginning at 10.4 to 10.0 Ma, a major increase in diatom sedimentation rates in both the EPB and in ODP 682 A signals an abrupt intensification of the Humboldt Current productivity, which may be related to the closure of the Central American Seaway to deeper water circulation. A fourth step coincides with the transgressive deposition of the diatom-rich Pisco-2 sequence at 8.4 Ma and is enhanced by the onset of the global Late Miocene Biogenic Bloom (8.0–4.5 Ma). Between 7 and 6 Ma, common subtropical diatoms in ODP 682 A contrast with high diatom deposition rates in the EPB, suggesting variable El Niño-La Niña conditions. During this time, fossil vegetation evidence from southern Peru supports the presence of seasonal periods of enhanced rainfall typical of El Niño conditions. The Miocene trend of increasing sea mammal diversity in the EPB parallels these steps, further supporting stepwise enhancement of primary productivity in the Humboldt Current during the Miocene.

Humboldt Current↗

USGS2025

Landslide-channel feedbacks amplify channel widening during floods

Channel widening is a major hazard during floods, particularly in confined mountainous catchments. However, channel widening during floods is not well understood and not always explained by hydraulic variables alone. Floods in mountainous regions often coincide with landslides triggered by heavy rainfall, yet landslide-channel interactions during a flood event are not well known or documented. Here we demonstrate with an example from the Great Colorado Flood in 2013, a 1000 year precipitation event, how landslide-channel feedbacks can substantially amplify channel widening and flood risk. We use a combination of DEM differencing, field analysis, and multiphase flow modeling to document landslide-channel interaction during the flood event in which sediment delivered by landslides temporarily dammed the channel before failing and generating substantial channel widening. We propose that such landslide-flood interactions will become increasingly important to account for in flood hazard assessment as flooding and landsliding both increase with extreme rainfall under climate change.

Colorado↗

USGS2024

Automated, near real-time ground-motion processing at the U.S. Geological Survey

We describe automated ground‐motion processing software named gmprocess that has been developed at the U.S. Geological Survey (USGS) in support of near‐real‐time earthquake hazard products. Because of the open‐source development process, this software has benefitted from the involvement and contributions of a broad community and has been used for a wider range of applications than was initially envisioned. Here, we give an overview and introduction to the software, including how it has leveraged other open‐source libraries. We highlight some key features that gmprocess provides, compare response spectra calculated with the automated processing approach of gmprocess to the response spectra provided by the Next Generation Attenuation projects, and summarize projects that have utilized gmprocess. These use‐cases demonstrate that this software development effort has been successfully leveraged in earthquake research activities both within and outside the USGS.

Seismological Research Letters↗

USGS2025

Bølling-Allerød productivity in the subarctic Pacific driven by seasonal upwelling

The Bølling-Allerød deglacial event is marked by high diatom productivity and opal deposition throughout the subarctic Pacific. This opal could either constitute a strengthened biological pump and thus carbon sequestration, or a weakened biological pump and release of marine-sequestered CO 2 to the atmosphere. We quantify silicic acid supply at IODP Site U1340 in the Bering Sea using biogenic opal and δ 30 Si of Coscinodiscus , a diatom genus. These records, along with diatom environmental indicators, suggest the Bølling-Allerød had high silicic acid availability related to a shift from stratification to seasonal upwelling dynamics. We thus propose the primary cause of the high productivity event was increased macronutrient supply from vertical exchange that injected old, nutrient-rich, CO 2 -rich waters into the surface. Enhanced CO 2 release from the subarctic Pacific may help explain critical intervals of CO 2 rise that occur at the onsets of the Bølling-Allerød and PreBoreal.

Geophysical Research Letters↗