USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geotechnical and Geoenvironmental Engineering”

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

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

1,722 records · Page 84Linked to original sources
USGS2026

Deep groundwater total dissolved solids mapping in the Dakota Group, Williston Basin, USA

Growing concern about the quantity of available freshwater around the world has led to interest in surveying groundwater total dissolved solids (TDS) below water well depths. Deep TDS has not been systematically mapped, and there is much to learn about the distribution and controls on deeper groundwater. In sedimentary basins across the United States, groundwater resources often overlie hydrocarbon resources, providing an opportunity to use borehole geophysical data collected for hydrocarbons to characterize groundwater and pore space resources. This study adapts a recently developed subsurface geostatistical and geophysical modeling approach to continuously map groundwater TDS, porosity, and temperature in the Dakota Group of the Williston Basin—an undercharacterized regional aquifer system overlying deeper hydrocarbon reservoirs. Groundwater TDS in the Dakota Group ranges from approximately 4800 to 26,900 mg/L. TDS patterns are stratified with higher TDS in the lower and upper Dakota Group, and relatively lower TDS in the middle Dakota Group. The lower TDS in the middle zone may represent a preferential regional flow path for lower-TDS meteoric recharge from the west. The alternating pattern of TDS may also be evidence of higher-TDS inflows into the Dakota Group from underlying and potentially from overlying aquifers. Porosity is lower near the center of the Williston Basin and tends to be higher to the east, which may be related to grain size distributions. The new regional TDS and porosity modeling serves as a quantitative reference for water users and provides supporting evidence for hypotheses on Dakota Group recharge.

Montana, North Dakota, South Dakota↗
USGS2024

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington↗
USGS2026

Subsurface vertical connectivity shapes solute transport to montane streams: Insights from watershed-scale geophysics-informed modeling

Headwater catchments provide essential water and nutrients to downstream ecosystems. The timing and evolution of their transport is shaped in part by the “invisible” subsurface structure but have been largely unexplored at the watershed scale. This raises a fundamental question: how and to what extent does subsurface structure influence streamflow and solute transport in montane watersheds? Here, we test controls of three-dimensionally resolved subsurface structure on solute transport, which reflects source waters and water residence time. We integrate the geophysically mapped architecture of Coal Creek, a headwater catchment in the Colorado River Basin, into numerical flow-transport models. Our results reveal greater vertical connectivity elongates flow paths and enhances deep groundwater contributions to streams, therefore substantially influencing the timing and magnitude of solute transport. In contrast, watershed-integrated outlet discharge can be predicted without high-resolution subsurface information. These findings demonstrate incorporating subsurface architecture is critical for predicting ecosystem health and contaminant transport.

Colorado↗
USGS2025

Earth Mapping Resources Initiative protocols—Sampling hard-rock mine waste and perpetual mine water sources

Supporting the overarching goal to evaluate critical minerals nationwide, the mine waste characterization effort in the U.S. Geological Survey (USGS) Earth Mapping Resources Initiative has created a series of protocols to standardize sampling carried out under this effort by the participating State geological surveys and their cooperators. The protocols are based on published, reviewed methods that can be deployed in the field. The protocols include (1) collecting and processing composite samples of mine and mill waste, including tailings, waste rock, gangue, heap leach piles, ore stockpiles, slag, or other mineralized and processed materials and (2) collecting and preserving water samples from perpetual or long-term mine water sources. The protocols also specify information to document on field sheets and detail the collection of geospatial data. The analytical methods used by the USGS and USGS contract laboratories are described in this report, including the data delivery pathway for USGS-derived data.

Scientific Investigations Report↗
USGS2025

Reservoir thermal energy storage pre-assessment for the United States

