USGS Science⌕ Search

SEARCH · USGS Science

Results for “Natural Resources Research”

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.

837 records · Page 23Linked to original sources
USGS2023

Sediment thickness map of United States Atlantic and Gulf Coastal Plain Strata, and their influence on earthquake ground motions

With the recent successful accounting of basin depth ground-motion adjustments in seismic hazard analyses for select areas of the western United States, we move toward implementing similar adjustments in the Atlantic and Gulf Coastal Plains by constructing a sediment thickness model and evaluating multiple relevant site amplification models for central and eastern United States seismic hazard analyses. We digitize and combine existing sediment thickness data sets into a composite surface that delineates the base of Cretaceous sediments under the Atlantic Coastal Plain and the base of Mesozoic sediments under the Gulf Coastal Plain. Amplification models dependent on sediment thickness, site natural period, and source-to-site path length are compared with data sets of observed ground motions to evaluate the ability of the new models to improve ground motion estimates. We find that the amplification models can account for observed trends in sediment-thickness and period-dependent residuals, but some tuning is required. For example, the model of Chapman and Guo requires a reference V S 30 , the time-averaged shear-wave velocity within 30 m of the Earth’s surface, for non-Coastal Plain sites, which we estimate to be between about 1 and 2 km/s. Along with our sediment thickness model, we estimate a velocity profile for application to the Harmon et al. site-natural-period-based model in order to best match the Chapman and Guo period dependence for a broad range of sediment thicknesses. The Next Generation of Attenuation models for the eastern United States Gulf Coast path-based adjustment models can also account for seismic attenuation in the Coastal Plain sediments and reduce the standard deviation of total residuals. If enacted in the U.S. Geological Survey National Seismic Hazard Model, these amplification models will reduce predicted short-period (<1 s) and increase predicted long-period (>1 s) ground motions in the Coastal Plains appreciably.

Earthquake Spectra↗
USGS2025

Water supply in the conterminous United States, Alaska, Hawaii, and Puerto Rico, water years 2010–20

We present an assessment of water supply across the conterminous United States (CONUS), Alaska, Hawaii, and Puerto Rico covering water years 2010–20. Our analysis drew on two national hydrologic models, the National Hydrologic Model Precipitation-Runoff Modeling System and the Weather Research and Forecasting model hydrologic modeling system. Both models produced estimates of streamflow, evapotranspiration, soil moisture, snow water equivalent, and other hydrologic states and fluxes. The models were driven by the bias-adjusted 4-kilometer-resolution, long-term regional hydroclimate simulation over the conterminous United States dataset (CONUS404). We assessed spatial and temporal error distributions by comparing monthly simulations at the 12-digit hydrologic unit code and regional scale from both models against external benchmarking datasets. Results showed that average annual rainfall across the CONUS was 857 millimeters per year for the period of analysis, with water year 2012 the driest year (729 millimeters) and water year 2019 the wettest year (995 millimeters). Key interannual variability results included the following: (1) the California–Nevada hydrologic region had the highest variability in precipitation and snow accumulation, and (2) the Texas hydrologic region was among hydrologic regions with the highest variability in precipitation. We related interannual variability in precipitation to storage volumes in soil moisture, snow water equivalent, and lakes and reservoirs to highlight areas with little storage and large year-to-year variability in precipitation. These areas included the Southern High Plains, Central High Plains, Texas, Souris–Red–Rainy, Mississippi Embayment, and Midwest regions. Our analysis of groundwater-level data showed that several of these areas overlap aquifers where groundwater levels were considerably lower than historical averages, including the Colorado Plateaus aquifers, the Rio Grande aquifer system, and the Central and Southern regions of the High Plains aquifer. Many of these lowered groundwater levels are continuations of decades-long declines from overpumping that started well before the assessment period. The resulting water budgets and their analyses provide a high-resolution foundational assessment of the mean state and variability of the terrestrial hydrologic cycle across the CONUS and Alaska, Hawaii, and Puerto Rico to support a wide range of water resource management applications.

Professional Paper↗
USGS2026

Continental-scale prediction of hydrologic signatures and processes

Understanding how dominant hydrologic processes and their drivers vary across diverse continental-scale landscapes is critical for hydrologic modeling and water management applications. Our research addresses this question by synthesizing large-sample watershed datasets, Caravan and GAGES-II, and developing random forest models to identify patterns in hydrologic function. We assessed dominant processes by examining hydrologic signatures – summary indicators of watershed function derived from hydroclimatic time series and random forest models across 14 146 gauged United States watersheds. The results reveal clear continental-scale gradients in hydrologic processes, including baseflow, overland flow, storage, and water balance losses. Our map of dominant processes highlights, for example, the transition from baseflow to fast responses and back to baseflow along the elevation gradient from the Appalachian spine, through the Piedmont, to the Eastern Coastal Plain; a distinct outer ring around the Great Lakes region; and sharp contrasts between coastal and inland processes in the West. Variable importance analysis from random forest models show that processes in the western U.S. are primarily controlled by climate, whereas in the eastern U.S., soil, geology, and topography play larger roles, with distinct human influences apparent in urban areas. Our approach of estimating dominant processes and their drivers facilitates extending process knowledge from research watersheds to the continental scale, assessing current hydrological understanding, and evaluating hydrological model structures.

