USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ore Geology Reviews”

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.

830 records · Page 45Linked to original sources
USGS2025

Structural setting and geothermal potential of northeastern Reese River Valley, north-central Nevada: Highly prospective detailed study site for the INGENIOUS project

The northeastern part of the Reese River basin situated ~15 km southeast of Battle Mountain, Nevada, scored highly in the Nevada geothermal play fairway analysis (PFA) for hosting potential hidden geothermal systems. This site (also referred to as Argenta Rise) was therefore chosen for detailed study in the INGENIOUS project (INnovative Geothermal Exploration through Novel Investigations Of Undiscovered Systems). The high PFA scores resulted primarily from favorable structural settings (e.g., fault intersections and pull aparts) with relatively high slip rates on Quaternary faults. The INGENIOUS project is utilizing additional parameters and more rigorous analytical techniques to further advance exploration at this site. This includes integration of geological (e.g., Quaternary fault mapping) and new geophysical datasets (e.g., gravity, magnetics, MT data, and five reprocessed seismic reflection profiles) to build a structural model and to identify specific favorable sites for potential geothermal upwellings. Two-meter temperature surveys were also conducted in the area (139 measurements). This part of north-central Nevada is characterized by systems of intersecting northerly and ENE-striking faults within the broader Humboldt structural zone, a poorly understood belt of ENE-striking faults and relatively high heat flow extending across northern Nevada. Kinematic analysis of exposed fault surfaces shows that ENE-striking faults have accommodated sinistral-normal slip, and normal slip characterizes N- to NNE-striking faults. Northeastern Reese River Valley lies within a broad left step between major ENE-striking fault zones on the northern flanks of the Argenta Rim and Shoshone Range and thus corresponds to a broad pull-apart in the ENE-striking sinistral-normal fault system. Notably, the nearby Beowawe geothermal system in Whirlwind Valley (with abundant sinter, hot springs, and a geothermal power plant) occupies a fault intersection in a relatively small left step in a major ENE-striking sinistral-normal fault and may serve as an analogue for a potential hidden system in northeastern Reese River Valley. Existing geological maps, high-resolution lidar, and seismic reflection data demonstrate that northeastern Reese River Valley is structurally complex with multiple intersections between the ENE- and N- to NNE-striking fault systems. Some of these fault intersections correspond to low resistivity anomalies, magnetic lows, and/or very subtle 2-m temperature anomalies, which may indicate hidden geothermal upwellings. Three-dimensional modeling and temperature-gradient drilling are planned to further evaluate these sites for geothermal activity.

Nevada↗
USGS2025

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences↗
USGS2025

ShakeAlert® version 3: Expected performance in large earthquakes

The ShakeAlert earthquake early warning (EEW) system partners along with U.S. Geological Survey (USGS) licensed operators deliver EEW alerts to the public and trigger automated systems when a significant earthquake is expected to impact California, Oregon, or Washington. ShakeAlert’s primary goal is to provide usable warning times before the arrival of damaging shaking. EEW is most likely to achieve this goal in large‐magnitude earthquakes. In recent years, ShakeAlert has gone through a series of upgrades to its underlying scientific algorithms aimed at improved performance during large earthquakes. Version 3 of this software recently went live in the production system and includes improvements to all algorithms. The main seismic algorithms that detect an earthquake and characterize its location, magnitude, and fault rupture orientation are faster than older versions. Other key changes include: using real‐time geodetic data to characterize the magnitude growth in large earthquakes; the introduction of an alert pause procedure to compromise between speed near the epicenter and improved accuracy at larger distances; and the inclusion of a nonergodic site‐response model in the ground‐motion predictions. ShakeAlert has achieved its primary goal of usable warning times before strong shaking at some locations in real‐time operations in recent M 6 earthquakes. Using offline tests, we demonstrate usable warning times are possible for many sites with peak shaking values of modified Mercalli intensity (MMI) 7–8 in M 7+ earthquakes and also for many MMI 8–9 sites in M 8+ earthquakes. ShakeAlert partners use a variety of MMI and magnitude thresholds in deciding when to alert their users within bounds set by the USGS. Our study shows that there is room to raise the magnitude thresholds up to about M 5.5 without adversely affecting performance in large earthquakes. The ground‐motion criteria are more complex owing to a significant drop‐off in warning times between the MMI 4 and 5 levels of predicted shaking. However, widely used ShakeAlert products, such as the MMI 3 and 4 contour products, can provide sufficiently long warning times before strong shaking in moderate‐to‐great earthquakes to enable a range of protective actions.

Bulletin of the Seismological Society of America↗
USGS2025

Base-flow sampling to enhance understanding of the groundwater flow component of nitrogen loading in small watersheds draining into Long Island Sound

