Craig Ulrich

Managed aquifer recharge in urban basins increasingly relies on small-footprint stormwater control measures (SCMs), yet their actual infiltration performance and contribution to groundwater recharge remain poorly constrained after installation. This study evaluated two contrasting SCMs in the Los Angeles Basin — a deep drywell beneath a paved parking lot and a shallow bioswale within a narrow right-of-way — using controlled recharge experiments monitored with borehole-based, time-lapse electrical resistivity tomography (ERT). Instrumented boreholes were installed around each SCM to enable 3D (drywell) and 2D (bioswale) imaging, and ERT data were collected autonomously over several days, during which known volumes of water were applied to each system for calibration. Soil moisture was independently measured with neutron probes and in situ sensors to validate resistivity-derived wetting patterns. Results showed that the drywell rapidly conveyed water from shallow pretreatment features and deeper screened intervals into laterally extensive and stratigraphically controlled pathways, with preferential flow observed along higher-permeability layers. By contrast, infiltration in the bioswale was limited to the upper 2 m, even under elevated water levels exceeding typical storm conditions, indicating shallow storage and limited deep percolation. In both settings, changes in resistivity closely matched changes in measured soil moisture, confirming the ability of borehole ERT to resolve hydrologic processes in constrained built environments. The findings demonstrated that SCM infiltration performance is highly site specific, influenced by engineering design and subsurface heterogeneity, and that geophysical monitoring provides critical information for accurate recharge estimation, SCM efficacy and siting, maintenance planning, and long-term water-resource management in urban basins.

Mountainous headwater streams are sustained by both snowmelt-driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model-data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface-subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed datasets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater level. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process-based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Groundwater interactions with mountain streams are often simplified in model projections, potentially leading to inaccurate estimates of streamflow response to climate change. Here, using a high-resolution, integrated hydrological model extending 400 m into the subsurface, we find groundwater an important and stable source of historical streamflow in a mountainous watershed of the Colorado River. In a warmer climate, increased forest water use is predicted to reduce groundwater recharge resulting in groundwater storage loss. Losses are expected to be most severe during dry years and cannot recover to historical levels even during simulated wet periods. Groundwater depletion substantially reduces annual streamflow with intermittent conditions predicted when precipitation is low. Expanding results across the region suggests groundwater declines will be highest in the Colorado Headwater and Gunnison basins. Our research highlights the tight coupling of vegetation and groundwater dynamics and that excluding explicit groundwater response to warming may underestimate future reductions in mountain streamflow.

We analyze results of high-resolution seismic imaging at The Geysers geothermal reservoir in northern California, USA, using a dense seismic network to image the spatial heterogeneity of the reservoir structure and flow paths. The project uses 92 seismic stations spaced at approximately 500 m over a 5 km × 5 km study area. Microseismic data for more than 17,000 earthquakes have been acquired over a period of 13 months and automatically processed for P- and S-wave phase arrival times. The data are subsequently inverted using a joint inversion approach to image the spatial heterogeneity of the reservoir including the 3D P- and S-wave velocity structure and VP/VS ratio, and to locate earthquake hypocenters. The resulting tomographic images are appraised by integration into The Geysers’ 3D reservoir model and by a spatial correlation to the injection and production wells. Spatial correlation of P-wave velocity images to water injection and steam production wells reveal higher velocities below injection wells, due to higher water saturation, and lower velocities in the vicinity of steam-producing wells, due to the presence of steam in the surrounding reservoir rocks. The spatial correlation of VP/VS to steam in the reservoir indicates decreased estimates in the vicinity of steam production wells. In contrast, the VP/VS ratio reveals high values in the reservoir for regions near water injection wells and along the potential flow path of water through the reservoir. The estimates of shear modulus indicate high values in a region of competent graywacke, which is known for a lack of fractures and steam production, and low values in a region that is dominated by water flow, suggesting that fractured rock created a pathway for the water through the reservoir. The heterogeneity observed in the S-wave velocities indicates a compartmentalized reservoir, which correlates spatially with fault projections in the 3D reservoir model.

Bedrock property quantification is critical for predicting the hydrological response of watersheds to climate disturbances. Estimating bedrock hydraulic properties over watershed scales is inherently difficult, particularly in fracture-dominated regions. Our analysis tests the covariability of above- and belowground features on a watershed scale, by linking borehole geophysical data, near-surface geophysics, and remote sensing data. We use machine learning to quantify the relationships between bedrock geophysical/hydrological properties and geomorphological/vegetation indices and show that machine learning relationships can estimate most of their covariability. Although we can predict the electrical resistivity variation across the watershed, regions of lower variability in the input parameters are shown to provide better estimates, indicating a limitation of commonly applied geomorphological models. Our results emphasize that such an integrated approach can be used to derive detailed bedrock characteristics, allowing for identification of small-scale variations across an entire watershed that may be critical to assess the impact of disturbances on hydrological systems.

