USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Waterway, Port, Coastal, and Ocean Engineering”

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,722 records · Page 44Linked to original sources
USGS2026

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the effects of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study, we constructed the Strawberry Creek integrated hydrological model to increase this understanding. The Strawberry Creek integrated hydrological model is an integrated surface runoff and groundwater model, GSFLOW to simulate surface runoff and infiltration and groundwater flow. The model has 263 rows and 176 columns in each of its three layers. The model area includes the Strawberry Creek and four adjacent watersheds. The precipitation-runoff modeling system was calibrated using two streamflow gaging stations, and the groundwater and surface-water flow model was calibrated to the reported spring diversion discharge and a sparse number of groundwater-level measurements. The Strawberry Creek integrated hydrological model was run with and without diversions active, and simulated streamflow was compared, indicating that in the headwaters of Strawberry Creek, about 35 percent of the diversion flow was captured from base flow.

California↗
USGS2026

Reconstructing the Quaternary depositional history using geologic mapping and three-dimensional modeling of the subsurface near Fort Morgan, northeastern Colorado

Centered on Fort Morgan, Colorado, this study is intended to build from previous work by adding a three-dimensional (3D) view of the subsurface to better understand the depositional history of Quaternary deposits. A 1:100,000 scale geologic map was made by combining previous geologic maps, regional soil maps, and recent field investigations. In addition to the geologic mapping, drill hole lithologic data from water wells and oil and gas exploration were compiled and lithologic units simplified to best represent the stratigraphy of the Quaternary deposits. From these subsurface data, a 3D subsurface model was constructed, trimmed at the surface by a digital elevation model, and a bedrock surface foundation gridded from drill hole data was added. The surface of the 3D model was then compared visually to the surficial geologic map. Cross sections were constructed from the 3D model and compared to site-specific drilling that was done as part of this project. Finally, the model was examined in detail to reconstruct the depositional history of the subsurface alluvial and eolian units. Alluvial and fluvial drainage basins exposed in the subsurface have a greater areal extent than the present-day narrow drainages. Older eolian sand in the subsurface tends to be interbedded with loess indicating coeval deposition. Holocene sand, both eroded from bedrock exposed at the surface north of the study area and reworked from the South Platte River, buries most of the interbedded older sand and loess.

Colorado↗
USGS2025

Microbial ecology of permafrost soils: Populations, processes, and perspectives

Permafrost microbial research has flourished in the past decades, due in part to improvements in sampling and molecular techniques, but also the increased focus on the permafrost greenhouse gas feedback to climate change and other ecological processes in high latitude and alpine permafrost soils. Permafrost microorganisms are adapted to these extreme environments and remain active at low temperatures and when resources are limited. They are also an important component of global elemental cycles as they regulate organic matter turnover and greenhouse gas production, particularly as permafrost thaws. Here we review the permafrost microbiology literature coupled with an exploration of its historical aspects, with a particular focus on a new understanding advanced by molecular biology techniques. We further identify knowledge gaps and ways forward to improve our understanding of microbial contributions to ecosystem biogeochemistry of permafrost-affected systems.

Permafrost and Periglacial Processes↗
USGS2025

Quantitative subsurface characterization illuminates the origin of the Quaternary Mississippi River Valley alluvial aquifer

The Mississippi River Valley alluvial aquifer (MRVA) is vital to U.S. food security and global agricultural supply. However, quantitative understanding of its Quaternary origin, architecture, and hydrologic function remains incomplete. Here we develop a three-dimensional hydrostratigraphic model to characterize the deposition of clay and silt, fine-medium sands, and graveliferous sands using lithologic data from 75,000 boreholes compiled across the Lower Mississippi Valley and a geostatistical method—interval kriging. We find that cyclic glacial entrenchments, evidenced by remnants of pre-Wisconsinan postglacial sediments, alongside geodynamic activities shaped the MRVA basal configuration. Stratal weakening from faulting and salt diapirism enhanced glacial incision and thereby produced abrupt aquifer thickening. We demarcate the top of graveliferous sands as the regional marker of the Pleistocene-Holocene transition. The MRVA hydrostratigraphy reveals hydrologic function and geologic controls on groundwater storage and quality, advancing the assessment of aquifer sustainability under a changing climate, with implications for alluvial aquifers globally.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗
USGS2025

