Geology Reports⌕ Search

Geology topics

Research about United States

Source-linked reports with geographic coverage including United States.

1,087 records · Page 18Linked to original sources
USGS2025

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana↗
USGS2025

Pre-eruptive characteristics of “suspect” silicic magmas in Carlin-type Au-forming systems

World-class Carlin-type Au deposits hosted in sedimentary rock were formed when profuse Eocene silicic magmatism swept across northern Nevada in response to arc migration. Carlin-type Au deposits formed along with porphyry/skarn Cu-Mo-W-Au deposits, epithermal Ag-Au deposits, and distal disseminated Ag-Au deposits. But unlike these other Au-bearing deposits that have clear associations with igneous intrusions, Carlin-type ore deposits appear to have formed distant from concealed plutons, and their origin remains controversial. Despite decades of abundant geophysical, geochronological, and geochemical studies suggesting the involvement of magmas, concrete evidence for magmatic involvement is still lacking. Consequently, the involvement of contemporaneous igneous systems remains inferred based on age, proximity, and variable isotopic, geochemical, and geophysical clues. A recent synthesis of deposit models postulates that Carlin-type Au deposits are intrusion-related, but that the causative magmas reside deeper (∼6–12 km) than in typical porphyry and peripheral systems (∼3–5 km), meaning that Carlin-type deposits are perhaps more distal expressions of igneous intrusions. We investigate a collection of “suspect” magmatic systems over a ∼7 m.y. timespan (∼41–34 Ma) that are contemporaneous with and near known Carlin-type ore deposits. We report results of a multifaceted array of in situ geochemical analyses (FTIR, EMP, SHRIMP-RG, LA-ICP-MS) of quartz-hosted melt inclusions, biotite, and quartz to better characterize the pre-eruptive characteristics of these magmas. We also report results of thermobarometry and thermodynamic phase equilibria modeling to help place constraints on magmatic reservoir depths and processes. Rather than a single “flavor” of silicic magma, we observe a surprisingly broad compositional spectrum of rhyolites, with one end of the spectrum exhibiting more arc-like (I-type) characteristics and the other end displaying more post-subduction, thick-crust extensional (A-type) characteristics. This broad compositional spectrum suggests a more complex picture of silicic crustal magmatism operating over a narrow span of time during slab rollback. Despite this spectrum, magmatic systems in this study are consistently ferroan and generally peraluminous, which we interpret as an expression of the relatively elevated geotherm at the time and incorporation of variable amounts of highly peraluminous metasedimentary crustal components. The silicic magma spectrum encompasses a range of mineralization associations, including subduction-related Cu-Mo-W-Au-Ag and post-subduction, thick-crust extensional rare-metal Mo-Sn-W-F-Be-Ag-Au, consistent with the prolific and diverse array of ore deposits that formed during this time. Carlin-type Au deposition appears to be associated with nearly the entire magmatic spectrum. This apparent indifference to silicic magma “flavor” would seem to imply that if magmas are involved in Carlin-type Au deposit genesis, they perhaps do not need to be compositionally specialized and/or possibly are only relevant as heat sources driving circulation to remobilize and redistribute metals.

Nevada↗
USGS2025

Geochemistry and Soils of the Big Smoky Valley Fens, Nevada

Fens are groundwater-fed wetlands that can provide habitat for plants and animals. Due to anthropogenic activities and climate change, many fens around the world are at risk. This paper presents the results of a study of the hydrology and geochemistry of fens in Big Smoky Valley, central Nevada to support the Bureau of Land Management’s activities in the area. A water sample from the largest fen in the study area was analyzed for its water chemistry and compared to a nearby alluvial aquifer and hot spring. The high SiO 2 concentration of the fen sample implies that the fen water may originate from geothermal water. A soil core was taken to analyze radiocarbon age and soil type. A majority of the core was composed of silt and clay interlayered with water-filled voids. Changes in the character of the clay with depth suggest that there may have been changes in the depositional environment over time. Radiocarbon dating of Ruppia seeds showed longevity of the fen, with the minimum 14 C age of the core as 4,375±40 years. This paper provides reconnaissance-level information on the Big Smoky Valley fens, but further information would be needed to better understand the source of water to the fens or how the fen environment has changed over time with climate.

