dMODELS: A MATLAB software package for modeling crustal deformation near active faults and volcanic centers

We have developed a MATLAB software package for the most common models used to interpret deformation measurements near faults and active volcanic centers. The emphasis is on analytical models of deformation that can be compared with data from the Global Positioning System (GPS), InSAR, tiltmeters and strainmeters. Source models include pressurized spherical, ellipsoidal and sill-like magma chambers in an elastic, homogeneous, flat half-space. Dikes and faults are described following the mathematical notation for rectangular dislocations in an elastic, homogeneous, flat half-space. All the expressions have been checked for typographical errors that might have been present in the original literature, extended to include deformation and strain within the Earth's crust (as opposed to only the Earth's surface) and verified against finite element models. A set of GPS measurements from the 2006 eruption at Augustine Volcano (Alaska) is used to test the software package. The results show that the best fit source to the GPS data is a spherical intrusion (ΔV=5×10−4ΔV=5×10−4 km3), about 880 m beneath the volcano's summit.

► The software is a practical tool to investigate volcano and earthquake deformation. ► Emphasis is on analytical models of crustal deformation. ► Formulas were checked for typographical errors existing in the original references. ► All analytical models were verified against numerical models. ► The software package is freely available and open source.