Assessment of density pattern retention of generalized data for 1:100,000-scale United States topographic maps

Cartographic generalization reduces the complexity of geographic data to produce legible, smaller-scale displays that retain essential information and logical geographic patterns. Generalization is a vital process in topographic map production. An important challenge in this process is managing and evaluating consistency across scale in the density and spatial distribution of map features such as buildings, roads, streams, water bodies, and elevation contours. Density patterns in these features reflect underlying physiographic conditions, which include factors such as bedrock geology, tectonics, climate, and landforms. Assessments of an acceptable level of change in feature density patterns are critical to ensuring the readability, usability, and accuracy of generalized maps and data. Preserving realistic density patterns across mapping scales also supports sustainable development goals in cartography, by helping to prioritize and communicate the relative reliability of geospatial data at specific scales.

Conference Paper↗
USGS2025

Colored shaded-relief bathymetry and acoustic backscatter of Lake Sammamish, Washington

Evidence of strong earthquakes (such as underwater landslides and associated deposits) may be recorded within the lacustrine sediments of Pacific Northwest lakes. The floor of Lake Sammamish, Wash., an approximately 11 kilometer (6.8 mile) long, 2 kilometer (1.2 mile) wide, and 35 meter (114.8 feet) deep lake located in a populated region just east of Seattle, was mapped by the U.S. Geological Survey in November of 2021 to search for evidence of past earthquakes. Mapping was conducted using a SWATHplus-M 234-kHz interferometric side-scan sonar system was pole-mounted on the U.S. Geological Survey research vessel Parke Snavely, and the system collected full-coverage bathymetric and acoustic backscatter data, which were processed to 2-meter spatial resolution. Two maps were created, a colored shaded-relief bathymetric map showing lake floor morphology (sheet 1), and an acoustic-backscatter map showing backscatter intensities (sheet 2). The results may then be utilized together to investigate past earthquake activity.

Washington↗
USGS2025

Methodology and technical input for the 2025 U.S. List of Critical Minerals—Assessing the potential effects of mineral commodity supply chain disruptions on the U.S. economy

The Secretary of the Interior, acting through the Director of the U.S. Geological Survey, is tasked by section 7002 (“Mineral Security”) of title VII (“Critical Minerals”) of the Energy Act of 2020 (Public Law 116–260, December 27, 2020, 116th Congress) with reviewing and revising the methodology used to evaluate mineral commodity supply risk and the U.S. List of Critical Minerals (LCM) no less than every 3 years. Following two previous LCM assessments, this analysis represents the latest technical input for evaluating each mineral commodity’s supply risk and determining their recommended status on the LCM. We evaluated mineral commodity supply risk using two criteria: (1) an economic effects assessment that quantified the potential effects of various trade disruption scenarios on the U.S. economy, and (2) an examination of whether the mineral commodity’s U.S. supply chain relied on a sole domestic producer that represented a single point of failure. For the first criterion, postdisruption equilibrium quantities and prices for each mineral commodity were calculated based on their price elasticities of supply and demand and the availability of excess production capacity for each yearlong foreign trade disruption scenario. Subsequently, a nonlinear optimization routine was used with detailed economic input-output tables to estimate the potential economic effects on the U.S. economy of over 1,200 scenarios for 84 mineral commodities. After accounting for the probability of each scenario’s occurrence, the overall results are presented in terms of changes in U.S. gross domestic product (GDP) by individual industry and the economy overall. The results, which ranged from a net decrease in U.S. GDP of nearly $4.5 billion to a net increase of $33 million, largely reflect U.S. import dependency and world production concentration. Using the Jenks natural breaks optimization method, a statistical classification technique, we categorized the mineral commodities into several classes based on this overall risk quantification. Mineral commodities with annualized probability-weighted net decreases in U.S. GDP greater than $2 million were recommended for inclusion on the LCM. If a mineral commodity did not meet the threshold for inclusion on the LCM under the first criterion, its domestic supply chain was examined under the second criterion, which recommended a mineral commodity for inclusion on the LCM if there was only a single domestic producer. Ultimately, the two criteria resulted in the recommendation of the addition of six mineral commodities (in descending risk order, potash, silicon, copper, silver, rhenium, and lead) to and the removal of two mineral commodities (arsenic and tellurium) from the LCM. By using an economic effects assessment, the results of this analysis provide a prioritization that can also be compared directly against other risk analyses and the cost of various risk mitigation strategies.