Nevada↗
USGS2025

Detecting hidden sedimentary geothermal systems in the Upper Colorado River Basin

Geothermal resources exist in sedimentary rock where circulation of water for efficient extraction or storage of heat is possible. Except in rare instances where hot water is expressed at the land surface, sedimentary geothermal resources are hidden, so the identification of these systems is optimally accomplished using predictive subsurface modeling. An integrated approach using detailed paleogeographic interpretations, subsurface geologic mapping, and numerical modeling has produced regional geologic and temperature models for the Upper Colorado River Basin, a large watershed in central North America that contains many sedimentary basins. These models identify areas of hidden sedimentary geothermal resource potential in low temperature (<90°C), moderate temperature (90–150°C), and high temperature (>150°C) fairways across the study area. These models incorporate maps of key horizons in outcrop and the subsurface to create a robust structural framework that can be used to target favorable geology for natural or engineered permeability. This framework is populated with lithologies derived from detailed palaeogeographical maps and over 40,000 bottom hole temperature (BHT) values were used to create a calibrated three-dimensional (3D) temperature model across the region. The resulting maps serve as a regional sedimentary geothermal play fairway screening tool for evaluating different grades of sedimentary geothermal resources and for identifying areas of interest where more detailed, prospect-scale studies can be undertaken.

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

Quality of groundwater used for domestic supply in the Gilroy-Hollister basin and surrounding areas, California, 2022

More than 2 million Californians rely on groundwater from domestic wells for drinking-water supply. This report summarizes a 2022 California Groundwater Ambient Monitoring and Assessment Priority Basin Project (GAMA-PBP) water-quality survey of 33 domestic and small-system drinking-water supply wells in the Gilroy-Hollister Valley groundwater basin and the surrounding areas, where more than 20,000 residents are estimated to utilize privately owned domestic wells. The study area includes the Llagas subbasin in the north, the North San Benito subbasin in the south, and the surrounding uplands. The study was focused on groundwater resources used for domestic drinking-water supply, which are mostly drawn from shallower parts of aquifer systems rather than those of groundwater resources used for public drinking-water supply in the same area. This assessment characterized the quality of ambient groundwater in the aquifer before filtration or treatment, rather than the quality of drinking water delivered to the tap. To provide context, the measured concentrations of constituents in groundwater were compared to Federal and California State regulatory and non-regulatory benchmarks for drinking-water quality. A grid-based method was used to estimate the areal proportions of groundwater resources used for domestic drinking wells that have water-quality constituents present at high concentrations (above the benchmark), moderate concentrations (between one-half of the benchmark and the benchmark for inorganic constituents, or between one-tenth of the benchmark and the benchmark for organic constituents), and low concentrations (less than one-half or one-tenth the benchmark for inorganic and organic constituents, respectively). This method provides statistically representative results at the study-area scale and permits comparisons to other GAMA-PBP study areas. In the study area, inorganic constituents in groundwater were greater than regulatory benchmarks (U.S. Environmental Protection Agency [EPA] or State of California maximum contaminant levels [MCLs]) for public drinking-water quality in 24 percent of domestic groundwater resources. The inorganic constituents present at concentrations greater than MCLs for drinking water were nitrate (as nitrogen), barium, chromium, and selenium. Total dissolved solids (TDS) or manganese were present at concentrations greater than the secondary maximum contaminant levels (SMCLs) that the State of California uses as aesthetic-based benchmarks in 48 percent of domestic groundwater resources. No volatile organic compounds or pesticide constituents were present at concentrations greater than regulatory benchmarks. Total coliform bacteria and enterococci were detected in 4 percent of domestic groundwater resources. Per- and polyfluoroalkyl substances (PFAS) were detected in 19 percent of domestic groundwater resources, and 10 percent had concentrations greater than recently enacted (April 2024) EPA MCLs. Physical and chemical factors from natural and anthropogenic sources that could affect the groundwater quality were evaluated using results from statistical testing of associations between constituent concentrations and potential explanatory variables. In this study, relevant physical factors include well construction characteristics, groundwater age, site proximity to groundwater recharge or discharge zones, and potential sources of contamination. Relevant chemical factors include the initial chemistry of the recharge water, the mineralogy of the aquifer sediments, and the subsequent shifts in chemistry as biologic and geologic reactions alter groundwater in the subsurface. Nitrate concentrations were correlated to agricultural land use, distance from the boundary of the Gilroy-Hollister Valley groundwater basin, and the proportion of modern (post-1950s) water captured by the well. Denitrification under anoxic redox conditions can mitigate some nitrate derived from fertilizer application. Total dissolved solids primarily were derived from water-rock interactions with soils and aquifer materials in the study area, but there were high concentrations where agricultural practices contributed additional TDS. Mineralogy of aquifer sediments and rocks also affect barium, selenium, boron, and chromium concentrations in the Gilroy-Hollister Valley groundwater basin. PFAS were positively correlated with urban land use and the proportion of modern water captured by the well.

