GRACE and GRACE Follow-On gravity observations of intermediate-depth earthquakes contrasted with those of shallow events

Earthquakes involve mass redistribution within the solid Earth and the ocean, and as a result, perturb the Earth's gravitational field. For most of the shallow (<60 km) earthquakes with M_w > 8.0, the GRACE satellite gravity measurements suggest considerable volumetric disturbance of rocks. At a spatial scale of hundreds of km, the effect of volumetric change exceeds gravity change by vertical deformation; for example, negative gravity anomalies associated with volumetric expansion are characteristic patterns after shallow thrust events. In this study, however, we report contrasting observations of gravity change from two intermediate-depth (100–150 km) earthquakes of 2016 & 2017 Mw 8.0 (two combined) Papua New Guinea thrust faulting events and 2019 Mw 8.0 Peru normal faulting and highlight the importance of compressibility in earthquake deformation. The combined 2016/17 thrust events resulted in a positive gravity anomaly of 5–6 microGal around the epicenter, while the 2019 normal faulting produced a negative gravity anomaly of 3–4 microGal. Our modeling found that these gravity changes are manifestation of vertical deformation with limited volumetric change, distinct from gravity changes after the shallow earthquakes. The stronger resistance of rocks to volume change at intermediate-depth results in largely incompressible deformation and thus in a gravity change dominated by vertical deformation. In addition, malleable rocks under high pressure and temperature at depth facilitated substantial afterslip and/or fast viscoelastic relaxation causing additional vertical deformation and gravity change equivalent to the coseismic change. For the Papua New Guinea events, this means that postseismic relaxation enhanced coseismic uplift and relative sea level decrease.