USGS Science⌕ Search

SEARCH · USGS Science

Search USGS Science

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.

At least 595 records · Page 33Linked to original sources
USGS2025

Chemostratigraphy of the Cretaceous Hue Shale in Arctic Alaska: Exploring paleoceanographic controls on trace element enrichment, organic matter accumulation, and source-rock evolution

We document chemostratigraphy in an outcrop of late Albian to early Campanian (∼103–82 Ma) marine source rocks to better understand paleoenvironmental controls on trace element (TE) enrichment and organic matter accumulation in the distal Colville foreland basin of Arctic Alaska and how those drivers are linked to arc volcanism and successions of Cretaceous oceanographic and climatic biogeochemical events. This unique, 113-m-thick section of Cretaceous Hue Shale deposited during a series of previously undocumented Arctic Cretaceous oceanic anoxic events (Lease et al., 2024) is the only known exposure of thermally immature (0.48–0.52% R o , random vitrinite reflectance) Hue Shale in Arctic Alaska. Strata comprise mainly clay-rich mudstone with elevated total organic carbon (TOC) and hydrogen index values reaching 26.3 wt% (mean = 7.5 wt%) and 689 mg hydrocarbon (HC)/g TOC (mean = 385 mg HC/g TOC), respectively. Maceral composition consists predominantly of fluorescent amorphous organic matter, with abundant brightly fluorescent alginite, including Tasmanites , acritarchs, and Leiosphaeridia . Discrete layers of volcanic ash (preserved as bentonite) are present throughout the section and provide quantitative age control based on U–Pb dates. Chemostratigraphic trends are investigated to advance our understanding of local oceanographic conditions and controls on stratigraphic and temporal heterogeneity of Brookian source rocks. Concurrent sedimentary enrichment in Mo, U, V, Pb, and Cu across the Albian–Cenomanian boundary of the exposed basal gamma-ray zone, may reflect anoxic to euxinic benthic redox conditions favoring organic matter accumulation and preservation. Fluctuating degrees of anoxia-euxinia are inferred throughout the overlying Hue Shale succession, reflected by varying patterns of TE enrichment and TE–TOC covariation. Elevated C org /P molar ratios (>250) across most of the section, with several values exceeding 690, signify that enhanced biological productivity is sustained throughout deposition. Enhanced productivity, recorded by both C org /P and excess Ba, also parallels increases in source rock richness (elevated TOC and S 2 values) during the late Albian–early Cenomanian and late Cenomanian – Turonian. Enhanced productivity and variations in oceanic circulation/stratification likely both drove changes in benthic redox conditions that favored organic carbon accumulation and preservation. Increased continental arc volcanism (e.g., Okhotsk–Chukotka volcanic belt) and High Arctic Large Igneous Province magmatic eruptions throughout the Cretaceous, inferred to have influenced nutrient cycling and local aqueous nutrient availability, also have been invoked as potential drivers of organic carbon burial and source-rock development across the sedimentary sequence. Results presented here document the organic-rich and oil-prone source-rock quality of the Hue Shale in the distal part of the Colville foreland basin and bolster the potential for a Cretaceous petroleum system beneath the eastern North Slope.

Alaska↗
USGS2025

Performance analysis of oil recovery and CO2 retention in a greenfield residual oil zone: CO2-EOR in Tall Cotton Field (Permian Basin, West Texas, USA)