California↗
USGS2025

Depth of magma crystallization and fluid exsolution beneath the porphyry-skarn Cu deposits at Santa Rita and Hanover-Fierro, New Mexico, USA

The depth level at which porphyry Cu–forming magmas fractionated and exsolved mineralizing fluids is actively debated. In the classic model, extensive magma fractionation occurs in large, upper crustal magma chambers, and concomitant fluid exsolution leads to forceful expulsion of residual magmas in the form of porphyry dikes, stocks, and breccia pipes, which subsequently serve as pathways for the mineralizing fluids. In contrast, some recent studies highlighting the role of deep crustal magma fractionation in the production of fertile magmas essentially deny the existence of upper crustal magma chambers at the time of mineralization. To address this, we conducted a detailed thermobarometric investigation of 13 intermediate to felsic, porphyritic intrusive rocks related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, United States, representing two premineralization magmas (61–60 Ma), seven synmineralization magmas (60–58 Ma), and four late- to postmineralization magmas (58–57 Ma). For each sample, the pressure of last magma crystallization before final magma ascent to the current exposure level was reconstructed based on Al-in-hornblende barometry of small hornblende inclusions trapped within quartz phenocrysts and through titanium-in-quartz (TitaniQ) thermobarometry of the host quartz phenocrysts themselves. Since quartz is one of the last crystallizing magmatic minerals, and no significant phenocryst growth could have occurred in small dikes and stocks after final magma emplacement, quartz phenocrysts and their contained hornblende inclusions record the depth of last magma crystallization before final magma ascent. When present, hornblende phenocrysts and hornblende inclusions within other major phenocrysts were also analyzed. Both quartz and hornblende barometers return consistent average pressures of 3.2 ± 0.4 kbar for the entire suite of pre- to postmineralization magmas, corresponding to depths of 11 to 14 km. The synmineralization magmas return even more consistent average pressures of 3.1 ± 0.2 kbar, corresponding to a depth of 12 ± 1 km. The volume of the mineralizing porphyry dikes and stocks at the emplacement level is far too small to have provided all the fluids and metals required to form the observed ore deposits. Therefore, the majority of the ore-forming fluids must have originated from the magmas that crystallized at 12 ± 1 km depth. The ore deposits, conversely, formed at ~5-km paleodepth. This implies that most of the mineralizing fluids traveled an average vertical distance of ~7 km from their magmatic source to the eventual site of ore precipitation. The relatively unaltered nature and low veining degree of deeper parts of mineralized porphyry dikes and stocks suggest that the fluid transport through these intrusive bodies occurred mostly at near-solidus conditions by means of fluid percolation along grain boundaries. In summary, our results suggest that (1) a large, upper crustal pluton exists ~7 km beneath the Santa Rita and Hanover-Fierro deposits; (2) abundant phenocryst crystallization occurred at this depth level; and (3) this pluton was the main source for the exsolution of ore-forming fluids. However, the investigated rocks have elevated whole-rock Sr/Y ratios, indicating magma fractionation at deep crustal levels. As a result, our preferred model is a combination of the two end-member models introduced above, with most magma fractionation having occurred in the deep crust and with residual, intermediate to felsic melts having ascended and accumulated at 11 to 14 km paleodepth, where they continued to crystallize with comparatively little crystal-liquid separation, before some of these magmas ascended further to shallow levels and quenched to porphyries.

