An official website of the United States government
Abstract We propose a new approach for performing drained and undrained loading of elastoplastic geomaterials over large deformations using smoothed particle hydrodynamics (SPH), a meshfree continuum particle method, combined with the modified Cam Clay (MCC) model of critical state soil mechanics. The numerical approach draws upon a novel one‐particle two‐phase penalty‐method based formulation for handling undrained loading in saturated soils, which allows tracking of the buildup of pore‐water pressures under combined shearing and compression. Large‐scale parallelized simulations are employed to accommodate a significant number of degrees of freedom in a three‐dimensional setting. After verification and benchmark testing, the SPH based formulation is used to analyze the propagation of reverse faults through fluid‐saturated clay deposits and the rupture of strike‐slip faults across earthen embankments. The computational methodology tests the robustness of the meshfree approach in situations where the soil tends to dilate on the ‘dry’ side of the critical state line and to compact on the ‘wet’ side, but cannot, because of the incompressibility constraint imposed by undrained loading. Our results extend the current understanding of fault rupture modeling and further demonstrate the potential of our framework together with the SPH method for large deformation analyses of complex problems in geotechnics.
more »
« less
Fault rupture propagation through stratified sand-clay deposits and engineered earth structures: a meshfree and critical-state modeling approach
Permanent deformation and uplift caused by fault rupture is one of the most significant hazards posed by earthquakes on the built environment. In this paper, we use Smoothed Particle Hydrodynamics (SPH) to explore the effects of soil layering or stratification on the trajectories and deformation patterns caused by rupturing reverse faults in bedrock, as well as in the foundations of engineered earth structures. SPH is a continuum meshfree numerical method highly adept at modeling large deformation problems in geotechnics. Through the use of constitutive models involving softening behavior as well as critical state type models, we isolate the effects of rigid body rotation from critical state behavior of soil in helping explain the frequently observed rotation of shear bands emanating from the bedrock fault. This analysis is facilitated by the fact that the SPH method allows us to track the propagation of shear bands over substantial amounts of vertical uplift (more than 50\% of the total height of the soil deposit), far beyond many previous computational studies employing the finite element method (FEM). We observe and characterize various emergent features including fault bifurcations, stunted faults, and tension cracking, while providing insights on practical guidelines regarding the potential surface distortion width, and the critical amount of fault displacement required for surface rupture depending on the multilayered constitution of the soil deposit. Finally, we predict the expected amount of surface distortion and internal damage to earthen embankments depending on varying fault location and soil makeup.
more »
« less
- Award ID(s):
- 1914780
- PAR ID:
- 10549176
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Acta geotechnica
- ISSN:
- 0374-0633
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Despite a lack of modern large earthquakes on shallowly dipping normal faults, Holocene M w > 7 low-angle normal fault (LANF; dip<30°) ruptures are preserved paleoseismically and inferred from historical earthquake and tsunami accounts. Even in well-recorded megathrust earthquakes, the effects of non-linear off-fault plasticity and dynamically reactivated splay faults on shallow deformation and surface displacements, and thus hazard, remain elusive. We develop data-constrained 3D dynamic rupture models of the active Mai’iu LANF that highlight how multiple dynamic shallow deformation mechanisms compete during large LANF earthquakes. We show that shallowly-dipping synthetic splays host more coseismic slip and limit shallow LANF rupture more than steeper antithetic splays. Inelastic hanging-wall yielding localizes into subplanar shear bands indicative of newly initiated splay faults, most prominently above LANFs with thick sedimentary basins. Dynamic splay faulting and sediment failure limit shallow LANF rupture, modulating coseismic subsidence patterns, near-shore slip velocities, and the seismic and tsunami hazards posed by LANF earthquakes.more » « less
-
Abstract The Eastern California shear zone is a complex set of dextral faults that accommodates significant plate motion and has produced large earthquakes. The evolution of this system and why it consists of closely spaced, irregular faults that fail in multi‐fault ruptures are not well understood. Here we analyze the geometry, spatial distribution, and Quaternary slip activity of right‐lateral faults in the southern Mojave block. We find these faults are oriented favorably for accommodating regional dextral plate motion and do not show evidence of replacement following counterclockwise rotation to unfavorable positions, although activity may be migrating westward as previously proposed. We also confirm that the shear zone is transpressive, with widespread restraining bends, distributed convergent deformation, and significant impact on near‐fault topography. Observations also show that faults are geometrically complex, as represented by along‐strike variability in fault strike. We document a correlation between strike variability and fault activity (slip rate or net slip), which is evident within the shear zone as well as for a control group of other faults. We suggest that strike variability represents a form of geometric roughness, which may inhibit fault slip and result in complex ruptures, slip‐strengthening behavior, and a prevalence of off‐fault deformation. Other factors, including preexisting crustal fabric, edge effects, and changes in the stress field, may further complicate kinematics. These results suggest that faults of the shear zone are still juvenile and somewhat unique, yet offer an important window into how broadly distributed shear may evolve into a through‐going continental transform system.more » « less
-
Abstract Determining conditions for earthquake slip on faults is a key goal of fault mechanics highly relevant to seismic hazard. Previous studies have demonstrated that enhanced dynamic weakening (EDW) can lead to dynamic rupture of faults with much lower shear stress than required for rupture nucleation. We study the stress conditions before earthquake ruptures of different sizes that spontaneously evolve in numerical simulations of earthquake sequences on rate‐and‐state faults with EDW due to thermal pressurization of pore fluids. We find that average shear stress right before dynamic rupture (aka shear prestress) systematically varies with the rupture size. The smallest ruptures have prestress comparable to the local shear stress required for nucleation. Larger ruptures weaken the fault more, propagate over increasingly under‐stressed areas due to dynamic stress concentration, and result in progressively lower average prestress over the entire rupture. The effect is more significant in fault models with more efficient EDW. We find that, as a result, fault models with more efficient weakening produce fewer small events and result in systematically lower b‐values of the frequency‐magnitude event distributions. The findings (a) illustrate that large earthquakes can occur on faults that appear not to be critically stressed compared to stresses required for slip nucleation; (b) highlight the importance of finite‐fault modeling in relating the local friction behavior determined in the lab to the field scale; and (c) suggest that paucity of small events or seismic quiescence may be the observational indication of mature faults that operate under low shear stress due to EDW.more » « less
-
ABSTRACT We seek to improve our understanding of the physical processes that control the style, distribution, and intensity of ground surface ruptures on thrust and reverse faults during large earthquakes. Our study combines insights from coseismic ground surface ruptures in historic earthquakes and patterns of deformation in analog sandbox fault experiments to inform the development of a suite of geomechanical models based on the distinct element method (DEM). We explore how model parameters related to fault geometry and sediment properties control ground deformation characteristics such as scarp height, width, dip, and patterns of secondary folding and fracturing. DEM is well suited to this investigation because it can effectively model the geologic processes of faulting at depth in cohesive rocks, as well as the granular mechanics of soil and sediment deformation in the shallow subsurface. Our results show that localized fault scarps are most prominent in cases with strong sediment on steeply dipping faults, whereas broader deformation is prominent in weaker sediment on shallowly dipping faults. Based on insights from 45 experiments, the key parameters that influence scarp morphology include the amount of accumulated slip on a fault, the fault dip, and the sediment strength. We propose a fault scarp classification system that describes the general patterns of surface deformation observed in natural settings and reproduced in our models, including monoclinal, pressure ridge, and simple scarps. Each fault scarp type is often modified by hanging-wall collapse. These results can help to guide both deterministic and probabilistic assessment in fault displacement hazard analysis.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
