Applied Engineering Using Schumann Resonance for Earthquakes Monitoring

For populations that may be affected, the risks of earthquakes and tsunamis are a major concern worldwide. Therefore, early detection of an event of this type in good time is of the highest priority. The observatories that are capable of detecting Extremely Low Frequency (ELF) waves (<300 Hz) today represent a breakthrough in the early detection and study of such phenomena. In this work, all earthquakes with tsunami associated in history and all existing ELF wave observatories currently located worldwide are represented. It was also noticed how the southern hemisphere lacks coverage in this matter. In this work, the most suitable locations are proposed to cover these geographical areas. Also, ELF data processed obtained from the observatory of the University of Almeria in Calar Alto, Spain are shown. This is a contribution to help protect against natural disasters such as those caused by earthquakes and tsunamis.


Introduction
The global belt of Earth's subduction zones is the region is where most of the most destructive earthquakes and tsunamis occur.The magnitude (M) of an earthquake is the most important factor that determines its violence and destructiveness, but that is not always the case; a relatively low magnitude earthquake may be devastating in a given zone, since construction criteria in the zone must be accounted for, among other factors.In addition, great tsunamis or large earthquakes (M > 8.0) are usually generated whenever extensive areas of the subduction megathrust face a rupture [1].There are examples of earthquake-produced tsunamis in the 2004 Sumatra earthquake (M-9.0)[2], accounting for 220,000 casualties, or in 2006, with the tsunami that developed from the Java earthquake (M-7.8)[3].The intensity and destruction capacity of earthquakes and tsunamis are estimated by considering the area of vertical uplift of the sea bed, and these magnitudes are related to the geometry of the slipping fault when the earthquake rupture is located near the sea floor [4].Seismic events produce variations in the magnetic field as well; the study of electromagnetic phenomena precursory of earthquakes and volcanic eruptions has been a very active subject during the last years [5].Examples of this being the observation of Ultra Low Frequency (ULF) anomalies (300 Hz-3000 Hz) prior to the development of the Loma Prieta earthquake in 1989 near San Francisco [6,7]; later a geomagnetic field variation of 7.2 nT was detected approximately 7 min before the Tohoku earthquake occurred in 2011 and detected by the Iwaki observatory, at 210 km from the epicenter [8][9][10].Another example is the possible detection of a series of pre-earthquakes magnetic anomalies that occurred on 25 April 2015, with epicentre in Nepal and a 7.8 magnitude.Detection was possible through the use of the Swarm magnetic satellites [11].
In addition, after the Fukushima Daiichi nuclear plant disaster and the Great East Japan Earthquake Tsunami, which occurred in 11 March 2011, the requirement of cover human-environment symbiotic points is very important [12].
Tsunamis also have a major impact on economic aspects.After a tsunami occurs the flooding of the surrounding lands by seawater is common, damaging those lands and destroying plot borders.This causes damage to property rights and to the land administration system in general [13].Therefore, after a Tsunami, it is necessary to re-establish the property lines again, using Global Positioning System (GPS) normally, as was the case in Indonesia [14].From a geodesic point of view, it is well known that the definitive displacement of the zones near the epicentre is the results of earthquakes [15], e.g., the postseismic deformation after 2008 Wencham earthquake in China [16], or postseismic changes in Thai geodetic network owing to the Sumatra-Andaman mega-thrust earthquake in 2004 [17] and the Nias earthquake in 2005 [18].
Mega-splay faults are caused when very extensive thrust faults rising from the plate frontier megathrust intersects the seafloor along the lower earring of the margin.It has been hypothesized that these mega-splay faults transfer displacement efficiently to the near surface, contributing to the genesis of seismic phenomena.Recently, these facts have been identified as first order characteristics in the Nankai Trough [19], and are usual in other subduction areas such as Alaska [20], Sunda [21], and Colombia [22].
To understand the development of tsunamis during the formation of large earthquakes, it is necessary to determine the exact location of the slip of the decomposition system, called splay-frontal [23].Instrumental measurements of earthquake and tsunamis largely come from tide seismological and gauge observatories and lately, bottom pressure recorders located in the deep ocean.The first time that the satellites Jason-1 and Topex-Poseidon collected transects of sea-height data using radar altimeter was in the 2004 Indian Ocean tsunami.These data showed a tsunami signal during the across the Indian Ocean [24].The number of publications based on earthquake monitoring increases significantly and becomes a topic of increasing interest in the last decade [25,26].
Some authors suggest that a probabilistic approach might be a method for assessing the risk posed by seismic phenomena for a wide range of magnitudes [27].Due to the relationship between earthquakes and tsunamis, both the recurrence allocation of events in time and the delivery of earthquake or landslide sizes can be used to calculate the tsunami probability.In coastal locations with a broad register of tsunami is the distribution scope that is similar to the ones from other natural damage such as landslides, forest fires, and earthquakes [1], which could be described as a power law [28].But apart from this approach, methods for early detection are needed.
Sensors can collect data to detect seismic phenomena that can lead to major natural disasters by means of a suitable seismic system and a correct treatment of the information.In addition, these systems could help with the early detection of a disaster, supporting the decisions to minimize human and environmental damage [12].In the literature there are several types of sensor systems that measure different environmental responses after an earthquake and could determine the material damage caused.An example is the sensor Structural Health Monitoring (SHM), which measures the acceleration of the response movement of a building after an earthquake [29]; other sensors measure post-seismic deformations in horizontal structures [30] or in vertical structures [31].
Early detection of this type of seismic events can help save many lives.One such method of early detection could be the presence of ELF wave's observatories, because of the relationship between Very Low Frequency/Extremely Low Frequency (VLF/ELF) waves and seismic phenomena [32,33].During the eruption of the in the Kelud volcano eruption, seismic signals showing peaks at 3.7 mHz, 4.8 mHz, 5.7 mHz, and 6.8 mHz were observed in Indonesia.This fact suggests that an ionosphere-atmosphere coupling phenomenon was present along with the lithosphere-atmosphere coupling [34].In addition, mentioning that such electromagnetic phenomena might be used to predict earthquakes, an example of this is shown in the precursors of the Loma-Prieta earthquake [35], or for the Great East Japan earthquake of 2011 (or Tohoku earthquake) [36].
The aim of this paper is to represent the major worldwide earthquakes and tsunamis of twentieth and twenty-first centuries, and to show the location of ELF wave's observatories existing in the world accounting for the relationship between the two phenomena.This allows us to visualize the zones high probability of seismic phenomena.Taking into consideration the presence of seismic phenomena and their possible detection, different areas of the world will be shown which are suited to this kind of observatories.

