skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Friday, April 12 until 2:00 AM ET on Saturday, April 13 due to maintenance. We apologize for the inconvenience.

Title: Benchmarking the Optimal Time Alignment of Tsunami Waveforms in Nonlinear Joint Inversions for the Mw 8.8 2010 Maule (Chile) Earthquake
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
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Frontiers in Earth Science
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract

    The 2018 Palu tsunami contributed significantly to the devastation caused by the associated7.5 earthquake. This began a debate about how the moderate size earthquake triggered such a large tsunami within Palu Bay, with runups of more than 10 m. The possibility of a large component of vertical coseismic deformation and submarine landslides have been considered as potential explanations. However, scarce instrumental data have made it difficult to resolve the potential contributions from either type of source. We use tsunami waveforms derived from social media videos in Palu Bay to model the possible sources of the tsunami. We invert InSAR data with different fault geometries and use the resulting seafloor displacements to simulate tsunamis. The coseismic sources alone cannot match both the video‐derived time histories and surveyed runups. Then we conduct a tsunami source inversion using the video‐derived time histories and a tide gauge record as inputs. We specify hypothetical landslide locations and solve for initial tsunami elevation. Our results, validated with surveyed runups, show that a limited number of landslides in southern Palu Bay are sufficient to explain the tsunami data. The Palu tsunami highlights the difficulty in accurately capturing with tide gauges the amplitude and timing of short period waves that can have large impacts at the coast. The proximity of landslides to locations of high fault slip also suggests that tsunami hazard assessment in strike‐slip environments should include triggered landslides, especially for locations where the coastline morphology is strongly linked to fault geometry.

    more » « less
  3. Abstract

    We use tsunami waveforms recorded on deep water absolute pressure gauges (Deep‐ocean Assessment and Reporting of Tsunamis), coastal tide gauges, and a temporary array of seafloor differential pressure gauges (DPG) to study the tsunami generated by the 15 July 2009 magnitude 7.8 Dusky Sound, New Zealand, earthquake. We first use tsunami waveform inversion applied to Deep‐ocean Assessment and Reporting of Tsunamis seafloor pressure gauge and coastal tide gauge data to estimate the fault slip distribution of the Dusky Sound earthquake. This fault slip estimate is then used to generate synthetic tsunami waveforms at each of the DPG sites. DPG instruments are unfortunately not well calibrated, but comparison of the synthetic tsunami waveforms to those observed at each DPG site allows us to determine an appropriate amplitude scaling to apply. We next use progressive data assimilation of the amplitude‐scaled DPG observations to retrospectively forecast the Dusky Sound tsunami wavefields and find a good match between forecast and observed tsunami wavefields at the Charleston tide gauge station on the west coast of New Zealand's South Island. While an advantage of the data assimilation method is that no initial condition is needed, we find that our forecast is improved by merging tsunami forward modeling from a rapid W‐phase earthquake source solution with the data assimilation method.

    more » « less
  4. Abstract On 19 September 2022, a major earthquake struck the northwestern Michoacán segment along the Mexican subduction zone. A slip model is obtained that satisfactorily explains geodetic, teleseismic, and tsunami observations of the 2022 event. The preferred model has a compact large-slip patch that extends up-dip and northwestward from the hypocenter and directly overlaps a 1973 Mw 7.6 rupture. Slip is concentrated offshore and below the coast at depths from 10 to 30 km with a peak value of ∼2.9 m, and there is no detected coseismic slip near the trench. The total seismic moment is 3.1×1020  N·m (Mw 7.6), 72% of which is concentrated in the first 30 s. Most aftershocks are distributed in an up-dip area of the mainshock that has small coseismic slip, suggesting near-complete strain release in the large-slip patch. Teleseismic P waveforms of the 2022 and 1973 earthquakes are similar in duration and complexity with high cross-correlation coefficients of 0.68–0.98 for long P to PP signal time windows, indicating that the 2022 earthquake is a quasi-repeat of the 1973 earthquake, possibly indicating persistent frictional properties. Both the events produced more complex P waveforms than comparable size events along Guerrero and Oaxaca, reflecting differences in patchy locking of the Mexican megathrust. 
    more » « less
  5. Abstract

    Nearly every type of slip behavior of which megathrusts are known to be capable of has been observed within a narrow ~75‐km‐wide strip of the Sunda megathrust—the Banyak Islands subsegment. The diverse list of recorded slip events starts with the great 1861MW ~ 8.5 earthquake and includes the 1907MW ~ 8.2 shallow tsunami earthquake, the 1966–1981 slow slip event, the great 2005MW8.6 Nias‐Simeulue earthquake and associated afterslip, many small (M < 6) to moderate (6 M < 7) earthquakes, and the large 6 April 2010MW7.8 Banyak Islands earthquake. In this paper we map out the spatial and temporal relationships between these slip events to search for a pattern and possible controls on slip behavior. We use GPS and coral geodetic data to derive a coseismic slip distribution for the relatively understudied 2010 event and find that this event fits like a puzzle piece into a low slip patch left by the 2005 event. We compare the recent slip sequence with the historical sequence associated with the 1861 event and find that the modern slip pattern could be a repeat of the historical slip pattern, possibly suggesting the importance of stationary fault properties in controlling slip. The diverse set of slip events mostly fit our expectations within depth‐dependent rupture zones; however, the overlaps between events of different slip types may complicate this conceptual model. Along strike, we suggest that a subducting fracture zone could act to diversify the slip behavior, and we explore possible mechanisms for this slip control.

    more » « less