Main Slip Patches Research Articles

Co-seismic and post-seismic slip associated with the 2021 Mw5.9 Arkalochori, Central Crete (Greece) earthquake constrained by geodetic data and aftershocks

The co-seismic and post-seismic deformation field associated with the Mw5.9 Arkalochori main shock that occurred in central Crete (Greece) on 27 September 2021 is analyzed using Copernicus Sentinel-1A & 1B images, GNSS measurements and seismological data. The fault geometry is constrained through the joint inversion of multiple datasets and the slip distribution for the co-seismic and post-seismic period is obtained using a homogeneous half-space elastic model and the Steepest Descent Method. The results indicate a blind normal fault striking 215° with a 55° dip to the northwest and the co-seismic slip model suggests a nearly circular main slip patch (8 × 6 km2) with a maximum slip of 0.98 m. Post-seismic displacements started rapidly after the main shock followed by a gradual decay as highlighted by the calculated InSAR time series. The temporal evolution of post-seismic slip is described by a simple logarithmic function, decaying faster at the southwest part of the fault. The cumulative afterslip model suggests that the maximum post-seismic slip of 0.23 m occurred within a similar depth range compared to the co-seismic one, yet with a shift towards the southwest. Post-seismic slip inside the main shock rupture area is sustained, highlighting the slow recovery of locking in the co-seismic slip region. Afterslip (seismic or aseismic) played a dominant role in the early post-seismic period acting complementarily to the main rupture. Indications suggest that the spatiotemporal evolution of the productive aftershock sequence may be driven afterslip, alongside other potential factors.

Revisiting the 2015Mw = 8.3 Illapel earthquake: unveiling complex fault slip properties using Bayesian inversion

SUMMARYThe 2015 moment magnitude Mw = 8.3 Illapel earthquake is the largest mega-thrust earthquake that has been recorded along the Chilean subduction zone since the 2010 Mw = 8.8 Maule earthquake. Previous studies indicate a rupture propagation from the hypocentre to shallower parts of the fault, with a maximum slip varying from 10 to 16 m. The amount of shallow slip differs dramatically between rupture models with some results showing almost no slip at the trench and other models with significant slip at shallow depth. In this work, we revisit this event by combining a comprehensive data set including continuous and survey GNSS data corrected for post-seismic and aftershock signals, ascending and descending InSAR images of the Sentinel-1A satellite, tsunami data along with high-rate GPS, and doubly integrated strong-motion waveforms. We follow a Bayesian approach, in which the solution is an ensemble of models. The kinematic inversion is done using the cascading capability of the AlTar algorithm, allowing us to first get a static solution before integrating seismic data in a joint model. In addition, we explore a new approach to account for forward problem uncertainties using a second-order perturbation approach. Results show a rupture with two main slip patches, with significant slip at shallow depth. During the rupture propagation, we observe two regions that are encircled by the rupture, with no significant slip, westward of the hypocentre. These encircling effects have been previously suggested by back-projection results but have not been observed in finite-fault slip models. We propose that the encircled regions correspond to zones where the yield stress largely exceeds the initial stress or where fracture energy is too large to be ruptured during the Illapel earthquake. These asperities may potentially break in the future and probably already broke in the past.

Joint Inversion of Geodetic and Strong Motion Data for the 2012, Mw 6.1–6.0, May 20th and May 29th, Northern Italy Earthquakes: Source Models and Seismotectonic Interpretation

