Postseismic Slip Research Articles

Co-Seismic and Post-Seismic Slip Properties Associated with the 2024 M 7.5 Noto Peninsula, Japan Earthquake Determined by GNSS Observations

Based on GNSS observations, the co-seismic and post-seismic slip of the 2024 Noto Peninsula earthquake and the spatio-temporal pattern of afterslip are investigated in this paper. The co-seismic slip is mainly distributed in the depth range of 2 to 15 km with the maximum value of 5.94 m. Compared with the co-seismic rupture pattern, a shallow afterslip can be observed after the earthquake, and the afterslip patch is formed northeast of the epicenter. The maximum value of afterslip during the post-seismic 180 days is 1.13 m, which is situated at the longitude of 137.53°, latitude of 37.75°, and epth of 5.43 km. The spatio-temporal evolution of afterslip indicates that the fault activity has continued throughout the post-seismic 180 days, and the coverage and magnitude of afterslip have gradually increased. As time goes on, the fault activity tends to weaken, as evidenced by a decrease in slip rate. The daily images of afterslip demonstrate that the fault activity is particularly strong in the early time period following the earthquake. The maximum value of afterslip in the first week accounts for about 18% of that in the post-seismic 180 days, and the maximum slip rate reaches 0.043 m/day. In addition, the Coulomb stress analysis indicates that afterslip and most aftershocks appear in the positive Coulomb stress region, suggesting that co-seismic Coulomb stress changes may be the driving mechanism of afterslip and aftershocks.

Afterslip and Creep in the Rate‐Dependent Framework: Joint Inversion of Borehole Strain and GNSS Displacements for the Mw 7.1 Ridgecrest Earthquake

AbstractThe elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub‐nanostrain precision and sub‐second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best‐fitting models adhere to mechanical rate‐and‐state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Surface ruptures of the 2022 Guanshan-Chihshang earthquakes in central Longitudinal Valley area, eastern Taiwan

AbstractThe Mw 6.4 and 6.8 Guanshan-Chihshang earthquakes occurred on 17 and 18 September 2022 resulted in prominent surface ruptures within the Longitudinal Valley in eastern Taiwan, particularly along the Yuli fault. Approximately 18 h after the mainshock, we began to document the surface rupture near Yuli Town. Our result suggests the surface rupture formed a confined single left-lateral trace in the town of Yuli, characterized by a series of en échelon right-stepping left-lateral faulting geometry. The rupture of 2022 roughly matches the locations of 1951 surface ruptures inside Yuli Town, with a similar amount of left-lateral cross-fault displacement. North and South of the Yuli residential area, we identified several sections of the surface rupture distributed in the water-saturated paddy fields. The maximum left-lateral displacement recorded across the rupture can reach 1.4 m just south of Yuli, with the fault scarp resembling a high-angle west-dipping fault geometry. In addition to the co-seismic surface ruptures, our repeating cross-fault measurements show significant post-seismic shallow after-slip along the Yuli fault. The amount of post-seismic deformation within 3 months after the mainshock is close to, or even higher than the co-seismic cross-fault displacement, consistent with local witness accounts and post-event field photos which showed continuous damage and displacement of building floors and roads after the earthquake. Such shallow post-seismic slips were also observed along the main fault trace in the 2014 South Napa earthquake, and likely represent the shallow elastoplastic behavior of the sub-vertical fault in the young alluvial sediments.

Early postseismic slip of the 21 November 2022 Mw 5.6 Cianjur, Indonesia, earthquake based on GPS measurements

ABSTRACT Global Positioning System (GPS) data after the 21 November 2022 Mw 5.6 Cianjur earthquake are crucial for understanding the physical processes during the postseismic phase. New continuous GPS stations were installed and recorded data for periods ranging from five to twenty-five days after the mainshock. The daily time series data at GPS stations shows an evolving logarithmic decay over time, as expected due to frictional afterslip. We inverted the postseismic displacement data to solve the postseismic slip distribution, utilising Akaike’s Bayesian information criterion. A cumulative aseismic slip moment of 1.54 × 1016 N·m, equivalent to Mw 4.7, is modelled during the period from five to twenty-five days after the mainshock. This aseismic slip moment is equivalent to approximately 7% of the seismic moment during the mainshock. Our analysis on the newly identified active fault using GPS displacements highlights the imperative need to identify active faults prior to future earthquake occurrences, especially when conducting hazard mitigation in populated regions in Indonesia.

