Breaking Up: Comminution Mechanisms in Sheared Simulated Fault Gouge

The microstructural state and evolution of fault gouge has important implications for the mechanical behaviour, and hence the seismic slip potential of faults. We use 3D discrete element (DEM) simulations to investigate the fragmentation processes operating in fault gouge during shear. Our granular fault gouge models consist of aggregate grains, each composed of several thousand spherical particles stuck together with breakable elastic bonds. The aggregate grains are confined between two blocks of solid material and sheared under a given normal stress. During shear, the grains can fragment in a somewhat realistic way leading to an evolution of grain size, grain shape and overall texture. The ‘breaking up’ of the fault gouge is driven by two distinct comminution mechanisms: grain abrasion and grain splitting. The relative importance of the two mechanisms depends on applied normal stress, boundary wall roughness and accumulated shear strain. If normal stress is sufficiently high, grain splitting contributes significantly to comminution, particularly in the initial stages of the simulations. In contrast, grain abrasion is the dominant mechanism operating in simulations carried out at lower normal stress and is also the main fragmentation mechanism during the later stages of all simulations. Rough boundaries promote relatively more grain splitting whereas smooth boundaries favor grain abrasion. Grain splitting (plus accompanying abrasion) appears to be an efficient mechanism for reducing the mean grain size of the gouge debris and leads rapidly to a power law size distribution with an exponent that increases with strain. Grain abrasion (acting alone) is an effective way to generate excess fine grains and leads to a bimodal distribution of grain sizes. We suggest that these two distinct mechanisms would operate at different stages of a fault’s history. The resulting distributions in grain size and grain shape may significantly affect frictional strength and stability. Our results therefore have implications for the earthquake potential of seismically active faults with accumulations of gouge. They may also be relevant to the susceptibility of rockslides since non-cohesive basal shear zones will evolve in a similar way and potentially control the dynamics of the slide.     

汶川地震断裂带粒度分布特征：对地震碎裂机制的约束

断层岩的粒度分布包含岩石破裂机制、摩擦性质和地震能量分配等重要信息。筛分-称重和激光测量是分析断层岩三维粒度的2种有效方法,但每一种方法的测量范围仅有3个量级,难以全面反映断层岩的粒度分布特征。利用上述2种方法对汶川地震断层滑动带上的断层岩（简称断层岩）的粒度分布进行了测量,粒径测量范围从0.2μm至16mm,跨度达到5个数量级。结果显示：1）存在临界粒径d_c（0.95～1.90μm）。粒度大于和小于d_c的颗粒满足不同的颗粒数（N_d）-粒径（d）分布关系,表明该断层岩的粒径分布不具有自相似性特征。2）利用粒度大于d_c的颗粒计算出的分形维数与断层岩类型有很好的相关性,即断层带边缘的角砾岩的平均分形维数为2.6,核部压碎角砾岩的平均分形维数约为3.0,中心断层泥的分形维数约为3.5。粒径小于d_c的颗粒的分形维数为1.7～2.1。分形维数的突变反映出断层破裂机制的复杂性,即在不同的粒度域,岩石的破裂机制不尽相同。3）依据粒度分析结果,估算出汶川地震断层泥的单位破裂能（E_s）为0.63MJ/m2。     

Multimodal particle size distributions of fine-grained sediments: mathematical modeling and field investigation

Multimodal particle size distributions (PSDs) of fine-grained cohesive sediments are common in marine and coastal environments. The curve-fitting software in this study decomposed such multimodal PSDs into subordinate log-normal PSDs. Four modal peaks, consisting of four-level ordered structures of primary particles, flocculi, microflocs, and macroflocs, were identified and found to alternately rise and sink in a flow-varying tidal cycle due to shear-dependent flocculation. The four modal PSD could be simplified further into two discrete size groups of flocculi and flocs. This allowed the development of a two-class population balance equation (TCPBE) model with flocculi and flocs to simulate flocculation involving multimodal PSDs. The one-dimensional vertical (1-DV) TCPBE model further incorporated the Navier-Stokes equation with the k-ε turbulence closure and the sediment mass balance equations. Multimodal flocculation as well as turbulent flow and sediment transport in a flow-varying tidal cycle could be simulated well using the proposed model. The 1-DV TCPBE was concluded to be the simplest model that is capable of simulating multimodal flocculation in the turbulent flow field of marine and coastal zones.     