New Mexico↗
USGS2025

Comparisons of shoreline positions from satellite-derived and traditional field- and remote-sensing techniques

Satellite-derived shorelines (SDS) have the potential to help researchers answer critical coastal science questions and support work to predict coastal change by filling in the spatial and temporal gaps present in current field-based and remote-sensing data collection methods. The U.S. Geological Survey conducted comparison analyses of traditionally sourced shorelines and SDS in diverse coastal landscapes to determine how SDS could be used in ongoing and future work across varied coastal environments and provided some initial findings that could be used for implementation. Using CoastSeg, a browser-based program for SDS detection and mapping, SDS for the period 1984–2023 for multiple locations across the United States were compared to shoreline positions from traditionally sourced shoreline data. In this report, the authors present these comparisons alongside lessons learned and challenges encountered when building SDS workflows in different coastal locations. Results show that individual SDS have larger uncertainty and yet produced similar linear trends to sparser, traditionally sourced shoreline data; because SDS methods provide orders of magnitude more data than traditional shoreline-detection methods, they can be used to evaluate shoreline behaviors. Refining average scalar slopes used in tidal corrections did not result in substantial decreases in uncertainty. Using lessons from this work to outline needs for regional implementation, initial setup time would be considerable, being on the order of weeks. However, once complete, shoreline detections and analyses are fast (on the order of minutes to hours) and achievable using a desktop computer.

Alaska, Florida, Massachusetts, Washington↗
USGS2025

Detection of viral, bacterial, and protozoan pathogens and microbial source tracking markers in paired large- and small-volume water samples

When sampling for waterborne microbes, researchers may need to diverge from recommended sample volumes due to logistical constraints, novel targets, or challenging matrices, with little guidance about the potential impact on results. In field studies, we measured bacteria, viruses, and protozoa (15 quantitative polymerase chain reaction assays) in paired large- and small-volume samples to evaluate method performance and relevant factors. Concordance between methods was low. Large-volume ultrafiltration yielded more detections than small-volume sampling, especially for pathogens in groundwater. Greater microbial concentrations were associated with more frequent detections in small-volume samples and greater concordance between paired samples. Large-volume samples appeared to be more susceptible to diminished sensitivity from complex sample matrices. In laboratory studies, recovery of microbes was poorer for large- than small-volume methods, although large-volume methods more reliably detected low-concentration targets. Large-volume samples were less stable than small-volume samples during storage. Overall, large-volume sampling was superior for detecting pathogens but may underestimate concentrations; small-volume sampling was more prone to false negatives but was adequate when concentrations were relatively high, like we observed for microbial source tracking in surface waters.

Wisconsin↗
USGS2025