Seismic Phenomena and ELF Waves
The importance of the study of precursors or indicators of seismic activity is a key to early detecting and preventing succession of major disasters [37].
Seismic events can be related to different types of events.There is evidence indicating that a link between the lithosphere-atmosphere-ionosphere occurs before the occurrence of seismic events [38].Understanding these processes' physics and possibility of it being used to develop an early detection system is the subject of great research and development in recent decades.It has demonstrated the influence of seismic events on this type of ELF signals, setting both types of signals artificially and analyzing the behavior of electromagnetic waves when they are motivated by earthquakes [32].One of the first attempts to standardize ELF wave monitoring on a global scale dates back to 1985, thus it is considered a fairly developed work topic in recent years [39].
Also, there exist other publications that detect anomalies in signals produced by the seismic phenomenon.One example is the correlation between the seismic phenomena and the Outgoing Long-wave Radiation (OLR) data.These signals can be observed by satellites [16], such as the case of the Chinese Earthquake (Wenchuan and Lushan earthquakes) that occurred between 2006 and 2013.
The detection of such anomalies can also be carried out by observatories or ELF/VLF stations [33].These observatories are located in very low-noise locations.In them, the ELF/VLF range (signals up to 300 Hz) and the orthogonal components of the magnetic field (B x , B y and B z ) can be measured.Three orthogonal induction coil magnetometers are normally employed.The B x component usually represents the North-South (NS) geomagnetic component of the magnetic field.It can be determined by the induction magnetometer whose axis is aligned with the direction NS, so that the component B x is sensitive to waves propagating in the East-West (EW) direction.The other component of interest is the B y component that is sensitive to waves propagating in the meridian plane [40].ELF/VLF signals propagate through the environment, and are detected by an appropriately located ELF/VLF receiver over a zone where enhanced radio emissions are happening prior to a great earthquake.Some studies about ionosphere-atmosphere-groundwater phenomena great earthquakes (M > 6.5) associate this occurrence with strong ELF noises or anomalies present in signals located in the ELF frequency band [41].
One of the most important ELF signals, which are affected by the presence of these seismic events is the Schumann Resonance (SR) [42].These resonances are produced by natural phenomena; they were predicted theoretically in the 50's, but no graphical representation of them was established until 1962 [43].The SR signals are very weak signals, on the order of PicoTesla, and occur in the space between the earth and the ionosphere, as indicated in Figure 1.This space acts as a resonant cavity or waveguide for these signals, since it has a finite or limited size [44].In regards to the relationship between the mean of these phenomena ELF and the existence of seismic events, there are relatively recent works showing the appearance of abnormal variations in the amplitude and phase of these ELF signals when crossing regions of certain seismic activity.These changes are due to displacement of the inner boundary of the ionosphere in a few kilometers when seismic activity occurs [48].
Other publications have shown the fact of the signals derived from this phenomenon in a period in which large earthquakes occurred, such as the case developed and understood in Taiwan between 1999 and 2004.It is noted that earthquakes produce a variation of about 1 Hz in the localization of the peak frequency.Another anomaly that is detected in this type of work is an atypical amplitude increase of all its modes, highlighting the variation in the fourth mode of resonance, as observed in large earthquakes M > 6 [49].Other publications indicate that the monitoring this signal allows the checking of the existence of depressions in the frequency of the fourth mode resonance occurring between two and six days before one of the study earthquakes [50].In most of the works, variations or anomalies in these signals can be observed with days or even weeks in advance before the occurrence of any earthquake.For example, variations have been detected three weeks before the earthquake that happened on the 11 of March 2011, with an 8.9 magnitude in front of the Japanese coast [51].Signals anomalies were detected for few days, as well ahead before the Sichuan earthquake that took place with a magnitude of 8.0 [52], and the Kobe earthquake of 6.0 magnitude in 2013 [53].
Therefore, the study and analysis of these signals ELF may be a method for detecting seismic events.The presence of ELF stations or observatories in regions with seismic activity or nearby is useful, so that detection is more effective, despite the resonance phenomenon being a global one.

