Earthquake Rupture Area Research Articles

On the use of tsunami-source data for high-resolution fault imaging of offshore earthquakes

The source imaging for offshore earthquakes using terrestrial geodetic data has a limited estimation performance due to the low data resolution. One approach to overcome this limitation is the use of seafloor geodetic data. In this study, we focus on tsunami-source data, which is the spatial distribution of vertical crustal displacements above the source area and can be derived from tsunami waveform records. We evaluate how the use of this spatial seafloor geodetic data improves the estimation of a rectangular fault model. Here, the fault model of the 2016 off-Fukushima earthquake in Japan, which was a shallow intraplate earthquake (Mw 7.0), was estimated by three inversions: terrestrial Global Navigation Satellite System (GNSS) data only, tsunami-source data only, and a combination of the GNSS data and tsunami-source data. A Bayesian inversion approach was used to understand the distribution of the estimated fault parameters and their relationship. The results indicated that the terrestrial GNSS data have a low resolution for the analysis of the offshore earthquake, which resulted in a biased solution with large uncertainty. Conversely, the use of tsunami-source data significantly improved the resolution and reliability of source imaging and reduced the dependency among fault parameters. These results suggested that the high-spatial-resolution information of tsunami source is a powerful tool in source imaging of offshore shallow earthquakes. Moreover, the combined use of the two different geodetic data leads to a more robust estimation of fault parameters. We believe that the use of tsunami-source data is useful, not only for the post hoc source analysis, but also for estimating an earthquake rupture area just after a large earthquake, where GNSS data are currently used.Graphical

Analysis of the coseismic slip behavior for the MW = 9.1 2011 Tohoku-Oki earthquake from satellite GOCE vertical gravity gradient

Over the past decade, the three largest and most destructive earthquakes in recent history with associated tsunamis occurred: the Mw = 9.2 Sumatra-Andamam in 2004, then the Mw = 8.8 Maule in 2010, and finally the Mw = 9.1 Tohoku- Oki in 2011. Due to the technological and scientific developments achieved in recent decades, it has been possible to study and model these phenomena with unprecedented resolution and precision. In addition to the coseismic slip models, for which joint inversions of data from various sources are carried out (e.g., teleseismic data, GNSS, INSAR, and Tsunami, among others), depicting the space-time evolution of the rupture, we have high-resolution models of the degree of interseismic coupling (based on GNSS) and also maps of seismic b-value changes. Among these advances, new Earth gravity field models allow mapping densities distribution homogeneously and with a resolution (in wavelengths) of approximately the large rupture areas of great megathrust earthquakes. In this regard, the maximum resolution of GOCE-derived static models is in the order of λ/2≈66 km, while GRACE monthly solutions are in the order of λ/2≈300 km. From the study of the static and dynamic gravitational field, it has been possible to infer mass displacements associated with these events, which have been modeled and compared to the deformation inferred using other methods, yielding very good results. In this work we study the kinematic behavior of the rupture process for one of these largest events, the Mw = 9.1 Tohoku-Oki 2011 earthquake, employing the vertical gradient of gravity derived from the GOCE satellite, finding that the maximum slip occurred close to a lobe of minimum Tzz, as was observed for other case-studies in other subduction-related settings studied in previous works (e.g., the Maule earthquake and the Sumatra-Andaman earthquake, among others). In addition, from the rupture propagation using kinematic models, it can be observed that the rupture is arrested when it approaches high-density structures and, it is enhanced when connecting with lobes of low vertical gravity gradient. We also mapped a block expressed as a low Tzz lobe, developed along the marine forearc, which is controlled by a parallel-to-the-trench normal fault that accommodates subsidence during the interseismic period, as it is coupled with the subducted slab. Then, after rupturing the plate interface, this block is decoupled promoting tectonic inversion and uplift. In this way, the hypothesis that the density structure along the forearc is the ultimate first-order factor that governs the rupture process is reinforced.