Storing thermal energy underground for later use in electricity production or direct-use heating/cooling is a promising, viable, and economical green energy option. Reservoir thermal energy storage (RTES) is one such option, which stores energy in underutilized permeable strata with low ambient groundwater flow rates and more geochemically evolved (e.g. brackish/saline) waters relative to overlying principal aquifer systems. The U.S. Geological Survey has begun assessing RTES potential nationally by focusing on five generalized geologic regions (Basin and Range, Coastal Plain, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Hydrogeologic reservoir models are developed for the following eight metropolitan area cities within those regions to evaluate RTES performance across different climates and subsurface conditions: Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon. Evaluated metrics include estimated required well spacing, thermal storage capacity, and thermal recovery efficiency through time. Also considered for each reservoir are potential complicating factors, including reservoir depth, thermally driven free convection, and groundwater salinity. This work focuses on direct-use cooling because the need for cooling modern office buildings greatly exceeds that for heating in most parts of the country (Falta and others, 2016); however, the evaluated metrics are also relevant to heating and electricity applications. Results indicate that favorable RTES conditions exist in each region, with the Coastal Plain and Basin and Range being especially favorable for thermal storage capacity, while the Pacific Northwest and Michigan Basin excel at energy recovery for the evaluated cooling application. The results underscore the utility of developing maps of thermal storage capacity, subsurface temperature models, and volumetric estimates of thermal storage capacity to serve as key RTES resource classification standards. Overall, this pre-assessment provides a basic understanding of RTES potential in several cities and geologic regions throughout the country and will aid ongoing thermal energy storage assessment efforts.

Arizona, Illinois, Michigan, New Mexico, Oregon, S↗
USGS2025

Megafaunal community structure on ferromanganese and phosphorite hardgrounds in the Southern California Borderland

The Southern California Borderland (SCB) is a topographically complex region on the active continental margin that hosts varied hardground habitats, including ferromanganese (FeMn) crusts and phosphorites, marine minerals being considered for resource extraction. The SCB is influenced by seasonal upwelling and terrestrial inputs, and has a well-defined oxygen-minimum zone (OMZ). We analyzed megafaunal community composition, density, and diversity across 41 ROV video transects at 10 SCB hard substrate sites spanning depth (378–2765 m), temperature (1.79–7.96 °C), and oxygen (3.91–105.67 μM) gradients, with varying rock types (FeMn crusts, phosphorites, other rocks). We counted 32,426 individuals representing 146 unique taxa over a total area of 21,935 m 2 . Echinodermata contributed 48 % of the total abundance; Cnidaria 24 %; Porifera 13 %; Annelida 6 %; Arthropoda 6 %; Chordata 2 %; Mollusca, Foraminifera, and Hemichordata <1 %. Megafauna communities showed clear heterogeneity, with density, diversity, and community composition varying among transects within sites and across sites. Rock type recorded no influence on megafaunal density, but habitats that included FeMn crusts had higher diversity and distinct taxonomic composition. Location, depth, and oxygen were the most important variables explaining variance among megafaunal communities, with distinct communities formed at deeper depths and within the OMZ. This study provides baseline information about the megafauna inhabiting SCB hardgrounds, particularly those associated with minerals considered for their resource potential. It extends existing documentation of such communities to greater depths, reveals increased representation of sponges and corals below 1000 m, and highlights the variability within and among deep-sea megafaunal communities, offering context for informed policy decisions.

Southern California borderland↗
USGS2026

Earthquake ground-motion model adjustments for the San Francisco Bay area

We develop adjustments to ergodic ground‐motion models (GMMs) to improve their performance in the San Francisco Bay Area (SFBA). GMMs are widely used in hazard assessments to estimate characteristics of ground shaking based on known properties of the source, path, and site. Such models are often developed using datasets containing records from various regions, resulting in models that represent median ground‐motion behavior, which may not adequately represent ground motions within subregions. This is true for the SFBA, where ground motions attenuate more rapidly with distance than in many other parts of California that dominate GMM databases. To support improved seismic hazard estimates in the SFBA, we calculate regional constants and anelastic attenuation coefficient adjustments relative to two commonly used ergodic GMMs: BSSA14 ( Boore et al. , 2014 ) and ASK14 ( Abrahamson et al. , 2014 ). These adjustments are obtained for a suite of ground‐motion intensity measures (peak ground acceleration, peak ground velocity, and 5%‐damped pseudospectral acceleration at oscillator periods ranging from 0.075 to 10 s) using mixed‐effects regression. Use of the regionally adjusted models reduces the overall bias by up to 0.5 natural log units for BSSA14 and up to 0.6 natural log units for ASK14. We demonstrate one application of our attenuation adjustments and their implications in an earthquake early warning case study of the 2014 M 6.0 South Napa earthquake. The predicted extent of shaking using the adjusted models better matches observed shaking at large source‐to‐site distances, especially for lower shaking intensities, thus potentially reducing overalerting. We encourage the use of our model adjustments when ergodic models are considered for seismic hazard studies in the SFBA.

California↗
USGS2024

Slip rate for the Rose Canyon fault through San Diego, California, based on analysis of GPS data: Evidence for a potential Rose Canyon–San Miguel-Vallecitos fault connection?

The Rose Canyon fault is the southern extension of the larger Newport–Inglewood–Rose Canyon fault system, which represents a major structural boundary in the Inner Continental Borderland (ICB) offshore of southern California. Ten to fifteen percent of total plate boundary motion in southern California is thought to be accommodated by the faults of the ICB, but the exact distribution of slip is uncertain. With an onshore segment, the Rose Canyon fault offers an opportunity to measure the slip rate using traditional geodetic methods. In this study, we use Global Positioning System (GPS) surface velocities from a combined campaign and continuous GPS network to constrain elastic models of the Rose Canyon fault. We then compare the observed surface velocities with proposed conceptual models of regional fault connections that facilitate the transfer of slip into the Rose Canyon fault to assess how well the observations are explained by the models. The results of elastic half‐space models suggest that the Rose Canyon fault may be slipping toward the higher end of geologic estimates, with the preferred model indicating a slip rate of 2.4 ± 0.5 mm/yr. Although limited in terms of near‐fault benchmarks, we find an improved model fit using an asymmetrical elastic half‐space model and a higher slip rate, suggesting a potential rheological contrast across the Rose Canyon fault, similar to observations from the northern Newport–Inglewood fault segments. Observed GPS surface velocities, background seismicity, and gravity anomalies south of San Diego Bay point toward a more easterly trace for the Rose Canyon fault, suggesting a possible connection with the San Miguel–Vallecitos fault system. Such a connection could increase the potential rupture lengths of future earthquakes and have important consequences for regional seismic hazards.

California↗
USGS2025

Fossil diatoms in Arctic marine surface sediments

Diatoms are one of the main constituents of marine phytoplankton in the Arctic, and thanks to their siliceous skeletons, diatom fossils are relatively well preserved in sediments. Due to their species-specific sensitivity to different ocean conditions, their abundance and assemblages in sediments are routinely used by paleoceanographers to reconstruct the state of the surface ocean in the past (Koç 2007). By using statistical methods, we can obtain valuable knowledge about their ecological preferences (Oksman et al. 2019) and generate quantitative reconstructions of various parameters, such as sea-surface temperature and sea-ice concentration through time (Krawczyk et al. 2021; Sha et al. 2014). The Marine Arctic Diatoms (MARDI) working group (WG) ( pastglobalchanges.org/mardi ) aims to advance knowledge on marine-diatom ecology and diatom-based reconstructions by compiling and harmonizing data from surface-sediment samples across the Arctic.

Pages Magazine↗
USGS2024

Parsimonious high-resolution landslide susceptibility modeling at continental scales

Landslide susceptibility maps are fundamental tools for risk reduction, but the coarse resolution of current continental-scale models is insufficient for local application. Complex relations between topographic and environmental attributes characterizing landslide susceptibility at local scales are not transferrable across areas without landslide data. Existing maps with multiple susceptibility classifications under-represent landslide potential in moderate and gently sloping terrain. We leverage an extensive landslide database ( N = 613,724), a high-resolution digital elevation model (10-m), and high-performance computing resources, to develop a new nationwide susceptibility map for the contiguous United States, Hawaii, Alaska, and Puerto Rico. We calculate four alternative linear and nonlinear thresholds of topographic slope and relief using an objective split-sample calibration. We down-sample our results to a 90-m grid to account for uncertainty in the digital elevation model and landslide position, and evaluate these thresholds' ability to differentiate areas of greater susceptibility. The less conservative nonlinear model optimally balances our priorities of capturing observed landslides (99%) while minimizing area covered by susceptible terrain (43%). Independent evaluation with four statewide landslide inventories ( N = 172,367) reinforces our model selection but highlights spatially variable performance. Therefore, we propose a novel approach to susceptibility classification using the concentration of landslide-prone terrain within each down-sampled grid. While landslides are possible within any cells containing susceptible terrain, those with the highest concentration capture the majority of observed landslides. Our new map characterizes landside susceptibility more consistently than prior models; our transparent classification approach also provides flexibility for accommodating different tolerances in risk reduction measures.

AGU Advances↗
USGS2026

Riverscape heterogeneity shapes population diversity for a migratory fish

Habitat patch dynamics can scale up to influence population demography and diversity with implications for resilience to environmental stochasticity. But how the spatial arrangement and size of habitat patches interact with other components of habitat heterogeneity to shape population diversity at larger spatial scales is not well understood. For riverine fishes, there is increasing evidence that tributary streams provide critical demographic support to main stem rivers. However, the extent to which main stem rivers rely on demographic contributions from tributaries, and the factors underlying this dependence, have not been assessed. Here, we used genetic stock identification to evaluate the effect of tributaries on population diversity of Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) occupying the main stem Snake River, Wyoming, USA. We found that the main stem relied almost entirely on tributaries for demographic support, but main stem composition varied spatially among river sections. Distance between habitat patches, catchment area, and groundwater availability acted in concert to determine the contribution of specific tributaries to the main stem, but contributions were ultimately modulated by habitat connectivity. We also found evidence for multi-scale spatial structure in tributary contributions, providing insight into untested drivers of main stem river population diversity. Our results demonstrate how spatially discrete and distributed riverscape attributes influence population diversity at broader spatial scales, illustrating how ecosystem resilience emerges from the dynamic, two-way exchange of individuals and energy across habitat networks. Management plans for large rivers that address the ecological contributions of tributaries may be needed to achieve optimal outcomes. Similarly, conservation strategies that exclusively focus on headwater streams may fail to capture the broader habitat requirements necessary to maintain robust cold-water fish populations and associated recreational fisheries, particularly under global environmental change.