ELF Sensor
ELF sensors are built by combining vertical electric field sensors and horizontal magnetic induction coils [54].This process is very complex due to the large number of factors involved, such as its low frequency, the weakness of the signal to be detected, the non-uniformity of the ionosphere, the manifestation of natural electromagnetic noise, and the interference caused by industrial zones, appearing as background noise in the process [55].These are some of the reasons behind the small number of ELF sensors compared to other observatories; although that number has increased fourfold since the five observatories accounted for in 1999.
These observatories are located throughout the world.However, the largest concentration of such observatories is in Japan, where the occurrence of seismic events is greater [56].Many of these observatories are attempting to establish a correlation between the ELF signals measured and seismic activity, and to compare with data from other stations.For this, a correlation is established within the The measurement of these signals allows developing various applications.This includes determining the global lightning activity or storms by diurnal variations in the amplitude of the signal continuous records [45], or the study of planetary electromagnetic environment [46].Another application of SR signals that has aroused great interest in recent years is as a global tropical thermometer.Since a frequency shift of resonance peaks has been observed related to the global temperature increase of the planet or global warming [47].
In regards to the relationship between the mean of these phenomena ELF and the existence of seismic events, there are relatively recent works showing the appearance of abnormal variations in the amplitude and phase of these ELF signals when crossing regions of certain seismic activity.These changes are due to displacement of the inner boundary of the ionosphere in a few kilometers when seismic activity occurs [48].
Other publications have shown the fact of the signals derived from this phenomenon in a period in which large earthquakes occurred, such as the case developed and understood in Taiwan between 1999 and 2004.It is noted that earthquakes produce a variation of about 1 Hz in the localization of the peak frequency.Another anomaly that is detected in this type of work is an atypical amplitude increase of all its modes, highlighting the variation in the fourth mode of resonance, as observed in large earthquakes M > 6 [49].Other publications indicate that the monitoring this signal allows the checking of the existence of depressions in the frequency of the fourth mode resonance occurring between two and six days before one of the study earthquakes [50].In most of the works, variations or anomalies in these signals can be observed with days or even weeks in advance before the occurrence of any earthquake.For example, variations have been detected three weeks before the earthquake that happened on the 11 of March 2011, with an 8.9 magnitude in front of the Japanese coast [51].Signals anomalies were detected for few days, as well ahead before the Sichuan earthquake that took place with a magnitude of 8.0 [52], and the Kobe earthquake of 6.0 magnitude in 2013 [53].
Therefore, the study and analysis of these signals ELF may be a method for detecting seismic events.The presence of ELF stations or observatories in regions with seismic activity or nearby is useful, so that detection is more effective, despite the resonance phenomenon being a global one.