Excessive nitrogen discharge is a major concern for the Long Island Sound. Programs have been implemented to reduce point sources of nitrogen to the sound, but little is known about the nonpoint sources. This study aims to better understand the current groundwater contributions of nitrogen from nonpoint sources in the Long Island Sound watershed. During the spring and summer of 2022, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, collected water-quality samples to analyze nutrients (nitrogen and phosphorus), chloride, and bromide at 45 stations in the Long Island Sound watershed in Connecticut, New York, and Rhode Island. The stations were in small drainage watersheds (5 to 30 square kilometers) in the southern part of the Long Island Sound watershed. During two separate synoptic sampling events, water-quality samples and instantaneous streamflow measurements were collected under base-flow conditions (where the streamflow is dominated by groundwater inputs rather than overland flow or runoff flow). One sampling event was in the nongrowing season (April 24–25, 2022), and the other was in the growing season (June 30–July 1, 2022). To calculate instantaneous nitrogen loads and yields, streamflow was measured at the time of sample collection. Nitrogen concentrations, loads, and yields varied among sampling stations and by season. Total filtered nitrogen concentrations were generally lower in the nongrowing season (from less than 0.14 to 1.9 milligrams per liter) than in the growing season (from less than 0.23 to 3.0 milligrams per liter). Nitrate plus nitrite concentrations showed little variation between the nongrowing and growing seasons. Unfiltered ammonia plus organic nitrogen concentrations were generally lower in the nongrowing season (from less than 0.07 to 0.83 milligram per liter) than in the growing season (from 0.11 to 0.98 milligram per liter). In contrast, total filtered and unfiltered nitrogen loads and yields were higher in the nongrowing season than during the growing season, likely because streamflows were higher during the nongrowing season. Total unfiltered nitrogen yields during the nongrowing season ranged from less than 0.15 to 5.0 kilograms per square kilometer per day. Total unfiltered nitrogen yields during the growing season ranged from less than 0.12 to 2.5 kilograms per square kilometer per day. Total filtered nitrogen yields during the nongrowing season ranged from less than 0.13 to 5.2 kilograms per square kilometer per day. Total filtered nitrogen yields during the growing season ranged from less than 0.06 to 2.5 kilograms per square kilometer per day.

Connecticut, New York, Rhode Island↗
USGS1999

Water quality of Rob Roy Reservoir and Lake Owen, Albany County, and Granite Springs and Crystal Lake Reservoirs, Laramie County, Wyoming, 1997-98

The water quality of four reservoirs was assessed during 1997 and 1998 as a cooperative project between the Cheyenne Board of Public Utilities and the U. S. Geological Survey. The four reservoirs, Rob Roy, Lake Owen, Granite Springs, and Crystal Lake, provide approximately 75 percent of the public water supply for Cheyenne, Wyoming. Samples of water and bottom sediment were collected and analyzed for selected physical, chemical, and biological characteristics to provide data about the reservoirs. Water flows between the reservoirs through a series of pipelines and stream channels. The reservoirs differ in physical characteristics such as elevation, volume, and depth. Profiles of temperature, dissolved oxygen, specific conductance, and pH were examined. Three of the four reservoirs exhibited stratification during the summer. The profiles indicate that stratification develops in all reservoirs except Lake Owen. Stratification developed in Rob Roy, Granite Springs, and Crystal Lake Reservoirs by mid-July in 1998 and continued until September, with the thickness of the epilimnion increasing during that time. Secchi disk readings indicated Rob Roy Reservoir had the clearest water of the four reservoirs studied. The composition of the phytoplankton community was different in the upper two reservoirs from that in the lower two reservoirs. Many of the species found in Rob Roy Reservoir and Lake Owen are associated with oligotrophic, nutrient-poor conditions. In contrast, many of the species found in Granite Springs and Crystal Lake Reservoirs are associated with mesotrophic or eutrophic conditions. The total number of taxa identified also increased downstream. The chemical water type in the reservoirs was similar, but dissolved-solids concentrations were greater in the downstream reservoirs. Water in all four reservoirs was a calcium-bicarbonate type. In the fall of 1997, Rob Roy Reservoir had the lowest dissolved-solids concentration (19 milligrams per liter), whereas Crystal Lake Reservoir had the highest concentration (63 milligrams per liter). Relatively little differences in the concentrations of major-ion species were noted between samples collected near the surface and near the bottom of the same reservoir. In contrast, iron and manganese concentrations generally were higher in samples collected near the bottom of a reservoir than in near-surface samples collected from the same reservoir. Composite bottom-sediment samples from all four reservoirs contained similar concentrations of bulk constituents such as aluminum, iron, phosphorus and titanium, but varied in concentrations of trace elements. Trace-element concentrations in Rob Roy Reservoir and Lake Owen were similar to the crustal average, whereas in Granite Springs and Crystal Lake Reservoirs the concentrations were similar to granitic rocks.

Wyoming↗
USGS2025

Gravity and magnetic surveys of the Skaergaard intrusion, East Greenland

