Gravity-driven infiltration of liquid water into unsaturated porous media can be a spatially heterogeneous process due to the gravity fingering instability. When such infiltration occurs in a subfreezing porous medium, liquid water can readily freeze, leading to both the removal of liquid water available for transport and a reduction in local permeability. As a result of the coupling between gravity fingering and freezing, macroscopic frozen structures can form that record the shape and history of the wetting front. These structures have been observed in the field in terrestrial snowpack and glacial firn layers and are believed to have profound impacts on how liquid water and its accompanying thermal content distribute during infiltration. However, a more detailed physics-based understanding of freezing infiltration has been missing. In this work, we use a thermodynamic nonequilibrium infiltration model to investigate the emergence of refrozen structures during water infiltration into an initially homogeneous and subfreezing porous medium. From scaling analysis, we recover the relevant nondimensional groups that govern the physics of the freezing infiltration process. We identify two key mechanisms caused by freezing that reduce the effective infiltration rate, calculated as the maximum depth of infiltration per elapsed time. In the first mechanism, the effective infiltrate rate decreases because a portion of the liquid water is consumed due to freezing, and such an effect can be well quantified by the freezing Damk & ouml;hler number. For the second mechanism, we report on a new phenomenon termed secondary fingering, where new flow paths are established in between the primary infiltration channels. We find that secondary fingering reduces the degree of flow channelization and thus weakens the effective rate of infiltration via flow field homogenization. Finally, we identify a regime under high Damk & ouml;hler number in which freezing-induced permeability reduction completely clogs flow.

We perform a direct comparison between experiment and simulation of density-driven convective mixing in three-dimensional (3D) porous media. We find excellent agreement between the experiment and the model in terms of both the convection fingering pattern and the average rate of fluid mixing. In particular, the experiment exhibits dynamic self-organization of columnar plumes into a reticular pattern, which, until now, had only been observed in 3D simulations. We also report good quantitative agreement between the experiment and simulation in the evolution of the state of mixing by comparing, over time, (i) the average concentration at depth, (ii) the variance of the concentration field, (iii) the scalar dissipation rate, and (iv) the dissolution flux. We derive a relation between the scalar dissipation rate and the dissolution flux in a closed system, and we show that the flux in a 3D system is approximately similar to 30% higher than in a 2D system, confirming previous numerical estimates.

The microstructure of snow determines its fundamental properties such as mechanical strength, reflectivity, or thermo-hydraulic properties. Snow undergoes continuous microstructural changes due to local gradients in temperature, humidity, or curvature, in a process known as snow metamorphism. In this work, we focus on wet snow metamorphism, which occurs when the temperature is close to the melting point and involves phase transitions among liquid water, water vapor, and solid ice. We propose a pore-scale phase-field model that simultaneously captures the three relevant phase change phenomena: sublimation (deposition), evaporation (condensation), and melting (solidification). The phase-field formulation allows one to track the temperature evolution among the three phases and the water vapor concentration in the air. Our three-phase model recovers the corresponding two-phase transition model when one phase is not present in the system. 2D simulations of the model unveil the impact of humidity and temperature on the dynamics of wet snow metamorphism at the pore scale. We also explore the role of liquid melt content in controlling the dynamics of snow metamorphism in contrast to the dry regime before percolation onsets. The model can be readily extended to incorporate two-phase flow and may be the basis for investigating other problems involving water phase transitions in a vapor-solid-liquid system, such as airplane icing or thermal spray coating.

Vadose zone soil moisture is often considered a pivotal intermediary water reservoir between surface and groundwater in semi-arid regions. Understanding its dynamics in response to changes in meteorologic forcing patterns is essential to enhance the climate resiliency of our ecological and agricultural system. However, the inability to observe high-resolution vadose zone soil moisture dynamics over large spatiotemporal scales hinders quantitative characterization. Here, utilizing pre-existing fiber-optic cables as seismic sensors, we demonstrate a fiber-optic seismic sensing principle to robustly capture vadose zone soil moisture dynamics. Our observations in Ridgecrest, California reveal sub-seasonal precipitation replenishments and a prolonged drought in the vadose zone, consistent with a zero-dimensional hydrological model. Our results suggest a significant water loss of 0.25 m/year through evapotranspiration at our field side, validated by nearby eddy-covariance based measurements. Yet, detailed discrepancies between our observations and modeling highlight the necessity for complementary in-situ validations. Given the escalated regional drought risk under climate change, our findings underscore the promise of fiber-optic seismic sensing to facilitate water resource management in semi-arid regions. A fiber-optic seismic sensing principle provides a spatiotemporal view of vadose zone soil moisture dynamics, offering a scalable solution for enhanced long-term, large-scale water resource management in face of climate variability.

