On 16 September 2015, the M_W = 8.2 Illapel megathrust earthquake ruptured the Central Chilean margin. Combining inversions of displacement measurements and seismic waveforms with high frequency (HF) teleseismic backprojection, we derive a comprehensive description of the rupture, which also predicts deep ocean tsunami wave heights. We further determine moment tensors and obtain accurate depth estimates for the aftershock sequence. The earthquake nucleated near the coast but then propagated to the north and updip, attaining a peak slip of 5–6 m. In contrast, HF seismic radiation is mostly emitted downdip of the region of intense slip and arrests earlier than the long period rupture, indicating smooth slip along the shallow plate interface in the final phase. A superficially similar earthquake in 1943 with a similar aftershock zone had a much shorter source time function, which matches the duration of HF seismic radiation in the recent event, indicating that the 1943 event lacked the shallow slip.

The Maule earthquake (2010 February 27, M_w 8.8, Chile) broke the subduction megathrust along a previously locked segment. Based on an international aftershock deployment, catalogues of precisely located aftershocks have become available. Using 23 well-located aftershocks, we calibrate the classic teleseismic backprojection procedure to map the high-frequency seismic radiation emitted during the earthquake. The calibration corrects traveltimes in a standard earth model both with a static term specific to each station, and a ‘dynamic’ term specific to each combination of grid point and station. The second term has been interpolated over the whole slipping area by kriging, and is about an order of magnitude smaller than the static term. This procedure ensures that the teleseismic images of rupture development are properly located with respect to aftershocks recorded with local networks and does not depend on accurate hypocentre location of the main shock. We track a bilateral rupture propagation lasting ∼160 s, with its dominant branch rupturing northeastwards at about 3 km s⁻¹. The area of maximum energy emission is offset from the maximum coseismic slip but matches the zone where most plate interface aftershocks occur. Along dip, energy is preferentially released from two disconnected interface belts, and a distinct jump from the shallower belt to the deeper one is visible after about 20 s from the onset. However, both belts keep on being active until the end of the rupture. These belts approximately match the position of the interface aftershocks, which are split into two clusters of events at different depths, thus suggesting the existence of a repeated transition from stick-slip to creeping frictional regime.

A new approach for estimating methane emission at mud volcanoes is here proposed based on measurements of the seismic tremor on their surface. Data obtained at the Dashgil mud volcano in Azerbaijan reveal the presence of energy bursts characterized by well-determined features (i.e. waveforms, spectra and polarization properties) that can be associated with bubbling at depth. Counting such events provides a possible tool for monitoring gas production in the reservoir, thus minimizing logistic troubles and representing a cheap and effective alternative to more complex approaches. Specifically, we model the energy bursts as the effect of resonant gas bubbles at depth. This modelling allows to estimate the dimension of the bubbles and, consequently, the gas outflow from the main conduit in the assumption that all emissions from depth occur by bubble uprising. The application of this model to seismic events detected at the Dashgil mud volcano during three sessions of measurements carried out in 2006 and 2007 provides gas flux estimates that are in line with those provided by independent measurements at the same structure. This encouraging result suggests that the one here proposed could be considered a new promising, cheap and easy to apply tool for gas flux measurements in bubbling gas seepage areas.