Aeromagnetic and gravity surveys of the Skaergaard intrusion in East Greenland were carried out in July–August 1971 as part of a grant to the University of Oregon Center for Volcanology to refine the models of crystallization and differentiation of the intrusion, specifically to test whether the intrusion is underlain by dense rocks of a reservoir 20 kilometers (km) thick (referred to as a “hidden zone”). The Skaergaard intrusion is a source of platinum group elements that are critical mineral resources for many technologies, and because no new data have been collected these legacy datasets remain a valuable asset. The total-intensity aeromagnetic survey was flown in early July 1971 with a proton precession magnetometer at a constant barometric altitude of 1.5 km (5,000 feet) with a nominal line spacing of 1 km. Two gravimeters were used to acquire 168 stations of which 86 were at known altitudes (mainly sea level) and 82 had altitudes measured by altimetry in late July–August 1971. Finally, a north-south ground vertical-intensity magnetic traverse was completed across the intrusion together with collection of oriented hand specimens. The hand specimens were measured for remnant magnetization and density, along with density measurements of more specimens collected by expedition geologists for other purposes. The intrusion is composed of layered gabbro with extensive crystal fractionation that is dense and strongly reversely polarized. After terrain correction and standard Bouguer gravity reduction, the gravity anomaly dataset was corrected for all rock above sea level using the density measurements of the various zones of the intrusion and the topographic and geologic maps (variable density Bouguer gravity reduction). A large regional gradient in the gravity anomaly data was removed using orthogonal polynomial fitting to the gridded data. The zonal volumes of rock below sea level were calculated from the dipping polygonal layer gravity model of the intrusion below sea level and combined with elliptic cross–section cylinders for the various zones above sea level to approximate the original zonal volumes of the intrusion. The residual gravity anomaly of 18–20 milligals (mGal) was only about half of the expected anomaly if a large hidden zone proposed from petrologic considerations were present, and both two-dimensional and three-dimensional models imply that the exposed series of intrusion zones explain the gravity anomaly by their down-dip extension below sea level together with a small hidden-zone volume. A three-dimensional model of the exposed rocks and their down-dip extension below sea level also can account for the aeromagnetic anomaly with little or no requirement for hidden-zone rock. The middle and upper zone units of the intrusion contain the most magnetite and account for most of the aeromagnetic anomaly.

Skaergaard intrusion↗
USGS2025

Performance analysis of oil recovery and CO2 retention in a greenfield residual oil zone: CO2-EOR in Tall Cotton Field (Permian Basin, West Texas, USA)

Residual oil zones (ROZs) can offer significant oil resources via enhanced oil recovery (EOR) as well as subsurface carbon dioxide (CO 2 ) retention during injection. If injected CO 2 is anthropogenic, the ROZs can offer a substantial geologic storage potential. The ROZs below the oil/water contact (OWC) of main pay zones (MPZ) in conventional reservoirs or brownfields, are more commonly developed for CO 2 injection and oil production and reported in the literature. However, CO 2 -EOR in greenfield ROZs, reservoirs without a MPZ present, have rarely been developed for CO 2 -EOR operation. The Tall Cotton Field of West Texas, Permian Basin, which started production in 2015 (Phase 1) and expanded in 2017 (Phase 2) from the San Andres Limestone, is one of the first examples of greenfield ROZs developed for EOR by injecting CO 2 . This paper analyses EOR and CO 2 retention performance of Tall Cotton Field using allocated injection and production data from inverted 5-spot well patterns of Phase-1 and -2 developments. Production and injection data allocated to each of the 28 identified patterns (nine 20-acre patterns for Phase-1, three 20-acre and sixteen 10-acre patterns for Phase-2) were analyzed for historical and forecasted oil recovery using ratio-trend decline analysis, and for CO 2 retention performance of the patterns. The allocated data were further used to calculate injected reservoir pore volume and void replacement ratios (VRR) for the analysis period. Quantitative results indicated that oil recovery factors of the 5-spot patterns varied between 4–10 %, and 5–30 % between the end of injection and the forecast periods, respectively. Storage of CO 2 , on the other hand, increased to a mean value of ∼7130 MMscf per pattern in Phase-1 and to a mean storage of 3700 MMscf per pattern in Phase-2 until the end of injection, followed by a decline after the end of injection and into the forecast period. Resulting CO 2 utilization factors ∼6–50 Mscf/bbl were estimated at the end of injection. Overall, presented results suggested that developing greenfield ROZs for CO 2 -EOR can be as promising as brownfield ROZs and mature MPZs for EOR and underground storage of injected CO 2 . For Tall Cotton Field, results suggest that Phase-2 patterns generally outperformed Phase-1 for oil recovery factors, while Phase-1 performed better in CO 2 retention performance metrics. This is the first study in the literature that reports a detailed CO 2 -EOR performance analysis of a greenfield ROZ in the Permian Basin, which can potentially allow for comparison with MPZs and brownfield ROZs.

Texas↗
USGS2025

Estimated hydrogeologic, spatial, and temporal distribution of self-supplied domestic groundwater withdrawals for aquifers of the Virginia Coastal Plain

