USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers for Young Minds”

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,406 records · Page 30Linked to original sources
USGS2024

Hydroacoustic observations reveal drivers of mixing and salinization of a karst subterranean estuary during intense precipitation

Karst subterranean estuaries within globally ubiquitous carbonate aquifers are coastal groundwater ecosystems that provide an essential water resource for human populations. To understand the drivers of salinization within a coastal aquifer in the Yucatan Peninsula (Mexico), we employed hydroacoustics in flooded caves to observe how oceanic and atmospheric events facilitate mixing between the meteoric lens (fresh-brackish groundwater) and the saline groundwater on tidal and episodic timescales. Precipitation during Tropical Storm Carlotta increased the flow and salinity of the meteoric lens without evidence for vertical mixing across the halocline. We postulate that vertical migration of haloclines in the conduit relative to those within the rock matrix during precipitation creates lateral density gradients that drive mixing, and ultimately creates a brackish layer within the meteoric lens. These results provide a mechanistic explanation for vertical and lateral exchange in a coastal carbonate aquifer, which has implications for groundwater response to future climatic change.

Ox Bel Ha cave system, Yucatan peninsula↗
USGS2026

Storage timescales and the crystal growth history recorded in the layered reservoir of the Unit 3 Kulanaokuaiki Tephra Member of the Uēkahuna Ash from Kīlauea volcano

The Kulanaokuaiki Tephra Unit 3 (900 C.E.) was a subplinian basaltic scoria eruption, the largest known to have occurred at Kīlauea’s summit. The initiation mechanism and cause(s) of this highly hazardous eruption remain poorly understood, particularly as this style of activity has not been observed in the modern period. We use olivine diffusion timescales from Fe-Mg and lithium zoning to assess magma storage timescales and evaluate evidence for mixing prior to eruption. Of the 200 crystals assessed, 55% are normally zoned for forsterite and 45% are non-zoned. Zoned crystals (> 1 mm) ( n = 96) have average compositions of Fo 89 cores and Fo 81 rims, while small crystals 0.5–1 mm ( n = 14) have Fo 85 cores and Fo 82 rims. Modeled Fo diffusion timescales range from < 1 to 12.6 years (median of 3.0 years). Lithium maps in 32 of the forsterite-zoned crystals reveal a variety of lithium zoning types, most of which we attribute to crystal growth. Many lithium profiles are coupled to phosphorus within the outer ~ 100 μm of the crystal ( n = 17), reflecting growth processes and precluding their use for diffusion studies. Overall, our Unit 3 olivine dataset displays no evidence of magmatic priming shortly prior to eruption, unlike prior observations for the basaltic subplinian eruption of the Keanakāko‘i Tephra Unit E (1650 C.E.). Thus, explosive eruption magmas may be stored for considerably longer than previously thought, and magmatic priming is not a necessary precursor for Kīlauea’s explosive eruptions. We hypothesize that gas-driven filter pressing from crystallization of microlites and macrocrysts in the reservoir slurry contributed to the explosivity of the Unit 3 eruption.

Hawaii↗
USGS2025

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming↗
USGS2025

Factors affecting the distribution of water-bearing fractures in the bedrock aquifers of West Virginia