AbstractWe present the first rupture models of the two mainshocks of the 2012 northern Italy sequence, determined by jointly inverting seismic and geodetic data. We aim at providing new insights into the mainshocks for which contrasting seismotectonic interpretations are proposed in literature. Sources' geometric parameters were constrained by seismic reflection profiles, 3‐D relocations and focal mechanisms of mainshocks/aftershocks. Site‐specific velocity profiles were used to model accelerograms affected by strong propagation effects related to the Po basin. Our source models differ significantly from previous ones relying on either seismic or geodetic data. Their comparison against geological sections and aftershock distribution provides new insights about the ruptured thrust faults. The May 20th Mw6.1 mainshock activated the Middle Ferrara thrust‐ramp dipping ∼45° SSW‐wards, breaking a main eastern slip patch 4–15 km deep in Mesozoic carbonates (maximum slip 0.7–0.8 m) and Paleozoic‐Triassic basement rocks, and a small western patch in the basement. The May 29th Mw6.0 mainshock featured two separated asperities along the Mirandola thrust‐ramp dipping ∼42° S‐wards: an eastern asperity 4–15 km deep in Mesozoic carbonates and basement rocks (maximum slip 0.7 m) and a deeper western one (7–16 km depth) mainly in the basement (slip peak 0.8 m). On‐fault aftershocks were concentrated within the basement and Mesozoic carbonates, devoiding high‐slip zones. Slip and aftershock distribution was controlled by the rheological transition between Mesozoic carbonates and Cenozoic sediments. Unlike previous thin‐skinned tectonic interpretations, our results point to a complex rupture process along moderately dipping (40°–45°) thrust‐ramps deeply rooted into the Paleozoic crystalline basement.

The Mechanisms of Tsunami Amplification and the Earthquake Source of the 2021 M 7 Acapulco, Mexico, Earthquake

ABSTRACT We present a slip model for the 2021 M 7 Acapulco, Mexico, earthquake produced by inversion of strong motion, Global Navigation Satellite Systems, tide gauge, and Interferometric Synthetic Aperture Radar data. The earthquake occurs within the Guerrero gap, identified as a region of concern for its seismogenic potential and paucity of large events. We find that rupture was compact, constrained to depths between 10 and 20 km, and consistent of two main slip patches. The slip model leaves a broad swath of the megathrust unbroken, and, whether the event signals a reactivation of large earthquakes in the region remains unknown. We find that tide gauge recordings inside Acapulco Bay for the M 7 1962 earthquake and the 2021 event are strikingly similar, thus we interpret this as weak evidence that 2021 is a repeat of 1962. We also produce a high-resolution hydrodynamic model of the resulting tsunami using the slip model as initial condition and place special emphasis in understanding the long duration (∼17 hr) of waves inside the bay. We find that simple bay resonance alone does not account for the features of the event. Rather it is a complex interaction with shelf modes and edge waves that continuously re-excite the bay resonance that leads to the protracted tsunami disturbances. Furthermore, we find that significant currents in excess of 1 m/s occur in localized portions of the bay even when wave amplitudes remain small.

An Atypical Shallow Mw 5.3, 2021 Earthquake in the Western Corinth Rift (Greece)

AbstractModerate‐to‐large earthquakes in rifts may occur on leading boundary faults or inner antithetic faults. Here we show a rare case of the 2020–2021 seismic sequence in the Corinth rift, that culminated in the shallow rupture of the antithetic fault, neither preceded nor followed by the leading fault rupture. The hypocenter of the largest shock (Mw 5.3 of 17 February 2021) was located at ∼8 km depth. However, seismic waveform data, supported by satellite‐geodetic and tide gauge measurements, pointed to rupture at shallow depth (∼3 km), where no earthquakes were previously observed. We show that the earthquake most probably ruptured two orthogonal, conjugate fault segments: a weak nucleation phase occurred in the microseismically highly active sub‐horizontal detachment layer, followed – a few seconds later – by a larger, shallow moment release on a high‐angle, south‐dipping normal fault. The latter is the Mornos offshore fault, antithetic to the leading, north‐dipping Psathopyrgos fault. Our study presents the first instrumental/observational evidence of a very shallow Mw 5+ event in this rift – and one of the few reported worldwide. The depth limit of the main shallow slip patch coincides with the expected crossing of the Mornos fault with the Psathopyrgos fault, stressing the importance of fault segmentation and rooting inherited from the rift history. This unusual shallow slip in a depth range with little background seismicity and few aftershocks needs to be further investigated by dynamic modeling as a possible prototype of hazardous events in rift environments.

