Microseismicity Appears to Outline Highly Coupled Regions on the Central Chile Megathrust

Fast relocking and afterslip-seismicity evolution following the 2015 Mw 8.3 Illapel earthquake in Chile

Large subduction earthquakes induce complex postseismic deformation, primarily driven by afterslip and viscoelastic relaxation, in addition to interplate relocking processes. However, these signals are intricately intertwined, posing challenges in determining the timing and nature of relocking. Here, we use six years of continuous GNSS measurements (2015-2021) to study the spatiotemporal evolution of afterslip, seismicity and locking after the 2015 Illapel earthquake ([Formula: see text] 8.3). Afterslip is inverted from postseismic displacements corrected for nonlinear viscoelastic relaxation modeled using a power-law rheology, and the distribution of locking is obtained from the linear trend of GNSS stations. Our results show that afterslip is mainly concentrated in two zones surrounding the region of largest coseismic slip. The accumulated afterslip (corresponding to [Formula: see text] 7.8) exceeds 1.5 m, with aftershocks mainly occurring at the boundaries of the afterslip patches. Our results reveal that the region experiencing the largest coseismic slip undergoes rapid relocking, exhibiting the behavior of a persistent velocity weakening asperity, with no observed aftershocks or afterslip within this region during the observed period. The rapid relocking of this asperity may explain the almost regular recurrence time of earthquakes in this region, as similar events occurred in 1880 and 1943.

B-value variations in the Central Chile seismic gap assessed by a Bayesian transdimensional approach

The b-value can be used to characterize the seismic activity for a given earthquake catalog and provide information on the stress level accumulated at active faults. Here we develop an algorithm to objectively estimate variations of b-value along one arbitrary dimension. To this end, we employ a Bayesian transdimensional approach where the seismic domains will be self-defined according to information in the seismic catalog. This makes it unnecessary to prescribe the location and extent of domains, as it is commonly done. We first show the algorithm's robustness by performing regressions from synthetic catalogs, recovering the target models with great accuracy. We also apply the algorithm to a microseismicity catalog for the Central Chile region. This segment is considered a seismic gap where the last major earthquake with shallow slip was in 1730. Our results illuminate the downdip limit of the seismogenic zone and the transition to intraslab seismicity. In the along-strike direction, low b-value coincides with the extent of locked asperities, suggesting a high-stress loading at the Central Chile seismic gap. Our results indicate the reliability of the Bayesian transdimensional method for capturing robust b-value variations, allowing us to characterize the mechanical behavior on the plate interface of subduction zones.

Measuring Coastal Subsidence after Recent Earthquakes in Chile Central Using SAR Interferometry and GNSS Data

Coastal areas concentrate a large portion of the country’s population around urban areas, which in subduction zones commonly are affected by drastic tectonic processes, such as the damage earthquakes have registered in recent decades. The seismic cycle of large earthquakes primarily controls changes in the coastal surface level in these zones. Therefore, quantifying temporal and spatial variations in land level after recent earthquakes is essential to understand shoreline variations better, and to assess their impacts on coastal urban areas. Here, we measure the coastal subsidence in central Chile using a multi-temporal differential interferometric synthetic aperture radar (MT-InSAR). This geographic zone corresponds to the northern limit of the 2010 Maule earthquake (Mw 8.8) rupture, an area affected by an aftershock of magnitude Mw 6.8 in 2019. The study is based on the exploitation of big data from SAR images of Sentinel-1 for comparison with data from continuous GNSS stations. We analyzed a coastline of ~300 km by SAR interferometry that provided high-resolution ground motion rates from between 2018 and 2021. Our results showed a wide range of subsidence rates at different scales, of analyses on a regional scale, and identified the area of subsidence on an urban scale. We identified an anomalous zone of subsidence of ~50 km, with a displacement <−20 mm/year. We discuss these results in the context of the impact of recent earthquakes and analyze the consequences of coastal subsidence. Our results allow us to identify stability in urban areas and quantify the vertical movement of the coast along the entire seismic cycle, in addition to the vertical movement of coast lands. Our results have implications for the planning of coastal infrastructure along subduction coasts in Chile.