Dr. Benavente-Bravo, Roberto

Geodetically constrained interseismic interplate coupling has been widely used to assess seismic potential in subduction zones. Modeling interseismic deformation is challenging, as it involves interplate coupling and often ignores continental internal deformation processes. We present a novel methodology to jointly estimate interplate coupling along with upper plate rigid motion and surface strain, constrained by GNSS‐derived velocities. We use a least squares inversion with a spatially variable Equal Posterior Information Condition Tikhonov regularization, accounting for observational and elastic structure uncertainties. Our modeling reveals three megathrust regions with high tsunamigenic earthquake potential located within the Atacama Seismic Gap (Chile). This study indicates the presence of a downdip segmentation located just above the 1995 (Mw8.0) Antofagasta earthquake rupture, raising concerns for the potential of tsunamigenic earthquake occurrence at shallower depths. Additionally, we show that surface motion is dominated by strain, with rather negligible rigid motion, dismissing the rigid Andean microplate model typically assumed in previous studies.

Inversions of geodetic data are regularly used to estimate interseismic locking in subduction zones. However, the ill‐posed nature of these problems motivates us to include prior information, physically consistent with processes of the subduction seismic cycle. To deal with model instabilities, we present an inversion method to estimate both plate‐locking and model uncertainties by inverting Global Navigation Satellite System derived velocities based on a Bayesian model selection scheme. Our method allows us to impose positivity constraints via a multivariate folded‐normal distribution, with a specified covariance matrix. Model spatial correlations are explored and ranked to find models that best explain the observed data and for a better understanding of locking models. This approach searches for hyperparameters of the prior joint multivariate probability density function (PDF) of model parameters that minimize the Akaike Bayesian Information Criterion (ABIC). To validate our approach, we invert synthetic displacements from analytic models, yielding satisfactory results. We then apply the method to estimate the plate‐locking in Central Chile (28°–39°S) and its relation to the coseismic slip distribution of earthquakes with magnitudes Mw > 8.0, on the subduction zone since 2010. We also search among different prior PDFs for a single ductile‐fragile limit depth. Our results confirms a spatial correlation between locked asperities and the 2010 Mw 8.8 Maule and 2015 Mw 8.3 Illapel earthquake rupture zones. The robustness of our locking model shows potential to improve future seismic and tsunami hazard estimations.

The detection of transient events related to slow earthquakes in GNSS positional time series is key to understanding seismogenic processes in subduction zones. Here, we present a novel Principal and Independent Components Correlation Analysis (PICCA) method that allows for the temporal and spatial detection of transient signals. The PICCA is based on an optimal combination of the principal (PCA) and independent component analysis (ICA) of positional time series of a GNSS network. We assume that the transient signal is mostly contained in one of the principal or independent components. To detect the transient, we applied a method where correlations between sliding windows of each PCA/ICA component and each time series are calculated, obtaining the stations affected by the slow slip event and the onset time from the resulting correlation peaks. We first tested and calibrated the method using synthetic signals from slow earthquakes of different magnitudes and durations and modelled their effect in the network of GNSS stations in Chile. Then, we analyzed three transient events related to slow earthquakes recorded in Chile, in the areas of Iquique, Copiapó, and Valparaíso. For synthetic data, a 150 days event was detected using the PCA-based method, while a 3 days event was detected using the ICA-based method. For the real data, a long-term transient was detected by PCA, while a 16 days transient was detected by ICA. It is concluded that simultaneous use of both signal separation methods (PICCA) is more effective when searching for transient events. The PCA method is more useful for long-term events, while the ICA method is better suited to recognize events of short duration. PICCA is a promising tool to detect transients of different characteristics in GNSS time series, which will be used in a next stage to generate a catalog of SSEs in Chile

We compiled a novel microseismicity catalog for the Central Chile megathrust (29°–35°S), comprising 8,750 earthquakes between April 2014 and December 2018. These events describe a pattern of three trenchward open half-ellipses, consisting of a continuous, coast-parallel seismicity band at 30–45 km depth, and narrow elongated seismicity clusters that protrude to the shallow megathrust and separate largely aseismic regions along strike. To test whether these shapes could outline highly coupled regions (“asperities”) on the megathrust, we invert GPS displacement data for interplate locking. The best-fit locking model does not show good correspondence to seismicity, possibly due to lacking resolution. When we prescribe high locking inside the half-ellipses, however, we obtain models with similar data fits that are preferred according to the Bayesian Information Criterion (BIC). We thus propose that seismicity on the Central Chile megathrust may outline three adjacent highly coupled regions, two of them located between the rupture areas of the 2010 Maule and the 2015 Illapel earthquakes, a segment of the Chilean margin that may be in a late interseismic stage of the seismic cycle.