The Gorgona island, Colombia, earthquake of 10 September 2007 (Mw 6.8); rupture process and implications on the seismic hazard in the region

The subduction zone in southern Colombia and northern Ecuador has been the site of significant events in the area that broke the great 1906 earthquake in Ecuador. The 2007 earthquake on Gorgona Island is located within this area, presenting a normal-type focal mechanism and a moment magnitude Mw 6.8. This event is due to the compressive stresses exerted on the continental plate of South America by the Nazca oceanic plate. We performed a kinematic inversion source process of the Gorgona earthquake. Our results show the main slip patch at the hypocenter and a secondary slip release deeper at SW on the fault plane, which agrees with the location of aftershocks; the source time function obtained reflects these energy releases with a total duration of 15 s. We also analyzed the stress state field in this region using the focal mechanism of its early aftershock sequence, obtaining an extensive regime with an azimuth Shmin of 108° and NW-SE direction. The average interseismic coupling in this convergent margin zone is not so high as in the southern areas. Nevertheless, observations of this kind of seismicity and GPS studies in the upper plate of a subduction zone are essential in understanding earthquake cycles and, perhaps, in anticipating slip distributions in future subduction events.

The Crete Isl. (Greece) Mw6.0 Earthquake of 27 September 2021: Expecting the Unexpected

The 27 September 2021 damaging mainshock (Mw6.0) is the first known strong earthquake that ruptured the Arkalochori area, Crete Isl., Greece, during the entire historical period, making it an unexpected event in the long-term sense. The area is characterized by the presence of the normal active Kastelli Fault (KF) striking NNE-SSW and dipping towards ~WNW. The KF, of surface exposure only ~6 km, at its southern tip is truncated by the nearly perpendicular active Nipiditos fault. The main shock was preceded by foreshock activity lasting for ~3.9 months, thus the mainshock turned out to be an expected event in the short-term sense. Maximum ground subsidence of ~20 cm was estimated from InSAR images, but this also incorporates deformation that may have been caused by the largest aftershock (Mw5.1) of 28 September 2021. The fault model produced from the inversion of InSAR observations indicated strike 216°, dip towards ~NW at angle 53°, rake −95°, and is consistent with fault-plane solutions obtained from routine moment tensor analysis. The geodetic seismic moment calculated from the Okada’s formalism is 1.14 × 1018 N·m (Mw6.0), while a maximum slip of 1.03 m was found at depths from 3.5 km to 5 km. The entire aftershock epicenters cloud strikes in a ~SW-NE direction but is distributed in two clusters, the southern and the northern ones. The foreshock cloud, the main slip patch, the deformation area, and the strongest aftershocks all fall within the southern cluster. The foreshocks concentration at the deepest edge of the main slip patch was a foreshadow of the mainshock nucleation area. The northern cluster, which is very likely due to the gradual expansion of aftershocks, is situated in the KF hanging wall block. To interpret the main seismic slip in the southern cluster area we propose the existence of a buried KF segment at the SSW-wards prolongation of the emerged at the surface segment. Assuming a rectangular seismic fault stress drop Δσ~7 bars was found. However, for a circular fault area, which in this case is more realistic, we get Δσ = 55 bars. This is a relatively large value for Greek earthquakes but is explainable by increased fault rigidity as a result of the long repeat time of strong earthquakes in KF.

Insights on fault reactivation during the 2019 November 11, Mw 4.9 Le Teil earthquake in southeastern France, from a joint 3-D geological model and InSAR time-series analysis

