USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers in Marine Science”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,606 records · Page 49Linked to original sources
USGS2026

Spatial variation in effectiveness and costs of sagebrush restoration treatments across the western United States

Natural resource managers strive to improve restoration efficiency across the vast network of public lands of the United States, seeking the best return on investment through innovation, geographic prioritization, and adaptive management. This challenge is exemplified by ecosystems dominated by sagebrush ( Artemisia spp.), which occur across the western United States and have been degraded by increasing wildfire frequency driven by invasive plant spread, among other factors. However, both costs and effectiveness of sagebrush restoration can vary spatially due to management practices in addition to biotic and abiotic factors, and characterizing this variation could inform broad-scale planning. We leveraged published models of sagebrush recovery and treatment costs to predict and evaluate cost-effectiveness of aerial and drill seeding Artemisia spp. across 429 718 km 2 that overlapped 12 US states. Compared to natural recovery, effectiveness for both methods was generally low, yet effectiveness was often greater, and costs were always greater, for drill than aerial seeding, resulting in slightly higher cost-effectiveness on average for drill seeding. Cost-effectiveness for both treatments increased near major roads and in areas with more repeated burns. Cost-effectiveness also increased with predisturbance cover of sagebrush but decreased with herbaceous cover. However, we also identified areas where aerial seeding was more cost-effective than drill seeding, particularly with greater slope. Our results provide spatially explicit estimates and potential mechanisms of cost-effectiveness for two common seeding methods for sagebrush, which can help prioritize limited resources, guide land use, and improve restoration efficiency and effectiveness across public rangelands of the western United States.

western United States↗
USGS2025

Mapping eelgrass (Zostera marina) cover and biomass at Izembek Lagoon, Alaska, using in-situ field data and Sentinel-2 satellite imagery

The U.S. Geological Survey and the U.S. Fish and Wildlife Service have developed a three-tiered strategy for monitoring eelgrass ( Zostera marina ) beds at Izembek Lagoon, Alaska, that targets different spatial and temporal scales. The broadest-scale monitoring (tier-1) uses satellite imagery about every 5 years to delineate the spatial extent of eelgrass beds throughout the lagoon. This report describes the most recent (mid-2020s) tier-1 eelgrass monitoring at Izembek Lagoon. The monitoring effort began by canvasing all satellite imagery collected during summer, under clear daytime skies and at low-tide, since the last tier-1 effort in 2006. Two eelgrass maps of Izembek Lagoon were generated by first creating maps of spectrally unique classes from two Sentinel-2 satellite images collected on July 1, 2016, and August 14, 2020, then attributing those spectral classes with information about eelgrass conditions based on field data. Specifically, maps depicting various eelgrass metrics, such as percentage of cover and modeled biomass, were generated using summaries of the ground data that spatially intersected each spectral class. Comparisons of the 2016 and 2020 Sentinel-2 maps showing eelgrass distributional extent, as well as a 2006 Landsat map, indicated that areas where eelgrass presence may have declined during 2006–20 were most prevalent in the central part of Izembek Lagoon. More recently, during 2016-20, areas of possible biomass decline were more prevalent in the southern part of the lagoon. Monitoring eelgrass conditions at Izembek Lagoon with satellite imagery and concurrent ground data allows conditions to be compared over time, but the influences of tide levels, growing season phenology, and spatiotemporal co-registration accuracy should be considered when designing and interpreting change detection analyses.

Alaska↗
USGS2025

