Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 44774.
Published online 2017 March 21. doi:  10.1038/srep44774
PMCID: PMC5359576

Causal mechanisms of seismo-EM phenomena during the 1965–1967 Matsushiro earthquake swarm


The 1965–1967 Matsushiro earthquake swarm in central Japan exhibited two unique characteristics. The first was a hydro-mechanical crust rupture resulting from degassing, volume expansion of CO2/water, and a crack opening within the critically stressed crust under a strike-slip stress. The other was, despite the lower total seismic energy, the occurrence of complexed seismo-electromagnetic (seismo-EM) phenomena of the geomagnetic intensity increase, unusual earthquake lights (EQLs) and atmospheric electric field (AEF) variations. Although the basic rupture process of this swarm of earthquakes is reasonably understood in terms of hydro-mechanical crust rupture, the associated seismo-EM processes remain largely unexplained. Here, we describe a series of seismo-EM mechanisms involved in the hydro-mechanical rupture process, as observed by coupling the electric interaction of rock rupture with CO2 gas and the dielectric-barrier discharge of the modelled fields in laboratory experiments. We found that CO2 gases passing through the newly created fracture surface of the rock were electrified to generate pressure-impressed current/electric dipoles, which could induce a magnetic field following Biot-Savart’s law, decrease the atmospheric electric field and generate dielectric-barrier discharge lightning affected by the coupling effect between the seismic and meteorological activities.