ELF Sensor
ELF sensors are built by combining vertical electric field sensors and horizontal magnetic induction coils [54].This process is very complex due to the large number of factors involved, such as its low frequency, the weakness of the signal to be detected, the non-uniformity of the ionosphere, the manifestation of natural electromagnetic noise, and the interference caused by industrial zones, appearing as background noise in the process [55].These are some of the reasons behind the small number of ELF sensors compared to other observatories; although that number has increased fourfold since the five observatories accounted for in 1999.
These observatories are located throughout the world.However, the largest concentration of such observatories is in Japan, where the occurrence of seismic events is greater [56].Many of these observatories are attempting to establish a correlation between the ELF signals measured and seismic activity, and to compare with data from other stations.For this, a correlation is established within the different components of the magnetic field [57], and the development of theoretical models for describing the intensity of the SR by disturbances in the ionosphere over the top of the epicenter of an earthquake.One of these theoretical models is based on ionosphere's vertical conductivity profile model, describing the regular Earth-ionosphere cavity.The possible disturbance is introduced with the modified model.The localized ionosphere modification shows a Gaussian radial dependence; it has a 1-Mm radius, with the decrease reaching 20 km in the lower ionosphere height over the epicenter of the earthquake (Taiwan).Also, the diffraction problem in the Earth-ionosphere cavity with a disturbance can be resolved by using the Stratton-Chu integral equation [58].
One of these ELF stations has been developed and steadily improved by the research group TIC019 of Almeria University and it is located in Calar Alto (Almeria, Spain), Figure 2.This observatory began operation in 2011, being the first observatory of its kind in Spain.
Appl.Sci.2017, 7, 1113 5 of 20 different components of the magnetic field [57], and the development of theoretical models for describing the intensity of the SR by disturbances in the ionosphere over the top of the epicenter of an earthquake.One of these theoretical models is based on ionosphere's vertical conductivity profile model, describing the regular Earth-ionosphere cavity.The possible disturbance is introduced with the modified model.The localized ionosphere modification shows a Gaussian radial dependence; it has a 1-Mm radius, with the decrease reaching 20 km in the lower ionosphere height over the epicenter of the earthquake (Taiwan).Also, the diffraction problem in the Earth-ionosphere cavity with a disturbance can be resolved by using the Stratton-Chu integral equation [58].
One of these ELF stations has been developed and steadily improved by the research group TIC019 of Almeria University and it is located in Calar Alto (Almeria, Spain), Figure 2.This observatory began operation in 2011, being the first observatory of its kind in Spain.The technology used in this type of observatories has evolved significantly in recent years.The general aspects of it are described in a few papers, but only a handful of articles explicitly set the methodology used to capture and measure ELF signals.Figure 3 shows the block diagram of Calar Alto ELF observatory (Almeria, Spain).As the diagram shows, the system consists of a set of stages from signal capture and local real-time processing to transmission to the University of Almeria and the insertion of post-processed data into an interactive database.The stages of the measuring system are as follows: (1) Signal Capture: formed by the ELF magnetic sensors.The technology used in this type of observatories has evolved significantly in recent years.The general aspects of it are described in a few papers, but only a handful of articles explicitly set the methodology used to capture and measure ELF signals.Figure 3 shows the block diagram of Calar Alto ELF observatory (Almeria, Spain).As the diagram shows, the system consists of a set of stages from signal capture and local real-time processing to transmission to the University of Almeria and the insertion of post-processed data into an interactive database.The stages of the measuring system are as follows: (1) Signal Capture: formed by the ELF magnetic sensors.The magnetic sensors are two magnetic antennas oriented along the main cardinal directions (north-south and east-west).Each sensor is made up of two symmetrical coils that provide balanced differential signals helping to reduce common mode effects [59,60].An optional notch filter was implemented as part of the measurement system using a double Active cell [61] to attenuate the power signals.
The A/D converter system consists of a 24-bit high resolution system that is based on oversampling techniques.The obtained data from the converter is subject to a decimating process through temporary averaging.This causes a reduction in the acquisition noise, as well as in the sampling frequency.By a high-performance embedded Digital Signal Processor (DSP) system is done by this pre-processing.Ultimately, data are sent in real-time by a high-speed link to the research group facilities in Almeria University, where is worked using Power Spectral Density (PSD) estimation techniques [62].Through the development of this system, the Schumann resonance signal capture efficiency will be enhanced, thereby facilitating their subsequent analysis and study.The data collected from these observatories are like those presented in Figure 4a,b.In these figures, a periodogram and a spectrogram are represented, respectively, recorded from the Calar Alto ELF station (Spain).These representations were established from the data recordings in normal conditions of the day 16 September 2015.The periodogram (Figure 4a) was obtained by 30 min averaging, the result of Fast Fourier Transform (FFT) of 20 s data segments.The 24 h spectrogram is composed from periodograms that are obtained by averaging each 10-min segments of 20 s.
In this periodogram, the first 7 resonance peaks can be observed, as well as the signal from the power lines.In the spectrogram, you can see beside the Resonance Schumann modes, variations in power line signals and other signals that affect all frequencies about the 500 min, which is due to The magnetic sensors are two magnetic antennas oriented along the main cardinal directions (north-south and east-west).Each sensor is made up of two symmetrical coils that provide balanced differential signals helping to reduce common mode effects [59,60].An optional notch filter was implemented as part of the measurement system using a double Active cell [61] to attenuate the power signals.
The A/D converter system consists of a 24-bit high resolution system that is based on oversampling techniques.The obtained data from the converter is subject to a decimating process through temporary averaging.This causes a reduction in the acquisition noise, as well as in the sampling frequency.By a high-performance embedded Digital Signal Processor (DSP) system is done by this pre-processing.Ultimately, data are sent in real-time by a high-speed link to the research group facilities in Almeria University, where is worked using Power Spectral Density (PSD) estimation techniques [62].Through the development of this system, the Schumann resonance signal capture efficiency will be enhanced, thereby facilitating their subsequent analysis and study.The data collected from these observatories are like those presented in Figure 4a