Hydrogeology, water budget, and simulated groundwater availability in the Salt Fork Arkansas River and Chikaskia River alluvial aquifers, northern Oklahoma, 1980–2020

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s aquifers to determine the maximum annual yield for each groundwater basin. The U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an updated hydrologic investigation of the Salt Fork Arkansas River and Chikaskia River alluvial aquifers in northern Oklahoma for the study period spanning 1980–2020 and evaluated the simulated effects of potential groundwater withdrawals on groundwater flow and availability in the Salt Fork Arkansas River alluvial aquifer. A hydrogeologic framework and conceptual model were developed to guide the development of a numerical model. Three groundwater-availability scenarios were evaluated by using the calibrated numerical model, which was focused on the Salt Fork Arkansas River alluvial aquifer. These scenarios were used to (1) estimate equal-proportionate-share groundwater withdrawal rates, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought. The 20-, 40-, and 50-year equal-proportionate-share groundwater withdrawal rates for the Salt Fork Arkansas River alluvial aquifer under normal recharge conditions were about 0.63, 0.58, and 0.57 acre-foot per acre per year, respectively. Projected 50-year groundwater withdrawal scenarios were used to simulate the effects of modified well withdrawal rates. Because well withdrawals were less than 2 percent of the calibrated numerical-model water budget, changes to the well groundwater withdrawal rates had little effect on simulated Salt Fork Arkansas River base flows and groundwater storage in the Salt Fork Arkansas River alluvial aquifer. A hypothetical 10-year drought scenario was used to simulate the potential effects of a prolonged period of reduced recharge on groundwater storage. Groundwater storage at the end of the hypothetical drought period was 14.5 percent less than the groundwater storage of the calibrated numerical model without the simulated drought.

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

Evaluation of nutrient, alkalinity, and acid-neutralizing capacity stabilities in water samples analyzed by the U.S. Geological Survey National Water Quality Laboratory, 2023–24

The U.S. Geological Survey evaluated the stability of water-sample chemical analysis of nutrient, alkalinity, and acid-neutralizing capacity constituents with respect to the duration between sample collection and laboratory analysis, also known as the sample holding time. A study began in the spring of 2023 to evaluate the sample stability, between 2 and 180 days after sample collection, of the chemical properties and chemical constituents of alkalinity as calcium carbonate, filtered; acid-neutralizing capacity as calcium carbonate, unfiltered; total ammonia as nitrogen, filtered; total ammonia plus organic nitrogen as nitrogen, filtered and unfiltered; nitrite as nitrogen, filtered; nitrate plus nitrite as nitrogen, filtered; total nitrogen, filtered and unfiltered; orthophosphate as phosphorous, filtered; and total phosphorus as phosphorus (filtered and unfiltered) in water. Both surface water and groundwater matrices were represented. Sample instability varied by observed property and matrix; therefore, providing general guidance for sample holding time is not possible based on matrices alone. No correlations between field measurements of sample characteristics and sample instability were observed. Although observations for some properties indicate sample stability that exceeds the recognized U.S. Geological Survey National Water Quality Laboratory method holding times, this is not necessarily the case for matrices and seasonal characteristics that were not investigated. Based on the limited number of six sample sources used in this study, some patterns emerge for the 12 observed properties studied. Five observed properties generally indicate stability for as many as 180 days after sampling (total nitrogen as nitrogen, both filtered and unfiltered; orthophosphate as phosphorus, filtered; and phosphorus as phosphorus, both filtered and unfiltered). Other observed properties indicate stability for as many as 180 days for some matrices, but not for others. Finally, some observed properties indicate instability well before 180 days.

Scientific Investigations Report↗
USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗
USGS2024

Deep structure of Siletzia in the Puget Lowland: Imaging an obducted plateau and accretionary thrust belt with potential fields

