We define the geometric and kinematic characteristics of the fault activated during the Mw = 6.3, 6 April 2009 L'Aquila earthquake, from the modeling of Envisat and COSMO‐SkyMed (the first ever X‐band interferogram inverted for a coseismic dislocation study) DInSAR interferograms. Our best‐fit solution for the main shock is represented by a normal fault ∼16 km long and ∼12 km wide, with a small right‐lateral component, dipping 47°SW with a maximum slip of ∼90 cm. Although the seismic dislocation probably ended at 1 km depth, the updip projection of the fault plane corresponds to the northern segment of the mapped Paganica–S. Demetrio fault, where alignment of surface breaks was observed in the field. The absence of this fault in existing seismic source catalogues suggests that an improved approach, involving detailed surface and subsurface geological and geophysical investigations, is needed for a better assessment of the seismic hazard at the local scale.

We investigate a large geodetic data set of interferometric synthetic aperture radar (InSAR) and GPS measurements to determine the source parameters for the three main shocks of the 2016 Central Italy earthquake sequence on 24 August and 26 and 30 October (Mw 6.1, 5.9, and 6.5, respectively). Our preferred model is consistent with the activation of four main coseismic asperities belonging to the SW dipping normal fault system associated with the Mount Gorzano‐Mount Vettore‐Mount Bove alignment. Additional slip, equivalent to a Mw ~ 6.1–6.2 earthquake, on a secondary (1) NE dipping antithetic fault and/or (2) on a WNW dipping low‐angle fault in the hanging wall of the main system is required to better reproduce the complex deformation pattern associated with the greatest seismic event (the Mw 6.5 earthquake). The recognition of ancillary faults involved in the sequence suggests a complex interaction in the activated crustal volume between the main normal faults and the secondary structures and a partitioning of strain release.

Several independent indicators imply a high probability of a great (M N 8) earthquake rupture of the subduction megathrust under the Mentawai Islands of West Sumatra. The human consequences of such an event depend crucially on its tsunamigenic potential, which in turn depends on unpredictable details of slip distribution on the megathrust and how resulting seafloor movements and the propagating tsunami waves interact with bathymetry. Here we address the forward problem by modelling about 1000 possible complex earthquake ruptures and calculating the seafloor displacements and tsunami wave height distributions that would result from the most likely 100 or so, as judged by reference to paleogeodetic data. Additionally we carry out a systematic study of the importance of the location of maximum slip with respect to the morphology of the fore-arc complex. Our results indicate a generally smaller regional tsunami hazard than was realised in Aceh during the December 2004 event, though more than 20% of simulations result in tsunami wave heights of more than 5 m for the southern Sumatran cities of Padang and Bengkulu. The extreme events in these simulations produce results which are consistent with recent deterministic studies. The study confirms the sensitivity of predicted wave heights to the distribution of slip even for events with similar moment and reproduces Plafker's rule of thumb. Additionally we show that the maximum wave height observed at a single location scales with the magnitude though data for all magnitudes exhibit extreme variability. Finally, we show that for any coastal location in the near field of the earthquake, despite the complexity of the earthquake rupture simulations and the large range of magnitudes modelled, the timing of inundation is constant to first order and the maximum height of the modelled waves is directly proportional to the vertical coseismic displacement experienced at that point. These results may assist in developing tsunami preparedness strategies around the Indian Ocean and in particular along the coasts of western Sumatra.

[1] We model the spatial and temporal evolution of seismicity during the 1997 Umbria-Marche seismic sequence in terms of subsequent failures promoted by fluid flow. The diffusion process of pore-pressure relaxation is represented as a pressure perturbation generated by coseismic stress changes and propagating through a fluid saturated medium

No abstract

No abstract

S U M M A R YWith the continued growth in availability of differential interferometry from synthetic aperture radar and GPS data, space based geodesy has been widely applied to image the coseismic displacement field and to retrieve the static dislocation over the fault plane for almost all the significant earthquakes of the past two decades. This is performed by linear data inversion over a set of subfaults, generally characterized by a constant and predefined or manually adjusted dimensions.In this paper, we propose a new algorithm to automatically retrieve an optimized fault subdivision in the linear inversion of coseismic geodetic data. The code iteratively keeps the parameter resolution close to a predefined high value. We first discuss the rationale supporting our algorithm and, after a detailed description of its implementation, we analyse the advantages of its introduction in the data inversion.The algorithm was tested against an exhaustive range of synthetic and real data sets and fault mechanisms. Among them, we present the results for the M w 6.2, 2009 L'Aquila (Central Italy) earthquake and compare the new and previously published slip distributions showing the disappearance of misleading slip pattern and the increased resolution for shallower zones.

[1] We analyze the coseismic stress perturbation during the 17 June 2000 south Iceland seismic sequence; the main shock (M s 6.6) was followed by three large events within a few tens of seconds (8, 26, and 30 s) located within 80 km. The aim of this paper is to investigate short-term fault interaction and instantaneous triggering. This happens when a fault perturbed by a stress change fails before the end of the transient stress perturbation. We compute the shear, normal, and Coulomb stress changes as functions of time in a stratified elastic half-space by using discrete wave number and reflectivity methods. We calculate dynamic stresses caused by the main shock at the hypocenters of these three subsequent events. Our numerical results show that the onset of the last two events is slightly delayed with respect to the arrival time of the second positive peak of Coulomb stress variation, while the first event occurred after the first positive stress peak. We have also analyzed the response of a spring slider system representing a fault governed by a rate-and state-dependent friction law, perturbed by shear and normal stress variations caused by the main shock. The fault response to the computed stress perturbations is always clock advanced. We have found suitable constitutive parameters of the modeled fault that allow the instantaneous dynamic triggering of these three earthquakes. If the initial sliding velocity is comparable with the tectonic loading velocity, we obtained failure times close to the observed origin times for low values of the initial effective normal stress.