Worldwide Map of ELF Observatories and Seismic Events
The number of stations or observatories ELF today is not too large, due to many factors, among which is the difficulty of the measurement process or the study of such signals.This study is an issue that is being widely developed in recent years.ELF/VLF stations or observatories exist today can be classified into several types: (a) Observatories that measure the disturbances caused by seismic activity on the ELF radiation, considering subELF observatories.These subELF observatories usually record the possible variations in the Earth's magnetic field, where some have very low frequencies (less than 0.03 Hz).An example of this is the station of Uchimura (Japan), which allows for studying the data collected in the range comprised between 0.01 to 0.033 Hz frequencies [52].In other publications, the signals of lower frequencies (0.01-0.02Hz) are analyzed [54] to establish anomalies or unaccounted perturbations of the earth's magnetic field before the occurrence of several earthquakes of interest [53].Considering the development of such stations of great interest for possible applications in the field of early detection is considered as an emerging theme in the study.In this periodogram, the first 7 resonance peaks can be observed, as well as the signal from the power lines.In the spectrogram, you can see beside the Resonance Schumann modes, variations in power line signals and other signals that affect all frequencies about the 500 min, which is due to strong wind.Therefore, by monitoring the evolution of these electromagnetic signals, we can establish a study of their variation and find out the mechanisms that are related these changes with the occurrence of seismic events.To relate ELF signals' alterations with strong earthquakes, several cases must be studied.The significance of the data would depend on the number of ELF stations around the epicenter, so common ELF alterations could be found and lay out the criteria that might allow early earthquake detection.

