Search

SUMMARYWe investigated the kinematic rupture model of the 2018 Mw 6.8 Zakynthos, Ionian Sea (Greece), earthquake by using a non-linear joint inversion of strong motion data, high-rate GPS time-series and static coseismic GPS displacements. We also tested inversion results against tide-gauge recordings of the small tsunami generated in the Ionian Sea. In order to constrain the fault geometry, we performed several preliminary kinematic inversions by assuming the parameter values resulting from different published moment tensor solutions. The lowest cost function values were obtained by using the geometry derived from the United States Geological Survey (USGS) focal solution. Between the two conjugate USGS planes, the rupture model which better fits the data is the one with the N9°E-striking 39°ESE-dipping plane. The rupture history of this model is characterized by a bilateral propagation, featuring two asperities; a main slip patch extending between 14 and 28 km in depth, 9 km northeast from the nucleation and a slightly shallower small patch located 27 km southwest from the nucleation. The maximum energy release occurs between 8 and 12 s, when both patches are breaking simultaneously. The maximum slip is 1.8 m and the total seismic moment is 2.4 × 1019 Nm, corresponding to a Mw value of 6.8. The slip angle shows a dominant right-lateral strike-slip mechanism, with a minor reverse component that increases on the deeper region of the fault. This result, in addition to the observed possibility of similar mechanisms for previous earthquakes occurred in 1959 and 1997, suggests that the tectonic deformation between the Cephalonia Transform Fault Zone and the northern tip of the Hellenic Arc Subduction zone may be accommodated by prevailing right lateral low-dipping faults, occurring on re-activated structures previously experiencing (until Pliocene) compressional regime. Comparison of predicted and observed tsunami data suggests the need of a better characterization of local harbour response for this type of relatively short-wavelength events, which is important in the context of tsunami early warning. However, the suggested dominantly strike-slip character would in turn imply a reduced tsunami hazard as compared to a dominant thrust faulting regime from this source region.

SUMMARY In this study, we investigate the rupture history of the 2009 April 6 (Mw 6.1) L’Aquila normal faulting earthquake by using a non-linear inversion of strong motion, GPS and DInSAR data. Both the separate and joint inversion solutions reveal a complex rupture process and a heterogeneous slip distribution. Slip is concentrated in two main asperities: a smaller shallow patch of slip located updip from the hypocentre and a second deeper and larger asperity located southeastwards along-strike direction. The key feature of the source process emerging from our inverted models concerns the rupture history, which is characterized by two distinct stages. The first stage begins with rupture nucleation and with updip propagation at relatively high (∼4.0 km s−1), but still subshear, rupture velocity. The second stage starts nearly 2.0–2.5 s after nucleation and it is characterized by the along-strike rupture propagation. The largest and deeper asperity fails during this stage of the rupture process. The rupture velocity is larger in the updip than in the along-strike direction. The updip and along-strike rupture propagation are separated in time and associated with a Mode II and a Mode III crack, respectively. The comparison between the source models inferred in this study with the Poisson ratio anomalies in the crustal volume containing the fault plane allows the interpretation of the delay in along-strike rupture propagation in terms of a structural control of the rupture history. Our results show that the L’Aquila earthquake featured a very complex rupture, with strong spatial and temporal heterogeneities suggesting a strong frictional and/or structural control of the rupture process.

SUMMARY In this study, we aim to understand the dependence of the critical slip weakening distance (D c) on the final slip (D tot) during the propagation of a dynamic rupture and the consistency of their inferred correlation. To achieve this goal we have performed a series of numerical tests suitably designed to validate the adopted numerical procedure and to verify the actual capability in measuring D c. We have retrieved two kinematic rupture histories from spontaneous dynamic rupture models governed by a slip weakening law in which a constant D c distribution on the fault plane as well as a constant D c/D tot ratio are assumed, respectively. The slip velocity and the shear traction time histories represent the synthetic ‘real’ target data which we aim to reproduce. We use a 3-D traction-at-split nodes numerical procedure to image the dynamic traction evolution by assuming our modelled slip velocity as a boundary condition on the fault plane. We assume a regularized Yoffe function as source time function in our modelling attempts and we measure the critical slip weakening distance from the inferred traction versus slip curves at each point on the fault. We compare the inferred values with those of the target dynamic models. Our numerical tests show that fitting the slip velocity functions of the target models at each point on the fault plane is not enough to retrieve good traction evolution curves and to obtain reliable measures of D c. We find that the estimation of D c is very sensitive to any small variation of the slip velocity function. An artificial correlation between D c/D tot is obtained when a fixed shape of slip velocity is assumed on the fault (i.e. constant rise time and constant time for positive acceleration) which differs from that of the target model. We point out that the estimation of fracture energy (breakdown work) on the fault is not affected by biases in measuring D c.

