Reviewing megathrust slip behavior for recent Mw > 8.0 earthquakes along the Peru-Chilean margin from satellite GOCE gravity field derivatives

The characterization of a seismogenic zone associated with the rupture process that occurs during great megathrust earthquakes has been approached from different perspectives. Different studies lightened the structural complexity and compositional heterogeneity of the interplate region. In this work, we address this relationship from the analysis of density contrasts provided by satellite gravity data and comparing these results with slip behavior along the rupture zone. We performed the calculation of gravity anomalies and gradients both corrected by the effect of topography and sediments. Then we analyzed this relation with the rupture zones of four megathrust earthquakes Mw>8 along the Peru-Chilean coast, associated with the largest earthquakes that occurred in the last years along this setting, finding that the maximum vertical displacements were located close to gravity minimums in the forearc zone. We finally obtained density models by inverting four trench-parallel profiles from the Gravity disturbance, finding low-density contrasts related to the maximum slips, and high density contrasts at its edges related to a decrease of the displacement or to the lateral ending of the rupture. The along-strike seismic segmentation observed previously in the gravity derivatives is also observed along dip, being slip increased at gravity lows and arrested at gravity highs, probably indicating that the forearc structure plays a key role in across strike rupture behavior in depth. Our results agree with the hypothesis that persistent tectonic features, modeled at a high degree by the oceanic plate morphology, may control strain accumulation and release along the megathrust.