conterminous United States↗
USGS2026

Bottom-up characterization of geologic methane emissions in the San Juan Basin in the southwestern USA

Methane is a potent greenhouse gas that plays an important role in atmospheric chemistry and global warming. The current global methane budget has large uncertainties, and a better understanding of the budget would help to guide strategies for reducing anthropogenic emissions to fight climate change. Natural geologic methane emissions are a particularly poorly constrained source, with top-down estimates from 14 C in ice cores suggesting much lower geologic emissions than bottom-up scaling of direct flux measurements. Our study aims to contribute to resolving this discrepancy through improved bottom-up characterization of geologic methane seepage in the San Juan Basin in southwestern Colorado and northwestern New Mexico, USA. We performed 983 new flux chamber measurements in this basin during summer 2022 and winter 2023 field campaigns. Our results, in combination with prior measurements, suggest that natural seepage in the San Juan Basin only occurs on or near the Fruitland coal outcrop. Specifically, our new measurements confirm previous measurements of seepage along the northwestern exposure of the Fruitland outcrop in Colorado (a known hydrodynamic overpressure region) and for the first time, identified seepage locations along the southernmost Fruitland outcrop exposure in New Mexico, in association with a coal cleat and a fault. Overall, seepage along the Fruitland coal outcrop is heterogeneously distributed, with both positive and negative (interpreted as microbial soil sink) methane fluxes. Features that are hypothesized to be predictive of seepage (e.g., faults) were not associated with positive methane fluxes in areas outside of the Fruitland outcrop. Our best estimate for total geologic methane seepage in the San Juan Basin from spatial interpolation and statistical upscaling is approximately 0.14 Tg CH 4 /yr, with a range from 0.029 to 0.48 Tg CH 4 /yr. This best-estimate value is lower than a previous bottom-up estimate from a gridded seepage inventory, but higher than a previous top-down estimate.

Colorado, New Mexico↗
USGS2025

Hydrogeologic characterization of the Cahuilla Valley and Terwilliger Valley Groundwater Basins, Riverside County, California

The relation between the groundwater and the amount of natural recharge to the Cahuilla Valley and Terwilliger Valley groundwater basins is not well understood. During the 20th century, the reliance on groundwater near Anza, California, used for agricultural, domestic, and municipal reasons has increased, and there is the potential for changes in groundwater availability related to climate change. Several types of existing data were evaluated, and new data were collected for this study, with the goal of characterizing the region’s hydrogeology. The study’s scope included constructing a geologic framework model to show where the groundwater-bearing units are present and their relation to each other, estimating the major components of the groundwater budget, and understanding local short-term and regional long-term groundwater flow and how that has changed since the early 1900s. Two electrical resistivity tomography surveys were done in the Durasno Valley about 2,150 feet apart to identify the thickness of the alluvium, its horizontal extent, and the depth-to-basement along two profiles perpendicular to Cahuilla Creek. The subsurface sediments were mostly horizontally layered and the transitional boundary between the alluvium and basement was thinner and shallower along the upgradient profile where the depth-to-basement was about 70 feet below land surface; the depth-to-basement at the downgradient profile was more than about 140 feet below land surface. The results from the surveys were used to place four monitoring wells at two sites along the survey profiles. Artesian flow from the deepest well at the downgradient site indicated that the decomposed and competent basement likely contributed some groundwater to the overlying alluvium, laterally, from below, or both. A digital three-dimensional geologic framework model was constructed using EarthVision software to represent the subsurface geometry of the alluvium, decomposed basement, and competent basement. Maps and cross sections of the modeled thicknesses of the alluvium and decomposed basement, and the modeled elevation of the top of the competent basement, were made to show the subsurface geometry of vertical faults, selected wells, and the groundwater-bearing units. Because natural recharge is related to the variable cycles of precipitation, estimates are difficult to quantify. Recharge and runoff have extreme interannual variability in the study area; recharge and runoff can be sporadic, and a substantive amount may not occur in some years. Estimates of recharge from a previous study and the regional-scale Basin Characterization Model for California for four different periods ranged from 3,800 acre-feet/year for 1897–1947 to 5,900 acre-feet/year for 1971–2000. Potential recharge from the disposal of domestic septic systems may have been as much as 500 acre-feet in 2020. It was estimated that between about 400 and 2,400 acre-feet/year of groundwater is lost through evapotranspiration by vegetation and evaporation from open water bodies, but the main source of discharge is through pumpage, mainly used for agriculture from the alluvium in the Cahuilla Valley and Terwilliger Valley groundwater basins. The estimated total pumpage for 1991–2021 ranged from about 1,140 acre-feet in 2019 to about 3,450 acre-feet in 1994. When summed, the cumulative amount of estimated pumpage between 1991 and 2021 was about 81,400 acre-feet. The general direction of groundwater flow is from the northeast along the San Jacinto fault zone at the headwaters of Cahuilla and Hamilton Creeks, to the surface-water outlets at the west and southeast parts of the study area. Groundwater-level data from the 1950s and earlier indicate that there was a natural groundwater divide between the Cahuilla Valley and Terwilliger Valley groundwater basins, but the changing magnitude and extent of the groundwater depressions caused by pumping since about 1950 indicate that the location of the natural groundwater boundary between the Cahuilla Valley and Terwilliger Valley groundwater basins has migrated over time. Flow from the upper to the lower parts of the Cahuilla Valley groundwater basin roughly follows the course of Cahuilla Creek through the narrow Durasno Valley where an estimated volume of flow in April 2019 was about 10–150 acre-feet/year. Short-term trends in groundwater levels, particularly in wells where groundwater is shallow and in the basement unit, show how some areas respond quickly to recharge and discharge. Wells located further to the east within the Cahuilla Valley groundwater basin in the alluvium show much less of a response to recharge events; areas of sustained pumpage from the alluvium, primarily for agriculture, show long-term declines in groundwater levels and generally do not show the effects of storm events or recent runoff. Groundwater levels in wells that are farthest from where most of the recharge occurs and where pumping has been the greatest, had some of the largest long-term groundwater-level declines at a rate of about 0.8 foot/year between 1971 and 2021.

