Subsurface Aquifer Research Articles

Dynamic Adaptive and Fully Unstructured Tetrahedral Gridding: Application to CO2 Sequestration with Consideration of Full Fluid Compressibility

Numerical simulations of complex subsurface flow problems and geomechanics will advance enormously by dynamic adaptive gridding. Only a small part of a large domain where sharp changes occur may require fine gridding. In this work, we introduce a methodology to carry out dynamic adaptive gridding for large-scale flow problems. The algorithm allows 2D and 3D unstructured gridding with consideration of full fluid compressibility in single-phase, and two-phase compositional flow. We divide a triangular element into four and a tetrahedron element into eight which creates hanging nodes at each level of refinement. The handing nodes are eliminated by splitting extra elements. As a result, a transition region exists between the fine and coarse grid regions of the domain. We have applied the method to CO2 sequestration in subsurface aquifers. The conditions are selected such that gravity fingers develop from density increase by the dissolution of CO2 in the aqueous phase. The selected examples include large domains where neither systematic studies nor dynamic adaptive gridding have been reported in the past. Results from comparison with uniform gridding reveal a speedup of up to three orders of magnitude in 3D.

Deep learning-based geological parameterization for history matching CO2 plume migration in complex aquifers

History matching is crucial for reliable numerical simulation of geological carbon storage (GCS) in deep subsurface aquifers. This study focuses on inferring highly complex aquifer permeability fields with multi- and intra-facies heterogeneity to improve the characterization of CO2 plume migration. We propose a deep learning (DL)-based parameterization strategy combined with the ensemble smoother with multiple data assimilation (ESMDA) algorithm to formulate an integrated inverse framework. The DL model is employed to parameterize non-Gaussian permeability fields using low-dimensional latent variables in a Gaussian distribution, thereby mitigating the non-Gaussianity issue faced by the ensemble-based ESMDA inverse method and simultaneously alleviating the computational burden of high-dimensional inversion. The efficacy of the integrated DL-ESMDA inverse framework is demonstrated using a 3-D GCS model, where it estimates the non-Gaussian permeability field characterized by multi- and intra-facies heterogeneity. Results show that the DL model is able to represent the highly complex and high-dimensional permeability fields using low-dimensional latent vectors. The DL-ESMDA framework sequentially updates these low-dimensional latent vectors instead of the original high-dimensional permeability field to obtain posterior estimations of the permeability field. The resulting CO2 plume migration closely matches historical measurements, suggesting a significantly improved model reliability after history matching. Additionally, a substantial reduction in uncertainty for future plume migration predictions beyond the history matching period is observed. The proposed framework provides an effective approach for reliable characterization of CO2 plume migration in highly heterogeneous aquifers, enhancing GCS project operation and risk analysis.

Geoelectrical resistivity mapping for sustainable groundwater management in Umuahia South: Insights from vertical electrical sounding

The study of geoelectrical resistivity provides critical insights into subsurface characteristics and aquifer dynamics, particularly in regions with varying geological formations. This research investigates the geoelectrical properties of the aquifer system in Umuahia South, highlighting its significance for sustainable groundwater management. The primary aim of this study is to characterize the geoelectrical layers, assess aquifer thickness and resistivity, and evaluate the implications for groundwater resource management. Vertical Electrical Sounding (VES) was conducted across ten locations to measure the resistivity and thickness of subsurface layers. The data collected were analyzed to identify distinct geoelectrical layers, their respective resistivity values, and thicknesses. The study reveals significant variations in aquifer properties across the region, with VES 1 showcasing the largest thickness of 69.2 m and VES 9 the smallest at 5.2 m. Aquifer thickness decreases from the Benin Formation toward the Ameki and Ogwashi-Asaba Formations. Hydraulic conductivity ranges from 0.207 to 0.954 m/day, indicating varying groundwater flow potential, while transmissivity values vary from 4.965 to 30.441 m²/day, with 60% of aquifers classified as intermediate. Groundwater Potential Index (GWPI) highlights approximately 70% of the area as high potential, underscoring the need for targeted management strategies to optimize water resource extraction and sustainability. The variability in resistivity and thickness across the study area suggests significant subsurface heterogeneity. High resistivity values indicate favorable conditions for groundwater flow, particularly in areas like VES 6. Conversely, lower resistivity in locations such as VES 9 indicates potential challenges for groundwater storage. This study provides vital information on the aquifer system in Umuahia South, emphasizing the need for targeted groundwater management strategies based on the geoelectrical characteristics identified. The research contributes to a deeper understanding of the aquifer dynamics in the region and offers a framework for sustainable groundwater resource management, highlighting the importance of geoelectrical assessments.