SUMMARY The 2019, Mw 4.9 Le Teil earthquake occurred in southeastern France, causing substantial damage in this slow deforming region. Field observations, remote sensing and seismological studies following the event revealed that coseismic slip concentrates at shallow depth along a ∼5 km long rupture associated with surface breaks and a thrusting mechanism. We further investigate this earthquake by combining geological field mapping, 3-D geology, InSAR time-series analysis and a coseismic slip inversion. From structural, stratigraphic and geological data collected around the epicentre, we first produce a 3-D geological model of the region surrounding the rupture using the GeoModeller software. Our model includes the geometry of the geological layers and the main faults, including the La Rouvière Fault, (LRF) the Oligocene normal fault that ruptured during the earthquake. We generate a time-series of surface displacement from Sentinel-1 SAR data ranging from early 2019 January to late 2020 January using the NSBAS processing chain. The spatio-temporal patterns of surface displacement for this time span show neither a clear pre-seismic signal nor significant post-seismic transient deformation. We extract the coseismic displacement pattern from the InSAR time-series, highlighting along-strike variations of coseismic surface slip. The maximum relative displacement along the line of sight is up to ∼16 cm and is located in the southwestern part of the rupture. We invert for the slip distribution on the fault from the InSAR coseismic surface displacement field. Constraining our fault geometry from the geological model, acceptable fault dip ranges between 55° and 60°. Our model confirms the reactivation of LRF, with reverse slip at very shallow depth and two main slip patches reaching, respectively, 30 and 24 cm of slip, both around 500 m depth. We finally discuss how the 3-D fault geometry and geological structure may have impacted the slip distribution and propagation during the earthquake. This study is a step to reassess the seismic hazard of the many faults similar to the La Rouvière one along the Cévennes fault system, in a densely populated area hosting several sensitive nuclear sites.

The 22 July 2020 MW 7.8 Shumagin seismic gap earthquake: Partial rupture of a weakly coupled megathrust

The earthquake potential of the Shumagin seismic gap along the Alaska Peninsula (∼162°W to ∼158.5°W) has been debated for more than 40 years. On 22 July 2020, the eastern half of the gap hosted an MW 7.8 earthquake involving a patchy rupture of the megathrust in the depth range of 20 to 45 km. The space-time slip distribution is determined by joint inversion of teleseismic P and SH waves and static displacements from regional GPS stations. The event initiated near the epicenter of the 10 November 1938 (MW 8.2) event, and ruptured westward, with little/no overlap with the 1938 rupture zone. The main slip patch has peak slip of ∼3.8 m below the Shumagin Islands, and produced ∼30 cm uplift and ∼25 cm SSE horizontal displacement on Chernabura Island. The slip model predicts well the small (<1 cm) tsunami signals persisting for more than ten hours observed at deep-water DART seafloor pressure recordings along the Alaska-Aleutian arc. Aftershocks with depths from 20 to 40 km fringe the large-slip patches, and show westward concentration during the first month after the mainshock. Aftershocks up-dip of the 1948 MW 7.1 event contribute to the high level of modest-size background seismicity extending to the trench in the region of very low seismic coupling (0.0–0.1) in the western Shumagin gap east of the 1 April 1946 (MW 8.6) rupture zone. The 31 May 1917 event is the last major earthquake to rupture the eastern half of the Shumagin gap, and has a lower surface wave magnitude (MSG-R 7.4, horizontal components) compared to the 2020 event (MSG-R 7.7, vertical components). Comparison of instrument-equalized waveforms for the 1917 and 2020 events indicates similar size contrast and differences in overall rupture duration and slip complexity. The 2020 rupture has average slip of ∼1.9 m over the 3600 km2 region with co-seismic slip ≥1 m. This is much less than the ∼6.7 m of potentially accumulated slip deficit since 1917, consistent with geodetic estimates of low average seismic coupling coefficient of 0.1–0.4. The megathrust seaward of the 2020 event has low seismicity and may either be aseismic or capable of comparable size ruptures. Comparisons are made with other subduction zones that have experienced relatively deep megathrust slip in regions with moderate seismic coupling.

Coseismic Slip Distribution of the 24 January 2020 Mw 6.7 Doganyol Earthquake and in Relation to the Foreshock and Aftershock Activities