Mechanically Coupled Areas on the Plate Interface in the Nankai Trough, Japan and a Possible Seismic and Aseismic Rupture Scenario for Megathrust Earthquakes

AbstractQuantifying the stress distribution or finding mechanically coupled areas on the plate interface is fundamentally important for conjecturing megathrust earthquakes that may occur in the future. Kinematically coupled areas, or slip deficit distributions, on plate interfaces were commonly estimated by geodetic‐data analyses. However, mechanically coupled areas are not identical to the kinematically coupled areas. The present study develops an inversion method to estimate the stress rate distribution as mechanically coupled areas. We apply this method to the Nankai trough subduction zone, southwestern Japan to detect mechanically coupled areas. Some of the estimated coupled areas correspond to the rupture areas of historical earthquakes. Others are in a deeper part, which may release the stress as aseismic slip. We then construct a rupture scenario that can occur in the Nankai trough in the future based on the estimated mechanical coupling, assuming that an effective stress accumulation period is 100 years. The scenario consists of a foreshock of MW 8.0 followed by an afterslip of MW 7.7 and a mainshock of MW 8.2. Although the moment magnitude of the afterslip is slightly smaller than the foreshock, the strain energy released by the foreshock is significantly larger than the afterslip. In order to compare events when the stress drops are significantly different, it is useful to use the minimum released strain energy, or available energy, to characterize the magnitude of the events.

Upper-plate structure and tsunamigenic faults near the Kodiak Islands, Alaska, USA

AbstractThe Kodiak Islands lie near the southern terminus of the 1964 Great Alaska earthquake rupture area and within the Kodiak subduction zone segment. Both local and trans-Pacific tsunamis were generated during this devastating megathrust event, but the local tsunami source region and the causative faults are poorly understood. We provide an updated view of the tsunami and earthquake hazard for the Kodiak Islands region through tsunami modeling and geophysical data analysis. Using seismic and bathymetric data, we characterize a regionally extensive seafloor lineament related to the Kodiak shelf fault zone, with focused uplift along a 50-km-long portion of the newly named Ugak fault as the most likely source of the local Kodiak Islands tsunami in 1964. We present evidence of Holocene motion along the Albatross Banks fault zone, but we suggest that this fault did not produce a tsunami in 1964. We relate major structural boundaries to active forearc splay faults, where tectonic uplift is collocated with gravity lineations. Differences in interseismic locking, seismicity rates, and potential field signatures argue for different stress conditions at depth near presumed segment boundaries. We find that the Kodiak segment boundaries have a clear geophysical expression and are linked to upper-plate structure and splay faulting. The tsunamigenic fault hazard is higher for the Kodiak shelf fault zone when compared to the nearby Albatross Banks fault zone, suggesting short wave travel paths and little tsunami warning time for nearby communities.

Marine forearc structure of eastern Java and its role in the 1994 Java tsunami earthquake

Abstract. We resolve a previously unrecognized shallow subducting seamount from a re-processed multichannel seismic profile crossing the 1994 Mw 7.8 Java tsunami earthquake rupture area. Seamount subduction occurs where the overriding plate experiences uplift by lateral shortening and vertical thickening. Pronounced back-thrusting at the landward slope of the forearc high and the formation of splay faults branching off the landward flank of the subducting seamount are observed. The location of the seamount in relation to the 1994 earthquake hypocentre and its co-seismic slip model suggests that the seamount acted as a seismic barrier to the up-dip co-seismic rupture propagation of this moderate-size earthquake.

Frictional Segmentation of the Chilean Megathrust From a Multivariate Analysis of Geophysical, Geological, and Geodetic Data