Wyoming↗
USGS2025

Carbonatite-hosted residual REE deposits

Rare earth elements (REEs) occur in magmatic rocks but are especially enriched in carbonatite and alkaline silicates. If these rocks are chemically weathered, then the REEs may become further enriched within the regolith developed from these rocks. Primary magmatic REE minerals, as well as the various carbonate minerals and apatite, provide the REEs which, under pervasive chemical weathering, are incorporated within low-temperature REE minerals forming within the regolith. Many of these minerals, as well as their textures, are characteristic of this mode of formation. Lateritic conditions of weathering are instrumental in producing a thick, weathered, or regolith, profile, and the roles of sulfide oxidation, fluctuating groundwater tables, and downward mass wasting due to carbonate dissolution are identified as the most important controls on REE enrichment in the regolith.

Book chapter↗
USGS2024

Predicting hydrothermal reservoir depth from chemical geothermometers using a three-dimensional temperature model in the Great Basin, USA

Recent work in the Great Basin region of the western United States has made it possible to predict the depth of hydrothermal reservoirs (i.e., the depth at which heat is accumulated prior to ascent via hydrothermal upflow) identified through geochemistry and to contextualize the spatial patterns of these reservoir depths. Chemical geothermometers use the chemical and mineral constituents of hydrothermal fluids to predict the temperature at which fluids equilibrated with the host rocks at depth. Assuming that most of the Great Basin is dominated by conductive conditions until a vertically connected hydrothermal flow path is created (e.g., by faulting), geothermometers reflect the chemical and thermal conditions at the depth interval that the fluid has conductively equilibrated over a long period before a vertical conduit allows convective upflow. By pairing geothermometer temperature estimates with our recent three-dimensional temperature model of conductive heat flow in the Great Basin, we estimate the corresponding reservoir depths and construct a map of circulation depths. The predicted depths from geothermometers have spatial patterns across the Great Basin that relate to patterns seen in other geologic and geophysical data. Deeper springs generally occur disproportionately in areas with higher strain rates and in basins. We posit that current elevated strain rates reflect patterns of historic deformation where ongoing tectonic activity maintains permeable pathways to deeper reservoirs, some of which are estimated to exceed 6 km depth. Basins, as expected, contain a disproportionate number of these deep systems, because the underlying aquifers are closer to the surface in basins, thus requiring less water pressure to reach the surface than in mountain ranges. Most springs estimated to have their source in a deep reservoir occur at places known to host a hydrothermal system; these refined depth estimates of the source reservoir can help to better constrain the source depth for many known hydrothermal systems across the Great Basin.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2025