AbstractOn 24 January 2020 (UTC), a destructive Mw 6.7 earthquake struck the east Anatolian fault of eastern Turkey after a series of foreshocks, causing many casualties and significant property damage. In this study, the rupture process of this earthquake is investigated with teleseismic broadband body-wave and surface-wave records. Results indicate that this earthquake is a left-lateral strike-slip event, and the rupture extends mainly to south. The main slip patch spreads ∼30 km along strike in the shallow above 14 km with a peak slip of ∼1.2 m, and the total seismic moment is 1.69×1019 N·m. The east–west component of horizontal surface displacement predicted with our slip model ranges from ∼0.4 to −0.3 m. The predicted displacements are consistent with the observed ones obtained from satellite images. We relocate 459 foreshocks and early aftershocks to explore the relationship between foreshock and aftershock sequences and coseismic slip. It is noted that there is an anticorrelation relationship between the distributions of early aftershocks and the coseismic slip. The strain energy in the large slip patch may have been sufficiently released by the mainshock; therefore, fewer early aftershocks occurred in that patch. Although we note a similar pattern between the relocated foreshock and coseismic slip, and a migration of foreshock, our dataset may not well resolve the correlation and migration due to the incomplete relocation foreshock catalog. Based on the slip model, we calculate the coulomb stress changes on the surrounding faults caused by the mainshock. The results reveal that the mainshock promoted stress accumulation on the northern and southern ends of the Elazig–Matalya segment and may reactivate the locked fault segment, leading to a high seismic risk in these regions. Although this earthquake does not significantly increase the coulomb stress change, the seismic risk of the Matalya–Kahraman Maras–Antakya segment should draw attention.

The 2018 Mw 6.8 Zakynthos (Ionian Sea, Greece) earthquake: seismic source and local tsunami characterization

SUMMARYWe investigated the kinematic rupture model of the 2018 Mw 6.8 Zakynthos, Ionian Sea (Greece), earthquake by using a non-linear joint inversion of strong motion data, high-rate GPS time-series and static coseismic GPS displacements. We also tested inversion results against tide-gauge recordings of the small tsunami generated in the Ionian Sea. In order to constrain the fault geometry, we performed several preliminary kinematic inversions by assuming the parameter values resulting from different published moment tensor solutions. The lowest cost function values were obtained by using the geometry derived from the United States Geological Survey (USGS) focal solution. Between the two conjugate USGS planes, the rupture model which better fits the data is the one with the N9°E-striking 39°ESE-dipping plane. The rupture history of this model is characterized by a bilateral propagation, featuring two asperities; a main slip patch extending between 14 and 28 km in depth, 9 km northeast from the nucleation and a slightly shallower small patch located 27 km southwest from the nucleation. The maximum energy release occurs between 8 and 12 s, when both patches are breaking simultaneously. The maximum slip is 1.8 m and the total seismic moment is 2.4 × 1019 Nm, corresponding to a Mw value of 6.8. The slip angle shows a dominant right-lateral strike-slip mechanism, with a minor reverse component that increases on the deeper region of the fault. This result, in addition to the observed possibility of similar mechanisms for previous earthquakes occurred in 1959 and 1997, suggests that the tectonic deformation between the Cephalonia Transform Fault Zone and the northern tip of the Hellenic Arc Subduction zone may be accommodated by prevailing right lateral low-dipping faults, occurring on re-activated structures previously experiencing (until Pliocene) compressional regime. Comparison of predicted and observed tsunami data suggests the need of a better characterization of local harbour response for this type of relatively short-wavelength events, which is important in the context of tsunami early warning. However, the suggested dominantly strike-slip character would in turn imply a reduced tsunami hazard as compared to a dominant thrust faulting regime from this source region.

Heterogeneous Behavior of the Campotosto Normal Fault (Central Italy) Imaged by InSAR GPS and Strong-Motion Data: Insights from the 18 January 2017 Events