Best Management Practices to Mitigate Contamination of Karstic Aquifers from Emergency Fire-Control Runoff

We propose best management practices to increase karstic aquifer drinking water supplies’ resilience to potential hazardous material impacts where first responders and public safety officials suspect emergency firefighting runoff has entered the subsurface in the aquifer recharge zone. Karstic aquifers are unique because of their direct openings to the land’s surface, which allows contaminants in runoff or other surface waters to rapidly enter the subsurface and impact aquifer water quality. In the United States, 20% of the land surface is karst, and about a third of the groundwater used for drinking comes from karstic aquifers. Proposed best management practices emphasize on-site, real-time evaluation of the transport and fate of HAZMAT that may enter the subsurface and focus on water quality sampling, runoff and groundwater flow modeling, nontoxic dye tracing, and related studies for use in planning before and during an event and after emergency response has ended. We recommend best management practices for evidence-based screening of high-risk sites to facilitate placement of hazardous material sampling devices and emergency responder deployment. The best management practices, tools, curricula, and training described combine with earlier work by the authors to provide a comprehensive menu of actions to (1) help first responders prevent or limit flushing of hazardous material into a karstic aquifer; (2) help emergency management officials understand the consequences of contamination and issue data-driven geographically focused health and safety risk communications; and (3) help provide health and safety officials with relevant science- and data-based information that can help mitigate human and environmental health risks.

Source localization in subsurface aquifers based on conservation data by learning a Gaussian kernel

Source localization in subsurface aquifers based on conservation data by learning a Gaussian kernel

Evaluation of subsurface geothermal groundwater aquifers at Southern Kirthar Range Sindh province Pakistan through the application of vertical electrical sounding

Pakistan has significant geothermal potential but due to poor energy management and low utilization of rich geothermal energy, the country is ranked among the energy deficient nations in the world. In this study, integrated geophysical techniques were used to assess the geothermal energy potential of various hot springs in the Sindh province. The study revealed that Naing Shareef and Gaji Shah hot springs have remarkable potential for geothermal energy having pool water temperature of 43 °C and 45 °C, respectively. Geophysical resistivity survey was carried out at both Naing Shareef and Gaji Shah hot spring to determine depth, thickness and quality of thermal water. At Naing Shareef hot spring, only a single thermal aquifer zone has been pinpointed within the shale composition of the Oligocene Nari Formation, and it is on average 55.07 m depth and 29.8 m thick. At an average depth of 100.66 m and an aquifer thickness of 77.03 m, the Gaji Shah hot springs exclusively possess a thermal aquifer region established within the sandstone composition of the Oligocene Nari Formation.

Investigation of Aquifer Zone using Central Loop Time-domain Electromagnetic in Institut Teknologi Bandung

Abstract Time-domain electromagnetic (TDEM) has been widely used to investigate the subsurface aquifer. In this research, the TDEM method with a central loop configuration was employed to define the aquifer zone in the southern area of the Institut Teknologi Bandung (ITB) campus in Bandung, Indonesia. Moreover, TDEM is effective for mapping aquifers due to its sensitivity to subsurface conductivity. The acquisition consists of five sounding points spaced 20 meters apart, along with a W-E profile in the southern area of ITB. The data were interpreted using 1D inversion software (1X1D, Interpex, USA), utilizing weights from a smoothing constraint scheme (Occam) to obtain the resistivity distribution of the subsurface structure. The interpretation results indicate three layers from top to bottom: (a) the topsoil was interpreted as a sand layer with resistivity of 100 to 350 Ωm; (b) beneath it lies a shale layer with resistivity ranging from 5 to 100 Ωm that has good porosity to allow water from the surface to flow into the ground; and (c) further down, there is a clay layer with low resistivity of 1 to 5 Ωm which act as impermeable layer.

Hydrogeological Resolution Using Composite Geophysical Methods: A Case Study of Inland Salinity in Part of Agra, Uttar Pradesh India

