Geology topics

Browse USGS reports

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

At least 1,531 records · Page 85Linked to original sources

USGS2025

Using visualization science to inform the design of environmental decision-support tools—A case study of the U.S. Geological Survey Waterwatch

Environmental decision-support tools are increasingly being used to serve both expert and non-expert audiences. Many existing tools are primarily expert-focused, and redesigning them can be challenging because experts and non-experts interact with tools differently, existing users may be resistant to changes, and there is little guidance on how to prioritize redesign efforts and demonstrate their efficacy. In this report, we present a case study of a user-centered redesign of an established environmental decision-support tool—the U.S. Geological Survey WaterWatch. WaterWatch supports flood, drought, and other water resource management decisions through the display of water levels at gages across the United States. Using a participatory process, we identified a functional change (replacing the existing rainbow colormap), created an alternative design, and tested the alternative’s usability through two general public surveys. The results showed that replacing the rainbow colormap with a more intuitive diverging colormap improves usability, regardless of the audience’s subjective preference for the rainbow color scheme. In addition, we demonstrated the importance of using legends to improve the audience’s understanding of the map symbols. This study demonstrates how user-centered design approaches can be used to inform the design of high-profile products and tools.

Scientific Investigations Report↗

USGS2025

A 10-year continuous daily simulation of chloride flux from a suburban watershed in Fairfax County, Virginia, USA

Increasing levels of chloride in surface water are associated with detrimental effects on water quality, aquatic ecosystems, infrastructure, and human health. Numerous mass-balance studies have inferred watershed transport processes by interpreting chloride inputs and outputs, but few represent internal dynamics explicitly. We constructed a coupled water/chloride mass balance model to gain insights into storage, residence time, and transport processes in a 10-km 2 urban watershed. The model, which operates over a 10-year period at a daily time scale, represents storage in a dynamic soil-moisture reservoir, quick-flow runoff from storm events, and slow-flow runoff that sustains streamflow in dry weather. The calibrated model accurately represented (a)the observed transition from a streamflow enrichment regime in cold months to a dilution regime in warmer months, (b) the observed tendency for late-summer concentrations to be higher after winters with heavy snowfall, and (c) a period-of-record downward trend in chloride concentration likely associated with a downward trend in annual snowfall. Estimated chloride inputs averaged 195 metric tons per year, while the average output was 270 metric tons per year. In contrast, estimated storage was only 107 metric tons. The estimated mean residence time in groundwater was 1.27 years. This short residence time indicates that efforts to reduce inputs will manifest as decreased concentrations in streamflow on a management-relevant time scale of several years. The coupled mass balance model yielded insights into internal watershed dynamics that would not be possible from simple input/output analysis; such models can be useful tools for gaining insight into small watershed hydrology and pollutant transport.

Virginia↗

USGS2025

MTAB 113, December 2025

This Memo to All Banders (MTAB 113) was released in December 2025. Subjects in this this memo are 1. The Chief’s Chirp – End-of-Year Message and BBL’s Year in Numbers; 2. Alerts –Highly Pathogenic Avian Influenza and Screw-worms; 3. News – GameBirds Data Release and Banders Without Borders Updates; 4. A note from the permitting shelves – project description template and double-check contact info reminder; 5. A note from the supply room – band order reminders; 6. Data management – NABBP Database Species Changes Updates and prepare your data reminders; 7. Frequently asked questions – What bird status extra info code should I use? How can I process repeat warnings faster during the data submission process?; 8. Auxiliary marker corner – reminder to get data in; 9. Message to the Flyways – Data submission, band orders, and Winter Flyway Council Meetings; 10. Moments in history – bird band pop quiz; 11. Upcoming events; 12. Recent literature; 13. Request for information; and 13. Appendix: NABBP Database Species Changes Update.

Memo to All Banders (MTAB)↗

USGS2025

The value of reducing uncertainties to support the management of a high‐elevation endemic salamander

