No abstract

On 16 September 2015, the MW = 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.

A GNSS water vapour tomography system developed to reconstruct spatially resolved humidity fields in the troposphere is described. The tomography system was designed to process the slant path delays of about 270 German GNSS stations in near real-time with a temporal resolution of 30 minutes, a horizontal resolution of 40 km and a vertical resolution of 500 m or better. After a short introduction to the GPS slant delay processing the framework of the GNSS tomography is described in detail. Different implementations of the iterative algebraic reconstruction techniques (ART) used to invert the linear inverse problem are discussed. It was found that the multiplicative techniques (MART) provide the best results with least processing time, i. e. a tomographic reconstruction of about 26000 slant delays on a 8280 cell grid can be obtained in less than 10 minutes. Different iterative reconstruction techniques are compared with respect to their convergence behaviour and some numerical parameters. The inversion can be considerably stabilized by using additional non-GNSS observations and implementing various constraints. Different strategies for initialising the tomography and utilizing extra information are discussed. At last an example of a reconstructed field of the wet refractivity is presented and compared to the corresponding distribution of the integrated water vapour, an analysis of a numerical weather model (COSMO-DE) and some radiosonde profiles.

Along a subduction zone, great megathrust earthquakes recur either after long seismic gaps lasting several decades to centuries or over much shorter periods lasting hours to a few years when cascading successions of earthquakes rupture nearby segments of the fault. We analyze a decade of continuous Global Positioning System observations along the South American continent to estimate changes in deformation rates between the 2010 Maule (M8.8) and 2015 Illapel (M8.3) Chilean earthquakes. We find that surface velocities increased after the 2010 earthquake, in response to continental‐scale viscoelastic mantle relaxation and to regional‐scale increased degree of interplate locking. We propose that increased locking occurs transiently during a super‐interseismic phase in segments adjacent to a megathrust rupture, responding to bending of both plates caused by coseismic slip and subsequent afterslip. Enhanced strain rates during a super‐interseismic phase may therefore bring a megathrust segment closer to failure and possibly triggered the 2015 event.

The potential temporal shifts in the integrated water vapor (IWV) time series obtained from reprocessed data acquired from global navigation satellite systems (GNSS) were comprehensively investigated. A statistical test, the penalized maximal t test modified to account for first-order autoregressive noise in time series (PMTred), was used to identify the possible mean shifts (changepoints) in the time series of the difference between the GPS IWV and the IWV obtained from the European Centre for Medium-Range Weather Forecasts (ECMWF) interim reanalysis (ERA-Interim) data. The approach allows for identification of the changepoints not only in the GPS IWV time series but also in ERA-Interim. The IWV difference time series formed for 101 GPS sites were tested, where 47 of them were found to contain in total 62 changepoints. The results indicate that 45 detected changepoints were due to the inconsistencies in the GPS IWV time series, and 16 were related to ERA-Interim, while one point was left unverified. After the correction of the mean shifts for the GPS data, an improved consistency in the IWV trends is evident between nearby sites, while a better agreement is seen between the trends from the GPS and ERA-Interim data on a global scale. In addition, the IWV trends estimated for 47 GPS sites were compared to the corresponding IWV trends obtained from nearby homogenized radiosonde data. The correlation coefficient of the trends increases significantly by 38% after using the homogenized GPS data.