On inverting gravity changes with the harmonic inversion method: Teide (Tenerife) case study (original) (raw)

Here we investigate the applicability of the harmonic inversion method to time-lapse gravity changes observed in volcanic areas. We carry out our study on gravity changes occured over the period of 2004–2005 during the unrest of the Central Volcanic Complex on Tenerife, Canary Islands. The harmonic inversion method is unique in that it calculates the solution of the form of compact homogeneous source bodies via the mediating 3-harmonic function called quasigravitation. The latter is defined in the whole subsurface domain and it is a linear integral transformation of the surface gravity field. At the beginning the seeds of the future source bodies are introduced: these are quasi-spherical bodies located at the extrema of the quasigravitation (calculated from the input gravity data) and their differential densities are free parameters preselected by the interpreter. In the following automatic iterative process the source bodies change their size and shape according to the local values...