Matlab scripts for calculating ground surface displacements due to thermoelastic contraction. The model was developed by Furuya (2005), who gave the analytical quasi-static solution of ground deformation in response to an instantaneous heat source in an elastic half-space. The model can be applied to infer post-eruptive deformation after magma intrusion or lava emplacement.
Please cite the following papers when using the code:
- Masato Furuya (2005), Quasi-static thermoelastic deformation in an elastic half-space: theory and application to InSAR observations at Izu-Oshima volcano, Japan, Geophysical Journal International, 161(1), 230–242. doi:10.1111/j.1365-246X.2005.02610.x.
- Xiaowen Wang & Yosuke Aoki (2019). Posteruptive thermoelastic deflation of intruded magma in Usu volcano, Japan, 1992–2017, Journal of Geophysical Research: Solid Earth,124. doi.org/10.1029/2018JB016729