AbstractWe present a formal seismic segmentation model of the Chile Margin derived from a multivariate analysis of gravity anomalies, basal friction and interplate locking from GPS velocities. These independent proxies of the frictional structure along the megathrust integrate over timescales of 101–107 years. We analyze their covariance using Principal Component Analysis (PCA) and extracting the Empirical Orthogonal Functions (EOFs) and Principal Components (PCs). We find that (a) the first and first + second PCA modes for the triple combination of the proxies explain (respectively) 41.4% and 69.4% of the total variance and (b) the first‐order spatial pattern of each isolated proxy is largely reproduced when it is combined with the other two proxies. This suggests that the three proxies share a first order along‐trench tendency that must be connected through a common physical process. Comparing our results with rupture areas of earthquakes and using the rate‐and‐state friction theory, we infer that trench‐perpendicular PCs reflect the limits of the velocity‐weakening (VW) seismogenic zone. Polarity changes in along‐trench EOFs (mostly for gravity) mark changes between unstable and conditionally stable frictional regimes inside the VW zone and we used them to define the boundaries of 17 unitary segments. These segments correlate with first‐order tectonic features of the Andes at >103‐km scales as well as with zones of multi‐segment ruptures at 102‐km scale. We analyze the combined influence of subducting and upper‐plate geologic features on the nature of persistent seismic barriers. Our results have implications for understanding seismotectonic processes along the Andean margin and elsewhere.

Scaling relations between seismic moment and rupture area of earthquakes in stable continental regions

This article describes the development of scaling relations between seismic moment and rupture area of earthquakes in stable continental regions (SCRs). The article reviews the relations developed by Somerville and compares them with relations developed by other investigators. It also compares the scaling relations of SCR earthquakes with those in tectonically active continental regions (TCRs). Three different methods of estimating rupture area, based on aftershocks, slip models, and duration methods were used by Somerville to analyze the relation between seismic moment and rupture area, using earthquake source parameters compiled from published literature. For each category of data, the relations obtained were not significantly different from those obtained by constraining them to be self-similar (scale-invariant), so self-similar relations were adopted. The stress drops corresponding to these scaling relations range from 51 to 86 bars, with an average of 65 bars. This value is comparable with the value of 58 bars obtained by Leonard, and it is recommended that the Leonard scaling relations for SCR earthquakes be used for the NGA East Project. To a first approximation, the results of Somerville and those of Somerville et al. indicate that the rupture areas of SCR earthquakes are about half those of TCR earthquakes, and their stress drops are about 2.8 times higher. Allmann and Shearer find less of a difference, presumably because their intraplate category includes some earthquakes that we would assign to TCR instead of SCR. Their study indicates that the rupture areas of intraplate earthquakes are about two-thirds those of TCR earthquakes, and their stress drops are about 2 times higher.

Source model for the Mw 6.0 earthquake in Jiashi, China on 19 January 2020 from Sentinel-1A InSAR data

On January 19, 2020, an Mw 6.0 earthquake occurred in Jiashi, Xinjiang Uygur Autonomous Region of China. The epicenter was located at the basin-mountain boundary between the southern Tian Shan and the Tarim Basin. Interferometric Synthetic Aperture Radar (InSAR) is used to obtain the coseismic deformation field from both ascending and descending Sentinel-1A satellite images of the European Space Agency. The results showed that the coseismic deformation is distributed between the Kalping fault and the Ozgertaou fault. The earthquake produced significant deformation over an area of approximately 40 km by 30 km. The maximum and minimal displacements along the line of sight (LOS) are 5.3 cm and − 4.2 cm for the ascending interferogram and are 7.2 cm and − 3.0 cm for the descending interferogram, respectively. The fault geometry from the Multi peak Particle Swarm Optimization computation indicates that the seismogenic fault is a shallow low-dipping planar fault that is 4.58 km depth underground. The finite slip model inverted by the Steepest Descent Method implies that the rupture is dominated by a thrust fault. The slips are concentrated in a depth of 5–7 km with a maximum slip of 0.29 m. The estimated total seismic moment is 1.688 × 1018 Nm, corresponding to a magnitude of Mw 6.1. The seismogenic fault is the Kalping fault which has a listric structure. The coseismic deformation only occurred on the décollement layer and did not involve the ramp segment. The coseismic Coulomb stress changes have enhanced the stress on the deep margin of the Jiashi earthquake rupture area, indicating that there is still the possibility of strong earthquakes in this region in the future.

