An official website of the United States government
Abstract Geophysical fluid‐granular flows, such as pyroclastic currents and debris flows, owe much of their runout and hazard behavior to the occurrence and time‐variant decay of a flow‐internal fluid pore pressure. However, modeling the effects of fluid pore pressure to forecast hazards is challenging because a unified method in Earth Sciences to quantitatively determine the permeability of these natural mixtures is currently missing. Here we combine experiments on fluidization and defluidization of pyroclastic materials, eolian sediments, and glass beads mixtures with numerical multiphase simulations to compare previous attempts to compute the permeability of complex natural particle‐fluid mixtures. In analogy to particle‐engineering studies on simple gas‐particle mixtures, we demonstrate that the effective length‐scale in the characterization of the fluid‐particle interaction of complex natural mixtures is the product of the Sauter mean diameter and the particle sphericity. Its use in the Kozeny‐Carman equation allows accurate prediction of mixture permeability, and we suggest the routine calculation of the Sauter mean from grain size distributions of the deposits of geophysical mass flows in Earth Sciences. We also show, through defluidization experiments, that the duration of gas retention in natural mixtures is well described when using the Sauter mean as the effective particle size. Further, we show through multiphase simulations that initial bed expansion extends the pore pressure diffusion timescale up to nine times. These results can be applied to small‐to‐large volume dense pyroclastic currents where the ranges of Sauter mean diameter predict gas retention for long duration and to debris flows and snow avalanches.
more »
« less
The fragmentation-induced fluidisation of pyroclastic density currents
Abstract Pyroclastic density currents (PDCs) are the most lethal volcanic process on Earth. Forecasting their inundation area is essential to mitigate their risk, but existing models are limited by our poor understanding of their dynamics. Here, we explore the role of evolving grain-size distribution in controlling the runout of the most common PDCs, known as block-and-ash flows (BAFs). Through a combination of theory, analysis of deposits and experiments of natural mixtures, we show that rapid changes of the grain-size distribution transported in BAFs result in the reduction of pore volume (compaction) within the first kilometres of their runout. We then use a multiphase flow model to show how the compressibility of granular mixtures leads to fragmentation-induced fluidisation (FIF) and excess pore-fluid pressure in BAFs. This process dominates the first ~2 km of their runout, where the effective friction coefficient is progressively reduced. Beyond that distance, transport is modulated by diffusion of the excess pore pressure. Fragmentation-induced fluidisation provides a physical basis to explain the decades-long use of low effective friction coefficients used in depth-averaged simulations required to match observed flow inundation.
more »
« less
- Award ID(s):
- 1751905
- PAR ID:
- 10498497
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 14
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Test of the Effective Stress Law for Semibrittle Deformation Using Isostatic and Triaxial Load PathsAbstract For brittle friction and rock deformation, the coefficientαin the general effective stress relationσe = σ − αPpcan be approximated as unity with sufficient accuracy. However, it is uncertain ifαdeviates from unity for semibrittle flow when both brittle and intracrystalline‐plastic deformation is involved. We conducted triaxial and isostatic compression experiments on synthetic salt‐rocks (∼300 ppm water) at room temperature to test the effective stress relation in the semibrittle regime using silicone oil and argon gas as pore fluids. Confining and pore pressures were cycled while their difference (differential pressure) was kept constant, such that changes in the mechanical behavior would indicate deviation ofαfrom unity. Microstructural observations were used to determine the dependence ofαon true area of grain contact from asperity yielding. In triaxial compression experiments, semibrittle flow involves grain boundary cracking and sliding, and intragranular dislocation glide and cracking. Flow strength remains constant for changes in pore fluid pressure of more than two orders of magnitude. In isostatic compression experiments, samples show combined processes of microcracking, grain boundary sliding, dislocation glide, and fluid‐assisted grain boundary migration recrystallization. Volumetric strain depends directly on the differential pressures (i.e.,αequals one). Analysis of grain‐contact area in both experiments indicates thatαis independent of the true area of contact defined by plastic yielding at grain boundaries. The observation ofαeffectively equals one may be explained by operation of pressure‐independent intracrystalline‐plastic mechanisms and transmission of pore pressure at grain boundaries through thin fluid films.more » « less
-
This paper derives a closed-form criterion to assess the risk of flowslide runout in loose frictional soil. The derivations rely on a recently proposed framework to simulate pre- and post-failure motion in infinite slopes. An analytical solution of the coupled differential equations capturing flowslide hydromechanics is obtained by specifying them for a perfectly plastic constitutive law. This result enables a comprehensive examination of the factors that control whether the landslide motion, once triggered, autonomously comes to rest (self-regulating behaviour with low mobility) or continues to propagate (self-feeding behaviour with high mobility). It is found that the time history of motion is regulated by non-dimensional property groups reflecting the timescale of excess pore pressure dissipation and the inertial properties of the liquefied zone, which are in turn governed by material (e.g. hydraulic conductivity, dilation coefficient, elastic moduli) and slope properties (e.g. thickness, inclination). The solution is used to build charts identifying the critical ranges of soil properties and triggering factors that differentiate between high-mobility and low-mobility flowslides. Most importantly, it is shown that the fate of flowslide motions is predicted by a critical ratio expressed in terms of excess pore pressure and flow velocity, here defined as the factor of mobility, FM, with values above 1 indicating a self-feeding runout.more » « less
-
Abstract Fault-zone fluids control effective normal stress and fault strength. While most earthquake models assume a fixed pore fluid pressure distribution, geologists have documented fault valving behavior, that is, cyclic changes in pressure and unsteady fluid migration along faults. Here we quantify fault valving through 2-D antiplane shear simulations of earthquake sequences on a strike-slip fault with rate-and-state friction, upward Darcy flow along a permeable fault zone, and permeability evolution. Fluid overpressure develops during the interseismic period, when healing/sealing reduces fault permeability, and is released after earthquakes enhance permeability. Coupling between fluid flow, permeability and pressure evolution, and slip produces fluid-driven aseismic slip near the base of the seismogenic zone and earthquake swarms within the seismogenic zone, as ascending fluids pressurize and weaken the fault. This model might explain observations of late interseismic fault unlocking, slow slip and creep transients, swarm seismicity, and rapid pressure/stress transmission in induced seismicity sequences.more » « less
-
This study outlines a probabilistic cyclic shear strain-based procedure for the determination of the minimum shear strain, γcl, required to initiate liquefaction in gravelly soils. The proposed formulation accounts for the influence of void ratio through the shear wave velocity and the grain size distribution through the coefficient of uniformity, Cu. Separate equations for γcl are derived considering four cyclic resistance models that rely on shear wave velocity as a measure of probabilistic liquefaction resistance. Similarities and differences in the resulting γcl for each of these models are identified. The accuracy and uncertainty of cyclic strain-based models in predicting liquefaction in gravelly soils are demonstrated using existing liquefaction case histories where grain size distributions are available. The excess pore pressure response of gravelly soils subjected to earthquake ground motions is evaluated using a subset of the available liquefaction case histories and the cyclic shear strain and energy-based frameworks and is compared to laboratory test specimens. Although the trends in excess pore pressure generation from critical layers in the case histories are comparable to laboratory-based responses, a greater rate of excess pore pressure generation is calculated for the field cases. The models presented in this study can help identify sites that have a high potential for ground failure when used together with other established models.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1038/s41467-023-37867-1