Surface Rupture and Fault Characteristics Associated With the 2020 Magnitude (MW) 6.6 Masbate Earthquake, Masbate Island, Philippines

AbstractOn 18 August 2020, Masbate Island was struck by a magnitude (MW) 6.6 earthquake. This seismic event represents the second occurrence of a strong earthquake (M > 6) in 17 years, emphasizing the necessity for further investigation into the characteristics of this event. In this study, we employ Interferometric Synthetic Aperture Radar, seismicity analysis, and field investigations to comprehensively characterize the coseismic and postseismic slip associated with the event. Our findings reveal a 50‐km‐long fault rupture along the Masbate segment of the Philippine Fault, with ∼23 km surface rupture mapped onshore, despite the occurrence of interseismic creep. The slip distribution demonstrates decreasing displacements northwestward toward the creeping section, with a maximum left‐lateral displacement of 0.97 m near the epicenter. Toward the southeast offshore, the rupture terminates at a left stepover of a fault. While the surface rupture appears relatively straight and narrowly concentrated, the secondary ruptures and mapped offshore faults reveal a more complex transtensional fault structure in the southeastern part of Masbate Island. This fault complexity represents an asperity that facilitates high‐stress accumulation and rupture initiation. Postseismic slip persists for several months along the onshore creeping segment. Based on comprehensive measurements of both cumulative and coseismic slip along the Masbate fault segment, we calculate a slip rate ranging between 2.8 and 3.8 cm/year and a recurrence interval of 16–41 years for earthquakes similar to the 2020 earthquake. Our study highlights how heterogeneity in fault properties, including geometry and coupling state, influences the distribution of slip and magnitude of earthquakes. The 2020 Masbate earthquake provides valuable insights into the rupture dynamics and fault behavior of the Philippine Fault in the Masbate region.

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.

Coseismic and Early Postseismic Deformation of the 2024 Mw7.45 Noto Peninsula Earthquake

AbstractAn unexpected Mw7.45 earthquake struck the Noto Peninsula on 1 January 2024, preceded by several long‐living earthquake swarms, providing a valuable opportunity to study seismic and aseismic slips, as well as their interactions. We derived coseismic and 19‐day postseismic slip distributions by inverting co‐ and post‐seismic displacements from Global Navigation Satellite System (GNSS) data. The inverted coseismic slip distribution shows two slip patches, with a maximum slip of ∼4 m. The early postseismic afterslip is 0.1–0.25 m within coseismic slip asperity and 0.1–0.6 m northward of the rupture area. The afterslip within the rupture area is accompanied by numerous aftershocks and coincides with a ∼6 MPa stress drop, suggesting that aftershocks are likely driven by the afterslip. The pattern of poroelastic rebound implies a potential effect of fluid flow on aftershock triggering. This study sheds lights on the intricate interplay between seismic and aseismic processes following large earthquakes.

Assessing Slip Rates on the Xianshuihe Fault Using InSAR with Emphasis on Phase Unwrapping Error and Atmospheric Delay Corrections

Located on the southeastern periphery of the Tibetan Plateau, the Xianshuihe fault (XSHF) is an active left-lateral strike-slip fault renowned for its frequent and intensive seismic activities. This highlights the necessity of employing advanced geodetic methodologies to precisely evaluate the fault kinematics and seismic hazard potential along this fault. Among these techniques, interferometric synthetic aperture radar (InSAR) stands out for its high spatial resolution and regular revisit intervals, enabling accurate mapping of interseismic deformation associated with fault motion. However, the precision of InSAR in measuring deformation encounters several challenges, particularly artifacts stemming from phase unwrapping errors and atmospheric phase delays. In this study, we utilize ascending and descending Sentinel-1 InSAR images spanning from January 2017 to January 2023 to drive the line-of-sight (LOS) mean crustal velocities associated with the XSHF with emphasis on phase unwrapping errors and atmospheric delay corrections. Then, the reliability of the derived LOS velocities is assessed using independent observations from the Global Navigation Satellite System (GNSS). The inferred fault slip rate along the XSHF shows significant along-strike variations, gradually decreasing from ~11.1 mm/yr at the Luhuo section to ~6.6 mm/yr at the Kangding section and then sharply increasing to ~13.0 mm/yr towards its eastern terminus at the Moxi section. The fault locking depth shows similar along-strike variations, decreasing from ~19.5 km in the northwestern part to ~4.8 km at the Kangding section, before increasing to 19.6 km at the Moxi segment. Notably, apparent surface fault creeping, characterized by a slip rate of ~2.7 mm/yr, is observed at the Kangding segment, likely resulting from postseismic slip following the 2014 Mw 6.3 Kangding earthquake.

