Slip-weakening Distance Research Articles

Effects of fault roughness on estimating critical slip-weakening distance from fault slip history: A laboratory study

Earthquakes are the dynamic rupture of faults governed by fault weakening processes. Critical slip-weakening distance (Dc) is a crucial source parameter of earthquakes, and the determination of Dc is of great concern to semiologists. However, determining Dc for natural earthquakes is challenging due to the trade-off in inversed source models. To solve this problem, Fukuyama and his coworkers proposed a simple method (denoted as the F&M method) to estimate Dc directly from slip-velocity functions. According to the F&M method, the fault slip at the peak slip velocity (Dc') can be used as an approximation of Dc. However, the feasibility of this method has not been completely resolved. Here, we performed laboratory earthquake rupture experiments to examine the validity of the F&M method. Experiments were conducted on faults with different roughness and stress level. We studied the effect of fault roughness on the validity of the F&M method. Our results show that the increase in fault roughness could complicate the fault weakening process, producing repeated weakening and strengthening phase, which leads to a prominent deviation between Dc' and Dc. Furthermore, we also observed a correlation between the critical slip-weakening distance Dc and the final fault slip D. Such correlation implies the scale-dependent nature of Dc, which is consistent with seismic observations in the field.

Fractal Spatial Distributions of Initial Shear Stress and Frictional Properties on Faults and Their Impact on Dynamic Earthquake Rupture

ABSTRACT We investigate the influence of the heterogeneous slip-weakening distance (DC) in dynamic rupture simulations, in which DC is proportional to the fault irregularities. Specifically, we compare a heterogeneous fractal DC distribution to a uniform DC over the entire fault when the initial shear stress is also heterogeneous. We find that even small changes in the average value of DC (<1 mm) can lead to significant differences in the rupture evolution; that is, the average DC and the way DC is distributed determines if the rupture is a runaway, self-arrested, or nonpropagating. We find that the self-arrested ruptures differ from runaway ruptures in the amount of area characterized by large slips (asperities). Self-arrested ruptures match the Somerville et al. (1999) asperity criteria in which ∼25% of ruptured area radiate ∼45% of the total seismic moment. This criterion is not satisfied for runaway ruptures. For runaway ruptures, ∼50% of the ruptured area radiates about 70% of the seismic moment, indicating that the ruptured area is not linearly proportional to the seismic moment. Self-arrested ruptures are characterized by dynamic shear stress drops (SDs) in the range ∼2.9–5.5 MPa, whereas for runaway ruptures the dynamic SDs increase to values between ∼12 and 20 MPa. Self-arrested ruptures generated by fractal distributed DC resemble the rupture properties of observed earthquakes. In addition, results show that the conditions for self-arrested ruptures are connected to the decrease of residual energy at rupture boundaries.

Dynamic Rupture Models of the 2016 ML 5.8 Gyeongju, South Korea, Earthquake, Constrained by a Kinematic Rupture Model and Seismic Waveform Data

ABSTRACT The ML 5.8 earthquake that jolted Gyeongju in southeastern Korea in 2016 was the country’s largest inland event since instrumental seismic monitoring began in 1978. We developed dynamic rupture models of the Gyeongju event constrained by near-source ground-motion data using full 3D spontaneous dynamic rupture modeling with the slip-weakening friction law. Based on our results, we propose two simple dynamic rupture models with constant strength excess (SE) and slip-weakening distance (Dc) that produce near-source ground-motion waveforms compatible with recorded ones in the low-frequency band. Both dynamic models exhibit relatively large stress-drop values, consistent with previous estimates constrained by source spectrum analyses. The fracture energy estimates were also larger than those predicted by a scaling relationship with the seismic moment. The dynamic features constrained in this study by spontaneous rupture modeling and waveform comparison may help understand the source and ground-motion characteristics of future large events in southeastern Korea and thus the seismic hazard of the region.

Hydro-mechanical earthquake cycles in a poro-visco-elasto-plastic fluid-bearing fault structure

A major goal in earthquake physics is to derive a constitutive framework for fault slip that captures the dependence of shear strength on fault rheology, sliding velocity, and pore-fluid pressure. In this study, we present H-MEC (Hydro-Mechanical Earthquake Cycles), a newly-developed two-phase flow numerical framework — which couples solid rock deformation and pervasive fluid flow — to simulate how crustal stress and fluid pressure evolve during the earthquake cycle. This unified, continuum-based model, incorporates a staggered finite difference–marker-in-cell method and accounts for inertial wave-mediated dynamics and fluid flow in poro-visco-elasto-plastic compressible medium. Global Picard-iterations and an adaptive time stepping allow the correct resolution of both long- and short-time scales, ranging from years to milliseconds. We present a comprehensive in-plane strike-slip setup in which we test analytical poroelastic benchmarks of pore-fluid pressure diffusion from an injection point along a finite fault width. We then investigate how pore-fluid pressure evolution and solid–fluid compressibility control sequences of seismic and aseismic fault slip. While the weakening phase is controlled by localized compaction of pores and dynamic self-pressurization of fluids inside the undrained fault zone, the subsequent propagation of dynamic ruptures is driven by pore-pressure waves. Furthermore, pore-fluid pressure conditions on the fault and shear strength weakening associated with rapid self-pressurization of fluids control the characteristic slip-weakening distance, the final size of seismic events, and the scaling between slip and fracture energy observed for large earthquakes. Our modeling results demonstrate that fault failure can occur due to poroelastic coupling on a finite-width shear zone, thus highlighting the importance of considering the realistic hydro-mechanical structure of faults to investigate fluid-driven seismic and aseismic slip, either as a natural process or induced by human activities.

Deciphering the Origins of Transient Seismic Moment Accelerations by Realistic Dynamic Rupture Simulations

ABSTRACTProperties of earthquake source physics can be inferred from the comparison between seismic observations and results of dynamic rupture models. Although simple self-similar rupture models naturally explain the space and time observations at the scale of the whole earthquake, several observational studies based on the analysis of source time functions (STFs) suggest that they are unable to reproduce the initial accelerating phases of the rupture. We here propose to reproduce the observed transient moment accelerations, without affecting the global self-similarity of the rupture, to constrain their possible physical origins. Simulated STFs are generated from dynamic simulations with heterogeneous slip-weakening distance Dc. Heterogeneity is introduced on the fault plane through a fractal number-size distribution of circular patches, in which Dc takes a value proportional to their radius. As a consequence of the stochastic spatial distribution of the patches, rupture development exhibits a large variability, and delays between initiation and main rupture activation frequently occur. This variability, together with the dynamic correlation between rupture velocity and slip velocity inside each broken patch, successfully perturbs the self-similar properties: rather than growing quadratically with time, STFs have an higher apparent time exponent, close to the observed value of 2.7. In a broader perspective, our simulations show that to respect observed STF shapes, realistic dynamic models should generate bursts of seismic moment, most likely by episodes where slip and rupture velocity are correlated. Such a behavior appears to emerge more naturally when considering heterogeneities in the friction parameters rather than in the initial stress.

Repeating caldera collapse events constrain fault friction at the kilometer scale

Fault friction is central to understanding earthquakes, yet laboratory rock mechanics experiments are restricted to, at most, meter scale. Questions thus remain as to the applicability of measured frictional properties to faulting in situ. In particular, the slip-weakening distance [Formula: see text] strongly influences precursory slip during earthquake nucleation, but scales with fault roughness and is challenging to extrapolate to nature. The 2018 eruption of K̄ılauea volcano, Hawaii, caused 62 repeatable collapse events in which the summit caldera dropped several meters, accompanied by [Formula: see text] 4.7 to 5.4 very long period (VLP) earthquakes. Collapses were exceptionally well recorded by global positioning system (GPS) and tilt instruments and represent unique natural kilometer-scale friction experiments. We model a piston collapsing into a magma reservoir. Pressure at the piston base and shear stress on its margin, governed by rate and state friction, balance its weight. Downward motion of the piston compresses the underlying magma, driving flow to the eruption. Monte Carlo estimation of unknowns validates laboratory friction parameters at the kilometer scale, including the magnitude of steady-state velocity weakening. The absence of accelerating precollapse deformation constrains [Formula: see text] to be [Formula: see text] mm, potentially much less. These results support the use of laboratory friction laws and parameters for modeling earthquakes. We identify initial conditions and material and magma-system parameters that lead to episodic caldera collapse, revealing that small differences in eruptive vent elevation can lead to major differences in eruption volume and duration. Most historical basaltic caldera collapses were, at least partly, episodic, implying that the conditions for stick-slip derived here are commonly met in nature.

Inferring Critical Slip-Weakening Distance from Near-Fault Accelerogram of the 2014 Mw 6.2 Ludian Earthquake

Abstract To better assess potential earthquake hazards requires a better understanding of fault friction and rupture dynamics. Critical slip-weakening distance (Dc) as one of the key friction parameters, however, is hard to determine on natural faults. For strike-slip earthquakes, we may directly estimate the Dc from Dc″—the double near-fault ground displacement at the time of the peak velocity (Fukuyama and Mikumo, 2007). Yet near-fault observations are very few, and, thus, there were only limited earthquakes with such Dc″ estimation. In 2014, an Mw 6.2 strike-slip event—the Ludian earthquake—occurred in southwest China. The strong-motion station (LLT) that is ∼0.45 km from the fault recorded the earthquake and enabled us to estimate Dc″ from the accelerograms. We inspect the polarity of the accelerometers and compare the integrated velocities with waveforms of nearby broadband stations. We also analyze the particle motion at the LLT station and retrieve the earthquake initiation at the intersection of the conjugated faults. We then apply the baseline correction to the seismograms, recover the ground velocities and displacements, and obtain the value of Dc″=0.1 m at the station. The recovered final displacements are compared with the predicted ground displacements of a finite-fault model. The discrepancy of fault-parallel displacements might imply limited underestimation of Dc″, and the estimated upper limit is 0.3 m. Comparison between the Dc″ and final slip on the fault patch follows the scaling of previous larger earthquakes. Analysis of the near-fault accelerometer data enhances our understanding on the earthquake source of the Ludian earthquake. This case extends the lower magnitude boundary of the Dc″ values obtained from natural faults and opens a window into the friction property in the seismically active region.

Spatiotemporal Evolution of Microseismicity Seismic Source Properties at the Irpinia Near-Fault Observatory, Southern Italy

ABSTRACTWe estimate the source parameters of small-magnitude earthquakes that occurred during 2008–2020 in the Irpinia faults area (southern Italy). We apply a spectral decomposition approach to isolate the source contribution from propagation and site effects for ∼3000 earthquakes in the local magnitude range between ML 0 and 4.2. We develop our analyses in three steps. First, we fit the Brune (1970) model to the nonparametric source spectra to estimate corner frequency and seismic moment, and we map the spatial distribution of stress drop across the Irpinia area. We found stress drops in the range 0.4–8.1 MPa, with earthquakes deeper than 7 km characterized by higher average stress drop (i.e., 3.2 MPa). Second, assuming a simple stress-release model (Kanamori and Heaton, 2000), we derive fracture energy and critical slip-weakening distance. The spatial variability of stress drop and fracture energy allows us to image the present stress conditions of fault segments activated during the 23 November 1980 Ms 6.9 earthquake. The variability of the source parameters shows clear patterns of the fault mechanical properties, suggesting that the Irpinia fault system can be divided into three main sectors, with the northern and southern ones showing different properties from the central one. Our results agree with previous studies indicating the presence of fluids with different composition in the different sectors of the Irpinia fault system. In the third step, we compare the time evolution of source parameters with a time series of geodetic displacement recorded near the fault system. Temporal trends in the correlation between geodetic displacement and different source parameters indicate that the poroelastic deformation perturbation generated by the karst aquifer recharge is modulating not only the occurrence rate of microseismicity (D’Agostino et al., 2018) but may lead to rupture asperities with different sizes and characteristics.

Dynamic Rupture Study of Near-Field Velocity Pulses during the 2016 Kumamoto Earthquake, Japan

ABSTRACT Simulating the ground motions of future earthquakes requires a proper understanding and modeling of source, path, and site effects. Ground motions recorded during recent earthquakes very close to their ruptured faults provide new evidence of the importance of source effects and suggest that physics-based rupture modeling is critical to account for them. Here, we develop dynamic rupture models to simulate the near-fault ground motions generated by the 2016 Kumamoto, Japan, earthquake (Mw 7.0) at Nishihara village, which feature a large-amplitude velocity pulse. Comparison of mainshock and foreshock waveforms suggests that the source of the velocity pulse is on the Futagawa fault segment located very close to the site. Our dynamic models use the spectral element method and are built upon a previous kinematic description of the event via a so-called “characterized source model,” with three strong-motion generation areas (SMGAs) on the assumed fault plane. We first develop a reference model that reproduces the main features of the rupture process in agreement with previous results of kinematic source inversion. We then examine the sensitivity of the simulated near-fault ground motions to the frictional parameters (critical slip-weakening distance and stress drop) in the shallow part of the fault and to the geometrical properties of the shallow SMGA. Even assuming drastically different frictional properties in the shallow part of the fault, the amplitude of the simulated ground motions was affected little. On the other hand, changes of geometrical properties of the shallow SMGA generated large differences in simulated ground motions. The results indicate that geometrical features of the shallow SMGA played a more important role in generating near-fault ground motions with velocity pulses as observed at Nishihara village during the 2016 Kumamoto earthquake.

Source Time Function Clustering Reveals Patterns in Earthquake Dynamics

AbstractWe cluster a global database of 3529 Mw>5.5 earthquakes in 1995–2018 based on a dynamic time warping distance between earthquake source time functions (STFs). The clustering exhibits different degrees of complexity of the STF shapes and suggests an association between STF complexity and earthquake source parameters. Most of the thrust events have simple STF shapes across all depths. In contrast, earthquakes with complex STF shapes tend to be located at shallow depths in complicated tectonic regions, exhibit long source duration compared with others of similar magnitude, and tend to have strike-slip mechanisms. With 2D dynamic modeling of dynamic ruptures on heterogeneous fault properties, we find a systematic variation of the simulated STF complexity with frictional properties. Comparison between the observed and synthetic clustering distributions provides useful constraints on frictional properties. In particular, the characteristic slip-weakening distance could be constrained to be short (<0.1 m) and depth dependent if stress drop is in general constant.

Application of a parameter-shifted grey wolf optimizer for earthquake dynamic rupture inversion

Global optimization is an essential approach to any inversion problem. Recently, the grey wolf optimizer (GWO) has been proposed to optimize the global minimum, which has been quickly used in a variety of inversion problems. In this study, we proposed a parameter-shifted grey wolf optimizer (psGWO) based on the conventional GWO algorithm to obtain the global minimum. Compared with GWO, the novel psGWO can effectively search targets toward objects without being trapped within the local minimum of the zero value. We confirmed the effectiveness of the new method in searching for uniform and random objectives by using mathematical functions released by the Congress on Evolutionary Computation. The psGWO algorithm was validated using up to 10,000 parameters to demonstrate its robustness in a large-scale optimization problem. We successfully applied psGWO in two-dimensional (2D) synthetic earthquake dynamic rupture inversion to obtain the frictional coefficients of the fault and critical slip-weakening distance using a homogeneous model. Furthermore, this algorithm was applied in inversions with heterogeneous distributions of dynamic rupture parameters. This implementation can be efficiently applied in 3D cases and even in actual earthquake inversion and would deepen the understanding of the physics of natural earthquakes in the future.

Parameter interdependence of dynamic self-similar crack with distance-weakening friction

SUMMARY In several analytical and numerical studies, the slip rate function and energy release rate for dynamic self-similar crack growth have been investigated, and the results obtained have contributed to a theoretical understanding and estimation of on-fault energetics. However, the relationships among physical parameters, including stress state, process zone size, rupture velocity, peak slip rate and energy release rate, are still unclear. Therefore, the aim of this study is to derive an analytical solution of the slip rate distribution of antiplane self-similar crack growth under distance-weakening friction that mimics slip-weakening friction. To satisfy the condition that the slip rate starts from zero at the rupture front, a trade-off relationship among rupture velocity, process zone size and breakdown stress-to-stress drop ratio is proposed. The peak slip rate, slip-weakening distance and fracture energy obtained using the proposed model provide a possible mechanism for the determination of the rupture velocity and the estimation of the fracture energy of the self-similar crack growth, based on the seismic observables.

Effects of seismogenic width and low-velocity zones on estimating slip-weakening distance from near-fault ground deformation

SUMMARYFault weakening process controls earthquake rupture propagation and is of great significance to impact the final earthquake size and seismic hazard. Critical slip-weakening distance (${D_c}$) is one of the key parameters, which however is of difficult endeavours to be determined on natural faults, mainly due to its strong trade-off with the fault strength drop. An estimation method of ${D_c}$ proposed by Fukuyama et al. provides a simple and direct reference of ${D_c}$ on real faults from the near-fault ground displacement at the peak of ground velocity (${D_c}^{\prime\prime}$). However, multiple factors may affect the observed near-fault ground velocity and thus need to be considered when estimating ${D_c}.$ In this work we conduct 3-D finite element numerical simulations to examine the effects of finite seismogenic width and near-fault low velocity zones (LVZs) on the results of ${D_c}^{\prime\prime}$. In uniform models with constant prescribed ${D_c}$, the derived ${D_c}^{\prime\prime}$ values increase with seismogenic width. Furthermore, the scaling between ${D_c}^{\prime\prime}$ and final slip in models with a constant ${D_c}$ indicates that the scale-dependent feature of ${D_c}^{\prime\prime}$ might not be related to variation in friction properties. With a near-fault LVZ, ${D_c}^{\prime\prime}$ values show significant magnification. The width of the LVZ plays a more important role in enlarging ${D_c}$ estimation compared to the depth of the LVZ. Complex wavefields and multiple wiggles introduced by the LVZ could lead to delay pick and then cause large deviation. The value of ${D_c}$ on the fault may be overestimated through ${D_c}^{\prime\prime}$ from limited stations only.

Effects of Initial Stress Field and Critical Slip-Weakening Distance on Rupture Selectivity of 3D Buried Branched Faults

ABSTRACTUsing the boundary integral equation method with the slip-weakening friction law, we investigate the effects of the initial stress field and the critical slip-weakening distance on the rupture selectivity on the 3D buried branched faults. The numerical results show that after reaching the bifurcation line between the main fault and the branched faults, rupture continues to propagate on one or both bifurcation planes (BPs), with rupture on one plane more favorable than the other. The initial stress distribution plays a decisive role in the selection of the favorable BP (FBP), and there is a critical status of stress distribution, around which rupture propagates on both planes, whereas the FBP switches between the two. For a given fault geometry, the critical status of initial stress, which is described by a ratio Fp between normal stresses, is related to the critical slip-weakening distance Dc.

Dynamic Modelling of Induced Seismicity by Using Seismic Efficiency Constraints and a New Scaling Law for Slip-Weakening Distance

In the numerical simulation of induced seismicity, much attention is generally paid to the calibration of the frictional resistance of the causative fault to obtain a seismic moment consistent with that of the actual event, whereas sufficient investigation is not made in the estimation of the slip-weakening distance Dc as well as in the calibration of seismically radiated energy. The present study addresses this problem by numerically and analytically investigating the relation between Dc and seismic source parameters. First, this study performs the dynamic simulation of an induced seismic event caused by a decrease in the effective normal stress. The analysis demonstrated that seismic efficiency η can be used to improve the accuracy of estimating the critical slip-weakening distance and the coefficient of kinetic friction µd whilst considering not only seismic moment but also radiated energy in the calibration. This gave insight into the development of the new calibration method for induced seismicity that considers energy-related seismic source parameters. Furthermore, a new scaling law of the slip-weakening distance was derived from the theoretical expression of seismic efficiency η, considering seismic moment Mo and scaled energy $$ \hat{e} $$. The proposed scaling law can yield the relation between Dc and Mo, which is shown to be similar to that obtained from a previous study, but additionally considers the relation between seismically radiated energy and Dc. The dependency of Dc on seismically radiated energy implied from the proposed scaling law has been verified from the dynamic analyses where η = 0.06 was used to place a constraint on Dc for seismic events with different magnitudes. The developed numerical simulation methodology of induced seismicity as well as the scaling law considering the energy indices significantly contributes to improving the accuracy of back-analysis, thus leading to a more accurate estimation of the mechanical properties of faults and/or shear zones in seismically active regions of deep underground mines or reservoirs composed of discontinuous hard rock masses.

Effect of Slip-Weakening Distance on Seismic\u2013Aseismic Slip Patterns

In order to explain the seismic–aseismic slip patterns observed on megathrust faults, numerical simulations were carried out using the quasi-dynamic asperity fault model with the slip-dependent friction and stress-dependent healing. Two friction law parameters, strength and slip-weakening distance, are interpreted in the subduction channel context. The parameters are treated as random fields with specified characteristics. Their distributions define heterogeneities of the interplate frictional coupling. The higher strength regions accumulate stresses, whereas the slip-weakening distance lengths control the stress release rates. The simulation results indicate that the slip-dependent asperity model reproduces key features of real megathrust fault behavior. First, stable and unstable slip movements can occur at the same locations, even if the friction parameters are fixed. Second, two rupture styles, single asperity breaks and wide, smooth, propagating rupture fronts, can be distinguished. The latter style is responsible for large slips near the free surface, where lower fault strengths are expected. The reason for these effects is that slip instabilities depend both on local friction and on the system stiffness, which is related to the slipping area size and distribution of slips. It is also shown that the high-strength interplate patches, such as subducted seamounts, can both promote and restrain large earthquakes, depending on the slip-weakening distance lengths.

Near-fault deformation and Dc″ during the 2016 Mw7.1 Kumamoto earthquake

An Mw7.1 Kumamoto earthquake occurred at 01:25:05 on April 16, 2016 (JST). The earthquake involved a rupture at a shallow depth along a strike-slip fault with surface breaks. Near-fault ground motion records, especially those of a strike-slip earthquake, can provide us with direct information on the earthquake source process. During the earthquake, near-fault seismograms were obtained at KMMH16 station located about 500 m off the fault. The ground displacements were well recovered from the double numerical integration of accelerograms at KMMH16 both on the surface and at the bottom of the 252-m-deep borehole. Fault-parallel static displacement was estimated to be about 1.1 m from the acceleration waveforms. The Dc″ value, which is defined as double the fault-parallel displacement at peak velocity time, was proposed as a proxy of the slip-weakening distance. Using both the velocity and displacement fault-parallel waveforms, the Dc″ value was estimated at about 1 m. This value was between 30 and 50% of the total slip on the fault, which is consistent with previous observations. Graphical abstract GRAPHICAL ABSTRACT CAPTION

On the scale dependence of earthquake stress drop

We discuss the debated issue of scale dependence in earthquake source mechanics with the goal of providing supporting evidence to foster the adoption of a coherent interpretative framework. We examine the heterogeneous distribution of source and constitutive parameters during individual ruptures and their scaling with earthquake size. We discuss evidence that slip, slip-weakening distance and breakdown work scale with seismic moment and are interpreted as scale dependent parameters. We integrate our estimates of earthquake stress drop, computed through a pseudo-dynamic approach, with many others available in the literature for both point sources and finite fault models. We obtain a picture of the earthquake stress drop scaling with seismic moment over an exceptional broad range of earthquake sizes (−8 < MW < 9). Our results confirm that stress drop values are scattered over three order of magnitude and emphasize the lack of corroborating evidence that stress drop scales with seismic moment. We discuss these results in terms of scale invariance of stress drop with source dimension to analyse the interpretation of this outcome in terms of self-similarity. Geophysicists are presently unable to provide physical explanations of dynamic self-similarity relying on deterministic descriptions of micro-scale processes. We conclude that the interpretation of the self-similar behaviour of stress drop scaling is strongly model dependent. We emphasize that it relies on a geometric description of source heterogeneity through the statistical properties of initial stress or fault-surface topography, in which only the latter is constrained by observations.

Dynamic rupture activation of backthrust fault branching

We perform dynamic rupture simulations to investigate the possible reactivation of backthrust branches triggered by ruptures along a main thrust fault. Simulations with slip-weakening fault friction and uniform initial stress show that fast propagation speed or long propagation distance of the main rupture promotes reactivation of backthrust over a range of branch angles. The latter condition may occur separately from the former if rupture speed is limited by an increasing slip-weakening distance towards the junction direction. The results suggest a trade-off between the amplitude and duration of the dynamic stress near the main rupture front for backthrust reactivation. Termination of the main rupture by a barrier can provide enhanced loading amplitude and duration along a backthrust rooted near the barrier, facilitating its reactivation especially with a high frictional resistance. The free surface and depth-dependent initial stress can have several additional effects. The sign of the triggered motion along the backthrust can be reversed from thrust to normal if a deeply nucleated main rupture breaks the free surface, while it is preserved as thrust if the main rupture is terminated by a barrier at depth. The numerical results are discussed in relation to several recent megathrust earthquakes in Sumatra, Chile, and Japan, and related topics such as branch feedbacks to the main fault. The dynamic view on backthrust fault branching provided by the study fills a gap not covered by quasi-static models or observations. A specific examined case of antithetic fault branching may be useful for indicating a barrier-like behavior along the main fault.

Effect of slip-weakening distance on selected seismic source parameters of mining-induced fault-slip

Fault-slip caused by mining activities in underground mines could cause severe damage to nearby mine openings due to seismic waves arising from the fault-slip. A better understanding of seismic source parameters, such as seismically radiated energy and the maximum slip rate, is of pivotal importance in estimating the damage induced by the fault-slip. It is widely recognized that dynamic behaviour of fault-slip is largely affected by slip-weakening behaviour. In the present study, the effect of the slip-weakening behaviour on selected seismic source parameters of fault-slip is investigated using a mine-wide numerical model encompassing a fault running parallel to an orebody. A linear slip-weakening law is employed in a parametrical study with respect to characteristic slip-weakening distance, d0, assuming two types of fault-slip: (a) zonal fault-slip taking place across the fault; (b) local fault-slip taking place within limited areas of the fault, i.e., it is assumed that fault-slip is arrested by surrounding rockmasses. Results obtained from the dynamic analysis indicate that seismically radiated energy and slip rates are significantly affected by d0 for both types of fault-slip. On the other hand, the seismic moment is susceptible to d0 only for zonal fault-slip, given that d0 is smaller than the maximum shear displacement increment during the fault-slip. These results imply that considering slip-weakening distance is critical in estimating seismically radiated energy and slip rates, whereas slip-weakening behaviour might not have a large influence on seismic moment of local fault-slip.