Many salamander populations are declining, and methods to determine how best to allocate limited resources to slow or reverse these declines could support land managers in their decision‐making processes. Multiple types of uncertainty may delay management decisions, including when (1) knowledge of a species' ecology is incomplete, (2) climate change effects on environmental covariates are uncertain, and (3) the efficacy of management alternatives is unknown. For management decisions, a value‐of‐information analysis can identify which uncertainties are critical to reduce in order to identify an optimal strategy from a set of possible management actions. If the same management action is optimal across the full range of uncertainties, then resources for research can be redirected toward active management. Using value‐of‐information analyses, we examine the effect of uncertainty on identifying optimal management to maximize the future expected occupancy of Plethodon shenandoah , a Federally Endangered high‐elevation endemic salamander that is threatened by climate change. Out of 11 management actions proposed by National Park Service managers, those that increase environmental moisture are expected to maximize occupancy, and we find that the selection of this action is robust to all the identified uncertainties. We show that, even in systems with multiple sources of large uncertainty, value of information analyses discriminate among investments in species management.

Ecosphere↗

USGS2025

Big River bluegill: Combining vital rates and long-term monitoring to understand population dynamics in large rivers

Long-term monitoring indicates bluegill catch rates are relatively stable in some reaches of the Upper Mississippi River and highly variable in others, whereas in the Illinois River, catch rates have decreased. A lack of age structure information precludes understanding population processes responsible for patterns in catch rates. To build a better understanding of why catch rates have changed over time, we integrated short-term age structure information with long-term monitoring data to quantify and assess spatial patterns in bluegill population dynamics across six study reaches of these two rivers. Specifically, we estimated and compared age and size structure, growth, maturity, mortality, and recruitment. We used quantile regressions to apply age estimates to long-term data for investigating trends in age-based catch rates reflective of recruitment (age-1 catch rates), mortality (using age-1 and age-2+ catch rates), and spawning stock (age-2+). Our findings indicate trends in bluegill age-2+ catch rates increased and then stabilized across upstream study reaches, but dynamic rates, size structure, and age at maturity varied among reaches. Bluegill populations in downstream study reaches had low maximum size, early maturation, low mean age, low proportional stock density, declining recruitment, and declining age-2+ catch rates. An insufficient number of bluegills were collected from the study reach furthest downstream to adequately quantify dynamic rates. Our results support life history theory in that bluegill respond to unstable environmental conditions through life history adaptations. These findings show how integrating periodic age structure information with long-term monitoring can enhance population assessments.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗

USGS2025

Evaluating uncertainties with sample-collection method and source selection in sediment fingerprinting: an example from a Great Lakes tributary

Purpose Identifying suspended-sediment (SS) sources, seasonal variability, and phosphorus (P) transported with SS is critical information for basin managers, although there may be concerns about comparability between flow-integrated SS samples used for sediment fingerprinting and discrete samples used for concentrations and loads in basins where SS is mostly silt + clay and(or) one land-use predominates. Objectives were to determine if (1) sample-collection method and (2) source consideration influenced apportionment of the largest SS source. Methods Concurrent-replicate, SS samples were collected during 2022 from the East River, Wisconsin using an automated sampler, commonly used for water-quality sampling, and passive samplers, frequently used for SS fingerprinting. Samples were evaluated for differences in physical and chemical characteristics that may affect source apportionment. Considered sources included three upland land-use (cropland, forest, and roads), two in-channel (streambank and streambed sediment), and one that connects uplands to the stream channel (gullies). Source apportionment used established methods in the SedSAT tool. Source scenarios included land-use + streambank (4src), 4src + gully, 4src + streambed, and 4src + gully + streambed (6src). Results There were no statistically significant differences in median grain size, organic carbon, or sediment-bound P as a function of collection method. In-channel sources were the largest proportional SS source, regardless of season, hydrologic condition, collection method, or source scenario. Source verification highlighted which source fingerprints were most accurately defined and implications for SS target apportionment. Conclusion Varying the source scenarios for sediment fingerprinting indicated that improved management of hydrologic connectivity between upland land use and the stream channel has the potential to mitigate SS loads.

Wisconsin↗

USGS2025

Assessing environmental drivers and protist community dynamics that shaped the historic August 2022 Heterosigma akashiwo bloom in San Francisco Bay, California

