A process-based approach to understanding and managing triggered seismicity.

There is growing concern about seismicity triggered by human activities, whereby small increases in stress bring tectonically loaded faults to failure. Examples of such activities include mining, impoundment of water, stimulation of geothermal fields, extraction of hydrocarbons and water, and the injection of water, CO2 and methane into subsurface reservoirs1. In the absence of sufficient information to understand and control the processes that trigger earthquakes, authorities have set up empirical regulatory monitoring-based frameworks with varying degrees of success2,3. Field experiments in the early 1970s at the Rangely, Colorado (USA) oil field4 suggested that seismicity might be turned on or off by cycling subsurface fluid pressure above or below a threshold. Here we report the development, testing and implementation of a multidisciplinary methodology for managing triggered seismicity using comprehensive and detailed information about the subsurface to calibrate geomechanical and earthquake source physics models. We then validate these models by comparing their predictions to subsequent observations made after calibration. We use our approach in the Val d'Agri oil field in seismically active southern Italy, demonstrating the successful management of triggered seismicity using a process-based method applied to a producing hydrocarbon field. Applying our approach elsewhere could help to manage and mitigate triggered seismicity.

Dynamic imaging of force chains in 3D granular media.

Granular media constitute the most abundant form of solid matter on Earth and beyond. When external forces are applied to a granular medium, the forces are transmitted through it via chains of contacts among grains-force chains. Understanding the spatial structure and temporal evolution of force chains constitutes a fundamental goal of granular mechanics. Here, we introduce an experimental technique, interference optical projection tomography, to study force chains in three-dimensional (3D) granular packs under triaxial shear loads and illustrate the technique with random assemblies of spheres and icosahedra. We find that, in response to an increasing vertical load, the pack of spheres forms intensifying vertical force chains, while the pack of icosahedra forms more interconnected force-chain networks. This provides microscopic insights into why particles with more angularity are more resistant to shear failure-the interconnected force-chain network is stronger (that is, more resilient to topological collapse) than the isolated force chains in round particles. The longer force chains with less branching in the pack of round particles are more likely to buckle, which leads to the macroscopic failure of the pack. This work paves the way for understanding the grain-scale underpinning of localized failure of 3D granular media, such as shear localization in landslides and stick-slip frictional motion in tectonic and induced earthquakes.

Wetting transition and fluid trapping in a microfluidic fracture.

Immiscible fluid-fluid displacement in confined geometries is a fundamental process occurring in many natural phenomena and technological applications, from geological CO2 sequestration to microfluidics. Due to the interactions between the fluids and the solid walls, fluid invasion undergoes a wetting transition from complete displacement at low displacement rates to leaving a film of the defending fluid on the confining surfaces at high displacement rates. While most real surfaces are rough, fundamental questions remain about the type of fluid-fluid displacement that can emerge in a confined, rough geometry. Here, we study immiscible displacement in a microfluidic device with a precisely controlled structured surface as an analogue for a rough fracture. We analyze the influence of the degree of surface roughness on the wetting transition and the formation of thin films of the defending liquid. We show experimentally, and rationalize theoretically, that roughness affects both the stability and dewetting dynamics of thin films, leading to distinct late-time morphologies of the undisplaced (trapped) fluid. Finally, we discuss the implications of our observations for geologic and technological applications.

Soft matter physics of the ground beneath our feet.

The soft part of the Earth's surface - the ground beneath our feet - constitutes the basis for life and natural resources, yet a general physical understanding of the ground is still lacking. In this critical time of climate change, cross-pollination of scientific approaches is urgently needed to better understand the behavior of our planet's surface. The major topics in current research in this area cross different disciplines, spanning geosciences, and various aspects of engineering, material sciences, physics, chemistry, and biology. Among these, soft matter physics has emerged as a fundamental nexus connecting and underpinning many research questions. This perspective article is a multi-voice effort to bring together different views and approaches, questions and insights, from researchers that work in this emerging area, the soft matter physics of the ground beneath our feet. In particular, we identify four major challenges concerned with the dynamics in and of the ground: (I) modeling from the grain scale, (II) near-criticality, (III) bridging scales, and (IV) life. For each challenge, we present a selection of topics by individual authors, providing specific context, recent advances, and open questions. Through this, we seek to provide an overview of the opportunities for the broad Soft Matter community to contribute to the fundamental understanding of the physics of the ground, strive towards a common language, and encourage new collaborations across the broad spectrum of scientists interested in the matter of the Earth's surface.

