scholarly journals Testing Dynamic Earthquake Rupture Models Generated With Stochastic Stress Drop

2018 ◽  
Author(s):  
A Pitarka
Author(s):  
Percy Galvez ◽  
Anatoly Petukhin ◽  
Paul Somerville ◽  
Jean-Paul Ampuero ◽  
Ken Miyakoshi ◽  
...  

ABSTRACT Realistic dynamic rupture modeling validated by observed earthquakes is necessary for estimating parameters that are poorly resolved by seismic source inversion, such as stress drop, rupture velocity, and slip rate function. Source inversions using forward dynamic modeling are increasingly used to obtain earthquake rupture models. In this study, to generate a large number of physically self-consistent rupture models, rupture process of which is consistent with the spatiotemporal heterogeneity of stress produced by previous earthquakes on the same fault, we use multicycle simulations under the rate and state (RS) friction law. We adopt a one-way coupling from multicycle simulations to dynamic rupture simulations; the quasidynamic solver QDYN is used to nucleate the seismic events and the spectral element dynamic solver SPECFEM3D to resolve their rupture process. To simulate realistic seismicity, with a wide range of magnitudes and irregular recurrence, several realizations of 2D-correlated heterogeneous random distributions of characteristic weakening distance (Dc) in RS friction are tested. Other important parameters are the normal stress, which controls the stress drop and rupture velocity during an earthquake, and the maximum value of Dc, which controls rupture velocity but not stress drop. We perform a parametric study on a vertical planar fault and generate a set of a hundred spontaneous rupture models in a wide magnitude range (Mw 5.5–7.4). We validate the rupture models by comparison of source scaling, ground motion (GM), and surface slip properties to observations. We compare the source-scaling relations between rupture area, average slip, and seismic moment of the modeled events with empirical ones derived from source inversions. Near-fault GMs are computed from the source models. Their peak ground velocities and peak ground accelerations agree well with the ground-motion prediction equation values. We also obtain good agreement of the surface fault displacements with observed values.


2018 ◽  
Vol 733 ◽  
pp. 148-158 ◽  
Author(s):  
Agnès Chounet ◽  
Martin Vallée ◽  
Mathieu Causse ◽  
Françoise Courboulex

2020 ◽  
Vol 224 (1) ◽  
pp. 487-516 ◽  
Author(s):  
E H Madden ◽  
M Bader ◽  
J Behrens ◽  
Y van Dinther ◽  
A-A Gabriel ◽  
...  

SUMMARY How does megathrust earthquake rupture govern tsunami behaviour? Recent modelling advances permit evaluation of the influence of 3-D earthquake dynamics on tsunami genesis, propagation, and coastal inundation. Here, we present and explore a virtual laboratory in which the tsunami source arises from 3-D coseismic seafloor displacements generated by a dynamic earthquake rupture model. This is achieved by linking open-source earthquake and tsunami computational models that follow discontinuous Galerkin schemes and are facilitated by highly optimized parallel algorithms and software. We present three scenarios demonstrating the flexibility and capabilities of linked modelling. In the first two scenarios, we use a dynamic earthquake source including time-dependent spontaneous failure along a 3-D planar fault surrounded by homogeneous rock and depth-dependent, near-lithostatic stresses. We investigate how slip to the trench influences tsunami behaviour by simulating one blind and one surface-breaching rupture. The blind rupture scenario exhibits distinct earthquake characteristics (lower slip, shorter rupture duration, lower stress drop, lower rupture speed), but the tsunami is similar to that from the surface-breaching rupture in run-up and length of impacted coastline. The higher tsunami-generating efficiency of the blind rupture may explain how there are differences in earthquake characteristics between the scenarios, but similarities in tsunami inundation patterns. However, the lower seafloor displacements in the blind rupture result in a smaller displaced volume of water leading to a narrower inundation corridor inland from the coast and a 15 per cent smaller inundation area overall. In the third scenario, the 3-D earthquake model is initialized using a seismo-thermo-mechanical geodynamic model simulating both subduction dynamics and seismic cycles. This ensures that the curved fault geometry, heterogeneous stresses and strength and material structure are consistent with each other and with millions of years of modelled deformation in the subduction channel. These conditions lead to a realistic rupture in terms of velocity and stress drop that is blind, but efficiently generates a tsunami. In all scenarios, comparison with the tsunamis sourced by the time-dependent seafloor displacements, using only the time-independent displacements alters tsunami temporal behaviour, resulting in later tsunami arrival at the coast, but faster coastal inundation. In the scenarios with the surface-breaching and subduction-initialized earthquakes, using the time-independent displacements also overpredicts run-up. In the future, the here presented scenarios may be useful for comparison of alternative dynamic earthquake-tsunami modelling approaches or linking choices, and can be readily developed into more complex applications to study how earthquake source dynamics influence tsunami genesis, propagation and inundation.


2019 ◽  
Vol 219 (2) ◽  
pp. 1237-1252
Author(s):  
J S Neely ◽  
Y Huang ◽  
W Fan

SUMMARY The 280-km-long San Cristobal Trough (SCT), created by the tearing of the Australia plate as it subducts under the Pacific Plate near the Solomon and Vanuatu subduction zones, has hosted strike-slip earthquake sequences in 1993 and 2015. Both sequences, which likely represent a complete seismic cycle, began along the oldest section of the SCT—the portion farthest from the tear that has experienced the most cumulative displacement—and migrated to the younger sections closer to the tear. The SCT's abundant seismicity allows us to study transform boundary development—a process rarely observed along a single fault system—through observations of earthquake rupture properties. Using the spectral ratio method based on empirical Green's functions (EGFs), we calculate the corner frequencies of three Mw ∼7 2015 earthquakes and colocated smaller earthquakes. We utilize two different spectral ratio stacking methods and fit both Brune and Boatwright models to the stacked spectral ratios. Regardless of stacking methods and spectral models, we find that the corner frequencies of the 2015 Mw ∼7 earthquakes decrease slightly with distance from the tear. Assuming a constant rupture velocity and an omega-square spectral model, this corner frequency decrease may be due to an increase in rupture length with distance from the tear. The spectrum of the 2015 earthquake farthest from the tear also deviates from the omega-square model, which may indicate rupture complexity. Stress drop estimates from the corner frequencies of the 2015 Mw ∼7 earthquakes range between 1 and 7 MPa, whereas stress drop estimates of their EGFs range from ∼0.05 to 10 MPa with most values between 0.1 and 1 MPa. Independent evidence from a second moments analysis of the 2015 earthquake sequence also indicates a possible increase in rupture length with distance from the tear, confirming the results from the spectral ratio analysis. We also observe an increase in normalized centroid time-delay values, a first-order proxy for rupture behaviour, with distance from the tear for the 2015 sequence. A similar trend for the 1993 sequence suggests that earthquake rupture varies systematically along the SCT. Since distance from the tear corresponds to cumulative fault displacement, these along-strike rupture variations may be due to a displacement-driven fault maturation process.


2016 ◽  
Author(s):  
Christie D. Rowe ◽  
◽  
W. Ashley Griffith ◽  
Catherine Ross ◽  
Benjamin Melosh ◽  
...  
Keyword(s):  

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.


Sign in / Sign up

Export Citation Format

Share Document