海原活动断裂中断层泥的特征、成因及其对断层滑动性能的影响

本文对海原断裂的断层泥的各类特征进行了系统的研究,分析了断层泥带内部的小型断裂类型,探讨了断层泥组分的有机地球化学意义。通过对比,确立了石英碎屑表面机械作用成因的电子显微构造特征组合为活断层的判别标志。本文将石英碎屑表面电子显微溶蚀构造分为五类,据此分析了断层活动时期和状态。研究中还发现了断层泥带内部电子显微结构类型,进而分析了断层活动与微结构、粒度分布间的关系。最后还讨论了断层泥形成演化过程以及断层泥对全新世以来断层滑动性能影响及其与地震活动的关系     

热水条件下角闪石断层泥的摩擦滑动性质——与斜长石断层泥的对比

基于速率与状态依赖性摩擦本构关系理论框架,在热水条件下研究了角闪石断层泥的摩擦滑动性质并与闪长岩的另一种主要矿物斜长石的摩擦滑动性质进行了对比.摩擦实验是在三轴实验系统上完成,有效正应力200 MPa,孔隙压力30 MPa,并将加载速率在1.22 μm/s和 0.122 μm/s之间实施了切换.结果表明角闪石的摩擦系数均值为0.70±0.01,随着温度增加没有系统性的变化,整体低于斜长石的摩擦系数(0.75±0.01);计算与实验表明,角闪石和斜长石的摩擦系数的体积分数加权平均值与闪长岩的摩擦系数基本一致;角闪石在实验温度范围内(100～614 ℃)显示速率强化(a-b>0),与斜长石在整个温度范围内的速率弱化(a-b<0)正好相反;角闪石的速率依赖性在整个实验温度范围内无系统性的变化.     

Triaxial testing of Lopez Fault gouge at 150 MPa mean effective stress

Triaxial compression experiments were performed on samples of natural granular fault gouge from the Lopez Fault in Southern California. This material consists primarily of quartz and has a self-similar grain size distribution thought to result from natural cataclasis. The experiments were performed at a constant mean effective stress of 150 MPa, to expose the volumetric strains associated with shear failure. The failure strength is parameterized by the coefficient of internal friction , based on the Mohr-Coulomb failure criterion.Samples of remoulded Lopez gouge have internal friction =0.6±0.02. In experiments where the ends of the sample are constrained to remain axially aligned, suppressing strain localisation, the sample compacts before failure and dilates persistently after failure. In experiments where one end of the sample is free to move laterally, the strain localises to a single oblique fault at around the point of failure; some dilation occurs but does not persist. A comparison of these experiments suggests that dilation is confined to the region of shear localisation in a sample. Overconsolidated samples have slightly larger failure strengths than normally consolidated samples, and smaller axial strains are required to cause failure. A large amount of dilation occurs after failure in heavily overconsolidated samples, suggesting that dilation is occurring throughout the sample. Undisturbed samples of Lopez gouge, cored from the outcrop, have internal friction in the range =0.4–0.6; the upper end of this range corresponds to the value established for remoulded Lopez gouge. Some kind of natural heterogeneity within the undisturbed samples is probably responsible for their low, variable strength. In samples of simulated gouge, with a more uniform grain size, active cataclasis during axial loading leads to large amounts of compaction. Larger axial strains are required to cause failure in simulated gouge, but the failure strength is similar to that of natural Lopez gouge.Use of the Mohr-Coulomb failure criterion to interpret the results from this study, and other recent studies on intact rock and granular gouge, leads to values of that depend on the loading configuration and the intact or granular state of the sample. Conceptual models are advanced to account for these descrepancies. The consequences for strain-weakening of natural faults are also discussed.     