The relationship between kinematics and fault geometry for surface coseismic ruptures on across-strike faults: New observations of slip vectors and displacements along the Pisia and Skinos faults from the 1981 Eastern Gulf of Corinth, Greece earthquakes

The relationships between kinematics and fault geometry for the coseismic ruptures from the 24th and February 25, 1981 earthquake sequence in the eastern Gulf of Corinth (Ms 6.7 and 6.4) are analysed. The two earthquakes ruptured faults located across strike rather than along strike as typifies other earthquake sequences. In detail, surface ruptures formed on the sub-parallel Pisia and Skinos Faults, with an 8 km along-strike overlap zone, separated across strike by < 2 km. The largest coseismic offsets occurred in the overlap zone. The 41-year-old ruptures are still well preserved as bedrock fault plane lichen-free stripes and colluvial ruptures, allowing detailed structural mapping at 213 rupture localities. A comparison between our measurements and Jackson et al. (1982) showed no overall consistent signal of post-seismic slip as some of our measurements were greater and some smaller than those recorded in 1981. The ruptures produced a single maximum asymmetric profile (Pisia: maximum throw of 223 cm) and a double maxima profile (Skinos: maximum throw of 109 cm and 130 cm). The shapes of the profiles differed in previous earthquakes on these faults, as evidenced by an older lichen-free stripe, implying non-characteristic earthquakes. Summing the two overlapping throw profiles across-strike reveals a single maximum symmetric bell-like profile. Using the above observations on coseismic offsets, kinematic information, and the geometry of faults, a rupture scenario has been proposed in terms of fault bends and corrugation orientations which suggests that parts of each fault may have ruptured in each earthquake.

GNSS imaging coseismic and postseismic slip associated with the 2021 M 8.2 Chignik, Alaska earthquake

The coseismic and postseismic slip associated with the 2021 M 8.2 Chignik, Alaska earthquake is imaged using the geodetic modeling method. Further, the spatio-temporal evolution of afterslip and aftershocks (0–60 days) is also investigated based on the postseismic deformations captured by Global Navigation Satellite System (GNSS) sites in this study. The coseismic slip is mainly concentrated at depths of 20–50 km with a maximum of 3.92 m. After this earthquake, a significant afterslip can be observed in the shallow region, which is distributed around the coseismic rupture region. Afterslip is primarily accumulated at depths of 10–30 km, and the maximum reaches 1.29 m, which is located at a longitude of 157.19°W, latitude of 55.09°N, and depth of 21.80 km. The spatio-temporal properties of afterslip indicate that the spatial extent and magnitude of afterslip continuously increase during the 60 days after the earthquake, and the fault activity in the first seven days is significantly strong. Meanwhile, the fault slip rate during the time span of 0–7 days is the largest with a maximum of 0.054 m/day. The fault slip rate gradually decreases with time, implying that the fault activity is slowly weakening. The Coulomb Failure Stress (CFS) distribution suggests that afterslip and most aftershocks are located at the southeastern region of the hypocenter, which is covered by positive CFS. Moreover, in addition to coseismic rupture effects, afterslip may be a possible triggering mechanism of early aftershocks as evidenced by the fact that the spatial distributions of afterslip and aftershocks are highly consistent.

Geomechanical modelling of injection-induced seismicity: the case study of the Muara Laboh geothermal plant

SUMMARY In this work, we study the induced seismicity recorded during an injection operation at the Muara Laboh geothermal plant (Indonesia). The swarm, consisting of three bursts activating a normal fault zone, is characterized by rapid earthquake (km d−1) migration. We use a 2-D rate-and-state asperity model to better understand the physical mechanisms controlling the evolution of this induced swarm. The model suggests that the observed rapid seismic migration can be explained by the interaction among asperities through the expansion of slow post-seismic slip fronts. Also, it shows that the amount of seismicity generated by the fluid injection is strongly controlled by the background seismicity of the system, that is by the seismicity determined by the tectonic load charging the fault. This close correlation between natural and induced seismicity suggests that the injection in Muara Laboh principally stimulates critically stressed faults, which release the seismicity determined by their natural seismic cycle.

