Several major earthquakes (Mw>7) have occurred in this gap since 1850 (Fig. 1); the largest until now was the Mw 7.7 Tocopilla earthquake in 2007, which broke the southern rim of this segment beneath and north of Mejillones Peninsula along a total length of 150 km. Only the downdip end of the locked zone slipped in this event, and the total slip in the rupture area was less than 2.6 m 6,7 leaving most of the past slip deficit of c. 8-9 m accumulated since 1877 3 approaches. First, we performed waveform modelling of local strong motion seismograms and teleseismic body waves to constrain the kinematic development of the rupture towards the final displacement in a joint inversion with continuous GPS data of static displacements (Fig. 1, 2a). Second, we use the backprojection technique applied to stations in North America to map the radiation of high frequency seismic waves (HFSR; 1-4 Hz) 9,10 . The latter technique is not sensitive to absolute slip amplitudes, but rather to changes in slip and rupture velocity.During the first 35-40s the rupture propagated downdip with increasing velocity, nearly reaching the coastline (Fig. 2a,b). Surprisingly, towards the end of the rupture, the area near the epicenter was reactivated. In spite of the relatively complicated kinematic history of the rupture the cumulative slip shows a simple 'bull's eye' pattern with a peak coseismic slip of (Fig. 3a). The Iquique main shock nucleated at the 4 northwestern border of a locked patch and ruptured towards its center (Fig. 2a, 3a). The downdip end of the main shock as well as for the large Mw 7.6 aftershock rupture mapped both by the HFSR and co-seismic slip agrees quite accurately with the downdip end interseismic coupling (Fig. 2a,c 3a). The accelerated downdip rupture propagation for both earthquakes closely followed the gradient towards higher locking. Therefore, the Iquique event and its largest aftershock appear to have broken the central, only partly locked segment of the Northern Chile Southern Peru seismic gap releasing part of the slip deficit accumulated here since 1877 (cf. Fig. 1).The seismicity before the Iquique earthquake also concentrates in this zone of intermediate locking at the fringe of the highly locked -high slip patch (Fig. 3a). Starting in July 2013, three foreshock clusters with increasingly larger peak magnitudes and cumulative seismic moment occurred here (Fig. 2c, 3a,c). The mainshock rupture started at the northern end of the foreshock zone, inside the region of intermediate locking (Fig. 2c, 3a). Interestingly, the second foreshock cluster (January 2014) is associated with a weak transient deformation, whereas the third cluster (March 2014) shows a very distinct transient signal. GPS displacement vectors calculated over the times spanning these foreshock clusters point towards the cluster epicentres (Extended Data Figure 4). Deformation for both transients is entirely explained by the cumulative coseismic displacement of the respective foreshock clusters (Fig. 3d inset, Extended Data Figure 4). The ar...

40Large volcanic eruptions on Earth commonly occur with collapse of the roof of a crustal magma 41 reservoir, forming a caldera. Only a few such collapses occur per century and lack of detailed 42 observations has obscured insight on mechanical interplay between collapse and eruption. We use Calderas are 1 -100 km diameter depressions found in volcanic regions of Earth and other planets. basaltic andesite) intrusive activity and eruptions (2,(9)(10)(11)(12). 59The consensus from field and modelling studies is that caldera collapse progresses from initial 60 surface downsag to fault-controlled subsidence (1, 8, 13, 14). The pre-collapse topography is obtained by subtracting the subsidence observed at the surface. As we recorded the caldera subsidence mainly on the ice (Fig. 1, Fig. S1), we made corrections and (Fig. 3A). We therefore conclude that suggestions of a large increase in ice flow out of the caldera 147 during these events (25) cannot be fitted with our data. 148Bedrock subsidence exceeding 1 m occurred within an area of 110 km 2 that extended beyond the 149 pre-existing caldera (Fig. 1, Fig. S1). After termination of collapse the total subsidence at the pre-150 existing caldera rims amounted to 3 to 11 meters ( Fig. 1D and 1E). Using subglacial radio-echo GPS station in the center of the caldera (Fig. 1A), including the rate of vertical rate of ice surface Cumulative number of M>4 caldera earthquakes, with magnitude evolution colored in red, blue and 176 grey representing clusters on the southern rim, the northern rim and smaller clusters, respectively 177 (see Fig. S5). E) Cumulative seismic moment for M>4 caldera earthquakes. from analysis of subaerial gas measurements (Fig. 4). This depth concurs with our regional on FTIR and Multi-GAS measurements (24). 194Seismicity and subsurface structure 195 We used seismic data and Distinct Element Method (DEM) numerical modelling (24), to 196 characterize the deeper collapse structure as the reactivation of a steeply-inclined ring fault (Fig. 5). 197We mostly observed seismicity at depths of 0-9 km beneath the northern and southern caldera rims 198( Fig. 5B), with earthquakes being more numerous on the northern rim. This spatial pattern of 199 seismicity is consistent with fracturing above a deflating magma reservoir that was elliptical in (Fig. 5C, D). Our best fitting models had preexisting faults dipping out at 80-85¡ from the caldera 207 center on the north side and at 85-90¡ toward the caldera center on the south side. The modeled pre- 208existing faults lay at 1-2 km below the surface on the north side and 3-4 km on the south side. 209Modeling of a more complex fault geometry or the inclusion of greater material heterogeneity may 210 further improve the data fit, but presently lacks robust geophysical constraints. components of the observed earthquakes at B ‡rdarbunga. We, however, narrowed down on 222 plausible solutions by using the micro-earthquakes (Fig. 5A). The moment tensor solutions are well 223 constrained, but the inferred d...

The moment magnitude () 5.5 earthquake that struck South Korea in November 2017 was one of the largest and most damaging events in that country over the past century. Its proximity to an enhanced geothermal system site, where high-pressure hydraulic injection had been performed during the previous 2 years, raises the possibility that this earthquake was anthropogenic. We have combined seismological and geodetic analyses to characterize the mainshock and its largest aftershocks, constrain the geometry of this seismic sequence, and shed light on its causal factors. According to our analysis, it seems plausible that the occurrence of this earthquake was influenced by the aforementioned industrial activities. Finally, we found that the earthquake transferred static stress to larger nearby faults, potentially increasing the seismic hazard in the area.

No abstract