A heuristic probabilistic approach to estimating size-dependent mobility of nonuniform sediment

This paper presents a heuristic probabilistic approach to estimating the size-dependent mobilities of nonuniform sediment based on the pre- and post-entrainment particle size distributions (PSDs), assuming that the PSDs are lognormally distributed. The approach fits a lognormal probability density function to the pre-entrainment PSD of bed sediment and uses the threshold particle size of incipient motion and the concept of sediment mixture to estimate the PSDs of the entrained sediment and post-entrainment bed sediment. The new approach is simple in physical sense and significantly reduces the complexity and computation time and resource required by detailed sediment mobility models. It is calibrated and validated with laboratory and field data by comparing to the size-dependent mobilities predicted with the existing empirical lognormal cumulative distribution function approach. The novel features of the current approach are: (1) separating the entrained and non-entrained sediments by a threshold particle size, which is a modified critical particle size of incipient motion by accounting for the mixed-size effects, and (2) using the mixture-based pre- and post-entrainment PSDs to provide a continuous estimate of the size-dependent sediment mobility.     

Development of shear localization in simulated quartz gouge: Effect of cumulative slip and gouge particle size

Frictional sliding experiments were conducted on two types of simulated quartz gouge (with median particle diameters 5 m and 25 m, respectively) at confining pressures ranging from 50 MPa to 190 MPa in a conventional triaxial configuration. To investigate the operative micromechanical processes, deformation texture developed in the gouge layer was studied in samples which had accumulated different amounts of frictional slip and undergone different stability modes of sliding. The spatial patterning of shear localization was characterized by a quantitative measurement of the shear band density and orientation. Shear localization in the ultrafine quartz gouge initiated very early before the onset of frictional sliding. Various modes of shear localization were evident, but within the gouge zoneR ₁-shears were predominant. The density of shear localization increased with cumulative slip, whereas the angle subtended at the rock-gouge interface decreased. Destabilization of the sliding behavior in the ultrafine quartz gouge corresponded to the extension ofR ₁-shears and formation of boundaryY-shear segments, whereas stabilization with cumulative slip was related to the coalescence ofY-shear segments to form a throughgoing boundary shear. In the coarse quartz gouge, the sliding behavior was relatively stable, probably because shear localization was inhibited by distributed comminution. Two different models were formulated to analyze the stress field within the gouge zone, with fundamentally different predictions on the orientations of the principal stresses. If the rock-gouge interface is assumed to be bonded without any displacement discontinuity, then the maximum principal stress in the gouge zone is predicted to subtend an angle greater than 45° at the interface. If no assumption on displacement or strain continuity is made and if the gouge has yielded as a Coulomb material, then the maximum principal stress in the gouge zone is predicted to subtend an angle less than 45°. If the apparent friction coefficient increases with overall slip (i.e., slip-hardening), then the Riedel shear angle progressively decreases with increasing shear strain within the gouge layer, possibly attaining a zero value which corresponds to a boundaryY-shear. Our quantitative data on shear localization orientation are in reasonable agreement with this second model, which implies the coefficient of internal friction to be about 0.75 for the ultrafine quartz gouge and 0.8 for the coarse gouge. The wide range of orientations for Riedel shear localization observed in natural faults suggests that the orientations of principal stresses vary as much as in an experimental gouge zone.     

Particle size distribution of cataclastic fault materials from Southern California: A 3-D study

The particle size distributions of fault gouge from the San Andreas, the San Gabriel, and the Lopez Canyon faults in Southern California were measured using sieving and Coulter-Counter techniques over a range of particle sizes from 2 m to 16 mm. The distributions were found to be power law (fractal) for the smaller fragments and log-normal by mass for sizes near and above the peak size. The apparent fractal dimensionD of the smaller particles in gouge samples from the San Andreas fault, the San Gabriel fault and the Lopez Canyon gouge were 2.4–3.6, 2.6–2.9 and 2.4–3.0, respectively. The averageD for the Lopez Canyon gouge was 2.7±0.2, which is in agreement with earlier studies of this gouge using planar 2-D sections. The fractal dimension of the finer fragments from all three faults is observed to be correlated with the peak fragment size, with finer gouges tending to have a largerD. A computer automaton is used to show that this observation may be explained as resulting from a fragmentation process which has a grinding limit at which particle reduction stops.     