Segmentation of the Main Himalayan Thrust Illuminated by Bayesian Inference of Interseismic Coupling

AbstractWe use a recent compilation of geodetic data of surface displacements in a fully Bayesian approach to derive a probabilistic estimate of interseismic coupling along the Main Himalayan Thrust (MHT). Our probabilistic estimate of interseismic coupling highlights four large, highly coupled patches separated by three potential barriers of low coupling. Locked patches overlap with estimated rupture areas of historical large earthquakes over the past centuries. The coincident spatial variability in coupling, seismicity, and prominent active topography suggests a structural segmentation of the MHT imposed by inherited tectonic structures from the India‐Eurasia collision. This correlation implies that inherited tectonic structures may affect how stress builds up along the MHT, thereby influencing the location and size of large Himalayan earthquakes and the growth of the mountain range.

Empirical relocation of large subduction-zone earthquakes via the teleseismic network correlation coefficient method

The network correlation coefficient method is a powerful tool for relocating events under low signal-to-noise ratio conditions, as well as small (M_{text{W}} < 5.0) near-field earthquakes, where the relocation is performed by maximizing the sum of the waveform correlation coefficients between pairs of events in the network. Here we extend the method to relocating large (6.0 le M_{text{W}} le 7.3) teleseismic subduction-zone events. We apply a source duration correction to better correlate different-magnitude events. Four case studies are presented where we precisely locate interplate events near the rupture areas of giant megathrust earthquakes, namely the 2011 M_{text{W}} ;9.1 , Tohoku, 2010 M_{text{W}} ;8.8 Maule, 2005 M_{text{W}} ; 8.6 Nias, and 2007 M_{text{W}} ;8.5 Bengkulu earthquakes. We demonstrate the extensive links between these interplate events, with a clear alignment of earthquakes observed along each subduction interface.

Tsunamigenic Splay Faults Imply a Long‐Term Asperity in Southern Prince William Sound, Alaska

AbstractCoseismic slip partitioning and uplift over multiple earthquake cycles is critical to understanding upper‐plate fault development. Bathymetric and seismic reflection data from the 1964 Mw9.2 Great Alaska earthquake rupture area reveal sea floor scarps along the tsunamigenic Patton Bay/Cape Cleare/Middleton Island fault system. The faults splay from a megathrust where duplexing and underplating produced rapid exhumation. Trenchward of the duplex region, the faults produce a complex deformation pattern from oblique, south‐directed shortening at the Yakutat‐Pacific plate boundary. Spatial and temporal fault patterns suggest that Holocene megathrust earthquakes had similar relative motions and thus similar tsunami sources as in 1964. Tsunamis during future earthquakes will likely produce similar run‐up patterns and travel times. Splay fault surface expressions thus relate to plate boundary conditions, indicating millennial‐scale persistence of this asperity. We suggest structure of the subducted slab directly influences splay fault and tsunami generation landward of the frontal subduction zone prism.

The December 2nd, 2015 Bingöl, Eastern Anatolia-TURKEY, earthquake (Mw = 5.3): A rupture on optimally oriented fault plane

