Analysis of Swarm Satellite Magnetic Field Data Before the 2016 Ecuador (Mw = 7.8) Earthquake Based on Non-negative Matrix Factorization

In this paper, based on non-negative matrix factorization (NMF), we analyzed the ionosphere magnetic field data of the Swarm Alpha satellite before the 2016 (Mw = 7. 8) Ecuador earthquake (April 16, 0.35°N, 79.93°W), including the whole data collected under quiet and disturbed geomagnetic conditions. The data from each track were decomposed into basis features and their corresponding weights. We found that the energy and entropy of one of the weight components were more concentrated inside the earthquake-sensitive area, which meant that this weight component was more likely to reflect the activity inside the earthquake-sensitive area. We focused on this weight component and used five times the root mean square (RMS) to extract the anomalies. We found that for this weight component, the cumulative number of tracks, which had anomalies inside the earthquake-sensitive area, showed accelerated growth before the Ecuador earthquake and recovered to linear growth after the earthquake. To verify that the accelerated cumulative anomaly was possibly associated with the earthquake, we excluded the influence of the geomagnetic activity and plasma bubble. Through the random earthquake study and low-seismicity period study, we found that the accelerated cumulative anomaly was not obtained by chance. Moreover, we observed that the cumulative Benioff strain S, which reflected the lithosphere activity, had acceleration behavior similar to the accelerated cumulative anomaly of the ionosphere magnetic field, which suggested that the anomaly that we obtained was possibly associated with the Ecuador earthquake and could be described by one of the Lithosphere–Atmosphere–Ionosphere Coupling (LAIC) models.

This research was supported by the National Natural Science Foundation of China under Grant No. 41974084 and the International Cooperation Project of Department of Science and Technology of Jilin Province No. 20200801036GH.

Agurto-Detzel, H., Font, Y., Charvis, P., Régnier, M., Rietbrock, A., Ambrois, D., et al. (2019). Ridgesubduction and afterslip control aftershock distribution of the 2016 Mw 7.8 Ecuador earthquake. Earthand Planet. Sci. Lett. 520, 63–76. doi: 10.1016/j.epsl.2019.05.029

CrossRef Full Text | Google Scholar

Akhoondzadeh, M., De Santis, A., Marchetti, D., Piscini, A., and Cianchini, G. (2018). Multi precursorsanalysis associated with the powerful Ecuador (Mw= 7.8) earthquake of 16 April 2016 using Swarmsatellites data in conjunction with other multi-platform satellite and ground data. Adv. Space Res. 61, 248–263. doi: 10.1016/j.asr.2017.07.014

CrossRef Full Text | Google Scholar

Akhoondzadeh, M., De Santis, A., Marchetti, D., Piscini, A., and Jin, S. (2019). Anomalous seismo-laivariations potentially associated with the 2017 Mw=7.3 Sarpol-e Zahab (Iran) earthquake from Swarmsatellites, GPS-TEC and climatological data. Adv. Space Res. 64, 143–158. doi: 10.1016/j.asr.2019.03.020

CrossRef Full Text | Google Scholar

Alken, P. (2016). Observations and modeling of the ionospheric gravity and diamagnetic current systemsfrom CHAMP and Swarm measurements. J. Geophys. Res. Space Phys. 121, 589–601. doi: 10.1002/2015ja022163

CrossRef Full Text | Google Scholar

Benioff, H. (1949). Seismic evidence for the fault origin of oceanic deeps. Bull. Geol. Soc. Am. 60, 1837–1856. doi: 10.1130/0016-7606(1949)60[1837:SEFTFO]2.0.CO;2

CrossRef Full Text | Google Scholar

Bolt, B. A. (1999). Earthquake. NewYork, NY: W. H. Freeman.

Google Scholar

Bouffard, J., Floberghagen, R., and Olsen, N. (2019). The Swarm Satellite Trio Studies Earth and Its Environment. EOS.

Google Scholar

Cichocki, A., Zdunek, R., and Amari, S. (2006). “New algorithm for non-negative matrix factorization inapplication to blind source separation,” in Proceedings of ICASSP'06 (Toulouse: IEEE Press), 621–624.

Google Scholar

De Santis, A., Abbattista, C., Alfonsi, L., Amoruso, L., Campuzano, S. A., Carbone, M., et al. (2019a). Geosystemics view of earthquakes. Entropy 21:412. doi: 10.3390/e21040412

PubMed Abstract | CrossRef Full Text | Google Scholar