Open-File Report↗
USGS2024

Agricultural tile drains increase the susceptibility of streams to longer and more intense streamflow droughts

Streamflow droughts are receiving increased attention worldwide due to their impact on the environment and economy. One region of concern is the Midwestern United States, whose agricultural productivity depends on subsurface pipes known as tile drains to improve trafficability and soil conditions for crop growth. Tile drains accomplish this by rapidly transporting surplus soil moisture and shallow groundwater from fields, resulting in reduced watershed storage. However, no work has previously examined the connection between tile drainage and streamflow drought. Here, we pose the question: does the extent of watershed-level tile drainage lead to an increased susceptibly and magnitude of streamflow droughts? To answer this, we use daily streamflow data for 122 watersheds throughout the Midwestern United States to quantify streamflow drought duration, frequency, and intensity. Using spatial multiple regression models, we find that agricultural tile drainage generates statistically significant ( p < 0.05) increases in streamflow drought duration and intensity while significantly reducing drought frequency. The magnitude of the effect of tile drainage on streamflow drought characteristics is similar to that of water table depth and precipitation seasonality, both of which are known to influence streamflow droughts. Furthermore, projected changes in regional precipitation characteristics will likely drive the installation of additional tile drainage. We find that for each 10% increase in tile-drained watershed area, streamflow drought duration and intensity increase by 0.03 d and 12%, respectively, while frequency decreases by 0.10 events/year. Such increases in tile drainage may lead to more severe streamflow droughts and have a detrimental effect on the socio-environmental usage of streams throughout the Midwest.

Illinois, Indiana, Iowa, Michigan, Minnesota, Ohio↗
USGS2025

Groundwater hydrology, groundwater and surface-water interactions, water quality, and groundwater-flow simulations for the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado, 2017–19