Crossover from viscous fingering to fracturing in cohesive wet granular media: a photoporomechanics study.

We study fluid-induced deformation and fracture of cohesive granular media, and apply photoporomechanics to uncover the underpinning grain-scale mechanics. We fabricate photoelastic spherical particles of diameter d = 2 mm, and make a monolayer granular pack with tunable intergranular cohesion in a circular Hele-Shaw cell that is initially filled with viscous silicone oil. We inject water into the oil-filled photoelastic granular pack, varying the injection flow rate, defending-fluid viscosity, and intergranular cohesion. We find two different modes of fluid invasion: viscous fingering, and fracturing with leak-off of the injection fluid. We directly visualize the evolving effective stress field through the particles' photoelastic response, and discover a hoop effective stress region behind the water invasion front, where we observe tensile force chains in the circumferential direction. Outside the invasion front, we observe compressive force chains aligning in the radial direction. We conceptualize the system's behavior by means of a two-phase poroelastic continuum model. The model captures granular pack dilation and compaction with the boundary delineated by the invasion front, which explains the observed distinct alignments of the force chains. Finally, we rationalize the crossover from viscous fingering to fracturing by comparing the competing forces behind the process: viscous force from fluid injection that drives fractures, and intergranular cohesion and friction that resist fractures.

Crustal fingering facilitates free-gas methane migration through the hydrate stability zone.

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.

Restoring universality to the pinch-off of a bubble.

The pinch-off of a bubble is an example of the formation of a singularity, exhibiting a characteristic separation of length and time scales. Because of this scale separation, one expects universal dynamics that collapse into self-similar behavior determined by the relative importance of viscous, inertial, and capillary forces. Surprisingly, however, the pinch-off of a bubble in a large tank of viscous liquid is known to be nonuniversal. Here, we show that the pinch-off dynamics of a bubble confined in a capillary tube undergo a sequence of two distinct self-similar regimes, even though the entire evolution is controlled by a balance between viscous and capillary forces. We demonstrate that the early-time self-similar regime restores universality to bubble pinch-off by erasing the system's memory of the initial conditions. Our findings have important implications for bubble/drop generation in microfluidic devices, with applications in inkjet printing, medical imaging, and synthesis of particulate materials.

Comprehensive comparison of pore-scale models for multiphase flow in porous media.

Multiphase flows in porous media are important in many natural and industrial processes. Pore-scale models for multiphase flows have seen rapid development in recent years and are becoming increasingly useful as predictive tools in both academic and industrial applications. However, quantitative comparisons between different pore-scale models, and between these models and experimental data, are lacking. Here, we perform an objective comparison of a variety of state-of-the-art pore-scale models, including lattice Boltzmann, stochastic rotation dynamics, volume-of-fluid, level-set, phase-field, and pore-network models. As the basis for this comparison, we use a dataset from recent microfluidic experiments with precisely controlled pore geometry and wettability conditions, which offers an unprecedented benchmarking opportunity. We compare the results of the 14 participating teams both qualitatively and quantitatively using several standard metrics, such as fractal dimension, finger width, and displacement efficiency. We find that no single method excels across all conditions and that thin films and corner flow present substantial modeling and computational challenges.

Rheology of bacterial superfluids in viscous environments.