Water use from private-domestic wells accounts for nearly 40 percent of total groundwater withdrawals in the Virginia Coastal Plain Physiographic Province (henceforth called the Virginia Coastal Plain). However, because self-supplied domestic water use generally falls below the Virginia Department of Environmental Quality (VDEQ) reporting and management threshold of 300,000 gallons per month, quantifying these withdrawals is challenging. This report builds upon the foundation of previous U.S. Geological Survey investigations by providing revised techniques to improve estimates of the aquifer source, spatial distribution, and monthly magnitude of these groundwater withdrawals. The aquifer sources of private-domestic wells in the Virginia Coastal Plain were estimated by cross-referencing 8,264 well records from the VDEQ and the Virginia Department of Health to a digital model of the Virginia Coastal Plain hydrogeologic framework. This analysis highlights the regional importance of the Yorktown-Eastover, Potomac, and surficial aquifers. Collectively, these three aquifers account for 80 percent of self-supplied domestic groundwater withdrawals. The population using self-supplied domestic water was estimated using census blocks, well-use ratios, building footprints, and land-use and land-cover data to produce a high-resolution, disaggregated, raster-based dataset. This approach improves upon previous models at the census-block or road-network scale by reducing the low-density spread of the self-supplied domestic population across undeveloped areas and concentrating the population and its corresponding water use in the areas where it is most likely to occur. Results show that an estimated 475,332 people comprise the 2020 self-supplied domestic population of the Virginia Coastal Plain, an increase of 5.7 percent since 2010, and the greatest concentrations of self-supplied domestic population surround large cities. Estimates could be further refined with the addition of current and complete spatial data on public water-system service areas. The quantity of water used by the self-supplied domestic population was estimated by modifying published state per-capita water-use coefficients with the corresponding monthly variability assessed from Virginia Coastal Plain public water-system withdrawal data. This analysis estimates an average increase of 12 percent from June through August and an average decrease of 8 percent from December through March from the baseline annual average of 80 gallons per day per capita, which generally matches similar studies in the eastern United States. The application of these revised methodologies for the estimation of private-domestic wells and the self-supplied domestic population improves understanding of domestic groundwater use in the Virginia Coastal Plain across hydrogeologic, spatial, and temporal scales. These revisions help better inform water-resource managers and decision makers and support higher resolution groundwater modeling. Furthermore, these methods are transferrable to other areas where self-supplied domestic water withdrawals are important to the overall water budget.

Virginia↗
USGS1999

Simulation of the effects of development of the ground-water flow system of Long Island, New York

Extensive development on Long Island since the late 19th century and projections of increased urbanization and ground-water use makes effective water-resource management essential for preservation of the island's hydrologic environment and maintenance of a reliable source of water supply. This report presents results of a ground-water flow simulation analysis of the effects of development on the Long Island ground-water system. It describes ground-water levels, stream-flow, and the ground-water budget for the predevelopment period (pre-1900), the 1960's drought, and a more recent (1968-83) period with significant hydrologic stress. The report also presents estimated effects of a proposed water-supply strategy for the year 2020. Long Island has three major aquifers-the upper glacial (water-table), the Magothy, and the Lloyd aquifers-that are separated to varying degrees by confining units. Before development, recharge from precipitation entered the ground-water system at a rate of more than 1.1 billion gallons per day. An equal amount discharged to streams (41 percent), the shore (52 percent), and subsea boundaries (7 percent) . Urbanization and withdrawal of more than 400 Mgal/d (million gallons per day) from wells have resulted in local effects that include declines in ground-water levels, drying up and burial of streams and wetlands, reduction of ground-water recharge by increased overland flow to the ocean, a general decrease in ground-water discharge, and salt water intrusion. In some areas, the reduction in recharge is mitigated by leakage from water-supply and wastewater disposal lines, and infiltration of storm water through recharge basins. During 1968-83, a net loss of 240 Mgal/d from the ground-water system caused a decrease in ground-water discharge to streams (135 Mgal/d), to the shore (82 Mgal/d), and to subsea boundaries (23Mgal/d).The greatest adverse effects have been in western Long Island, where the most severe development has occurred. This analysis shows stream base flow to be highly sensitive to water-table fluctuations, and long streams to be more sensitive than short ones. A water-supply scenario for the year 2020 was simulated that employs redistribution of pumping centers to mitigate extreme local effects . Although the net stress on the ground-water system was projected to increase 57 Mgal/d (24 percent) above that of 1968-83, redistribution of ground-water withdrawals across the island would allow recovery of cones of depression in western Long Island, thereby reducing the threat of salt water intrusion and increasing base flow of some streams . The increased stress would cause a net decrease in base flow island wide of 44 Mgal/d; total base flow would be 281 Mgal/d - 39 percent below predevelopment levels or 14 percent below 1968-83 levels. The most severe effects would be in Nassau and western Suffolk Counties.

New York↗
USGS2026

Determining Volcanic Risk in Auckland (DEVORA) Research Programme—A transdisciplinary approach to address the challenge of distributed volcanism in an urban environment

The Determining Volcanic Risk in Auckland (DEVORA) Research Programme was launched in 2008 to address the challenges associated with monogenetic volcanism in an urban setting and to enhance volcanic risk management in Tāmaki Makaurau Auckland in Aotearoa New Zealand. It is a multi-agency, increasingly transdisciplinary (defined here as research that transcends traditional disciplinary boundaries by integrating diverse types of knowledge, perspectives, and methods from academic and non-academic participants to create novel solutions to complex problems), and collaborative research program jointly led by Waipapa Taumata Rau University of Auckland and Earth Sciences New Zealand (ESNZ; formerly GNS Science), with core funding from Natural Hazards Commission Toka Tū Ake (NHC; formerly the Earthquake Commission, EQC) and Te Kaunihera o Tāmaki Makaurau Auckland Council (AC). The primary research focus of DEVORA is to investigate the geologic history, volcanic hazards, and risk posed by the basaltic intraplate Auckland Volcanic Field. Disruption from ash fall and gas from other Aotearoa New Zealand volcanoes is also considered. DEVORA’s work to explore exposure and vulnerability in Tāmaki Makaurau Auckland is also useful for assessing risks from other non-volcanic natural hazards, such as seismic and tsunami hazards. The greater Tāmaki Makaurau Auckland region has an ethnically and socio-economically diverse population of approximately 1.7 million, representing about one-third of the Aotearoa New Zealand population, and hosts critical infrastructure of national significance. The size and nature of the populace, consequential economic base, and important infrastructure within Tāmaki Makaurau Auckland mean that the effects of a volcanic eruption would be felt nationally, including through the disruption of air travel to Aotearoa New Zealand. The hazards from such an eruption could potentially affect hundreds of thousands of people, businesses, and lifelines (critical infrastructure). A considerable challenge for emergency and risk managers is the monogenetic nature of the volcanic field. It is not known where or when the next eruption will occur, how much warning we may get before an eruption, nor how an eruption and its effects might unfold. In this contribution, we highlight the concept and collaborative intent of the DEVORA Programme and show how it has evolved over the 16 years since its inception. We describe how DEVORA has unified more than 100 researchers (including more than 50 graduate students) and numerous stakeholders to address key issues facing Tāmaki Makaurau Auckland and describe how research findings are being implemented into policy and communicated to stakeholder agencies and the public. We also illustrate the broader influence of the DEVORA Programme and provide some learnings that might benefit others embarking on similar integrated projects, especially those focused on distributed volcanism in and near populated areas.