Residual oil zones (ROZs) can offer significant oil resources via enhanced oil recovery (EOR) as well as subsurface carbon dioxide (CO 2 ) retention during injection. If injected CO 2 is anthropogenic, the ROZs can offer a substantial geologic storage potential. The ROZs below the oil/water contact (OWC) of main pay zones (MPZ) in conventional reservoirs or brownfields, are more commonly developed for CO 2 injection and oil production and reported in the literature. However, CO 2 -EOR in greenfield ROZs, reservoirs without a MPZ present, have rarely been developed for CO 2 -EOR operation. The Tall Cotton Field of West Texas, Permian Basin, which started production in 2015 (Phase 1) and expanded in 2017 (Phase 2) from the San Andres Limestone, is one of the first examples of greenfield ROZs developed for EOR by injecting CO 2 . This paper analyses EOR and CO 2 retention performance of Tall Cotton Field using allocated injection and production data from inverted 5-spot well patterns of Phase-1 and -2 developments. Production and injection data allocated to each of the 28 identified patterns (nine 20-acre patterns for Phase-1, three 20-acre and sixteen 10-acre patterns for Phase-2) were analyzed for historical and forecasted oil recovery using ratio-trend decline analysis, and for CO 2 retention performance of the patterns. The allocated data were further used to calculate injected reservoir pore volume and void replacement ratios (VRR) for the analysis period. Quantitative results indicated that oil recovery factors of the 5-spot patterns varied between 4–10 %, and 5–30 % between the end of injection and the forecast periods, respectively. Storage of CO 2 , on the other hand, increased to a mean value of ∼7130 MMscf per pattern in Phase-1 and to a mean storage of 3700 MMscf per pattern in Phase-2 until the end of injection, followed by a decline after the end of injection and into the forecast period. Resulting CO 2 utilization factors ∼6–50 Mscf/bbl were estimated at the end of injection. Overall, presented results suggested that developing greenfield ROZs for CO 2 -EOR can be as promising as brownfield ROZs and mature MPZs for EOR and underground storage of injected CO 2 . For Tall Cotton Field, results suggest that Phase-2 patterns generally outperformed Phase-1 for oil recovery factors, while Phase-1 performed better in CO 2 retention performance metrics. This is the first study in the literature that reports a detailed CO 2 -EOR performance analysis of a greenfield ROZ in the Permian Basin, which can potentially allow for comparison with MPZs and brownfield ROZs.

Texas↗
USGS2025

Characterizing sedimentary organic carbon in a hydrothermal spreading center, the Escanaba Trough

Sediments in critical marine mineral environments are of wide importance due to their preservation of both marine minerals and organic carbon (OC) stocks. However, OC storage and cycling is often overlooked in mineral system studies. This work characterizes sedimentary OC within the Escanaba Trough, a hydrothermal sulfide system off the coast of northern California. By utilizing ROV-based push coring, we collected sediments near and far from high temperature, low temperature, and inactive vents. We applied a multipronged organic geochemical approach, measuring bulk sediment, OC fractions of varying labilities, and biomarkers to tease apart the storage, source, and cycling of OC within this complex system. In contrast to past work indicating a primarily terrestrial source to deeper, Pleistocene sediments, our results suggest a primarily marine source in surface sediments at least 50 m away from active venting. Near active venting, we see evidence of locally produced chemosynthetic OC in addition to the background marine source. This chemosynthetic OC appears to be rapidly remineralized and supports more complex deep-sea food webs compared to hemipelagic sites. Still, the greatest contribution of labile OC was observed at inactive vent sites, which we suggest is due to hydrothermal fluid alteration of more labile OC at actively venting sites. Depending on fluid temperature, hydrothermally altered OC is either preserved in the sediments as petroleum or migrated with fluid flow. Considering inactive sites are those most likely to be targeted by potential mining, we suggest additional studies could verify if these results apply at other sedimented seafloor massive sulfide systems.

California↗
USGS2025

The impact of burial diagenesis on soil-formed minerals in paleosols using stable isotopes of phyllosilicates and carbonate clumped isotopes

To understand the effects of burial diagenesis on the stable isotope geochemistry of soil-formed clay and carbonate minerals in paleosols, samples were collected from seven cores, spanning middle- to upper-Pennsylvanian strata of the Illinois Basin, with varied maximum burial depths of 1–3 km. Mixed-layer illite-smectite and kaolinite mixtures give δ 2 H and δ 18 O values of −83 ‰ to −36 ‰ and 11.9 ‰ to 21.1 ‰ (VSMOW), respectively. After carbonates were screened petrographically for diagenetic textures using transmitted light and cathodoluminescence, measured clumped isotope Δ 47 values range from 0.504 to 0.563 ‰ (I-CDES). Resulting mineral formation temperatures for phyllosilicate mineral mixtures are 28 to 66 °C (mean = 47 °C), whereas T(Δ 47 ) estimates for calcites are 36 to 61 °C (mean = 45 °C). Calculated δ 18 O water values from which phyllosilicate minerals and calcites precipitated under isotopic equilibrium ranges from −7.1 to −1.2 ‰ and − 1.4 to +4.9 ‰, respectively. Closed and open-system phyllosilicate-fluid exchange modeling indicates that phyllosilicate alteration occurred in the presence of a low temperature brine or meteoric water and is interpreted to occur in a layer-by-layer illitization transformation. Due to the lack of diagenetic textures and positively correlated T(Δ 47 ) and δ 18 O water , calcites are interpreted to have undergone solid-state bond reordering. Despite low to moderate temperatures (<125 °C) and varying depths of shallow burial (1–3 km), solid-state transformation of phyllosilicates and calcites indicates paleosols had prolonged exposure to burial conditions which has implications for the use of paleosol minerals for paleoenvironmental reconstructions.

