On April 6, 2009, 01:32:39 GMT, the city of L'Aquila was struck by a Mw 6.3 earthquake that killed 307 people, causing severe destruction and ground cracks in a wide area around the epicenter. Four days before the main shock we augmented the existing permanent GPS network with five GPS stations of the Central Apennine Geodetic Network (CaGeoNet) bordering the L'Aquila basin. The maximum horizontal and vertical coseismic surface displacements detected at these stations was 10.39 ± 0.45 cm and −15.64 ± 1.55 cm, respectively. Fixing the strike direction according to focal mechanism estimates, we estimated the source geometry with a non linear inversion of the geodetic data. Our best fitting fault model is a 13 × 15.7 km2 rectangular fault, SW‐dipping at 55.3 ± 1.8°, consistent with the position of observed surface ruptures. The estimated slip (495 ± 29 mm) corresponds to a 6.3 moment magnitude, in excellent agreement with seismological data.

[1] We (re)analyzed the source of the 26 December 2004 Sumatra-Andaman earthquake and tsunami through a nonlinear joint inversion of an inhomogeneous data set made up of tide gauges, satellite altimetry, and far-field GPS recordings. The purpose is twofold: (1) the retrieval of the main kinematics rupture parameters (slip, rake, and rupture velocity) and (2) the inference of the rigidity of the source zone. We independently estimate the slip from tsunami data and the seismic moment from geodetic data to derive the rigidity. Our results confirm that the source of the 2004 Sumatra-Andaman earthquake has a complex geometry, constituted by three main slip patches, with slip peaking at $30 m in the southern part of the source. The rake direction rotates counterclockwise at the northern part of the source, according to the direction of convergence along the trench. The rupture velocity is higher in the deeper than in the shallower part of the source, consistent with the expected increase of rigidity with depth. It is also lower in the northern part, consistent with known variations of the incoming plate properties and shear velocity. Our model features a rigidity (20-30 GPa) that is lower than the preliminary reference Earth model (PREM) average for the seismogenic volume. The source rigidity is one of the factors controlling the tsunami genesis: for a given seismic moment, the lower the rigidity, the higher the induced seafloor displacement. The general consistence between our source model and previous studies supports the effectiveness of our approach to the joint inversion of geodetic and tsunami data for the rigidity estimation.

S U M M A R YThe post-seismic response of a viscoelastic Earth to a seismic dislocation can be computed analytically within the framework of normal-modes, based on the application of propagator methods. This technique, widely documented in the literature, suffers from several shortcomings; the main drawback is related to the numerical solution of the secular equation, whose degree increases linearly with the number of viscoelastic layers so that only coarse-layered models are practically solvable. Recently, a viable alternative to the standard normal-mode approach, based on the Post-Widder Laplace inversion formula, has been proposed in the realm of postglacial rebound models. The main advantage of this method is to bypass the explicit solution of the secular equation, while retaining the analytical structure of the propagator formalism. At the same time, the numerical computation is much simplified so that additional features such as linear non-Maxwell rheologies can be simply implemented. In this work, for the first time, we apply the Post-Widder Laplace inversion formula to a post-seismic rebound model. We test the method against the standard normal-mode solution and we perform various benchmarks aimed to tune the algorithm and to optimize computation performance while ensuring the stability of the solution. As an application, we address the issue of finding the minimum number of layers with distinct elastic properties needed to accurately describe the post-seismic relaxation of a realistic Earth model. Finally, we demonstrate the potentialities of our code by modelling the post-seismic relaxation after the 2004 Sumatra-Andaman earthquake comparing results based upon Maxwell and Burgers rheologies.

We adopted a multidisciplinary approach to investigate the seismotectonic scenario of the 30 October 2016, Mw 6.5, Norcia earthquake, the largest shock of the 2016–2017 central Italy earthquake sequence. First, we used seismological and geodetic data to infer the dip of the main slip patch of the seismogenic fault that turned out to be rather low‐angle (~37°). To evaluate whether this is an acceptable dip for the main seismogenic source, we modeled earthquake deformation using single‐ and multiple‐fault models deduced from aftershock pattern analyses. These models show that the coseismic deformation generated by the Norcia earthquake is coherent with slip along a rather shallow‐dipping plane. To understand the geological significance of this solution, we reconstructed the subsurface architecture of the epicentral area. As the available data are not robust enough to converge on a single fault model, we built three different models encompassing all major geological evidence and the associated uncertainties, including the tectonic style and the location of major décollement levels. In all models the structures derived from the contractional phase play a significant role: from controlling segmentation to partially reusing inherited faults, to fully reactivating in extension a regional thrust, geometrically compatible with the source of the Norcia earthquake. Based on our conclusions, some additional seismogenic sources falling in the eastern, external portions of the Apennines may coincide with inherited structures. This may be a common occurrence in this region of the chain, where the inception of extension is as recent as Middle‐Upper Pleistocene.

The dynamics governing the movement of the radon are complex and dependent on many factors. In the present study, we characterise the nature of temporal variations of 2-hourly and daily radon measurements in several monitoring sites of the Italian Radon mOnitoring Network (IRON) in Italy. By means of continuous wavelet transformation, a spectral analysis in time-frequency domain is performed. The results reveal that there are sub-daily, daily and yearly persistent periodicities that are common for all the stations. We observe structural seasonal breaks, that occur at the same frequency but at distinct time. Variations in radon concentration and local temperature are studied in terms of frequency contents and synchronicity. When analysing several long time series together, it is evident that the phase difference at low frequency movements (365-day period) between the radon and local temperature time series is depending on the sites’ location and therefore strongly controlled by local factors. This could at least partially explain the apparently contrasting results available in the literature obtained investigating smaller dataset about the relationships between temperature and radon variations. On the other hand, results show that all radon time series are characterised by marked cycles at 1 and 365-days and less evident cycles at 0.5-day and 180-days. They would be all ascribable to environmental-climatic factors: the short-period cycles to temperature and pressure variations, the long-period cycles also to seasonal rainfall variations.