Relocating earthquakes in clusters based on variations in the intervals between their first P- and S-waves

Abstract

Довжину інтервалу між першими вступами Р- та S-хвиль часто використовують для приблизного оцінювання епіцентральної відстані. Ми пропонуємо алгоритм визначення координат землетрусів шляхом одночасного порівнювання великої кількості таких інтервалів. Для визначення різниці між інтервалами на кожній зі станцій обчислюється функція взаємної кореляції між відповідними хвильовими формами - без визначення абсолютної довжини інтервалів. У поточній версії алгоритму припускаємо, що розміри кластера набагато менші за відстань до станцій; азимути станцій і кути виходу перших Р- та S-хвиль доволі точно відомі принаймні для одного (опорного) землетрусу; промені перших хвиль лежать у вертикальній площині, що містить вогнище і станцію. Зв'язок між координатами і різницею довжини інтервалів у цьому разі стає суто геометричний і лінійний, а відповідну систему рівнянь легко розв'язати. Результати серії модельних експериментів з використанням різної кількості станцій та їхньої конфігурації, рівня шуму у спостережених даних, ступеня неповноти даних, неточних азимутів та кутів виходу довели стійкість і надійність алгоритму і можливість його застосування у подальшому до реальних даних. З огляду на велику кількість обмежень на координати кожного із землетрусів алгоритм найдоцільніше використовувати у разі дуже слабких землетрусів, або малої кількості станцій, коли істотна частина даних відсутня. Алгоритм можна використовувати окремо, або з метою підтвердження правильності координат, визначених іншими методами, або вбудувати його в інші методи, аби підвищити їх надійність завдяки істотному збільшенню кількості обмежень.The length of the interval between the first P- and S-waves is routinely used as a rough estimator of epicentral distance. We propose an algorithm for the relocation of earthquakes occurring in clusters, based on the simultaneous comparison of a large number of intervals. Variations in the intervals at each station are measured by cross-correlation between the respective portions of records directly and without a reference to any absolute times. In the current version of the algorithm, it is assumed that the size of the cluster is much smaller than the distance to the stations; the azimuths of the stations, as well as the angles of the emergence of the first P- and S-waves, are more or less accurately known for at least one (reference) earthquake; and the rays of the first waves lie in the vertical plane that contains the earthquake and the station. Under these assumptions, the relationship between the locations and the variations in the intervals becomes purely geometrical and linear, and the corresponding system can easily be solved. A series of synthetic experiments with different numbers and configurations of stations, levels of noise in the observed data, sparse data, and inaccuracies in azimuths and angles of emergence have demonstrated the stable and reliable performance of the algorithm and its potential applicability to real data. Due to the large number of constraints on each location, the algorithm can be used primarily in the case of small earthquakes or sparse networks when a large portion of data is missing. It can be used independently, to validate the locations determined by other methods, or be integrated into them, thereby improving their reliability by providing a large number of additional constraints.