Detailed understanding of crustal components and tectonic history of forearcs is important due to their geological complexity and high seismic hazard. The principal component of the Cascadia forearc is Siletzia, a composite basaltic terrane of oceanic origin. Much is known about the lithology and age of the province. However, glacial sediments blanketing the Puget Lowland obscure its lateral extent and internal structure, hindering our ability to fully understand its tectonic history and its influence on modern deformation. In this study, we apply map-view interpretation and two-dimensional modeling of aeromagnetic and gravity data to the magnetically stratified Siletzia terrane revealing its internal structure and characterizing its eastern boundary. These analyses suggest the contact between Siletzia (Crescent Formation) and the Eocene accretionary prism trends northward under Lake Washington. North of Seattle, this boundary dips east where it crosses the Kingston arch, whereas south of Seattle the contact dips west where it crosses the Seattle uplift (SU). This westward dip is opposite the dip of the Eocene subduction interface, implying obduction of Siletzia upper crust at this southern location. Elongate pairs of high and low magnetic anomalies over the SU suggest imbrication of steeply-dipping, deeply rooted slices of Crescent Formation within Siletzia. We hypothesize these features result from duplication of Crescent Formation in an accretionary fold-thrust belt during the Eocene. The active Seattle fault divides this Eocene fold-thrust belt into two zones with different structural trends and opposite frontal ramp dips, suggesting the Seattle fault may have originated as a tear fault during accretion.

Washington↗
USGS2024

Generalized potentiometric maps of the Fort Union, Hell Creek, and Fox Hills aquifers within the Standing Rock Reservation

Generalized potentiometric surfaces of the Fort Union, Hell Creek, and Fox Hills aquifers were constructed to assess the groundwater resources of the Standing Rock Reservation. Additionally, this information can provide water managers with tools and data to effectively manage water resources in the future. Previous studies that mapped the geology and hydrogeology of the area at differing scales were used to confirm in which aquifer the study wells were completed. Water-level data from wells are provided by the U.S. Geological Survey Groundwater System Inventory database, the South Dakota Department of Agriculture and Natural Resources, and the North Dakota Department of Water Resources. Hydrographs were constructed for five selected observation wells to evaluate historical water-level fluctuations and trends. Hydrographs for the Hell Creek aquifer showed a flat trend with a rise in 2020. Hydrographs for the deeper Fox Hills aquifer showed that water levels fluctuated in response to climatic conditions and demonstrated an increasing trend in water-level elevations starting in 2010. Hydrographs were not constructed for any wells completed in the Fort Union Formation because none of the wells had continuous long-term measurements. Generalized potentiometric surfaces, constructed from interpolating water-level elevations, gave insight into groundwater flow directions. Groundwater in the Fort Union aquifer likely flows radially outward from the northwest part of the study area to the northeast and south-southeast parts. Groundwater in the Hell Creek aquifer generally flows from higher elevations in the northwest towards lower areas, where surface-water tributaries have incised into the aquifer. Groundwater in the Fox Hills aquifer likely flows from higher elevations in the west, southwest, and central parts of the study area towards the valleys of the Grand River and Missouri River. Most wells used for constructing potentiometric maps had only one recorded water-level measurement from drillers at the time of well construction. These measurements are often subject to error because the well is still recovering and because of spatial limitations of data availability. Also, because single water-level measurements were recorded at different points in time, additional uncertainty is introduced by fluctuating climatic conditions effect on water levels. Potentiometric map interpretation limitations are a result of areas with sparse data. Limitations also arise from potentiometric surfaces generalizing a complex and dynamic hydrogeologic system; however, the generalized potentiometric surface maps can be used to assist water managers and can help prioritize locations for future monitoring in areas with high uncertainty from sparse existing data.

North Dakota, South Dakota↗
USGS2024

Post-glacial stratigraphy and late Holocene record of great Cascadia earthquakes in Ozette Lake, Washington, USA

