Modeling the Contribution of Poroelastic Deformation to Postseismic Geodetic Signals

To constrain the poroelastic component of postseismic deformation, we model the subsurface hydrologic response to the Mw 7.6 subduction zone earthquake that occurred on the plate interface beneath the Nicoya peninsula in Costa Rica on 5 September 2012. The model shows that poroelastic relaxation occurs on multiple time scales and the associated deformation can be up to 2 cm for the trench‐perpendicular component. By modeling the time‐dependent deformation associated with poroelastic relaxation, we can begin to remove its contribution from the observed geodetic signal. Inversions for after slip that ignore poroelastic deformation have errors of 10–20% overall and up to 50% locally. Poroelastic effects can both mute and amplify the inferred afterslip.

[1]  Y. Fukuda,et al.  Probing the Poisson's ratio of poroelastic rebound following the 2011 Mw 9.0 Tohoku earthquake , 2018, Geophysical Journal International.

[2]  A. Gualandi,et al.  Poroelasticity and Fluid Flow Modeling for the 2012 Emilia-Romagna Earthquakes: Hints from GPS and InSAR Data , 2018, Geofluids.

[3]  M. Hesse,et al.  Modeling the poroelastic response to megathrust earthquakes: A look at the 2012 Mw 7.6 Costa Rican event , 2017 .

[4]  T. Dixon,et al.  Slow slip events in the early part of the earthquake cycle , 2017 .

[5]  C. Kyriakopoulos,et al.  Large and primarily updip afterslip following the 2012 Mw 7.6 Nicoya, Costa Rica, earthquake , 2017 .

[6]  C. Kyriakopoulos,et al.  Large and primarily updip afterslip following the 2012 M w 7.6 Nicoya, Costa Rica, earthquake: Afterslip From 2012 Nicoya Earthquake , 2017 .

[7]  C. Bignami,et al.  Aftershocks, groundwater changes and postseismic ground displacements related to pore pressure gradients: Insights from the 2012 Emilia‐Romagna earthquake , 2017 .

[8]  T. Dixon,et al.  Strain release at the trench during shallow slow slip: The example of Nicoya Peninsula, Costa Rica , 2017 .

[9]  Hiroshi P. Sato,et al.  Resolving depth-dependent subduction zone viscosity and afterslip from postseismic displacements following the 2011 Tohoku-oki, Japan earthquake , 2017 .

[10]  Andrew T. T. McRae,et al.  Firedrake , 2016, Mathematics of Planet Earth.

[11]  David A. Ham,et al.  An Algorithm for the Optimization of Finite Element Integration Loops , 2016, ACM Trans. Math. Softw..

[12]  S. Schwartz,et al.  Monitoring transient changes within overpressured regions of subduction zones using ambient seismic noise , 2016, Science Advances.

[13]  T. Dixon,et al.  Multiscale postseismic behavior on a megathrust: The 2012 Nicoya earthquake, Costa Rica , 2015 .

[14]  H. Villinger,et al.  Slow and delayed deformation and uplift of the outermost subduction prism following ETS and seismogenic slip events beneath Nicoya Peninsula, Costa Rica , 2015 .

[15]  T. Dixon,et al.  Earthquake and tsunami forecasts: Relation of slow slip events to subsequent earthquake rupture , 2014, Proceedings of the National Academy of Sciences.

[16]  X. Kuang,et al.  An integrated permeability‐depth model for Earth's crust , 2014 .

[17]  Yan Hu,et al.  Contributions of poroelastic rebound and a weak volcanic arc to the postseismic deformation of the 2011 Tohoku earthquake , 2014, Earth, Planets and Space.

[18]  Agustan,et al.  A comprehensive model of postseismic deformation of the 2004 Sumatra–Andaman earthquake deduced from GPS observations in northern Sumatra , 2014 .

[19]  S. Owen,et al.  Nicoya earthquake rupture anticipated by geodetic measurement of the locked plate interface , 2014 .

[20]  S. Owen,et al.  The 5 September 2012 Nicoya, Costa Rica Mw 7.6 earthquake rupture process from joint inversion of high‐rate GPS, strong‐motion, and teleseismic P wave data and its relationship to adjacent plate boundary interface properties , 2013 .

[21]  Anders Logg,et al.  Unified form language: A domain-specific language for weak formulations of partial differential equations , 2012, TOMS.

[22]  Eric J. Fielding,et al.  Fault slip models of the 2010–2011 Canterbury, New Zealand, earthquakes from geodetic data and observations of postseismic ground deformation , 2012 .

[23]  Michael Manga,et al.  Hydrological effects of the M W 7.1 Darfield (Canterbury) earthquake, 4 September 2010, New Zealand , 2012 .

[24]  Timothy H. Dixon,et al.  Active deformation near the Nicoya Peninsula, northwestern Costa Rica, between 1996 and 2010: Interseismic megathrust coupling , 2012 .

[25]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[26]  Sylvain Barbot,et al.  A unified continuum representation of post-seismic relaxation mechanisms: semi-analytic models of afterslip, poroelastic rebound and viscoelastic flow , 2010 .

[27]  Walter D. Mooney,et al.  Poroelastic stress-triggering of the 2005 M8.7 Nias earthquake by the 2004 M9.2 Sumatra–Andaman earthquake , 2010 .

[28]  Anders Logg,et al.  DOLFIN: Automated finite element computing , 2010, TOMS.

[29]  Paul Segall,et al.  Earthquake and Volcano Deformation , 2010 .

[30]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[31]  Y. Fialko,et al.  Hydrologic detection and finite element modeling of a slow slip event in the Costa Rica prism toe , 2009 .

[32]  B. E. Shaw,et al.  Afterslip and aftershocks in the rate‐and‐state friction law , 2007, physics/0703249.

[33]  S. Ge,et al.  Hydrodynamic response of subduction zones to seismic activity: A case study for the Costa Rica margin , 2006 .

[34]  Eric J. Fielding,et al.  Coseismic and Postseismic Slip of the 2004 Parkfield Earthquake from Space-Geodetic Data , 2006 .

[35]  Yuri Fialko,et al.  Evidence of fluid-filled upper crust from observations of postseismic deformation due to the 1992 Mw7.3 Landers earthquake , 2004 .

[36]  Hugo Perfettini,et al.  Postseismic relaxation driven by brittle creep: A possible mechanism to reconcile geodetic measurements and the decay rate of aftershocks, application to the Chi-Chi earthquake, Taiwan , 2004 .

[37]  Paul Lundgren,et al.  Geodetic and Seismic Constraints on some Seismogenic Zone Processes in Costa Rica , 2003 .

[38]  Paul Segall,et al.  Post-earthquake ground movements correlated to pore-pressure transients , 2003, Nature.

[39]  Herbert F. Wang.  Theory of Linear Poroelasticity with Applications to Geomechanics and Hydrogeology , 2000 .

[40]  Kenneth W. Hudnut,et al.  Poroelastic rebound along the Landers 1992 earthquake surface rupture , 1998 .

[41]  T. Dixon,et al.  Seismogenic zone structure beneath the Nicoya Peninsula, Costa Rica, from three‐dimensional local earthquake P‐ and S‐wave tomography , 2006 .

[42]  Large and primarily updip afterslip following the 2012 M w 7.6 Nicoya, Costa Rica, earthquake Geophysical Research: Solid Earth , 2022 .