De Santis, A., Balasis, G., Pavón-Carrasco, F. J., Cianchini, G., and Mandea, M. (2017). Potential earthquake precursory pattern from space: the 2015 Nepal event as seen by magnetic Swarm satellites. Earth Planet. Sci. Lett. 461, 119–126. doi: 10.1016/j.epsl.2016.12.037

CrossRef Full Text | Google Scholar

De Santis, A., De, F. G., Spogli, L., Perrone, L., Alfonsi, L., Qamili, E., et al. (2015). Geospaceperturbations induced by the earth: the state of the art and future trends. Phys. Chem. Earth 85–86, 17–33. doi: 10.1016/j.pce.2015.05.004

CrossRef Full Text | Google Scholar

De Santis, A., Marchetti, D., Pavon-Carrasco, F. J., Cianchini, G., Perrone, L., Abbattista, C., et al. (2019b). Precursory worldwide signatures of earthquake occurrences on Swarm satellite data. Sci. Rep. 9:20287. doi: 10.1038/s41598-019-56599-1

PubMed Abstract | CrossRef Full Text | Google Scholar

De Santis, A., Marchetti, D., Spogli, L., Cianchini, G., Pavón-Carrasco, F. J., Franceschi, G. D., et al. (2019c). Magnetic field and electron density data analysis from Swarm satellites searching for ionosphericeffects by great earthquakes: 12 case studies from 2014 to 2016. Atmosphere 10:371. doi: 10.3390/atmos10070371

CrossRef Full Text | Google Scholar

D'Errico, J. (2021). SLM - Shape Language Modeling, MATLAB Central File Exchange. Available online at: https://www.mathworks.com/matlabcentral/fileexchange/24443-slm-shape-language-modeling (accessed February 23, 2021).

Dobrovolsky, I. R., Zubkov, S. I., and Myachkin, V. I. (1979). Estimation of the size of earthquakepreparation zones. Pure Appl. Geophys. 117, 1025–1044. doi: 10.1007/BF00876083

CrossRef Full Text | Google Scholar

Finlay, C. C., Olsen, N., Kotsiaros, S., Gillet, N., and Tøffner-Clausen, L. (2016). Recent geomagneticsecular variation from Swarm and ground observatories as estimated in the CHAOS-6 geomagnetic fieldmodel. Earth Planets Space 68:112. doi: 10.1186/s40623-016-0486-1

CrossRef Full Text | Google Scholar

Finlay, C. C., Olsen, N., and Tøffner-Clausen, L. (2015). DTU candidate field models for IGRF-12 and theCHAOS-5 geomagnetic field model. Earth Planets Space 67:114. doi: 10.1186/s40623-015-0274-3

CrossRef Full Text | Google Scholar

Freund, F. (2000). Time-resolved study of charge generation and propagation in igneous rocks. J. Geophys. Res. Solid Earth 105, 11001–11019. doi: 10.1029/1999JB900423

CrossRef Full Text | Google Scholar

Freund, F. (2011). Pre-earthquake signals: underlying physical processes. J. Asian Earth Sci. 41, 383–400. doi: 10.1016/j.jseaes.2010.03.009

CrossRef Full Text | Google Scholar

Friis-Christensen, E., Lühr, H., and Hulot, G. (2006). Swarm: a constellation to study the Earth's magneticfield. Earth Planets Space 58, 351–358. doi: 10.1186/BF03351933

CrossRef Full Text | Google Scholar

Friis-Christensen, E., Lühr, H., Knudsen, D., and Haagmans, R. (2008). Swarm - an earth observationmissioninvestigatinggeospace. Adv. Space Res. 41, 210–216. doi: 10.1016/j.asr.2006.10.008

CrossRef Full Text | Google Scholar

Fu, C., Lee, L., Ouzounov, D., and Jan, J. (2020). Earth's outgoing longwave radiation variability prior to M ≥ 6.0earthquakes in the Taiwan area during 2009-2019. Front. Earth Sci. 8:364. doi: 10.3389/feart.2020.00364

CrossRef Full Text | Google Scholar

Han, P., Hattori, K., Hirokawa, M., Zhuang, J., Chen, C., Febriani, F., et al. (2014). Statistical analysis ofULF seismomagnetic phenomena at Kakioka, Japan, during 2001-2010. J. Geophys. Res. Space Phys. 119, 4998–5011. doi: 10.1002/2014JA019789

CrossRef Full Text | Google Scholar