Asperity Size and Neighboring Segments Can Change the Frictional Response and Fault Slip Behavior: Insights From Laboratory Experiments and Numerical Simulations

AbstractAccurate assessment of the rate and state friction parameters of rocks is essential for producing realistic earthquake rupture scenarios and, in turn, for seismic hazard analysis. Those parameters can be directly measured on samples, or indirectly based on inversion of coseismic or postseismic slip evolution. However, both direct and indirect approaches require assumptions that might bias the results. Aiming to reduce the potential sources of bias, we take advantage of a downscaled analog model reproducing megathrust earthquakes. We couple the simulated annealing algorithm with quasi‐dynamic numerical models to retrieve rate and state parameters reproducing the recurrence time, rupture duration and slip of the analog model, in the ensemble. Then, we focus on how the asperity size and the neighboring segments' properties control the seismic cycle characteristics and the corresponding variability of rate and state parameters. We identify a tradeoff between (a–b) of the asperity and (a–b) of neighboring creeping segments, with multiple parameter combinations that allow mimicking the analog model behavior. Tuning of rate and state parameters is required to fit laboratory experiments with different asperity lengths. Poorly constrained frictional properties of neighboring segments are responsible for uncertainties of (a–b) of the asperity in the order of per mille. Roughly one order of magnitude larger uncertainties derive from asperity size. Those results provide a glimpse of the variability that rate and state friction estimates might have when used as a constraint to model fault slip behavior in nature.

GEODETIC EVIDENCE FOR POST-SEISMIC DEFORMATIONS FOLLOWING THE 2014 NORTH AEGEAN MW 6.9 EARTHQUAKE

The 2014 North Aegean Sea earthquake was a strong (Mw 6.9) event that caused significant crustal deformations. In the present study we investigate the long-term impact of the earthquake on the kinematics of the North Aegean Trough (NAT). For this purpose, we analyzed GPS observations collected from May 2010 to April 2022 at five permanent GPS reference stations. Two of these stations are located close to the epicenter(s) on the Islands of Lemnos and Samothrace. We processed the data using the Precise Point Positioning (PPP) technique. The analysis of the obtained coordinate time-series revealed a post-seismic deformation (PSD) period lasting for more than two years leading to cumulative 2D post-seismic displacement of 22 mm and 27 mm for Samothrace and Lemnos, respectively. The magnitudes of these post-seismic slip vectors correspond to 23% and 49% of the co-seismic vectors at Samothrace and Lemnos, respectively. The long-term analysis showed that after the end of the PSD period the stations are characterized by stable velocities that are noticeably different compared to the velocities prior to the event. We observed a change in the velocity in the order of 2 mm/yr for both Samothrace and Lemnos. It is the first time that PSD and velocity changes have been reported for the 2014 North Aegean Sea earthquake shedding light on the characteristics and the impact of this important earthquake on the kinematics of NAT.

Transient postseismic slip and aftershock triggering: A case study of the 2008 MW7.9 Wenchuan Earthquake, China

In this study, we investigate how a stress variation generated by a fault that experiences transient postseismic slip (TPS) affects the rate of aftershocks. First, we show that the postseismic slip from Rubin-Ampuero model is a TPS that can occur on the main fault with a velocity-weakening frictional motion, that the resultant slip function is similar to the generalized Jeffreys-Lomnitz creep law, and that the TPS can be explained by a continuous creep process undergoing reloading. Second, we obtain an approximate solution based on the Helmstetter-Shaw seismicity model relating the rate of aftershocks to such TPS. For the Wenchuan sequence, we perform a numerical fitting of the cumulative number of aftershocks using the Modified Omori Law (MOL), the Dieterich model, and the specific TPS model. The fitting curves indicate that the data can be better explained by the TPS model with a B/A ratio of approximately 1.12, where A and B are the parameters in the rate- and state-dependent friction law respectively. Moreover, the p and c that appear in the MOL can be interpreted by the B/A and the critical slip distance, respectively. Because the B/A ratio in the current model is always larger than 1, the model could become a possible candidate to explain aftershock rate commonly decay as a power law with a p-value larger than 1. Finally, the influence of the background seismicity rate r on parameters is studied; the results show that except for the apparent aftershock duration, other parameters are insensitive to r.