Ozette Lake is an ~100-m-deep coastal lake located along the outer coast of the Olympic Peninsula (Washington, USA); it is situated above the locked portion of the northern Cascadia megathrust but also relatively isolated from active crustal faults and intraslab earthquakes. Here we present a suite of geophysical and geological evidence for earthquake-triggered mass transport deposits (MTDs) and related turbidite deposition in Ozette Lake since ca. 14 ka. Comprehensive high-resolution bathymetry data, seismic reflection profiles, and sediment cores are used to characterize the post-glacial stratigraphic framework and examine paleoseismic evidence in the lacustrine sediments. Stacked sequences of MTDs along the steep eastern flanks of the lake appear to grade basin-ward from thick, chaotic, blocky masses to thin, parallel-bedded turbidite beds. The discrete turbidite event layers are separated by fine-grained (silt and clay) lake sedimentation. The event layers are observed throughout the lake, but the physical characteristics of the deposits vary considerably depending on proximity to primary depocenters, steep slopes, and subaqueous deltas. A total of 30–34 event deposits are observed in the post-glacial record. Radiometric dating was used to reconstruct a detailed sedimentation history over the last ~5.5 k.y., develop an age model, and estimate the recurrence (365–405 yr) for the most recent 12 event layers. Based on sedimentological characteristics, temporal overlap with other regional paleoseismic chronologies, and recurrence estimates, at least 10 of the dated event layers appear to be sourced from slope failures triggered by intense shaking during megathrust ruptures; the recurrence interval for these 10 events is 440–560 yr. Thus, Ozette Lake contains one of the longest and most robust geological records of repeated shaking along the northern Cascadia subduction zone.

Washington↗
USGS2024

Time varying crustal anisotropy at Whakaari/White Island volcano

Whakaari/White Island has been the most active New Zealand volcano in the 21st century, producing small phreatic and phreatomagmatic eruptions, which are hard to predict. The most recent eruption occurred in 2019, tragically claiming the lives of 22 individuals and causing numerous injuries. We employed shear-wave splitting analyses to investigate variations in anisotropy during quiescence, unrest, and an eruption between 2018 and 2020. We examined spatial and temporal variations in 7972 shear-wave splitting measurements and derived 1402 V p /IV s ratio measurements. Comparing shear-wave splitting parameters from similar earthquake paths across different times indicates that the observed temporal changes are unlikely to result from variations in earthquake paths through media with spatial variability. Instead, these changes may stem from variations in anisotropy over time, likely caused by changes in crack alignment due to stress or varying fluid content.

Whakaari/White Island volcano↗
USGS2025

Identifying precursors and tracking pulses of magma ascent in multidisciplinary data during the 2018–2023 phreatomagmatic eruption at Semisopochnoi Island, Alaska

The 2018–2023 phreatomagmatic eruptions at Semisopochnoi Island, Alaska produced abundant long-period (LP) seismicity, harmonic and broadband tremor, and explosion signals over several well-monitored periods of eruption and quiescence. The corresponding dataset provides an excellent opportunity to investigate precursory and syn -eruptive geophysical signals of long-lived phreatomagmatic eruptions using multiparameter observations. We generated explosion and LP event catalogs through novel implementations of the REDPy ( Hotovec-Ellis, 2024 ) repeating event detector in mid-2021 following a network upgrade and the onset of a new phase of the eruption. The hundreds of detected explosions show a high degree of infrasound waveform similarity over more than a year, indicating a repeating source mechanism likely associated with explosive magma-water interaction. The seismic LP catalog shows that events began over a month prior to renewed explosive activity at the beginning of August 2021, and that lower frequency index (FI) LPs were generated in the week prior to the onset of explosions. We applied a recently developed machine learning tool (VOISS-Net, Tan et al., 2024 ) to catalog abundant broadband and harmonic seismic tremor recorded before and during the renewed explosive activity, along with LPs and explosions. The tremor catalogs complement the LP and explosion catalogs by filling out the seismic sequence with the dominant signal types. Together, these catalogs reveal a seismic sequence of renewed unrest that started with several weeks of LP events, followed by LPs with lower FI values and harmonic tremor in the days prior to explosive activity, and finally the onset of discrete explosions and broadband eruption tremor. We interpret this sequence as the ascent of a new pulse of magma that first interacted with the hydrothermal/groundwater system to produce LPs, followed by harmonic tremor, and that ultimately drove explosive magma-water interactions and periods of continuous ash emissions. The 2021 seismic sequence, in combination with long-term records of satellite SO 2 emissions, deformation from interferometric synthetic aperture radar (InSAR) analysis, ash sample analysis, infrasound, and volcano tectonic seismicity, allows us to interpret the entire 9-year period of unrest and eruption that began with an intrusion and earthquake swarm in 2014.

