Time-dependent dike propagation from joint inversion of seismicity and deformation data

• Published in: Journal of geophysical research. Solid earth 2013-11, Vol.118 (11), p.5785-5804
• Main Authors: Segall, P., Llenos, A. L., Yun, S.-H., Bradley, A. M., Syracuse, E. M.
• Format: Article
• Language: English
• Subjects: Decay; Deformation; dike; Dikes; Earth sciences; Earth, ocean, space; Earthquakes; Exact sciences and technology; Geophysics; Global positioning systems; GPS; Intrusion; inversion; Inversions; Mathematical models; Propagation; Seismic activity; Seismic phenomena; Seismicity; Seismology; Volcanoes
• ISSN: 2169-9313; 2169-9356
• DOI: 10.1002/2013JB010251

Dike intrusions both deform Earth's surface and induce propagating earthquake swarms. We develop methods to utilize both deformation and seismicity from brittle, volcano‐tectonic earthquakes to image time‐dependent dike propagation. Dieterich's (1994) seismicity‐rate theory is used to relate dike‐induced stress changes to seismicity rate and is combined with elastic Green's functions relating dike opening to deformation. Different space‐time patterns of seismicity develop if earthquakes occur at the same depth as the dike compared to above/below the dike. In the former, seismicity initiates near the dike's leading edges but shuts off as the dike tips pass and seismogenic volumes fall into stress shadows. In the latter, seismicity continues at a decaying rate after the tips pass. We focus on lateral propagation and develop a nonlinear inversion method that estimates dike length and pressure as a function of time. The method is applied to the 2007 Father's Day intrusion in Kilauea Volcano. Seismicity is concentrated at ∼3 km depth, comparable to geodetic estimates of dike depth, and decays rapidly in time. With lateral propagation only and a vertical dike‐tip line, it is difficult to fit both GPS data and the rapid down‐rift jump in seismicity, suggesting significant vertical propagation. For the events to have occurred below the dike requires a very short aftershock decay time, hence unreasonably high background stressing rate. The rapid decay is better explained if the dike extends somewhat below the seismicity. We suggest that joint inversion is useful for studying the diking process and may allow for improved short‐term eruption forecasts. Key Points Forward model is developed to predict dike induced deformation and seismicity. Inversion for dike length and pressure from deformation and seismicity data. Inversion method applied to 2007 Father's Day intrusion at Kilauea volcano.