Simultaneous Study of VLF/ULF Anomalies Associated with Earthquakes in Japan

We carried out a simultaneous study of ground-based magnetic field and lower ionospheric anomalies during major earthquakes occurring around Japan in 2010 and 2012. Ultra Low Frequency (ULF) geomagnetic field wave-forms of Esashi station and Very Low Frequency (VLF) Japanese transmitter (JJY) electric signal amplitude received in Moshiri station Hokkaido during nighttime (22:00-02:00 LT) were used to minimize the local interference. Twenty earthquakes having magnitude greater than 5.5 were considered for the data analysis for two years. Nighttime amplitude fluctuations and polarization from the received VLF transmitter signal amplitude and ULF magnetic field respectively were calculated to identify anomalous signatures in relation to every earthquake. We found most earthquakes analyzed indicating VLF amplitude anomalies simultaneously occurred with ULF magnetic field anomalies within a week prior to the earthquakes. Stronger anomalies were observed for larger magnitude and shallower earthquakes. Focal mechanism of earthquakes was also examined to identify the effectiveness of generating anomaly. Both VLF and ULF anomalies were observed for reverse fault type earthquakes occurring under the strong pressure in the crust. Obtained results may indicate the common anomaly source both for VLF and ULF in the lithosphere and are consistent with currently proposed Lithos-phere-Atmosphere-Ionosphere (LAI) coupling scenarios during the earthquake preparation period.


Introduction
Considering the global concern of the effect of earthquake disaster, studies of earthquake prediction are in progress globally using different methods and technologies. Earthquake prediction study is still very complex and needs many serious efforts not only to achieve its high accuracy but also to understand mechanism of precursory signatures. Among various research fields, monitoring electromagnetic emissions and ionospheric perturbations prior to the earthquakes has proved its worth by showing some positive results for the earthquake precursors. Electromagnetic emissions in relation to major seismic activities from DC to HF range have been detected using numerous ground and space-based experiments [1]- [15]. Although individual seismo-electromagnetic parameters have shown their success to identify precursors, a considerable amount of false seismogenic anomalies are presumably due to other geophysical phenomena and local interferences. Therefore monitoring of different electromagnetic parameters has been proposed recently to enhance the success rate of identifying seismogenic anomalies [16] [17] [18].
Among many seismo-electromagnetic parameters, anomalous signatures of ULF magnetic field and VLF transmitter amplitude during the seismic activities have been found much promising for earthquake prediction purpose. Studies of ULF magnetic field anomalies during the seismic activities have been reported by using different methods such as amplitude anomaly method, polarization parameter method, fractal analysis method etc. [5] [19] [20] [21] [22] [23]. These methods are based upon the assumption that the anomalous ULF field is generated by the current source around the earthquake preparation region. In recent years, a new method has been proposed so-called depression method [24]. In this method, amplitudes of ULF waves penetrated from magnetosphere are monitored. Depressions in ULF magnetic field in the frequency range from 0.3 to 0.5 Hz were found before earthquakes presumably due to the attenuation of ULF waves during the propagation through the perturbed ionosphere [25].
Seismo-ionospheric perturbations are monitored by subionospheric VLF transmitter signals. Amplitude reduction method and terminator time method are in use but producing ambiguous results [13] [17] [26]. The reduction method was further refined by trend and nighttime fluctuation methods [27]. In spite of using rather sophisticated data analysis, there are still false seismogenic anomalies because of many perturbation sources of the lower ionosphere other than seismogenic ones such as space weather and atmospheric parameters [28] [29]. Therefore, multi-parameter observations are promising to identify the seismo-electromagnetic signals more accurately rather than observations of single parameter.
In this study we carry out the simultaneous study of anomalies observed by ULF magnetic field and VLF transmitter signals in relation to major earthquakes occurring in Japan in the years of 2010 and 2012. We also try to study the physical mechanism of anomalies for analyzed earthquakes inferred from types of earthquakes determined by their focal mechanism. titude 44.4˚N, longitude 142.3˚E), Hokkaido emitted from JJY transmitter in Fukushima (latitude 37.37˚N, longitude 140.85˚E) to monitor the lower ionospheric conditions. Geomagnetic field data in ULF range from Esashi geomagnetic observatory (latitude 39.23˚N, longitude 141.35˚E) in Iwate are used because Esashi station is situated around the great circle path (GCP) between JJY and Moshiri (JJY-MSR path). Then we study lower ionospheric and geomagnetic anomalies associated with major earthquakes occurred around JJY-MSR path.
VLF subionospheric perturbations have been continuously monitored by the network observation operated by the University of Electro-Communications, Tokyo, Japan. The years of data analysis are 2010 and 2012. We excluded the data of 2011 because of several long data gaps due to the power failure associated with the 2011 earthquake off the Pacific coast of Tohoku. We analyzed the VLF transmitter amplitude data using nighttime fluctuation method [27]. In this method, the trend value of nighttime amplitude was calculated by subtracting the running average of the last 15 days of averaged nighttime amplitude from the daily nighttime average amplitude from 22:00 to 02:00 LT.
The ULF data of Esashi station was obtained from its website www.vldb.gsi.go.jp.
We used the 3-component (H x : North-South, H y : East-West, and H z : Vertical) magnetic field waveform data of 1 Hz sampling rate. The power spectral density for three components (P x , P y , and P z ) were calculated by using Fast Fourier Transform (FFT) between 0.0017 and 0.5 Hz. Furthermore, we calculated polarizations P z /P x and P z /P y using its horizontal and vertical components for different frequency bands. This parameter has been found to increase before earthquakes (started from a few weeks before an earthquake) as reported by [20] [23]. Magnetic Indices data and earthquake data were obtained from www.wdc.kugi.kyoto-u.ac.jp and www.earthquake.usgs.gov respectively. We consider the earthquakes with magnitude greater than 5.5 and depth less than 200 km for our data analysis. Furthermore, we took a spatial limitation by choosing epicenters located within the first Fresnel zone of JJY-Moshiri for VLF data analysis or within the circle of radius 300 km around the Esashi station for ULF data analysis shown in Figure 1 and