ABSTRACT This manuscript is an attempt at utilizing a combination of varied Geophysical methods (n = 4) – Vertical Electrical Sounding (VES), Electrical Resistivity Tomography (ERT), Time Domain Electromagnetic (TEM), and Borehole logging for finer resolution of potential fresh groundwater and saline groundwater aquifer zones in marginal alluvial plains of Agra district, Uttar Pradesh state. A conceptual model of sub-surface aquifer disposition displaying zones of fresh and saline groundwater zones was prepared using an interpreted and processed multi-parameter dataset. A positive correlation was found between TEM (44.6 m), Well log data (43.0 m), VES data (47.2), and ERT data (45.0 m) demarcating the boundary up to which fresh groundwater potential lies from the ground surface. Conventional methods rely on resistivity values to identify groundwater potential zones, and a rough estimation of salinity can be carried out. Instead of the conventional VES method, the combination of aforesaid methods results in finer resolution and easier demarcation of fresh groundwater zones from those of saline zones. This method can be successfully reproduced to demarcate the extent of saline water ingress in coastal areas and finding freshwater lenses within inland saline aquifers.

Analysis of facies proportions as a tool to quantify reservoir heterogeneity

The geometry, orientation and stacking patterns of sedimentary geobodies play a key control on fluid migration within subsurface reservoirs, aquifers, and repositories. Outcrops and more recently virtual outcrops are commonly used as analogues to characterise geobody architecture and stacking. Traditionally, both reservoirs and their analogues are described with respect to the bulk proportion of different facies and architectural elements within the system. These are combined with geometric measurements of the different geobodies to describe heterogeneity.This study presents a new methodology for quantifying the range of facies proportions within architectural element panels taken from outcrop analogues. This method measures the facies proportions within a series of vertical profiles and considers the statistical distribution, rather than the simple mean across the entire panel. The basic premise being that within a layer cake system, the spread in facies proportions will be very narrow whereas in a more “jigsaw” type system the spread will be much wider.To test this hypothesis, a series of nine virtual outcrops from a range of clastic depositional systems including fluvial, shallow marine and deep marine settings, have been analysed. Depositional facies (architectural elements) were interpreted in LIME and a series of depositional strike-orientated orthorectified panels were extracted. The panels were then analysed, and the extracted facies proportions were translated to “net sand” distribution. Results were then compared with conventional analysis of reservoir heterogeneity.Different depositional systems show distinct patterns. When arranged in depositional dip order, there is a systematic decrease in sand distribution spread down dip with significant resets at the boundaries between the major gross depositional environments (continental, shallow marine/shelf, deep marine).

Application of 2-dimensional tomography for subsurface lithology imaging and aquifer determination in Etim Ekpo, Southern Nigeria.

Application of 2-dimensional tomography for subsurface lithology imaging and aquifer determination in Etim Ekpo, Southern Nigeria.

Investigation on radon levels in soil and water associated with Mount Tampomas geothermal activity in Indonesia

ObjectiveTo investigate radon emanation in geothermal manifestations around Mount Tampomas, West Java, Indonesia, and assess radon concentrations in soil and water samples. MethodsRadon measurements were conducted using the Durridge Rad7 instrument, supplemented with a soil gas probe for in-situ soil radon measurements at a depth of 80 ​cm. In-situ water radon measurements were performed using the Rad Aqua instrument, while radon measurements for hot water samples were conducted separately. Radon measurements for hot water samples were corrected for decay using a radon decay correction factor. ResultsThe analysis of radon measurements revealed a wide range of concentrations in soil and water samples. Soil radon concentrations ranged from 15 Bq/m³ to 4,660 Bq/m³, with localized hotspots exhibiting exceptionally high concentrations. Water radon measurements showed elevated levels, ranging from 0.2 Bq/L to 13.4 Bq/L in-situ, particularly in hot springs. In collected water samples, radon concentrations ranged from 1 Bq/L to 6 Bq/L. These combined results highlight significant variability in radon levels across different water sources influenced by geothermal activity. ConclusionsThese findings indicate active emanation processes influenced by geological factors and underscore the role of subsurface geology and aquifer characteristics in radon transport mechanisms. The presence of localized radon hotspots suggests the need for comprehensive monitoring and proactive management strategies to mitigate environmental and public health risks associated with radon exposure.

Application of modified drastic model for oil spills pollution affecting water quality system in part of Niger delta region of Nigeria