California↗
USGS2025

Variable partitioning of lithium in rhyolitic melt during decompression and ascent

The partitioning behavior of Li in magmatic systems is increasingly being investigated due to the economic importance of Li in the transition to sustainable energy resources (e.g., batteries). However, at upper crustal pressures, it remains uncertain whether Li preferentially partitions into the vapor or liquid (brine) phase or remains in the silicate melt. This complicates our ability to determine where Li resides—silicate melt, minerals, or fluid phase—upon eruption, a crucial factor for understanding its postdepositional movement and concentration into a brine or volcano-sedimentary deposit. Here, we present a novel investigation into the behavior of Li within natural evolved melts during continuous magma decompression and ascent using melt embayments (open melt inclusions). Mineral-hosted melt embayments preserve records of the evolving composition of the exterior melt, including degassing pathways and ascent timescales, when paired with appropriate diffusion coefficients. Lithium concentration profiles were measured in quartz-hosted melt embayments from the rapidly quenched eruptive phases of five rhyolitic, caldera-forming eruptions to investigate the behavior of Li during magma decompression and ascent, where vapor partitioning and ascent dynamics were previously established by investigating H 2 O and CO 2 profiles. We find that in four systems, embayments contain lower interior Li concentrations than the coerupted melt inclusions; the fifth system contains the same Li concentrations in embayments and melt inclusions. However, many of these embayments contain gradients, with 84 % preserving Li enrichment near the melt-bubble interface, as compared to their interior concentration. We interpret these characteristics to represent two distinct stages of Li partitioning during magma decompression and ascent, in contrast to existing literature that proposes only one type of partitioning behavior. The first stage is interpreted as melt depletion of Li, likely driven by partitioning into an exsolved supercritical fluid phase, supported by the strong correlation between the extent of Li depletion and Cl concentration in the melt, as well as the decompression rate. This behavior then fundamentally shifts, where Li reenriches in the melt, postulated to be driven by the unmixing of the supercritical fluid phase at shallow pressures. For the one system that did not develop Li gradients through decompression, we attribute this to the lower values of Na and Cl in the melt, potentially inhibiting the partitioning of Li into a fluid phase. Importantly, the behavior of Li during decompression is not consistent within or between volcanic centers, highlighting the need for systematic experimental investigation in variable composition melts at pressures relevant to conduit dynamics. This knowledge would improve our ability to model Li profiles to understand magma decompression, and predict where Li resides (e.g., stored in volcanic glass, gas, or crystals) upon eruption prior to any later extraction.

Economic Geology↗
USGS2000

Water quality of selected springs and public-supply wells, Pine Ridge Indian Reservation, South Dakota, 1992-97