Results and Discussion
In 2010, 5 earthquakes occurred within the first Fresnel zone (events 1, 2, 4, 6, and 8) of JJY-MSR VLF transmitter-receiver path. Detailed information of 8 earthquakes analyzed are given in Table 1. In Figure 2 Vertical lines in Figure 2 show the occurrence days of the earthquakes, and the event numbers of the earthquakes are given at the top of each line. Centroid Moment Tensor (CMT) solution of each earthquake is also shown in the third panel. The CMT solution indicates the relation between the six independent parameters of moment tensor representation of earthquake and ground motion  that earthquake generates [30] [31]. This solution provides the information about the motion of plates causing the earthquake, which leads to the type of  earthquakes indicated in Table 1. The trend value of VLF shows a significant decrease within 7 days prior to all 5 earthquakes. In particular, the trend decreases of more than 2σ (i.e. the value-2 in the top panel in Figure 2) before the occurrence of the earthquakes 2 and 4, and the trend decreases around 1.5σ (i.e. −1.5 in the vertical axis) before the earthquakes 1, 6 and 8. It can be seen from Table 1 that magnitudes of the earthquakes 2 and 4 are greater than 6 and depth less than 35 km, while the magnitudes of earthquakes 1, 6 and 8 are less than 6 and depth greater than 35 km. However, there are cases when amplitude anomalies (a significant decrease in trend value) occurred without seismic activities as shown in Figure 2, the reasons of these anomalies may be due to the factors other than seismic activities such as the effect of global geomagnetic disturbance originated from solar activity and atmospheric phenomena toward the lower ionosphere as reported by [28] [29].
To examine the seismo-electromagnetic effects different from VLF anomaly for these earthquakes, ULF magnetic field data observed at Esashi station (located within the first Fresnel zone of the JJY-MSR path) were analyzed using polarization method. We calculated power spectral density (P x , P y , and P z ) and then obtained two polarizations P z /P x and P z /P y for two frequency ranges MP3   within about one week before the earthquake). Earthquake event 4 does not show any clear anomaly despite its proximity to the ULF observation station and relatively large magnitude. The earthquake event 7 has a large polarization but it occurred 8 days before the earthquake, which is the longest lead time among 7 earthquakes with an anomaly. In this figure, epicenters of four earthquakes (events 3, 5, 7 and 9) do not lie in the Fresnel zone of JJY-Moshiri, but the ULF anomaly was identified within the circle of radius of 300 km around Esashi station as shown in Figure 1.
Now we compare the results between VLF and ULF data analysis. All 5 earthquakes within the first Fresnel zone of the VLF path indicate anomalies. Among them, 4 earthquakes show anomalies in the ULF magnetic field. Anomalies both in VLF and ULF occurred within about a week before earthquakes despite lead time varies between earthquakes but occurred the similar timings between VLF and ULF. This simultaneous occurrence both in time and spatial coordinates may indicate that observed anomalies in VLF and ULF are from the same seismogenic origin and also indicate LAI-coupling if the source of ULF field anomalies is originated in the lithosphere. However, many ULF anomalies were identified during 2010 other than earthquake timings and occurred more frequently than those of VLF. ULF magnetic field is perturbed by many external forcing effects including global geomagnetic activity, interferences from human activities in addition to the local seismicity. Furthermore, magnitude of VLF anomaly (trend) has a qualitative relationship with earthquake parameters. For example, two large earthquakes (evets 2 (M = 6.5) and 4 (M = 6.3)) show large anomalies greater than −3σ with relatively small depth, and other events have much smaller anomalies with smaller magnitude and large depth. Whilst, the magnitude of ULF anomaly (polarization) does not clearly show the qualitative relation with earthquake parameters such as magnitude, distance from epicenter. For example, earthquake event 4 is a reverse type (we explain later why this earthquake type is important) with a large magnitude and shallow epifocus. Moreover this event occurred at the shortest distance from the ULF station during moderate ΣKp, so that a large anomaly is expected but the magnetic field does not show any clear anomaly in contrast to VLF. We have not found a reasonable explanation yet and we have to investigate this.
Finally, we focus on the CMT solution for an earthquake, almost all earthquakes analyzed in 2010 were reverse fault type except event 1. The event 1 is a normal type earthquake and occurred one day before the event 2 of a reverse type with much larger magnitude. Therefore, the observed anomaly would be due to the reverse type; however we are not able to clearly identify which earthquake is responsible for the observed anomalies. Nevertheless, both VLF and ULF anomalies occur for reverse type earthquakes for majority of earthquake events.
Same type of data analysis was carried out for the year of 2012. In 2012, 8 earthquakes occurred within the first Fresnel zone of JJY-MSR path (events 1, 2, 6, 7, 8, 9, 10, 11) (Figure 4). And earthquake events 3, 4, and 5 occurred beyond the first Fresnel zone but were located the distance within 300 km from Esashi ULF observation station. Details of the earthquake parameters are summarized in Table 2.
In Figure 5, we present the results of VLF data analysis for 2012. The top panel shows daily variation of nighttime amplitude deviation from running mean normalized by running standard deviations, and ΣKp index is also plotted in the second panel of the figure. In the bottom panel, CMT solution is shown. Earthquakes events 2 shows a significant decrease (i.e. anomaly) in the trend greater than 2σ, while events 6, 8, 9, and 10 show the decrease less than 2σ but greater than 1.5σ. And these anomalies occurred up to 7 days before the earthquakes. The earthquake event 7 does not show any clear anomaly. From Figure  6 we can see that many earthquakes show significant anomaly on ULF magnetic field. The earthquake events 1, 2, 3, 4, 5, 8, and 11 show the enhancement in polarization greater than 2.5σ, while events 6, 7 and 10 do not show clear ULF anomalies.
Comparing the results between VLF and ULF, earthquake event 2 is the 2nd largest magnitude of M = 6.1 and also shallow (D = 15 km) in Table 2    is understandable that a large anomaly was observed both in VLF and ULF because large energy was accumulated close to the ground within the earthquake preparation region. On the other hand, earthquake event 7 does not show significant anomaly both in VLF and ULF. Possible reason for this may be due to the epicenter located under the sea away from the coast and to depth being not shallow (46 km), which make it difficult for anomalous ULF field generated     around epifocus penetrating through the deep water over the epicenter. Although some other earthquakes also occurred underwater, they are very close to coastal area and so an anomaly appeared on the ground without significant attenuation. In the event 6, the anomaly in VLF occurred continuously for a few days around 1 week before the earthquake was identified but not in ULF. This VLF anomaly may be due to global geomagnetic disturbance inferred from rather high ΣKp as previously reported [28] and the local seismo-magnetic effect was masked. Overall, half of the earthquakes have a correlation between VLF and ULF anomalies except high ΣKp period (event 6). Three events (1, 2 and 8) show anomaly and one event (7) does not show anomaly both in VLF and ULF.
Either VLF or ULF anomaly is observed for two earthquakes (10 and 11). Regarding event 10, observed anomaly in VLF can be due to event 9 rather than event 10, because the epicenter of the event 9 is in the VLF path but far from the ULF station (>300 km), which occurred 4 days before event 10 and ULF does not show anomaly for event 10. The reason why no anomaly observed for the event 11 in VLF can be due to the epicenter being located at the edge of the first Fresnel zone, while ULF anomaly was identified because of proximity of the ULF observation station.
We also examined CMT solution for these earthquakes and found these two earthquake events 6 and 9 are normal type earthquakes while other events indicate reverse type(s). The earthquake event 6 occurred during high ΣKp time period and VLF and ULF variations were affected. For the event 9, VLF anomaly was identified probably for the normal type of earthquake. Because of lack of numbers of normal type earthquakes and, we are not able clarify if normal type earthquake is effective to generate anomalies as the reverse type earthquake.
In this study we have made an effort to examine anomalies occurred prior to major earthquakes based upon the ULF geomagnetic field and remotely received VLF transmitter electric amplitude indicating lower ionospheric perturbations.
During the about one week period before earthquakes, both ULF/VLF anomalies were found simultaneously for majority of earthquakes. Earthquakes with larger magnitude and small depth tend to show a stronger anomaly. However, ULF has many anomalies unrelated to earthquakes than VLF during the time period of analysis. These non-seismogenic anomalies are possibly due to other external forcing factors such as both natural and artificial noises affecting the geomagnetic field. Hence anomalies in ULF seem to be less correlated with seismic activities than those in VLF. On the other hand some earthquakes do not show clear anomalies, which may be due to other parameters such as global geomagnetic disturbance causing the anomaly in VLF amplitude but reducing ULF polarization (i.e. horizontal magnetic components tend to increase during high ΣKp). There are other possible reasons why no significant anomalies were observed for some earthquakes during geomagnetically quiet time period. Then we use the focal mechanism of the occurrence of all the earthquakes using CMT solution. Now for explaining the LAI coupling during the seismic activities, some physical scenarios are proposed; first possibility is that the radon gas emanated from the seismic regions propagated upward and created inhomogeneity in the electron density at the base of the ionosphere in the form of localized enhancement.
The observations of intensified radioactive gas emanations before the earthquakes have already been reported e.g. [35]. This mechanism is also applicable for the ionospheric disturbances during the earthquake occurring under the sea.
In addition to the above mentioned mechanism, there is another possibility in which the ionosphere perturbation is caused by seismic activities. This is related to the electric fields generated from the highly pressurized EQ preparation region which may penetrate through the lower ionosphere and modify the ionospheric layer by a few kilometers so that the phase and amplitude of the subionospheric VLF signals are influenced [26]. But this effect may not work effectively for the earthquakes occurred under the deep water.
The other possibility could be the generation of long period gravity waves during the earthquake preparation process and their intensification at the base of the ionosphere. This possibility has been discussed by [15] [16] in great detail.
Above mentioned coupling scenarios are almost verified from our experimental results in which most of the reverse type earthquake are showing anomalies both in VLF and ULF simultaneously, while the normal type may not show clearly as seen from the 5-year statistical study [36]. Furthermore, most VLF anomalies (i.e. lower ionospheric perturbations) are observed almost simultaneously with ULF magnetic anomalies, which indicates rather clear LAI coupling process for seismic preparation process. Finally, multi-parameter approach such as simultaneous observations of VLF amplitude and ULF polarization having taken in this paper contributes not only V. Singh, Y. Hobara to shed light on the generation and propagation mechanism of the precursory anomaly signatures of earthquakes but also enhances the accuracy of precursor identification because these two parameters are physically independent in general, and so spatial extent of anomaly detection and source of interferences (except global geomagnetic activity) are different. Therefore, we are able to validate the anomalies by comparing both VLF and ULF results.

Conclusion
We have presented the simultaneously analyzed results of ULF magnetic field data and VLF subionospheric propagation data for the years of 2010 and 2012 in relation to major earthquakes around Japan. We found that most earthquakes occurring within the first Fresnel zone show a decrease in VLF amplitude (i.e. ionospheric perturbations) within a week before. ULF magnetic field has more non-seismogenic anomalies than VLF. However the majority of the earthquakes located <300 km from the ULF observation station are associated with ULF geomagnetic anomalies as well. Both VLF and ULF anomalies get stronger for larger magnitude and shallower earthquakes. Furthermore, many of these earthquakes with anomalies occurred under strong pressure from reverse fault type earthquakes that were determined by the CMT solution. These results indicate sharing the causative sources generating anomalies in the lithosphere during the earthquake preparation process and are consistent with existing LAI coupling scenarios. Finally, simultaneous observations of ULF and VLF data may enhance the accuracy of precursor identification.