Bedrock fracture systems facilitate weathering, allowing fresh mineral surfaces to interact with corrosive waters and biota from Earth's surface, while simultaneously promoting drainage of chemically equilibrated fluids. We show that topographic perturbations to regional stress fields explain bedrock fracture distributions, as revealed by seismic velocity and electrical resistivity surveys from three landscapes. The base of the fracture-rich zone mirrors surface topography where the ratio of horizontal compressive tectonic stresses to near-surface gravitational stresses is relatively large, and it parallels the surface topography where the ratio is relatively small. Three-dimensional stress calculations predict these results, suggesting that tectonic stresses interact with topography to influence bedrock disaggregation, groundwater flow, chemical weathering, and the depth of the "critical zone" in which many biogeochemical processes occur.

Geophysics provides a multidimensional suite of investigative methods that are transforming our ability to see into the very fabric of the subsurface environment, and monitor the dynamics of its fluids and the biogeochemical reactions that occur within it. Here we document how geophysical methods have emerged as valuable tools for investigating shallow subsurface processes over the past two decades and offer a vision for future developments relevant to hydrology and also ecosystem science. The field of “hydrogeophysics” arose in the late 1990s, prompted, in part, by the wealth of studies on stochastic subsurface hydrology that argued for better field‐based investigative techniques. These new hydrogeophysical approaches benefited from the emergence of practical and robust data inversion techniques, in many cases with a view to quantify shallow subsurface heterogeneity and the associated dynamics of subsurface fluids. Furthermore, the need for quantitative characterization stimulated a wealth of new investigations into petrophysical relationships that link hydrologically relevant properties to measurable geophysical parameters. Development of time‐lapse approaches provided a new suite of tools for hydrological investigation, enhanced further with the realization that some geophysical properties may be sensitive to biogeochemical transformations in the subsurface environment, thus opening up the new field of “biogeophysics.” Early hydrogeophysical studies often concentrated on relatively small “plot‐scale” experiments. More recently, however, the translation to larger‐scale characterization has been the focus of a number of studies. Geophysical technologies continue to develop, driven, in part, by the increasing need to understand and quantify key processes controlling sustainable water resources and ecosystem services.

[1] Cross-well electrical resistivity tomography (ERT) was used to monitor the migration of a saline tracer in a two-well pumping-injection experiment conducted at the Massachusetts Military Reservation in Cape Cod, Massachusetts. After injecting 2200 mg/L of sodium chloride for 9 hours, ERT data sets were collected from four wells every 6 hours for 20 days. More than 180,000 resistance measurements were collected during the tracer test. Each ERT data set was inverted to produce a sequence of 3-D snapshot maps that track the plume. In addition to the ERT experiment a pumping test and an infiltration test were conducted to estimate horizontal and vertical hydraulic conductivity values. Using modified moment analysis of the electrical conductivity tomograms, the mass, center of mass, and spatial variance of the imaged tracer plume were estimated. Although the tomograms provide valuable insights into field-scale tracer migration behavior and aquifer heterogeneity, standard tomographic inversion and application of Archie's law to convert electrical conductivities to solute concentration results in underestimation of tracer mass. Such underestimation is attributed to (1) reduced measurement sensitivity to electrical conductivity values with distance from the electrodes and (2) spatial smoothing (regularization) from tomographic inversion. The center of mass estimated from the ERT inversions coincided with that given by migration of the tracer plume using 3-D advective-dispersion simulation. The 3-D plumes seen using ERT exhibit greater apparent dispersion than the simulated plumes and greater temporal spreading than observed in field data of concentration breakthrough at the pumping well.Citation: Singha, K., and S. M. Gorelick (2005), Saline tracer visualized with three-dimensional electrical resistivity tomography: Field-scale spatial moment analysis, Water Resour. Res., 41, W05023,

No abstract

Abstract:We want to develop a dialogue between geophysicists and hydrologists interested in synergistically advancing process based watershed research. We identify recent advances in geophysical instrumentation, and provide a vision for the use of electrical and magnetic geophysical instrumentation in watershed scale hydrology. The focus of the paper is to identify instrumentation that could significantly advance this vision for geophysics and hydrology during the next 3-5 years. We acknowledge that this is one of a number of possible ways forward and seek only to offer a relatively narrow and achievable vision. The vision focuses on the measurement of geological structure and identification of flow paths using electrical and magnetic methods. The paper identifies instruments, provides examples of their use, and describes how synergy between measurement and modelling could be achieved. Of specific interest are the airborne systems that can cover large areas and are appropriate for watershed studies. Although airborne geophysics has been around for some time, only in the last few years have systems designed exclusively for hydrological applications begun to emerge. These systems, such as airborne electromagnetic (EM) and transient electromagnetic (TEM), could revolutionize hydrogeological interpretations. Our vision centers on developing nested and cross scale electrical and magnetic measurements that can be used to construct a three-dimensional (3D) electrical or magnetic model of the subsurface in watersheds. The methodological framework assumes a 'top down' approach using airborne methods to identify the large scale, dominant architecture of the subsurface. We recognize that the integration of geophysical measurement methods, and data, into watershed process characterization and modelling can only be achieved through dialogue. Especially, through the development of partnerships between geophysicists and hydrologists, partnerships that explore how the application of geophysics can answer critical hydrological science questions, and conversely provide an understanding of the limitations of geophysical measurements and interpretation.