This report presents results of a water-quality study for the Pine Ridge Indian Reservation, South Dakota. The study was a cooperative effort between the U.S. Geological Survey and the Water Resources Department of the Oglala Sioux Tribe. Discharge and water-quality data were collected during 1992-97 for 14 contact springs located in the northwestern part of the Reservation. Data were collected to evaluate potential alternative sources of water supply for the village of Red Shirt, which currently obtains water of marginal quality from a well completed in the Inyan Kara aquifer. During 1995-97, water-quality data also were collected for 44 public-supply wells that serve about one-half of the Reservation's population. Quality-assurance sampling was used to evaluate the precision and accuracy of environmental samples. Ten of the springs sampled contact the White River Group, and four contact the Pierre Shale. Springs contacting the White River Group range from calcium bicarbonate to sodium bicarbonate water types. Two springs contacting the Pierre Shale have water types similar to this; however, sulfate is the dominant anion for the other two springs. In general, springs contacting the White River Group are shown to have better potential as alternative sources of water supply for the village of Red Shirt than springs contacting the Pierre Shale. Nine of the springs with better water quality were sampled repeatedly; however, only minor variability in water quality was identified. Six of these nine springs, of which five contact the White River Group, probably have the best potential for use as water supplies. Discharge from any of these six springs probably would provide adequate water supply for Red Shirt during most periods, based on a limited number of discharge measurements collected. Concentrations of lead exceeded the U.S. Environmental Protection Agency (USEPA) action level of 15 µg/L for three of these six springs. Five of these six springs also had arsenic concentrations that exceeded 10 µg/L, which could be problematic if the current maximum contaminant level (MCL) is lowered. Blending of water from one or more springs with water from the existing Inyan Kara well may be an option to address concerns regarding both quantity and quality of existing and potential sources. All nine springs that were sampled for indicator bacteria had positive detections on one or more occasions during presumptive tests. Although USEPA standards for bacteria apply only to public-water supplies, local residents using spring water for domestic purposes need to be aware of the potential health risks associated with consuming untreated water. One spring contacting the White River Group and two springs contacting the Pierre Shale exceeded 15 pCi/L for gross alpha; these values do not necessarily constitute exceedances of the MCL, which excludes radioactivity contributed by uranium and radon. Additional sampling using different analysis techniques would be needed to conclusively determine if any samples exceeded this MCL. Nine springs were sampled for selected pesticides and tritium. The pesticides atrazine, carbaryl, and 2,4-D were not detected in any of the samples. The nine springs were analyzed for tritium in order to generally assess the age of the water and to determine if concentrations exceeded the MCL established for gross beta-particle activity. Tritium results indicated two springs are composed primarily of water recharged prior to atmospheric testing of nuclear bombs and two other springs have a relatively large percentage of test-era water. The remaining five springs had tritium values that indicated some percentage of test-era water; however, additional sampling would be needed to determine whether water is predominantly pre- or post-bomb age. Of the 44 public-supply wells sampled, 42 are completed in the Arikaree aquifer, one is completed in an alluvial aquifer, and one is completed in the Inyan Kara aquifer. Water from the alluvial well is a sodium bicarbonate water type, water from Arikaree aquifer ranges from calcium bicarbonate to sodium bicarbonate types, and water from the Inyan Kara well is a calcium sulfate bicarbonate type. Of the 44 wells sampled, 28 (64 percent) tested positive for indicator bacteria in presumptive tests. Because these were single samples that generally were collected upstream from chemical treatment feeders, positive detections do not necessarily constitute exceedances of drinking-water standards. A single sample from an Arikaree well exceeded the MCL for arsenic of 50 µg/L. Arsenic exceeded 10 µg/L for six additional Arikaree wells and for the alluvial well and the Inyan Kara well, which could be problematic if the current MCL is lowered. The alluvial well also exceeded the secondary maximum contaminant level (SMCL) for dissolved solids, which is non-enforceable, and the action level for lead. The Inyan Kara well exceeded the SMCL's for iron and for manganese and the MCL of 5 pCi/L for radium-226 and 228 combined. Several Arikaree wells exceeded SMCL's for either pH, sulfate, dissolved solids, iron, or manganese. One Arikaree well exceeded the MCL of 4.0 mg/L for fluoride and another exceeded the MCL of 10 mg/L for nitrite plus nitrate. Ten Arikaree wells equalled or exceeded 15 pCi/L for gross alpha; however, these values do not necessarily constitute exceedances of the MCL, which excludes radioactivity contributed by uranium and radon. Additional sampling using different analysis techniques would be needed to conclusively determine if any samples exceeded this MCL. Eight wells, all from the Arikaree aquifer, equalled or exceeded the proposed MCL of 20 µg/L for uranium and 33 wells (75 percent) equalled or exceeded one-half of the proposed MCL. Although this standard has only been proposed, additional information regarding the extent of elevated uranium concentrations in the Arikaree aquifer, and the geochemical processes involved, may be beneficial. It was determined from analyses of uranium isotope data for five wells that the source of elevated uranium concentrations is naturally occurring, rather than anthropogenic.

South Dakota↗
USGS2026

Decadal trends in the quality of groundwater used for public drinking-water supply in California, 2004–2023, California groundwater ambient monitoring and assessment program, priority basin project

This study provides a comprehensive assessment of decadal changes in the quality of groundwater used for public drinking-water supply at 444 monitoring sites across California during 2004–2023. We assessed decadal step trends in groundwater quality for 145 water-quality constituents and geochemical indicators statewide and across geographic and land-use based network groups. We evaluated the statistical significance of directional changes (predominant increase or decrease of constituent concentrations) and the magnitude of those changes across all network groups. Uranium showed the most widespread directional and high-magnitude increases of all constituents with regulatory benchmarks statewide, particularly in the agriculture-dominated Central Valley as well as urban- and desert-dominated regions of Southern California. Fluoride and perchlorate showed the most widespread directional and high-magnitude decreases of all constituents with regulatory benchmarks statewide, which were also most pronounced in Southern California. Although arsenic and nitrate did not often register significant directional changes across network groups, they showed widespread, high-magnitude changes in both directions (increase and decrease) at levels often exceeding 10 percent of respective regulatory benchmarks statewide. Triazine herbicides (atrazine and simazine) and the gasoline oxygenate methyl tert-butyl ether (MTBE) showed significant directional decreases statewide, but not at levels considered to be of high magnitude compared to respective regulatory benchmarks. We observed significant directional and high-magnitude increases of total dissolved solids (TDS) statewide, which were most pronounced in agricultural areas. Analysis of explanatory geochemical indicators indicated that prevalent statewide increases of alkalinity and calcium were the predominant components of the observed statewide increases in TDS by mass. Widespread increases in groundwater alkalinity and calcium across agricultural and urban areas may be related, in part, to warm-season irrigation and other anthropogenic factors that have shifted soil weathering dynamics over the long term. Increasing alkalinity concentrations were related to increasing uranium concentrations, particularly in areas with aquifer materials derived from granitic rocks. Conversely, increasing calcium concentrations were related to decreasing fluoride concentrations, particularly in areas where fluoride occurred naturally at elevated concentrations. Decrease of perchlorate, triazine herbicides, and MTBE are likely related to decreased anthropogenic source inputs over time and natural attenuation in aquifers.

