skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: The M w = 6.6 earthquake and tsunami of south Crete on 2020 May 2
SUMMARY On 2020 May 2, an Mw = 6.6 earthquake struck about 63 km south of Ierapetra in Crete, Greece. The earthquake generated a small tsunami which agitated local harbours. We studied this event in the context of earthquakes with seismic records in 1908, 1910, 1923, 1952, 2009 and 2013, all of similar magnitudes located south of Crete. Based on an energy-to-moment ratio, our analysis suggests that this event was neither slow nor fast, hence appropriate for using scaling laws to infer seafloor deformations. We also performed a field survey, three days after the event and present field observations from seven locations, including the island of Chrisi, where our highest measurement of 0.95 m was located. Runup along the coast of southern Crete ranged from 0.24 to 0.87 m. One tide gauge record is available for this event, and we did image analysis to obtain accurately timed water surface elevations from eyewitness videos and images. We undertook high-resolution hydrodynamic simulations using published moment tensor solutions to identify the source of the tsunami. Simulations were performed with two models, MOST (a nonlinear shallow water model) and COULWAVE (a Boussinesq-type model), to infer how different approximations of the parent equations of motion affect predictions for tsunamis of this size, which are fairly common in the Eastern Mediterranean and routinely trigger Tsunami Service Providers to issue warning messages. Based on the inter-model comparison, we conclude that the shallow-water equations are adequate in modelling this event at the distances considered, suggesting that such codes can be used to infer the tsunami source and to estimate tsunami impacts. Last, our field work revealed lack of knowledge of tsunami hazards, as most eyewitnesses remained near the waterfront, filming the associated unusual water motions instead of taking shelter on high ground.  more » « less
Award ID(s):
1830056 1940351
PAR ID:
10326815
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
230
Issue:
1
ISSN:
0956-540X
Page Range / eLocation ID:
480 to 506
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract On 18 November 2022, a large earthquake struck offshore southern Sumatra, generating a tsunami with 25 cm peak amplitude recorded at tide gauge station SBLT. OurW‐phase solution indicates a shallow dip of 6.2°, compatible with long‐period surface wave radiation patterns. Inversion of teleseismic body waves indicates a shallow slip distribution extending from about 10 km deep to near the trench with maximum slip of ∼4.1 m and seismic moment of  Nm (MW7.3). Joint modeling of seismic and tsunami data indicates a shallow rigidity of ∼23 GPa. We find a low moment‐scaled radiated energy of , similar to that of the 2010MW7.8 Mentawai event () and other tsunami earthquakes. These characteristics indicate that the 2022 event should be designated as a smaller moment magnitude tsunami earthquake compared to the other 12 well‐documented global occurrences since 1896. The 2022 event ruptured up‐dip of the 2007MW8.4 Bengkulu earthquake, demonstrating shallow seismogenic capability of a megathrust that had experienced both a deeper seismic event and adjacent shallow aseismic afterslip. We consider seismogenic behavior of shallow megathrusts and concern for future tsunami earthquakes in subduction zones globally, noting a correlation between tsunami earthquake occurrence and subducting seafloor covered with siliceous pelagic sediments. We suggest that the combination of pelagic clay and siliceous sediments and rough seafloor topography near the trench play important roles in controlling the genesis of tsunami earthquakes along Sumatra and other regions, rather than the subduction tectonic framework of accretionary or erosive margin. 
    more » « less
  2. ; ; ; ; ; ; ; ; ; (, Geophysical Journal International)
    SUMMARY Hazardous tsunamis are known to be generated predominantly at subduction zones. However, the 2018 Mw 7.5 Palu (Indonesia) earthquake on a strike-slip fault generated a tsunami that devastated the city of Palu. The mechanism by which this tsunami originated from such an earthquake is being debated. Here we present near-field ground motion (GPS) data confirming that the earthquake attained supershear speed, i.e. a rupture speed greater than the shear wave speed of the host medium. We subsequently study the effect of this supershear rupture on tsunami generation by coupling the ground motion to a 1-D non-linear shallow-water wave model accounting for both time-dependent bathymetric displacement and velocity. With the local bathymetric profile of Palu bay around a tidal station, our simulations reproduce the tsunami arrival and motions observed by CCTV cameras. We conclude that Mach (shock) fronts, generated by the supershear speed, interacted with the bathymetry and contributed to the tsunami. 
    more » « less
  3. ; ; ; ; ; ; (, Frontiers in Earth Science)
    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
  4. ; ; ; (, ESS Open Archive)
    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
  5. ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract Most large tsunamis are generated by earthquakes on offshore plate boundary megathrusts. The primary factors influencing tsunami excitation are the seismic moment, faulting geometry, and depth of the faulting. Efforts to provide rapid tsunami warning have emphasized seismic and geodetic methods for quickly determining the event size and faulting geometry. It remains difficult to evaluate the updip extent of rupture, which has significant impact on tsunami excitation. TeleseismicPwaves can constrain this issue; slip under deep water generates strongpwPwater reverberations that persist as ringingPcodaafter the directPphases from the faulting have arrived. Event‐averagedPcoda/Pamplitude measures at large epicentral distances (>80°), tuned to the dominant periods of deep waterpwP(~12–15 s), correlate well with independent models of whether slip extends to near the trench or not. Data at closer ranges (30° to 80°) reduce the time lag needed for inferring the updip extent of rupture to <15 min. Arrival ofPPandPPPphases contaminates closer distancePcodameasures, but this can be suppressed by azimuthal or distance binning of the measures. Narrowband spectral ratio measures and differential magnitude measures ofPcodaand directP(mB) perform comparably to broader band root‐mean‐square (RMS) measures.Pcoda/Plevels for large nonmegathrust events are also documented. Rapid measurement ofPcoda/Pmetrics after a large earthquake can supplement quick moment tensor determinations to enhance tsunami warnings; observation of largePcodalevels indicates that shallow submarine rupture occurred and larger than typical tsunami (for givenMW) can be expected. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email