In 2017, the U.S. Geological Survey, in cooperation with the Upper Arkansas Water Conservancy District, began a study to provide a comprehensive analysis of the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado. The study included collection of data pertaining to groundwater hydrology, groundwater and surface-water interactions, and water quality in the alluvial aquifer. In addition to providing foundational information on the hydrology of the alluvial aquifer, a numerical groundwater-flow model was developed to estimate the potential effects of additional storage of groundwater in the alluvial aquifer. Groundwater-level elevation data from 30 wells were used to estimate groundwater-flow directions in the alluvial aquifer, which were generally from the southwest to northeast, away from the Sangre de Cristo Mountains and towards perennial streams in the center of the valley. Although some seasonal variation was apparent in groundwater-level elevation records, no statistically significant seasonal trends were indicated. Statistically significant long-term trends were indicated in groundwater-level elevation records for 8 of the 30 wells, and of these wells with statistically significant trends, all but 1 indicated a negative trend of groundwater-level elevations. Spatial evaluation of wells with statistically significant negative groundwater-level elevation trends showed many are in areas of denser well drilling for domestic or other uses, indicating increasing groundwater use could potentially be causing groundwater-level elevation declines. There were instances of wells with no statistically significant groundwater-level elevation trends also located in areas of greater density of well completions. Additional investigations may be necessary to more fully characterize the processes responsible for negative groundwater-level elevation trends. Streamflow gain or loss calculations were completed for low flow in 2017–19 and for high flow in 2018 in nine reaches of streams within the study area. Stream reaches of the upper Texas Creek, upper Grape Creek, upper-middle Grape Creek, and Taylor Creek displayed consistent streamflow loss in each period from 2017 to 2019. These stream reaches represent long-term sources of recharge to the alluvial aquifer. Streamflow gain or loss varies through time in other stream reaches (lower Texas Creek, lower-middle Grape Creek, lower Grape Creek below Westcliffe, and lower Grape Creek above DeWeese Reservoir). The temporally variable behavior indicates these stream reaches may be sources of groundwater recharge or areas of groundwater discharge, likely depending on temporal dynamics between the elevation of the water table and the stream. Water-quality samples were collected from 10 groundwater wells and 10 stream sites during September through November 2019. All groundwater and stream samples were analyzed for major and trace elements and stable isotopes of water. A subset of groundwater samples was also analyzed for the environmental tracers sulfur hexafluoride, tritium, and noble gases. Comparison of water-quality results to U.S. Environmental Protection Agency drinking water-quality standards indicated no constituents exceeded primary standards for human health. Spatial evaluation of water quality indicated the concentrations of various constituents are likely controlled by groundwater and surface-water interactions and by spatial variability in bedrock geology underlying the alluvial aquifer. Specifically, streams shown to gain from groundwater had water chemistry constituent compositions similar to groundwater, whereas streams exiting the Sangre de Cristo Mountains tended to have compositions consistent with snowmelt. Groundwater geochemistry appeared to be partially controlled by oxidation-reduction processes and by proximity to igneous rocks in the Wet Mountains. Environmental tracers used to estimate groundwater age indicated all sampled groundwater contained tracers representing modern recharge (approximately less than 65 years old) but mixing of premodern recharge (approximately more than 65 years old) also occurs. Spatial evaluation of environmental tracers indicated large faults may be conduits for upwelling of older groundwater. No trends were observed in groundwater age with well depth, indicating all sampled wells are located within the zone of active groundwater flow. The presence of modern groundwater in wells with statistically significant negative groundwater-level elevation trends indicates groundwater storage depletions may be partially offset by capture of modern recharge. Repeated sampling of groundwater age would be necessary, however, to determine if any trends in groundwater age exist, which may indicate changing groundwater recharge, storage, or discharge. Additional investigations could also consider quantifying groundwater age in deeper wells to more fully define the depth of active groundwater flow. A numerical groundwater-flow model was developed to estimate components of the water budget, simulate groundwater and surface-water interactions, and evaluate the potential effects of aquifer storage and recovery. Simulated groundwater-level elevations from the calibrated groundwater-flow model are similar to the observed pattern of groundwater-level elevations with higher elevations in the western part of the study area along the Sangre de Cristo Mountains. Simulated water-budget components indicate most of the recharge to the alluvial aquifer is derived from streamflow losses, which is consistent with observations of losing streams along the mountain front. The largest groundwater discharge component of the alluvial aquifer was to streams in the center of the valley, where observations of stream gain or loss indicated the predominance of gaining conditions. Comparison of groundwater and surface-water interactions between the calibrated groundwater-flow model for 2000-19 (the base-case model) and a simulation including additional recharge, representing potential aquifer storage and recovery operations, indicated the additional recharge distributed throughout the area had minimal effects on streamflow in the nearby Grape Creek. An analysis of subregional groundwater budgets showed approximately 54 percent of the additional recharge flowed back to nearby Grape Creek, and the other 46 percent was distributed laterally into adjacent cells in the alluvial aquifer. The comparison of simulations and subregional water budget show the additional recharge did not substantially alter groundwater-level elevations or basin wide groundwater storage. Although the analysis of additional recharge provided in the numerical groundwater-flow model considers only one of many possible recharge scenarios, the model provides a useful tool that could be modified for various scenarios to understand potential effects of managed aquifer recharge.

Colorado↗
USGS2024

Static and dynamic strain in the 1886 Charleston, South Carolina, earthquake

During the 1886 Mw 7.3 Charleston, South Carolina, earthquake, three railroads emanating from the city were exposed to severe shaking. Expansion joints in segmented railroad tracks are designed to allow railroad infrastructure to withstand a few parts in 10,000 of thermoelastic strain. We show that, in 1886, transient contractions exceeding this limiting value buckled rails, and transient extensions pulled rails apart. Calculated values for dynamic strain in the meizoseismal region are in reasonable agreement with those anticipated from the relation between strain and moment magnitude proposed by Barbour et al. (2021) and exceed estimated tectonic strain released by the earthquake by an order of magnitude. Almost all of the documented disturbances of railroad lines, including evidence for shortening of the rails, can thus be ascribed to the effects of dynamic strain changes, not static strain. Little or no damage to railroads was reported outside the estimated 10 −4 dynamic strain contour. The correspondence between 10 −3 and 2×10 −4 contours of dynamic strain and Mercalli intensity 9 and 8, anticipated from the dependence of each quantity on peak ground velocity, suggests it may be possible to use railroad damage to quantitatively estimate shaking intensity. At one location, near Rantowles, ≈20 km west of Charleston, a photograph of buckled track taken one day after the earthquake has been cited as evidence for shallow dextral slip and has long focused a search for a causal fault in this region. Photogrammetric analysis reveals that the buckle was caused by transient contraction of <10 cm with no dextral offset. Our results further weaken the evidence for faulting in the swamps and forests south of the Ashley River in 1886, hitherto motivated by the photograph and limited macroseismic evidence for high‐intensity shaking.

