PlatEMotion: a Matlab® Tool for Geodetic Reference Frame Definition

Istituto Nazionale di Geofisica e Vulcanologia

An important concept related to plate tectonics theory is the internal rigidity of tectonic plates together with Euler’s theorem, which allows modelling the relative motion of plates quantitatively. Roughly speaking, plates can be viewed as rigid caps on the surface of a sphere. The motion of a plate on a sphere can be described by a rotation about a virtual axis that passes through the centre of the sphere (Euler’s theorem). This implies that an angular velocity vector originating at the centre of the Earth can describe motions of plates. The most widespread parameterization of such a vector is using latitude and longitude to describe the location where the rotation axis cuts the surface of the Earth, and a rotation rate that corresponds to the magnitude of the angular velocity. Formally, the latitude and longitude of the angular velocity vector constitute the so called ‘Euler pole’.
Here, we describe the main features of the developed software tool, called PlatE-Motion (PEM), which allows calculating both the expected velocity value for any point located on the Earth providing a Euler pole (direct problem) and inferring the Euler pole parameters by inverting the observed velocities at a set of sites located on a rigid block (inverse problem). PEM is developed in Matlab® framework and is designed for easy-to-use file exchange with the GAMIT/GLOBK software package [Herring et al., 2006], and has a user-friendly graphical user interface (GUI). The file formats of the data used by PEM are not only restricted to the GAMIT/GLOBK file formats, but they can also be simple ASCII files generated from other kind of software. The tool is available for the scientific community. The interested user can freely obtain it by simply contacting the authors

Herring, T.A., King, R.W. and McKlusky, S.C., (2006). Introduction to GAMIT/GLOBK Release 10.3. Massachusetts Institute of Technology, Cambridge, M.A., pp. 45.

Mulargia, F. and Gasperini P., (1995). Evaluation of the applicability of the time and slip-predictable earthquake recurrence models to Italian seismicity. Geophys. J. Int., 120, 453-473.

Nocquet, J.-M., Calais, E, Altamimi, Z, Sillard P. and Boucher C., (2001). Intraplate deformation in the western Europe deduced from an analysis of the International Terestrial Reference Frame 1997 (ITRF97) velocity field. J. Geophys. Res., 106, B6, 11239-11257.

Palano, M., Puglisi, G. and Gresta, S., (2008). Ground deformation patterns at Mt. Etna from 1993 to 2000 from joint use of InSAR and GPS techniques. J. Volcanol. Geotherm. Res., 169, doi:10.1016/ j.jvolgeores.2007.08.014.

Palano, M., Rossi M., Cannavò F., Bruno, V., Aloisi M., Pellegrino D., Pulvirenti M., Siligato G. and Mattia, M., 2010. Etn@ref: a geodetic reference frame for Mt. Etna GPS networks. Annals of geophysics, 53, 4, doi: 10.4401/ag-4879.

Tarantola, A. and Vallette, B., (1982). Generalized nonlinear inverse problems solved using the least squares criterion. Rev. Geophys. Space Phys., 20, 219-232.