Ultralong, supershear rupture of the 2025 Mw 7.7 Mandalay earthquake reveals unaccounted risk

The 28 March 2025 moment magnitude ( M w ) 7.7 earthquake in Mandalay, Burma (Myanmar), ruptured 475 kilometers of the Sagaing Fault, which was more than twice the length predicted by magnitude scaling relationships. Kinematic slip models and observation of a Rayleigh Mach wave that passed through parts of Thailand confirmed that rupture occurred at supershear velocities of greater than 5 kilometers per second. The anomalous length exposed a vast population to violent near-fault shaking. The Mandalay earthquake is a modern analog for the M w 7.9 1906 San Francisco earthquake, another atypically long and fast rupture. Probabilistic seismic hazard analyses use scaling relations that do not account for such long ruptures at moderate magnitudes. This limitation, in conjunction with a likely increased population and infrastructure exposure for atypically long ruptures, contributes to a potential mischaracterization of seismic risk.

Mandalay↗
USGS2025

Ambient field seismology in critical zone hydrological sciences

Passive ambient noise monitoring is an emerging tool in environmental seismology, leveraging the ambient seismic field to assess temporal variations in shallow subsurface properties. This review focuses on the potential and challenges of using scattered coda waves from noise correlation functions to monitor critical zone dynamics. The sensitivity of seismic velocities to various environmental factors, including precipitation, snowmelt, atmospheric pressure, and groundwater fluctuations, underscores the method’s versatility. While coda waves excel in detecting subtle changes due to their scattered nature, ballistic waves provide higher spatial resolution, albeit with challenges in source stability. Advances in seismic sensing, including distributed acoustic sensing and low-cost geophone networks, have enabled high-resolution monitoring of hydrological processes, subsurface deformation, and seismic hazards. Integrating seismic data with hydrological models provides insights into water storage, pore pressure changes, and soil moisture dynamics. However, limitations in spatial resolution, calibration with ground truth data, and coupled effects between environmental factors remain key challenges. This review emphasizes the importance of interdisciplinary approaches in refining methodologies, enhancing sensor deployments, and addressing data gaps. Passive seismic monitoring offers opportunities to understand critical zone processes and their broader impacts on seismic hazards and environmental sustainability.