Worldwide Map of ELF Observatories and Seismic Events
The number of stations or observatories ELF today is not too large, due to many factors, among which is the difficulty of the measurement process or the study of such signals.This study is an issue that is being widely developed in recent years.ELF/VLF stations or observatories exist today can be classified into several types: (a) Observatories that measure the disturbances caused by seismic activity on the ELF radiation, considering subELF observatories.These subELF observatories usually record the possible variations in the Earth's magnetic field, where some have very low frequencies (less than 0.03 Hz).An example of this is the station of Uchimura (Japan), which allows for studying the data collected in the range comprised between 0.01 to 0.033 Hz frequencies [52].In other publications, the signals of lower frequencies (0.01-0.02Hz) are analyzed [54] to establish anomalies or unaccounted perturbations of the earth's magnetic field before the occurrence of several earthquakes of interest [53].Considering the development of such stations of great interest for possible applications in the field of early detection is considered as an emerging theme in the study.(b) Other observatories had recorded signals in the 1 Hz to 50/60 Hz spectrum.The latter are important frequencies due to the presence of strong interference signals in this band [63].It might be interesting to consider the different sub-bands within this range [64].(c) There are also observatories that monitor the exact value of frequencies.Mainly, these values happen to be multiples of the principal frequency depending on the country.These frequencies are usually 223 Hz (in countries having 60 Hz as principal frequency) and 233 Hz (in countries having 50 Hz as principal frequency) [65].Other observatories use a frequency of 17 Hz [66].
All of the observatories mentioned earlier measure signals of original nature, but there are other stations that monitor artificial signals that are generated by human activity, such as signals in the VLF band that are issued for the navigation aid [67].
Furthermore, not all of the stations have the same measuring devices.The most common are seismometers, acoustic meters, or telluric current meters.These instruments are of great importance to correlate the different parameters of interest with the detection of natural phenomena.The data above can be related to other meteorological data, such as wind speed, humidity, temperature, and atmospheric pressure [48].
Table 1 shows the World-Wide ELF/VLF observatories.This list is in descending order, showing the estimated station range according to the references listed in the table itself.The observatories with range above 1000 km are usually used for communications and shorter ranges are aimed for the detection of natural disasters.The station's sensitivity range is the maximum distance at which it can detect an atmospheric phenomenon.Therefore, this range depends on the magnitude of the different phenomena it can detect.In order to build the worldwide map of seismic events, records were used from the Significant Earthquake Database of National Oceanic and Atmospheric Administration (NOAA), containing information on destructive earthquakes from 2150 B.C. to the present day.All of the featured earthquakes meet at least one of the following: Moderate damage (at least 1 million dollars), 10 or more casualties, Magnitude 7.5 or greater, or the earthquake generated a tsunami [78].More than 1400 earthquakes have been represented in this work, also considering those with tsunami associated.
These particular earthquakes that are accompanied with tsunami occurrence were selected among all others due to the fact that these events have produced huge disasters throughout the history.Thus, their detection is of great interest.The presentation of other earthquakes in the map would have made it much complex and difficult to read due to the high number of seismic events that have occurred in the considered period of time (the last two centuries).They wouldn't add relevant information as well, since the seismic zones are still active.
The previous data was laid out over a map, displaying recent seismic events of interest and the stations indicated in Table 1, with their corresponding detection range (as set by relevant literature).Figure 5 shows the final results.It is worth mentioning that annually, there are about one million earthquakes of M = 2.0 [83].
Appl.Sci.2017, 7, 1113 9 of 20 earthquakes meet at least one of the following: Moderate damage (at least 1 million dollars), 10 or more casualties, Magnitude 7.5 or greater, or the earthquake generated a tsunami [78].More than 1400 earthquakes have been represented in this work, also considering those with tsunami associated.These particular earthquakes that are accompanied with tsunami occurrence were selected among all others due to the fact that these events have produced huge disasters throughout the history.Thus, their detection is of great interest.The presentation of other earthquakes in the map would have made it much complex and difficult to read due to the high number of seismic events that have occurred in the considered period of time (the last two centuries).They wouldn't add relevant information as well, since the seismic zones are still active.
The previous data was laid out over a map, displaying recent seismic events of interest and the stations indicated in Table 1, with their corresponding detection range (as set by relevant literature).Figure 5 shows the final results.It is worth mentioning that annually, there are about one million earthquakes of M = 2.0 [83].The WinkelTripel projection has been used to build the map.The name Tripel (German for "triple") is a reference to Winkel's goal of reducing three kinds of distortion: area, direction, and distance.The WinkelTripel has the smallest skewness, so it is suited for mapping the entire world [84].Nowadays it is found that this projection has been used to show the location of earthquakes recorded all over the earth, mainly due to the good representation of highest latitudes.It is used from The WinkelTripel projection has been used to build the map.The name Tripel (German for "triple") is a reference to Winkel's goal of reducing three kinds of distortion: area, direction, and distance.The WinkelTripel has the smallest skewness, so it is suited for mapping the entire world [84].Nowadays it is found that this projection has been used to show the location of earthquakes recorded all over the earth, mainly due to the good representation of highest latitudes.It is used from 1998 by the National Geographic Society for worldwide maps [85], replacing the Robinson projection as the standard projection.