Viscous environments are ubiquitous in nature and in engineering applications, from mucus in lungs to oil recovery strategies in the earth's subsurface - and in all these environments, bacteria also thrive. The behavior of bacteria in viscous environments has been investigated for a single bacterium, but not for active suspensions. Dense populations of pusher-type bacteria are known to create superfluidic regimes where the effective viscosity of the entire suspension is reduced through collective motion, and the main purpose of this study is to investigate how a viscous environment will affect this behavior. Using a Couette rheometer, we measure shear stress as a function of the applied shear rate to define the effective viscosity of suspensions of Escherichia coli (E. coli), while varying both the bacterial density within the suspension and the viscosity of the suspending fluid. We document the remarkable observation that E. coli decreases the effective suspension viscosity to near-zero (superfluidic regime) for all solvent viscosities tested (1-17 mPa s). Specifically, we observe that the bacterial density needed to trigger this superfluidic regime and the maximum shear rate under which this regime can be sustained both decrease with increasing solvent viscosity. We find that the resulting rheograms can be well approximated by the Carreau-Yasuda law. Using this, we propose a constitutive model as a function of the solvent viscosity and the bacterial concentration only. This model captures the onset of the superfluidic regime and offers promising avenues for the modelling of flow of bacterial suspensions in viscous environments.

Hand-Hygiene Mitigation Strategies Against Global Disease Spreading through the Air Transportation Network.

The risk for a global transmission of flu-type viruses is strengthened by the physical contact between humans and accelerated through individual mobility patterns. The Air Transportation System plays a critical role in such transmissions because it is responsible for fast and long-range human travel, while its building components-the airports-are crowded, confined areas with usually poor hygiene. Centers for Disease Control and Prevention (CDC) and World Health Organization (WHO) consider hand hygiene as the most efficient and cost-effective way to limit disease propagation. Results from clinical studies reveal the effect of hand washing on individual transmissibility of infectious diseases. However, its potential as a mitigation strategy against the global risk for a pandemic has not been fully explored. Here, we use epidemiological modeling and data-driven simulations to elucidate the role of individual engagement with hand hygiene inside airports in conjunction with human travel on the global spread of epidemics. We find that, by increasing travelers engagement with hand hygiene at all airports, a potential pandemic can be inhibited by 24% to 69%. In addition, we identify 10 airports at the core of a cost-optimal deployment of the hand-washing mitigation strategy. Increasing hand-washing rate at only those 10 influential locations, the risk of a pandemic could potentially drop by up to 37%. Our results provide evidence for the effectiveness of hand hygiene in airports on the global spread of infections that could shape the way public-health policy is implemented with respect to the overall objective of mitigating potential population health crises.

Wettability control on multiphase flow in patterned microfluidics.

Multiphase flow in porous media is important in many natural and industrial processes, including geologic CO2 sequestration, enhanced oil recovery, and water infiltration into soil. Although it is well known that the wetting properties of porous media can vary drastically depending on the type of media and pore fluids, the effect of wettability on multiphase flow continues to challenge our microscopic and macroscopic descriptions. Here, we study the impact of wettability on viscously unfavorable fluid-fluid displacement in disordered media by means of high-resolution imaging in microfluidic flow cells patterned with vertical posts. By systematically varying the wettability of the flow cell over a wide range of contact angles, we find that increasing the substrate's affinity to the invading fluid results in more efficient displacement of the defending fluid up to a critical wetting transition, beyond which the trend is reversed. We identify the pore-scale mechanisms-cooperative pore filling (increasing displacement efficiency) and corner flow (decreasing displacement efficiency)-responsible for this macroscale behavior, and show that they rely on the inherent 3D nature of interfacial flows, even in quasi-2D media. Our results demonstrate the powerful control of wettability on multiphase flow in porous media, and show that the markedly different invasion protocols that emerge-from pore filling to postbridging-are determined by physical mechanisms that are missing from current pore-scale and continuum-scale descriptions.

Nonequilibrium Thermodynamics of Hydrate Growth on a Gas-Liquid Interface.