EarthArXiv↗
USGS2026

Deaf, deafblind, and hard of hearing university student experiences with earthquake early warning in the United States: Evaluating language planning and technology access

The growing literature on deaf and hard of hearing (DHH+) populations and disasters demonstrates that emergency communication (including alerts) is not reaching global DHH + individuals with dangerous impacts for morbidity and mortality. This is the first research study in the U.S. to qualitatively explore the experiences of DHH + persons with earthquake early warning (EEW) through group-based dialogue sessions. The study investigates eight DHH + university students'past earthquake experiences, access to EEW alerts, and perceptions of ShakeAlert Ⓡ , an EEW system for detecting earthquakes and alerting residents of California, Oregon, and Washington. Findings highlight key gaps in disaster alert usability within four thematic areas: lack of messaging in participants' language(s), unclear alert messaging, deficient message delivery mechanisms for deafblind persons, and insufficient access to earthquake information and training that leads to dependence on informal information networks. Weaknesses identified in these four themes reduce DHH + trust in EEW systems and compromise the capacity of alert recipients to take swift protective action or to mentally prepare before shaking starts. The study also underscores structural factors such as insufficient linguistic representation in disaster language planning and technology design, which ignores the linguistic and sensory access needs of DHH + individuals. Building on disaster language planning frameworks, we recommend involving DHH + populations to co-develop EEW alerts. By centering DHH + perspectives, this research contributes to ongoing efforts to ensure that EEW systems reach everyone.

California↗
USGS2026

Ungulate migrations of the Western United States, volume 6

This report, volume 6 in the “Ungulate Migrations of the Western United States” report series, showcases the migrations of 23 ungulate herds in the Western United States. The report series is produced by the Corridor Mapping Team (CMT). Led by the U.S. Geological Survey, the CMT is a collaboration among 11 State agencies, as well as regional and Federal partners, and an expanding number of Tribal wildlife agencies. The CMT was initiated in response to the U.S. Department of the Interior Secretarial Order 3362, which was signed in 2018 and provided Federal support to expand existing research efforts to study ungulate populations and conserve their migrations throughout the Western United States. Including this volume, the report series has detailed the migrations of 237 unique ungulate herds throughout the Western United States and continues to serve as a valuable resource to guide local and regional management, policy, and on-the-ground work necessary to maintain intact and functional ungulate migrations. This report highlights several guiding principles of the CMT that facilitate collaboration among the diverse set of partners and contribute to the program’s continued successes. Notably, raw global positioning system data are not shared among participating agencies and the U.S. Geological Survey, delineating migration corridors and seasonal ranges relies on empirical data, the CMT provides flexible approaches to participating State and Tribal partners, and regular CMT meetings create a framework for open communication among agency partners that supports transboundary mapping of migrations. The 237 ungulate migrations that have been included in the report series are an expanding inventory, which can help maintain ungulate migrations in perpetuity.

Arizona, California, Colorado, Idaho, Montana, Nev↗
USGS2023

Hydrogeology, karst, and groundwater availability of Monroe County, West Virginia