San Francisco Bay, California, typically has chlorophyll a (chl- a ) concentrations below 10 µg L −1 , despite nutrient loadings exceeding those in many estuaries with recurring harmful algal blooms (HABs). However, in August 2022 there was a Heterosigma akashiwo (raphidophyte) bloom with chl- a concentrations exceeding 450 µg L −1 , resulting in widespread hypoxia and fish die-off. We used protist community (18S ribosomal ribonucleic acid gene metabarcoding) and environmental data collected 2–3 times per month from 2015–2023 to assess differences between 2022 and other years, and capture dynamics pre-, during, and post-bloom. The 2022 protist community clustered separately from other years in non-metric multidimensional scaling and hierarchical cluster analyses. A diatom bloom in June generated the highest median June chl- a concentration (15 µg L −1 ) recorded between 1977–2023 and drew down nutrients (nitrogen, silicate), leading to the establishment of a heterotrophic community in July with low dinoflagellate abundance (0.9% of protists versus 5.9% average in other years). The absence of dinoflagellates, potentially including H. akashiwo predators , plus suitable habitat (13–14 daylight hours, temperatures 17–21˚C, nitrate concentrations 17–24 µmol L −1 , neap tides) likely provided opportunity for H. akashiwo to proliferate from 0.03% to 48% of protists between July 20 and August 25. As the bloom collapsed, the genus Gyrodinium (dinoflagellate predator of H. akashiwo ) increased from 3.0% to 52% of protists between August 25 and August 31, likely exerting grazing pressure. Results show that antecedent protist community dynamics play important roles in shaping HABs, and interactions between biotic and abiotic factors govern the development of blooms.

California↗

USGS2025

Where to restore and conserve? A regional benefit cost analysis of coral reef protection and restoration for coastal flood resilience

Momentum is growing for the management of coral reefs as a strategy to reduce climate risks in tropical coastlines. Yet, quantification of the life-time costs, impacts, and benefits remains limited. This study provides one of the first rigorous, spatially explicit, regional-scale Benefit:Cost Analyses (BCA) for coral reef restoration and was designed to meet the BCA requirements of major hazard mitigation funding and programs. This study simulates coastal flooding using a hydrodynamic model under different scenarios representing current coral reef conditions, reef degradation, and reef restoration. These coastal flood maps are used to estimate socioeconomic damages, which are included in a BCA to assess cost-effectiveness and priority areas for coral reef conservation and restoration. The United States Virgin Islands is used as a case study given recent impacts from storms and their new policy that declares reefs as natural infrastructure. The results show that flood risk across the islands of Saint Croix, Saint John, and Saint Thomas is $51.4 million USD per year. Annually, coral reefs prevent flood damages to 481 people and $43.6 million USD of infrastructure, which represents 87 % of the flood risk. These results identify the communities that could receive the greatest benefits from coral conservation and restoration, which overlap significantly. Coral reef restoration is found to be a cost-effective strategy for flood mitigation although how reefs are restored matters. Hybrid coral reef restoration provides greater flood mitigation benefits than ecological restoration, but its cost-effectiveness is overall lower given its greater costs. Many coastal areas have benefit:cost ratios exceeding 1.5, which complies with typical government agencies’ requirements for accessing hazard mitigation funds. This valuation framework helps provide rigorous regional-scale quantification of nature-based coastal protection solutions for coastal risk management decisions.

Journal of Environmental Management↗

USGS2025

Speciation genomics in the tiger whiptail lizards ( Aspidoscelis tigris complex)

The transition from small genetic to genome-scale datasets for studying biodiversity has revealed that genetic exchange through introgressive hybridization is a widespread phenomenon in nature. Despite this, a lack of high-quality reference genomes for most non-model species limits our understanding of the impact of this process for many taxonomic groups. This restricts the range of insights that genomic tools can provide for conservation biologists, who often hope to employ genomic datasets to accurately identify historically isolated lineages to protect and to predict their evolutionary fate in the face of environmental change. Tiger whiptail lizards ( Aspidoscelis tigris complex) are an abundant and important ecological component of ecosystems across the southwestern United States. In this study, we assembled and annotated a chromosome-level reference genome for A. t. stejnegeri from coastal California. We then used this reference genome to reconstruct patterns of speciation and admixture within the larger species complex, finding evidence that gene flow is widespread both geographically and across the genome.

northern Mexico, southwestern United States↗

USGS2025

An analysis of the linked decisions in the confiscation of illegally traded turtles