An evaluation of the effects of different deicing salt application rates on three watersheds in Essex County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Transportation, evaluated the effects of different deicing salt application rates on surface water, groundwater, and highway runoff quality near State highways in northern New York. Three reaches of State highways were tested with different deicing treatments between October 2019 and November 2022: a salt-sand mixture (Treatment A), a salt mixture applied at a lower rate (Treatment B), and a control mixture consistent with typical deicing salt amounts and application rates. Data on pavement conditions and the quality of surface water, highway runoff, and groundwater were collected. Surface electromagnetic data were also collected. Surface-water and groundwater quality downgradient from the State highways were compared with water quality at upgradient locations. The percentage of snow or ice coverage was used to evaluate the effectiveness of the salt applications. This report provides an overview of the transport of deicing salt. The Treatment B watershed had deicing mixture applied more frequently than other highway reaches, which caused it to have the highest annual total chloride application. Despite differences in chloride application, flow-weighted mean chloride concentrations in highway runoff were comparable across treatments. Chloride concentrations were elevated in surface water and groundwater downgradient from highways relative to chloride concentrations upgradient from highways. A chloride mass balance, calculated for one treatment watershed, indicated that groundwater affected by legacy deicing practices may be contributing additional chloride to surface water. Spatial patterns from electromagnetic surveys show a shallow saline plume alongside the highway in that area. Differences in winter severity and pavement-surface conditions drove deicing salt applications in the treatment areas. This study found that several factors affect chloride loads in the watersheds, including variable winter conditions, adaptive snow and ice management, legacy management practices, and area-specific aquifer and groundwater conditions.

New York↗
USGS2025

Groundwater structures fish growth and production across a riverscape

Landscapes are composed of habitat patches and conditions that vary across space and time. While habitat variability and complexity can support important ecological processes and ecosystem services, the dynamic nature of habitats can also constrain organismal growth and production as optimal conditions are fleeting. In riverine ecosystems, groundwater discharge to streams stabilises water temperature and flow regimes, thus mediating how habitat complexity is expressed. Yet, how stable habitats structure growth and production within the broader landscape matrix is not well understood. In this study, we explored the effects of groundwater on spatiotemporal variation in growth and production for juvenile Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) across the upper Snake River catchment, Wyoming, USA. We combined machine learning techniques and remotely sensed landscape data to estimate groundwater availability across the river network, which we linked to stream temperature regimes and conspecific density. We then used Bayesian hierarchical models to quantify the effects of temperature, density and groundwater on spatiotemporal variation in fish growth and production in 52 focal reaches. Finally, we predicted body size trajectories and trends in total production continuously over both space and time to understand the effect of groundwater at the riverscape scale. Groundwater discharged to streams where topography changes abruptly in valley-bottom areas underlain by coarse glacial deposits. Groundwater stabilised temperature regimes and was associated with high trout densities. Temperature and density, in turn, interacted to influence growth rates: growth increased strongly with temperature, but this effect was reduced when density was high. Accordingly, variation in groundwater availability among stream reaches diversified growth and production regimes. In reaches with low groundwater availability, growth and production declined over time from summer maxima. In contrast, in reaches with high groundwater availability, temporal trends in growth and production were hump-shaped—peaking in autumn—and mean production was greater. At the riverscape scale, temporal asynchrony in growth rates generated convergent spatial variation in growth capacity, but—when combined with density—led to the formation of distinct hotspots of production. Our results demonstrate how groundwater, an important driver of aquatic ecosystem heterogeneity, structures trout growth and production across space and time. Importantly, rare, but stable habitats may disproportionately affect ecological processes and serve as key sources of population diversity at larger spatial scales.

Wyoming↗
USGS2025

Long‐period ground motions from dynamic rupture simulations of large earthquakes on the creeping Hayward–Calaveras–Rodgers Creek fault system