Since the seismic activities of the Matsushiro earthquake swarm occurred over a long period, i.e., 1965–1967, a considerable amount of reliable data has been reported on the seismic activities1,2,3,4,5,6 and the complexed seismo-EM activities7,8,9,10,11,12,13,14,15,16 as given in Figs 1,,22 and and3,3, where the source of the map in Fig. 1 is Geospatial Information Authority of Japan website ( created from the map data source: Landsat8 mosaic image (GSI, TSIC, GEO Grid/AIST), Landsat 8 image (courtesy of the U.S. Geological Survey), Geological Information Authority of Japan. The total number of felt earthquakes was ~600,000, and the number of large earthquakes with magnitude M  4 was 881,2. The highest magnitude reached M = 5.4. The focal depths ranged from 2–5 km1. The spatiotemporal variation in the seismic activities was divided into four stages1: The 1st–3rd periods corresponded to three peaks of the felt earthquakes (c.f. Fig. 2a). As the epicentral region developed stage by stage in the northeast region, the vertical and horizontal deformation of the crust accelerated rapidly in the 3rd stage (c.f. Fig. 2b)1. A fissure zone, 0.3 km wide and 1.8 km long, resulting from an underlying strike-slip fault, then developed near the northeast of Mt. Minakami where the open-surface cracks of a regular en echelon arrangement were successively constituted (c.f. Fig. 1 right-hand-insert)3,4,5. A large amount of ground water was expelled from the surface cracks accompanying the CO2 gas as shown in Fig. 27. In the 4th stage, the activities were diminished to the lower level1,2.

Figure 1
Locations of typical observation sites in Matsushiro area map.
Figure 2
Spatiotemporal variation of seismic activities in the 1965–1967 Matushiro earthquake swarm.
Figure 3
Temporal variation of seismo-EM phenomena.

The rupture mechanism is reasonably understood as follows18,19,20,21,22. When the magma that caused the eruption (0.35 Ma) of Mt. Minakami, a Pleistocene andesite volcano, cooled, large amount of andesitic and CO2-rich water was dehydrated and became highly pressurized beneath the impermeable and seismic-wave reflective layer lying at ~10 km deep20,22. The highly pressurized water then broke the impermeable layer and upwelled, passing through vertical flaws/cracks in the crust18,19,21, possibly caused by the shrinkage of the crust during the cooling process of the volcanic vent of Mt. Minakami. During upwelling of the water, the pressure decreased which resulted in degassing of CO2, volume-expansion and then the horizontal expansion of the opening of the vertical flaws/cracks (so called CO2/water eruption)18,19. The following supplementary information is added to the above-mentioned previous studies18,19,20,21,22: CO2 behaves as a supercritical fluid (sc-CO2) at a temperature T and pressure P above its critical point of Tc = 31 °C and Pc = 7.4 MPa, at which point sc-CO2 can dissolve various minerals in the crust causing corrosion cracking and local weakness in the surrounding crusts. Pc might be equivalent to the pressure at a depth of 3.3 km with a specific rock gravity sg of 2.3, where the depth corresponded to that of the seismic zone. A phase change from sc-CO2 to gas occurred, leaving dissolved mineral deposits in the crack channels (referred to as sc-CO2 eruption). The sudden rock ruptures likely generated remarkable earthquake pop-like sounds as often noted during the active seismic periods2. Furthermore, as the critical point of water is T = 374 °C and P = 22.1 MPa, the dehydrated water/CO2 might have been in a supercritical state (sc-water/CO2) below the estimated depth of 12.4 km with an sg of 2.3. Beneath the impermeable layer at ~10 km deep the phase change of sc-water and the resulting volume expansion of the super-heated steam might have broken the impermeable layer.

The geomagnetic observation at the site C in Fig. 1 showed an anomalous increase in the total geomagnetic field on the order of 10 nT associated with the seismic activities in the 3rd stage13 as shown in Fig. 3b. The geomagnetic variation exceeded a noise level of 0.2 nT, but could not be observed at the time of shallow the earthquakes of M  3.012. The apparent absence of correlation between the earthquake occurrence and the geomagnetic variation indicated that the observed geomagnetic variation should not simply be attributed to the coseismic stress change.

The AEF observation from October to November 1966 demonstrated that the electric fields were often decreased by the felt-earthquakes16 (c.f. Fig. 3a). A similar decrease in the AEF may occur during an approach of a storm cloud with an electric dipole in the upper positive and lower negative range23. However, the AEF decreased with rapid alternative variation in the positive and negative directions on 22 Oct. 1966 (c.f. Fig. 3a), coinciding with favourable weather throughout the day as well as before and after. This conditions might indicate that electric dipoles generated in the underground seismic zone affected the ground electrification in a similar manner as the storm-cloud dipole. The AEF alternative variation became much more enhanced on the night of 14 Nov. 1966 (c.f. Fig. 3a). Although the weather recovered during the observed period, a cold front passed through during the day. The coupled interaction of seismic activities with meteorological variation, therefore, may have influenced the enhanced AEF variation.

The EQLs were photographed for the first time by T. Kuribayashi, an amateur photographer living in Matsushiro, as shown typical photos 1–3 in Fig. 1. Yasui14 evaluated the EQL reports as follows: Among 34 events recorded in the EQL candidates, 18 events could be identified; these events occurred in the period from the 1st stage until the early 3rd stage. As for the other 16 events, there were some fears of misjudgments due to other luminous phenomena, such as lightning or the sun or the moon, artificial lights and others (c.f. Fig. 3b). As seen in photos 1–3 in Fig. 1, the EQLs appeared in the mountainous areas within the seismic zone. Photo 3 was taken using a camera with a fish-eye lens on 26 September 1966, and it shows the scale of luminous zone, i.e., approximately 8 km (the yellow-coloured dotted circle in Fig. 1). The luminosity continued for 96 seconds. The felt earthquakes occurred at 03:11 and 03:37 L.T. but not at the time of the EQLs. The lack of direct coincidence between the EQLs and the earthquake occurrence indicated that the EQLs could not be explained by any model due to the coseismic stress change mechanism. Unlike the geomagnetic variation, which was notable in the 3rd stage, the EQLs were commonly observed through the 1st–3rd stages (c.f. Fig. 3b). This finding indicates that even if the source mechanism is identical, the way in which the source mechanism affects both seismo-EM phenomena differs. The EQLs are likely related to a certain mechanism involving the upwelling of CO2/water and the atmospheric condition, as the EQLs often occurred when a cold-front passed through14.

In terms of the stress-change mechanisms, the piezomagnetic model24 and the modified piezoelectric model25 have been proposed to explain either the geomagnetic variation or the EQLs. These models appear to provide a plausible explanation regarding whether the observed seismo-EMs were coseismic: however, this was not the case. A conceptual model was recently proposed that the electric fields at the ground-to-air interface due to an influx of stress activated charge carriers (positive holes) become so steep as to trigger corona discharges26,27. Corona discharges generate audible noise, but the Matsushiro EQLs did inaudible14.

An electrokinetic flow model28, where the electric current is generated due to groundwater flow passing through open cracks, seemed to be consistent with the observations in the 3rd stage when the outflow rate of water at the fissure zone increased abruptly (c.f. the blue-coloured line in Fig. 2a). The estimated geomagnetic change of 4 nT28 was still lower than the observed change of 5–15 nT. The electrokinetic flow model would require rapid and implausibly continuous fluid flow29. This model, furthermore, could not explain why the EQLs occurred not only at the 3rd stage, when the electrokinetic effect might be significant, but also at the 1st and 2nd stages when it might be less important.

In this context, a new source mechanism for the complexed seismo-EMs is needed to elucidate the underlying causal relationship between the hydro-mechanical and the electromagnetic processes. We consider that exoeletrons are emitted transiently from various trapped sites in lattice defects on fresh fractured surfaces30,31,32, and any gases that interacts with fractured fresh surfaces might be negatively electrified due to electron attachment reactions32. Taking into account the coupled electromagnetic interaction of the cracks with the gas flowing in as a working hypothesis for the seismo-EM phenomena33, we conducted laboratory experiments of uniaxial rock rupture coupled with high-pressure CO2 flow as described in the Methods: Fracture tests. The typical well-measured signals are shown in Fig. 4a for quartz diorite, the main rock constitution of the Matsushiro epicentre crust, with and without CO2 flow. The electric current was successfully measured for the rocks with CO2 flow, followed for approximately 2 msec after fully development of the crack. After the signal peak, the small alternating signals followed for ~10 msec. The signal variations, observed both with and without CO2 flow, should be attributed to the vibrations excited by the electric field fluctuation between the charge-separated mating fracture surfaces. The observed peak currents I(labo) are shown in Fig. 4b as a function of the gas-fracture interacting area S(labo), defined as [thickness of the rock sample] × [width of collecting electrode]. The relationship between I(labo) and S(labo) can be expressed by

Figure 4
Electric currents due to the fracturing rocks (left) and the induced geomagnetic effects (right).
An external file that holds a picture, illustration, etc.
Object name is srep44774-m1.jpg

The present results suggest that as many water/CO2 gas-bearing pores are distributed in the shallow seismic zone, the dipole generation gave rise to short-term transient and temporally decaying electrification at the ground level. This process may provide a plausible explanation for the AEF variation (c.f. Fig. 3a) during the active earthquake periods.

During the 3rd stage, as the geomagnetic intensity level gradually increased13, the successive opening of en echelon-type cracks, resulting from an underlying strike-slip fault, was observed in the fissure zone3,4,5. The laboratory experiments suggest that whenever an en echelon crack array extended from the deep seismic zone formed as a result of coupled left-lateral movement interactions and sc-CO2 eruptions, large vertical dipole currents I (field) were generated with induced observable geomagnetic variation following Biot-Savart’s law. Based on the model shown in Fig. 4(c), we then estimated the geomagnetic variation |ΔB| assuming a line dipole current I(field) element as follows.

An external file that holds a picture, illustration, etc.
Object name is srep44774-m2.jpg

where μ0 is the permittivity of free space μ0 =  × 10−7 WbA−1m−1, θ1and θ 2 denote the angle shown in Fig. 4c, and R is the distance between the observation site and the fissure zone. We assumed that Eq. (1) holds even during the field-level rupture processes. The degree of geomagnetic variation |ΔB| is shown in Fig. 4d as a function of R for dipole lengths of 1.5 km and 3 km, where the dipole currents were estimated to be 450 A and 900 A using Eq. (1) for rupture areas of 0.3 km × 1.5 km = 0.45 km2 and 0.3 km × 3 km = 0.9 km2, respectively. These rupture areas are comparable to those associated with earthquakes of M = 3.6 and 3.9, respectively34. The estimated |ΔB| values are in reasonable agreement with those of the observed total geomagnetic variation level13. The present results suggest that because the fissure zone with en echelon crack arrays extended 1 km to 2.3 km from the geomagnetic observation site C in Fig. 1, it was located within a sensitive region that was not too distant or close for the detection of seismo-geomagnetic variation.

Assuming the rapid alternative AEF variations, such as those observed on 14 Nov. 1966 (c.f. Fig. 3a), which might be the source mechanism for the EQLs, we conducted the discharge experiments of the model fields as described in the Methods: Discharge tests. Figure 5a and b shows the model field of the embankments with and without twigs of an evergreen needle-shaped tree at an AC 15 kHz, respectively. The case in which twigs were included showed pale bluish-purple-coloured lights, similar to those of the eye-witnessed EQLs, but less overall lights than the case without the twigs. The twigs clearly act as a dielectric material for dielectric barrier discharge (DBD); so-called silent (inaudible) discharge, with lower discharge current than the spark discharge (SD) with audible noises, as seen in Fig. 6. The spectrum lines of the lights mainly consist of 2nd positive nitrogen molecules within the wavelength range of 380–440 nm (c.f. Fig. 5c). The DBD due to AEF variations resulting from the coupled interaction of thunderstorm and seismic activities thus might be responsible for the EQLs, where trees covering the mountainous area played an important role in the silent EQLs.

Figure 5
Laboratory discharge experimental results using the model-field.
Figure 6
Laboratory discharge experimental results using the model-field (continued).

In summary, a set of hydro-mechanical processes and the associated seismo-EMs are consistently explained with the present model as illustrated in Fig. 7, which improved the previous studies20,22. This is the first reported comprehensive modelling study to describe the causal relationships among a diverse set of seismic anomalies associated with the 1965–1967 Matsushiro earthquake swarm.

Figure 7
Causal mechanisms of seismo-EM phenomena during the 1965–1967 Matsushiro earthquake swarm.


Fracture tests

The experimental set up involving to a universal testing machine is shown in Fig. 8a. A flat-ended chisel, made of hardened carbon steel S45C with Brinell hardness of ~230, was loaded against an as-received rock block (normally 50 × 50 × 20 mm in size); quartz diorite, crustal rock of the Matsushiro epicentral zone. CO2 gas at the pressure ranging from 0.3 to 0.5 MPa and at room temperature was introduced into a flow-channel inside the chisel. At an instant when the rock was subjected to guillotine-type fracture at a load of ~12 kN, the pressurized CO2 gas flowed into the crack gap of the rock passing through the open slit (1 mm × 20 mm) equipped with the flat-ended chisel. To this end, a seal made of a 1-mm-thick PTFE sheet was set between the chisel of the flat-ended area around the open slit and the rock surface (c.f. Fig. 8a: upper insert). Undesirable three-point-type fractures of the rock were often observed, where the crack initiates from the bottom side to the top. Therefore, to ensure the crack could initiate from the upper side to allow CO2 gas to flow, the four corners of the upper rock surface were suppressed loosely with clamps (not shown in Fig. 8a). A clean copper mesh #20 electrode was set on the bottom of the rock (see the lower insert of Fig. 8a and biased at +77 dc volts to collect negatively charged gas flow. Two-ply sheets consisting of PTFE, 1 mm thick each, were set between the sample rock and the basement rocks of fine-grained granite in order to assist the horizontal displacement when the crack was opened. After the tests, work clay was pressed on the cracked rock surface to measure the gap size which typically ranged from 0.4 mm to 1 mm.

Figure 8
Experimental set up (a) Universal tester for measurement of electrified gas flow associated with rock fracture: the upper inset shows the PTFE sealed slit at the flat ended chisel and the lower insert shows a copper-mesh electrode equipped on the back ...

Discharge tests

An experimental set up for discharge tests in an open-air room is shown in Fig. 8b, where a model field of embankment filled with wet soil (the resistivity was ~260 ohm. m) with or without twigs of an evergreen needle-shaped tree (Monterey cypress Goldcrest) was set between the electrodes with separation distance of 10 cm and subjected to high voltage up to 20 kV DC or AC 1 kHz and 15 kHz. Discharge lights were photographed using a CMOS camera with ISO sensitivity of 25,600 at an exposure time of 5 sec. The spectrum was analyzed using a high-sensitivity spectral radiance meter. The discharge currents were measured in terms of voltage drop across 2.8 k resistor as shown in the block diagram of Fig. 8b9.

Additional Information

How to cite this article: Enomoto, Y. et al. Causal mechanisms of seismo-EM phenomena during the 1965–1967 Matsushiro earthquake swarm. Sci. Rep. 7, 44774; doi: 10.1038/srep44774 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:


The authors thank Genesis Research Institute, Inc. for the financial support; Konica Minolta, Inc. and Tamaoki Electronic Co. Ltd for lending of the spectrometer and the high power supply, respectively and Hideki Kuribayashi, Matsushiro Earthquake Center, Earthquake Research Institute, The University of Tokyo & Geospatial Information Authority of Japan for providing various information.


The authors declare no competing financial interests.

Author Contributions Y.E. planned the research and prepared the manuscript. Y.E. and T.Y. performed the experiments. N.O. was responsible for launching the present project. Y.E., T.Y. and N.O. discussed on the results.


  • Hagiwara T. & Iwata T. Summary of the seismographic observation of Matsushiro swarm earthquakes. Bull. Earthquake Res. Inst. Univ. Tokyo 60, 485–515 (1968).
  • Japan Metrological Agency. Report on Matsushiro earthquake swarm Augast 1965-Decenber 1967 (in Japanese with English Abstract). Tech. Report JMA 62, 1–556 (1968).
  • Nakamura K. & Tsuneishi Y. Ground cracks at Matsushiro probably of understanding strike-slip fault origin, I-Preliminary Report. Bull. Earthquake. Instit. Univ. Tokyo 44, 1371–1384 (1966).
  • Nakamura K. & Tsuneishi Y. ibid II-The Matsushiro earthquake fault. Bull. Earthquake. Instit. Univ. Tokyo 45, 417–471 (1967).
  • Tsuneishi Y. & Nakamura K. Faulting associated with the Matsushiro Swarm earthquakes. Bull. Earthquake. Instit. Univ. Tokyo 48, 29–51 (1970).
  • Izutuya S. Revisited results of levelling surveys during the Matsushiro earthquake swarm (in Japanese with English Abstract). Bull. Earthquake Res. Inst. Univ. Tokyo 50, 273–280 (1975).
  • Rikitake T. et al. . Geomagnetic and geoelectric studies of the Matsushiro earthquake swarm (1). Bull. Earthquake Res. Inst., Univ. Tokyo 44, 363–408 (1966).
  • Rikitake T. et al. . Geomagnetic and geoelectric studies of the Matsushiro earthquake swarm (2). Bull. Earthquake Res. Inst., Univ. Tokyo 44, 409–418 (1966).
  • Rikitake T. et al. . Geomagnetic and geoelectric studies of the Matsushiro earthquake swarm (3). Bull. Earthquake Res. Inst., Univ. Tokyo 44, 1335–1370 (1966).
  • Rikitake T. et al. . Geomagnetic and geoelectric studies of the Matsushiro earthquake swarm (4). Bull. Earthquake Res. Inst., Univ. Tokyo 44, 1735–1758 (1966).
  • Rikitake T. et al. . Geomagnetic and geoelectric studies of the Matsushiro earthquake swarm (5). Bull. Earthquake Res. Inst., Univ. Tokyo 45, 395–416 (1967).
  • Rikitake T. et al. . Geomagnetic and geoelectric studies of the Matsushiro earthquake swarm (6). Bull. Earthquake Res. Inst., Univ. Tokyo 45, 919–944 (1967).
  • Yamazaki Y. & Rikitake T. Local anomalous changes in the geomagnetic field at Matsushiro. Bull. Earthquake Res. Inst., Univ. Tokyo 48, 637–643 (1970).
  • Yasui Y. A study on the luminous phenomena accompanied with earthquake (Part I) (in Japanese with English Abstract). Mem. Kakioka Mag. Obs. 13, 25–61 (1968).
  • Derr J. S. Earthquake lights: a review of observations and present theories. Bull. Seismo. Soc. Ameri. 63, 2177–2187 (1973).
  • Kondo G. The variation of the atmospheric electric field at the time of earthquake. Mem. Kakioka Mag. Obs. 13, 11–23 (1968)
  • Kasuga I. Aspects of the relation of thermal water and Matsushiro earthquakes in Kagai hot spring area, Nagano Prefecture (in Japanese with English Abstract). J. Geogr. 76, 16–26(1967).
  • Tsuneishi Y. Chishitsu-Kozo no Kagaku (in Japanese); Science of Geological Structure: translated title), (ed. Kimura T.) Ch.4, 370 (Asakura Syoten, 1984).
  • Sasai Y. A generating mechanism of Matsushiro, Central Japan, earthquake swarm: spontaneous hydraulic fracturing (in Japanese). Proc. Conductivity Anomaly Res. Group, 181–195 (1994).
  • Yoshida N., Okuzawa T. & Tsukahara H. Origin of deep Matsushiro earthquake swarm fluid interred from isotope ratios (in Japanese with English Abstract). Jisin 55, 207–216 (2002).
  • Cappa F., Rutqvist J. & Yamamoto K. Modelling crustal deformation and rupture processes related to upwelling of deep CO2-rich fluids during the 1965–1967 Matsushiro earthquake swarm in Japan. J. Geophys. Res. 114, B10304–B10304 (2009).
  • Tsukahara H. Understanding the meaning of water discharge at the Matsushiro earthquake swarm area (in Japanese with a figure noted in English). Abstract, Japn. Geosci. Union Meeting, May 23–28, Makuhari, Chiba, U004–08 (2010).
  • Uman M. A. Lightning Discharge, Academic Press, p112 (1987).
  • Sasai Y. The piezomagnetic field associated with the Mogi model. Bull. Earthquake Inst. 54, 1–29 (1979).
  • Ikeya M. & Takaki S. Electromagnetic fault for earthquake lightning. Japan. J. Appl. Phys. 35, L355–L357 (1996).
  • Freund F. T. Toward a unified solid state theory for pre-earthquake signals. Acta Geophysica 58, 719–766 (2010).
  • Thèriault R., St-Laurent F., Freund F. T. & Derr J. S. Prevalence of earthquake lights associated with rift environments. Seis. Res. Lett. 85, 159–178 (2014).
  • Mizutani H. & Ishido T. A new interpretation of magnetic field variation associated with the Matsushiro earthquakes. J. Geomag. Geolectr. 28, 179–188 (1976).
  • Johnston M. J. S., Mueller R. J. & Sasai Y. Magnetic field observations in the near-field the 28 June 1992 Mw 7.3 Landers, California, earthquakes. Bull. Seis. Soc. Amer. 84, 792–798 (1994).
  • Lewis D. R. Exoelectron emission phenomena and geological applications. Geolog. Soc. Amer. Bull. 77, 761–770 (1966).
  • Enomoto Y. & Hashimoto H. Emission of charged particles from indentation fracture of rocks. Nature 346, 641–643 (1990).
  • Scuriero L., Dickinson J. T. & Enomoto Y. The electrification of flowing gases by mechanical abrasion of mineral surfaces. Phys. Chem. Minerals 25, 566–573(1998).
  • Enomoto Y. Coupled interaction of earthquake nucleation with deep Earth gases: a possible mechanism for seismo-electromagnetic phenomena. Geophys. J. Inter. 191, 1210–1214 (2012).
  • Hanks T. C. & Bakun W. H. A bilinear source-scaling model for M-log A observation of continental earthquakes. Bull. Seismol. Soc. Amer. 92, 1841–1846 (2002).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group