Over the last few decades, freshwater turtles have become more common in the illegal wildlife trade because of growing global demand. Illegally traded turtles may be intercepted by several different agencies with separate jurisdictions. When turtles are confiscated, uncertainties may make releasing them back to the wild difficult. We used tools from decision analysis to achieve the following three objectives: (1) map elements of the decision process and their relationships in the illegal turtle trade using conceptual models, (2) outline the linked decisions for turtle confiscation and repatriation using decision trees, and (3) evaluate the decision trees for two example scenarios, one with moderate uncertainty and one with high uncertainty. We used the wood turtle ( Glyptemys insculpta ) as a case study, which is a species of conservation concern in part due to illegal wildlife trafficking. We conducted 23 semi-structured interviews of decision makers in law enforcement, biologists, land managers, and zoo staff. Interviews revealed that decisions regarding the disposition of confiscated turtles are complicated by uncertainty in disease status and origin. Decision makers that handle confiscated turtles also recognize that their decisions are often made in sequence and dependent on the outcome of antecedent decisions. In evaluating our decision trees, we found that the optimal decisions for example scenarios were similar and insensitive to uncertainty. Future applications of the decision trees by decision makers would involve a decision analyst to parameterize and interpret the choices and consequences involved in working through these decision trees. Collectively, our work shows how the use of decision trees can help structure and evaluate risky decisions for repatriating confiscated wildlife.

Conservation Science and Practice↗

USGS2025

How do hydrological variability and human activities control the spatiotemporal changes of riverine nitrogen export in the Upper Mississippi River Basin?

Excessive nitrogen export from agricultural watersheds remains a critical water quality challenge, with the Upper Mississippi River Basin (UMRB) significantly contributing to downstream eutrophication and hypoxia in the Gulf. This study investigates the spatiotemporal dynamics of riverine nitrate plus nitrite (NO 3 – + NO 2 – -N) export across the UMRB at high spatial resolution (12-digit Hydrologic Unit Codes or HUC12 subwatershed scale) during 2001–2020 and quantifies the effects of anthropogenic activities and hydrological variability on riverine NO 3 – + NO 2 – -N export changes in the region between 2001–2005 and 2016–2020. Our results revealed hotspots of substantial increases in NO 3 – + NO 2 – -N yields across the UMRB, with distinct regional patterns in driving factors. Over the entire UMRB, NO 3 – + NO 2 – -N yields increased by 9.7 kg/ha/yr on average from 2001–2005 to 2016–2020, with anthropogenic activities contributing 4.8 kg/ha/yr and hydrological variability contributing 4.9 kg/ha/yr. The northern and western UMRB had combined influences from both anthropogenic activities and hydrological variability, while the east-central regions had predominantly hydrologically driven changes. Agricultural sources, including fertilizer, manure, and biological nitrogen fixation, collectively contributed over 80% of NO 3 – + NO 2 – -N loading throughout the basin. This framework for disentangling human and hydrological impacts provides critical insights for developing effective and targeted watershed management strategies to reduce nutrient losses and improve water quality.

Upper Mississippi River Basin↗

USGS2025

Source(s) of the smooth Caloris exterior plains on Mercury: Mapping, remote analyses, and scenarios for future testing with BepiColombo data

Mercury hosts widespread smooth plains that are concentrated in the Caloris impact basin, in an annulus surrounding the Caloris basin, and in the adjacent northern smooth plains. The origins of these smooth plains are uncertain, although prior work suggests these plains in the northwestern Caloris annulus might reflect volcanic activity, impact ejecta, or a combination of the two. Deciphering the timing and mode of emplacement of these plains would provide a critical constraint on regional late-stage volcanism or impact effects. In this work, the region northwest of Caloris was investigated using geomorphological and color-based mapping, crater counting techniques, and spectral analyses with the goal of placing constraints on the source of the observed units and identifying the primary emplacement mechanism. Mapping and spectral analyses confirm previous findings of two distinct, yet intermingled, units within these plains, each with similar crater count model ages that postdate the formation of the Caloris impact basin. Mapping, spectra analysis, ages, and the identification of potential flow pathways are more consistent with a predominantly volcanic origin for the smooth plains materials, although these data do not rule out contributions from impact ejecta or impact melt. We propose several hypothetical scenarios, including post-emplacement modification by near-surface volatiles, to explain these observations and clarify the emplacement mechanism for these specific smooth plains regions. Further observations from the BepiColombo mission should provide data to potentially address the outstanding questions from this work.

