Stress Drop during Earthquakes: Effect of Fault Roughness Scaling

2011 ◽  
Vol 101 (5) ◽  
pp. 2369-2387 ◽  
Author(s):  
T. Candela ◽  
F. Renard ◽  
M. Bouchon ◽  
J. Schmittbuhl ◽  
E. E. Brodsky
Keyword(s):  
2020 ◽  
Vol 224 (2) ◽  
pp. 1371-1380
Author(s):  
Aglaja Blanke ◽  
Grzegorz Kwiatek ◽  
Thomas H W Goebel ◽  
Marco Bohnhoff ◽  
Georg Dresen

SUMMARY Earthquake source parameters such as seismic stress drop and corner frequency are observed to vary widely, leading to persistent discussion on potential scaling of stress drop and event size. Physical mechanisms that govern stress drop variations are difficult to evaluate in nature and are more readily studied in controlled laboratory experiments. We perform two stick-slip experiments on fractured (rough) and cut (smooth) Westerly granite samples to explore fault roughness effects on acoustic emission (AE) source parameters. We separate large stick-slip events that generally saturate the seismic recording system from populations of smaller AE events which are sensitive to fault stresses prior to slip. AE event populations show many similarities to natural seismicity and may be interpreted as laboratory equivalent of natural microseismic events. We then compare the temporal evolution of mechanical data such as measured stress release during slip to temporal changes in stress drops derived from AEs using the spectral ratio technique. We report on two primary observations: (1) In contrast to most case studies for natural earthquakes, we observe a strong increase in seismic stress drop with AE size. (2) The scaling of stress drop with magnitude is governed by fault roughness, whereby the rough fault shows a more rapid increase of the stress drop–magnitude relation with progressing large stick-slip events than the smooth fault. The overall range of AE sizes on the rough surface is influenced by both the average grain size and the width of the fault core. The magnitudes of the smallest AE events on smooth faults may also be governed by grain size. However, AEs significantly grow beyond peak roughness and the width of the fault core. Our laboratory tests highlight that source parameters vary substantially in the presence of fault zone heterogeneity (i.e. roughness and narrow grain size distribution), which may affect seismic energy partitioning and static stress drops of small and large AE events.


2017 ◽  
Vol 44 (2) ◽  
pp. 777-783 ◽  
Author(s):  
O. Zielke ◽  
M. Galis ◽  
P. M. Mai

Author(s):  
Serge A. Shapiro ◽  
Carsten Dinske

AbstractSometimes, a rather high stress drop characterizes earthquakes induced by underground fluid injections or productions. In addition, long-term fluid operations in the underground can influence a seismogenic reaction of the rock per unit volume of the fluid involved. The seismogenic index is a quantitative characteristic of such a reaction. We derive a relationship between the seismogenic index and stress drop. This relationship shows that the seismogenic index increases with the average stress drop of induced seismicity. Further, we formulate a simple and rather general phenomenological model of stress drop of induced earthquakes. This model shows that both a decrease of fault cohesion during the earthquake rupture process and an enhanced level of effective stresses could lead to high stress drop. Using these two formulations, we propose the following mechanism of increasing induced seismicity rates observed, e.g., by long-term gas production at Groningen. Pore pressure depletion can lead to a systematic increase of the average stress drop (and thus, of magnitudes) due to gradually destabilizing cohesive faults and due to a general increase of effective stresses. Consequently, elevated average stress drop increases seismogenic index. This can lead to seismic risk increasing with the operation time of an underground reservoir.


1999 ◽  
Vol 122 (2) ◽  
pp. 203-208 ◽  
Author(s):  
Chuwei Zhou ◽  
Wei Yang ◽  
Daining Fang

Mechanical properties and damage evolution of short-fiber-reinforced metal matrix composites (MMC) are studied under a micromechanics model accounting for the history of cooling and thermal cycling. A cohesive interface is formulated in conjunction with the Gurson-Tvergaard matrix damage model. Attention is focused on the residual stresses and damages by the thermal mismatch. Substantial stress drop in the uniaxial tensile response is found for a computational cell that experienced a cooling process. The stress drop is caused by debonding along the fiber ends. Subsequent thermal cycling lowers the debonding stress and the debonding strain. Micromechanics analysis reveals three failure modes. When the thermal histories are ignored, the cell fails by matrix damage outside the fiber ends. With the incorporation of cooling, the cell fails by fiber end debonding and the subsequent transverse matrix damage. When thermal cycling is also included, the cell fails by jagged debonding around the fiber tops followed by necking instability of matrix ligaments. [S0094-4289(00)01202-0]


2021 ◽  
Author(s):  
Itzhak Lior ◽  
Anthony Sladen ◽  
Diego Mercerat ◽  
Jean-Paul Ampuero ◽  
Diane Rivet ◽  
...  

