1. Introduction
There are many countries or regions suffered by strong earthquakes in the world as Japan, China, Indonesia, Turkey, Italia, USA, Central America, South America, and so on. One of the most efficient disaster mitigations has been assumed to predict them and take appropriate actions beforehand to save human lives before occurrence of the events. We have only one successful example of earthquake prediction, the Haicheng Earthquake (1975 Feb. 04) in China. Many countries, especially earthquake-prone countries have been involved to find practical methods more than half of a century. But there are still no practical methods despite great efforts. There spreads general pessimistic atmosphere against positive stance to the earthquake prediction with the result of decrease of number of scientists and engineers interested in the investigation. Notwithstanding, there remain still not a small number of experts engaged in the earthquake prediction and try to find a new breakthrough.
We have been investigating underground electric field anomalies since 1989 using special antenna. Analyses of the large amount of data at different geological sites and for crustal activities from 1989 to 2002 have revealed that there are plausible anomalies in both ULF and VLF bands. In addition, we found the electric signals propagated some ten times faster than the seismic P wave, i.e. 100 km/s, suggesting the possibility of much earlier alarms compared to EEW based on seismic waves. The most interesting and profitable evidence for the prediction was detected at the time of the great Tohoku earthquake in 2011. The project was planned to reveal the source of the high-speed phenomena propagating such a high speed several times of seismic P-wave by using measurement instrument of wider frequency band. The special pulse-like variations in the ELF band were found to occur only in the nucleation stage just a week before the earthquake, suggesting that the practical method of imminent prediction is possible. Analyses suggested that the events were caused by micro-cracks occurring in the nucleation stage, having features consistent with recent theoretical and experimental results in seismo-electromagnetism. Results of research in seismo-electromagnetics have a wide scope of applications in variety field of technology.
2. Method
2.1. Sensor
We use a special sensor, made of a long vertical steel casing pipe that is installed into a borehole 200 m to 1800 m deep. The antenna invented by [1] is of the mono-pole type and is immersed underground (Figure 1) in the conductive crust analogous to the towing tail cable antenna used for submarines underwater communication. The steel casing pipe functions as a monopole of the antenna in the conductive substance and the copper wire buried underground, surrounding the casing pipe, as the reflector [2] [3] [4] [5]. The boreholes were originally for different purposes such as for the measurements of seismic motion, ground water level and chemical contents, ground deformation, pumping of natural gas and belonged to several organizations (Fujinawa et al. [6] [7]). The recording frequencies are DC band (0 - 0.7 Hz), ULF (0.001 - 0.7 Hz) and ELF/VLF (1 kHz - 9 kHz). The electromagnetic field observation stations have been constructed in central Japan since 1988, and thirteen observation sites were used (Figure 2).
2.2. Observation
Prior observation confirmed that the system has high robustness against both meteorological and urban noise, enabling us to find signals related to crustal activities [4] [5] [6] [7]. Exceptions are the large noises induced by nearby
Figure 1. Schematic diagram of multiple component electromagnetic field observations. The electric field is measured by a pair of horizontal dipoles of 8 m in length with a span of 10 ~ 50 m and using a vertical mono-pole of 100 - 1200 m length surrounded by wire with a radius about 5 - 50 m buried at a depth of about l m.
Figure 2. Locations of the electro-magnetic field observation sites in central Japan utilizing borehole antenna 150 - 1200 m long as the principal sensor at eleven sites: Hasaki (HAS), Katsuura (KTU), Tsukuba (TKB), Chikura (CKR), Awano (AWN), Izu-Oshima (OHS), Miyakejima (MIY), Nagaoka (NGO), Koufu (KOF), Sagara (SAG) and Hodaka (HOD).
lightning discharges, which are easily discriminated because of large difference of strength and waveform. After confirming the performance of the borehole antenna our observation was conducted only with the borehole antenna without magnetic sensors and electric dipoles (Figure 1).
The observation network was built and has operated since 1989, and the sites are situated at several kinds of geological conditions, which are nearly uniform in central Japan (Figure 2). Several kinds of anomalies were picked up and investigated from the points of view of finding precursors related to earthquakes and volcanic eruptions. The ULF band anomalies [3] have been most extensively investigated because they were detected in almost all cases of nearby volcanic eruptions and seismic swarms [5] [6] [7].
2.3. Temporal Experiment
Sometimes temporal experiments were focused on topics. For example, we temporally made high-speed recordings with a time resolution of 1 sec for about 27 months (from June 1994 to December 2002) to investigate the characteristics of precursory anomalous signals. These experiments suggested that there are electromagnetic variations starting from the origin times of moderate earthquakes. Moreover, the evolution of signals correlated with seismic motion at the site, some 100 km distance from the epicenter.
We could not analyze higher-frequency phenomena in detail because of the system’s slow recording speed. Our new observation, which started on March 3, 2011 at Hasaki (now Kamisu City, Ibaraki Prefecture) enabled us to investigate particular co-rupture signals using an updated dynamic range. We could detect very interesting electric variations, suggested to be induced by microcracks in the nucleation stage of earthquakes [8].
3. ULF Band Anomalies
3.1. Observational Result
Electromagnetic approaches have been providing many interesting results, shedding light on the preparatory process for earthquakes [9] [10] [11]. There are many observations of electromagnetic anomalies to be asserted as precursory phenomena. However, there are few well-organized observations that confirm intrinsic relationships between the crustal fracture process and the candidates of precursors. We need a framework of investigation containing generally approved evaluation criteria [12]. Here, we present the electric field observation from a sensitive special sensor referring to the guidelines as far as possible.
Moreover, there are many generation mechanisms that induce electromagnetic anomalies. The fluid-phase has been noted to understand the complex physio-chemical process of earthquakes. Fluids are suggested to play important roles in long-term structural evolution, fault creeping, and fault ruptures (e.g., [13] [14] [15]. Ground water levels and chemical content have been observed in quantities that suggest the existence of changes in the ground water regime before an earthquake (e.g., [9] [16] [17] [18]). A new generation model was presented to explain the characteristic features of the GUV [5] [6].
3.2. Characteristics of Anomalies
The ULF band characteristic waveform was first identified at the time of volcanic eruption activities in 1990 at Izu-Oshima [3]. The amplitude increased suddenly before gradual decreasing to about half the peak amplitude, then decreasing sharply to one-tenth of the peak, and finally returning to the gradual trend level (Figure 3). The signals have never been detected in a normal state, but have been observed in almost all volcanic eruption activities and seismic swarms occurring near the observation sites: the earthquake swarm around Mt. Hodaka in 1998 [7] and Mt. Nikko-Shirane in 1999, and the volcanic eruption at the Miyake Island in 2000. The phenomena have been suggested to be induced by rapid ground motions through the electro-kinetic effect [19] [20] [21] [22] [23]. The ULF variation is called: “geyser like ULF signals “(GUV) because of the similarity of the temporal change of the electric strength with the water height change of the geyser. We thought that the whole of the pulse-like variations detected by the slow speed recordings are of this type.
3.3. Generation Mechanism
Electromagnetic field variations, as self-potential change, have been attributed to piezo-effects, MHD effects, and electro-kinetic effects [19] [20]. Quantitative comparison of these effects suggests that the electro-kinetic effects are about one-order larger than other effects [24]. [22] first attempted to explain the magnetic anomalies associated with the Matsushiro seismic swarm based on electro-kinetic effects. In addition, a self-potential (SP) anomaly has been interpreted on basis of the streaming potential model, based on the Onsager relation [17].
Previous descriptions of the electromagnetic field, due to the electro-kinetic effect, have been limited to a stationary case, and the source region is assumed to be at the boundaries of different substances (e. g. [23]). An extended estimation
Figure 3. GUV type anomaly. a) the observed waveform. b) simulated waveform assuming the source function.
of the induced field in the case of transient generating forces was first attempted by [25] using the Maxwell equation. The transient electric field can be estimated if the pressure distributions are given in a specified conductivity distribution. The time evolution of the electrical field is highly similar to that at the zone near the source, because of the diffusion property i.e. the phase 1 and phase 2 of a sharp change are conserved in the near field. However, at several kilometers distance from the source, the signal is sufficiently diffused to lose the original sharpness at the driving point. The property could be also used to infer the source distance.
The analysis of GUV indicates the following points,
1) The electric field anomaly is well correlated with the seismic swarm and volcanic eruptions around the sites.
2) The electric anomaly is suggested to be induced by variations in the confined water circulation regime through electric-kinetic effects.
3) The pore water pressure anomalies, preceding rapid crustal deformation from earthquakes, are suggested to be detected by the borehole antenna.
4. Co-Seismic EM Field
4.1. Co-Rupture Electric Variation
There are many reports indicating electromagnetic signals are associated with seismic waves, but there is little field observation of waves which propagate faster than P-wave. The first observation was reported using the magnetic measurement [26]. The electromagnetic wave before P wave arrival, but not starting from the origin time, has been reported frequently [27] [28] [29]. A model based on the piezoelectric effect and the dislocation theory of fault was developed by [27] to explain co-seismic electric signals.
Field detection of the co-rupture electromagnetic anomalies is very rare. One was detected at the time of the great Hokkaido Earthquake in 1994, and the other was the electric field detection starting from the origin time of two moderate earthquakes [30] as presented in this paper.
4.2. Observed Result
The anomalies B (Figure 4) observed at Hasaki associated with the major earthquake 4.9 on October 16, 2002 are quite different from the well-known co-seismic electric variation: the co-seismic electric variation has good correlation with the seismic waveform. At this temporary observation we made a continuous high-speed recording of time resolving power of 0.1 s using the radio time signal adjustment. The recording is shown in the upper part of Figure 4(a). The initial phase of B is clearly seen from 42 m 50 s (T0 in the Fig.) to 43 m 5 s (P) in the record. T0 is the origin time of the earthquake, determined by the National Research Institute for Earth Science and Disaster Prevention (NIED) [31]. The lower part of Figure 4(b) is the seismic record at Choshi of the strong motion seismograph network (K-NET).
Figure 4. Electric field evolved in good correlation with seismic wave propagation. 1) A high-speed monitoring record is shown. In the recording the time resolution is 1 sec and the clock is calibrated to of 1 ms. T0 shows the origin time of earthquake, and P, S are P, S waves arrival, respectively at Choshi. The starting time of anomaly is 11 h 42 m 49.0 s (±1 s), very close to the origin time of the earthquake, T0 = 11 h 42 m 50.0 s. 2) The strong motion record of the earthquake at Choshi near Hasaki (Figure 2). The time axis of the strong motion record is shifted corresponding to the position difference between the electric field observation site HASAKI and the K-net site Choshi.
Comparison of the electric field variation and the observed seismic waveform leads us to the following points,
1) the anomaly started from the origin time of the earthquake,
2) the electric field strength increased gradually until P wave arrival,
3) the strength decreased after P-wave arrival,
4) the strength increased sharply from S-wave arrival, and decreased slowly as the seismic wave diminished,
5) the propagation speed was estimated to be larger than 100 km/s.
Similar high-speed recordings were taken by Dr. Takahashi for about 27 months, from June 1994 to December 2002: 20 months at Hasaki and 7 months at Chikura. The records were checked by eye to find similar co-seismic signals induced by moderately large earthquakes. Seismic events with a magnitude M ≥ 4, and epicentral distance less than 100 km are selected by use of the official seismic catalog of Japan Meteorological Agency (JMA). Total of 55 earthquakes are picked up. Careful analysis found that there were three cases of co-rupture anomalies in 42 events at Hasaki, and 2 cases in 24 events at Chikura. Those signals were more or less similar, i.e. the biggest signal at S-wave arrival, but the correlation with the seismic wave form is not so clear compared with that of October 16, 2002. Anyhow it is important to know there are another 5 examples of co-seismic electric signals travel with speed larger than 100 km/s.
4.3. Mechanism
The strong correlation of the electric field variation with the P-wave leads us to infer that the generation zone is limited to the area around the seismic wave front. [24] indicated that Seismic Electric Signals (SES) are generated, first by active sources such as crustal deformations, crack generation, and pour water movement, second by passive sources such as changes in the physical parameters, i.e. electrical conductivity, magnetic susceptibility and electric permittivity, and third by apparent sources i.e. observation and environment noises. The active sources are of piezo-magnetic, piezo-electric, non-classical piezo-electric, electro-kinetic, and induction effects. The present phenomena do not seem to have been induced by the passive causes because of the large difference of the time constants between the two phenomena. External causes are also discarded because of the high S/N ratio of the sensor.
Order estimation of the SES strength [32] suggests that the dominant causes are the seismo-dynamic effect [28] and/or the electro-kinetic effect. The unstable detection of the present phenomena excludes the induction effect. We suspect that the electro-kinetic effect induces an increase in the confined pressure at the wave front [30] [32]. The SES for the case of infinite half space was estimated to indicate that the far field changes as r−2 (r: focal distance) [32] and the signal reaches the point up to more than 160 km from the wave front. These results agree with the present observations.
The phenomenon suggests a new measurement method in seismology using electromagnetic sensors. We can approximate that the first detection time is the origin time, with the result to reduce one of the unknowns in the focal parameters determination. One of most interesting application is for earthquake early warning (EEW). The most developed EEW has been working since 2007 [33] [34] [35] [36] [37]. This system uses the focal parameters determined by P-wave arrival data at a small number of observation sites. The practical use of this information started from October 2007 in Japan with the implementation of counter-measure systems, such as the emergency shutdown of gas fittings and home electrical devices (28). But it generally provides insufficient time (seconds) in cases of shallow depth just under the city. A new electromagnetic sensor will afford additional seconds, along with an accompanying decrease in the extent of earthquake disaster [36] [37].
5. Microcrack
The identification of anomalous phenomena is a prerequisite for understanding the physical processes to be applied to earthquake forecasting (e.g., [16] [38]). For short-term prediction, numerous observations by seismic, geodetic, hydrologic and electromagnetic approaches have been conducted to find precursory phenomena. Of them, pre-shocks have been the most direct and plausible phenomenon for identifying the earthquake nucleation stage. Previous investigations for pre-shock activity of natural earthquakes or the micro-cracks activity of laboratory experiments (e.g., [39] [40] [41] [42] [43]) have provided essential evidence for the building of the physical and statistical models of the nucleation process.
As for natural earthquakes, however, the evidence is limited to an average image over many samples. There is considerable variability in the results with few plausible clues to be applied for practical prediction (e.g., [38]). On the other hand, with the rock fracture experiment we can always observe the characteristic rapid increase of acoustic emission activity in the nucleation period [42] [43], especially in the accelerated increase of micro-cracks known to follow the modified Ohmori formula. However, there are only a few pioneering field studies (e.g., [44] [45] [46]). Those researchers found that the acoustic emission intensity is high during a period spanning about a half-day before and after an earthquake, without means to define the nucleation phase.
5.1. Results
An updated data logger was installed at Hasaki on March 3, 2011, 8 days before the Tohoku earthquake. We checked the AC and DC daily records obtained for March 3 to March 31, 2011. A daily record shows nothing except for background noises, seen as a belt on the calmest state and as a short time-period pulse-like activity on active days. All the pulse-like variations were picked up. There is a total of 224 events, including the large event of the co-seismic variation corresponding to the main shock, being consistent with one of the guidelines to approve the candidate of precursor [12].
Whole events are classified into three groups, depending on their signal duration: type A (Figure 5(a)), type B (Figures 5(b)-(d)), and type C (Figure 5(e)). The type A event (Figure 5(a)) is simply a so-called ULF geyser-like variation (GUV). The type B and C events have a much smaller duration, from 10 ms to several tens of ms for type B, and 0.3 - 0.5 ms for type C. Type B pulses are subdivided into three types, depending on their waveform (Figures 5(b)-(d)). Type B-1 (Figure 5(b)) is similar to type A in its abrupt rise and maintaining polarity, similar to the relaxation waveform, but is different in decaying stepwise. Similar electric variation is reported to be associated with acoustic emissions in rock deformation experiments, in a hydraulic fracturing test, and in a simulation to investigate the Loma Prieta Earthquake. Type B-2 (Figure 5(c)) is a wave packet similar to seismic waveform. Type C (Figure 5(e)) is a wave packet with a carrier frequency of some 5 kHz. The crack theory [38] suggests that type A and type B1 correspond to tensile cracks, types B-2 and C to shear modes, and type B-3 to the compounded mode. Figure 6 shows the time evolution of the number of the pulses of types A (green), B (red), C (blue). Just before the great Tohoku Earthquake (“EQ3.11” in the Figure), the type B variations were dominant, and type A rare. After the earthquake, types A and C appeared without type B. The three kinds of B band anomalies evolved similarly, but the wave packet type (B-2) tended to be dominant approaching the main shock, in agreement with the rock fracture experiment results.
The time evolution of the type B phenomena had a strong similarity to that of the acoustic emissions just before the rupture in the rock experiment, and to that of the pre-shocks before the main shock. Fitting the inverse power law (tf - t)−n to the rate of occurrence number N(t) of cracks at arbitrary time t (day starting March 1, 2011) from the final rupture time tf results in the index number n = 2.0 and tf = 10.5. The index number 2.0 is nearly the same as 1.7 for the foreshock,
Figure 5. Waveform of the pulse-like signals.
Figure 6. The Tohoku Earthquake, magnitude 9.0 occurred on March 11, 2011. (a) Time evolution pulse-like signal of Type A (green). Almost all of the type A variation occurred after the main shock. (b) Time evolution of Type B (red) including B-1, B-2, and B-3. The type B variation is suggested to be particularly useful for imminent prediction. The continuous curve shows the cumulative number of the whole B-type variations. (c) Time evolution of Type C (blue). This type occurred only after the after-effect stage of the seismic cycle. The Tohoku Earthquake, magnitude 9.0 occurred on March 11, 2011. There are several time periods (grey horizontal line) without data due to logger problems and malfunction caused by the severe main shock.
but about 4 times larger than n = 0.5 for the acoustic emission in the laboratory experiment. The assumed time 10.5 agrees satisfactorily with the actual value of 11.6. Enough of many microcracks enabled us to know the time evolution of activity and to suggest the final rupture time.
5.2. Generation Mechanism (MC)
The study of seismo-electric variation had a new era after the work of [47]. The physics of wave generation and propagation in porous media, saturated with fluid, was formulated from the microscopic point of view considering the individual motion of main rock substance and the fluid phase. He derived equations for the coupled poro-elastic and EM field on the basis of Biot’s poro-elastodynamic equations and from Maxwell’s equations of electrodynamics. The wave mode equation for a fixed frequency has been shown to have four modes, fast and slow compressional P-wave (Pf and Ps) and two transversal modes (S –wave and electromagnetic modes: S, EM). The former three modes are co-seismic electric field changes accompanying seismic waves. The EM mode is an electromagnetic wave in the fully saturated porous media with dynamic electric permittivity and dynamic conductivity, depending on the frequency. The analytical form was given by [47] [48]. Those seismo-electric waves have been confirmed in the numerical and laboratory experiments [49] - [55] indicating three phases Pf, S, EM modes are present though the EM mode is extremely small to be observed in the field [49] [50] [51].
The electrokinetic coupling model was analyzed by [49] for the double couple sources providing a useful means for multi-disciplinary investigation of earthquake precursor through electromagnetic and seismological approaches. Here the observed anomalous electric field variations are analyzed based on recent results on the investigation of seismo-electromagnetism started by [48].
5.3. Electric Strength of Each Phases
The electric field variations associated with microcracks sometimes have a clear P-phase besides a dominant S-phase and some events have a small phase, which was assumed to be EM mode. We made quantitative and qualitative comparisons of the observed and estimated amplitude of each mode using electric field description induced by double couple point source in the isotropic medium [48] [49]. We used the basic formulation of [49] and the subsequent formulation for the double couple source of [49] to compare the observed amplitude with those calculated. [49] showed that the electric field E, induced by double couple of the dipole moment M0, is described by the far field, intermediate, and near field parts as per the seismic waves in the elastic medium. In the far field approximation, the electric field Epf of fast P-wave (Pf) and Es of S-wave(S) field are compared [55] for events of shear mode. We used only shear mode events.
The ratio Ap/As is about 30 times: the strength of the Pf-wave is much larger than the S phase. On the other hand, the observed ratio shows that the strength of the EM mode to that of the S mode is 0.030 ± 0.02. The value depends on the wave frequency parameters of the medium. In the case of frequency 300 Hz, the ratio remains almost the same for Apf/As, and one third for Aem/As, respectively, compared with those of 500 Hz. The permeability difference of 2 order (ko = 10−14 m2, ko = 10−12 m2) results in almost the same values for Apf/As, and about a 10% decrease for Eem/Es.
We can see that there is considerable agreement between the observed and estimated electric strength. At the same time, it is noted that S and EM modes have observed values an order larger than those by estimation; the Pf mode has about 1/50 of the estimated value. The exceptionally large estimated values are discussed in the following.
The electric strength of fast-P, S, and EM modes of the shear type events are compared. The averages of observed amplitude orders as As > Apf > Aem, are partially different from the present and previous numerical simulation results. Analytical estimation is conducted on the basis of the seismo-electromagnetic theory of [47], assuming the isotropic homogeneous medium is filled with confined solvent. The contradictory relation between the observed and numerical simulation results may be attributed to the assumed whole space model rather than the half space model. Overall agreement between the observed and estimated results support that, 1) Pulse-like events of B-type of electrogram before the Tohoku Earthquake in 2011 were induced by the rupture of microcracks in the nucleation period of the main shock, 2) Observed phases of waveform, P, S, and EM are confirmed to be co-seismic Pf, S and free electromagnetic wave assumed by [47], 3) the SES phenomena in the field can be analyzed from the formation of [47] and from successive progress by successors (e.g., Revil).
6. Conclusions
We have been observing electric field changes at 13 ground electric field monitoring stations in the central parts of Japan, in relation to the prediction on earthquakes and volcanic eruptions since 1989. The detection sensor is a special antenna made of a vertical pipe casing placed in boreholes. The system is very sensitive to crustal activities, such as earthquakes and volcanic eruptions, and is highly robust to both meteorological and urban noise. Three stages of observation have led to the following results,
1) GUV (Geyser type ULF Variation)
a) The electric field anomaly is well correlated with seismic swarm or volcanic eruptions around sites.
b) It is suggested to be induced by changes of the confined water circulation through electric-kinetic effects.
c) The GUV can be profitable signals to detect precursors of earthquakes and volcanic eruptions.
2) Co-rupture Signal
We presented the observational data showing that there is electromagnetic field variation at the seismic wave front with the effective propagation speed of order 100 km/s. This phenomenon is suggested to be more useful for EEW technology, producing more lead-time before arrival of the “big shock”, when compared to the usual monitoring that relies on elastic wave detection. Further studies are needed before we can develop instruments able to detect electric field variations. Future efforts may need to concentrate on increasing the dynamic range of the measurements as well as gaining understanding of the generation mechanism of the phenomenon.
3) Micro-crack
We have observed that the electromagnetic method can detect microcracks preceding natural earthquakes by use of a special underground antenna. In addition, the detection distance for microcracks is suspected to be much larger than by the acoustic emission method. The long distance detection of microcracks, in comparison to acoustic emissions, is caused by the small decay rate of the seismo-electric signals trapped on the conductivity contrast. The nucleation stage of a crustal rupture can be identified through monitoring of the electric field variations, focusing on particular waveforms in the selected frequency ranges.
These pulse-like phenomena are suggested to be one of the most probable candidates of precursors for an imminent prediction, satisfying many of the conditions required for the precursor by the IAESPI committee. Our present findings could provide a breakthrough for the imminent prediction of earthquakes. This method can also be applied to volcanic eruptions and large-scale ground slides.
4) Remarks
Precursory phenomena are investigated by means of seismic activity, crustal deformation, groundwater anomalies, and electromagnetic anomalies. Present findings of the overall agreement of the observation and estimation for the seismo-elecromagnetic phenomena suggest that seismological, crustal, groundwater, and electromagnetic phenomena that appear in the preparatory process of earthquake occurrence can be analyzed and discussed on the theory of behavior of the seismic activity, hydro dynamical and electromagnetic, and phenomena induced by variety of cracks taking account of the electro kinetic effect. Many of the observational electric anomalies can be more satisfactory interpreted based on the unified scheme. A multidisciplinary approach, based on the unified theory induced by varieties of crack, is expected to open a new window to building up practical earthquake prediction.