Eastern Anatolia is a ∼2 km high plateau shaped by the continent-continent collision of the Arabian and Eurasian plates. The left lateral East Anatolian Fault Zone, the right lateral North Anatolian Fault Zone and the Bitlis fold-thrust belt are the major tectonic boundaries of this convergence zone. The Arabian-Eurasian collision has resulted in high volcanism and well recorded seismic activity in Eastern Anatolia and its surroundings. Karlıova Junction is located at the intersection of these major fault systems and contains secondary faults such as the Karakoçan fault and the Sancak-Uzunpınar fault. We processed high quality waveform data collected from a recent seismic activity in the proximity of the city of Bingöl where the largest event is the moderate-size earthquake (Mw = 5.3) occurred on 2 December 2015. This event is located to the west of Karlıova Junction and to the northwest of Sancak-Uzunpınar fault. The spatial distribution of the aftershocks points out unmapped faults with NNW-SSE alignment towards the west of Sancak-Uzunpınar fault. Aftershock depth distribution indicates a nearly 15 km deep brittle seismogenic zone. The relocated aftershock distributions and seismic moment calculations yield a rupture area 9 km in length and 5 km in width with an average 8 cm of slip. The latest four earthquakes preceding the 2015 Eq with magnitudes larger than Mw > 6 enhanced the Coulomb stress failure in the 2015 Bingöl Earthquake rupture area. The ruptured fault plane is in-line with the optimally oriented right-lateral strike-slip faults.

Spatial and Temporal Variations in Earthquake Stress Drop on Gofar Transform Fault, East Pacific Rise: Implications for Fault Strength

AbstractOn Gofar Transform Fault on the East Pacific Rise, the largest earthquakes (6.0 ≤ MW ≤ 6.2) have repeatedly ruptured the same portion of the fault, while intervening fault segments host swarms of microearthquakes. These long‐term patterns in earthquake occurrence suggest that heterogeneous fault zone properties control earthquake behavior. Using waveforms from ocean bottom seismometers that recorded seismicity before and after an anticipated 2008 MW 6.0 mainshock, we investigate the role that differences in material properties have on earthquake rupture at Gofar. We determine stress drop for 138 earthquakes (2.3 ≤ MW ≤ 4.0) that occurred within and between the rupture areas of large earthquakes. Stress drops are calculated from corner frequencies derived using an empirical Green's function spectral ratio method, and seismic moments are obtained by fitting the omega‐square source model to the low frequency amplitude of the displacement spectrum. Our analysis yields stress drops from 0.04 to 3.2 MPa with statistically significant spatial variation, including ~2 times higher average stress drop in fault segments where large earthquakes also occur compared to fault segments that host earthquake swarms. We find an inverse correlation between stress drop and P wave velocity reduction, which we interpret as the effect of fault zone damage on the ability of the fault to store strain energy that leads to our spatial variations in stress drop. Additionally, we observe lower stress drops following the MW 6.0 mainshock, consistent with increased damage and decreased fault strength after a large earthquake.

Development and application of a tsunami fragility curve of the 2015 tsunami in Coquimbo, Chile

Abstract. The last earthquake that affected the city of Coquimbo took place in September 2015 and had a magnitude of Mw=8.3, resulting in localized damage in low-lying areas of the city. In addition, another seismic gap north of the 2015 earthquake rupture area has been identified; therefore, a significant earthquake (Mw=8.2 to 8.5) and tsunami could occur in the near future. The present paper develops a tsunami fragility curve for the city of Coquimbo based on field survey data and tsunami numerical simulations. The inundation depth of the 2015 Chile tsunami in Coquimbo was estimated by means of numerical simulation with the Non-hydrostatic Evolution of Ocean WAVEs (NEOWAVE) model and five nested grids with a maximum grid resolution of 10 m. The fragility curve exhibited behavior similar to that of other curves in flat areas in Japan, where little damage was observed at relatively high inundation depths. In addition, it was observed that Coquimbo experienced less damage than Dichato (Chile); in fact, at an inundation depth of 2 m, Dichato had a ∼75 % probability of damage, while Coquimbo proved to have only a 20 % probability. The new fragility curve was used to estimate the damage by possible future tsunamis in the area. The damage assessment showed that ∼50 % of the structures in the low-lying area of Coquimbo have a high probability of damage in the case of a tsunami generated off the coast of the study area if the city is rebuilt with the same types of structures.