Particulate precipitation, deposition, and accumulation, including the formation of salt and mineral crystals, frequently occur in a wide range of subsurface applications involving multiphase flow through porous media. Consequently, there has been a considerable emphasis on researching and understanding these phenomena. However, modeling particle dynamics in flows through porous media with low Reynolds numbers has always been a challenging problem as it requires resolving fluid flow around the moving solid particles, the solid–solid contact mechanics, and the solid–fluid coupling. The discrete element method coupled with fluid solvers has been widely used to study particle-laden flow. Most fluid-solid numerical schemes involve solving the full or generalized Navier–Stokes equations, which often yields relatively accurate fluid-solid interactions at the cost of computation time and particle shape limitations. In this paper, we present a novel method to study mono-layered particle-laden flow by coupling the level set discrete element method (LS-DEM) with Hele-Shaw flow model. Utilizing the Hele-Shaw flow model allows us to simplify flow computation, while incorporating LS-DEM enables the simulation of arbitrarily shaped particles. Cases of mono-layered particle flow through a simplified micromodel geometry are studied and validated against published experimental results. Moreover, the effects of particle friction and shape on clogging statistics are investigated.

The transport of meltwater through porous snow is a fundamental process in hydrology that remains poorly understood but essential for more robust predictions of how the cryosphere will respond under climate change. Here, we propose a continuum model that resolves the nonlinear coupling of preferential melt flow and the nonequilibrium thermodynamics of ice‐melt phase change at the Darcy scale. We assume that the commonly observed unstable melt infiltration is due to the gravity fingering instability and capture it using the modified Richards equation, which is extended with a higher‐order term in saturation. Our model accounts for changes in porosity and the thermal budget of the snowpack caused by melt refreezing at the continuum scale, based on a mechanistic estimate of the ice‐water phase change kinetics formulated at the pore scale. We validate the model in 1D against field data and laboratory experiments of infiltration in snow and find generally good agreement. Compared to existing theory of stable melt infiltration, our 2D simulation results show that preferential infiltration delivers melt faster to deeper depths, and as a result, changes in porosity and temperature can occur at deeper parts of the snow. The simulations also capture the formation of vertical low porosity annulus known as ice pipes, which have been observed in the field but lack mechanistic understanding to date. Our results demonstrate how melt refreezing and unstable infiltration reshape the porosity structure of snow and impacts thermal and mass transport in highly nonlinear ways that are not captured by simpler models.

In this study, we experimentally investigate fluid–fluid displacement in a Hele–Shaw cell where the two fluids react, upon mixing, to form solid precipitates. Under the conditions of our experiments, we observe that precipitation reaction along the moving fluid–fluid interface generates solids in the form of mineral particle suspensions. We find that both electrostatic and hydrodynamic forces control the extent of particle–particle agglomeration during the suspension flow. Such particle suspension decreases the overall mobility of the multiphase mixture, thus altering the overall displacement. Although the injected fluids are viscously stable, the precipitation band that forms between the fluids becomes unstable to form finger-like flow channels compartmentalized by solid-deposited walls and clusters. We show that the emergence, growth and decay of the fingering pattern are strongly influenced by the injection rate and the initial fluid chemical concentrations. In addition, we show that precipitation-induced fingering has a strong feedback on fluid–fluid mixing and the subsequent precipitation rate. Lastly, we find that, counter to intuition, a higher injection rate results in a larger amount of precipitates that are securely deposited in the Hele–Shaw cell.

This data deposit contains the data necessary to reproduce results presented in the paper "Hydrate Formation on Marine Seep Bubbles and the Implications for Water Column Methane Dissolution", to appear in JGR Oceans. The files "Fu_AGUSupplementary_S11.xlsx" and "Fu_AGUSupplementary_S12.xlsx" are referenced in the supplement of the original paper. The file "Fu_AGUSupplementary_additionalDATA.xlsx" contains all the data and inputs used to create Figure 14, 15 and 16 of the original manuscript.

