Main Article Content
We propose a MATLAB toolbox for the computation of the strain rate field from the coordinate time series of some continuous GNSS stations. It consists of several functions, also compatible with GNU Octave, implementing the following steps: (i) time series download from a data repository (e.g., the Nevada Geodetic Laboratory database); (ii) calculation of velocities of the selected stations by means of the Maximum Likelihood Estimation (MLE) method implemented in the external package Hector, including modeling of offsets, outliers, noise and periodic components; (iii) (optional) filtering of Common Mode Errors; (iv) calculation of the strain rate field with the modified least squares method, in which a scale factor can be introduced to define the locality of the deformation analysis and, besides uncertainty estimation, a geometric evaluation of the significance of the results is provided; (v) visualization of the results for immediate use and easy interpretation for scientific purposes. The toolbox is divided into two components: the first one, named StaVel, performs the steps (i)-(iii) and the second component, GridStrain, performs the steps (iv) and (v). The potential of the toolbox is demonstrated on a real dataset. Time series from several continuous GNSS stations in South-Eastern Sicily (Southern Italy) are processed by means of StaVel and GridStrain in order to provide the strain rate field.
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.
No Permission Required
Istituto Nazionale di Geofisica e Vulcanologia applies the Creative Commons Attribution License (CCAL) to all works we publish.
Under the CCAL, authors retain ownership of the copyright for their article, but authors allow anyone to download, reuse, reprint, modify, distribute, so long as the original authors and source are cited. No permission is required from the authors or the publishers.