he Hayward, Calaveras, and Rodgers Creek faults in the San Francisco Bay region of California have a high probability of producing a large earthquake in the next decades. Although these faults creep, the creep is insufficient to keep up with their relatively rapid slip rates on their deepest sections, so they have been storing tectonic strain since their last large earthquakes, with the Hayward’s and Rodgers Creek’s more than 150 yr ago. We do not know what the next large Hayward–Calaveras–Rodgers Creek earthquakes will look like or how strongly they will shake the San Francisco Bay region. Harris et al. (2021) used the 3D dynamic (spontaneous) rupture method to simulate large earthquakes on these creeping faults. In this article, we examine the resulting simulated long‐period ( T > 1 s) ground shaking from 0 to 50 km distance, for earthquakes nucleating on the Hayward fault and earthquakes nucleating on the Rodgers Creek fault. We compare these simulated long‐period ground motions with the Boore et al. (2014) well‐established empirically based ground‐motion model suitable for the slowest material velocity in our 3D velocity structure. We find that the simulated long‐period ground motions from the creeping‐fault earthquake scenarios produce a reasonable agreement with the empirical expectations if frictional cohesion is included only where it is appropriate.

California↗
USGS2025

The 1912 Ms 7.2 earthquake in the Denali region of central Alaska

The 2002 M w 7.9 earthquake in central Alaska ruptured 340 km on three faults—Susitna Glacier thrust fault, Denali fault, Totschunda fault—crossing both the Richardson Highway and the Alaska Pipeline. Its occurrence prompted renewed interest in historical large earthquakes that possibly originated on the Denali fault. One of these earthquakes was a M s 7.2 event on July 7, 1912, which we revisit with two approaches: (1) probabilistic relocation of the epicenter using globally recorded arrival times, and (2) compilation and reassessment of shaking intensity reports to estimate a macroseismic epicenter. Our preferred instrumental epicenter is west of the Parks Highway and in agreement with the maximum‐reported shaking, which was from the Parker–Browne expedition of Denali. We also relocated a M s 6.4 aftershock, whose epicenter is 11 km from the mainshock. Candidate faults for the 1912 earthquake include the Denali fault, the McLeod Creek thrust fault, and the Kantishna Hills thrust fault. Future analysis of active faults, paleoseismic results, 1912 instrumental data, and 1912 felt reports may help in interpreting the fault and mechanism of the 1912 earthquake.

Alaska↗
USGS2025

Observations of tear-drinking by lepidopterans on moose ( Alces alces americana ) in northeastern North America

Lepidoptera have long been known to feed on the tears of vertebrates as a presumed source of minerals or nutrients. While this unusual behavior has been observed in a variety of species, only a single previous record has been documented outside of the tropics. Here, we present the first documentation of moths visiting the eyes of a bull moose ( Alces americanus americanus ), captured via trail camera in Green Mountain National Forest, Vermont, United States. We discuss the biogeography of this behavior, how it may differ between tropical and temperate climates, and its potential impact on moose health.

Vermont↗
USGS2025

Searching for seismic precursors - The Barry Landslide hazard clean up

The Barry Landslide, located in Barry Arm of Prince William Sound, Alaska, poses a major hazard due to its steep, unstable slopes and the potential for a massive landslide-generated tsunami. With an estimated volume of 500–700 million cubic meters, the Barry Landslide could trigger highly destructive waves. In this study, we focus on seismic signals from the Barry Landslide, which are critical for providing timely tsunami warnings. Since the summer of 2020, the region has been instrumented to monitor the landslide, but the seismic record is complicated by the presence of nearby glaciers and frequent regional earthquakes. Among these signals, we analyze a specific class of short-duration, high-frequency seismic events that exhibit strong seasonal variability, increasing in rate from late summer to mid-winter before ceasing abruptly in late winter or early spring. Our analysis suggests that the source of these signals is likely near or beneath Cascade Glacier, adjacent to the landslide, rather than within the landslide mass itself. We apply detection algorithms to construct a time history for this signal type, which we then compare with environmental factors like precipitation, temperature, and slope displacement data from ground-based radar and remote sensing. Correlations indicate that these seismic events may be driven by seasonal hydrological changes, particularly the freeze-up of subglacial water pathways. While these events are not directly linked to landslide motion, they serve as indirect markers of subsurface hydrological conditions that influence slope stability. Our findings highlight the complex interplay between glaciers, groundwater, and landslide dynamics, emphasizing the need for multi-parameter monitoring to assess evolving geohazards in the region.