Abstract Fluid–fluid displacement in porous media has been viewed through the lens of Lenormand's phase diagram since the late 1980s. This diagram suggests that the character of the flow is controlled by two dimensionless parameters: the capillary number and the viscosity ratio. It is by now well known, however, that the wettability of the system plays a key role in determining the pore-scale displacement mechanisms and macroscopic invasion patterns. Here, we endow Lenormand's diagram with the impact of wettability using dynamic and quasi-static pore-network models. By using the fractal dimension and the ratio of characteristic viscous and capillary pressures we delineate the five principal displacement regimes within the extended phase diagram: stable displacement, viscous fingering, invasion percolation, cooperative pore filling and corner flow. We discuss the results in the context of pattern formation, displacement-front dynamics, pore-scale disorder and displacement efficiency.

Modeling water flow and solute transport in the vadose zone is essential to understanding the fate of soil pollutants and their travel times towards groundwater bodies. It also helps design better irrigation strategies to control solute concentrations and fluxes in semiarid and arid regions. Heterogeneity, soil texture and wetting front instabilities determine the flow patterns and solute transport mechanisms in dry soils. When water is already present in the soil, the flow of an infiltration pulse depends on the spatial distribution of soil water and on its mobility. We present numerical simulations of passive solute transport during unstable infiltration of water into sandy soils that are prone to wetting front instability. We study the impact of the initial soil state, in terms of spatial distribution of water content, on the infiltration of a solute-rich water pulse. We generate random fields of initial moisture content with spatial structure, through multigaussian fields with prescribed correlation lengths. We characterize the patterns of water flow and solute transport, as well as the mass fluxes through the soil column. Our results indicate a strong interplay between preferential flow and channeling due to fingering and the spatial distribution of soil water at the beginning of infiltration. Fingering and initial water saturation fields have a strong effect on solute diffusion and dilution into the ambient water during infiltration, suggesting an effective separation between mobile and inmobile transport domains that are controlled by the preferential flow paths due to fingering.

Earth and Space Science Open Archive PosterOpen AccessYou are viewing the latest version by default [v1]LIP flows may not have been as thick as they appearAuthorsJonasKatonaiDXiaojingFuTusharMittaliDMichaelMangaStephenSelfiDSee all authors Jonas KatonaiDCorresponding Author• Submitting AuthorYale UniversityiDhttps://orcid.org/0000-0002-2127-8912view email addressThe email was not providedcopy email addressXiaojing FuUniversity of California Berkeleyview email addressThe email was not providedcopy email addressTushar MittaliDUniversity of California BerkeleyiDhttps://orcid.org/0000-0002-8026-0018view email addressThe email was not providedcopy email addressMichael MangaUniversity of California Berkeleyview email addressThe email was not providedcopy email addressStephen SelfiDUniversity of California BerkeleyiDhttps://orcid.org/0000-0002-8081-3532view email addressThe email was not providedcopy email address

Widespread seafloor methane venting has been reported in many regions of the world oceans in the past decade. Identifying and quantifying where and how much methane is being released into the ocean remains a major challenge and a critical gap in assessing the global carbon budget and predicting future climate [C. Ruppel, J. D. Kessler. Rev. Geophys. 55, 126-168 (2017)]. Methane hydrate ([Formula: see text]) is an ice-like solid that forms from methane-water mixture under elevated-pressure and low-temperature conditions typical of the deep marine settings (>600-m depth), often referred to as the hydrate stability zone (HSZ). Wide-ranging field evidence indicates that methane seepage often coexists with hydrate-bearing sediments within the HSZ, suggesting that hydrate formation may play an important role during the gas-migration process. At a depth that is too shallow for hydrate formation, existing theories suggest that gas migration occurs via capillary invasion and/or initiation and propagation of fractures (Fig. 1). Within the HSZ, however, a theoretical mechanism that addresses the way in which hydrate formation participates in the gas-percolation process is missing. Here, we study, experimentally and computationally, the mechanics of gas percolation under hydrate-forming conditions. We uncover a phenomenon-crustal fingering-and demonstrate how it may control methane-gas migration in ocean sediments within the HSZ.