South Carolina↗
USGS2026

Global pegmatite-hosted lithium, cesium, and rubidium resources: A dataset for grade and tonnage modeling

Quantitative mineral resource assessments of potential undiscovered deposits can inform future mineral supply scenarios, but their accuracy is conditional on building robust grade and tonnage models of known deposits. This study presents an up-to-date global compilation and analysis of recently discovered and original, in-situ pegmatite-hosted Li, Cs, and Rb resources prior to historic production. Our analysis yields a median tonnage of 21.2 million tons (Mt) and grade of 1.12% Li 2 O, respectively, for global Li pegmatite deposits ( n = 73). The grades and tonnages of Li pegmatite resources vary depending on the age of the bedrock host domain, pegmatite crystallization age, and primary ore mineralogy. Lithium pegmatite resources hosted in Archean to transitional Archean-Paleoproterozoic domains have the largest median tonnage (29.8 Mt; n = 38), and those hosted in Paleoproterozoic to Mesoproterozoic domains have smaller median tonnages (6.5 Mt; n = 16). Cesium deposits where pollucite is the primary ore mineral have a bimodal grade distribution, with modes of 2.40 and 0.035 wt% Cs 2 O for high- and low-grade deposits, respectively, while Rb deposits are more unimodal with a median grade of 0.247 wt% Rb 2 O. Pegmatite-hosted Cs and Rb resources have median tonnages of 7.6 and 6.3 Mt, respectively. Covariation between ore mineralogy and the degree of crustal enrichment in pegmatite-hosted deposits is diagnostic of petrogenetic differences, including melt source characteristics, magma evolution, or variable degrees of volatile solubility. The Li pegmatite compilation is suitable for fitting robust numerical models to support quantitative assessments. More well-defined Rb and Cs pegmatite resources are required for quantitative assessments, but these data provide useful information about original in-place resources for framing supply discussions.

Ore Geology Reviews↗
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↗
USGS1998

Water resources of Mellette and Todd counties, South Dakota