Bedrock aquifers cover 23,601 square miles within the State of West Virginia and comprise 97.4 percent of the surficial area within the State; the remaining 2.6 percent (621 square miles) consists of alluvial sand-and-gravel and glacial outwash aquifers bordering the State’s major rivers. While West Virginia’s alluvial aquifers have been studied extensively, bedrock aquifers have only been characterized for studies completed in a few areas in Jefferson, McDowell, and Monroe Counties. Bedrock aquifers are water supplies for public supply, agriculture, industry, and residential homeowner use. In this study, the U.S. Geological Survey, in cooperation with the West Virginia Department of Environmental Protection Division of Water and Waste Management, provides a statewide assessment of the occurrence and distribution of fractures within bedrock aquifers of the State and the various topographic, physiographic, and lithologic influences controlling the occurrence and distribution of bedrock fractures. The results of this study provide an increased understanding of the distribution of fractures in bedrock aquifers in West Virginia and help to verify trends that have been suspected for many years but were never well documented or verified by data. The results confirmed that the density of fractures and those that were determined to be water bearing decrease significantly with depth. A statistically significant difference in the density of fractures was observed at a depth of 215 feet for wells in the Appalachian Plateaus Physiographic Province’s and in the Valley and Ridge Physiographic Province’s aquifers; a higher density of fractures and water-bearing fractures were above a depth of 215 feet than below that depth. This is an important consideration when drilling wells for residential, commercial, industrial, or agricultural water supply. Abandoned underground coal mines are commonly believed to form large pools of water in the interconnected mine entries in abandoned room and pillar coal mines. Such pools of water can and do exist in abandoned underground coal mines, but many mines lack open entries and are held up by overburden strata and pillars that can collapse and form aquifers comprised of vast interconnected rubble zones (gob), especially in older mines. Data assessed for this study showed that shale-corrected values of effective porosity for limestone aquifers in West Virginia had a median value of 2 percent and an average value of 4 percent and generally are mineralized with low effective porosity. Argillaceous or sandy limestone has a median shale-corrected porosity of 4 percent and an average shale-corrected porosity of 5 percent. The median and average shale-corrected porosity of sandstone aquifers was estimated to be 14 percent, but the median shale-corrected porosity for argillaceous or calcareous sandstone was 5 percent and the average shale-corrected porosity for argillaceous or calcareous sandstone was 6 percent. Even though shale has a relatively high total sonic porosity compared to other lithologies, shale and siltstone had relatively low shale-corrected porosity, ranging from 0 to 2 percent. Well yields were previously documented to be highest in valley settings, lowest on hilltops, and intermediate on hillsides. Transmissivity data provided by this study confirm this general pattern within the Appalachian Plateaus Province; however, the Valley and Ridge Province does not follow this pattern. While still lowest on hilltop settings, the highest well yields were in hillside settings. The trend for the Valley and Ridge Province was likely skewed because of 9 high-yield wells specifically targeting deeper thin limestone units, such as the Tonoloway and Helderberg Limestones, at depths with transmissivity in excess of 2,000 feet squared per day in Mineral County, West Virginia, or targeting karst aquifers in Berkeley, Jefferson, or Greenbrier Counties, West Virginia. Finally, water-bearing fractures have been hypothesized to comprise a small number of all fractures within a typical bedrock well in West Virginia. Data collected for this study support this theory. A total of 3,403 fractures were identified during this study; 3,151 (92.6 percent) of those fractures are low-transmissive, and only 252 (7.4 percent) fractures are water-bearing. Even though a well may contain many fractures, less than 8 percent are considered water-bearing fractures.

West Virginia↗
USGS1996

Magnitude and frequency of floods in the Suwannee River Water Management District, Florida

Flood-frequency statistics for 2-, 5-, 10-, 25-, 50-. 100-, 200-, and 500-year recurrence intervals, based on three methods of analysis, are presented for 25 continuous-record and seven peak flow partial-record gaging stations in the Suwannee River Water Management District. The first method, for gaged stations, utilizes station records; the second method, for ungaged sites, utilizes regional regression analysis; and the third method uses a weighted combination of the station and regional values. Because the weighted values utilize two more or less independent estimates of the peak flow statistic, they are considered more accurate than the station estimates or the regression estimates alone. Also, the use of another weighting scheme to improve estimates of flood frequency statistics at ungaged sites is demonstrated. The karstic nature of much of the Suwannee River Water Management District significantly attenuates flood peaks in some streams by providing substantial subsurface storage when river stages are high. At such times, springs discharging into rivers may reverse flow temporarily and become sinks.

Florida↗
USGS2024

Elastic stress coupling between supraglacial lakes