Acoustic and vibration precursors from shear-slip characteristics of simulated zigzag-type gouge fault

Earthquake precursors are vitally important for warning and risk evaluation, but shedding light on them is a complex task for seismologists and geologists. We attempted to illuminate the precursory rules of fault failure by traditional biaxial-direction shear tests for a zigzag-type gouge with different size and arrangement of simulated faults manufactured from polymethyl methacrylate materials. Shear stress, acoustic emissions, and vibration signals were recorded and comparatively analyzed. Two important findings were obtained: (1) the linearly increasing trend of shear stress peaks generated by foreshocks may indicate the foreshock effect before a larger earthquake triggered by natural fault gouge failure; and (2) four kinds of earthquake precursors were discovered, and earthquake intensity and precursors may be dependent on the time of macro-fracture formation and the quantity of micro-fractures initiated before mainshock. The study contributes to interpreting earthquake shear-slip characteristics and may even help provide warnings of failure and instability.     

断层泥极限粒度存在的可能机理及其意义

通过对济南辉长岩的双剪摩擦实验发现,在一定的正压力下,存在着相应的断层泥极限粒度。应用颗粒破碎理论及破裂表面能密度进行分析,并结合对模拟断层带的显微观察,作者认为,摩擦功的消耗由颗粒磨细向形成R_1等剪裂面的转变是断层泥极限粒度存在的可能机理。另外,文中还对应用断层泥极限粒度反演其主要形成期围压值的可行性作了初步论证。     

Fault Wear by Damage Evolution During Steady-State Slip

Slip along faults generates wear products such as gouge layers and cataclasite zones that range in thickness from sub-millimeter to tens of meters. The properties of these zones apparently control fault strength and slip stability. Here we present a new model of wear in a three-body configuration that utilizes the damage rheology approach and considers the process as a microfracturing or damage front propagating from the gouge zone into the solid rock. The derivations for steady-state conditions lead to a scaling relation for the damage front velocity considered as the wear-rate. The model predicts that the wear-rate is a function of the shear-stress and may vanish when the shear-stress drops below the microfracturing strength of the fault host rock. The simulated results successfully fit the measured friction and wear during shear experiments along faults made of carbonate and tonalite. The model is also valid for relatively large confining pressures, small damage-induced change of the bulk modulus and significant degradation of the shear modulus, which are assumed for seismogenic zones of earthquake faults. The presented formulation indicates that wear dynamics in brittle materials in general and in natural faults in particular can be understood by the concept of a “propagating damage front” and the evolution of a third-body layer.     

Response of suspended sediment particle size distributions to changes in water chemistry at an Andean mountain stream confluence receiving arsenic rich acid drainage

