Determination of the horizontal strain rates tensor in Western Ukraine

Abstract

Дані GNSS спостережень (CORS) з 37 станцій, розташованих у районі Західної України, обробленО за допомогою модуля Bernese Processing Engine (BPE) Бернського програмного забезпечення GNSS версії 5.2 протягом періоду часу близько 2,5 року. Щоб досягти кращої згоди, вибрано станції IGS, найближчі до навколишнього району дослідження, з фіксованими координатами ITRF2008 в епоху 2005.0. Східна та північна складові швидкості спостережень GNSS з цих 37 постійних станцій, обчислені за результатами вимірювань GNSS, використані для побудови двовимірної моделі поля горизонтальних деформацій цієї місцевості. Це дослідження складається з трьох частин. По-перше, проаналізовано два точні рішення для компонентів 2D тензора швидкостей деформацій, отримані на геосфері на основі розв’язання власних величин – задачі власних векторів, ураховуючи симетричний тензор швидкості обертання. По-друге, на основі найпростіших і найкорисніших формул з першого етапу виконано строге оцінювання точності компонентів 2D тензора швидкостей деформацій на основі правила поширення коваріацій. Нарешті, обчислено компоненти 2D тензора швидкості деформації, швидкості дилатації та компоненти тензора рівних швидкостей в області. Для описаної області побудовано модель тензора швидкості обертання. Це привело до висновку, що область дослідження слід інтерпретувати як деформовану територію. На основі обчислень з GNSS-моделі цих компонентів горизонтальних деформацій встановлено норми основних значень та швидкості основних осей деформації земної кори. Основні тектонічні утворення показано як фонову інтенсивність різних компонентів швидкостей, швидкість обертання та тензори швидкості деформації. Топографічні особливості регіону ґрунтувались на моделі SRTM-3 (місія з топографії Shuttle) з роздільною здатністю 3²¥3². На перший погляд, найбільші значення отримано в районах, розташованих навколо Українських Карпат. Швидкість дилатації також має подібний розподіл. Тим не менше, оскільки в роботі обчислено лише власні числа та власні вектори без оцінки точності, це може призвести до сумнівних висновків щодо інтерпретації результатів і вимагає додаткового розв’язання суто математичної задачі. Потрібно знайти коваріаційну матрицю тензора деформації на основі заданої коваріаційної матриці компонентів швидкості, одержаних програмним забезпеченням Bernese. Оскільки досліджуваний регіон є дуже складним, то за отриманими результатами необхідне подальше ущільнення перманентних станцій GNSS.Doppler Orbitography and Radio-positioning Integrated by Satellite (CORS) observations from 37 Global Navigation Satellite System (GNSS) stations located in the Western Ukraine area were processed using Bernese Processing Engine module (BPE) of Bernese GNSS Software version 5.2 for a time span of about 2.5 years. To get a better agreement for constrains, the IGS stations closest to the surrounding area of study were chosen with fixed coordinates of ITRF2008 at epoch 2005.0. Eastern and Northern components of velocities of GNSS observations from these 37 permanent stations, calculated from GNSS measurements, were used to construct a 2D model of horizontal strain rates field for the area. This study is presented in three parts. Firstly, two exact solutions for the components of the 2D strain rate tensor derived on the geosphere based on solving the eigenvalues – eigenvectors problem were analyzed, including skew symmetric rotational rate tensor. Secondly, based on the most simple and useful formulas from the first stage, a rigorous estimation of the accuracy of components of the 2D strain rate tensor were obtained based on the covariance propagation rule. Finally, the components of the 2D strain rate tensor, dilatation rate and components of the sheer rate tensor in the region were computed. A model of the rotation rate tensor was constructed for the described area, which led to the conclusion that the region of study should be interpreted as a deformed territory. Based on the computations from the GNSS-data model of components of horizontal deformations, the rates of principal values and rates of principal axes of the Earth’s crust deformation were found. To be consistent, the main tectonic formations are shown as the background intensity of different components of velocities, the rotation rate and strain rate tensors. Topographic features of the region were based on the SRTM-3 model (Shuttle Radar Topography Mission) with resolution 3²¥3². At the first sight, the maximum sheer rates have greatest values in the areas located around the Ukrainian Carpathians. The dilatation rate has also a similar distribution. Nevertheless, because in the paper only eigenvalue – eigenvector problem without accuracy estimation has been considered, which possibly leads to doubtful conclusions regarding interpretation and requires an additional solution of a purely mathematical problem. The full covariance matrix of the strain rate tensor should be found based on given full covariance matrix of the velocity components obtained by Bernese software. As a matter of fact, the study region is very complex in terms of crustal movements, which, according to the results obtained, require further densification of permanent GNSS stations.