Alaska↗
USGS2025

Evolution of permeability and strength recovery of shear fracture under hydrothermal conditions

Geothermal energy is a clean and renewable resource that depends on the ability to move water through hot rock. In many locations, the ability to move water through rock requires the presence of extensive natural or human-made systems of fractures. However, these fracture systems are influenced by a variety of complex processes that occur at the temperature and pressure conditions found in geothermal reservoirs. To make geothermal systems more efficient, it is important to understand how these fractures evolve over time in response to these processes. This study explored how fractures in rock change under high-temperature and pressure conditions, like those found in geothermal reservoirs, through laboratory experiments and numerical simulations. In general fractures, tend to close due to pressure from the surrounding rock preventing the flow of water. One set of laboratory experiments revealed that fractures in granitic rock can weaken over time at elevated temperatures. These weak fractures can then slip slowly potentially keeping the fracture open for a longer period. This study also showed that at the highest temperature examined (250 °C) large differences between the chemical composition of the rock and the fluid moving through it could improve water flow temporarily likely due to the water dissolving parts of the rock. Flow rate predicted by models shows overall good agreement with the flow rate measured in our experiments but with some long time-scale variations that may be due to mechanical wear of the fracture surface. A second set of experiments found that, over the temperature range examined, the closure of fractures was caused by brittle failure driven by high stresses at points where fracture surfaces were in contact. This result suggests that fracture closure in geothermal reservoirs could be counteracted by maintaining high pore pressures. Additionally, we also tested a new method using electrical resistance to measure fracture closure in controlled systems with high precision.

Final Technical Report↗
USGS2025

Rapid emplacement of the Keaiwa Lava Flow of 1823 from the Great Crack in the Southwest Rift Zone of Kilauea volcano

The Keaīwa Lava Flow of 1823 in the Southwest Rift Zone of Kīlauea volcano is unusual for its expansive pāhoehoe sheet flow morphology and lack of constructive vent topography, despite having a similar tholeiitic basalt composition to other lavas erupted from Kīlauea. This lava flow issued from a ∼10-km-long continuous fissure now known as the Great Crack, and has an unusually thin sheet flow morphology with margin thicknesses of ∼15–110 cm (average of 42 cm). Based on field observations of the lava flow at its fissure vent (e.g., drain-back features), we propose that the Great Crack formed, or at least significantly widened, just prior to and syn-eruptively with this 1823 eruption. The absence of pyroclastic cones or spatter ramparts indicates that the eruption consisted of a rapid outpouring of relatively degassed lava as the fissure unzipped. The rapidly moving lava flow overtopped pre-existing tumuli and scoria cones (e.g., Lava Plastered Cones) up to ∼10 m tall. Glass and whole-rock chemistry yield homogeneous compositions for the lavas erupted from the Great Crack, with glass compositions of 6.40 ± 0.10 wt% MgO and whole-rock compositions of 7.39 ± 0.07 wt% MgO. Lava pads erupted from a short western fissure system are richer in mafic minerals (e.g., olivine and clinopyroxene), and show slightly more MgO-rich whole-rock compositions (7.79 ± 0.05 wt%). MgO-in-glass thermometry on juvenile spatter yield eruption temperatures of 1153 ± 13°C that are typical of Kīlauea lavas. Thus, the extensive sheet-like lava flow morphology is not a direct consequence of unusual magmatic or rheological conditions (i.e., low viscosity). Instead, the flow morphology is associated with high effusion rates caused by sudden drainage of uprift magma as it erupted from the Great Crack. Lava flow modeling on a 2-m-resolution digital elevation model indicates that a minimum bulk effusion rate of ∼5800 m 3 /s (∼3500 m 3 /s dense rock equivalent) and a minimum flow velocity of ∼11 m/s are required for the lava flow to overcome the topography of the Lava Plastered Cones. This effusion rate is among the highest inferred for eruptions in Hawaiʻi and around the world. This study highlights a less frequent eruption style at Hawaiian volcanoes characterized by a sudden outpouring of lava from an unusual fissure system. Local eyewitness accounts indicate that the 1823 eruption was preceded by seismicity. Given the complex magmatic-volcanic-tectonic relations across Kīlauea, we speculate that the south flank could have slipped over one or more events that ultimately triggered unzipping of the Great Crack and passive release of briefly stored uprift magma. An eruption similar to 1823 at Kīlauea or Mauna Loa, with an eruptive timeframe that could be as short as an hour, with high effusion rates and rapid flow front velocities, would not easily allow for a timely response.

