An official website of the United States government
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
INTER-COUPLED TSUNAMI MODELLING THROUGH AN ABSORBING-GENERATING BOUNDARY
For many practical and theoretical purposes, various types of tsunami wave models have been developed and utilized so far. Some distinction among them can be drawn based on governing equations used by the model. Shallow water equations and Boussinesq equations are probably most typical ones among others since those are computationally efficient and relatively accurate compared to 3D Navier-Stokes models. From this idea, some coupling effort between Boussinesq model and shallow water equation model have been made (e.g., Son et al. (2011)). In the present study, we couple two different types of tsunami models, i.e., nondispersive shallow water model of characteristic form(MOST ver.4) and dispersive Boussinesq model of non-characteristic form(Son and Lynett (2014)) in an attempt to improve modelling accuracy and efficiency.Recorded Presentation from the vICCE (YouTube Link): https://youtu.be/cTXybDEnfsQ
more »
« less
- PAR ID:
- 10221513
- Date Published:
- Journal Name:
- Coastal Engineering Proceedings
- Issue:
- 36v
- ISSN:
- 0589-087X
- Page Range / eLocation ID:
- 38
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY Tsunami generation by offshore earthquakes is a problem of scientific interest and practical relevance, and one that requires numerical modelling for data interpretation and hazard assessment. Most numerical models utilize two-step methods with one-way coupling between separate earthquake and tsunami models, based on approximations that might limit the applicability and accuracy of the resulting solution. In particular, standard methods focus exclusively on tsunami wave modelling, neglecting larger amplitude ocean acoustic and seismic waves that are superimposed on tsunami waves in the source region. In this study, we compare four earthquake-tsunami modelling methods. We identify dimensionless parameters to quantitatively approximate dominant wave modes in the earthquake-tsunami source region, highlighting how the method assumptions affect the results and discuss which methods are appropriate for various applications such as interpretation of data from offshore instruments in the source region. Most methods couple a 3-D solid earth model, which provides the seismic wavefield or at least the static elastic displacements, with a 2-D depth-averaged shallow water tsunami model. Assuming the ocean is incompressible and tsunami propagation is negligible over the earthquake duration leads to the instantaneous source method, which equates the static earthquake seafloor uplift with the initial tsunami sea surface height. For longer duration earthquakes, it is appropriate to follow the time-dependent source method, which uses time-dependent earthquake seafloor velocity as a forcing term in the tsunami mass balance. Neither method captures ocean acoustic or seismic waves, motivating more advanced methods that capture the full wavefield. The superposition method of Saito et al. solves the 3-D elastic and acoustic equations to model the seismic wavefield and response of a compressible ocean without gravity. Then, changes in sea surface height from the zero-gravity solution are used as a forcing term in a separate tsunami simulation, typically run with a shallow water solver. A superposition of the earthquake and tsunami solutions provides an approximation to the complete wavefield. This method is algorithmically a two-step method. The complete wavefield is captured in the fully coupled method, which utilizes a coupled solid Earth and compressible ocean model with gravity. The fully coupled method, recently incorporated into the 3-D open-source code SeisSol, simultaneously solves earthquake rupture, seismic waves and ocean response (including gravity). We show that the superposition method emerges as an approximation to the fully coupled method subject to often well-justified assumptions. Furthermore, using the fully coupled method, we examine how the source spectrum and ocean depth influence the expression of oceanic Rayleigh waves. Understanding the range of validity of each method, as well as its computational expense, facilitates the selection of modelling methods for the accurate assessment of earthquake and tsunami hazards and the interpretation of data from offshore instruments.more » « less
-
We present an algorithm to solve the dispersive depth-averaged Serre--Green--Naghdi equations using patch-based adaptive mesh refinement. These equations require adding additional higher derivative terms to the nonlinear shallow water equations. This has been implemented as a new component of the open source GeoClaw software that is widely used for modeling tsunamis, storm surge, and related hazards, improving its accuracy on shorter wavelength phenomena. We use a formulation that requires solving an elliptic system of equations at each time step, making the method implicit. The adaptive algorithm allows different time steps on different refinement levels and solves the implicit equations level by level. Computational examples are presented to illustrate the stability and accuracy on a radially symmetric test case and two realistic tsunami modeling problems, including a hypothetical asteroid impact creating a short wavelength tsunami for which dispersive terms are necessary.more » « less
-
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
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.9753/icce.v36v.currents.38