Extensive environmental pollution has affected the water resource, agricultural lands, economy and health of the inhabitants of Ahoada communities of the Niger Delta region, Nigeria. The main sources of pollution are oil spills emanating from pipeline interdiction and production-related spillages. Oil spills in the Niger Delta adversely affects the physical, chemical and biological properties of soil, which subsequently results in food shortages due to the reduction of nutrient contents of soils. Therefore, it is crucial to understand the patterns and the types of contaminants for appropriate remediation of the polluted areas. This study applied a novel DRASTIC model approach, incorporating a data-driven method Weight of Evidence (WoE) to examine subsurface aquifer vulnerability to oil spills contamination and present same as in the maps. Geographic Information System (GIS) in combination with DRASTIC and a modified DRASTIC model called DRAPSTIC were used in assessing the extent and vulnerability of the aquifers of the area to contamination. The maps of the DRAPSTIC show similarity in terms of the distribution and classification of the vulnerability index. Very high and high vulnerability classes covered larger areas in Ahoada east under DRAPSTIC compared to DRASTIC. However, under the same DRAPSTIC, larger areas were covered by low vulnerability class in Ahoada West.

Reporting nanoparticle tracers: Validation of performance in flow-through experiments simulating reservoir conditions

Characterization and monitoring of reservoir properties and conditions are key problems in the exploitation of subsurface aquifers and reservoirs, with tracer tests being an important tool that provides valuable insights into groundwater dynamics. The Reporting Nanoparticle Tracers (RNTs) approach was recently presented with the aim of expanding the scope of measureable parameters and the potential benefits of tracer tests in comparison to tests performed with traditional tracers. However, successful implementation of the concept depends on the ability of the nanoparticle tracers to exhibit a stable dispersion that facilitates sufficiently high recovery rates to allow for meaningful analysis. As a step forward toward the implementation of the approach, we used flow‑through studies with packed-bed quartz sand columns to compare the transport and retention properties of the RNTs with those of prevalent conservative molecular tracers and to show the feasibility of temperature detection. In the main experiments, the RNTs showed recovery rates around 80%, outperformed only by uranine (95 %) and surpassing eosin and sulforhodamine G (15 % and 50 % respectively). Thermally-activated RNTs could also be clearly differentiated from non-activated ones by a signal ratio difference of 50 %. These results validate the feasibility of application, especially considering the modularity of the nanoparticles, which enables adjustment of the RNTs to the hydrochemical conditions prevalent in the tested aquifer.

The dominance and growth of shallow groundwater resources in continuous permafrost environments

Water is a limited resource in Arctic watersheds with continuous permafrost because freezing conditions in winter and the impermeability of permafrost limit storage and connectivity between surface water and deep groundwater. However, groundwater can still be an important source of surface water in such settings, feeding springs and large aufeis fields that are abundant in cold regions and generating runoff when precipitation is rare. Whether groundwater is sourced from suprapermafrost taliks or deeper regional aquifers will impact water availability as the Arctic continues to warm and thaw. Previous research is ambiguous about the role of deep groundwater, leading to uncertainty regarding Arctic water availability and changing water resources. We analyzed chemistry and residence times of spring, stream, and river waters in the continuous permafrost zone of Alaska, spanning the mountains to the coastal plain. Water chemistry and age tracers show that surface waters are predominately sourced from recent precipitation and have short (<50 y) subsurface residence times. Remote sensing indicates trends in the areal extent of aufeis over the last 37 y, and correlations between aufeis extent and previous year summer temperature. Together, these data indicate that surface waters in continuous permafrost regions may be impacted by short flow paths and shallow suprapermafrost aquifers that are highly sensitive to climatic and hydrologic change over annual timescales. Despite the lack of connection to regional aquifers, continued warming and permafrost thaw may promote deepening of the shallow subsurface aquifers and creation of shallow taliks, providing some resilience to Arctic freshwater ecosystems.

Dynamics of hydrogen storage in subsurface saline aquifers: A computational and experimental pore-scale displacement study