Illinois, Indiana, Kentucky↗
USGS2025

Confocal laser-scanning microscopy (CLSM)-based thermal maturity of Tasmanites and progress in standardization of fluorescence microspectrometry

Evaluation of thermal maturity in vitrinite-free or vitrinite-deficient sediments via fluorescence microspectrometry can provide relevant information related to petroleum exploration and thermal history assessment. However, variation in spectral fluorescence properties of alginite macerals with increasing thermal maturity is largely underexplored. Here, authors of this study have applied confocal laser-scanning microscopy (CLSM) in conjunction with fluorescence microspectrometry to a maturity series of marine Upper Devonian Tasmanites algae from the Ohio Shale (Huron Member) and a single sample from the Marcellus Formation of the Appalachian Basin. Spectral fluorescence properties of Tasmanites were evaluated in relation to orientation, measurement location, and the number of measurements per sample, and were compared to published literature. Emission spectra of Tasmanites from continuous wave laser excitation (405 nm) were acquired from sections perpendicular and parallel to bedding and at different positions within individual Tasmanites bodies. The results showed a progressive red-shift in emission maxima (λ max ) in a large sample sized maturity series ( N = 19), e.g., 493 to 578 nm for the perpendicular section at middle position. Further, blue-shifted apex and mineral-adjacent positions within sections perpendicular to bedding were observed, with the latter being reported here for the first time. While blue-shift at apex positions can be attributed to mechanical deformation-induced reorientation of photoselected fluorophores, the blue-shifted mineral-adjacent positions could result from strain loading and development of a plastic deformation region at the mineral contact zone with Tasmanites . A decrease in standard deviation with increasing number of measured emission maxima is well-observed, and 15 to 20 individual measurements per sample appears sufficient for low standard deviation and coefficient of variance. CLSM-derived thermal maturity parameters indicated that a moderate positive correlation of red/green quotient (Q 650/500 ; R 2 = 0.67) with solid bitumen reflectance (BR o in %) exists. For reproducible results, the determination of λ max and Q 650/500 should be conducted exclusively in the middle position at perpendicular and parallel sections of the polished whole-rock pellets, where the lowest standard deviation in measurement was observed. These results strengthen the suitability and relevance of the CLSM technique in thermal maturity studies of dispersed organic matter (DOM) and contribute to the standardization of fluorescence microspectrometry methods in organic petrology investigation.

Kentucky, Ohio, Virginia, West Virginia↗
USGS2025

Refining the earthquake history of south-central Alaska through lake records

The Alaska–Aleutian subduction zone (AASZ) is one of the world's most seismically active plate boundaries and the source of the 1964 M w 9.2 Great Alaska earthquake–the second largest instrumentally recorded earthquake in the world. Understanding the nature and frequency of such earthquakes is necessary for seismic and tsunami hazard assessment, but instrumental and historical records that span less than 150 years are too short to allow a statistically reliable analysis of earthquake recurrence times. This calls for studies of evidence of past earthquakes, extending the earthquake catalog further back in time. Subduction-zone paleoseismology in south-central Alaska is predominantly based on coastal evidence of land-level changes and tsunamis generated by megathrust earthquakes and preserved in the geological record. A complementary approach is lacustrine paleoseismology, which is still a relatively young discipline in Alaska. However, globally, lake basins are well-established high-resolution and continuous recorders of paleoseismic activity along subduction zones, relying on the identification of underwater landslide deposits and turbidites generated by seismic shaking. As a result, lake basins not only register ground shaking from megathrust earthquakes, but also from intraslab and crustal earthquakes, which are typically not accompanied by significant land-level changes. In this review paper, we combine coastal and lacustrine paleoseismology approaches to refine the south-central Alaskan earthquake history by comparing the paleoseismic records from two lakes (i.e., Eklutna Lake, located in the Chugach Mountain Range, and Skilak Lake, situated on the Kenai Peninsula) with the coastal and crustal earthquake catalog in Alaska. The resulting age ranges of all known megathrust earthquakes involving the Alaskan megathrust between the Kodiak and Prince William Sound (PWS) sections are more precise and accurate for the last 1.3 kyrs BP than the previously published age ranges from coastal records. As a result, this study supports the following key conclusions: (1) The 1964 CE earthquake was an exceptionally strong and unique event in the last 2000 years, rupturing the PWS, Kenai, Barren Islands, and Kodiak sections simultaneously. (2) The high-resolution and seasonal markings of the varved lake records now disentangle for the first time closely timed earthquakes, which was not possible based on the coastal evidence alone. (3) No persistent megathrust rupture boundaries exist. So, the possibility of a full rupture of the entire eastern AASZ, from PWS to Semidi cannot be excluded. (4) The rupture pattern in the eastern AASZ reveals superimposed cycles of multi-asperity ruptures (1964 earthquake) and clustered complementary partial ruptures, or rupture cascades. (5) The PWS section hosts the largest asperity in the eastern AASZ. (6) The shaking record of megathrust earthquakes indicates a time-dependent (quasiperiodic) behavior for the study area, but the observation of complementary clusters means that the hazards will not drop to zero but instead may even increase for a neighboring section. (7) The time-independent behavior of intraplate earthquakes implies that the intraslab hazard did not decrease following the 2016 and 2018 earthquakes. This study utilizes an integrated approach for subduction zone paleoseismology as a solution for unraveling recurrence and rupture patterns in Alaska, which can be applied worldwide.