SUMMARY We introduce a rapid and robust, energy-duration procedure, based on the Haskell, extended-source model, to obtain an earthquake moment and a moment magnitude, MED. Using seismograms at teleseismic distances (30°–90°), this procedure combines radiated seismic energy measures on the P to S interval of broadband signals and source duration measures on high-frequency, P-wave signals. The MED energy-duration magnitude is scaled to correspond to the Global Centroid-Moment Tensor (CMT) moment-magnitude, MCMTw, and can be calculated within about 20 min or less after origin time (OT). The measured energy and duration values also provide the energy-to-moment ratio, Θ, used for identification of tsunami earthquakes. The MED magnitudes for a set of recent, large earthquakes match closely MCMTw, even for the largest, great earthquakes; these results imply that the MED measure is accurate and does not saturate. After the 2004 December 26 Sumatra-Andaman mega-thrust earthquake, magnitude estimates available within 1 hr of OT ranged from M= 8.0 to 8.5, the CMT magnitude, available about 3 hr after OT, was MCMTw= 9.0, and, several months after the event, Mw= 9.1–9.3 was obtained from analysis of the earth normal modes. The energy-duration magnitude for this event is MED= 9.2, a measure that is potentially available within 20 min after OT. After the 2006 July 17, Java earthquake, the magnitude was evaluated at M= 7.2 at 17 min after OT, the CMT magnitude, available about 1 hr after OT, was MCMTw= 7.7; the energy-duration results for this event give MED= 7.8, with a very long source duration of about 160 s, and a very low Θ value, indicating a possible tsunami earthquake.

Summary This paper presents a study of the tsunami following the disastrous earthquake that occurred on the 8th of September 1905 in Calabria, southern Italy. The shock caused devastation in many towns and villages leading to more than 500 victims. According to coeval sources the tsunami was not catastrophic, but was large enough to inundate low lying lands in several coastal segments and to affect boats. No specific study was ever devoted to this tsunami in the literature. Our analysis is carried out by means of numerical modelling and aims at filling this gap. Since the source fault of the earthquake has not yet been identified, we study three tsunamis produced by likely potential sources, according to macroseismic data and the seismotectonic knowledge of the region, namely the Capo Vaticano (CV) fault, the Vibo Valentia (VV) fault and the Lamezia (LA) fault. We compute tsunami waves under the hypotheses that the sea floor deformation is caused by a shear double-couple dislocation taking place over a rectangular fault and that sea waves propagate according to the shallow water approximation. The three tsunamis cause either a pure initial depression (VV and LA) or a predominant initial depression (CV) of the sea surface in the gulf of St. Eufemia. One very relevant finding is that tsunami energy is trapped within a narrow channel along the Calabrian coast to the north of the source with the consequence that wave amplitudes computed here are much larger than the values expected in the case of an ordinary energy decay law with distance; energy trapping along southern coastal segments is less efficient. Only a small fraction of tsunami energy penetrates the Straits of Messina to the south. Propagation to the west is affected by the Aeolian islands that produce a quasi-shadow zone on the side opposite to the tsunami impact. The LA tsunami is found to be too weak to be compatible with observations. A sensitivity analysis concerning some parameters of the faults (namely length, width and slip magnitude) has been performed in order to explore their influence on the tsunami features. The main result of this paper is that we cannot discriminate between the CV and VV sources: both are in agreement with some historical observations, whilst they do not match others. We believe that a better tuning of the source mechanism over these faults can improve the agreement between model computations and historical observations.

Summary Deep and large submarine slumps may generate tsunamis as disastrous as tsunamis of tectonic origin. Such a landslide is likely to be the origin of the 1998 July 17 tsunami of Papua New Guinea, the deadliest tsunami in the last 50 years. Water waves devastated a 20 km stretch of coastline, wiping out three villages and killing more than 2200 people. A numerical model has been developed to study the efficiency of deep slumps in producing tsunamis and has been applied to the Papua New Guinea event. The landslide is treated as the flow of a homogeneous gravity-driven continuum governed by a rheological law. Water waves are generated by sea-bottom displacements induced by the landslide. The shallow-water approximation is adopted for both the landslide and the associated water waves. The resulting differential equations are solved by a finite difference method based on shock-capturing. The shallow-water hypothesis is tested by comparison with a model solving Navier–Stokes equations for a mixture of water and sediments. Sensitivity tests carried out for a 2-D simplified geometry show that the water surface profile depends strongly on the constitutive law of the landslide. The 1998 event is simulated numerically by the shallow-water model, testing different friction laws. The observed inundation height distribution is well reproduced by the model for a volume of 4 km3, with its top located at a water depth of 550 m, and sliding with a Coulomb-type friction law over a distance of 5 km.