Remote Sensing↗

USGS2025

Correction to “Stable isotope reference materials and scale definitions – outcomes of the 2024 IAEA experts meeting”

This article corrects the following: Camin, F., Besic, D., Brewer, P., Allison, C., Coplen, T., Dunn, P., Gehre, M., Gröning, M., Meijer, H.-A.-J, Hélie, J.-F., Iacumin, P., Kraft, R., Krajnc, B., Kümmel, S., Lee, S., Meija, J., Mester, Z., Mohn, J., Moossen, H., Qi, H., Skrzypek, G., Sperlich, P., Viallon, J., Wassenaar, L. and Wielgosz, R. (2025). Stable Isotope Reference Materials and Scale Definitions—Outcomes of the 2024 IAEA Experts Meeting. Rapid Commun Mass Spectrom, 39: e10018, https://doi.org/10.1002/rcm.10018. In the above article, the text of note 3 in Section 4.2 is incorrect. The correct text is shown below: 3. Whilst limited amounts of the original VSMOW and SLAP scale-defining reference materials may still be available, to enable continuous wide distribution of RMs for VSMOW-SLAP scale realization (normalization, calibration) the scale-maintaining reference materials VSMOW2 and SLAP2 are available. VSMOW2 and SLAP2 have δ(18O) values identical to VSMOW and SLAP, respectively, but with an uncertainty with respect to VSMOW and SLAP. We apologize for this error. The CIAAW-IUPAC endorsed all isotope delta scale definitions appearing in Camin et al. (2025) and agreed to publish the current and all future updates and corrections to these definitions on its website (www.ciaaw.org).

Rapid Communications in Mass Spectrometry↗

USGS2025

Insights into widespread landsliding in southern Appalachia from Hurricane Helene

Between 23 and 27 September 2024, antecedent rain followed by Hurricane Helene produced one of the most damaging weather events in southern Appalachia history. The back-to-back storm events resulted in a maximum cumulative rainfall of 848 mm and hurricane-force wind gusts over 170 km/h in western North Carolina, eastern Tennessee, and southwestern Virginia. The resulting regional flooding, landslides, and tree blowdown caused over 100 fatalities, damaged or destroyed critical infrastructure and thousands of structures, and severed connectivity across the region. Over the next several weeks, a multi-agency landslide response produced a rapid hazard assessment and mapped 2217 landslides, 55% of which damaged infrastructure or property. Orographic uplift enhanced rainfall, resulting in concentrated landsliding along the ~250 km swath of the Blue Ridge escarpment in western North Carolina. Landslides initiated predominantly on windward-facing (southeast-facing) slopes, and localized clustering of initiation points indicated a strong influence of hillslope-scale meteorological and geomorphic factors. Many shallow landslides mobilized into larger, highly mobile, and damaging debris flows that graded into floods. Here, we put our preliminary observations in the context of historical storm-driven landslide events and open new avenues for investigating the nature and extent of landslides and their effects in southern Appalachia and similar environments.

Georgia, North Carolina, South Carolina, Tennessee↗

USGS2025

Regional characterization of coal resources in the U.S. Gulf Coast

There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The U.S. Gulf Coast and the Williston basins are the two major lignite-bearing basins within the country. Recent REE and CM studies of the lignite in these basins have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of the U.S. Gulf Coast lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from 31,181 drill holes and other data points. The estimated total coal resource in the Gulf Coast is about 83 billion metric tons in the upper 90 m (~ 300 ft) of the subsurface. Texas accounted for 40 percent (32 billion metric tons) of the total resource, followed by Mississippi (24 %, 20 billion metric tons), Louisiana (14 %, 12 billion metric tons), Tennessee (10 %, 8.5 billion metric tons), and Arkansas (6 %, 5.1 billion metric tons). The remaining states each accounted for less than 5 percent of the total resource. Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known state-wide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast and Williston Basin resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses. Coal-power electric generation has sharply decreased in past decades and mining of these coals for CM and REE could provide additional co-products such as activated carbon and other uses such as fertilizer (soil enhancer).

