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 km3), about 880 m beneath the volcanoʹs summit.