Geology ReportsSearch

SEARCH · Geology Reports

Results for “Natural Hazards and Earth Systems Sciences (NHESS)”

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,668 records · Page 88Linked to original sources

Unscrambling the Proterozoic supercontinent record of northeastern Washington State, USA

The time interval from Supercontinent Nuna assembly in the late Paleoproterozoic to Supercontinent Rodinia breakup in the Neoproterozoic is considered by some geologists to comprise the “Boring Billion,” an interval possibly marked by a slowdown in plate tectonic processes. In northeastern Washington State, USA, similar to much of western Laurentia, early workers generally thought the tectonostratigraphic framework of this interval of geologic time consisted of two major sequences, the (ca. 1480–1380 Ma) Mesoproterozoic Belt Supergroup and unconformably overlying (<720 Ma) Neoproterozoic Windermere Supergroup. However, recent research indicates that strata considered by early workers as Belt Supergroup equivalents are actually younger, and a post-Belt, pre-Windermere record is present within the <1360 Ma Deer Trail Group and <760 Ma Buffalo Hump Formation. Thus, the northeastern Washington region perhaps comprises the most complete stratigraphic record of the “Boring Billion” time interval in the northwestern United States and holds important insights into global Proterozoic supercontinent tectonic processes. In light of these exciting developments, this field guide will address the early historic economic geology and original mapping of these Proterozoic sequences in the northeastern Washington region, and from that foundation explore more recent isotopic provenance data and their regional to global context. Finally, the guide will end with a discussion of remaining questions with a goal of stimulating interest in these relatively understudied, yet important, rocks.

Washington

Updated three-dimensional temperature maps for the Great Basin, USA

As part of the periodic update of the geothermal energy assessments for the USA (e.g., last update by Williams and others, 2008), a new three-dimensional temperature map has been constructed for the Great Basin, USA. Williams and DeAngelo (2011) identified uncertainty in estimates of conductive heat flow near land surface as the largest contributor to uncertainty in previously published temperature maps. The new temperature maps incorporate new conductive heat flow estimates developed by DeAngelo and others (2023). Predicted temperatures at depth are compared with representative measurements (for conductively dominated conditions), showing good agreement under relatively simple uniform conditions. Inputs included radiogenic heat production for all layers of 1.89 μW/m 3 , effective bulk thermal conductivity of 2.7 W/m/°C for all rocks underlying sedimentary basins, and a previously published (Williams and DeAngelo, 2011) empirically driven estimate of increasing thermal conductivity with depth in sedimentary sequences. The resulting three-dimensional temperature model is published in a USGS data release associated with this manuscript (Burns and others, 2023).

Arizona, California, Idaho, Nevada, Oregon, Utah

Thickness of the Saudi Arabian crust

As part of a joint Saudi Geological Survey (SGS) and U.S. Geological Survey (USGS) project, we analyzed P-wave receiver functions from seismic stations covering most of the Kingdom of Saudi Arabia to map the thickness of the crust across the Arabia Plate. We present an update of crustal-thickness estimates and fill in gaps for the western Arabian Shield and the rifted margin at the Red Sea (the coastal plain), as well as the eastern Arabian Platform. We applied a conventional H-k stacking algorithm and included careful attention to stacking weights, two forms of sedimentary corrections for stations located on the Arabian Platform, and additional processing for noisy stations. We obtained useful results at 154 stations from 898 teleseismic events over a 2-year period from 1995–1997 (for non-SGS stations) and a 6-year period from 2008–2014 (for SGS stations). Average crustal thickness (that is, depth to the Mohorovičić discontinuity [Moho] below the surface) beneath the Red Sea coastal plain (the rift margin) is 29 kilometers (km), beneath the volcanic fields (known in Arabic as harra [plural] or harrat [singular]) is 35 km, beneath the Arabian Shield (excluding harrats) is 37 km, and beneath the Arabian Platform is 38 km. Crustal thinning appears not to extend east of the rift escarpment, suggesting uniform extension that is no broader at depth than at the surface. In contrast to some previous claims that the Arabian Platform crust is thicker than that of the Arabian Shield, we find no statistically significant difference between their whole crustal thicknesses. However, the average subsedimentary crustal thickness (that is, the crystalline crust) for stations on the Arabian Platform is 34 km, 3 km thinner than the crust of the Arabian Shield. Individual station P-wave (pressure) velocity and S-wave (shear) velocity ratios ( V P / V S ) are highly variable for the Arabia Plate, ranging from 1.60 to 1.97 and averaging 1.75, with a standard deviation of 0.07. There are no statistically significant differences between V P / V S ratios of the different geologic regions of Saudi Arabia. Similar V P / V S ratios, coupled with similar crustal thicknesses for harrats and the Arabian Shield, indicate that Cenozoic magmatism has contributed negligibly to crustal growth.