Details of Earth's shallow subsurface-a key component of the critical zone (CZ)-are largely obscured because making direct observations with sufficient density to capture natural characteristic spatial variability in physical properties is difficult. Yet this inaccessible region of the CZ is fundamental to processes that support ecosystems, society, and the environment. Geophysical methods provide a means for remotely examining CZ form and function over length scales that span centimeters to kilometers. Here we present a review highlighting the application of geophysical methods to CZ science research questions. In particular, we consider the application of geophysical methods to map the geometry of structural features such as regolith thickness, lithological boundaries, permafrost extent, snow thickness, or shallow root zones. Combined with knowledge of structure, we discuss how geophysical observations are used to understand CZ processes. Fluxes between snow, surface water, and groundwater affect weathering, groundwater resources, and chemical and nutrient exports to rivers. The exchange of gas between soil and the atmosphere have been studied using geophysical methods in wetland areas. Indirect geophysical methods are a natural and necessary complement to direct observations obtained by drilling or field mapping. Direct measurements should be used to calibrate geophysical estimates, which can then be used to extrapolate interpretations over larger areas or to monitor changing processes over time. Advances in geophysical instrumentation and computational approaches for integrating different types of data have great potential to fill gaps in our understanding of the shallow subsurface portion of the CZ and should be integrated where possible in future CZ research.

Both vertical and lateral flows of rock and water occur within eroding hills. Specifically, when considered over geological timeframes, rock advects vertically upward under hilltops in landscapes experiencing uplift and erosion. Once rock particles reach the land surface, they move laterally and down the hillslope because of erosion. At much shorter timescales, meteoric water moves vertically downward until it reaches the regional water table and then moves laterally as groundwater flow. Water can also flow laterally in the shallow subsurface as interflow in zones of permeability contrast. Interflow can be perched or can occur during periods of a high regional water table. The depths of these deep and shallow water tables in hills fluctuate over time. The fluctuations drive biogeochemical reactions between water, CO 2 , O 2 , and minerals and these in turn drive fracturing. The depth intervals of water table fluctuation for interflow and groundwater flow are thus reaction fronts characterized by changes in composition, fracture density, porosity, and permeability. The shallow and deep reaction zones can separate over meters in felsic rocks. The zones act like valves that reorient downward unsaturated water flow into lateral saturated flow. The valves also reorient the upward advection of rock into lateral flow through solubilization. In particular, groundwater removes highly soluble, and interflow removes moderately soluble minerals. As rock and water moves through the system, hills may evolve toward a condition where the weathering advance rate, W, approaches the erosion rate, E. If W = E, the slopes of the deep and shallow reaction zones and the hillsides must allow removal of the most soluble, moderately soluble, and least soluble minerals respectively. A permeability architecture thus emerges to partition each evolving hill into dissolved and particulate material fluxes as it approaches steady state.

[1] Hyporheic hydrodynamics are a control on stream ecosystems, yet we lack a thorough understanding of catchment controls on these flow paths, including valley constraint and hydraulic gradients in the valley bottom. We performed four whole-stream solute tracer injections under steady state flow conditions at the H. J. Andrews Experimental Forest (Oregon, United States) and collected electrical resistivity (ER) imaging to directly quantify the 2-D spatial extent of hyporheic exchange through seasonal base flow recession. ER images provide spatially distributed information that is unavailable for stream solute transport modeling studies from monitoring wells alone. The lateral and vertical extent of the hyporheic zone was quantified using both ER images and spatial moment analysis. Results oppose the common conceptual model of hyporheic ''compression'' by increased lateral hydraulic gradients toward the stream. We found that the extent of the hyporheic zone increased with decreasing vertical gradients away from the stream, in contrast to expectations from conceptual models. Increasing hyporheic extent was observed with both increasing and decreasing down-valley (i.e., parallel to the valley gradient) and cross-valley (i.e., from the hillslope to the stream, perpendicular to the valley gradient) hydraulic gradients. We conclude that neither cross-valley nor down-valley hydraulic gradients are sufficient predictors of hyporheic exchange flux nor flow path network extent. Increased knowledge of the controls on hyporheic exchange, the temporal dynamics of exchange flow paths, and their the spatial distribution is the first step toward predicting hyporheic exchange at the scale of individual flow paths. Future studies need to more carefully consider interactions between spatiotemporally dynamic hydraulic gradients and subsurface architecture as controls on hyporheic exchange.