Acid drainage is an important water quality issue in Andean watersheds, affecting the sustainability of urban, agricultural, and industrial activities. Mixing zones receiving acid drainage are critical sites where changes in pH and chemical environment promote the formation and dissolution of iron and aluminum oxy/hydroxides. These particles can significantly change the speciation of toxic metals and metalloids throughout drainage networks via sorption, desorption, and settling processes. However, little is known about the behavior of particle size distributions (PSDs) in streams affected by acid drainage and their relationship to metal speciation. This work studied: (a) the PSDs for a wide range of mixing ratios found at a fluvial confluence affected by acid drainage, and (b) the response of PSDs and arsenic speciation to environmental changes found when the particles approach complete mixing conditions. The confluence between the Azufre River (pH ~ 2, high concentration of dissolved metals) and Caracarani River (pH = 8.6, low concentration of dissolved metals) was used as a representative model for study. Field measurements show a bimodal PSD with modal diameters of ~50 and 300 μm. At shorter distances from the junction, the smaller modes with smaller particle volumes were dominant across the stream cross‐sections. A systematic shift towards larger particle sizes and larger particle volumes occurred downstream. The analysis of laboratory PSDs for Azufre/Caracarani mixing ratios between 0.01 and 0.5 (pH from 6.2 to 2.3) showed a bimodal trend with ~15 and 50 μm characteristic diameters; larger particles formed at pH>4. When particle suspensions were transferred in laboratory experiments from very low pH to full mixing conditions (pH ~ 2.8 and mixing ratio ~ 0.25) particle sizes varied, and the dissolved arsenic concentration decreased. The observed reaction kinetics were slow compared to the time scale of advective transport, creating opportunities for engineered controls for arsenic. This work contributes to a better understanding of the chemical‐hydrodynamic interactions in watersheds affected by mining, and identifying opportunities to improve water quality at points of use.     

Implications of Coulomb plasticity for the velocity dependence of experimental faults

Simulated fault gouges often deform more stably than initially bare surfaces of the same composition. It is important to understand why the sliding stability is enhanced because the presence of gouge on natural faults may have the same effect as seen in experiments, and thus explain the absence of earthquakes at shallow depths. Gouge stabilization in experiments has been attributed to positive contributions to velocity dependence within gouge layers from either dilation (Marone et al., 1990) or grain fracture (Biegel et al., 1989). In this study we test the hypothesis that some aspects of gouge and initially bare surface velocity dependence are identical by measuring the time-dependent constitutive parameterb. An important result follows however from stress analysis: if both sample configurations are frictional in the Mohr-Coulomb sense, each configuration is required to deform on planes of distinctly different orientation. The measured strength and velocity dependence will reflect this geometric difference. Our observed values ofb for simulated granite and quartz gouge are two to two and a half times smaller thanb for initially bare surfaces. This difference is completely accounted for if gouge is represented as a cohesionless-Coulomb plastic material. The analysis demonstrates the following points: 1) gouge deformation is fully consistent with Coulomb plasticity, 2) observed gouge velocity dependence is a function of observed strength and 3) the constitutive parameterb is the same for both bare surfaces and gouge. Furthermore, the results suggest that there is no time-dependent strengthening associated with stabilizing effects in gouge. These observations provide a framework for understanding how slip on initially bare surfaces and gouge deformation are related.     

Self-similar cataclasis in the formation of fault gouge

Particle-size distributions have been determined for gouge formed by the fresh fracture of granodiorite from the Sierra Nevada batholith, for Pelona schist from the San Andreas fault zone in southern California, and for Berea sandstone from Berea, Ohio, under a variety of triaxial stress states. The finer fractions of the gouge derived from granodiorite and schist are consistent with either a self-similar or a logarithmic normal distribution, whereas the gouge from sandstone is not. Sandstone gouges are texturally similar to the disaggregated protolith, with comminution limited to the polycrystalline fragments and dominantly calcite cement. All three rock types produced significantly less gouge at higher confining pressures, but only the granodiorite showed a significant reduction in particle size with increased confining pressure. Comparison with natural gouges showed that gouges in crystalline rocks from the San Andreas fault zone also tend to be described by either a self-similar or log-normal particle distribution, with a significant reduction in particle size with increased confining pressure (depth). Natural gouges formed in porous sandstone do not follow either a self-similar or a log-normal distribution. Rather, these are represented by mixed log-normal distributions. These textural characteristics are interpreted in terms of the suppression of axial microfracturing by confining pressure and the accommodation of finite strain by scale-independent comminution.     

Evolution of Wear and Friction Along Experimental Faults