Insights into Very Early Afterslip Associated with the 2021 M 8.2 Chignik, Alaska Earthquake Using Subdaily GNSS Solutions

Based on subdaily kinematic GNSS solutions, the fault slip properties during the very early postseismic phase after the 2021 M 8.2 Chignik earthquake are investigated in this paper. The very early postseismic deformations captured by near-field GNSS sites can be well depicted by the power model. The comparison of afterslip determined by daily and subdaily GNSS solutions suggests that neglecting very early afterslip can result in the underestimation of postseismic slip. Compared with coseismic slip, the cumulative afterslip of the first 24 h is mainly focused in the southeast of the hypocenter, and the shallow updip afterslip appears after this earthquake. The spatio-temporal evolution of the afterslip reveals that the patch of afterslip is immediately generated after the earthquake, and then the postseismic slip gradually grows along the afterslip patch. The magnitude of the afterslip patch varies remarkably within the 24 h following the earthquake, especially in the first several hours. Meanwhile, the spatio-temporal patterns of aftershocks and afterslip exhibit strong similarity during the first 24 h, suggesting that very early afterslip may be a possible driving factor of aftershocks. Moreover, most of the afterslip patches and aftershocks occurring immediately after this earthquake are situated in the area covered by positive Coulomb Stress Change (CSC), which implies that the immediate afterslip and aftershock activities can be influenced by the coseismic CSC. The following afterslip process further releases coseismic CSC and then influences the spatio-temporal variations of aftershock activities. Thus, the afterslip may be a possible triggering mechanism of very early aftershocks for this earthquake, alongside the effects of the CSC generated by coseismic rupture.

Coseismic and Early Postseismic Slip of the 2021 Mw 7.2 Nippes, Haiti, Earthquake: Transpressional Rupture of a Nonplanar Dipping Fault System

Abstract On 14 August 2021, an Mw 7.2 earthquake struck Nippes, Haiti, 11 yr after the devastating 2010 Mw 7.0 Port-au-Prince earthquake. This earthquake occurred in a remote region where the structure at the depth of the main boundary Enriquillo Plantain Garden fault (EPGF) is less known. Using Synthetic Aperture Radar imagery, we retrieve the coseismic and early postseismic deformation of the 2021 Haiti earthquake to constrain its fault geometry and slip distribution. Our modeling results show that the 2021 earthquake ruptured the high-angle Ravine du Sud fault and a bend fault ∼64° dipping to the north at depth. Although not only conclusive, the combination of coseismic and postseismic deformation, along with geomorphic features, and relocated aftershocks, suggest a nonplanar fault structure with significant variations in dip angles along both the depth and track of the EPGF in this region. East of the epicenter, we document a 25 km section along the EPGF that crept for ∼15 days. This distribution of aseismic slip utilizing stacked deformation indicates that only a small fraction of the accumulated strain near the surface was released during the earthquake, suggesting a high potential for seismic hazard in the region along the EPGF from the ruptured segment to the east, before reaching the 2010 rupture.

Deep Postseismic Creep Following Large Earthquakes Revealed by Repeating Aftershocks in the Southeastern Tibetan Plateau

Abstract Large earthquake occurrence and the subsequent postseismic period are the most dramatic part of a seismic cycle that usually lasts months to years. However, the fault dynamics that account for the postseismic events are yet to be fully understood. Here, we use the repeating aftershock sequences (RASs) to investigate postseismic slips following the Mw 6.6 Lushan, Mw 6.5 Jiuzhaigou, Mw 6.1 Jinggu, and Mw 6.2 Ludian earthquakes in the southeastern Tibetan Plateau and find 135 RASs following the mainshocks. The RAS seismicity suggests that seismogenic faults began to creep in depth within a few hours after the Lushan, Jiuzhaigou, and Jinggu mainshocks. The deep creeps mainly follow a velocity-strengthening friction mode and decay with an Omori law p-value of ∼1. The results suggest that the combination of fault healing and geometry together controls deep fault behaviors. We develop two conceptual models to explain our observations. Our results provide new insights into spatiotemporal fault evolution after large earthquakes.

Extraction of transient signal from GPS position time series by employing ICA