Comptes Rendus. Géoscience↗
USGS2025

Unlocking ecological insights from sub-seasonal visible-to-shortwave infrared imaging spectroscopy: The SHIFT campaign

We stand at the threshold of a transformative era in Earth observation, marked by space-borne visible-to-shortwave infrared (VSWIR) imaging spectrometers that promise consistent global observations of ecosystem function, phenology, and inter- and intra-annual change. However, the full value of repeat spectroscopy, the information embedded within different temporal scales, and the reliability of existing algorithms across diverse ecosystem types and vegetation phenophases have remained elusive due to the absence of suitable sub-seasonal spectroscopy data. In response, the Surface Biology and Geology (SBG) High-Frequency Time Series (SHIFT) campaign was initiated during late February 2022 in Santa Barbara County, California. SHIFT, designed to support NASA's SBG mission, addressed mission scoping, scientific advancement, applications development, and community building. This ambitious endeavor included weekly Airborne Visible InfraRed Imaging Spectrometer-Next Generation (AVIRIS-NG) imagery acquisitions for 13 weeks (spanning February 24 to May 29, 2022), accompanied by coordinated terrestrial vegetation and coastal aquatic data collection. We describe the rich datasets collected and illustrate how the complex sub-seasonal patterns of change can be linked to biological science and applications, surpassing insights from multispectral observations. Leveraging open-source processing methods and cloud-based analysis tools, the SHIFT campaign showcases the readiness of the scientific community to harness ecological insights from remotely sensed hyperspectral time series. We provide an overview of SHIFT's goals, data collections, preliminary results, and the collaborative efforts of early career scientists committed to unlocking the transformative potential of high-frequency time series data from space-borne VSWIR imaging spectrometers.