Water infiltration and unsaturated flow through heterogeneous soil control the distribution of soil moisture in the vadose zone and the dynamics of groundwater recharge, providing the link between climate, biogeochemical soil processes and vegetation dynamics. Infiltration into dry soil is hydrodynamically unstable, leading to preferential flow through narrow wet regions (fingers). In this paper we use numerical simulation to study the interplay between fingering instabilities and soil heterogeneity during water infiltration. We consider soil with heterogeneous intrinsic permeability. Permeabilities are random, with point Gaussian statistics, and vary smoothly in space due to spatial correlation. The key research question is whether the presence of moderate or strong heterogeneity overwhelms the fingering instability, recovering the simple stable displacement patterns predicted by most simplified model of infiltration currently used in hydrological models from the Darcy to the basin scales. We perform detailed simulations of constant-rate infiltration into soils with isotropic and anisotropic intrinsic permeability fields. Our results demonstrate that soil heterogeneity does not suppress fingering instabilities, but it rather enhances its effect of preferential flow and channeling. Fingering patterns strongly depend on soil structure, in particular the correlation length and anisotropy of the permeability field. While the finger size and flow dynamics are only slightly controlled by correlation length in isotropic fields, layering leads to significant finger meandering and bulging, changing arrival times and wetting efficiencies. Fingering and soil heterogeneity need to be considered when upscaling the constitutive relationships of multiphase flow in porous media (relative permeability and water retention curve) from the finger to field and basin scales. While relative permeabilities remain unchanged upon upscaling for stable displacements, the inefficient wetting due to fingering leads to relative permeabilities at the field scale that are significantly different from those at the Darcy scale. These effective relative permeability functions also depend, although less strongly, on heterogeneity and soil structure.

We develop a novel `moving capacitor' dynamic network model to simulate immiscible fluid-fluid displacement in porous media. Traditional network models approximate the pore geometry as a network of fixed resistors, directly analogous to an electrical circuit. Our model additionally captures the motion of individual fluid-fluid interfaces through the pore geometry by completing this analogy, representing interfaces as a set of moving capacitors. By incorporating pore-scale invasion events, the model reproduces, for the first time, both the displacement pattern and the injection pressure signal under a wide range of capillary numbers and substrate wettabilities. We show that at high capillary numbers the invading patterns advance symmetrically through viscous fingers. In contrast, at low capillary numbers the flow is governed by the wettability-dependent fluid-fluid interactions with the pore structure. The signature of the transition between the two regimes manifests itself in the fluctuations of the injection pressure signal.

Methane hydrate occurs naturally under pressure and temperature conditions that are not straightforward to replicate experimentally. Xenon has emerged as an attractive laboratory alternative to methane for studying hydrate formation and dissociation in multiphase systems, given that it forms hydrates under milder conditions. However, building reliable analogies between the two hydrates requires systematic comparisons, which are currently lacking. We address this gap by developing a theoretical and computational model of gas hydrates under equilibrium and nonequilibrium conditions. We first compare equilibrium phase behaviors of the Xe·H 2 O and CH 4 ·H 2 O systems by calculating their isobaric phase diagram, and then study the nonequilibrium kinetics of interfacial hydrate growth using a phase field model. Our results show that Xe·H 2 O is a good experimental analog to CH 4 ·H 2 O, but there are key differences to consider. In particular, the aqueous solubility of xenon is altered by the presence of hydrate, similar to what is observed for methane; but xenon is consistently less soluble than methane. Xenon hydrate has a wider nonstoichiometry region, which could lead to a thicker hydrate layer at the gas‐liquid interface when grown under similar kinetic forcing conditions. For both systems, our numerical calculations reveal that hydrate nonstoichiometry coupled with hydrate formation dynamics leads to a compositional gradient across the hydrate layer, where the stoichiometric ratio increases from the gas‐facing side to the liquid‐facing side. Our analysis suggests that accurate composition measurements could be used to infer the kinetic history of hydrate formation in natural settings where gas is abundant.

Rocks are heterogeneous multiscale porous media: two rock samples with identical bulk properties can vary widely in microstructure. The advent of digital rock technology and modern 3‐D printing provides new opportunities to replicate rocks. However, the inherent trade‐off between imaging resolution and sample size limits the scales over which microstructure and macrostructure can be identified and related to each other. Here, we develop a multiscale digital rock construction strategy by combining X‐ray computed microtomography and focused‐ion beam (FIB)‐scanning electron microscope (SEM) images, and we apply the technique to a tight sandstone. The computed tomography (CT) scanning images characterize macroscale pore structures, while the FIB‐SEM images capture microscale pore textures. The FIB‐SEM images are then coupled to CT images via a template‐matching algorithm and superposition. Bulk properties, including porosity and pore and throat size distribution, can be recovered with this approach. Permeability prediction with a pore network model for the largest connected pore network are 3 orders and 1 order of magnitude greater than the bulk rock measured value using the CT‐only and the SEM‐CT coupled images, respectively.

Liquid-vapor phase change is simulated in a porous medium, using a diffuse-interface model and the van der Waals equation of state. The role of wettability and pore space geometry on phase separation and on the propagation of evaporation fronts is investigated.