Hattori, K., and Han, P. (2018). “Statistical analysis and assessment of ultralow frequency magnetic signals inJapan as potential earthquake precursors,” in Pre-earthquake Processes: A Multidisciplinary Approach toEarthquake Prediction Studies, eds D. Ouzounov, S. Pulinets, K. Hattori, and P. Taylor (Hoboken, NJ: John Wiley & Sons), 229–240.

Google Scholar

Hattori, K., Serita, A., Gotoh, K., Yoshino, C., Harada, M., Isezaki, N., et al. (2004). ULF geomagneticanomaly associated with 2000 Izu Islands earthquake swarm, Japan. Phys. Chem. Earth Parts A B C 29, 425–435. doi: 10.1016/j.pce.2003.11.014

CrossRef Full Text | Google Scholar

Hayakawa, M. (2011). On the fluctuation spectra of seismo-electromagnetic phenomena. Nat. Hazards Earth Syst. Sci. 11, 301–308. doi: 10.5194/nhess-11-301-2011

CrossRef Full Text | Google Scholar

Hayakawa, M., and Molchanov, O. A. (2002). Seismo Electromagnetics Lithosphere-Atmosphere-Ionosphere Coupling. Tokyo: TERRA-PUB.

Google Scholar

Ho, Y., Jhuang, H., Lee, L., and Liu, J. (2018). Ionospheric density and velocity anomalies before M ≥ 6.5 earthquakes observed by DEMETER satellite. J. Asian Earth Sci. 166, 210–222. doi: 10.1016/j.jseaes.2018.07.022

CrossRef Full Text | Google Scholar

Kuo, C. L., Lee, L. C., and Huba, J. D. (2014). An improved coupling model for the lithosphere-atmosphere-ionosphere system. J. Geophys. Res.:Space Phys. 119, 3189–3205, doi: 10.1002/2013JA019392

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, D. D., and Seung, H. S. (1999). Learning the parts of objects by non-negative matrix factorization. Nature 401, 788–791.

PubMed Abstract | Google Scholar

Lee, D. D., and Seung, H. S. (2000). “Algorithms for non-negative matrix Factorization,” in International Conference on Neural Information Processing Systems, NIPS (Denver, CO: MIT Press), 556–562.

Li, M., Lu, J., Zhang, X., and Shen, X. (2019). Indications of ground-based electromagnetic observationsto a possible lithosphere-atmosphere-ionosphere electromagnetic coupling before the 12 May 2008 Wenchuan Ms 8.0 earthquake. Atmosphere 10:355. doi: 10.3390/atmos10070355

CrossRef Full Text | Google Scholar

Liu, J. Y., Chen, Y. I., Huang, C. C., Parrot, M., Shen, X. H., Pulinets, S. A., et al. (2015). A spatial analysison seismo-ionospheric anomalies observed by DEMETER during the 2008 M8.0 Wenchuan earthquake. J. Asian Earth Sci. 114, 414–419. doi: 10.1016/j.jseaes.2015.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Marchetti, D., and Akhoondzadeh, M. (2018). Analysis of Swarm satellites data showing seismo-ionosphericanomalies around the time of the strong Mexico (mw=8.2) earthquake of 08 September 2017. Adv. Space Res. 62, 614–623. doi: 10.1016/j.asr.2018.04.043

CrossRef Full Text | Google Scholar

Marchetti, D., De Santis, A., D'Arcangelo, S., Poggio, F., Jin, S., Piscini, A., et al. (2019a). Magneticfield and electron density anomalies from Swarm satellites preceding the major earthquakes of the2016-2017 Amatrice-Norcia (central italy) seismic sequence. Pure Appl. Geophys. 177, 305–319. doi: 10.1007/s00024-019-02138-y

CrossRef Full Text | Google Scholar

Marchetti, D., De Santis, A., D'Arcangelo, S., Poggio, F., Piscini, A. A., et al. (2019b). Pre-earthquake chain processes detected from ground to satellite altitude in preparation of the 2016-2017seismic sequence in central Italy. Remote Sens. Environ. 229, 93–99. doi: 10.1016/j.rse.2019.04.033

CrossRef Full Text | Google Scholar

Marchetti, D., De Santis, A., Shen, X., Campuzano, S. A., Perrone, L., Piscini, A., et al. (2019c). Possiblelithosphere-atmosphere-ionosphere coupling effects prior to the 2018 Mw=7.5 Indonesia earthquakefrom seismic, atmospheric and ionospheric data. J. Asian Earth Sci. 188:104097. doi: 10.1016/j.jseaes.2019.104097

CrossRef Full Text | Google Scholar