Mellette and Todd Counties are located in south-central South Dakota and have a combined area of 2,694 square miles. The White River and its tributaries, which include the Little White River, drain Mellette County and about one-half of Todd County. Tributaries to the Niobrara River, which include the Keya Paha River, drain the other one-half of Todd County. The average discharge of the Little White River is about 56 cubic feet per second as the river enters Todd County and is about 131 cubic feet per second as it discharges to the White River in northern Mellette County. The average discharge of the Keya Paha River just outside Todd County is about 39 cubic feet per second. The average annual runoff for Mellette and Todd Counties ranges from 0.94 to 2.36 inches based on records from nine streamflow-gaging stations in and near the counties. The average annual runoff is 1.62 inches, which compares with the average annual precipitation of about 19 inches. In Todd County, shallow wells completed in the alluvial, Ogallala, Arikaree, and White River aquifers generally can supply water that has low concentrations of dissolved solids, is fresh, and is soft to moderately hard. Ground water from shallow aquifers is limited in Mellette County; therefore, deep wells, often greater than 1,000 feet, are sometimes installed. The Pierre Shale often is used to supply rural domestic and stock wells in Mellette County even though well yields are low and the water has high dissolved solids, is moderately saline, and is very hard. Alluvial aquifers are present in both counties and store an estimated 1.6 million acre-feet of water. The water quality of the alluvial aquifers is dependent on the underlying deposits, and generally the water has low concentrations of dissolved solids, is fresh, and is soft to moderately hard where underlain by the Ogallala and Arikaree Formations; has moderate concentrations of dissolved solids, is slightly saline, and is hard where underlain by the White River Group; and has high concentrations of dissolved solids, is saline, and is very hard where underlain by the Pierre Shale. Also, yields often are lower where the alluvial aquifers are underlain by the Pierre Shale. The Ogallala aquifer is present in only Todd County, and the Arikaree aquifer is present throughout most of Todd County and southwestern and south-central Mellette County. The Ogallala aquifer contains an estimated 17 million acre-feet of water in storage, and the Arikaree aquifer contains an estimated 50 million acre-feet of water in storage. Both aquifers generally are suitable for irrigation, and yields from these aquifers are sometimes greater than 1,000 gallons per minute. Nitrate concentrations in 13 out of 92 water samples collected from the Ogallala aquifer exceeded the Primary Drinking Water Maximum Contaminant Level (MCL) of 10 milligrams per liter. In 11 out of 46 samples collected from the Arikaree aquifer, arsenic concentrations exceeded the MCL of 50 micrograms per liter. The White River aquifer, where present, is usually the shallowest source of ground water in Mellette County. The White River aquifer also is used in northern Todd County where the Ogallala and Arikaree aquifers are not present. The White River aquifer contains an estimated 50 million acre-feet of water in storage. Reported yields from the aquifer range from 1 to 30 gallons per minute, which generally is insufficient to support irrigation in most areas. However, yields are sufficient for livestock-watering and rural-domestic purposes. In both counties, the Pierre Shale is the shallowest bedrock aquifer and is exposed at the land surface throughout most of Mellette County. This aquifer is used primarily in Mellette County. Although the aquifer contains an estimated maximum of 1.5 million acre-feet of water in storage, it is not a viable source of ground water because the aquifer is relatively impermeable, yields are low, and water usually can be obtained from shallower sources, especially in Todd County. Reported yields from the Pierre Shale aquifer range from 1 to 8 gallons per minute. Because few test holes and wells penetrate below the Pierre Shale, little is known about the extent of the deeper bedrock aquifers. All wells completed in the Dakota Sandstone, Inyan Kara, and Minnelusa and Madison aquifers in the counties are used for stock-watering purposes. High concentrations of dissolved solids and hard water are characteristic of the water quality in the bedrock aquifers. Depths to the top of the deeper bedrock aquifers range from 1,270 feet to greater than 2,000 feet below land surface.

South Dakota↗
USGS2026

Computing flow-field distortion coefficients from well-construction and formation properties

Direct measurements of groundwater velocity made with borehole flowmeters in screened wells must be compensated for the effects of flow-field distortion (also known as borehole acceleration). A theoretical equation developed by Drost et al. (1968) and simple inputs describing hydraulic properties of well construction and geologic formation were programmed into an Excel workbook to facilitate computation by groundwater-flowmeter users. Tables describing the physical and hydraulic properties for well constructions and gravel pack media are provided with an example to facilitate use of the workbook. Groundwater flowlines converge or diverge as they pass from a geologic formation, through a gravel pack and well screen. The extent of flowline convergence or divergence and the value of the flow-field distortion coefficient is related to the relative changes in hydraulic conductivity of the well screen, gravel pack, and geologic formation. Convergence or divergence is accompanied by acceleration or deceleration of groundwater. Direct measurements of groundwater velocity at the center of the monitoring well can be adjusted to provide a more accurate estimate of velocity in the formation by applying a correction for flow-field distortion. Variables required to compute the flow-field distortion coefficient include the hydraulic conductivity of the gravel pack, well screen, and the geologic formation surrounding the well screen; the borehole radius, and the inside radius and outside radius of the well screen.

Groundwater↗
USGS2026

Distributed volcanic fields—Patterns, products, processes, and plumbing