Alaska↗
USGS2025

Simulation of the impacts of projected climate change on groundwater resources in the urban, semiarid Yucaipa Valley watershed, southern California using an integrated hydrologic model

Managing water resources in semiarid watersheds is challenging due to limited supply and uncertain future climate conditions. This paper examines the impact of future climate changes on an urban watershed in southern California using an integrated hydrologic model. GSFLOW modeling software is used to simulate the nonlinear relationships between climate trends and precipitation partitioning into ET, runoff, and subsurface storage. Four global circulation models (GCMs), each with two greenhouse-gas scenarios, RCP45 and RCP85 are used to project future climate conditions. GCMs include the CanESM2, CNRM-CM5, HadGEM2-ES, and MIROC5 models. The model's simulated hydrologic conditions are compared with historical data to assess changes in water budgets and groundwater supply. Results indicate decreased groundwater storage in most scenarios due to increased natural evapotranspiration, vegetation consumptive use, and streamflow out of the watershed. Only scenarios with substantially increased future precipitation show increased groundwater storage. The study also highlights increased future aridity despite the rise in precipitation and large precipitation events forecast by GCMs, which increase the risk of urban floods and decrease stream leakage and water available to vegetation.

California↗
USGS2025

Hydrologic response of groundwater and streamflow to natural and anthropogenic drivers of change in headwaters of the upper Colorado River basin during recent wet (1982–1999) and drought (2000–2022) conditions

Study region: Headwaters of the upper Colorado River basin (UCOL), USA Study focus: Surface-water and groundwater numerical models incorporating water-use information were used to investigate changes in climate, water use, and simulated hydrologic responses of snow processes, evapotranspiration, groundwater, and streamflow during recent wet (1982–1999) and drought (2000–2022) periods in the headwater subregions of the upper Colorado River basin. New hydrologic insights for the region: Decreases in average streamflow between wet and drought periods ranged from 20 % in the Colorado River headwaters subregion to 23 % in the Gunnison River headwaters subregion. Like streamflow, average surface runoff was statistically less during the drought than the wet period, with decreases from 24–31 % in the headwaters. On a volume basis, runoff decreases were greater than streamflow decreases in both the Colorado River and Gunnison River headwaters. Although the amount of water-year groundwater discharge to streams remained nearly the same between the wet and drought periods, groundwater as a percentage of streamflow increased between the wet and drought periods, highlighting the importance of groundwater in sustaining streamflow during drought conditions. Multiple linear regression analyses revealed that snowmelt-only models were better than the best precipitation and temperature models at explaining streamflow variability from all headwater subregions for both the wet and drought periods.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2025

Linking stream-reach nitrogen loads and groundwater “reachsheds” to inform wastewater-nitrogen management actions, Cape Cod, Massachusetts