Hawaii↗
USGS2024

Groundwater flow model for the Des Moines River alluvial aquifer near Des Moines, Iowa

Des Moines Water Works (DMWW) is a regional municipal water utility that provides residential and commercial water resources to about 600,000 customers in Des Moines, Iowa, and surrounding municipalities in central Iowa. DMWW has identified a need for increased water supply and is exploring the potential for expanding groundwater production capabilities in the Des Moines River alluvial aquifer, where it operates two radial collector wells (RCWs). The U.S. Geological Survey, in cooperation with DMWW, completed a study of the Des Moines River alluvial aquifer and interactions of the RCWs with the aquifer; no previously published model has included the existing well locations, which is the focus of this model. A conceptual and numerical groundwater flow model have been developed to characterize the Des Moines River alluvial aquifer under existing conditions, to simulate water levels observed in the RCWs, and to provide publicly accessible hydrologic data and research that advance understanding of the regional hydrologic system and can potentially be used in the future to evaluate groundwater production scenarios. Model performance was assessed by comparing observed and simulated groundwater levels that included water level elevations, water level changes, water level inequality observations, surface water streamflow, and change in surface water volume from upstream to downstream. Water table elevation in the aquifer layers is on average slightly overestimated with average absolute value error less than 1.5 meters at both RCWs and less than 2.5 meters for all observation wells in the alluvial aquifer layers. The model also accurately simulated water tables greater than the RCW design minimum (a water level threshold at which RCW pumping is reduced) in all timesteps for which water level observation data existed. Water table elevation error was higher in other model layers that were not the focus of the study, and the model did not accurately match streamflow targets.

Iowa↗
USGS2025

Metal-rich lacustrine sediments from legacy mining perpetuate copper exposure to aquatic-riparian food webs

Historic copper mining left a legacy of metal-rich tailings resulting in ecological impacts along and within Torch Lake, an area of concern in the Keweenaw Peninsula, Michigan, USA. Given the toxicity of copper to invertebrates, this study assessed the influence of this legacy on present day nearshore aquatic and terrestrial ecosystems. We measured the metal (Co, Cu, Ni, Zn, Cd) and metalloid (As) concentrations in sediment, pore water, surface water, larval and adult insects, and two riparian spider taxa collected from Torch Lake and a nearby reference lake. Overall, elevated metal and metalloid concentrations, particularly Cu, were measured in all sediment samples and some surface and pore water samples collected from Torch Lake. For instance, Cu concentrations in the Torch Lake sediment were ∼200% higher than the reference lake and all measured concentrations exceeded predicted effects concentrations by at least ninefold. Within larval insect tissues, we observed 160% higher Cu concentrations than measured in the reference lake, and Cu was the only measured element above predicted effects concentrations in Torch Lake. Adult insects collected at both lakes had similar metal concentrations irrespective of exposure levels. Yet we found 100% higher copper concentrations in Torch Lake riparian spiders, demonstrating elevated exposure risk to insectivores across the aquatic-terrestrial boundary. Our results highlight that other metals in the mixture may not be as concerning to adjacent riparian ecosystems, but copper remains a contaminant of concern in Torch Lake 60 years after mining ceased.