We develop a continuum-scale phase-field model to study gas-liquid-hydrate systems far from thermodynamic equilibrium. We design a Gibbs free energy functional for methane-water mixtures that recovers the isobaric temperature-composition phase diagram under thermodynamic equilibrium conditions. The proposed free energy is incorporated into a phase-field model to study the dynamics of hydrate formation on a gas-liquid interface. We elucidate the role of initial aqueous concentration in determining the direction of hydrate growth at the interface, in agreement with experimental observations. Our model also reveals two stages of hydrate growth at an interface-controlled by a crossover in how methane is supplied from the gas and liquid phases-which could explain the persistence of gas conduits in hydrate-bearing sediments and other nonequilibrium phenomena commonly observed in natural methane hydrate systems.

Forced Wetting Transition and Bubble Pinch-Off in a Capillary Tube.

Immiscible fluid-fluid displacement in partial wetting continues to challenge our microscopic and macroscopic descriptions. Here, we study the displacement of a viscous fluid by a less viscous fluid in a circular capillary tube in the partial wetting regime. In contrast with the classic results for complete wetting, we show that the presence of a moving contact line induces a wetting transition at a critical capillary number that is contact angle dependent. At small displacement rates, the fluid-fluid interface deforms slightly from its equilibrium state and moves downstream at a constant velocity, without changing its shape. As the displacement rate increases, however, a wetting transition occurs: the interface becomes unstable and forms a finger that advances along the axis of the tube, leaving the contact line behind, separated from the meniscus by a macroscopic film of the viscous fluid on the tube wall. We describe the dewetting of the entrained film, and show that it universally leads to bubble pinch-off, therefore demonstrating that the hydrodynamics of contact line motion generate bubbles in microfluidic devices, even in the absence of geometric constraints.

Wettability control of droplet durotaxis.

Durotaxis refers to cell motion directed by stiffness gradients of an underlying substrate. Recent work has shown that droplets also move spontaneously along stiffness gradients through a process reminiscent of durotaxis. Wetting droplets, however, move toward softer substrates, an observation seemingly at odds with cell motion. Here, we extend our understanding of this phenomenon, and show that wettability of the substrate plays a critical role: while wetting droplets move in the direction of lower stiffness, nonwetting liquids reverse droplet durotaxis. Our numerical experiments also reveal that Laplace pressure can be used to determine the direction of motion of liquid slugs in confined environments. Our results suggest new ways of controlling droplet dynamics at small scales, which can open the door to enhanced bubble and droplet logic in microfluidic platforms.

Thin Films in Partial Wetting: Internal Selection of Contact-Line Dynamics.

When a liquid touches a solid surface, it spreads to minimize the system's energy. The classic thin-film model describes the spreading as an interplay between gravity, capillarity, and viscous forces, but it cannot see an end to this process as it does not account for the nonhydrodynamic liquid-solid interactions. While these interactions are important only close to the contact line, where the liquid, solid, and gas meet, they have macroscopic implications: in the partial-wetting regime, a liquid puddle ultimately stops spreading. We show that by incorporating these intermolecular interactions, the free energy of the system at equilibrium can be cast in a Cahn-Hilliard framework with a height-dependent interfacial tension. Using this free energy, we derive a mesoscopic thin-film model that describes the statics and dynamics of liquid spreading in the partial-wetting regime. The height dependence of the interfacial tension introduces a localized apparent slip in the contact-line region and leads to compactly supported spreading states. In our model, the contact-line dynamics emerge naturally as part of the solution and are therefore nonlocally coupled to the bulk flow. Surprisingly, we find that even in the gravity-dominated regime, the dynamic contact angle follows the Cox-Voinov law.

Lifetime of carbon capture and storage as a climate-change mitigation technology.