Study Region Cape Cod, Massachusetts, U.S.A. Study Focus Anthropogenic nitrogen (N) is a key factor in degrading groundwater and surface-water quality, particularly in coastal New England where onsite wastewater systems are prevalent. This study evaluated whether direct N-load measurements in streams on Cape Cod, Massachusetts, coupled with flow-path information from groundwater-flow models, can effectively identify potential land areas where nitrogen mitigation could substantially reduce loads to receiving waters. Nitrogen fluxes were measured along 63 stream reaches during winter and summer and paired with simulated groundwater recharge areas to identify and rank potential areas for reduction of nonpoint-source N inputs. New Hydrologic Insights for the Region Reach-scale nitrate-N loads ranged from −39.1–1182 kg-N/yr per 100 m of stream, indicating spatially variable groundwater inputs across seasons. “Reachsheds” — areas contributing groundwater recharge to specific stream reaches — were delineated using a regional groundwater-flow model. Strong correlations were found between observed N loads and land-use characteristics, especially the number of septic systems and total N inputs from the sum of considered sources. Observed N loads were moderately correlated with recharge area size and wastewater flow estimates. Correlating reach-specific groundwater N loads with land use and parcel-scale nitrogen-yield data identified reachsheds with the highest potential for N load reduction. This approach enables targeted implementation of restoration efforts to optimize nutrient management and support regional load reduction.

Massachusetts↗
USGS2025

Systematic approach to prioritize wells for effective groundwater monitoring and management in the Arkansas Headwaters Basin, Colorado, USA

Study region The Arkansas Headwaters Basin, an intermountain basin in the Southern Rocky Mountains of North America. Study focus Our specific focus is choosing a set of wells to support a possible future regional groundwater-surface water model that would support water management. We present a three-step process using multiple criteria to score, predict, and choose prioritized wells that capture the full distribution of data including extremes. The three-step process provides accessible visualizations, fiscally efficient well prioritization, and screening useful for subsequent groundwater modeling. The novelty of the proposed methodology is the systematic approach integrating a scoring and a predictive approach to support a selection path. The systematic approach may be broadly adapted for other basins. New hydrological insights for the region Understanding regional hydrology hinges on efficient collection of hydrologic data that captures the relevant dynamics including extremes. The present study, a case study for a particular basin in the Southern Rocky Mountains, is the first use of a scripted (R software) strategy to select an economical and representative set of monitoring wells. Our findings suggest caution when using proximity as a proxy for correlation, because proximal wells in the same geologic formation and similar depths are not always correlated. In the Arkansas Headwaters Basin, subsurface geology may be less influential on groundwater elevations than broader hydrologic influences, such as regional drought.

Colorado↗
USGS2025

Statewide cumulative human health risk assessment of inorganics-contaminated groundwater wells, Montana, USA

Across the United States, rural residents rely on unregulated and generally unmonitored private wells for drinking water, which may pose serious health risks due to unrecognized contaminants. We assessed the nature, degree, and spatial distribution of cumulative health risks from inorganic contaminants in groundwater. Our analysis included nearly 84,000 data points from 6500+ wells, across 51 of Montana's 98 watersheds, using a public groundwater database. We compared a drinking water screening level cumulative risk assessment (CRA) for inorganics based on the U.S. Environmental Protection Agency (EPA) protective health thresholds (Maximum Contaminant Level Goals, Health Advisories [MCLG-HAs]) to a CRA based on EPA public supply enforceable Maximum Contaminant Levels (MCLs). Based on median concentrations of 19 inorganics (antimony, arsenic, barium, beryllium, boron, cadmium, chromium, copper, fluoride, manganese, molybdenum, nickel, nitrate, lead, selenium, strontium, thallium, uranium, zinc), 75% of watersheds had MCLG-HA-based cumulative risk values > 1.0; arsenic and uranium contributed the most risk, followed by strontium, fluoride, manganese and boron. Hence, this screening level (Tier I) CRA indicated widespread potential for unrecognized human health risk to private well users from inorganic contaminants considering both carcinogenic and non-carcinogenic risks. Sensitivity analysis showed that benchmarks applied (MCLG-HAs versus MCLs) exerted the largest control on results. Our findings identify priority regions for Tier 2 risk assessments to elucidate local sources and distributions of geogenic versus anthropomorphic contaminants. Our study is the first statewide assessment of cumulative health risk from groundwater that we are aware of, and results support increased statewide drinking water education and testing to reduce human health risks from contaminated private well water.

Montana↗
USGS2025

Evaluation of daily stream temperature predictions (1979-2021) across the contiguous United States using a spatiotemporal aware machine learning algorithm