On 18 January 2017, the 2016–2017 central Italy seismic sequence reached the Campotosto area with four events with magnitude larger than 5 in three hours (major event MW 5.5). To study the slip behavior on the causative fault/faults we followed two different methodologies: (1) we use Interferometric Synthetic Aperture Radar (InSAR) interferograms (Sentinel-1 satellites) and Global Positioning System (GPS) coseismic displacements to constrain the fault geometry and the cumulative slip distribution; (2) we invert near-source strong-motion, high-sampling-rate GPS waveforms, and high-rate GPS-derived static offsets to retrieve the rupture history of the two largest events. The geodetic inversion shows that the earthquake sequence occurred along the southern segment of the SW-dipping Mts. Laga normal fault system with an average slip of about 40 cm and an estimated cumulative geodetic moment of 9.29 × 1017 Nm (equivalent to a MW~6). This latter estimate is larger than the cumulative seismic moment of all the events, with MW &gt; 4 which occurred in the corresponding time interval, suggesting that a fraction (~35%) of the overall deformation imaged by InSAR and GPS may have been released aseismically. Geodetic and seismological data agree with the geological information pointing out the Campotosto fault segment as the causative structure of the main shocks. The position of the hypocenters supports the evidence of an up-dip and northwestward rupture directivity during the major shocks of the sequence for both static and kinematic inferred slip models. The activated two main slip patches are characterized by rise time and peak slip velocity in the ranges 0.7–1.1 s and 2.3–3.2 km/s, respectively, and by ~35–50 cm of slip mainly concentrated in the shallower northern part of causative fault. Our results show that shallow slip (depth &lt; 5 km) is required by the geodetic and seismological observations and that the inferred slip distribution is complementary with respect to the previous April 2009 seismic sequence affecting the southern half of the Campotosto fault. The recent moderate strain-release episodes (multiple M~5–5.5 earthquakes) and the paleoseismological evidence of surface-rupturing events (M~6.5) suggests therefore a heterogeneous behavior of the Campotosto fault.

Testing Different Tectonic Models for the Source of the Mw 6.5, 30 October 2016, Norcia Earthquake (Central Italy): A Youthful Normal Fault, or Negative Inversion of an Old Thrust?

AbstractWe 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.

Static Stress Increase in the Outer Forearc Produced by MW 8.2 September 8, 2017 Mexico Earthquake and its Relation to the Gravity Signal

An Mw 8.2 earthquake affected on September 8, 2017 the Cocos plate beneath the Caribbean plate next to the transform contact with the North American plate across the Tehuantepec gap, one of the two gaps that exist along the Mexican subduction zone offshore Chiapas. The epicenter occurred at a depth of 47 km and was associated with a highly steep normal fault parallel to the trench. The origin can be related to slab pull forces, dehydration processes and reactivation of outer rise faults. This earthquake produced positive changes in the Coulomb static stress, which would favor the development of thrust displacements along the interplate contact across the Tehuantepec gap instead of liberating tensions accumulated on it for decades. In the case that stresses were not released steadily in the future, a thrusting-related event with a magnitude of Mw 7.9 or higher could be expected after the September 8, 2017 Chiapas earthquake. However, the differential weight of the forearc offshore (that contains a region with a high positive gravity gradient) could inhibit rupture propagation to the coast precluding higher magnitudes. Combined and satellite only earth gravity field models show a highly heterogeneous density structure along the forearc of the North America-Caribbean plates and in particular at the studied zone. Maximum displacements in the rupture zone correlate to minimum gravity-derived anomalies showing a probable relationship between the differential weight of the forearc structure and the extensional faulting affecting the bending of the Cocos plate at depth. In this sense, variable sediment thicknesses over the forearc, and a high gravity gradient zone parallel to the trench in the outer forearc, could be reflecting physical heterogeneities that influence the propagation of the rupture zone. Additionally, to the north, a highly positive gradient signal related to the subducted Tehuantepec Fz projection beneath the North American plate marks the ending of the main slip patch acting as a barrier to the seismic rupture propagation. This event exemplifies (1) how extensional ruptures in the downgoing slab of a subduction zone can increment accumulated elastic strain across a seismic gap instead of liberating tensions, (2) that rupture propagation for this intraplate event was probably controlled (at shallow depths) by the density structure of the forearc, similarly to recent interplate earthquakes along the Chilean margin; (3) how a low gravity gradient signal and low density structures identified from the inversion model in the region of the outer forearc, where the normal stress increased after the main seismic event, suggest a higher risk of occurrence of a great megathrust earthquake.