Supraglacial lakes have been observed to drain within hours of each other, leading to the hypothesis that stress transmission following one drainage may be sufficient to induce hydro-fracture-driven drainages of other nearby lakes. However, available observations characterizing drainage-induced stress perturbations have been insufficient to evaluate this hypothesis. Here, we use ice-sheet surface-displacement observations from a dense global positioning system array deployed in the Greenland Ice Sheet ablation zone to investigate elastic stress transmission between three neighboring supraglacial lake basins. We find that drainage of a central lake can place neighboring basins in either tensional or compressional stress relative to their hydro-fracture scarp orientations, either promoting or inhibiting hydro-fracture initiation beneath those lakes. For two lakes located within our array that drain close in time, we identify tensional surface stresses caused by ice-sheet uplift due to basal-cavity opening as the physical explanation for these lakes' temporally clustered hydro-fracture-driven drainages and frequent triggering behavior. However, lake-drainage-induced stresses in the up-flowline direction remain low beyond the margins of the drained lakes. This short stress-coupling length scale is consistent with idealized lake-drainage scenarios for a range of lake volumes and ice-sheet thicknesses. Thus, on elastic timescales, our observations and idealized-model results support a stress-transmission hypothesis for inducing hydro-fracture-driven drainage of lakes located within the region of basal cavity opening produced by the initial drainage, but refute this hypothesis for distal lakes.

JGR Earth Surface↗
USGS2024

Temperature

The “National Field Manual for the Collection of Water-Quality Data” (NFM) provides guidelines and procedures for U.S. Geological Survey (USGS) personnel who collect data used to assess the quality of the Nation’s surface-water and groundwater resources. This chapter, NFM A6.1, provides guidance and protocols for the measurement of temperature of air, of a surface-water body or in groundwater, which include the scientific basis of the measurement, selection and maintenance of equipment, calibration verification, troubleshooting, and procedures for measurement and reporting. It updates and supersedes USGS Techniques of Water-Resources Investigations, book 9, chapter A6.1, version 2.0, by Franceska D. Wilde. Temperature of air and water is routinely measured when water samples are collected, is often measured continually at USGS streamgages, and is a parameter regularly measured during laboratory and field experiments. The field method for measuring temperature described in this chapter is applicable to air and most natural waters. Before 2017, the NFM chapters were released in the USGS Techniques of Water-Resources Investigations series. Effective in 2018, new and revised NFM chapters are being released in the USGS Techniques and Methods series; this series change does not affect the content and format of the NFM. More information is in the general introduction to the NFM (USGS Techniques and Methods, book 9, chapter A0) at https://doi.org/10.3133/tm9A0 . The authoritative current versions of NFM chapters are available in the USGS Publications Warehouse at https://pubs.usgs.gov/ . Comments, questions, and suggestions related to the NFM can be addressed to nfm@usgs.gov .

Techniques and Methods↗
USGS2024

Three-dimensional hydrogeologic framework model of the Russian River watershed, California

The Russian River watershed is in northern Sonoma County and southern Mendocino County, California, in the northern part of the California Coast Ranges. The Russian River serves as a supply for agricultural irrigation and for municipal, domestic, and commercial uses. Through a cooperative agreement with the California State Water Resources Control Board and Sonoma County Water Agency, the U.S. Geological Survey has completed studies to better understand the hydrogeologic system and develop numerical hydrologic modeling tools to evaluate and aid in managing groundwater resources. This report focuses on the development of a digital three-dimensional hydrogeologic framework model of the Russian River watershed for use in groundwater resource assessment and numerical models. The digital three-dimensional hydrogeologic framework model of the Russian River watershed portrays the altitude, thickness, and extent of five hydrogeologic units. These five hydrogeologic units include (1) a basement unit, (2) the Sonoma Volcanics, (3) a consolidated sedimentary rock unit, (4) an unconsolidated sediment unit, and (5) channel alluvium. Model input data were compiled from published geologic maps, interpreted well data, and a model of the top of basement derived from gravity data. These data were used to construct surfaces that represent the upper and lower subsurface boundaries of each hydrogeologic unit. Top surfaces were created for the five hydrogeologic units and then stacked in three dimensions to create a solid-volume digital model. The digital three-dimensional hydrogeologic framework model described in this report and the corresponding data represent the generalized geometry of the subsurface geologic units; the model reproduces the input geologic data with reasonable accuracy and is consistent with previously published subsurface conceptualizations of the region. The model indicates the overall geometry of the basement within the watershed and the spatial extent, altitude, and thickness of the basin-filling units. The hydrogeologic framework model is at a scale and resolution appropriate for use as the foundation for a numerical hydrologic model of the study area.