Monroe County is in southeastern West Virginia, encompassing an area of 474 square miles. The area consists of karst and siliciclastic aquifers of Ordovician, Silurian, Devonian, and Mississippian age and is in parts of two physiographic provinces: the Valley and Ridge Province to the east of Peters Mountain, and the Appalachian Plateau Province to the west of Peters Mountain. This study was developed in response to inquiries from the Monroe County Commission requesting assessment of the water resources of the county to better understand the quantity of the county’s groundwater resources, for both current [2023] and future demand, and to provide information to support protection and management of the county’s valuable groundwater resources. Various products were developed for this study that provide knowledge with respect to water availability and contamination susceptibility of the karst aquifers within the county. U.S. Geological Survey (USGS) geologists conducted extensive geologic mapping in support of the project, producing (1) a countywide bedrock geologic map, (2) a countywide hydrogeologic map, and (3) a light detection and ranging (lidar)-derived countywide digital elevation model and associated sinkhole map. A significant part of this work was to map in detail the Greenbrier Group at the formation level, which prior to this study had only partially been completed. The report also includes (4) a description of the lithologic units identified as part of the geologic mapping process. U.S. Geological Survey hydrologists completed several additional products for the hydrology part of the effort, including development of (1) a countywide potentiometric surface (water-table) map, (2) a countywide base-flow stream assessment, (3) countywide water-budget estimates, (4) well log surveys for 15 wells to better understand subsurface controls on groundwater flow within the study area, (5) two groundwater tracer tests to better refine the groundwater divide from the northern and southern parts of the karst aquifer in Monroe County; and finally, based on all available data collected for the study including the potentiometric surface map, geologic map, current [2023] and legacy fluorometric groundwater tracer tests, and base-flow stream assessments, (6) groundwater-basin delineations were reassessed for principal groundwater basins within the Greenbrier aquifer. In Monroe County, four principal hydrogeologic settings produce large yields of water for residential, agricultural, and other uses. The most relied upon water-bearing zone with respect to current [2023] public water supply is from springs along Peters Mountain. These springs are derived from intervals of fractured sandstone and resultant alluvial deposits. Groundwater flows downslope through these permeable alluvial deposits and discharges at the contact with less permeable strata, such as the Reedsville Shale. The second most relied upon water-bearing zone in Monroe County is within the karstic Greenbrier Group aquifer, in which the basal Hillsdale Limestone overlies the less permeable Maccrady Shale. This geologic contact between the Hillsdale Limestone and Maccrady Shale is not only targeted as a source of water for agricultural supply but also is targeted as a source of water for residential supply. The third most relied upon water-bearing zone is composed of shallow perched aquifers within the Greenbrier Group. The discontinuous nature of these perched aquifers makes mapping their extent impossible, but they are related to permeable geologic strata, such as karstified limestones with solutionally enhanced permeability that overlies less permeable shale or chert bedrock. During geologic mapping of the county, several of these perched aquifers were documented in the Pickaway, Union, and Alderson Limestones. A fourth zone consists of springs from Ordovician carbonates at the base of Peters Mountain, which are influenced by sinking streams as well as upwelling along faults. In terms of water quantity, the most sustainable springs are those having deeper-sourced flows. Public supplies are a principal source of water used for residential and commercial supply in the region, accounting for 0.49 million gallons per day (Mgal/d) of fresh-water withdrawals (0.14 Mgal/d of groundwater and 0.35 Mgal/d of surface water) for residential and commercial use and serving 6,645 individuals (49.2 percent of the population). An estimated 6,861 people, (50.8 percent of the population) primarily rely on private wells or other unregulated sources, such as springs, and withdraw 0.55 Mgal/d of groundwater for their residential use. Public water supply in the region is primarily (71.4 percent) derived from springs and augmented by stream withdrawals (backup sources mainly during low-flow periods), with the remaining portion (28.6 percent) derived from groundwater withdrawals from wells. For rural residents, however, 100 percent of their withdrawals are derived from groundwater (wells or springs).

West Virginia↗
USGS2026

Decadal trends in the quality of groundwater used for public drinking-water supply in California, 2004–23, California Groundwater Ambient Monitoring and Assessment Program, Priority Basin Project

This study provides a comprehensive assessment of decadal changes in the quality of groundwater used for public drinking-water supply at 444 monitoring sites across California during 2004–23. We assessed decadal step trends in groundwater quality for 145 water-quality constituents and geochemical indicators statewide and across geographic and land-use based network groups. We evaluated the statistical significance of directional changes (predominant increase or decrease of constituent concentrations) and the magnitude of those changes across all network groups. Uranium showed the most widespread directional and high-magnitude increases of all constituents with regulatory benchmarks statewide, particularly in the agriculture-dominated Central Valley as well as urban-and desert-dominated regions of Southern California. Fluoride and perchlorate showed the most widespread directional and high-magnitude decreases of all constituents with regulatory benchmarks statewide, which were also most pronounced in Southern California. Although arsenic and nitrate did not often register significant directional changes across network groups, they showed widespread, high-magnitude changes in both directions (increase and decrease) at levels often exceeding 10 percent of respective regulatory benchmarks statewide. Triazine herbicides (atrazine and simazine) and the gasoline oxygenate methyl tert -butyl ether showed significant directional decreases statewide, but not at levels considered to be of high magnitude compared to respective regulatory benchmarks. We observed significant directional and high-magnitude increases of total dissolved solids statewide, which were most pronounced in agricultural areas. Analysis of explanatory geochemical indicators indicated that prevalent statewide increases of alkalinity and calcium were the predominant components of the observed statewide increases in total dissolved solids by mass. Widespread increases in groundwater alkalinity and calcium across agricultural and urban areas may be related, in part, to warm-season irrigation and other anthropogenic factors that have shifted soil weathering dynamics over the long term. Increasing alkalinity concentrations were related to increasing uranium concentrations, particularly in areas with aquifer materials derived from granitic rocks. Conversely, increasing calcium concentrations were related to decreasing fluoride concentrations, particularly in areas where fluoride occurred naturally at elevated concentrations. Decrease of perchlorate, triazine herbicides, and methyl tert -butyl ether are likely related to decreased anthropogenic source inputs over time and natural attenuation in aquifers.