Stream temperature controls a variety of physical and biological processes that affect ecosystems, human health, and economic activities. We used 42 years (1979–2021) of data to predict daily summary statistics of stream temperature across >50,000 stream reaches in the contiguous United States using a recurrent graph convolution network. We comprehensively documented the performance – both across all reaches and by stream type (e.g., reservoir or groundwater influence) – as a baseline for future improvement. The model showed reach-level RMSE of <2 °C with 90 % prediction intervals that contain 90.7 % of observations. We also assessed how the model captured variability in ecologically relevant metrics (e.g., R 2 for annual 7-day maximum = 0.76; R 2 for days exceeding 25 °C = 0.75). This model does not outperform state-of-the-art machine learning efforts (e.g., RMSE ≤1.5 °C) due to a limited input set but does provide the most spatially complete modeling to date to support water availability assessments.

contiguous United States↗
USGS2025

Probing lower-crustal fault properties with frequency-dependent tidal tremor triggering

The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.

California↗
USGS2025

Smectite-rich horizons in Inceptisols trigger shallow landslides in tropical granitic terranes

Puerto Rico was affected by >70,000 landslides in the wake of 2017 Hurricane Maria, and landslide prevalence was especially high in the Utuado region in the Cordillera Central. Landslide density was highest where soil parent material is granodiorite; landslide slip surfaces tended to be shallow (<60 cm), and often were mobilized rapidly and with long runout distances. This study combines field observations with soil mineralogy (bulk and clay fractions), soil geochemistry (bulk fraction), and soil strength as determined by field cone penetrometer testing (CPT) to test the hypothesis that clay-rich subsoil horizons function as slip planes when water-logged. Soil pits were excavated to depths of ∼200 cm in Ultisols on an undulating plateau and to ∼100 cm in Inceptisols on steep slopes (36-43 o ) that flank the plateau and cone penetrometer tests (CPT) were done within 2 m of the soil pit. Six pits were located adjacent to scarps from previous landslides, enabling analysis of soil profiles downward through extrapolated slip surfaces. Results from X-ray diffraction (XRD), X-ray fluorescence (XRF) and thermogravimetric analysis (TGA) indicate that soils are heterogeneous, often with subsoil horizons enriched in clay minerals and immobile elements (Al, Fe, Ti). Inceptisols on steep slopes often contain smectite-rich horizons at 30–60 cm depth that appear to function as slip surfaces; in other Inceptisols, such horizons are not present and landslide susceptibility is potentially lower. In Ultisols, soil mineralogy is dominated by kaolinite with minor halloysite, and increased kaolinite content at soil depths ≥80 cm at some sites suggests potential slip surfaces enhancing probability of landslides. The origin of clay-rich horizons appears to be (1) fractures in granodiorite that facilitate water flow and leaching, accelerating mineral dissolution during early weathering stages, and (2) smectite-rich buried soils under permeable colluvium likely deposited by a prior mass wasting event. Where clay-rich layers occur beneath more-permeable horizons, rapid infiltration then absorption of water in clay-rich subsoil horizons causes decreased shear strength and increased landslide susceptibility.

Puerto Rico↗
USGS2025

Reservoir thermal energy storage pre-assessment for the United States