California↗
USGS2025

What 25+ years of "Did You Feel It" intensities tell us about shaking in California

“When will the Big One happen?” is a question that people often have for earthquake scientists. But while waiting for the “Big One” to occur, people will usually experience frightening or damaging shaking from multiple relatively smaller‐magnitude earthquakes. Given this context, it raises the question: “Where does most of the damage come from?” Could smaller, yet more frequent, earthquakes account for the majority of reported impactful shaking? To explore this question, we consider reports of earthquake damage and felt shaking experiences from a catalog of community‐collected intensity values from the U.S. Geological Survey’s “Did You Feel It?” system. Comparing these intensities to expectations from a ground‐motion model, we find that earthquakes of magnitudes smaller than expected are responsible for most reported intensities of community decimal intensities (CDI) 4.5 and above (moderate and higher shaking intensity levels). (Here “expected value” is meant in its mathematical sense of the mean or equivalently median shaking intensity.) We also present a regional analysis of observed earthquake shaking for specific areas, detailing the maximum intensity experienced within a predetermined area. We identify several instances of M < 4.5 events that generated maximum intensities of CDI > 5 in regions in California surrounding the cities of Eureka, San Francisco, Los Angeles, and San Diego. Our results motivate the need to include smaller‐magnitude earthquakes in communications about earthquake hazard and risk reduction.

California↗
USGS2025

A time-space model of graphite mineral systems

Increasing demand for graphite in energy storage systems warrants review of graphite ore genesis in a mineral systems framework. Orogenic graphite encompasses the metamorphic and orogenic mineral systems that produce flake graphite and hydrothermal vein (lump and chip) graphite deposits, respectively. A common feature of orogenic graphite deposits is an association with upper amphibolite- to granulite-facies metasedimentary rocks in continent-continent or continent-island arc collisional orogens. Orogenic flake graphite deposits form primarily through graphitization of organic carbon during regional metamorphism, but strain localization and partial melting of pelitic protoliths are likely important processes for graphite grade and quality enrichment. Orogenic vein graphite deposits precipitate from hydrothermal fluids, possibly derived from metamorphism or anatexis at depth. Decarbonation reactions in mixed calcareous-carbonaceous metasedimentary sequences are the most likely carbon sources for the veins. In contrast, intrusion-related graphite includes magmatic-hydrothermal and metamorphic mineral systems that form primarily in continental arc settings via the interaction of magmas with carbonaceous sedimentary packages. Magmatic-hydrothermal flake graphite deposits are hosted in plutonic and volcanic rocks, and result from the exsolution of CO 2 -CH 4 -rich fluids from melts contaminated by such packages. Contact metamorphism of carbonaceous sedimentary rocks by plutons produces some microcrystalline (amorphous) graphite deposits, including many in China. Compilation of geologic data from known graphite deposits globally suggests that pulses of carbon deposition in the Paleoproterozoic, Mesoproterozoic, and Neoproterozoic provided source carbon material. Subsequent supercontinent orogenesis at ca. 2,100 to 1,700 Ma (Columbia), ca. 1,300 to 1,000 Ma (Rodinia), and ca. 650 to 500 Ma (Gondwana) resulted in the genesis of orogenic flake and vein graphite deposits, where favorable geologic components overlapped with organic ± carbonate carbon-rich strata. Cryogenian deposition of graphite protoliths and Cryogenian – Cambrian metamorphic mineralization account for nearly 75% of all known resources globally and coincide with profound carbon isotope excursions and climate variability, implying a link with the global carbon budget. Comparatively few graphite deposits are associated with Pangea-forming orogens, attributed to less exhumation and/or denudation. High-temperature metasedimentary belts containing organic carbon-rich protoliths are most favorable for hosting orogenic flake graphite deposits, whereas sequences that also contain carbonate protoliths are favorable for orogenic graphite veins. Continent-scale orogenic belts may host both deposit types along with vanadium deposits. Use of a time-space mineral systems framework for graphite deposits can improve exploration models needed to ensure future supply of this critical mineral and provide insights into Earth’s long-term carbon cycle.