The catastrophic effects of climate change can be mitigated by transitioning to renewable energy sources such as wind, solar, and hydropower while utilizing hydrogen (H2) as an energy carrier. As renewable fuel sources are intermittent and location-specific, large-scale, long-term storage options for H2 must be explored with high necessity. Subsurface hydrogen storage in saline aquifers provides a vital scope to store H2 gas in available pore voids with minimal environmental risk. In this work, a highly heterogeneous porous micromodel was used to study the immiscible two-phase dynamics at high-temperature, high-pressure saline aquifer conditions using computational fluid dynamics (CFD) simulation based on volume of fluid (VOF) method. A set of microfluidic experiments were carried out under atmospheric conditions to validate the numerical model. The VOF model was observed to show good agreement with microfluidic experiments. An unstable displacement pattern with viscous fingering was observed during various flow conditions that captured high pressure (15 and 25 MPa) and temperature conditions (323 and 343 K). It was observed that fingering displacement of brine limits the storage spaces for hydrogen, which consequently follows the high permeable path. As a result, only 0.272 fraction of brine was invaded by the H2 at a temperature of 323 K and pressure of 15 MPa. A comparison between H2-brine and nitrogen (N2)-brine flow dynamics revealed a 46% less storage capacity of porous media for H2. Formations bearing high temperature (343 K) showed an 11.27% increase in H2 storage capacity, while pressure increase had negligible effect. At low capillary number (Ca), more snap-off effects resulted in higher trapping of H2 gas. This study aims to provide meaningful insights into the complexities of flow dynamics and displacement patterns that are crucial for optimizing hydrogen storage efficiency in saline aquifers and for potential ‘white’ hydrogen production.

Unveiling subsurface heterogeneity in porous aquifers: Insights from hydrogeophysics and derivative analysis

Groundwater is a crucial resource, particularly in urban areas where the need for alternative water sources is rising. Securing groundwater resources is paramount, requiring a concentrated effort to assess these resources and mitigate the increasing water shortages in urban regions. Despite its significance, the application of hydrogeophysics and derivative analysis in understanding aquifer dynamics is often overlooked and underused. Consequently, this study argues that neglecting the integration of hydrogeophysics dataset and derivative analysis in aquifer characterization leads to developing models that lack solution-oriented approaches, hindering effective groundwater management. This study aimed to enhance understanding of aquifer dynamics for solution-based modeling while emphasizing the importance of integrating hydrogeophysics dataset and derivative analysis to highlight aquifer system heterogeneities. The electrical resistivity tomography and derivative analysis of pumping test were used in this study to image the subsurface and amplify aquifer flow regimes respectively. The results of the electrical resistivity survey revealed a distinct layer of fine to medium-grain sand at depths of approximately 60 m in certain areas intercalated with medium-to-coarse-grain sand and thin layers of peat. Derivative analysis plots indicated that the predominant flow regime is linear and bilinear, with evidence of fracture dewatering during pumping cycles suggesting an unconfined aquifer. Additionally, aquifer heterogeneity and a no-flow boundary were evident during the pumping cycle. This study underscores the efficacy of combining hydrogeophysics and derivative analysis of pumping tests as a robust approach for imaging and validating subsurface conditions. The implications of these findings extend beyond the specific case study, offering valuable insights for groundwater utilization, monitoring, and management on a broader scale. By elucidating effective methodologies, this research enhances groundwater security and meets the escalating demand for sustainable water sources in urban areas. In conclusion, electrical resistivity tomography is an effective tool for characterizing subsurface aquifer systems, while derivative analysis of pumping tests is crucial for highlighting aquifer system heterogeneities. These methods offer a more practical interpretation compared to traditional drawdown versus time curve approaches, making the results invaluable for groundwater usage monitoring and management.

Quantifying genome-specific carbon fixation in a 750-meter deep subsurface hydrothermal microbial community.

Dissolved inorganic carbon has been hypothesized to stimulate microbial chemoautotrophic activity as a biological sink in the carbon cycle of deep subsurface environments. Here, we tested this hypothesis using quantitative DNA stable isotope probing of metagenome-assembled genomes (MAGs) at multiple 13C-labeled bicarbonate concentrations in hydrothermal fluids from a 750-m deep subsurface aquifer in the Biga Peninsula (Turkey). The diversity of microbial populations assimilating 13C-labeled bicarbonate was significantly different at higher bicarbonate concentrations, and could be linked to four separate carbon-fixation pathways encoded within 13C-labeled MAGs. Microbial populations encoding the Calvin-Benson-Bassham cycle had the highest contribution to carbon fixation across all bicarbonate concentrations tested, spanning 1-10mM. However, out of all the active carbon-fixation pathways detected, MAGs affiliated with the phylum Aquificae encoding the reverse tricarboxylic acid (rTCA) pathway were the only microbial populations that exhibited an increased 13C-bicarbonate assimilation under increasing bicarbonate concentrations. Our study provides the first experimental data supporting predictions that increased bicarbonate concentrations may promote chemoautotrophy via the rTCA cycle and its biological sink for deep subsurface inorganic carbon.

Unveiling the inhibitory mechanisms of chromium exposure on microbial reductive dechlorination: Kinetics and microbial responses