Professional Paper

The U.S. Geological Survey National Water Quality Network—Groundwater, 2024

The U.S. Geological Survey (USGS) operates a National Water Quality Network (NWQN) to monitor trends in groundwater quality and assess emerging contaminants of concern. It is a “network of networks” with 81 subnetworks being sampled on a decadal time scale. Each year, 8 of the subnetworks are sampled. Subnetworks have 20–30 wells each and include studies of domestic supply wells or shallow groundwater (20–50 feet deep) underlying urban land use or agricultural land use. Currently there are 2,089 wells in the network. All wells are sampled for physical properties, nutrients, major ions, trace elements, per- and polyfluoroalkyl substances (PFAS), and a subset of wells are sampled for pesticides, volatile organic compounds, radionuclides, and microbiological contaminants.

Contiguous United States

Rapid hydrothermal triggering of induced seismicity at the Coso Geothermal Field

The long-term producing Coso Geothermal Field (CGF) in California operates over 100 wells tapping into a reservoir characterized by an extensive fracture network, complex fluid pathways, and regular seismic activity. Understanding the interaction between seismicity and injection can shed important light on the hydrothermal characteristics of the field. Here, we analyze 15 years of local seismic and daily operational data from the CGF, identifying a strong correlation between short-term increase in seismicity rate and seasonal volumetric and temperature variations in the reinjected fluid. Furthermore, the seismic footprint during peak injection of colder fluids reveals a near-instantaneous response up to 2 km away from the injection well, too rapid for pore pressure diffusion alone. This short-term and distant response is observed to have directional preference, indicating structural or permeability anisotropy within the reservoir. Additionally, the seismic response correlates with the initial volumetric increase of colder fluids, but also with temperature decrease during stable injection periods, suggesting thermal effects alone can play an important role in triggering distant seismicity.

California

3D Converted wave reverse time migration imaging

We describe a newly developed method for recovering high-resolution images of seismic discontinuities, such as subducting slabs, in 3D. Our method makes use of converted S → P or P → S waves observed by dense arrays of seismometers to infer the locations and relative strengths of seismic discontinuities at depth in a target region. Observed direct and converted waves are backpropagated to their times of origin. The time-reversed wavefield is then separated into its constituent P and S components via the Helmholtz decomposition, and those separated wavefields are used to compute imaging functions that characterize the locations and relative strengths of seismic discontinuities. Imaging functions may be designed to use either S → P or P → S waves, so that users can target those arrivals expected to be most dominant in a given dataset. We have previously demonstrated the efficacy of our method in two dimensions, and we now present a 3D implementation of our technique which addresses the significant computational challenges posed by the size of volumetric wavefield data in three dimensions. Through a series of synthetic examples, we demonstrate that our method is capable of recovering the fine scale structure of a subducting slab given realistic station coverage and earthquake sources. We investigate optimal seismic station geometries for our technique and explore image interpretability in regions with poor data coverage. We find that linear station geometries yield more optimal, interpretable imaging functions than collections of small arrays can. We also show that our method can successfully recover bothS → P or P → S images when realistic shear earthquake sources are used, and we explore the additional computational challenges presented by the high frequency content of S waves. Our results demonstrate the potential for our technique to recover high-resolution information about subducting slabs in real-world regions, given that relatively sparse seismic arrays with only approximately 100 stations are capable of recovering interpretable imaging functions from just a few realistic earthquake sources for multiple discontinuities at significant depth in an area of approximately 400~sq~km.