Auckland↗
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↗
USGS2026

Precipitation-based flood-inundation maps for the East Fork Little Blue River and tributaries at Lee’s Summit, Missouri, 2024

The U.S. Geological Survey, in cooperation with the City of Lee’s Summit, Missouri, assessed flooding of the East Fork Little Blue River and tributaries for varying precipitation magnitudes and durations, varying antecedent runoff conditions, and projected climate-change conditions. The precipitation scenarios were used to develop a library of flood-inundation maps for a 2.95-mile reach of the East Fork Little Blue River and tributaries within the city. A two-dimensional U.S. Army Corps of Engineers Hydrologic Engineering Center–River Analysis System (HEC–RAS; ver. 6.5) rain-on-grid model was calibrated to selected runoff events representing a range of antecedent runoff conditions and hydrologic responses. Lowest adjacent grades for structures within the nearby study area were incorporated into the terrain, and depth grids and water-surface elevation grids were developed for the study area. Simulated velocities at selected bridge locations were also developed from the model. The model was calibrated using water-surface elevation data collected from water-level loggers (pressure transducers) and streamflow measurements and water-surface elevation measurements made at a reference point during runoff events. The calibrated HEC–RAS model was used to simulate streamflows from design rainfall events of 15-minute to 24-hour durations and ranging from a 100- to 0.1-percent annual exceedance probability (1-year to 1,000-year recurrence intervals). Flood-inundation maps were produced for depths at a reference location of 3 to 16 feet, or a depth exceeding the 0.1-percent annual exceedance probability interval precipitation. The results of each precipitation duration-frequency value were represented by a 1-foot-increment inundation map based on the generated peak streamflow from that rainfall event and the corresponding water-surface elevation at the East Fork Little Blue River reference location. Within the HEC–RAS model, 240 scenarios were developed from the design rainfall events with each of 3 antecedent conditions. Additional scenarios were created to simulate the effects of projected precipitation scenarios on the 100-year recurrence interval, 24-hour storm and the 100-year recurrence interval, 6-hour storm. All simulation results were assigned to a flood-inundation map condition based on the generated peak flow and corresponding water-surface elevation at the East Fork Little Blue River reference location. The flood-inundation maps are shown on a web mapping application made available to the public through the City of Lee’s Summit (hyperlink will be added when available). The flood-inundation maps are tied to real-time precipitation data obtained from the Automated Surface Observing System weather station at the Lee’s Summit Municipal Airport, accessible at https://mesonet.agron.iastate.edu/request/download.phtml?network=MO_ASOS . The availability of these maps, along with information regarding observed rainfall, could help provide emergency management personnel and residents with information that is critical for flood-response activities, such as evacuations and road closures, and for postflood recovery efforts.

Missouri↗
USGS2024

Three-dimensional temperature maps of the Williston Basin, USA: Implications for deep hot sedimentary and enhanced geothermal resources