Distributed volcanic fields—also referred to as “monogenetic fields” or simply “volcanic fields”—are found on every continent, on the seafloors, and in every tectonic setting on Earth. Yet, our understanding of what constitutes a distributed volcanic field, and our knowledge of the common products, spatiotemporal patterns in activity, and plumbing systems of these fields, are commonly poor and (or) debated. Here we provide an overview of the current understanding of distributed volcanic fields (DVFs) through review and compilation of published data, with an emphasis on understanding the characteristics of fields that are essential to prepare for future hazardous volcanic activity. Distributed fields can erupt restricted to variable magma compositions. Activity in these fields produces nearly all types of volcanoes, including scoria and spatter cones, maars, low shields, lava fields, stratovolcanoes, and calderas, but monogenetic volcanoes are the defining characteristic of DVFs. The number of vents and lifespans of DVFs vary greatly; vent numbers range from 2 to greater than 1,000, and fields can remain active for 150,000 to 10 million years. Some characteristics of DVFs appear field dependent. Although the volumes of single monogenetic eruptions within DVFs globally are variable (from less than 0.1 to a few cubic kilometers), each field may have a characteristic range in eruption volumes. Eruptions within distributed volcanic fields can have a random distribution and timing, but spatial and (or) temporal vent clustering and alignment are common, and eruptions may be time or volume predictable. Average eruption recurrence intervals in DVFs span over an order of magnitude, from less than 3 to greater than 100 thousand years. Exposures of the shallow portions of older systems reveal that many magma batches end their ascent before reaching the surface, with some dikes stopping within only meters of the surface. Petrologic insights into the deeper plumbing systems of DVF volcanoes indicate variable depths of magma crystallization and (or) stalling, and that magma transport and storage systems may be unique to a particular DVF and its underlying crustal structure. Both seismic and petrologic data for mafic DVF eruptions corroborate intrusions of magma into storage systems in the months to years before eruption, as well as variable ascent rates through the crust (days to weeks). Few distributed fields are closely monitored, as they erupt infrequently, but these observations imply that the amount of time available to install monitoring instrumentation after the first felt seismicity may be short, as precursory activity is frequently on the order of days to months for modern DVF eruptions. Satellite-based geodesy may provide useful warning in such cases. This chapter concludes with outstanding questions and areas for future research that will aid in our understanding and preparedness for future DVF eruptions.

Professional Paper↗
USGS2026

High‐resolution surface deformation and slip distribution observations for the 2023 Kahramanmaraş, Türkiye, earthquake sequence help constrain the rupture process

Splay, or branch, faults are a common geometric feature of earthquake surface ruptures and may provide constraints on the rupture behavior of an earthquake. The 2023 M w 7.8 Pazarcık and M w 7.5 Elbistan, Türkiye, earthquakes are examples of ruptures with multiple small splays, and the Pazarcık earthquake nucleated on a splay fault, the Narlı fault, before rupturing bilaterally on the East Anatolian fault (EAF). Here, we present 3‐m‐resolution surface displacement from subpixel correlation of Planet Dove optical images for the entirety of both ruptures with corresponding surface slip distributions. For a 30‐km‐long study region spanning the Narlı‐EAF intersection, we compare surface slip derived from five data sets with different resolutions (on‐the‐ground, WorldView, Planet Dove, Sentinel‐2, and Sentinel‐1) to elucidate complementary information. In addition, we integrate information from the surface expression of faulting with published dynamic rupture simulations and rupture process studies to constrain a rupture evolution for the Pazarcık earthquake that is consistent across data sets. This work highlights the complementary nature of disparate surface slip data sets and the role that high‐resolution surface displacement information, including from fault splays, can play in constraining nonunique rupture models and refining understanding of the earthquake rupture process.

Bulletin of the Seismological Society of America↗
USGS2023

A comprehensive fault system inversion approach: Methods and application to NSHM23

We present updated inversion‐based fault‐system solutions for the 2023 update to the National Seismic Hazard Model (NSHM23), standardizing earthquake rate model calculations on crustal faults across the western United States. We build upon the inversion methodology used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to solve for time‐independent rates of earthquakes in an interconnected fault system. The updated model explicitly maps out a wide range of fault recurrence and segmentation behavior (epistemic uncertainty), more completely exploring the solution space of viable models beyond those of UCERF3. We also improve the simulated annealing implementation, greatly increasing computational efficiency (and thus inversion convergence), and introduce an adaptive constraint weight calculation algorithm that helps to mediate between competing constraints. Hazard calculations show that ingredient changes (especially fault and deformation models) are the primary driver of hazard changes between NSHM23 and UCERF3. Updates to the inversion methodology are also consequential near faults in which the slip rate in UCERF3 was poorly fit or was satisfied primarily using large multifault ruptures that are now restricted by explicit b ‐value and segmentation constraints.

