An official website of the United States government
Abstract Shallow slow-slip events (SSEs) contribute to strain release near the shallow portions of subduction interfaces and may contribute to promoting shallow subduction earthquakes. Recent efforts in offshore monitoring of shallow SSEs have provided evidence of possible interactions between shallow SSEs and megathrust earthquakes. In this study, we use a dynamic earthquake simulator that captures both quasi-static (for SSEs) and dynamic (for megathrust earthquakes) slip to explore their interactions and implications for seismic and tsunami hazards. We model slip behaviors of a shallow-dipping subduction interface on which two locally locked patches (asperities) with different strengths are embedded within a conditionally stable zone. We find that both SSEs and earthquakes can occur, and they interact over multiple earthquake cycles in the model. Dynamic ruptures can nucleate on the asperities and propagate into the surrounding conditionally stable zone at slow speeds, generating tsunami earthquakes. A clear correlation emerges between the size of an earthquake and SSE activities preceding it. Small earthquakes rupture only the low-strength asperity, whereas large earthquakes rupture both. Before a large earthquake, periodic SSEs occur around the high-strength asperity, gradually loading stress into its interior. The critically stressed high-strength asperity can be ruptured together with the low-strength one in the large earthquake, followed by a relatively quiet interseismic period with very few SSEs and then a small earthquake. An SSE may or may not directly lead to nucleation of an earthquake, depending on whether a nearby asperity is ready for spontaneously dynamic failure. In addition, because of different SSE activities, the coupling degree may change dramatically between different interseismic periods, suggesting that its estimate based on a short period of observation may be biased.
more »
« less
Using a dynamic earthquake simulator to explore tsunami earthquake generation
SUMMARY Observations of historical tsunami earthquakes reveal that ruptures of these earthquakes propagate slowly at shallow depth with longer duration, depletion in high-frequency radiation and larger discrepancy of Mw–Ms than ordinary megathrust earthquakes. They can effectively generate tsunami and lead to huge damage to regional populated areas near the coast. In this study, we use a recently developed dynamic earthquake simulator to explore tsunami earthquake generation from a physics-based modelling point of view. We build a shallow-dipping subduction zone model in which locally locked, unstable patches (asperities) are distributed on a conditionally stable subduction interface at shallow depth. The dynamic earthquake simulator captures both quasi-static and dynamic processes of earthquake cycles. We find that earthquakes can nucleate on these asperities and propagate into the surrounding conditionally stable zone at slow speeds, generating tsunami earthquakes. A high normal stress asperity, representing a subducted seamount, can act as an asperity in some events but as a barrier in other events over multiple earthquake cycles. Low normal stress asperities typically act as asperities in tsunami earthquakes. The degree of velocity-weakening in the conditionally stable zone, which may sustain rupture at different speeds or stop rupture, is critical for tsunami earthquake generation and affects its recurrence interval. Distributed asperities may rupture in isolated events separated by tens of years, or in a sequence of events separated by hours to days, or in one large event in a cascade fashion, demonstrating complex interactions among them. The recurrence interval on a high normal stress asperity is much larger than that on low normal stress asperities. These modelling results shed lights on the observations from historical tsunami earthquakes, including the 1994 and 2006 Java tsunami earthquakes and 2010 Mentawai tsunami earthquake.
more »
« less
- Award ID(s):
- 2013695
- PAR ID:
- 10380711
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 229
- Issue:
- 1
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 255-273
- Size(s):
- p. 255-273
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The Mediterranean Hellenic Arc subduction zone (HASZ) has generated several 8 earthquakes and tsunamis. Seismic‐probabilistic tsunami hazard assessment typically utilizes uniform or stochastic earthquake models, which may not represent dynamic rupture and tsunami generation complexity. We present an ensemble of ten 3D dynamic rupture earthquake scenarios for the HASZ, utilizing a realistic slab geometry. Our simplest models use uniform along‐arc pre‐stresses or a single circular initial stress asperity. We then introduce progressively more complex models varying initial shear stress along‐arc, multiple asperities based on scale‐dependent critical slip weakening distance, and a most complex model blending all aforementioned heterogeneities. Thereby, regional initial conditions are constrained without relying on detailed geodetic locking models. Varying epicentral locations in the simplest, homogeneous model leads to different rupture speeds and moment magnitudes. We observe dynamic fault slip penetrating the shallow slip‐strengthening region and affecting seafloor uplift. Off‐fault plastic deformation can double vertical seafloor uplift. A single‐asperity model generates a 8 scenario resembling the 1303 Crete earthquake. Using along‐strike varying initial stresses results in 8.0–8.5 dynamic rupture scenarios with diverse slip rates and uplift patterns. The model with the most heterogeneous initial conditions yields a 7.5 scenario. Dynamic rupture complexity in prestress and fracture energy tends to lower earthquake magnitude but enhances tsunamigenic displacements. Our results offer insights into the dynamics of potential large Hellenic Arc megathrust earthquakes and may inform future physics‐based joint seismic and tsunami hazard assessments.more » « less
-
The Mediterranean Hellenic Arc subduction zone (HASZ) has generatedseveral Mw>=8 earthquakes and tsunamis.Seismic-probabilistic tsunami hazard assessment typically utilizesuniform or stochastic earthquake models, which may not represent dynamicrupture and tsunami generation complexity. We present an ensemble of ten3D dynamic rupture earthquake scenarios for the HASZ, utilizing arealistic slab geometry. Our simplest models use uniform along-arcpre-stresses or a single circular initial stress asperity. We thenintroduce progressively more complex models varying initial shear stressalong-arc, multiple asperities based on scale-dependent critical slipweakening distance, and a most complex model blending all aforementionedheterogeneities. Thereby, regional initial conditions are constrainedwithout relying on detailed geodetic locking models. Varying hypocenterlocations in the simplest, homogeneous model leads to different rupturespeeds and moment magnitudes. We observe dynamic fault slip penetratingthe shallow slip-strengthening region and affecting seafloor uplift.Off-fault plastic deformation can double vertical seafloor uplift. Asingle-asperity model generates a Mw~8 scenarioresembling the 1303 Crete earthquake. Using along-strike varying initialstresses results in Mw~8.0-8.5 dynamic rupture scenarioswith diverse slip rates and uplift patterns. The model with the mostheterogeneous initial conditions yields a Mw~7.5scenario. Dynamic rupture complexity in prestress and fracture energytends to lower earthquake magnitude but enhances tsunamigenicdisplacements. Our results offer insights into the dynamics of potentiallarge Hellenic Arc megathrust earthquakes and may inform futurephysics-based joint seismic and tsunami hazard assessments.more » « less
-
Faults are usually surrounded by damage zones associated with localized deformation. Here we use fully dynamic earthquake cycle simulations to quantify the behaviors of earthquakes in fault damage zones. We show that fault damage zones can make a significant contribution to the spatial and temporal seismicity distribution. Fault stress heterogeneities generated by fault zone waves persist over multiple earthquake cycles that, in turn, produce small earthquakes that are absent in homogeneous simulations with the same friction conditions. Shallow fault zones can produce a bimodal depth distribution of earthquakes with clustering of seismicity at both shallower and deeper depths. Fault zone healing during the interseismic period also promotes the penetration of aseismic slip into the locked region and reduces the sizes of fault asperities that host earthquakes. Hence, small and moderate subsurface earthquakes with irregular recurrence intervals are commonly observed in immature fault zone simulations with interseismic healing. To link our simulation results to geological observations, we will use simulated fault slip at different depths to infer the timing and recurrence intervals of earthquakes and discuss how such measurements can affect our understanding of earthquake behaviors. We will also show that the maturity and material properties of fault damage zones have strong influence on whether long-term earthquake characteristics are represented by single events.more » « less
-
Abstract The additional work of ploughing makes seamounts more resistant to subduction and more strongly coupled than smoother areas. Nevertheless, the idea that subducted seamounts are weakly coupled and slip aseismically has become dominant in the last decade. This idea is primarily based on the claim that a seamount being subducted in the southern Japan Trench behaves this way. The key element in this assertion is that largeM ∼ 7 earthquakes that abut the leading edge of the seamount require that the seamount be aseismically sliding to initiate them. More recent observations show instead that the surrounding region is aseismically sliding while the seamount acts as a stationary buttress. Here we re‐examine this case and model it with both weak and strong asperity assumptions. Our modeling results show that only a strong asperity model can produce this type of earthquake. Strong asperities also rupture the seamount in great earthquakes with long recurrence times. This provides the previously unknown source for a series of great tsunami earthquakes that have occurred along the southern Japan Trench, the most recent being the 1677 M8.3–8.6 Enpo Boso‐oki tsunami earthquake. The “weak asperity” hypothesis is thus found to be false in this foundational example.more » « less