As part of U.S. Geological Survey's (USGS) efforts to identify and assess geothermal energy resources of the US, a three-dimensional (3D) geologic and thermal model has been constructed for the Williston Basin, USA. The geologic model consists of all sedimentary units above the Proterozoic and Archean crystalline rock (called basement herein), with a total sedimentary thickness of up to 5 km near the basin center. Twenty-nine geologic units were mapped from interpreted formation tops from 16,465 wells. A 3D temperature model was constructed to a depth of 7 km by constructing a 3D heat flow model for the sedimentary units, followed by estimating underlying temperature using a one-dimensional (1D) analytic solution for heat flow within the underlying crystalline basement. Using the sedimentary basin model, heat flow was simulated in 3D and was calibrated using three temperature datasets: 1) 24 high-confidence static temperature logs (equilibrium thermal profiles), 2) more than15,000 drill stem test (DST) measurements from >7,000 wells, and 3) more than 45,000 bottomhole temperature (BHT) measurements from >14,000 wells. The DST and BHT datasets provide broad spatial coverage, but are lower confidence, primarily because measurements were made prior to attaining thermal equilibrium. DST and BHT measurements were binned regionally to develop representative thermal profiles that generally agree with these lower quality data (hereafter called pseudowell temperature profiles). Layer properties (primarily thermal conductivity and compaction curves) were set to best estimate values, then the heat flow model was calibrated to fit pseudowell and static temperature logs primarily by adjusting basal heat flow to approximate the overall temperature profile. Minor adjustments to thermal conductivity allowed adjusting changes in slope at lithologic contacts. Resulting maps include 3D temperature and basal (bottom of sedimentary units) heat flow estimates, which are used as input for the temperature model of the basement. The crystalline basement temperature model uses an analytic 1D solution to the heat flow equation that requires estimates of heat flow and temperature at the upper boundary (i.e., the sediment/basement contact), radiogenic heat production within the crystalline basement, and reference thermal conductivity (i.e., uncorrected for temperature). Two regions of high heat flow are identified: 1) in western North Dakota along the North American Central Plains Conductivity Anomaly and 2) in eastern Montana near the Poplar dome. Within the sedimentary column in the center of the basin of the basin, an area of approximately 100,000 km2 is predicted to have moderate- to high-temperature geothermal resources (>90 °C) under the thickest sequences of sediments. Where thick insulation and high heat flow coincide, electric-grade resources can be less than 4 km deep. Assuming a maximum feasible drilling depth of 7 km, temperatures are predicted to be as high as 175 °C. The geologic model may be used to identify strata at sufficient temperatures that may have natural permeability or that may have conditions that favor development of enhanced/engineered geothermal systems resources.

Montana, North Dakota, South Dakota↗
USGS2024

Quality of groundwater used for domestic supply in the eastern Sacramento Valley and adjacent foothills, California

Summary More than 2 million Californians rely on groundwater from privately owned domestic wells for drinking-water supply. This report summarizes a water-quality survey of domestic and small-system drinking-water supply wells in the eastern Sacramento Valley and adjacent foothills where more than 25,000 residents are estimated to use privately owned domestic wells. Study results show that inorganic and organic constituents in groundwater were present above regulatory (maximum contaminant level, MCL) benchmarks for public drinking-water quality in 8 and 3 percent, respectively, of the aquifer area used for domestic drinking-water supply (herein, “domestic groundwater resources”; fig. 1). The only inorganic constituent detected above regulatory benchmarks was arsenic. The only organic constituent exceeding regulatory benchmarks was the fumigant 1,2,3-trichloropropane (1,2,3-TCP). Three additional organic constituents—the disinfection by-product chloroform, the gasoline oxygenate methyl tert -butyl ether (MTBE), and the solvent tetrachloroethene (PCE)—were detected at low concentrations below one-tenth of regulatory benchmarks in 34, 10, and 10 percent of domestic groundwater resources, respectively. Total dissolved solids (TDS), iron, and manganese exceeded non-regulatory aesthetic guidelines for drinking water in 5, 10, and 26 percent of domestic groundwater resources, respectively. Per- and polyfluoroalkyl substances (PFASs) were detected in 29 percent of domestic groundwater resources,with 5 percent exceeding the recently enacted (April 2024) U.S. Environmental Protection Agency MCLs. Total coliform and enterococci bacteria were detected in 13 and 8 percent of domestic groundwater resources, respectively. Redox sensitive constituents in this study included arsenic, manganese, nitrate, and iron. In the lower elevation portions of the eastern Sacramento Valley study area, reducing conditions in groundwater aquifers promote elevated arsenic, iron, and manganese, and conversely lower concentrations of nitrate. The presence of the volatile organic compound (VOC) 1,2,3-TCP was related to its past history in select agricultural land uses (on orchards or vineyards) in the Sacramento Valley; however, unlike in the San Joaquin Valley where orchards and vineyards are more common, its detection frequency was low (only detected in one well in this study). Chloroform was frequently detected in this study at low levels. Chloroform is a disinfection byproduct commonly found in domestic wells treated by shock chlorination. The solvent PCE is among the most frequently detected VOCs in groundwater, which is primarily related to its long history of use and its persistence in groundwater in oxic conditions. The gasoline oxygenate MTBE was a contaminant introduced to groundwater through atmospheric exchange when it was used as a fuel additive to decrease smog inducing emissions from vehicles. Its occurrence in groundwater at low levels is expected and makes it a potentially useful tracer of relatively recent recharge water being withdrawn from wells. The PFASs are anthropogenic chemicals with hundreds of uses, and they have been incorporated into many different products, processes, and applications worldwide. Like MTBE, the occurrence of PFASs in groundwater may be in part due to atmospheric exchange, but there are several other pathways that contribute PFASs to the environment.

California↗
USGS1990

Recharge rates and aquifer hydraulic characteristics for selected drainage basins in middle and east Tennessee