Results and Discussion
When considering Figure 5, we have developed several maps that are raised or propose new ELF sensors as stations to provide coverage to potentially vulnerable regions.
The label of "high seismic risk" could be tagged to the zones of Central America, Andes, the most oriental part of Europe, Philippines, Indonesia, and Japan, with the three last regions having a high occurrence of earthquakes and tsunami.There are frequent seismic events of great magnitude that involve high damages: the Sumatra Tsunami in 2010 [18] or the Tsunami Japan 2011 [51], and others.This is one of the reasons for which the study and analysis of this type of event is of great importance for the population.Figure 6 shows the current ELF stations along with the ELF stations, as proposed in Table 2.It can be observed that in this figure the coverage of the areas with seismic risk is greater than the one established in Figure 5, where there are only represented the existing ELF stations.
high occurrence of earthquakes and tsunami.There are frequent seismic events of great magnitude that involve high damages: the Sumatra Tsunami in 2010 [18] or the Tsunami Japan 2011 [51], and others.This is one of the reasons for which the study and analysis of this type of event is of great importance for the population.Figure 6 shows the current ELF stations along with the ELF stations, as proposed in Table 2.It can be observed that in this figure the coverage of the areas with seismic risk is greater than the one established in Figure 5, where there are only represented the existing ELF stations.Table 2 proposes with a minimum of 10 new ELF observatories to cover the seismic risk of all these areas of the earth, with a standard range of 1000 km.With this range, all of the zones of interest could have been covered.So, South America can be covered by three observatories that located in Chile and Peru, while Oceania needs five observatories, and the south of Asia needs at least two of them.
If a smaller range is considered, the numbers of suggested stations should be higher, an example shown in Figure 7.In this figure the data in Table 3, where possible ELF stations are established to cover the vulnerable areas seismic coverage with a radius of 500 km are represented.If a smaller range is considered, the numbers of suggested stations should be higher, an example shown in Figure 7.In this figure the data in Table 3, where possible ELF stations are established to cover the vulnerable areas seismic coverage with a radius of 500 km are represented.The minimum number of stations to be considered in this case is 19.So, South America can be covered by six observatories located in Chile, Peru, and Ecuador, while Oceania needs 11 observatories, and the south of Asia needs at least two of them.
There is no theoretical model that establishes a specific range for the detection of this type of phenomena, and appointed error in these measures.The ranges described in Table 1 were estimated empirically according to the sensitivity of the phenomena captured, with usual distances of between 500 km and 1000 km.Consistent with these data, possible ELF stations with a range of 1000 km (Table 2) and 500 km (Table 3) have been proposed.The technology improvements allow for a greater sensitivity for the capture of these phenomena.
Recently, the largest number of the ELF observatories resides in the region of Japan; they cover the major part of its territory.The rest of observatories are distributed on the North American coast and Europe most commonly in its oriental zone.
In the Figure 5, it should be considered that the African continent has less seismic activity that can produce tsunami due to the localization of the tectonic plates, which makes it a stable zone.We can observe that there exist numerous zones in southern hemisphere, where the occurrence of seismic events is elevated.Currently, there are no ELF observatories in these zones, which make the prediction of seismic events possible and thus they are considered as vulnerable areas.
So, in Africa, there is no evidence of significant seismic events recorded and therefore installing ELF observatories would not be a priority, although it would be recommended.However, there are other areas of the world with high earthquake intensity and risk to human life as consequence to the high population density [86].
The low atmospheric attenuation with some latitude dependence of ELF signals (with an average of 0.64 dB/km at 40 Hz) allows for wave propagations at global scale [87].Despite the global character of these signals, the signals anomalies related to the occurrence of seismic events cannot be detected whenever the magnitude is less than 6 (M < 6).It should be considered that seismic phenomena of such magnitudes can also cause disasters and should be taken into account for their early detection [25].This is one of the reasons behind the necessity of establishing new ELF stations in the potentially active regions, such as the Asian zone or the south of America.
The development of new ELF stations in areas of frequent seismic activity could be considered as a very useful tool for the possible detection of seismic phenomena.There are researchers with great experience in the correlation of seismograms with previous anomalies of the natural patterns of ELF signals.The author Hayakawa M. is one of the main advocates of this theory, as indicated by the large number of articles he has developed in this topic, although the list of researchers is rather extensive.Table 4 shows some of the precursor phenomena of seismic events, which have been observed in different earthquakes in recent years.