We investigate the evolution of wear and friction along experimental faults composed of solid rock blocks. This evolution is analyzed through shear experiments along five rock types, and the experiments were conducted in a rotary apparatus at slip velocities of 0.002–0.97 m/s, slip distances from a few millimeters to tens of meters, and normal stress of 0.25–6.9 MPa. The wear and friction measurements and fault surface observations revealed three evolution phases: A) An initial stage (slip distances <50 mm) of wear by failure of isolated asperities associated with roughening of the fault surface; B) a running-in stage of slip distances of 1–3 m with intense wear-rate, failure of many asperities, and simultaneous reduction of the friction coefficient and wear-rate; and C) a steady-state stage that initiates when the fault surface is covered by a gouge layer, and during which both wear-rate and friction coefficient maintain quasi-constant, low levels. While these evolution stages are clearly recognizable for experimental faults made from bare rock blocks, our analysis suggests that natural faults “bypass” the first two stages and slip at gouge-controlled steady-state conditions.     

A prediction method for the size distribution of saltating particles above a mixed-size bed

In aeolian saltation, the sand bed is a mixture of sand particle with a wide range of particle sizes. Generally, the particle size distribution (PSD) of saltating particles is ignored by previous aeolian transport models, which will result in differences between predictions and observations. To better understand the saltation process, a prediction method of the PSD of saltating particles was proposed in this article. The probability of contact between incident sand and bed sand was introduced into the particle-bed collision process. An artificial PSD of the incident saltating particles was set as the initial condition. A stochastic particle-bed collision model considering contact probability was then used in each iteration step to calculate a new PSD of saltating particles. Finally, the PSD of saltating particles can be determined when aeolian saltation reaches a steady state (saltation is in a steady state when its primary characteristics, such as horizontal mass flux and the concentration of saltating particles, remain approximately constant over time and distance). Meanwhile, according to the experimental results, a calculation formula for the contact parameter n is given, which characterizes the shielding effect of particles on each other. That is, if soil PSD and friction velocity were given, the PSD of saltating particles can be determined. Our results do not depend on the initial conditions, and the predicted results are consistent with the experimental results. It indicated that our method can be used to determine the PSD of saltating particles. © 2020 John Wiley & Sons, Ltd.     

A Continuum Damage–Breakage Faulting Model and Solid-Granular Transitions

We present a thermodynamically-based formulation for mechanical modeling of faulting processes in the seismogenic brittle crust using a continuum damage–breakage rheology. The model combines previous results of a continuum damage framework for brittle solids with continuum breakage mechanics for granular flow. The formulation accounts for the density of distributed cracking and other internal flaws in damaged rocks with a scalar damage parameter, and addresses the grain size distribution of a granular phase in a failure slip zone with a breakage parameter. The stress–strain relation and kinetics of the damage and breakage processes are governed by the total energy function of the system, which combines the energy of the damaged solid with the energy of the granular material. A dynamic brittle instability is associated with a critical level of damage in the solid, leading to loss of convexity of the solid energy function and transition to a granular phase associated with lower energy level. A non-local formulation provides an intrinsic length scale associated with the internal damage structure, which leads to a finite length scale for damage localization that eliminates the unrealistic singular localization of local models. Shear heating during deformation can lead to a secondary finite-width internal localization. The formulation provides a framework for studying multiple aspects of brittle deformation, including potential feedback between evolving elastic moduli and properties of the slip localization zone and subsequent rupture behavior. The model has a more general transition from slow deformation to dynamic rupture than that associated with frictional sliding on a single pre-existing failure zone, and gives time and length scales for the onset of the dynamic fracturing process. Several features including the existence of finite localization width and transition from slow to rapid dynamic slip are illustrated using numerical simulations. A configuration having an existing narrow slip zone with localized damage produces for appropriate loading conditions an overall cyclic stick–slip motion. The simulated frictional response includes transitions from friction coefficient of ~0.7 at low slip velocity to dynamic friction below 0.4 at slip rates above ~0.1 m/s, followed by rapidly increasing friction for slip rates above ~1 m/s, consistent with laboratory observations.     

Simulation of the Influence of Rate- and State-dependent Friction on the Macroscopic Behavior of Complex Fault Zones with the Lattice Solid Model

-- In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by Dieterich and Kilgore (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D_c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.