Transient deformation, such as post-seismic slip, slow slip and pre-seismic slip events, is a limited low-frequency deformation that can last for hours to months, in contrast to a sudden slip on a fault caused by earthquakes. Continuous Global Positioning System (CGPS), one of the most common geodetic techniques for continuously monitoring crustal deformation, is capable of capturing transient deformation signals. A critical point in characterizing transient deformation signals is the development of extracting and deciphering transient deformation signals from a huge and messy data set of position time series. Principal Component Analysis (PCA), one of the data-driven methods, has been employed to derive transient deformation signals from position time series combing with Kalman filtering. Independent Component Analysis (ICA) performs well in recovering and separating the sources of observed data, however, it is rarely used in extracting transient deformation signals. We aim to decompose the transient deformation signals from the daily GPS observation deployed in Akutan Island from 2007 to 2015 with the ICA method and obtain the spatiotemporal responses to the source signals of transient deformation. Our results indicate that ICA method can also characterize effectively transient deformation signals spatially and temporally. Additionally, the independent relationship between sources obtained by ICA allows for flexibility in linearly combining different sources.

Seismic and aseismic slip during the 2006 Copiapó swarm in North-Central Chile

Earthquake swarms commonly occur along the Chilean subduction zone, witnessing fast seismic and slow aseismic slip behavior at the plate interface. However, the largest seismic swarms observed in Chile, particularly in the Copiapó-Atacama region, remain poorly documented, and the underlying processes have yet to be understood. Here, we perform seismological and geodetic analyses to investigate the 2006 Copiapó swarm, which developed in April and May 2006. The swarm began on April 19, with a magnitude Ml 5.3 earthquake. During the nine following days, we observe a migration of seismicity along the plate interface, the occurrence of doublets events, and a potential slow slip event in the GPS time series at site Copiapó. Then, on April 30, a first earthquake with Mw 6.6 occurred at 15 km depth at the plate contact. It likely triggered a second earthquake of magnitude Mw 6.5, which occurred 144 min later, 10 km northwest of the first earthquake. Using InSAR, we determined the slip distribution associated with these two earthquakes and detailed the postseismic slip they triggered in the next days and weeks. This “postseismic” phase appears to be predominantly aseismic, while the moment released during the “coseismic” phase is comparable to other seismic crises that occurred in Atacama. Although we did not find a larger seismic and aseismic ratio than in other swarms in South America, we suggest a similar mechanism of slow deformation as a driver of seismicity during seismic swarms. Finally, we propose that the slow and fast behavior of the 2006 Copiapó swarm is a consequence of the subduction of the Copiapó Ridge seamounts, which affects both the plate interface and the overriding plate by inducing complex interactions between seismic and aseismic processes.

Finite-Source Model of the 8 July 2021 M 6.0 Antelope Valley, California, Earthquake

Abstract We present a finite-source coseismic slip model of the 2021 Mw 6.0 Antelope Valley earthquake based on the joint inversion of regional seismic waveforms, Interferometric Synthetic Aperture Radar (InSAR), and Global Navigation Satellite Systems data. The results show that the mainshock rupture was dominated by normal slip along a nearly north–south-trending fault dipping to the east. The rupture lasted for ∼10 s, with primarily unilateral propagation toward the south. Most coseismic slip is found to be within a depth range between 6 and 10 km, with apparently no slip reaching the surface. Surface projection of the modeled fault plane matches well with the southern extension of the previously mapped Slinkard Valley fault (SVF). Aftershocks one year after the mainshock are mostly distributed within a relatively narrow band of 2–3 km thickness around the up-dip portion of the inferred coseismic rupture plane. There is little aftershock activity below 10 km, suggesting a relatively shallow brittle-to-ductile transition in this area. Aftershocks are also clustered at shallow depth beneath several branches of the Antelope Valley faults to the east of the mainshock rupture, including the Mw 4.4 event on 27 August 2021, which produced clear coseismic surface deformation observed by InSAR. Most aftershocks, immediately up-dip of the coseismic rupture and to the east beneath the Antelope Valley faults, are in areas of substantial coseismic Coulomb stress increase, particularly when assuming that all faults in this area dip to the east. This suggests that like the SVF that hosted the mainshock, the Antelope Valley faults in this area also dip to the east. There is little to no postseismic deformation seen from InSAR observations ∼2 months after the mainshock. The lack of clear coseismic and postseismic slip on the shallow portion of the fault suggests the potential for future shallow seismic activity.