Alaska↗
USGS2025

Groundwater response to managed aquifer recharge at the Southeast Houghton Artificial Recharge Project in Tucson, Arizona

Managed aquifer recharge is a widespread practice for storing water in the subsurface as groundwater. At a managed aquifer recharge facility in southern Arizona, groundwater-level and repeat microgravity data were collected to monitor aquifer response. These data were used to inform parameter identification for an unsaturated-zone flow model used to simulate the recharge process. The facility, the Southeast Houghton Artificial Recharge Project (SHARP), consists of 3 surface basins (about 27,600 square meters [6.8 acres] total surface area) where recycled water is distributed in recharge cycles lasting several months, with dry periods in between. During the study period, December 2020–December 2022, Tucson Water (the City of Tucson’s water utility) reported 6.56×10 6 cubic meters of water (5,320 acre-feet) recharged. Monitoring included groundwater-level observations at 3 monitoring wells and repeat microgravity measurements at as many as 22 locations (some stations were destroyed between surveys). Six gravity surveys were carried out using absolute- and relative-gravity meters. Large gravity increases, more than 250 microgals, were observed during the first repeat survey, 3.5 months after the start of recharge, but only in the immediate vicinity of the recharge basins. Data show that water moved downward to the water table, and storage changes in the unsaturated zone away from the facility were likely minimal. Gravity decreased at stations more than 1 kilometer from the facility, consistent with regional groundwater-level changes. Groundwater-level increases in wells adjacent to the recharge basins began 2 months after the second repeat gravity survey, and 5.5 months after recharge began. Unsaturated-zone flow modeling was carried out using software that simulates water movement and parameter estimation. Model calibration was carried out by minimizing an objective function calculated from the differences between simulated and observed groundwater levels, and between simulated and observed repeat microgravity data. Including repeat microgravity data in the objective function reduced the uncertainty in estimated parameter values for saturated hydraulic conductivity and saturated water content. Modeling indicated that the unsaturated zone between the recharge basins and the water table does not become saturated even after 685 days of simulated infiltration. This gradual wetting may account for increasing infiltration rates over time, as hydraulic conductivity increases with increasing water content. Unsaturated-zone water content decreased rapidly between recharge cycles. Model-simulated groundwater mounding extended about 1 kilometer from the center of SHARP after the 685-day period following the onset of recharge.

Arizona↗
USGS2025

An inset groundwater-flow model to evaluate the effects of layering configuration on model calibration and assess managed aquifer recharge near Shellmound, Mississippi

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model. The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Mississippi↗
USGS2025

The 3D Elevation Program—Supporting the economy of Massachusetts

Introduction Massachusetts extends from the mountains of the Appalachian system in the west of the State to the sandy beaches and rocky shorelines of the Atlantic coast in the east. Inland topographic data support a wide range of important activities, including geologic mapping, transportation planning, forest and wildlife management, quantifying ecological services, water supply protection, commonwealth-wide infrastructure planning, local site planning, and flood-plain management. Nearshore bathymetry can be used to support coastal portions of the Commonwealth by addressing the combined threats of ocean warming, strong storm surge, and rising sea levels. The maintenance and (or) expansion of Massachusetts ports (for instance, Boston, New Bedford) and Cape Cod sediment management depends upon the accurate mapping of bathymetry and the frequent influx of sediment and redeposition. Critical applications that address the broad range of requirements depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Massachusetts. The status of available and in-progress 3DEP baseline lidar data in Massachusetts is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $1.23 million in new benefits annually to Massachusetts. The top 10 Massachusetts business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Massachusetts↗