Seismica

The impact of source time function complexity on stress drop estimates

Earthquake stress drop—a key parameter for describing the energetics of earthquake rupture—can be estimated in several different, but theoretically equivalent, ways. However, independent estimates for the same earthquakes sometimes differ significantly. We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different estimates. We apply time‐ and frequency‐domain methods to estimate stress drops for real earthquakes in the SCARDEC (Seismic source ChAracteristics Retrieved from DEConvolving teleseismic body waves, Vallée and Douet, 2016 ) source time function (STF) database and analyze how rupture complexity drives stress‐drop estimate discrepancies. Specifically, we identify two complexity metrics—Brune relative energy (BRE) and spectral decay—that parameterize an earthquake’s complexity relative to the standard Brune model and strongly correlate with the estimate discrepancies. We find that the observed systematic magnitude–stress‐drop trends may reflect underlying changes in STF complexity, not necessarily trends in actual stress drop. Both the decay and BRE parameters vary systematically with magnitude, but whether this magnitude–complexity relationship is real remains unresolved.

Bulletin of the Seismological Society of America

The role of groundwater in contributing to surface water salinization in the Upper Colorado River Basin

Freshwater salinization impacts the availability of water for human use and ecosystem needs worldwide. It has been estimated that total dissolved solids (TDS) in the Colorado River Basin cause $350 million/year in damages and substantial resources are devoted to reducing TDS loading to streams. This study describes the development and application of coupled watershed models that enable TDS source tracking through the subsurface and across the landscape at a seasonal timestep for 35 years in the Upper Colorado River Basin. Results indicate that, on average, 75% of TDS loading to streams originates as baseflow, and 50% of loading is lagged in delivery by longer than one season. Snowmelt was identified as a dominant process controlling the transport of lagged TDS to streams. This approach informs when and where TDS mitigation efforts may be effective in a watershed that serves as a critical water supply for the southwestern United States.

Arizona, Colorado, New Mexico, Utah, Wyoming

Effects of snowmobile use on snowpack chemistry in Yellowstone National Park, 1998

Snowmobile use in Yellowstone National Park has increased substantially in the past three decades. In areas of greatest snowmobile use, elevated levels of by-products of gasoline combustion such as ammonium and benzene have been detected in snowpack samples. Annual snowpacks and snow-covered roadways trap deposition from local and regional atmospheric emissions. Snowpack samples representing most of the winter precipitation were collected at about the time of maximum annual snow accumulation at a variety of locations in the park to observe effects of a range of snowmobile traffic levels. Concentrations of organic and inorganic compounds in snow samples from pairs of sites located directly in and off snow-packed roadways used by snowmobiles were compared to concentrations in samples collected at nearby off-road sites. Concentrations of ammonium were 2 to 5 times higher for the in-road snow compared to off-road snow for each pair of sites. Thus, concentrations decreased rapidly with distance from roadways. In addition, concentrations of ammonium, nitrate, sulfate, benzene, and toluene in snow were positively correlated with snowmobile use.

Wyoming

Trimming the UCERF3-TD logic tree: Model order reduction for an earthquake rupture forecast considering loss exceedance

The Uniform California Earthquake Rupture Forecast version 3-Time Dependent depicts California’s seismic faults and their activity. Its logic tree has 5760 leaves. Considering 30 more model combinations related to ground motion produces 172,800 distinct models representing so-called epistemic uncertainties. To calculate risk to a portfolio of buildings, one also considers millions of earthquakes and spatially correlated ground-motion variability. We offer a tree-trimming technique that retains the probability distribution of portfolio loss and identifies the leading sources of uncertainty for further study. We applied it to a California statewide building portfolio and various levels of nonexceedance probability between one in 100 and one in 2500. We trimmed the logic tree from 172,800 leaves to as few as 15. The result: a supercomputer that would otherwise run 24 h to estimate the distribution of one-in-250-year loss can calculate it in moments with the reduced-order model. Others can use the reduced-order model to calculate risk to different California portfolios, and scientists can prioritize study to reduce the remaining epistemic uncertainty.