InSAR-Constrained Interseismic Deformation and Potential Seismogenic Asperities on the Altyn Tagh Fault at 91.5–95°E, Northern Tibetan Plateau

The present-day kinematic features of the different segments of the Altyn Tagh Fault (ATF) have been well-studied using geodetic data. However, on the eastern segment of the ATF at 91.5–95°E, high spatial resolution deformation has not been previously reported. Here, we processed 185 interferometric synthetic aperture radar (InSAR) images from three descending tracks of the C band ERS-1/2 and Envisat satellites spanning 1995–2011 and obtained the average deformation velocity field. Results show a left-lateral motion of ~4 mm/year along the fault-parallel direction across the ATF at 91.5–95°E, which is consistent with Global Positioning System (GPS) observations. The slip deficit rate distribution at shallow depths was resolved through the InSAR deformation velocity using a discretized fault plane. The slip deficit is capable of an Mw 7.9 earthquake, considering the elapsed time of the latest M 7.0 event. Two potential asperities that could be nucleation sites or rupture areas of future earthquakes were delineated based on the coupling coefficient and seismicity distributions along the fault plane. The larger asperity extends more than 100 km along the ATF at depths of 8–12 km. Our InSAR observations support the undeformed blocks model of the Indo-Eurasian collisional mechanism at the northern margin of the Tibetan plateau.

Directly estimating earthquake rupture area using second moments to reduce the uncertainty in stress drop

Author Posting. © The Author(s), 2018. This article is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 214 (2018): 2224–2235, doi:10.1093/gji/ggy201.

Dynamic triggering and earthquake swarms on East Pacific Rise transform faults

AbstractWhile dynamic earthquake triggering has been reported in several continental settings, offshore observations are rare. Oceanic transform faults share properties with continental geothermal areas known for dynamic triggering: high geothermal gradients, high seismicity rates, and frequent swarms. We study dynamic triggering along the East Pacific Rise by analyzing 1 year of seismicity recorded by Ocean Bottom Seismographs. By comparing the response to teleseismic waves from global earthquakes, we find triggering to be most sensitive to changes in normal stress and to preferentially occur above 0.25 kPa. The clearest example of triggering occurs on the Quebrada and Gofar faults after the Mw8.0 Wenchuan earthquake. On Gofar, triggered seismicity occurs between the rupture areas of large earthquakes, within a zone characterized by aseismic slip, abundant microseismicity, frequent swarms, and low Vp. We infer that lithological properties inhibiting rupture propagation, such as high porosity and fluid content, also favor dynamic triggering.

Afterslip following the 2007 Mw 8.4 Bengkulu earthquake in Sumatra loaded the 2010 Mw 7.8 Mentawai tsunami earthquake rupture zone

AbstractThe 12 September 2007 Mw 8.4 Bengkulu earthquake in Sumatra marked the first in a modern series of large earthquakes along the Mentawai section of the Sunda megathrust. Understanding the spatial distribution of coseismic slip and ensuing afterslip is important for assessing seismic hazard in neighboring unruptured regions of the megathrust. We reestimate the spatial distribution of coseismic slip during this earthquake with improved coseismic offsets from the Sumatran GPS Array (SuGAr) and estimate afterslip following this earthquake with SuGAr postseismic time series spanning ∼6.3 years after the earthquake. We invert for the spatiotemporal distribution of afterslip with the principal component analysis‐based inversion method (PCAIM), and we take into account viscoelastic deformation by incorporating into the inversion the estimation of strain within ductile deforming blocks located at asthenospheric depths. Our results suggest cumulative afterslip concentrated within, updip, and downdip of the 2007 coseismic rupture area and shallow afterslip that borders and overlaps the 2010 Mw 7.8 Mentawai earthquake rupture zone. The cumulative contribution of stress changes due to the coseismic event and the ensuing afterslip likely increased strain rates in the shallow portion of the megathrust adjacent to the Mentawai earthquake rupture area, potentially promoting its rupture in 2010.

