skip to main content


Title: Relationships Among Forearc Structure, Fault Slip, and Earthquake Magnitude: Numerical Simulations With Applications to the Central Chilean Margin
Abstract

Two adjacent segments of the Chile margin exhibit significant differences in earthquake magnitude and rupture extents during the 1960 Valdivia and 2010 Maule earthquakes. We use the discrete element method (DEM) to simulate the upper plate as having an inner and outer wedge defined by different frictional domains along the décollement. We find that outer wedge width strongly influences coseismic slip distributions. We use the published peak slip magnitudes to pick best fit slip distributions and compare our models to geophysical constraints on outer wedge widths for the margins. We obtain reasonable fits to published slip distributions for the 2010 Maule rupture. Our best‐fit slip distribution for the 1960 Valdivia earthquake suggests that peak slip occurred close to the trench, differing from published models but being supported by new seismic interpretations along this margin. Finally, we also demonstrate that frictional conditions beneath the outer wedge can affect the coseismic slip distributions.

 
more » « less
Award ID(s):
1723249
NSF-PAR ID:
10359900
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
48
Issue:
13
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We study the deep structure of the continental wedge along the south‐central Chilean convergent margin based on a joint interpretation of wide‐angle reflection and refraction seismic phases acquired between the southern part of the Maule 2010 (36°S) and the northern part of the Valdivia 1960 (41°S) earthquake rupture areas. When combined with results of previous seismic studies, our results provide new insights into the deep structure (to 15–20 km) of the overriding, continental South American plate. We observe a latitudinal variation ofPwave velocity, interpreted in terms of the distribution of Permo‐Triassic metamorphic complex near the coast, and deep reflectors at the base of these metamorphic units. A change in the observed deep structural style near 38°S suggests a segmentation in the interplate frictional properties between the Maule and Valdivia earthquakes rupture zones. This regional change in the continental wedge structure indicates a spatial correlation between the coseismic ruptures of large earthquakes (e.g., the 1960 and 2010 megathrust earthquakes) and the rheology/lithology along the interplate boundary.

     
    more » « less
  2. Abstract

    The south central Chilean margin regularly produces many of the world's largest earthquakes and tsunami, including the 2010Mw8.8 Maule and 1960Mw9.5 Valdivia events. In 2017, we acquired seismic reflection data along ~1,000 km of the margin using the R/VLangseth's 15 km long receiver array and 108.2 l (6,600 in3) seismic source to image structures associated with these ruptures. We focus on the Valdivia segment with the largest coseismic slip (~40 m). The outer 40 km of the forearc is an accretionary wedge constructed primarily of stacked sedimentary packages with irregular lengths and thicknesses and little along‐strike continuity. Forearc structures indicate that the accretionary wedge grows primarily through basal accretion of the downgoing trench fill. The décollement propagates along a weak boundary near the top of the trench fill but occasionally branches downward into the underthrust sediment along bedding horizons, peeling off slices that are underplated to the forearc. The shallow décollement level and the rarity of underplating events allow most of the trench sediment to subduct. As a result, only ~30% of the incoming sediment has been accreted since the Early Pliocene. This implies that, on average, ~1 km of sediment must subduct beyond the outer forearc, an inference that is supported by our seismic images. We propose that the thickness and great downdip and along‐strike extent of the underthrust layer, which separates the megathrust from the underlying roughness of the igneous ocean crust, ensures a smooth broad zone of strong coupling that generates the world's largest earthquakes and tsunami.

     
    more » « less
  3. Abstract

    The Shumagin seismic gap along the Alaska Peninsula experienced a major,MW7.8, interplate thrust earthquake on 22 July 2020. Several available finite‐fault inversions indicate patchy slip of up to 4 m at 8–48 km depth. There are differences among the models in peak slip and absolute placement of slip on the plate boundary, resulting from differences in data distributions, model parameterizations, and inversion algorithms. Two representative slip models obtained from inversions of large seismic and geodetic data sets produce very different tsunami predictions at tide gauges and deep‐water pressure sensors (DART stations), despite having only secondary differences in slip distribution. This is found to be the result of the acute sensitivity of the tsunami excitation for rupture below the continental shelf in proximity to an abrupt shelf break. Iteratively perturbing seismic and geodetic inversions by constraining fault model extent along dip and strike, we obtain an optimal rupture model compatible with teleseismicPandSHwaves, regional three‐component broadband and strong‐motion seismic recordings, hr‐GNSS time series and static offsets, as well as tsunami recordings at DART stations and regional and remote tide gauges. Slip is tightly bounded between 25 and 40 km depth, the up‐dip limit of slip in the earthquake is resolved to be well‐inland of the shelf break, and the rupture extent along strike is well‐constrained. The coseismic slip increased Coulomb stress on the shallow plate boundary extending to the trench, but the frictional behavior of the megathrust below the continental slope remains uncertain.

     
    more » « less
  4. null (Ed.)
    Finite-fault models for the 2010 M w 8.8 Maule, Chile earthquake indicate bilateral rupture with large-slip patches located north and south of the epicenter. Previous studies also show that this event features significant slip in the shallow part of the megathrust, which is revealed through correction of the forward tsunami modeling scheme used in tsunami inversions. The presence of shallow slip is consistent with the coseismic seafloor deformation measured off the Maule region adjacent to the trench and confirms that tsunami observations are particularly important for constraining far-offshore slip. Here, we benchmark the method of Optimal Time Alignment (OTA) of the tsunami waveforms in the joint inversion of tsunami (DART and tide-gauges) and geodetic (GPS, InSAR, land-leveling) observations for this event. We test the application of OTA to the tsunami Green’s functions used in a previous inversion. Through a suite of synthetic tests we show that if the bias in the forward model is comprised only of delays in the tsunami signals, the OTA can correct them precisely, independently of the sensors (DART or coastal tide-gauges) and, to the first-order, of the bathymetric model used. The same suite of experiments is repeated for the real case of the 2010 Maule earthquake where, despite the results of the synthetic tests, DARTs are shown to outperform tide-gauges. This gives an indication of the relative weights to be assigned when jointly inverting the two types of data. Moreover, we show that using OTA is preferable to subjectively correcting possible time mismatch of the tsunami waveforms. The results for the source model of the Maule earthquake show that using just the first-order modeling correction introduced by OTA confirms the bilateral rupture pattern around the epicenter, and, most importantly, shifts the inferred northern patch of slip to a shallower position consistent with the slip models obtained by applying more complex physics-based corrections to the tsunami waveforms. This is confirmed by a slip model refined by inverting geodetic and tsunami data complemented with a denser distribution of GPS data nearby the source area. The models obtained with the OTA method are finally benchmarked against the observed seafloor deformation off the Maule region. We find that all of the models using the OTA well predict this offshore coseismic deformation, thus overall, this benchmarking of the OTA method can be considered successful. 
    more » « less
  5. null (Ed.)
    The 2004 Mw 9.2 earthquake and tsunami that struck North Sumatra and the Andaman-Nicobar Islands devastated coastal communities around the Indian Ocean and was the first earthquake to be analyzed by modern techniques. This earthquake and the Tohoku-Oki Mw 9.0 earthquake and tsunami in 2011 showed unexpectedly shallow megathrust slip. In the case of North Sumatra, this shallow slip was focused beneath a distinctive plateau of the accretionary prism. This intriguing seismogenic behavior and forearc structure are not well explained by existing models or by relationships observed at margins where seismogenic slip typically occurs farther landward. The input materials of the North Sumatran subduction zone are a distinctive, thick (up to 4–5 km) sequence of primarily Bengal-Nicobar Fan–related sediments. This sequence shows strong evidence for induration and dewatering and has probably reached the temperatures required for sediment-strengthening diagenetic reactions prior to accretion. The correspondence between the 2004 rupture location and the overlying prism plateau, as well as evidence for a strengthened input section, suggests the input materials are key to driving the distinctive slip behavior and long-term forearc structure. The aim of Expedition 362 is to begin to understand the nature of seismogenesis in North Sumatra through sampling these input materials and assessing their evolution, en route to understanding such processes on related convergent margins. Properties of the incoming section affect the strength of the wedge interior and base, likely promoting the observed plateau development. In turn, properties of deeper input sediment control décollement position and properties, and hence hold the key to shallow coseismic slip. During Expedition 362, two primary, riserless sites (proposed Sites SUMA-11C and SUMA-12A) will be drilled on the oceanic plate to analyze the properties of the input materials. Coring, downhole pressure and temperature measurements, and wireline logging at these sites will constrain sediment deposition rates, diagenesis, thermal and physical properties, and fluid composition. Postexpedition experimental analyses and numerical models will be employed to investigate the mechanical and frictional behavior of the input section sediments/sedimentary rocks as they thicken, accrete, and become involved in plate boundary slip system and prism development. These samples and downhole measurements will augment the internationally collected site survey bathymetric, seismic, and shallow core data that provide the regional geological framework of the margin. 
    more » « less