Source Complexity of an Injection Induced Event: The 2016Mw5.1 Fairview, Oklahoma Earthquake

AbstractComplex rupture processes are occasionally resolved for weak earthquakes and can reveal a dominant direction of the rupture propagation and the presence and geometry of main slip patches. Finding and characterizing such properties could be important for understanding the nucleation and growth of induced earthquakes. One of the largest earthquakes linked to wastewater injection, the 2016 Mw5.1 Fairview, Oklahoma earthquake, is analyzed using empirical Green's function techniques to reveal its source complexity. Two subevents are clearly identified and located using a new approach based on relative hypocenter‐centroid location. The first subevent has a magnitude of Mw5.0 and shows the main rupture propagated toward the NE, in the direction of higher pore pressure perturbations due to wastewater injection. The second subevent appears as an early aftershock with lower magnitude, Mw4.7. It is located SW of the mainshock in a region of increased Coulomb stress, where most aftershocks relocated.

Complex Fault Geometry and Rupture Dynamics of the MW 6.5, 30 October 2016, Central Italy Earthquake

AbstractWe study the 30 October 2016 Norcia earthquake (MW 6.5) to retrieve the rupture history by jointly inverting seismograms and coseismic Global Positioning System displacements obtained by dense local networks. The adopted fault geometry consists of a main normal fault striking N155° and dipping 47° belonging to the Mt. Vettore‐Mt. Bove fault system (VBFS) and a secondary fault plane striking N210° and dipping 36° to the NW. The coseismic rupture initiated on the VBFS and propagated with similar rupture velocity on both fault planes. Updip from the nucleation point, two main slip patches have been imaged on these fault segments, both characterized by similar peak‐slip values (~3 m) and rupture times (~3 s). After the breakage of the two main slip patches, coseismic rupture further propagated southeastward along the VBFS, rupturing again the same fault portion that slipped during the 24 August earthquake. The retrieved coseismic slip distribution is consistent with the observed surface breakages and the deformation pattern inferred from interferometric synthetic aperture radar measurements. Our results show that three different fault systems were activated during the 30 October earthquake. The composite rupture model inferred in this study provides evidences that also a deep portion of the NNE trending section of the Olevano‐Antrodoco‐Sibillini thrust broke coseismically, implying the kinematic inversion of a thrust ramp. The obtained rupture history indicates that in this sector of the Apennines, compressional structures inherited from past tectonics can alternatively segment boundaries of NW trending active normal faults or break coseismically during moderate‐to‐large magnitude earthquakes.

Source Parameters of the 2016–2017 Central Italy Earthquake Sequence from the Sentinel-1, ALOS-2 and GPS Data

In this study, joint inversions of Synthetic Aperture Radar (SAR) and Global Position System (GPS) measurements are used to investigate the source parameters of four Mw &gt; 5 events of the 2016–2017 Central Italy earthquake sequence. The results show that the four events are all associated with a normal fault striking northwest–southeast and dipping southwest. The observations, in all cases, are consistent with slip on a rupture plane, with strike in the range of 157° to 164° and dip in the range of 39° to 44° that penetrates the uppermost crust to a depth of 0 to 8 km. The primary characteristics of these four events are that the 24 August 2016 Mw 6.2 Amatrice earthquake had pronounced heterogeneity of the slip distribution marked by two main slip patches, the 26 October 2016 Mw 6.1 Visso earthquake had a concentrated slip at 3–6 km, and the predominant slip of the 30 October 2016 Mw 6.6 Norcia earthquake occurred on the fault with a peak magnitude of 2.5 m at a depth of 0–6 km, suggesting that the rupture may have reached the surface, and the 18 January 2017 Mw 5.7 Campotosto earthquake had a large area of sliding at depth 3–9 km. The positive static stress changes on the fault planes of the latter three events demonstrate that the 24 August 2016 Amatrice earthquake may have triggered a cascading failure of earthquakes along the complex normal fault system in Central Italy.