EarthArXiv↗

USGS2025

Geologic map of the Guinevere Planitia quadrangle (V–30), Venus

This 1:5,000,000-scale geologic map of the Guinevere Planitia quadrangle divides the region into 15 geologic material units, defined using Magellan synthetic aperture radar (12.6-centimeter-wavelength radar system; 75 meters per pixel) datasets and including upland terrain units (2.4 percent of the surface area), plains materials units (59 percent), flow materials associated with named and unnamed eruptive centers (37.2 percent), small volcanic edifices, and impact crater materials (1.4 percent). Upland terrain units consist of tessera and lineated upland material, plains materials consist of Guinevere regional plains and Guinevere lineated and mottled plains, and flow materials consist of lobate flow material and plains-forming flow material. Specific lobate flows associated with Atanua Mons, Tuli Mons, Var Mons, and Uilata Fluctus are mapped separately. Other mapped units are impact crater material and small volcanic edifice. In addition to geologic units, we mapped linear features that show patterns of deformation or flow across the quadrangle. These consist of faults, wrinkle ridges, broad arches, channels, troughs, and flow direction indicators. The map region also contains several small volcanic features: shields, depressions, and craters. These, in combination with the plains, large volcanoes, and coronae, show the pervasive influence of volcanism across Venusian lowlands. The rims of nine identified impact features are delineated; large bright and dark haloes, which in some cases are associated with individual impact craters, are mapped as surficial mantling deposits. We documented spatial relationships using the stratigraphic and cross-cutting relationships of the quadrangle’s geologic units and features to provide a synthesis of the region’s geologic history. The upland terrain of the quadrangle indicates intense tectonic deformation and uplift. It is exposed as embayed remnants, typically within the plains, and represents the oldest geologic materials locally and across the region. Guinevere plains and the plains-forming flow unit appear to be assemblages of volcanic flows from multiple sources, including distinct coronae and corona-like structures. The temporal evolution of Guinevere lineated and mottled plains was likely protracted, with continued formation of small volcanic edifices over a long period. The morphologic and radar brightness characteristics of volcanoes in the region indicate their growth may have involved (1) multiple large-scale eruptive centers with recognizable spatial and temporal sequences, (2) extensive lava flow fields with a multitude of flows producing complex, overlapping patterns, and (3) numerous small volcanic edifices, including shields, domes, and cones. Although geologic patterns common to other regions of Venus are evident in the Guinevere Planitia quadrangle, local relative age relationships are inconsistent or unclear, preventing robust stratigraphic correlation. The mapping results do, however, indicate complicated local sequences of volcanic and tectonic activity.

Scientific Investigations Map↗

USGS2025

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

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

California↗

USGS2025

A probabilistic assessment methodology for the evaluation of geologic energy storage capacity—Natural gas storage in depleted hydrocarbon reservoirs

The need for energy storage, particularly underground, where capacity and duration may far exceed battery storage technologies, is especially relevant given the increasing demands for reliable power alongside the development of intermittent renewable electricity sources. Geologic energy storage facilities already exist, and expanded use would enable storing gases such as methane and hydrogen. In 2018, a National Academies of Sciences, Engineering, and Medicine report, “Future Directions for the U.S. Geological Survey's Energy Resources Program,” recommended that the U.S. Geological Survey (USGS) prioritize assessing underground energy storage in geologic formations in the United States. The U.S. Geological Survey has since developed a methodology for assessing natural gas storage capacities in depleted hydrocarbon reservoirs on a national scale. The methodology introduced in this report prescribes three approaches for calculating gas storage capacity. This methodology relies on the availability of input data, including cumulative hydrocarbon production records, reservoir petrophysical properties, and reservoir pressure data. Assessment inputs can be obtained from public, State-level databases and propriety national-scale databases, although the use of analogs could be warranted for estimating input parameters. Probabilistic assessment results are aggregated to play, petroleum province, regional, and national scales. The steps defined in this report are demonstrated on the Michigan Basin Province, which includes the Mississippian Sandstone Gas Play and the Clinton Structural Play. This methodology could be used to systematically and consistently assess hydrocarbon plays and provinces for natural gas storage capacity across the United States.

Michigan↗