Mineralium Deposita↗
USGS2024

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

Critical minerals in orogenic (gold) and Coeur d’Alene-type mineral systems of the United States

Orogenic and Coeur d’Alene-type mineral systems are produced by metamorphic devolatilization of thick volcanic or siliciclastic sedimentary rock sequences and the focused flow of hydrothermal fluids upwards along crustal-scale faults. Most orogenic systems are found along the Cordilleran orogen, stretching from California northwards into Alaska, whereas most Coeur d’Alene-type systems occur in the Proterozoic Belt Basin in Idaho and Montana. Although the deposit types in these systems are exploited for precious and base metals, potential exists for the production of critical minerals necessary for current (2025) societal needs in the United States. Publicly available geochemical data compiled for these mineral systems, coupled with mineralogical characteristics, indicate that several critical minerals could potentially be recovered from unmined resources and processed mine waste: arsenic, antimony, tellurium, cobalt, and tungsten from orogenic gold deposits and zinc, antimony, arsenic, and manganese from Coeur d’Alene-type systems. These critical minerals reside primarily in arsenopyrite (arsenic), scheelite (tungsten), siderite (manganese), sphalerite (zinc), tetrahedrite (antimony and arsenic), stibnite (antimony), and telluride (tellurium) minerals.

continental United States↗
USGS2025

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2022

Potentiometric-surface and depth-to-water maps for spring 2022 were created for the Mississippi River Valley alluvial aquifer (MRVA) using groundwater-altitude data from 1,136 wells completed in the MRVA and from the altitude of the top of the water surface in area rivers from 160 streamgages. The potentiometric-surface and depth-to-water maps for 2022 were created to support investigations to characterize the MRVA as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the MRVA for spring 2022 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 340 feet (ft) above the North American Vertical Datum of 1988. The regional direction of groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient was into the depression, and near rivers, where the groundwater gradient can be from aquifer to the river or from the river into the aquifer. There are large depressions in the potentiometric-surface map in the lower one-half of the Cache region and in much of the Grand Prairie and Delta regions. Depth to water in the MRVA, spring 2022, by well ranged from 5.00 ft above land surface to 145.66 ft below land surface.

Mississippi River Valley alluvial aquifer↗
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗
USGS2025

Hydrogeologic investigation, framework, and conceptual flow model of the Antlers aquifer, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. Every 20 years, the Oklahoma Water Resources Board is required to update the hydrologic investigation on which the maximum annual yield determinations were based. The maximum annual yield allocated per acre of land is used to set the equal-proportionate share pumping rate. The maximum annual yield of 5,913,600 acre-feet per year and equal-proportionate-share of 2.1 acre-feet per acre per year currently (2025) in place for the Antlers aquifer were issued by the Oklahoma Water Resources Board on February 14, 1995. Because more than 20 years have elapsed since the 1995 final order for the Antlers aquifer was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, completed an in-depth hydrologic study that included a hydrogeologic framework and conceptual groundwater-flow model for the 1980–2022 study period. The results of an analysis of land use, long-term climate patterns, streamflow and base-flow patterns, historical groundwater use, as well as groundwater-level fluctuations across the Antlers aquifer are described. In addition, groundwater quality was analyzed for total dissolved solids concentrations and major ions for the Antlers aquifer. An updated hydrogeologic framework was developed that included refining the aquifer boundary in Oklahoma, the creation of new potentiometric surface and saturated thickness of fresh groundwater maps, one multiple-well aquifer test, slug tests, and an analysis of lithologic logs across the aquifer. A conceptual groundwater flow model and water budget were developed by incorporating estimates of recharge from precipitation, saturated-zone evapotranspiration, streambed seepage, lateral groundwater flows, vertical leakage, and withdrawals from groundwater wells.