Bulletin of the Seismological Society of America↗
USGS2025

Timing of hydrothermal alteration and Au-Sb-W mineralization, Stibnite-Yellow Pine district, Idaho

The Stibnite-Yellow Pine district of central Idaho was mined from the early 1900s until the 1990s, extracting gold, antimony, tungsten, and mercury from veins and disseminated and replacement ores in mountainous terrain along the headwaters of the Salmon River. Mining during the two World Wars supplied critical antimony and tungsten to the war efforts. Recent exploration has delineated mineral resources of over 187 metric tons Au, 274 metric tons Ag, and 93,000 metric tons Sb. Mineralization is hosted in Cretaceous Idaho batholith granitic rocks and a sequence of Neoproterozoic to Paleozoic metasedimentary strata of carbonate and siliciclastic compositions. Historical studies outlined some of the complex paragenesis but debated the absolute age of mineralization. New petrographic and geochronologic work documents a sequence of five hydrothermal events in the Stibnite-Yellow Pine district. Event 1 is related to Cretaceous magmatic and hydrothermal activity and includes events ranging in age from 86 to 75 Ma, including sparse quartz-molybdenite veins dated at 86 Ma. Disseminated gold mineralization of event 2 is associated with sericitic alteration and sulfidation of igneous biotite and replacement of plagioclase by potassium feldspar, largely in granodiorite. Gold is present in zoned arsenian pyrite in both disseminated ores and in crosscutting carbonate-quartz veins containing pyrite and arsenopyrite. The large Yellow Pine deposit, localized at a dilatant bend in the Meadow Creek fault, hosts such disseminated and vein gold. Event 2 is interpreted as the major gold-forming event; 40 Ar/ 39 Ar ages of sericite and potassium feldspar alteration range, respectively, from 70 to 59 and 66 to 56 Ma. The long span is interpreted to reflect the age of gold mineralization and local overprinting by event 3. A narrower range from 66 to 61 Ma is interpreted to date the peak of gold mineralization and alteration. Event 3, tungsten mineralization with scheelite, is texturally later than event 2 gold and localized along the Meadow Creek structure. Event 3 scheelite has been dated by isotope dilution-thermal ionization mass spectrometry (ID-TIMS) and laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) U-Pb methods at 57 Ma. Event 4, best developed in the West End area, includes gold-silver–bearing quartz-carbonate-pyrite veins and breccias with epithermal textures and potassium feldspar alteration envelopes. Adularia from event 4 yields 40 Ar/ 39 Ar plateau ages of 52 to 51 Ma. Event 5 antimony and mercury mineralization consists of stibnite veins and breccia cements at the Yellow Pine and Hangar Flat deposits as well as cinnabar veins and replacements at the peripheral Fern and Hermes deposits; it is constrained by an LA-ICP-MS U-Pb date on scheelite (ca. 47 Ma) intergrown with stibnite. Minor propylitic and argillic alteration is evident in 47 Ma igneous dikes, which do not contain economic mineralization. The Au-Sb-W ores in the Stibnite-Yellow Pine mining district formed over an extended time period from about 70 to 45 Ma in multiple pulses that were localized along the Meadow Creek fault zone. Each event corresponds to episodes of magmatism and/or hydrothermal activity in the region. Insignificant event 1 skarn and molybdenum mineralization is similar in age to the Thompson Creek porphyry molybdenum deposit in central Idaho. Event 2 gold mineralization occurred during a magmatic gap in central Idaho but was synchronous with magmatism in the Bitterroot lobe further north; event 2 is similar in age to orogenic gold-arsenic mineralization at the Beartrack mine in eastern Idaho. Event 3 scheelite mineralization coincides with tungsten mineralization at the Quartz Creek deposit, late magmatism in the Bitterroot lobe, and rapid exhumation of the Atlanta lobe of the Idaho batholith. Event 4 gold mineralization is coincident with the onset of regional Challis magmatism and extension. Event 5 antimony and mercury mineralization is time-equivalent to epithermal gold mineralization in the nearby Thunder Mountain volcanic field and the peak of Challis magmatism.

Idaho↗