“…Heterogeneous crustal and thermal properties in the vicinity of a magmatic system are known to partition the resulting deformation field (Del Negro et al, ; Gottsmann & Odbert, ; Gregg et al, ; Hickey et al, , ), and so a natural continuation of this work is to consider the comparison of viscoelastic rheologies and deformation modes in models that more closely represent reality. This can be carried out with models specific to volcanic systems containing heterogeneous crustal properties, provided by seismic tomography, spatially variable rheological effects through the use of viscoelastic shells within an elastic medium (e.g., Currenti, ; Currenti & Williams, ; Delgado et al, ; Newman et al, ; Segall, ), or a temperature‐dependent viscosity distribution that accounts for crustal geotherms and the perturbation owing to a modeled magmatic source (Del Negro et al, ; Gottsmann et al, ; Gottsmann & Odbert, ; Gregg et al, ; Hickey et al, ). Further to this, an important consideration is the influence of regional stresses and strain fields (Costa et al, ; Currenti & Williams, ), such as active extension within the Taupo Volcanic Zone (e.g., Cabaniss et al, ), on the observed viscoelastic behaviors and the resultant deformation time series.…”