California↗
USGS2025

Cenozoic stratigraphy of Colorado

As a successor to previous Colorado stratigraphy charts (MS-53 Colorado Stratigraphic Chart and MS-54 Cretaceous Stratigraphy of Colorado), this Colorado Geological Survey (CGS) publication resulted from a collaboration between the CGS, USGS, and the Denver Museum of Nature and Science (DMNS). The chart was designed to illustrate Cenozoic stratigraphy spanning the state’s many sedimentary basins. It builds upon the work of dozens of colleagues and updates Richard Pearl’s seminal 1974 stratigraphy chart. The chart leverages the community’s stratigraphic work in both the subsurface and outcrop, and depicts new geochronologic constraints for many units. To facilitate comparison of strata to external forcing factors, the chart employs a linear timescale. Each unit’s dominant depositional environment is depicted, as are major mountain building events, erosional events, and regional unconformities.

Colorado↗
USGS2024

ShakeAlert® and schools: Incorporating earthquake early warning in school districts in Alaska, California, Oregon, and Washington

The U.S. Geological Survey-managed ShakeAlert® earthquake early warning system is the first public alerting system in the United States to provide rapid mass notification when an earthquake is detected. Although public alert delivery via mobile phones began in California in 2019 followed by Oregon and Washington in 2021, little is known about what might drive widespread implementation in at-risk institutional settings such as schools. For example, there has been limited research on how to best integrate earthquake early warning into existing emergency plans, alert systems, and drills to keep school children and personnel safe in an earthquake. To address this gap, in the spring of 2022, every school district superintendent in Alaska, California, Oregon, and Washington was sent a 15-min online survey. The survey assessed superintendent knowledge of ShakeAlert, preferences for alert messaging, and perceived opportunities and barriers to incorporating the system in schools. The results showed that superintendents had low awareness of ShakeAlert but held positive perceptions of the system's potential to enable life-saving protective actions. A major barrier to adoption included the perceived financial cost of implementing and maintaining the system. There were some statistically significant differences in state responses, and future research could investigate the specific needs of each state based on school district size and composition, hazard exposure, and earthquake experience. Together these findings can help inform targeted strategies to increase ShakeAlert adoption in schools and ultimately improve the safety of school children and staff during earthquakes.

Alaska, California, Oregon, Washington↗
USGS2024

A benchmark dataset and workflow for landslide susceptibility zonation

Landslide susceptibility shows the spatial likelihood of landslide occurrence in a specific geographical area and is a relevant tool for mitigating the impact of landslides worldwide. As such, it is the subject of countless scientific studies. Many methods exist for generating a susceptibility map, mostly falling under the definition of statistical or machine learning. These models try to solve a classification problem: given a collection of spatial variables, and their combination associated with landslide presence or absence, a model should be trained, tested to reproduce the target outcome, and eventually applied to unseen data. Contrary to many fields of science that use machine learning for specific tasks, no reference data exist to assess the performance of a given method for landslide susceptibility. Here, we propose a benchmark dataset consisting of 7360 slope units encompassing an area of about 4,100 km 2 "> 4,100 km 2 in Central Italy. Using the dataset, we tried to answer two open questions in landslide research: (1) what effect does the human variability have in creating susceptibility models; (2) how can we develop a reproducible workflow for allowing meaningful model comparisons within the landslide susceptibility research community. With these questions in mind, we released a preliminary version of the dataset, along with a “call for collaboration,” aimed at collecting different calculations using the proposed data, and leaving the freedom of implementation to the respondents. Contributions were different in many respects, including classification methods, use of predictors, implementation of training/validation, and performance assessment. That feedback suggested refining the initial dataset, and constraining the implementation workflow. This resulted in a final benchmark dataset and landslide susceptibility maps obtained with many classification methods. Values of area under the receiver operating characteristic curve obtained with the final benchmark dataset were rather similar, as an effect of constraints on training, cross–validation, and use of data. Brier score results show larger variability, instead, ascribed to different model predictive abilities. Correlation plots show similarities between results of different methods applied by the same group, ascribed to a residual implementation dependence. We stress that the experiment did not intend to select the “best” method but only to establish a first benchmark dataset and workflow, that may be useful as a standard reference for calculations by other scholars. The experiment, to our knowledge, is the first of its kind for landslide susceptibility modeling. The data and workflow presented here comparatively assess the performance of independent methods for landslide susceptibility and we suggest the benchmark approach as a best practice for quantitative research in geosciences.

Earth-Science Reviews↗
USGS2024

The spatially adaptable filter for error reduction (SAFER) process: Remote sensing-based LANDFIRE disturbance mapping updates