<p>The use of Distributed Acoustic Sensing (DAS) presents unique advantages for earthquake monitoring compared with standard seismic networks: spatially dense measurements adapted for harsh environments and designed for remote operation. However, the ability to determine earthquake source parameters using DAS is yet to be fully established. In particular, resolving the magnitude and stress drop, is a fundamental objective for seismic monitoring and earthquake early warning. To apply existing methods for source parameter estimation to DAS signals, they must first be converted from strain to ground motions. This conversion can be achieved using the waves’ apparent phase velocity, which varies for different seismic phases ranging from fast body-waves to slow surface- and scattered-waves. To facilitate this conversion and improve its reliability, an algorithm for slowness determination is presented, based on the local slant-stack transform. This approach yields a unique slowness value at each time instance of a DAS time-series. The ability to convert strain-rate signals to ground accelerations is validated using simulated data and applied to several earthquakes recorded by dark fibers of three ocean-bottom telecommunication cables in the Mediterranean Sea. The conversion emphasizes fast body-waves compared to slow scattered-waves and ambient noise, and is robust even in the presence of correlated noise and varying wave propagation directions. Good agreement is found between source parameters determined using converted DAS waveforms and on-land seismometers for both P- and S-wave records. The demonstrated ability to resolve source parameters using P-waves on horizontal ocean-bottom fibers is key for the implementation of DAS based earthquake early warning, which will significantly improve hazard mitigation capabilities for offshore and tsunami earthquakes.</p>


2021 ◽  
Author(s):  
Jagdish Chandra Vyas ◽  
Martin Galis ◽  
Paul Martin Mai

<p>Geological observations show variations in fault-surface topography not only at large scale (segmentation) but also at small scale (roughness). These geometrical complexities strongly affect the stress distribution and frictional strength of the fault, and therefore control the earthquake rupture process and resulting ground-shaking. Previous studies examined fault-segmentation effects on ground-shaking, but our understanding of fault-roughness effects on seismic wavefield radiation and earthquake ground-motion is still limited.  </p><p>In this study we examine the effects of fault roughness on ground-shaking variability as a function of distance based on 3D dynamic rupture simulations. We consider linear slip-weakening friction, variations of fault-roughness parametrizations, and alternative nucleation positions (unilateral and bilateral ruptures). We use generalized finite difference method to compute synthetic waveforms (max. resolved frequency 5.75 Hz) at numerous surface sites  to carry out statistical analysis.  </p><p>Our simulations reveal that ground-motion variability from unilateral ruptures is almost independent of  distance from the fault, with comparable or higher values than estimates from ground-motion prediction equations (e.g., Boore and Atkinson, 2008; Campbell and Bozornia, 2008). However, ground-motion variability from bilateral ruptures decreases with increasing distance, in contrast to previous studies (e.g., Imtiaz et. al., 2015) who observe an increasing trend with distance. Ground-shaking variability from unilateral ruptures is higher than for bilateral ruptures, a feature due to intricate seismic radiation patterns related to fault roughness and hypocenter location. Moreover, ground-shaking variability for rougher faults is lower than for smoother faults. As fault roughness increases the difference in ground-shaking variabilities between unilateral and bilateral ruptures increases. In summary, our simulations help develop a fundamental understanding of ground-motion variability at high frequencies (~ 6 Hz) due small-scale geometrical fault-surface variations.</p>


1990 ◽  
Vol 80 (6A) ◽  
pp. 1553-1570 ◽  
Author(s):  
R. B. Horner ◽  
R. J. Wetmiller ◽  
M. Lamontagne ◽  
M. Plouffe

Abstract Relative locations of 323 large aftershocks (M 3.0 or greater) in the period from 5 October 1985 to 25 March 1988 show that the Ms 6.6 event on 5 October 1985 initiated at 62.208°N, 124.217°W, about 2.5 km northeast of the Ms 6.9 main shock on 23 December 1985. The overall aftershock distribution suggests the October rupture was primarily a west-dipping, low-angle thrust. In subsequent aftershock activity, the main rupture plane was marked by a distinct quiescent area of about 200 km2 that persisted until the 23 December event. Most of the stress drop and slip occurred in this area. Following the 23 December rupture, a similar sized quiescent zone was also observed; however, it was only evident during the first 24 hr of the aftershock sequence, and the area was about 50 per cent too small to yield the overall stress drop. The additional area appeared to come from secondary rupture zones that developed coincident with the main shock rupture. Precise locations of 182 small (M 3.0 or less) aftershocks recorded during a third field survey from 12 to 21 September 1986 indicated at least one and probably three high-angle faults. Composite mechanism solutions showed thrust faulting except in a region directly south of the main shock rupture areas where there is a bend in one of the secondary fault zones and a concentration of aftershock activity. Mechanism solutions calculated for five of the largest aftershocks in the same region also indicated a similar variability. Development of secondary fault zones explained the increased complexity of the December event and may also provide an explanation for the vertical peak acceleration exceeding 2 g that was recorded about 10 sec after the December rupture initiated.


1983 ◽  
Vol 73 (3) ◽  
pp. 853-862
Author(s):  
J. Lomnitz-Adler

abstract A model is presented for the simulation of the statistical features of the earthquake process. An analytical solution is given for a simple case, and a numerical calculation of spatial and stress drop frequency distribution has been carried out. Extensions of the model are discussed.


Sign in / Sign up

Export Citation Format

Share Document