Quantitative information concerning aquifer hydrologic and hydraulic characteristics is needed to manage the development of ground-water resources. These characteristics are poorly defined for the bedrock aquifers in Middle and East Tennessee where demand for water is increasing. This report presents estimates of recharge rate, storage coefficient, diffusivity, and transmissivity for representative drainage basins in Middle and East Tennessee, as determined from analyses of stream-aquifer interactions. The drainage basins have been grouped according to the underlying major aquifer, then statistical descriptions applied to each group, in order to define area1 distribution of these characteristics. Aquifer recharge rates are estimated for representative low, average, and high flow years for 63 drainage basins using hydrograph analysis techniques. Net annual recharge during average flow years for all basins ranges from 4.1 to 16.8 in/yr (inches per year), with a mean value of 7.3 in. In general, recharge rates are highest for basins underlain by the Blue Ridge aquifer (mean value11.7 in/yr) and lowest for basins underlain by the Central Basin aquifer (mean value 5.6 in/yr). Mean recharge values for the Cumberland Plateau, Highland Rim, and Valley and Ridge aquifers are 6.5, 7.4, and 6.6 in/yr, respectively. Gravity drainage characterizes ground-water flow in most surficial bedrock aquifer in Tennessee. Accordingly, a gravity yield analysis, which compares concurrent water-level and streamflow hydrographs, was used to estimate aquifer storage coefficient for nine study basins. The basin estimates range from 0.002 to 0.140; however, most estimates are within a narrow range of values, from 0.01 to 0.025. Accordingly, storage coefficient is estimated to be 0.01 for all aquifers in Middle and East Tennessee, with the exception of the aquifer in the inner part of the Central Basin, for which storage coefficient is estimated to be 0.002. Estimates of aquifer hydraulic diffusivity are derived from estimates of the streamflow recession index and drainage density for 75 drainage basins; values range from 3,300 to 130,000 ft^2/d (feet squared per day). Basin-specific and site-specific estimates of transmissivity are computed from estimates of hydraulic diffusivity and specific-capacity test data, respectively. Basin-specific, or areal, estimates of transmissivity range from 22 to 1,300 ft^2/d, with a mean of 240 ft^2/d In general, areal transmissivity is highest for basins underlain by the Cumberland Plateau aquifer (mean value 480 ft^2/d) and lowest for basins underlain by the Central Basin aquifer (mean value 79 ft^2/d). Mean transmissivity values for the Highland Rim, Valley and Ridge, and Blue Ridge aquifer are 320,140, and 120 ft^2/d respectively. Site-specific estimates of transmissivity, computed from specific-capacity data from 118 test wells in Middle and East Tennessee range from 2 to 93,000 ft^2/d with a mean of 2,600 ft^2/d Mean transmissivity values for the Cumberland Plateau, Highland Rim, Central Basin, Valley and Ridge, and Blue Ridge aquifers are 2,800,1,200, 7,800, 390, and 65Oft Id, respectively.

Tennessee↗
USGS1998

Ground-water flow in the surficial aquifer system and potential movement of contaminants from selected waste-disposal sites at Naval Station Mayport, Florida

Ground-water flow through the surficial aquifer system at Naval Station Mayport near Jacksonville, Florida, was simulated with a two-layer finite-difference model as part of an investigation conducted by the U.S. Geological Survey. The model was calibrated to 229 water-level measurements from 181 wells during three synoptic surveys (July 17, 1995; July 31, 1996; and October 24, 1996). A quantifiable understanding of ground-water flow through the surficial aquifer was needed to evaluate remedial-action alternatives under consideration by the Naval Station Mayport to control the possible movement of contaminants from sites on the station. Multi-well aquifer tests, single-well tests, and slug tests were conducted to estimate the hydraulic properties of the surficial aquifer system, which was divided into three geohydrologic units�an S-zone and an I-zone separated by a marsh-muck confining unit. The recharge rate was estimated to range from 4 to 15 inches per year (95 percent confidence limits), based on a chloride-ratio method. Most of the simulations following model calibration were based on a recharge rate of 8 inches per year to unirrigated pervious areas. The advective displacement of saline pore water during the last 200 years was simulated using a particle-tracking routine, MODPATH, applied to calibrated steady-state and transient models of the Mayport peninsula. The surficial aquifer system at Naval Station Mayport has been modified greatly by natural and anthropogenic forces so that the freshwater flow system is expanding and saltwater is being flushed from the system. A new MODFLOW package (VAR1) was written to simulate the temporal variation of hydraulic properties caused by construction activities at Naval Station Mayport. The transiently simulated saltwater distribution after 200 years of displacement described the chloride distribution in the I-zone (determined from measurements made during 1993 and 1996) better than the steady-state simulation. The advective movement of contaminants from selected sites within the solid waste management units to discharge points was simulated using MODPATH. Most of the particles were discharged to the nearest surface-water feature after traveling less than 1,000 feet in the ground-water system. Most areas within 1,000 feet of a surface-water feature or storm sewer had traveltimes of less than 50 years, based on an effective porosity of 40 percent. Contributing areas, traveltimes, and pathlines were identified for 224 wells at Naval Station Mayport under steady-state and transient conditions by back-tracking a particle from the midpoint of the wetted screen of each well. Traveltimes to contributing areas that ranged between 15 and 50 years, estimated by the steady-state model, differed most from the transient traveltime estimates. Estimates of traveltimes and pathlines based on steady-state model results typically were 10 to 20 years more and about twice as long as corresponding estimates from the transient model. The models differed because the steady-state model simulated 1996 conditions when Naval Station Mayport had more impervious surfaces than at any earlier time. The expansion of the impervious surfaces increased the average distance between contributing areas and observation wells.

Florida↗
USGS2025

Updating and recalibrating the integrated Santa Rosa Plain Hydrologic Model to assess stream depletion and to simulate future climate and management scenarios in Santa Rosa, Sonoma County, California