Mouri, M., Funase, A., Cichocki, A., Takumi, I., Yasukawa, H., and Hata, M. (2009). “Implementation of matrix factorization based onminimizing quasi-absolute distance for electromagnetic global signal elimination,” in 17th Europe SignalProcessing Conference (Glasgow: EUSIPCO), 24–28.

Google Scholar

Natarajan, V., and Philipoff, P. (2018). Observation of surface and atmospheric parameters using “NOAA 18”satellite: a study on earthquakes of Sumatra and Nicobar Is regions for the year 2014 (M ≥ 6.0). Nat. Hazards 92, 1097–1112. doi: 10.1007/s11069-018-3242-y

CrossRef Full Text | Google Scholar

Olsen, N., Friis-Christensen, E., Floberghagen, R., Alken, P., Beggan, C. D., Chulliat, A., et al. (2013). The Swarm Satellite Constellation Application and Research Facility (SCARF) and Swarm data products. Earth Planets Space 65, 1189–1200. doi: 10.5047/eps.2013.07.001

CrossRef Full Text | Google Scholar

Olsen, N., Haagmans, R., Sabaka, T. J., Kuvshinov, A., Maus, S., Purucker, M. E., et al. (2006). The SwarmEnd-to-End mission simulator study: a demonstration of separating the various contributions to Earth'smagnetic field using synthetic data. Earth Planets Space 60, 359–370. doi: 10.1186/BF03351934

CrossRef Full Text | Google Scholar

Paatero, P., and Tapper, U. (1994). Positive matrix factorization: a non-negative factor model with optimalutilization of error estimates of data values. Environmetrics 5, 111–126. doi: 10.1002/env.3170050203

CrossRef Full Text | Google Scholar

Park, J., Noja, M., Stolle, C., and Lühr, H. (2013). The ionospheric Bubble Index deduced from magnetic field and plasma observations onboard Swarm. Earth Planets Space 65, 1333–1344. doi: 10.5047/eps.2013.08.005

CrossRef Full Text | Google Scholar

Parrot, M. (1995). Use of satellites to detect seismo-electromagnetic effects. Adv. Space Res. 15, 27–35. doi: 10.1016/0273-1177(95)00072-M

CrossRef Full Text | Google Scholar

Parrot, M. (2011). Statistical analysis of the ion density measured by the satellite DEMETER in relation withthe seismic activity. Earthquake Sci. 24, 513–521. doi: 10.1007/s11589-011-0813-3

CrossRef Full Text | Google Scholar

Perrone, L., De Santis, A., Abbattista, C., Alfonsi, L., Amoruso, L., Carbone, M., et al. (2018). Ionosphericanomalies detected by ionosonde and possibly related to crustal earthquakes in Greece. Ann. Geophys. 36, 361–371. doi: 10.5194/angeo-36-361-2018

CrossRef Full Text | Google Scholar

Pinheiro, K. J., Jackson, A., and Finlay, C. C. (2011). Measurements and uncertainties of the occurrencetime of the 1969, 1978, 1991, and 1999 geomagnetic jerks. Geochem. Geophys. Geosyst. 12:Q10015. doi: 10.1029/2011GC003706

CrossRef Full Text | Google Scholar

Piscini, A., Marchetti, D., and De Santis, A. (2019). Multi-parametric climatological analysis associatedwith global significant volcanic eruptions during 2002-2017. Pure Appl. Geophys. 176, 3629–3647. doi: 10.1007/s00024-019-02147-x

CrossRef Full Text | Google Scholar

Pulinets, S., and Boyarchuk, K. (2005). Ionospheric Precursors of Earthquakes. New York, NY: Springer-Verlag Berlin Heidelberg.

Google Scholar

Pulinets, S., and Ouzounov, D. (2011). Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model –anunified concept for earthquake precursors validation. J. Asian Earth Sci. 41, 371–382. doi: 10.1016/j.jseaes.2010.03.005

CrossRef Full Text | Google Scholar

Rikitake, T., and Yamazaki, Y. (1985). “The nature of resistivity precursor,” in Practical Approaches to Earthquake Prediction and Warning, eds C. Kisslinger and T. Rikitake (Dordrecht: Springer), 559–570.

Google Scholar

Schachtner, R., Pöppel, G., Tomé, A. M., and Lang, E. W. (2009). “Minimum determinant constraint for non-negative matrix factorization,” in Independent Component Analysis and Signal Separation, ICA (Paraty), 106–113.

Google Scholar

Shannon, C. (2001). A mathematical theory of communication. ACM SIGMOBILE Mobile Comput. Commun. Rev. 5, 3–55. doi: 10.1145/584091.584093