Chromium and organochlorine solvents, particularly trichloroethene (TCE), are pervasive co-existing contaminants in subsurface aquifers due to their extensive industrial use and improper disposal practices. In this study, we investigated the microbial dechlorination kinetics under different TCE-Cr(Ⅲ/VI) composite pollution conditions and elucidated microbial response mechanisms based on community shift patterns and metagenomic analysis. Our results revealed that the reductive dechlorinating consortium had high resistance to Cr(III) but extreme sensitivity to Cr(VI) disturbance, resulting in a persistent inhibitory effect on subsequent dechlorination. Interestingly, the vinyl chloride-respiring organohalide-respiring bacteria (OHRB) was notably more susceptible to Cr(III/VI) exposure than the trichloroethene-respiring one, possibly due to inferior competition for growth substrates, such as electron donors. In terms of synergistic non-OHRB populations, Cr(III/VI) exposure had limited impacts on lactate fermentation but significantly interfered with H2-producing acetogenesis, leading to inhibited microbial dechlorination due to electron donor deficiencies. However, this inhibition can be effectively mitigated by the amendment of exogenous H2 supply. Furthermore, being the predominant OHRB, Dehalococcoides have inherent Cr(VI) resistance defects and collaborate with synergistic non-OHRB populations to achieve concurrent bio-detoxication of Cr(VI) and TCE. Our findings expand the understanding of the response patterns of different functional populations towards Cr(III/VI) stress, and provide valuable insights for the development of in situ bioremediation strategies for sites co-contaminated with chloroethene and chromium.

Combined effects of aquifer heterogeneity and subsurface dam on nitrate contamination in coastal aquifers

Subsurface dams are effective for seawater intrusion mitigation, yet they can cause upstream nitrate accumulation. This research examines the interplay between subsurface dam construction and aquifer layering on nitrate pollution in coastal settings, employing numerical models to simulate density-driven flow and reactive transport. The study reveals that while subsurface dams are adept at curbing seawater intrusion, they inadvertently broaden the nitrate accumulation zone, especially when a low-K layer is present. Heterogeneous aquifers see more pronounced nitrate accumulation from subsurface dams. This effect is pronounced as it influences dissolved organic carbon dynamics, with a notable retreat inland correlating with the expansion of the nitrate pollution plume. A critical finding is that controlling seawater intrusion via dam height adjustment within the Effective Damming Region effectively reduces nitrate levels and bolsters freshwater output. However, exceeding the critical threshold—where the dam surpasses the low-K layer's bottom—results in a substantial shift in nitrate concentration, underscoring the need for precise dam height calibration to avoid aggravating nitrate pollution. This study's innovative contribution lies in its quantification of the nuanced effects of subsurface dams in stratified aquifers, providing an empirical basis for dam design that considers the layered complexities of coastal aquifers. The insights offer a valuable framework for managing nitrate contamination, thus informing sustainable coastal groundwater management and protection strategies.

The social connectivity of subsurface flows: Towards a better integration of the vertical dimension in socio‐hydrosystem studies

AbstractThis contribution points out that while the importance of hydrologic, geomorphic, ecological, temporal, and socio‐cultural connectivity in the functioning of hydrosystems has been acknowledged in three dimensions (longitudinal, lateral, and vertical), vertical connectivity has often been overlooked. Drawing on a multidisciplinary literature review, the authors aim to highlight the socio‐cultural connectivity of subsurface flows and aquifers as a crucial factor for socio‐hydrosystem understanding and management. The piece builds on emergent literature which underscores how groundwater, shallow groundwater, and the hyporheic zone are coproduced by nature and society through time. Furthermore, the review explores how verticality has become an important heuristic dimension at the intersection of the environmental and social sciences, and there has been a particular focus on the hyporheic zone to look at how notions of interstitiality and (in)visibility can be better integrated with socio‐hydrosystem science and management. Finally, the paper calls for further research to integrate the vertical dimension of hydrosystems into more comprehensive socio‐hydrological frameworks, which remain, at times, empirically and theoretically weak on questions of social power, even if they do incorporate aspects of political systems. Especially as societies' relationships to groundwater may be at the heart of climate change adaptation strategies, greater consideration of the social connectivity to subflows is a necessary direction for sustainable water resource management and scholarship.This article is categorized under: Human Water &gt; Water Governance Science of Water &gt; Hydrological Processes Science of Water &gt; Water and Environmental Change