The Santa Rosa Plain Hydrologic Model (SRPHM) was developed and published in 2014 through a collaboration between the U.S. Geological Survey (USGS) and Sonoma Water to analyze the hydrologic system in the Santa Rosa Plain watershed, help meet the increasing demand for fresh water, and prepare for future uncertainties in water resources. The original model simulated hydrological conditions and water use from water years 1975 to 2010. Recently (2023), the USGS, in cooperation with Sonoma Water and the California State Water Resources Control Board, updated the SRPHM model to extend its simulation period to the end of the 2018 calendar year, incorporate new estimates of rural and agricultural water use, and use efficient input format for climate variables. The updated model was recalibrated, and evaluation of the new model calibration is included in this report. This report presents the results of comparing the hydraulic heads, streamflow, and groundwater budget simulated by the updated model with those generated by the original model and observed data. The main difference in the simulated budget between the original and updated SRPHM is the estimates of agricultural pumping, rural domestic pumping, and return flow generated from rural water use that was not simulated in the original model. The revised agricultural pumping is simulated using the agricultural package, which constrains pumping to available groundwater. The use of the agricultural package leads to a more realistic estimation of agricultural water use, with revised agricultural pumping being one-third less than that in the original model. The revised rural pumping is about half of the pumping in the original model because of using detailed parcel data to estimate population density in rural areas instead of coarse census tracts. Overall, average total inflows for water years 2006–10 simulated by the updated model were about 2 percent less than the original model, and the average total updated outflows were nearly 5 percent less than the original model. The updated model was then used to generate stream depletion maps, simulate climate change scenarios during 2019–99, and simulate water rights allocation using the Model for Decision Support in Integrated River Basin Management (MODSIM). The results from simulating eight future climate scenarios indicated either an increase in groundwater storage or no significant change in the next 80 years, along with an increase in recharge, an increase in actual evapotranspiration in six out of eight climate projections, and an increase in surface runoff. The increases in the simulated future groundwater storage, recharge, evapotranspiration, and runoff in most climate projections are mainly driven by the projected increase in precipitation in most of the future climate scenarios. The updated model also was used to test a pilot case study demonstrating water-resource allocation among different users with different water rights using the integrated MODSIM-Groundwater and Surface-Water Flow Model (GSFLOW) platform. The updated SRPHM serves as a valuable tool for analyzing historical and future hydrologic conditions in the Santa Rosa Plain watershed and preparing for future uncertainties.

California↗
USGS2025

Cancer risk and estimated lithium exposure in drinking groundwater in the US

Importance Lithium is a naturally occurring element in drinking water and is commonly used as a mood-stabilizing medication. Although clinical studies have reported associations between receiving lithium treatment and reduced cancer risk among patients with bipolar disorder, to our knowledge, the association between environmental lithium exposure and cancer risk has never been studied in the general population. Objectives To evaluate the association between exposure to lithium in drinking groundwater and cancer risk in the general population. Design, Setting, and Participants This cohort study included participants with electronic health record and residential address information but without cancer history at baseline from the All of Us Research Program between May 31, 2017, and June 30, 2022. Participants were followed up until February 15, 2023. Statistical analysis was performed from September 2023 through October 2024. Exposure Lithium concentration in groundwater, based on kriging interpolation of publicly available US Geological Survey data on lithium concentration for 4700 wells across the contiguous US between May 12, 1999, and November 6, 2018. Main Outcome and Measures The main outcome was cancer diagnosis or condition, obtained from electronic health records. Stratified Cox proportional hazards regression models were used to estimate the hazard ratios (HRs) and 95% CIs for risk of cancer overall and individual cancer types for increasing quintiles of the estimated lithium exposure in drinking groundwater, adjusting for socioeconomic, behavioral, and neighborhood-level variables. The analysis was further conducted in the western and eastern halves of the US and restricted to long-term residents living at their current address for at least 3 years. Results A total of 252 178 participants were included (median age, 52 years [IQR, 36-64 years]; 60.1% female). The median follow-up time was 3.6 years (IQR, 3.0-4.3 years), and 7573 incident cancer cases were identified. Higher estimated lithium exposure was consistently associated with reduced cancer risk. Compared with the first (lowest) quintile of lithium exposure, the HR for all cancers was 0.49 (95% CI, 0.31-0.78) for the fourth quintile and 0.29 (95% CI, 0.15-0.55) for the fifth quintile. These associations were found for all cancer types investigated in both females and males, among long-term residents, and in both western and eastern states. For example, for the fifth vs first quintile of lithium exposure for all cancers, the HR was 0.17 (95% CI, 0.07-0.42) in females and 0.13 (95% CI, 0.04-0.38) in males; for long-term residents, the HR was 0.32 (95% CI, 0.15-0.66) in females and 0.24 (95% CI, 0.11-0.52) in males; and the HR was 0.01 (95% CI, 0.00-0.09) in western states and 0.34 (95% CI, 0.21-0.57) in eastern states. Conclusions and Relevance In this cohort study of 252 178 participants, estimated lithium exposure in drinking groundwater was associated with reduced cancer risk. Given the sparse evidence and unknown mechanisms of this association, follow-up investigation is warranted.

contiguous United States↗