Earthquake Spectra

Effect of deicing chemicals on the hydrologic environment in Massachusetts; evaluation of surface-water data collection

The objective of this study is to develop predictive relationships that can be used to describe the impact of highway deicing salts on the hydrologic environment. A method, presented in an earlier report, for estimating yearly mean chloride concentrations from estimated or actual runoff and salt-application data was applied to the period 1971 through1974 and the results compared with yearly mean chloride concentrations computed from records of specific conductance. Estimates were found to be from minus 63 percent to plus 56 percent in error when the error was expressed as the difference between estimated and computed values as a percentage of the computed values. Chloride concentrations used in this study are computed from specific conductance/chloride concentration relationships and records of specific conductance. Preliminary results from graphic analyses and least-squares regression analysis show that relationships between measured values of specific conductance and chloride concentration have correlation coefficients ranging from 0.30 to 0.97. Analysis of streamflow for all major dissolved constituents is recommended with the purpose of attempting to describe the variations in the specific conductance/chloride concentration relationships.

Massachusetts

Estimating domestic self-supplied water use in Rhode Island, 2014–21

Water withdrawal from private groundwater wells is often unaccounted for in water planning studies, and water from private wells can be a source of exposure to environmental contaminants. The sizes of populations that depend on private wells for domestic water use and the amounts of water that are withdrawn from these wells are generally poorly represented in data collection efforts because of the challenges of locating, metering, or gathering withdrawal information from individual property owners. To address this problem, the U.S. Geological Survey, in cooperation with the Rhode Island Water Resources Board, estimated the volume of water withdrawn from domestic self-supply wells and the populations who use them for the State of Rhode Island at a 30-meter pixel spatial resolution and one-month temporal resolution between July 2014 and June 2021. The number of people reliant on domestic self-supply wells has increased in Rhode Island over the study period; however, the statewide estimate of total water withdrawal has not statistically increased. Withdrawals from private wells are largest in the inland areas of the western part of the State, and the towns of Scituate and Charlestown have the highest estimated withdrawals. Statewide monthly withdrawals ranged from 3.987 million gallons per day in March 2018 to 7.767 million gallons per day in September 2016. The median per capita domestic water use rate was 46.0 gallons per capita per day.

Rhode Island

Cruise summary—Samoa Basin abyssal mapping—Box coring leg, Pago Pago, Territory of American Samoa to Pago Pago, American Samoa, April 11– May 1, 2026

Expedition Summary U.S. Geological Survey scientists led a box coring effort to the Samoa Basin to characterize minerals and the surrounding abyssal sediments and fauna. Thirty-eight box cores were deployed between April 13, 2026, and April 28, 2026. Thirty-six box cores recovered sufficient material to determine nodule density, and 35 recovered sufficient material for subcores to be collected. The purpose of this Data Report is to provide a summary of samples collected, initial results that were obtained shipboard, and briefly mention planned future analyses from this expedition.

Data Report

Footprints of past mining in Alaska (USA) derived from high-resolution satellite imagery

Mapping the land area used for mining in the past is essential for guiding the remediation of affected landscapes and assessing the resource potential of related waste products. Despite significant recent progress delineating footprints of active and inactive mining globally, the known inventory of such mine lands remains incomplete. Here, I describe a new map dataset of footprints of land surface disturbance and waste at sites of past mining in Alaska (USA) based on visual interpretation of satellite imagery. This dataset maps 6–14 times the area of previous regional and global mine footprint maps in Alaska and is the first in the region to explicitly delineate mine waste landforms (e.g., tailings piles). The data are publicly available from the U.S. Geological Survey under a “no rights reserved” Creative Commons (CC0) license agreement.

Alaska

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

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

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

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