MATLAB Tools for Earth’s Surface Deformation Studies
MATLAB Tools for Earth’s Surface Deformation Studies
French Title: Des Outils MATLAB pour les Études de Déformation de la Croûte Terrestre
Previously Published Material: Some parts of this poster were already published in GPS Solutions with DOI 10.1007/s10291-012-0296-2, and 10.1007/s10291-013-0354-4, and some parts are presented here for the first time. In this poster, however, we present them as a software suite and focus on their computational aspects rather than the mathematical background, as well as a case study for eastern Canada.
Abstract ID#: 33493
English Abstract:
We developed three software programs for research on Earth's surface deformation studies, including time series analysis, Euler pole estimation, and strain analysis. The first software named GPS Interactive Time Series Analysis (GITSA) is capable of visualizing the input time series, determining and removing jumps and outliers, data interpolation, calculating quality indicators and producing publication quality graphical outputs. Along with computing basic statistics such as mean and variance, it can accomplish the bivariate statistical analysis, residual analysis, and spectral analysis. The second application named Euler Pole Calculator (EPC), estimates the expected velocities for any points on Earth’s surface given the relevant Euler pole parameters. It also allows for calculating the Euler pole parameters from observed velocities of a set of sites on the same tectonic plate. The last program named GeoStrain, calculates the crustal strain rates using the least-squares collocation method. It is able to simultaneously determine the signal and noise of the velocities at the observation points with the best possible removal of the observation errors, or any other position with no velocity observation. Furthermore, it allows for calculating strain and rotation rate tensors at any points of interest. The advantage of the software is the fact that it can optionally consider the effect of the vertical velocities on the strain rates.
Although our programs are independent, they make a software suite: the output of each program is the input of another one. While inputs of the first program are the position time series of continuous GNSS sites obtained from a GNSS data processing software, the outputs of the last one are parameters of deformation, e.g., in the form of strain tensors. We developed these programs using MATLAB programming language with a graphical user interface, and released them freely as open source software. MATLAB Parallel Computing Toolbox was used to increase the performance of computationally intensive algorithms on multi-core processors.
Although our programs are independent, they make a software suite: the output of each program is the input of another one. While inputs of the first program are the position time series of continuous GNSS sites obtained from a GNSS data processing software, the outputs of the last one are parameters of deformation, e.g., in the form of strain tensors. We developed these programs using MATLAB programming language with a graphical user interface, and released them freely as open source software. MATLAB Parallel Computing Toolbox was used to increase the performance of computationally intensive algorithms on multi-core processors.