Oklahoma, Texas↗
USGS2026

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years. Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr. Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued. Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Nevada↗
USGS2025

A spatial analysis of the groundwater emergence flood hazard in Long Island, New York and near coastal areas surrounding Long Island Sound in New York, Connecticut, and Rhode Island

Long Island, New York and near coastal areas surrounding Long Island Sound are densely populated and, like other coastal areas, are susceptible to flooding from several potential sources, including stormwater from precipitation events, tidal flooding and storm surge, and groundwater inundation or groundwater emergence flooding. The latter refers to the intersection of a rising water table with land surface or critical infrastructure. Many studies of flood drivers either neglect or only briefly discuss how shallow groundwater conditions may contribute to or exacerbate flood conditions. As part of a comprehensive study of compound flood hazards in the near coastal areas surrounding Long Island and Long Island Sound, a spatial analysis was completed, in cooperation with the Environmental Protection Agency’s Long Island Sound Study, using available regional datasets to characterize the potential hazard for groundwater emergence flooding. The approximately 3,100 square mile study area was subdivided into 11,407 900-meter by 900-meter (approximately 3,000-feet by 3,000-feet) grid cells, for the purposes of integrating the spatial datasets to calculate and map the groundwater emergence flood hazard. The depth to the water table, hydrologic soil groups, and National Land Cover Database were harmonized to the common grid. A groundwater emergence flood hazard rank was calculated for each grid cell for current average conditions following a set of rules accounting for the depth to the water table and the percent of area within each cell with slow infiltrating soils. A higher sea level position scenario was also calculated for the Long Island part of the study area. The calculated groundwater emergence flood hazard rank was reviewed in concert with the National Land Cover Data Base to identify developed areas and associated infrastructure that may be at risk to groundwater emergence flooding. Study results indicate that the groundwater emergence flood hazard is highest in coastal areas and near surface water where the water table is close to ground surface. Inland areas away from surface water bodies are not likely to be exposed to groundwater emergence flooding. For Long Island, under a scenario with higher sea level position, a greater groundwater emergence flood hazard is calculated in some locations closer to the coast and where land is submerged. Away from the coast and surface-water drainage, the groundwater emergence flood hazard is similar between the current average sea level condition and a higher sea level position scenario.

EarthArXiv↗
USGS2024

Parsimonious high-resolution landslide susceptibility modeling at continental scales

Landslide susceptibility maps are fundamental tools for risk reduction, but the coarse resolution of current continental-scale models is insufficient for local application. Complex relations between topographic and environmental attributes characterizing landslide susceptibility at local scales are not transferrable across areas without landslide data. Existing maps with multiple susceptibility classifications under-represent landslide potential in moderate and gently sloping terrain. We leverage an extensive landslide database ( N = 613,724), a high-resolution digital elevation model (10-m), and high-performance computing resources, to develop a new nationwide susceptibility map for the contiguous United States, Hawaii, Alaska, and Puerto Rico. We calculate four alternative linear and nonlinear thresholds of topographic slope and relief using an objective split-sample calibration. We down-sample our results to a 90-m grid to account for uncertainty in the digital elevation model and landslide position, and evaluate these thresholds' ability to differentiate areas of greater susceptibility. The less conservative nonlinear model optimally balances our priorities of capturing observed landslides (99%) while minimizing area covered by susceptible terrain (43%). Independent evaluation with four statewide landslide inventories ( N = 172,367) reinforces our model selection but highlights spatially variable performance. Therefore, we propose a novel approach to susceptibility classification using the concentration of landslide-prone terrain within each down-sampled grid. While landslides are possible within any cells containing susceptible terrain, those with the highest concentration capture the majority of observed landslides. Our new map characterizes landside susceptibility more consistently than prior models; our transparent classification approach also provides flexibility for accommodating different tolerances in risk reduction measures.

AGU Advances↗