2017 Valparaíso earthquake sequence and the megathrust patchwork of central Chile

AbstractIn April 2017, a sequence of earthquakes offshore Valparaíso, Chile, raised concerns of a potential megathrust earthquake in the near future. The largest event in the 2017 sequence was a M6.9 on 24 April, seemingly colocated with the last great‐sized earthquake in the region—a M8.0 in March 1985. The history of large earthquakes in this region shows significant variation in rupture size and extent, typically highlighted by a juxtaposition of large ruptures interspersed with smaller magnitude sequences. We show that the 2017 sequence ruptured an area between the two main slip patches during the 1985 earthquake, rerupturing a patch that had previously slipped during the October 1973 M6.5 earthquake sequence. A significant gap in historic ruptures exists directly to the south of the 2017 sequence, with large enough moment deficit to host a great‐sized earthquake in the near future, if it is locked.

Near-source high-rate GPS, strong motion and InSAR observations to image the 2015 Lefkada (Greece) Earthquake rupture history

The 2015/11/17 Lefkada (Greece) earthquake ruptured a segment of the Cephalonia Transform Fault (CTF) where probably the penultimate major event was in 1948. Using near-source strong motion and high sampling rate GPS data and Sentinel-1A SAR images on two tracks, we performed the inversion for the geometry, slip distribution and rupture history of the causative fault with a three-step self-consistent procedure, in which every step provided input parameters for the next one. Our preferred model results in a ~70° ESE-dipping and ~13° N-striking fault plane, with a strike-slip mechanism (rake ~169°) in agreement with the CTF tectonic regime. This model shows a bilateral propagation spanning ~9 s with the activation of three main slip patches, characterized by rise time and peak slip velocity in the ranges 2.5–3.5 s and 1.4–2.4 m/s, respectively, corresponding to 1.2–1.8 m of slip which is mainly concentrated in the shallower (<10 km) southern half of the causative fault. The inferred slip distribution and the resulting seismic moment (M0 = 1.05 × 1019 N m) suggest a magnitude of Mw 6.6. Our best solution suggests that the occurrence of large (Mw > 6) earthquakes to the northern and to the southern boundaries of the 2015 causative fault cannot be excluded.

Application and evaluation of a rapid response earthquake-triggered landslide model to the 25 April 2015 Mw 7.8 Gorkha earthquake, Nepal

The 25 April 2015 Mw 7.8 Gorkha earthquake produced strong ground motions across an approximately 250km by 100km swath in central Nepal. To assist disaster response activities, we modified an existing earthquake-triggered landslide model based on a Newmark sliding block analysis to estimate the extent and intensity of landsliding and landslide dam hazard. Landslide hazard maps were produced using Shuttle Radar Topography Mission (SRTM) digital topography, peak ground acceleration (PGA) information from the U.S. Geological Survey (USGS) ShakeMap program, and assumptions about the regional rock strength based on end-member values from previous studies. The instrumental record of seismicity in Nepal is poor, so PGA estimates were based on empirical Ground Motion Prediction Equations (GMPEs) constrained by teleseismic data and felt reports. We demonstrate a non-linear dependence of modeled landsliding on aggregate rock strength, where the number of landslides decreases exponentially with increasing rock strength. Model estimates are less sensitive to PGA at steep slopes (>60°) compared to moderate slopes (30–60°). We compare forward model results to an inventory of landslides triggered by the Gorkha earthquake. We show that moderate rock strength inputs over estimate landsliding in regions beyond the main slip patch, which may in part be related to poorly constrained PGA estimates for this event at far distances from the source area. Directly above the main slip patch, however, the moderate strength model accurately estimates the total number of landslides within the resolution of the model (landslides≥0.0162km2; observed n=2214, modeled n=2987), but the pattern of landsliding differs from observations. This discrepancy is likely due to the unaccounted for effects of variable material strength and local topographic amplification of strong ground motion, as well as other simplifying assumptions about source characteristics and their relationship to landsliding.