In carbon capture and storage (CCS), CO(2) is captured at power plants and then injected underground into reservoirs like deep saline aquifers for long-term storage. While CCS may be critical for the continued use of fossil fuels in a carbon-constrained world, the deployment of CCS has been hindered by uncertainty in geologic storage capacities and sustainable injection rates, which has contributed to the absence of concerted government policy. Here, we clarify the potential of CCS to mitigate emissions in the United States by developing a storage-capacity supply curve that, unlike current large-scale capacity estimates, is derived from the fluid mechanics of CO(2) injection and trapping and incorporates injection-rate constraints. We show that storage supply is a dynamic quantity that grows with the duration of CCS, and we interpret the lifetime of CCS as the time for which the storage supply curve exceeds the storage demand curve from CO(2) production. We show that in the United States, if CO(2) production from power generation continues to rise at recent rates, then CCS can store enough CO(2) to stabilize emissions at current levels for at least 100 y. This result suggests that the large-scale implementation of CCS is a geologically viable climate-change mitigation option in the United States over the next century.

Synergetic fluid mixing from viscous fingering and alternating injection.

We study mixing of two fluids of different viscosity in a microfluidic channel or porous medium. We show that the synergetic action of alternating injection and viscous fingering leads to a dramatic increase in mixing efficiency at high Péclet numbers. Based on observations from high-resolution simulations, we develop a theoretical model of mixing efficiency that combines a hyperbolic mixing model of the channelized region ahead and a mixing-dissipation model of the pseudosteady region behind. Our macroscopic model quantitatively reproduces the evolution of the average degree of mixing along the flow direction and can be used as a design tool to optimize mixing from viscous fingering in a microfluidic channel.

Lab on a chip for a low-carbon future.

Transitioning our society to a sustainable future, with low or net-zero carbon emissions to the atmosphere, will require a wide-spread transformation of energy and environmental technologies. In this perspective article, we describe how lab-on-a-chip (LoC) systems can help address this challenge by providing insight into the fundamental physical and geochemical processes underlying new technologies critical to this transition, and developing the new processes and materials required. We focus on six areas: (I) subsurface carbon sequestration, (II) subsurface hydrogen storage, (III) geothermal energy extraction, (IV) bioenergy, (V) recovering critical materials, and (VI) water filtration and remediation. We hope to engage the LoC community in the many opportunities within the transition ahead, and highlight the potential of LoC approaches to the broader community of researchers, industry experts, and policy makers working toward a low-carbon future.

Macroscopic phase-field model of partial wetting: bubbles in a capillary tube.

Drops and bubbles are nonspreading, local, compactly supported features. They are also equilibrium configurations in partial wetting phenomena. Yet, current macroscopic theories of capillary-dominated flow are unable to describe these systems. We propose a framework to model multiphase flow in porous media with nonspreading equilibrium configurations. We illustrate our approach with a one-dimensional model of two-phase flow in a capillary tube. Our model allows for the presence of compactons: nonspreading steady-state solutions in the absence of external forces. We show that local rate dependency is not needed to explain globally rate-dependent displacement patterns, and we interpret dynamic wetting transitions as the route from equilibrium, capillary-dominated configurations, towards viscous-dominated flow. Mathematically, these transitions are possible due to nonclassical shock solutions and the role of bistability and higher-order terms in our model.

Capillary fracturing in granular media.

We study the displacement of immiscible fluids in deformable, noncohesive granular media. Experimentally, we inject air into a thin bed of water-saturated glass beads and observe the invasion morphology. The control parameters are the injection rate, the bead size, and the confining stress. We identify three invasion regimes: capillary fingering, viscous fingering, and "capillary fracturing," where capillary forces overcome frictional resistance and induce the opening of conduits. We derive two dimensionless numbers that govern the transition among the different regimes: a modified capillary number and a fracturing number. The experiments and analysis predict the emergence of fracturing in fine-grained media under low confining stress, a phenomenon that likely plays a fundamental role in many natural processes such as primary oil migration, methane venting from lake sediments, and the formation of desiccation cracks.