Precursor Phenomenon Earthquake (Year) Reference
Increase of Very Low Frequency/Extremely Low Frequency (VLF/ELF) electromagnetic noise before and after an earthquake.
Hyogo-Ken Nanbu earthquake (1995) [88] Anomalous increases up to 10 Hz in Ultra Low Frequency (ULF) signals were detected at Shigaraki, 90 km of the epicentre and at Kokubunji, 500 km east of the epicentre.
Hyogoken-Nanbu earthquake (1995) [91] Observation of ULF anomalies prior to the development of an earthquake.Loma Prieta earthquake (1989) [6,7] A geomagnetic field variation of 7.2 nT was detected approximately 7 min before an earthquake.Tohoku earthquake (2011) [8][9][10] Detection of a series of pre-earthquakes magnetic anomalies.Nepal earthquake (2015) [11] Abnormal variations in the amplitude and phase of these ELF signals when crossing regions of certain seismic activity.
earthquakes in Taiwan (1999-2004) [48] An atypical amplitude increase of all its modes, highlighting the variation in the fourth mode of resonance.Chi-chi earthquake (1999) [49] The existence of depressions in the frequency of the fourth mode resonance occurring between 2 and 6 days before one of earthquakes.
The Kamchatka region (during the last 30 years) [50] ELF Signals anomalies before an earthquake.Sichuan earthquake (2008) [52] ELF Signals anomalies before an earthquake.Kobe earthquake (2013) [53] Interference Signals in the ELF band.Japan Earthquake (2011) [64] The table above shows that the number of earthquake precursor electromagnetic phenomena is high, and that the majority are due to anomalies in the signals that were captured in different ELF stations.Such anomalies comprise from the increase or decrease of the usual intensity of the signal, a shift in the frequency of some of its modes of resonance or the existence of signals considered as ELF noise, signals not usually be in this band.
It should be also considered that there is no exact methodology for the temporary detection of these precursory phenomena.Some have been observed for several minutes before an earthquake, others a few days before or even weeks.Therefore, it is an empirical methodology that depends on several factors, such as the magnitude of the seismic event, the proximity to the ELF station, its instrumentation, and even the existing environment.In spite of this, there is evidence that a continuous study of the patterns of ELF band signals (such as RS signal) can be a powerful methodology that allows for the possible early detection of important seismic phenomena.
A seismograph is capable of collecting seismograms from any earthquake (M ≈ 7) at any point on the earth.The filtering effect of the high seismic frequencies allows us to determine the remoteness of the earthquake and to center our study of correlation with the earthquakes located within a radius of between 500 km and 1000 km, if it is located within range of a the ELF station.These are the estimated distances for the detection and correlation of the previous phenomena, as well as a possible explanation for the ELF stations performance ranges that are proposed in this paper.
Correlation studies of the ELF signals captured at the Calar Alto station with data from the Alboran earthquake (M = 6.3) [92], which occurred in January 2016, are being carried out.The results of this study may be adequate to extract quantitative conclusions of this methodology and to further strengthen its potential for the early detection of possible seismic events.

Conclusions
This paper shows the existing ELF stations and their experimental range of coverage, for the detection of seismic events.On the other hand, the earthquakes with related tsunamis that occurred throughout history were located.When considering the relationship between the detection of precursory electromagnetic phenomena and seismic events, this work highlights the lack of coverage for almost the entire southern hemisphere, warning the community about the seismic risk for the area of South America and South Asia for their high seismic activity.The location of ELF observatories is proposed to provide a minimum coverage of seismic risk alert worldwide, although this methodology is largely experimental and widely studied today.
The ELF stations proposed in this paper are set in the areas of the greatest potential risk in the development of high magnitude seismic events or damage to populated areas, depending on the experimental distances for detecting this kind of phenomena (500 km and 1000 km).The proposed increase in the number of stations allows for greater potential comparative study and better behavior of these events.With this procedure, it should have carried out a greater number of works and the distribution scope is similar to the ones from other natural damage such as landslides, forest fires, and earthquakes studies to confirm and continue with the work that reinforces the use of these observing systems as a tool for early detection of seismic events of interest.So, the goal is to avoid possible global catastrophes with the related study of ELF signals of the observatories.
Despite that the relation between the detection of these seismic events and the ELF signals measurements (highlighting the Schumann resonance signal) is not totally verified; many studies are being contributed into investigation in the last few years.Works that showed anomalies or variations in these ELF signals a few weeks or days before the occurrence of large seismic events occur.These investigations are of interest due the major implications that can soon be developed through them soon.Therefore, this work opens new perspectives in the early detection of earthquakes and tsunamis in the world.

( 2 )
Stage of conditioning: of signal, constituted by a high gain differential amplifier and a level scaled adapter.Through this process a suitable signal is obtained at the input of the next stage, with a maximum dynamic range input and the lowest output clipping.A passband filter ranging from 3 to 100 Hz is also included in this stage.This filter attenuates signals outside the desired band, preventing undesired intermodulation or cross-modulation phenomena.(3) The stage of the Analog/Digital converter (A/D converter): This stage has two channels, one for each sensor.It converts the analog signal into 24-bit samples per channel with a sample rate of 196 samples/s.(4) The last stage is the data transmission and storage.It is constituted by a local information storage

Figure 5 .
Figure 5. Tsunami risk alerts worldwide through ELF Sensors.

Figure 5 .
Figure 5. Tsunami risk alerts worldwide through ELF Sensors.

Figure 6 .
Figure 6.Tsunami risk alerts worldwide through proposed ELF observatories with a 1000 km of ratio.

Figure 6 .
Figure 6.Tsunami risk alerts worldwide through proposed ELF observatories with a 1000 km of ratio.

Figure 7 .
Figure 7. Tsunami risk alerts worldwide through proposed ELF observatories with a 500 km of radio.

Figure 7 .
Figure 7. Tsunami risk alerts worldwide through proposed ELF observatories with a 500 km of radio.