Interseismic deformation at subduction zones investigated by 2D numerical modeling: case study before the 2010 Maule earthquake

We study interseismic deformation preceding the Mw8.8 2010 Maule earthquake by means of two-dimensional finite-element modeling. Our goal is to gain insight into the fundamental factors controlling elastic strain build-up and release in subduction zones, and to evaluate different modeling approaches of surface displacement as observed by GPS. We developed a linear elasticity solver that allows us to implement a realistic subducting plate geometry constrained by geophysical data. We test the influence of subducting plate thickness, variations in the updip and downdip limit of a 100% locked interplate zone, elastic parameters, and velocity reduction at the base of the subducted slab. We compared our modeled predictions with interseismic GPS observations along an EW profile crossing the Maule earthquake rupture area, in order to determine best fitting parameters. Our results indicate little influence of the subducting plate thickness at a given downdip limit, which itself has a strong influence on surface deformation. However, the fit to observations is achieved only after reducing the velocity at the base of the subducted slab below the trench region to 10% of the far-field convergence rate. We link this novel result to complementary numerical models that gradually evolve toward considering longer time-scales and complex rheology in order to evaluate the mechanical meaning of the above mentioned inferred kinematic conditions. This allowed us to link the velocity reduction at the base of subducting slabs with a long-term state of high flexural stress resulting from the mechanical interaction of the slab with the underlying mantle. Even a small amount of theses high deviatoric stresses may transfer towards the upper portion of the slab as strain energy that could participate into the mechanical loading of the megathrust and therefore in triggering large earthquakes there.

Localised Effects of a Mega-Disturbance: Spatiotemporal Responses of Intertidal Sandy Shore Communities to the 2010 Chilean Earthquake.

Determining the effects of unpredictable disturbances on dynamic ecological systems is challenged by the paucity of appropriate temporal and spatial coverage of data. On 27 February 2010, an 8.8 Mw mega-earthquake and tsunami struck central Chile and caused coastal land-level changes, massive damage to coastal infrastructure, and widespread mortality of coastal organisms. Wave-exposed sandy beaches showed significant changes of species abundances from before to after the earthquake, but the highly dynamic biotic and abiotic conditions of these habitats make difficult to draw clear-cut conclusions from these patterns. Here, we analysed a beyond-BACI (Before-After Control-Impact) sampling design to test whether the effects of the Maule earthquake on sandy-shore species diversity, abundance, and structure were heterogeneous along the shore. Invertebrate species abundances were quantified before (i.e. February 2010) and after (i.e. March 2010, September 2010, and March 2011) the earthquake at three sandy shores randomly located within the earthquake rupture area and three sites within a “control” area located >400 km southward from epicentre. Immediately after the earthquake took place, the three sites located in the rupture area showed anomalous beach-profile uplifts that did not comply with the erosion (i.e. “negative” uplifts) that regularly occurs during late summer in the region. Species richness, abundance, and community structure significantly varied from before to after the strike, but these patterns of change varied among sites within both areas. Only the site with the strongest and persistent beach-profile uplift within the rupture area showed significant concomitant changes in species richness and community structure; after 13 months, this community showed a similar multivariate structure to the before-disturbance state. This site, in particular, was located in the section of the rupture area that received most of the impact of the after-earthquake tsunami. Therefore, our results suggest that the effects of the Maule mega-earthquake on the ecological communities were spatially heterogeneous and highly localised. We suggest that high mobility and other species’ adaptations to the dynamic environmental conditions of sandy beaches might explain the comparatively high resilience of these assemblages. With this work we hope to motivate further experimental research on the role of individual- and population-level properties in the response of sandy-beach communities to interacting sources of disturbances.