LANDFIRE (LF) has been producing periodic spatially explicit vegetation change maps (i.e., LF disturbance products) across the entire United States since 1999 at a 30 m spatial resolution. These disturbance products include data products produced by various fire programs, field-mapped vegetation and fuel treatment activity (i.e., events) submissions from various agencies, and disturbances detected by the U.S. Geological Survey Earth Resources Observation and Science (EROS)-based Remote Sensing of Landscape Change (RSLC) process. The RSLC process applies a bi-temporal change detection algorithm to Landsat satellite-based seasonal composites to generate the interim disturbances that are subsequently reviewed by analysts to reduce omission and commission errors before ingestion them into LF’s disturbance products. The latency of the disturbance product is contingent on timely data availability and analyst review. This work describes the development and integration of the Spatially Adaptable Filter for Error Reduction (SAFER) process and other error and latency reduction improvements to the RSLC process. SAFER is a random forest-based supervised classifier and uses predictor variables that are derived from multiple years of pre- and post-disturbance Landsat band observations. Predictor variables include reflectance, indices, and spatial contextual information. Spatial contextual information that is unique to each contiguous disturbance region is parameterized as Z scores using differential observations of the disturbed regions with its undisturbed neighbors. The SAFER process was prototyped for inclusion in the RSLC process over five regions within the conterminous United States (CONUS) and regional model performance, evaluated using 2016 data. Results show that the inclusion of the SAFER process increased the accuracies of the interim disturbance detections and thus has potential to reduce the time needed for analyst review. LF does not track the time taken by each analyst for each tile, and hence, the relative effort saved was parameterized as the percentage of 30 m pixels that are correctly classified in the SAFER outputs to the total number of pixels that are incorrectly classified in the interim disturbance and are presented. The SAFER prototype outputs showed that the relative analysts’ effort saved could be over 95%. The regional model performance evaluation showed that SAFER’s performance depended on the nature of disturbances and availability of cloud-free images relative to the time of disturbances. The accuracy estimates for CONUS were inferred by comparing the 2017 SAFER outputs to the 2017 analyst-reviewed data. As expected, the SAFER outputs had higher accuracies compared to the interim disturbances, and CONUS-wide relative effort saved was over 92%. The regional variation in the accuracies and effort saved are discussed in relation to the vegetation and disturbance type in each region. SAFER is now operationally integrated into the RSLC process, and LANDFIRE is well poised for annual updates, contingent on the availability of data.

Fire↗
USGS2026

Assessment of extreme subsurface hydrologic conditions captured during atmospheric river storms in the San Francisco Bay area (California, USA) with applications to shallow landslide early warning

An increase in soil pore water pressure is the typical trigger for the majority of landslides caused by rainfall. Atmospheric river storms, common to the west coast of North America during the winter season, can deliver landslide triggering rainfall resulting in severe impacts to coastal communities. Using a network of hydrologic monitoring stations situated within landslide-prone terrain in the San Francisco Bay area of California (USA), we assess the meteorologic conditions and resulting hydrologic and landslide response resulting from eight consecutive storm events that caused thousands of shallow landslides during the winter of 2022–2023. We find disparate hydrological responses and resultant degrees of observed landsliding ranging from < 1 landslide/km2 to 18 landslides/km2 at the monitoring sites that reflect the interplay and differences between rainfall delivery, subsurface hydrological characteristics, and geotechnical properties at each site. Antecedent soil moisture from both early season rainfall and the first storm in the sequence played a critical role in setting up some hillslopes for failure. Subsequent storms then generated elevated pore water pressures for several hours with associated landsliding. However, we find that the occurrence of widespread landsliding required not only sufficient pore water pressure magnitude in susceptible hillslopes, but also full and prolonged effective soil saturation throughout hillslope profiles. Landslides may still occur at lower values and durations of effective saturation but are likely to be less extensive regionally. We present these findings within the context of research directions and improvements to landslide early warning systems first suggested by researchers 40 years ago.

California↗
USGS2025

The 3D National Topography Model Call for Action—Part 2: The Next Generation 3D Elevation Program

The three-dimensional (3D) National Topography Model initiative to integrate elevation and hydrography data includes the next generation of hydrography data from the 3D Hydrography Program and the next generation of elevation data from the 3D Elevation Program (3DEP). The first-ever collection of light detection and ranging (lidar) data for the nation (IfSAR for Alaska) provides a critical baseline reference, and the addition of multiple repeat elevation mapping projects as part of the next generation of 3DEP would substantially expand analysis capabilities. As the U.S. Geological Survey (USGS) is closing in on our goal of complete coverage with 3DEP data available or in progress for 98.3 percent of the Nation at the end of fiscal year 2024, the USGS is already transitioning to the next generation of 3DEP. Based on the 3D Nation Study results and input from a broad range of stakeholders, the USGS National Geospatial Program has finalized a new design for 3DEP that provides increased lidar quality levels and refresh rates. The new program is designed with more flexibility to meet changing user needs and take advantage of improvements in mapping technologies. The program will aim to expand the level of interagency coordination for topobathymetric lidar acquisition for inland rivers. The next generation of 3DEP will also aim to emphasize research, including advancing program design, products, and services and engaging and leveraging the evolving 3D industry. Research goals also include becoming more flexible in meeting user needs and taking advantage of evolving remote-sensing technologies. The program also plans to move from focusing on producing standard products to producing a concept of a 3D Nation Ecosystem with a variety of inputs, products, and services.

Circular↗