Storing thermal energy underground for later use in electricity production or direct-use heating/cooling is a promising, viable, and economical green energy option. Reservoir thermal energy storage (RTES) is one such option, which stores energy in underutilized permeable strata with low ambient groundwater flow rates and more geochemically evolved (e.g. brackish/saline) waters relative to overlying principal aquifer systems. The U.S. Geological Survey has begun assessing RTES potential nationally by focusing on five generalized geologic regions (Basin and Range, Coastal Plain, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Hydrogeologic reservoir models are developed for the following eight metropolitan area cities within those regions to evaluate RTES performance across different climates and subsurface conditions: Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon. Evaluated metrics include estimated required well spacing, thermal storage capacity, and thermal recovery efficiency through time. Also considered for each reservoir are potential complicating factors, including reservoir depth, thermally driven free convection, and groundwater salinity. This work focuses on direct-use cooling because the need for cooling modern office buildings greatly exceeds that for heating in most parts of the country (Falta and others, 2016); however, the evaluated metrics are also relevant to heating and electricity applications. Results indicate that favorable RTES conditions exist in each region, with the Coastal Plain and Basin and Range being especially favorable for thermal storage capacity, while the Pacific Northwest and Michigan Basin excel at energy recovery for the evaluated cooling application. The results underscore the utility of developing maps of thermal storage capacity, subsurface temperature models, and volumetric estimates of thermal storage capacity to serve as key RTES resource classification standards. Overall, this pre-assessment provides a basic understanding of RTES potential in several cities and geologic regions throughout the country and will aid ongoing thermal energy storage assessment efforts.

Arizona, Illinois, Michigan, New Mexico, Oregon, S↗
USGS2025

Favorability mapping for hydrothermal power resource assessments of the Great Basin, USA

The U.S. Geological Survey (USGS) is updating the 2008 assessment of conventional hydrothermal resources for the Great Basin in the western United States. As part of this work, the workflow for hydrothermal resource favorability maps is being modified to integrate modern data-driven machine learning (ML) methods. Improvements include: [1] using new and refined evidence layers (features); [2] using an order of magnitude more training sites (labeled examples); [3] utilizing simple but non-linear supervised ML algorithms; [4] representing positive training sites (wells with measured heat flow) with their ordinal value proportional to the magnitude of convective upflow (i.e., low, high, or very high convective signals instead of past strategies using positive-negative labels); [5] supplementing training sites with additional sites with low convective signals to represent diverse under-sampled areas where hydrothermal systems are unlikely to exist; [6] comparing with competing approaches; and [7] utilizing Monte Carlo cross-validation to estimate and evaluate prediction uncertainty. For the new favorability map, over half of the power-producing systems (i.e., 15 of 28) are predicted in the 99th percentile of most favorable locations (i.e., the highest 1 % of favorability, corresponding to 1 % of the map area), exceeding the performance of past models that have explicitly used power plants as training sites. Previous favorability maps predicted approximately half of the power-producing hydrothermal systems above the 80th percentile (i.e., 20 % of the map area). For the new favorability map, 93 % of power-producing systems (i.e., 26 of 28) are above the 80th percentile. The power-producing systems for which the new model does not perform well are either comparatively small, low-temperature systems or systems also not predicted well by prior modeling approaches, suggesting that these few systems are unusual when compared with most power-producing systems. Focusing research on these known, seemingly different systems may yield new insights and subsequent discovery of new prospects.

California, Idaho, Nevada, Oregon, Utah↗
USGS2025

Simulating human behavior under earthquake early warning

Earthquakes are a rapid-onset hazard where advance planning and learning plays a key role in mitigating injuries and death to individuals. Recent advances in earthquake detection have resulted in the development of earthquake early warning (EEW) systems. These systems can provide advance warning to predetermined geographic regions that an earthquake is in progress, which may result in individuals receiving warning seconds before significant shaking is felt at their location. This additional time could allow individuals to take more effective protective actions during the immediate disaster. To demonstrate this, we created an agent-based simulation of a basic apartment that allowed us to randomly and repeatedly simulate an individual receiving and responding to an EEW message. The results of our preliminary simulation show that, in our study environment, earthquake early warning alerts have the potential to allow for sufficient time for individuals to take protective actions.

Washington↗
USGS2025

Estimated average annualized losses from potential building damage and fatalities due to earthquake-generated tsunamis in the United States

Earthquake-generated tsunamis represent substantial economic threats to states and territories in the United States (U.S.), but we are unaware of any effort to quantify potential impacts at the national level. This gap is partially due to the lack of nationally consistent data on tsunamigenic sources and associated return periods. This study addresses this issue and provides estimates of average annualized losses (AAL) for potential residential fatalities and capital stock losses associated with building damage (i.e., structural, non-structural, contents, and inventory damage) in the U.S. by curating tsunami-hazard information based on deterministic scenarios and probabilistic approaches, calculating potential losses, and estimating return periods where necessary. This assessment was done for the U.S. West Coast, Alaska, Hawaii, U.S. Pacific Territories, and U.S. Atlantic Territories. We estimate that earthquake-generated tsunamis that could affect these states and territories collectively represent $1 billion in potential AAL with 79 % of losses due to residential fatalities and 21 % of losses due to capital stock losses from building damage. We identify AAL variations based on county and county equivalents, states and territories, geographic regions, return periods, and departure-delay assumptions for evacuating residents. Results include high AAL values for potential fatalities in Puerto Rico and the U.S. Pacific Northwest region, high AAL values for potential building-related damage in Hawaii and California, and high building- and population-loss ratios for county equivalents in Alaska and U.S. territories.

International Journal of Disaster Risk Reduction↗