Multiple analysis of GPS signal fluctuation parameters before and after the megaearthquake in Japan on 11 March 2011

Category: 16-1
A.A. Lyubushin, P.V. Yakovlev, E.A. Rodionov


UDC 550.334



 A.A. Lyubushin (1), P.V. Yakovlev (2), E.A. Rodionov (2)


(1) Schmidt Institute of Physics of the Earth, Russian Academy of Sciences, Moscow, Russia

(2) Ordzhonikidze Russian State Geological Prospecting University, Moscow, Russia


Abstract. The field of GPS signals recorded at the net of 1203 stations in the Japanese archipelago from 30 January up to 26 March of 2011 is regarded. This time interval includes over 40 days before the mega-earthquake on 11 March 2011 (M=9.0) and almost 16 days after it. Signals from each station represent 3-component time series with a 30 minutes sampling time interval. Statistical properties of random fluctuations of GPS signals after linear trend removal and coming to increments before and after a seismic event are investigated. The maps of three noise waveform properties are regarded: maximum normalized eigenvalue of the correlation matrix and 2 smoothness indexes built by using orthogonal wavelets of Daubechies and Lang. The maps are constructed as values within nodes of a regular grid of the size 50´50 nodes which covers the investigated region by using information from 10 stations which are nearest to each node. Thus, 18 matrices of the noise properties within nodes of the grid are obtained: by 9 matrices before and after an event for three components of GPS signals (N – displacement northward, E – displacements eastward and Z – displacements upward). The final analysis consists in computing of the 1st principal components of investigated noise properties and constructing the maps of values of mP – the maximum normalized eigenvalue of the correlation matrix of principal components within a sliding spatial rectangular window. These operations are performed separately using data before and after the earthquake.

                The results of analysis reveal that before the earthquake the values of mP exceeds significantly those after the event and the region of maximum values of mP includes the epicenter of the future earthquake whereas after it the vicinity of the epicenter becomes a region of small correlations. These results confirm the hypothesis about increased correlations of geophysical field noises within zones of preparation of larger earthquakes and they could be used for seeking earthquakes precursors by GPS data.


Keywords: GPS signals, correlation analysis, principal component, wavelets, earthquake precursors.


Aivazyan S.A., Bukhshtaber V.M., Enyukov I.S., and Meshalkin L.D. Prikladnaya statistika. Klassifikatsiya i snizhenie razmernosti (Applied statistics. Classification and dimension reduction), Moscow: Finansy i statistika, 1989.

Farkov Yu. A. Orthogonal wavelets with compact supports on locally compact Abelian groups, Izv. Ross. Akad. Nauk Ser. Mat. 2005, vol. 69, no. 3, pp. 193–220

Lang W.C. Fractal multiwavelets related to the Cantor dyadic group, Intern. J. Math. and Math. Sci. 1998, vol. 21, pp. 307–317.

Lyubushin A.A., Seismic catastrophe in Japan on March 11, 2011: Long-term forecast based on low-frequency microseisms, Geofiz. Protsessy Biosfera, 2011b, vol. 10, no. 1, pp. 9–35.

Lyubushin A.A., Forecast of the Great Japan earthquake, Priroda (Moscow, Russ. Fed.), 2012a, no. 8, pp. 23–33.

Lyubushin A.A., Mapping the properties of low-frequency microseisms for seismic hazard assessment, Izv., Phys. Solid Earth, 2013a, vol. 49, no. 1, pp. 9–18.

Lyubushin A.A., Maps of linear synchronization for the properties of global low-frequency seismic noise, Geofiz. Issled., 2013b, vol. 14, no. 1, pp. 41–53.

Lyubushin A. Prognostic properties of low-frequency seismic noise, Natural Science. 2012, vol. 4, no. 8A, pp. 659–666.

Lyubushin A. How soon would the next mega-earthquake occur in Japan? Natural Science. 2013, vol. 5, no. 8A1, pp. 1–7.

Lyubushin A.A. Dynamic estimate of seismic danger based on multifractal properties of low-frequency seismic noise, Natural Hazards. 2014, vol. 70, Is. 1, pp. 471–483.

Mallat S. A wavelet tour of signal processing. San Diego; London; Boston; N.Y.; Sydney; Tokyo; Toronto: Academic Press, 1998. 577 p.

Rodionov E.A. and Farkov Yu.A. Estimates of the smoothness of dyadic orthogonal wavelets of Daubechies type, Math. Notes, 2009, vol. 86, no. 3, pp. 429–444.

Stroganov S.A. Estimate of smoothness of low-frequency microseismic oscillations using dyadic wavelets, Geofiz. Issled., 2012, vol. 13, no. 1., pp. 17–22.