We numerically study how the interactions between surface, subsurface, and rainfall create complex runoff behavior in a headwater-dominated high elevation, mountainous catchment in East River Watershed, Colorado, USA. In order to understand the effect of surface--subsurface interactions on the hydrological response, we compare model variations with different soil permeability. We compute hydrographs, hydrological indices, and spatio-temporal distributions of hydrological variables. The indices and distributions are then linked to the hydrograph at the outlet of the domain. Our results show that changing surface--subsurface interactions trigger different flow regimes, connectivity dynamics, and runoff generation mechanisms inside the catchment, and hence, affect the distributed hydrological response. Increased surface--subsurface interaction leads to a higher degree of connectivity in the catchment. Temporal dynamics of both distributed and aggregated hydrological response did not change significantly, which may be because the structural connectivity of the catchment remained untouched. Our results indicate that in higher order catchments, the effect of changes in distributed hydrological response may not always be observable through aggregated hydrological signatures.

Measuring hydro-mechanical properties of natural fractures is a prerequisite for optimizing hydraulic stimulation design and well placement. We completed experiments to characterize shear on natural fractures in schist, amphibolite, and rhyolite specimens drilled from EGS Collab Project’s field sites at the Sanford Underground Research Facility (SURF) in South Dakota. A triaxial direct shear method and coupled x-ray imaging were used to perform hydroshearing and mechanical shearing at the site’s in-situ stress conditions. This produced simultaneous measurements of fracture and matrix strength, permeability, stress-dependent aperture, dilation, and friction strength. Our results identified that only a subset of the natural fractures was weak enough for hydroshearing. Generally, hydroshearing increases fracture permeability by a factor of 10 or more and the enhancement is retainable over time. However, the shear slip does not always result in permeability enhancement. High content of phyllosilicates was found to associate with exceptionally weak fractures that also exhibited poor or even negative enhancement after stimulation. Combining our measurements with site data, we can predict that most observable fractures at the two EGS Collab sites do not meet the criteria for hydroshearing before tensile opening. In some cases, the visible fractures are low permeability and as strong as the adjacent rock. To induce hydroshearing before tensile opening, injection must target known weak and favorably oriented fractures with confirmed pre-existing permeability.

Understanding the initiation and arrest of earthquakes is one of the long-standing challenges of seismology. Here we report on direct observations of borehole displacement by a meter-sized shear rupture induced by pressurization of metamorphic rock at 1.5 km depth. We observed the acceleration of sliding, followed by fast co-seismic slip and transient afterslip. Total displacements were about 7, 5.5 and 9.5 micrometers, respectively for the observed pre-slip, co-seismic slip and afterslip. The observed pre-slip lasted about 0.4 seconds. Co-seismic slip was recorded by the 1 kHz displacement recording and a 12-component array of 3-C accelerometers sampled at 100 kHz. The observed afterslip is consistent with analytical models of arrest in a velocity-strengthening region and subsequent stress relaxation. The observed slip vector agrees with the activation of a bedding plane within the phyllite, which is corroborated by relocated seismic events that were observed during the later stages of the injection experiment.

Agricultural managed aquifer recharge (AgMAR) is a proposed management strategy whereby surface water flows are used to intentionally flood croplands with the purpose of recharging underlying aquifers. However, legacy nitrate (NO ₃ ^- ) contamination in agriculturally-intensive regions poses a threat to groundwater resources under AgMAR. To address these concerns, we use a reactive transport modeling framework to better understand the effects of AgMAR management strategies (i.e., by varying the frequency, duration between flooding events, and amount of water) on N leaching to groundwater under different stratigraphic configurations and antecedent moisture conditions. In particular, we examine the potential of denitrification and nitrogen retention in deep vadose zone sediments using variable AgMAR application rates on two-dimensional representations of differently textured soils, soils with discontinuous bands/channels, and soils with preferential flow paths characteristic of typical agricultural field sites. Our results indicate that finer textured sediments, such as silt loams, alone or embedded within high flow regions, are important reducing zones providing conditions needed for denitrification. Simulation results further suggest that applying recharge water all-at-once, rather than in increments, increases denitrification within the vadose zone, but transports higher concentrations of NO ₃ ^- deeper into the profile. This transport into deeper depths can be aggravated by wetter antecedent soil moisture conditions. We conclude that ideal AgMAR management strategies can be designed to enhance denitrification in the subsurface and reduce N leaching to groundwater, while specifically accounting for lithologic heterogeneity, antecedent soil moisture conditions, and depth to the water table.