Michigan↗
USGS2026

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

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

Hawaii↗
USGS2026

Geochemical disequilibrium at the brittle-ductile transition

We investigate the microtextural, microchemical, and isotopic effects of late-stage ductile deformation in quartzite mylonites and kyanite–muscovite–quartz veins from the Raft River shear zone (Utah). Quartz microstructures record pervasive disequilibrium, expressed by unannealed features including undulatory extinction, deformation lamellae, and poorly defined fabrics, typical of waning deformation in shear zones. Microchemical disequilibrium is best preserved in kyanite quartzites, where CL-zoned kyanite records repeated fractures, overgrowth, and mineral precipitation, and in muscovite from muscovite-poor quartzite mylonites that shows minor-element zoning consistent with syn-deformational overgrowth on detrital cores. In contrast, muscovite from muscovite-rich kyanite quartzites exhibits minimal chemical zoning. These microchemical variations correlate with 40 Ar/ 39 Ar age systematics. Chemically zoned muscovite preserves variable single-step ages, including ∼150 Ma ages, reflecting retention of detrital cores. In contrast, syndeformational muscovite consistently yields Miocene ages, indicating recrystallization and new growth below argon closure temperatures that reset inherited isotopic signatures. Similar trends are observed in quartzite mylonites, where increasing quartz recrystallization and stronger crystallographic preferred orientations occur toward deeper structural levels. Together, these observations indicate increasing retrograde deformation and recrystallization with depth in the Raft River shear zone and demonstrate that strain-driven recrystallization exerts a first-order control on muscovite ages. We suggest that apparent thermochronologic gradients in retrograde shear zones may reflect recrystallization gradients rather than temperature gradients. Where cooling limits the thermal driving force for recrystallization, isotopic relics are preserved, and local deformation and fluid availability control re-equilibration. Consequently, isotopic disequilibrium—particularly in the 40 Ar/ 39 Ar system—may be the rule rather than the exception of retrograde tectonic environments.

Idaho, Utah↗
USGS2026

Potential impacts of groundwater pumping on stream temperature are greatest in streams with substantial cold groundwater inflows

Groundwater pumping-induced reductions in streamflow (known as ‘streamflow depletion’) have been documented worldwide, but potential impacts of streamflow depletion on stream temperature are not well understood. Here, we use two types of models to identify potential impacts of pumping on stream temperature across the conterminous United States (CONUS) to determine which aspects of a stream's annual thermograph (thermal signatures) can be used to monitor and manage streamflow depletion impacts on stream temperature. We used long-term streamflow and stream temperature data from 30 streamgages across CONUS and surrogate models of streamflow depletion to analyse potential stream temperature impacts at each site. We compared two different stream temperature modelling approaches: (i) a process-based energy balance model and (ii) statistical regression models based on air temperature and stream discharge. We calculated a suite of thermal signatures under depleted and non-depleted conditions for each stream and found that maximum annual 7-day temperature and annual temperature range are potentially the most sensitive to streamflow depletion, with potential changes of at least 2°C at > 70% of the sites when using the process-based model. We also found that the regression-based models predicted much less sensitivity of stream temperature to streamflow depletion than the process-based model. This work provides an initial evaluation and sensitivity analysis of the potential impacts of streamflow depletion on stream temperature. We demonstrate that stream temperature may be most sensitive to pumping in streams with a high proportion of flow sourced from relatively cold groundwater inputs, and that regression-based stream temperature models may underpredict stream temperature changes caused by streamflow depletion.

conterminous United States↗