Ecosphere↗
USGS2026

Microtextural characteristics of adularia in banded quartz veins from the Midas low-sulfidation epithermal deposit, Nevada

High-grade ores at the Miocene Midas low-sulfidation epithermal deposit in northern Nevada are confined to crustiform quartz veins containing abundant adularia. Micro-X-ray fluorescence elemental mapping reveals that adularia is a common gangue mineral occurring in colloform bands, bands showing bladed textures, and bands with dendritic terminations. The adularia aggregates have delicate shapes and are comprised of stacked, submillimeter crystals hosted by fine-grained quartz. The textural evidence suggests that the adularia aggregates originally formed within a gel-like, noncrystalline silica matrix, which subsequently transformed into quartz. This indicates that the adularia did not precipitate in open space along the vein walls. Correlative microscopy, involving scanning electron microscopy-based automated mineralogy and optical petrography, demonstrates that bands containing abundant adularia are not the primary host to ore minerals. The ore minerals occur in different bands within the crustiform veins, implying that adularia and ore mineral precipitation did not always occur simultaneously. It is hypothesized here that fluid flow at Midas involved intermittent short-lived events of fluid flashing, causing rapid solute supersaturation in the liquid. During each flashing event, different amounts of vapor were produced along a given vein. Compositional differences between adjacent bands in the crustiform quartz veins may, therefore, be linked to variations in the amount of vapor formed during each flash event.

Nevada↗
USGS2025

Disparate groundwater responses to wildfire

Post-wildfire investigations of groundwater response reveal a range of outcomes, varying from substantial increases to notable decreases in recharge and baseflow, with some studies indicating negligible or short-lived effects. This review assesses these varied responses within five critical categories: climate, vegetation, hydrogeology, fire characteristics, and the cryosphere, examining both short-term (within 2 years) and intermediate (2–10 years post-fire) effects. Despite considerable variability, some consistent patterns emerge. For instance, in hydroclimatic settings where water input and evaporative demand cycles are out of sync, post-wildfire groundwater responses tend to be positive (i.e., increased flux or storage), whereas under low fire severity conditions or in vegetation types that quickly recover, groundwater responses tend to be negative (i.e., decreased flux or storage). We synthesize relevant findings into a compendium of testable hypotheses aimed at explaining the spatiotemporal variability in observed post-wildfire groundwater responses. A recurring theme is the critical influence of the pre-wildfire groundwater regime on expected response and recovery. We identify opportunities for specific improvements in post-wildfire monitoring and modeling that would further advance capabilities to predict groundwater response. A key area for further research is understanding how wildfire effects on snow dynamics and other cryospheric processes translate to changes in groundwater.

WIREs Water↗