CC BY
31
RUSSIAN JOURNAL OF EARTH SCIENCES, VOL. 16, ES3001, doi:10.2205/2016ES000568, 2016
Deformation precursors of large earthquakes derived from long term GNSS observation data
V. Kaftan1,2 and A. Melnikov2
Received 8 August 2016; accepted 20 August 2016; published 27 August 2016.
The results of observation of Global Navigation Satellite Systems (GNSS) in areas of large earthquakes are analyzed in this article. The Earth's surface deformation characteristics before, during and after the earthquake are researched, the results indicate the presence of abnormal deformation near their epicenters. The statistical evaluation of abnormal deformations with their root mean squares (rms) are presented. Conclusions about the possibilities of using local GNSS observation networks for evaluation of risk of strong seismic events are performed. KEYWORDS: Earthquake precursors; Earth's surface deformation; GNSS observation; San Andreas fault; least square method; seismic cycle related deformations; earthquake interaction, forecasting and prediction.
Citation: Kaftan, V. and A. Melnikov (2016), Deformation precursors of large earthquakes derived from long term GNSS observation data, Russ. J. Earth. Sci., 16, ES3001, doi:10.2205/2016ES000568.
Introduction
Researches of deformation earthquake precursors were on the front burner from the middle to the end of the previous century. The repeated conventional geodetic measurements such as precise leveling [Kaftan and Ostach, 1996] and linear-angular networks have been used for the study. Many examples of studies referenced to strong seismic events using conventional geodetic techniques are presented in [Rikitake, 1976]. One of the first case studies of geodetic earthquake precursors was done by Mescherikov [1968].
Rare repetitions, insufficient densities and locations of control geodetic networks made difficult predicting future places and times of earthquakes occurring.
Intensive development of Global Navigation Satellite Systems (GNSS) during the latest decades allows doing the research in a more effective level.
Investigation Area and Networks
Today permanent GNSS stations are being installed widely all over the world. It is possible now to study the Earth's surface deformation on a scale never possible before.
1 Geophysical Center of the Russian Academy of Sciences, Moscow, Russia
2Agro-Technological Institute, Peoples' Friendship University of Russia, Moscow, Russia
Copyright 2016 by the Geophysical Center RAS. http://elpub.wdcb.ru/journals/rjes/doi/2016ES000568-res.html
Some of permanent GNSS networks are covering the seismo-generating zones. One of the more investigated seismic areas is San Andreas Fault zone of California, USA [Wallace, 1990]. Two of GNSS networks of this zone are well placed to study Earth's surface deformation just near the epicenters of the strong Parkfield (September 28, 2004, Mw = 6.0) and El Mayor Cucapah (April 4, 2010, Mw = 7.2) earthquakes. The epicenters of the earthquakes are conveniently located several kilometers far from the permanent GNSS networks.
The Parkfield permanent GPS network of the Plate Boundary Observatory, USA, was used for the study of the Earth's surface deformation in relation to the Parkfield earthquake (Figure 1). The usage of the network for the study of seismic displacements is described in [Langbein and Bock, 2004].
The California Real Time Network (CRTN) of permanent GPS stations (Figure 2) is located not far from the epicenter of El Mayor Cucapah main shock and covers epicenters of some aftershocks of it. The Southern California Integrated GPS network is described particularly in [Hudnut et al., 2001].
The only shortest baseline vectors formed the Delauney triangulation were used in the processing as it is recommended in [Dokukin et al., 2010] and shown at Figure 1 and Figure 2.
Observation Data and Preprocessing
The block flow diagram of the analysis is shown at Figure 3.
Observation data used for the deformation analysis was received from the archive of the Scripps Orbit and Perma-
ES3001
1 of 9
■10000 -8000 -6000 -4000 -2000 0 2000 4000 6000 8000 km
Figure 1. Parkfield permanent GPS network used in the study.
nent Array Center (SOPAC) [http://sopac.ucsd.edu/]. Data sampling is equal to 30 s. Daily measurements were processed using MAGNET Tools software for the determination of baseline vectors and its co-variation matrices. Such processing was performed for the 6 measurement epochs for both GNSS networks. Then the obtained baseline vectors were adjusted using the special technique. The two dimensional approach of the deformation analysis has been used in the study. It is recommended as a proper and unified example in [Dermanis and Kotsakis, 2005].
terms of the normal equations, where Q+ = P is the matrix of the weights of the measurements.
This solution satisfies not only the condition v Pv = min but the xTQ+x = min too.
The adjustment procedure makes possible to calculate the plan deformation characteristics. Principal deformations , ^2 and dilatation A have been used in this study.
11 = [x2(dy3 - dyi) + y2(dx3 - dxi) - x^(dyi - dyi)-
Measurement Adjustment and Graphical Representation
The procedure of the adjustment is described as follows. Observation equations were presented as
v = Adx + I
(1)
where v - vector of the derived corrections to the differences of repeatedly measured baseline components of the order (3n — 3) X 1 for n baselines; A - matrix of the coefficients of the observation equation (1); dx - vector of point displacements of the order 3k X 1 in case of k determinate points; I - vector of the differences of the measured network elements of the order (3n — 3) X 1.
The Least Square solution of the observation equations (1)
dx = — N + L = — Qdx L
where N = AT Q+ A and L = AT Q+l, the so-called matrix of the normal equations coefficients and vector of the free
y3 (dX2 - dxi)]j (X2 y3 - X3 y2 ) 12 = [X2(dx3 - dxi) + y2(dy3 — dyi) - x3(dx2 - dxi)-y3(dy2 - dyi)}j(X2V3 - X3V2)
A = [xi(dy3 - dyi) - yi(dx3 - dxi) - X3(dyi - dyi)+
y3(dx2 - dxi)\j(x2y3 - X3yi)
where Xi, yi are plan coordinates; dxi, dyi - plan displacements; i - indices of vertex of triangles.
Displacements and deformations were calculated for the three epoch differences before and three epoch differences after the El Mayor Cucapah earthquake. Corresponding characteristics referenced to the Parkfield earthquake were received from the study [Kaftan et al., 2010]. The plan displacement vectors and vertical displacement isolines are presented in Figure 4 and Figure 5. The visual comparison of
Figure 2. CRTN network used in the study.
these characteristics does not allow to see some abnormal earth surface behavior before the earthquakes.
The graphical contour maps of the dilatation patterns and the main strain axes as a result of the "quick look" analysis are presented in Figure 6 and Figure 7 in comparison to each other. Figure 6 shows the existence of substantial Earth's surface deformations prior to the both earthquakes rising in time approaching to the moments of the main shocks. The
main extremes of deformations are placed near the earthquake epicenters. It is possible to consider that these features have to be considered as the earthquake precursors. The postseismic deformations are shown at Figure 7. It accelerates in the same locations in relation to the preseismic patterns.
The quantitative estimates of the dilatation A and root mean square errors aA of every observation interval for both study cases are represented in the Table 1. Preseismic and postseismic dilatation values are written in the top and bottom parts of the Table 1 respectively. The values of preseis-mic deformations vary from 0.1 to 1.8 x 10-5. The world practice demonstrates that so high level of deformation corresponds to seismic processes. These values of deformations can be considered as alarms of approaching to seismic event occurrences near its localizations.
In the case of the adopted adjustment technique it is possible to test the general statistical efficiency of displacement vectors determination. For this purpose we compute the ratio [Kaftan, 2003]
F =
dxT Ndx vT Pv
and set a confidence level a according to a freedom degree n — k with the use of F--distribution tables setting up the
Figure 3. Block flow diagram of the analysis.
critical value F(a,, shown in Table 2.
-fe).
The results of the testing are
fe
x 106
x 106
a)
I
Date 2007 1 31
- CS PftIATÎ , ppcv:
IBVHL
■cms -OAES
-W«T-WBC.KYVW
I
•®SSC -COTD ^rtjps
------MAJCACI
D) / I . ...... -P491
' X .pfcr \AZRY
^4®4 -««-G
•Surf4» ^^
>P474
'P4E2 'P4M
-RAAP
I
kesc
•P4«
-PS02
>p4sapjoi
3.986
3.75
x 10
xlo6
3.8
3.7
3.65
3.6
3.55
i 0.01
/ / (Date 2010 4 3
I
I 0.2
0.15
m
I 0.05
^POTR
10 cm
Alarm
4.5
5.5
6.5
C7"
7.5 x105
/ / /
' P47 2 NufKJP
e)
f)
SP4S6
^SIOS
/ Alarm
-C \
^POTR *-.PO«1
Alarm
<P494>4«
1. •
«
I
x 10s
Clear
* 10s
Figure 4. Preseismic horizontal (arrows) and vertical (contour lines) displacements before the Parkfield ( a, c, e ) and the El Mayor Cucapah ( b, d, f) earthquakes. Solid red dots are strongest foreshocks with M > 4. To zoom selected panel click on the corresponding number in red square. To remove zoomed image click on the "Clear" button.
Figure 5. Postseismic horizontal (arrows) and vertical (contour lines) displacements before the Parkfield (a, c, e) and the El Mayor Cucapah (b, d, f) earthquakes. Solid red dots are strongest shocks with M > 4. Stars are the main shock epicenters. To zoom selected panel click on the corresponding number in red square. To remove zoomed image click on the "Clear" button.
Figure 6. Preseismic principal strains (crosses of solid and dashed lines) dilatation patterns (contour lines) before the Parkfield ( a, c, e ) and the El Mayor Cucapah ( b, d, f) earthquakes. Solid red dots are strongest foreshocks with M > 4. To zoom selected panel click on the corresponding number in red square. To remove zoomed image click on the "Clear" button.
X 106
^XlO6
a)
X 10"
x 10 3.986 p
x106 3.8 -
b)
3.75
3.7
3.65
3.6
-OPCL
OPBL
*P623
-BLYT
GMPK
*POTR -P066
■ 10
*
0.5
0 X105
-0.5
1-1.5
Figure 7. Postseismic horizontal (arrows) and vertical (contour lines) displacements before the Parkfield (a, c, e) and the El Mayor Cucapah (b, d, f) earthquakes. Solid red dots are strongest shocks with M > 4. Stars are the main shock epicenters. To zoom selected panel click on the corresponding number in red square. To remove zoomed image click on the "Clear" button.
Table 1. Dilatation Extreme and Accuracy Estimations
Parkfield, September 2004, M = 6.0 El Mayor Cucapah, April 2010, M = 7.2
Time interval before Time interval from Time interval before Time interval from
and after the the first epoch A ± <7A(10-B) and after the the first epoch A ± <7A(10-B)
earthquake (years) (years) earthquake (years) (years)
-1.74 1.00 0.50 ± 0.04
-0.80 1.94 1.00 ± 0.07
-0.01 2.74 1.77 ± 0.06
+0.01 2.74 3.41 ± 0.04
+0.24 3.00 10.1 ± 0.06
+1.24 3.99 12.5 ± 0.06
-3.17 0.91 0.08 ± 0.01
-1.98 2.09 0.11 ± 0.01
-0.01 4.09 0.26 ± 0.01
+0.01 4.09 +1.00 ± 0.03
+1.74 5.83 +1.10 ± 0.01
+5.58 9.66 +1.40 ± 0.01
Table 2. Statistical Testing of Deformation by F and Range Criteria
Parkfield, September 2004, M = 6.0
El Mayor Cucapah, April 2010, M = 7.2
Time interval before the earthquake (years)
Time interval from the first epoch (years)
F
(fcrit =1.2,
P > 0.01, n - k = 288)
ÄA ± &AR
(10-5)
Time interval before the earthquake (years)
Time interval from the first epoch
F
(fcrit = 1.2,
P > 0.01,
ÄA ± &AR
(10-5)
(years) n - k = 288)
-1.74 1.00
-0.80 1.94
-0.01 2.74
+0.01 2.74
+0.24 3.00
+1.24 3.99
32.7 0.91 ± 0.06
106.5 1.83 ± 0.09
158.3 2.86 ± 0.11
2538.7 4976.6 8646.2 8.77 ± 0.08 15.5 ± 0.10 22.6 0.10
-3.17 0.91
-1.98 2.09
-0.01 4.09
+0.01 4.09
+1.74 5.83
+5.58 9.66
35.7 0.16 ± 0.02
38.6 0.35 ± 0.01
72.5 0.65 ± 0.02
22.5 2.51 ± 0.09
646.9 2.82 ± 0.02
415.0 3.16 0.03
As it is seen from the Table 2 all F ratios confidently exceed the critical values. It attests to the presumption that the level of displacements is much higher than the level of the measurement accuracy.
The range of extremal dilatation values _RA = Amax — Amin was used as other criterion of the abnormal deformation of the studied territory. It presented in the Table 2. It is possible to consider that the range values can be thought of as deformation precursors of the strong earthquake occurrence.
The further research will be continued to reveal the more detailed tendencies and regularities of Earth's surface deformation behavior in seismo-generating zones.
Acknowledgments. Authors are very grateful to GNSS webarchive and earthquake catalog creators for opened data access.
References
Conclusion
Two studied cases of horizontal deformation behavior before and after strong earthquakes show the contrast deformation changes from years by days before the events near their epicenters.
The estimated extremes can be considered as temporal precursors of the earthquake occurring.
It is possible that the places of the earthquake occurring can be predicted using the extremes too.
The permanent GNSS networks covering seismic generating fault zones can be effective tools for the earthquake prediction.
Dermanis, A., Ch. Kotsakis (2005), Estimating Crustal Deformation Parameters from Geodetic Data: Review of Existing Methodologies. Open Problems and New Challenges, International Association of Geodesy Symposia, 131, 7—18.
Dokukin, P. A., V. I. Kaftan, R. I. Krasnoperov (2010), The influence of the form of triangles in geodetic network on determining of Earth's surface deformation, Proceedings of the Univers. Geodesy and Aerophotosurveying, 5, 6—11. (in Russian)
Hudnut, K. W., et al. (2001), The Southern Califor-
nia Integrated GPS Network, Proceedings of the 10th FIG International Symposium on Deformation Measurements, 19-22 March 2001 p. 129-148, Orange, California, USA.
Kaftan, V. I. (2003), Temporal analysis of geospatial data: Kinematic models, Thesis for a doctor's degree, p. 285, Moscow
State University of Transportation, Moscow. (in Russian) Kaftan, V. I., O. M. Ostach (1996), Vertical land deformation in Caucasus region, Earthquake Prediction Research, 5, 235— 245.
Kaftan, V. I., R. I. Krasnoperov, P. P. Yurowsky (2010), Graphical representation of results of the earth's surface movements and deformations determination by means of global navigation satellite systems, Geodesy and Cartography, No. 11, 2—7. (in Russian)
Langbein, J., Y. Bock (2004), High-rate real-time GPS
network at Parkfield: Utility for detecting fault slip and seismic displacements, Geophysical Research Letters, 31, L15S20, doi:10.1029/2003gl019408 Mescherikov, J. A. (1968), Resent crustal movements in seismic
regions: Geodetic and geomorphic data, Tectonophysics, 6, 29, doi:10.1016/0040-1951(68)90024-3 Rikitake, T. (1976), Earthquake Prediction, 357 pp., Elsevier, New York.
Wallace, R. E., (ed.) (1990), The San Andreas fault system, California. U.S. Geological Survey professional paper 1515, 405 pp., Government Printing Office, Washington, D.C., U.S.
V. Kaftan, Organization of the Russian Academy of Sciences Geophysical Center, 3 Molodezhnaya Str., 119296 Moscow, Russia. ([email protected])
A. Melnikov, Agro-Technological Institute, Peoples' Friendship University of Russia, Moscow, Russia.
