An official website of the United States government
Abstract Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably‐dipping planar and variably‐curved nonplanar megathrusts using the volumetric, high‐order accurate codetandemto account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly‐dipping megathrusts to host larger earthquakes than steeply‐dipping ones. Under elevated pore pressure and less strongly velocity‐weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly‐dipping and sharply‐curved megathrusts host multi‐period supercycles of slow‐to‐fast, small‐to‐large slip events under higher effective stresses and more strongly velocity‐weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second‐order effects from structurally‐derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co‐ and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat faults.
more »
« less
Sensitivity Analysis of the Thermal Structure Within Subduction Zones Using Reduced‐Order Modeling
Abstract Megathrust earthquakes are the largest on Earth, capable of causing strong ground shaking and generating tsunamis. Physical models used to understand megathrust earthquake hazard are limited by existing uncertainties about material properties and governing processes in subduction zones. A key quantity in megathrust hazard assessment is the distance between the updip and downdip rupture limits. The thermal structure of a subduction zone exerts a first‐order control on the extent of rupture. We simulate temperature for profiles of the Cascadia, Nankai and Hikurangi subduction zones using a 2D coupled kinematic‐dynamic thermal model. We then build reduced‐order models (ROMs) for temperature using the interpolated Proper Orthogonal Decomposition (iPOD). The resulting ROMs are data‐driven, model agnostic, and computationally cheap to evaluate. Using the ROMs, we can efficiently investigate the sensitivity of temperature to input parameters, physical processes, and modeling choices. We find that temperature, and by extension the potential rupture extent, is most sensitive to variability in parameters that describe shear heating on the slab interface, followed by parameters controlling the thermal structure of the incoming lithosphere and coupling between the slab and the mantle. We quantify the effect of using steady‐state versus time‐dependent models, and of uncertainty in the choice of isotherm representing the downdip rupture limit. We show that variability in input parameters translates to significant differences in estimated moment magnitude. Our analysis highlights the strong effect of variability in the apparent coefficient of friction, with previously published ranges resulting in pronounced variability in estimated rupture limit depths.
more »
« less
- PAR ID:
- 10632619
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 26
- Issue:
- 5
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In subduction zones, along‐strike and downdip variations in megathrust slip behavior are linked to changes in the properties of the subducting and overriding plates. Although marine geophysical methods provide insights into subduction zone structures, most surveys consist of sparse 2D profiles, limiting our understanding of first‐order controls. Here, we use active‐source seismic data to derive a 3D crustal‐scale P‐wave velocity model of the Alaska Peninsula subduction zone that encompasses both plates and spans the Semidi segment and SW Kodiak asperity. Our results reveal modest variations within the incoming plate, attributed to a series of fracture zones, seamounts and their associated basement swell, collectively contributing to plate hydration. Basement swell appears to modulate the distribution and type of sediment entering the trench, likely impacting observed variations in slip behavior. The overriding plate exhibits significant heterogeneity. The updip limit and width of the dynamic backstop are similar between the SW Kodiak asperity and eastern Semidi segment but differ significantly from the Western Semidi segment. These distinctions likely account for differences in earthquake rupture patterns and interseismic coupling among these segments. Additionally, high velocities in the mid‐lower forearc crust coincide with the location of the megathrust slip during the Mw 8.2 2021 Chignik event. We interpret these velocities as intracrustal intrusions that contributed to the deep rupture of the 2021 event. Our findings suggest that the contrasting structural and material properties of both the incoming and overriding plates influence the spatially complex and semi‐persistent segmentation of the megathrust offshore the Alaska Peninsula.more » « less
-
Abstract Great subduction earthquakes (Mw ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward‐ or landward‐vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward‐vergent splays host more slip than those with steeper, landward‐vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip‐ and vergence‐dependence of splay rupture: static (mis)alignment relative to far‐field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip‐propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward‐vergent splays and inhibit backward branching onto landward‐vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.more » « less
-
Tsunamigenic megathrust earthquakes along the Cascadia subduction zone present a major hazard concern. We can better prepare to model the earthquake source in a rapid manner by imbuing fault geometry constraints based on prior knowledge and by evaluating the capabilities of using existing GNSS sensors. Near-field GNSS waveforms have shown promise in providing rapid coarse finite-fault model approximations of the earthquake rupture that can improve tsunami modeling and response time. In this study, we explore the performance of GNSS derived finite-fault inversions and tsunami forecasting predictions in Cascadia that highlights the impact and potential of geodetic techniques and data in operational earthquake and tsunami monitoring. We utilized 1300 Cascadia earthquake simulations (FakeQuakes) that provide realistic (M7.5-9.3) rupture scenarios to assess how feasibly finite-fault models can be obtained in a rapid earthquake early warning and tsunami response context. A series of fault models with rectangular dislocation patches spanning the Cascadia megathrust area is added to the GFAST inversion algorithm to calculate slip for each earthquake scenario. Another method used to constrain the finite-fault geometry is from the GNSS-derived CMT fault plane solution. For the Cascadia region, we show that fault discretization using two rectangular segments approximating the megathrust portion of the subduction zone leads to improvements in modeling magnitude, fault slip, tsunami amplitude, and inundation. In relation to tsunami forecasting capabilities, we compare coastal amplitude predictions spanning from Vancouver Island (Canada) to Northern California (USA). Generally, the coastal amplitudes derived using fault parameters from the CMT solutions show an overestimation bias compared to amplitudes derived from the fixed slab model. We also see improved prediction values of the run-up height and maximum amplitude at 10 tide gauge stations using the fixed slab model as well.more » « less
-
Hodges, K (Ed.)We develop a linear viscous constitutive relationship for pressure solution constrained by models of deformed metasedimentary rocks and observations of exposed rocks from ancient subduction zones. We include pressure and temperature dependence on the solubility of silica in fluid by parameterizing a practical van’t Hoff relationship. This general flow law is well suited for making predictions about interseismic behavior of subduction zones. We apply the flow law to Cascadia, where thermal structure, geometry, relative plate velocity, and Global Positioning System velocity field are well constrained. Results are consistent with the temperature conditions at which resolvable ductile strain is recorded in subducted mudstones (at depths near the updip limit of the seismogenic zone) and with relative plate motion accommodated completely by viscous deformation (at depths near the downdip limit of the seismogenic zone). The flow law also predicts the observed forearc tapering of slip rate deficit with depth.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1029/2024GC011937