CrossRef Full Text | Google Scholar

Shcherbakov, R., Zhuang, J., Zöller, G., and Ogata, Y. (2019). Forecasting the magnitude of the largestexpected earthquake. Nat. Commun. 10:4051. doi: 10.1038/s41467-019-11958-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, X., Zhang, X., Yuan, S., Wang, L., Cao, J., Huang, J., et al. (2018). The state-of-the-art of theChina Seismo-Electromagnetic Satellite mission. Scie. China Technol. Sci. 61, 634–642. doi: 10.1007/s11431-018-9242-0

CrossRef Full Text | Google Scholar

Smaragdis, P. (2004). “Non-negative matrix factor deconvolution; extraction of multiple sound sourcesfrom monophonic inputs,” in International Conference on Independent Component Analysis and Signal Separation (Granada), 494–499.

Google Scholar

Smaragdis, P., and Brown, J. (2003). “Non-negative matrix factorization for polyphonic music transcription,” in IEEE Workshop on Applications of Signal Processing to Audio and Acoustics (New Paltz).

Google Scholar

Wang, H., He, Y., Lühr, H., Kistler, L., Saikin, A., Lund, E., et al. (2019). Storm time EMIC waves observedby Swarm and Van Allen Probe satellites. J. Geophys. Res. Space Phys. 124, 293–312. doi: 10.1029/2018JA026299

CrossRef Full Text | Google Scholar

Wyss, M. (1991). Evaluation of proposed earthquake precursors. Eos Transac. Am. Geophys. Union 72:411. doi: 10.1029/90EO10300

CrossRef Full Text | Google Scholar

Wyss, M. (1997). Second round of earthquake of proposed earthquake precursors. Pure Appl. Geophys. 149, 3–16. doi: 10.1007/BF00945158

CrossRef Full Text | Google Scholar

Wyss, M., and Booth, D. C. (1997). The IASPEI procedure for the evaluation of earthquake precursors. Geophys. J. Int. 131, 423–424. doi: 10.1111/j.1365-246X.1997.tb06587.x

CrossRef Full Text | Google Scholar

Xie, W., Hattori, K., and Han, P. (2019). Temporal variation and statistical assessment of the b value offthe pacific coast of Tokachi, Hokkaido, Japan. Entropy 21:249. doi: 10.3390/e21030249

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, R., Wang, L., Hu, Z., Liu, D., Zhang, X., and Zhang, Y. (2013). Ionospheric disturbances before and afterstrong earthquakes based on DEMETER data. Acta Seismol. Sin. 35, 498–511. doi: 10.1155/2013/530865

CrossRef Full Text

Yizengaw, E., and Groves, K. M. (2018). Longitudinal and seasonal variability of equatorial ionosphericirregularities and electrodynamics. Space Weather 16, 946–968. doi: 10.1029/2018SW001980

CrossRef Full Text | Google Scholar

Zhang, X., Qian, J., Ouyang, X., Shen, X., Cai, J., and Zhao, S. (2009). Ionospheric electromagneticperturbations observed on DEMETER satellite before Chile M7.9 earthquake. Earthquake Sci. 22, 251–255. doi: 10.1007/s11589-009-0251-7

CrossRef Full Text | Google Scholar

Zhang, X., Zhao, S., Song, R., and Zhai, D. (2019). The propagation features of LF radio waves at topsideionosphere and their variations possibly related to Wenchuan earthquake in 2008. Adv. Space Res. 63, 3536–3544. doi: 10.1016/j.asr.2019.02.008

CrossRef Full Text | Google Scholar

Zhima, Z., Cao, J., Liu, W., Fu, H., Wang, T., Zhang, X., et al. (2014). Storm time evolution of ELF/VLFwaves observed by DEMETER satellite. J. Geophys. Res. Space Phys. 119, 2612–2622. doi: 10.1002/2013JA019237

CrossRef Full Text | Google Scholar

Zhima, Z., Hu, Y., Piersanti, M., Shen, X., De Santis, A., Yan, R., et al. (2020). The seismic electromagnetic emissions during the 2010 Mw 7.8 Northern Sumatra earthquake revealed by DEMETER satellite. Front. Earth Sci. 8. doi: 10.3389/feart.2020.572393

CrossRef Full Text | Google Scholar

Zhu, K., Li, K., Fan, M., Chi, C., and Yu, Z. (2019). Precursor analysis associated with the Ecuadorearthquake using swarm a and c satellite magnetic data based on PCA. IEEE Access 7, 93927–93936. doi: 10.1109/ACCESS.2019.2928015