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.