*Kosuke HEKI1
(1.Toyama Industrial Technology Center)
Keywords:Tohoku-Oki earthquake, Precursor phenomena, Geomgnetic, Declination, GPS-TEC, Modelling
The GPS-TEC enhancement stating ?40 minutes before the 2011 Tohoku-Oki earthquake has been observed (Heki, GRL, 2011). The geomagnetic declination change was confirmed nearby the fracture zone, at Easashi (ESA), Mizusawa (MIZ) by GSI and at Kakioka (KAK), Kanozan (KNZ) by JMA in respect to Kanoya (KNY) in synchronization with the GPS-TEC anomaly(Heki & Enomoto, JGR, 2013).These anomalies satisfy the criteria of earthquake precursor candidate (Wyss, AGU, 1991).The magnetic declination; the difference between the direction of the horizontal components H of the Earth's magnetic field and the magnetic north is normally ?6.9 degree westward (=-415.7 arc min) from true north at the ESA site, but, as seen in the Figure, ∆D ([ESA]-[KNY]) gradually changes to the positive direction (eastward) starting from 40 minutes to the maximum ∆D value of ∆D? +0.32 arc min (= 9.31x10-5 rad) just before the main shock. This change should be affected by generation of preseismic magnetic field ?∆B?. As ∆D is small, we can approximate the relationship between ∆H, ∆D and ?∆B? as shown in an insert of the Figure; i.e.?∆B?[nT] ? ?H? [nT]?∆D[rad] The amount of H is normally 29037nT at ESA. The ?∆B? is then 29037 nTx9.31x10-5 rad = 2.70 nT, which is in agreement with the observed preseismic variation of ∆H (?2.38 nT) as seen in the Figure. The preseismic geomagnetic field |∆B|, resulting from the time-varying current at the earthquake nucleus zone by Biot-Savalt's law, is expressed, by assuming the time-varying source current element of the length wc t*, as:?∆B?=μ0 ?I? sinθ wct*/4πR2, where μ0 is the permittivity of free space, ?I? is the pressure-impressed current, which is 170kA in the Tohoku-Oki earthquake (Enomoto & Heki, GJI, submitted), is an angle shown in the Figure, wc is the earthquake nucleation size, t* is the normalized preseismic time duration, R is the distance between observation site from the epicenter. The present model of the above equation at t*=1 gives ?∆B?=1.78 nT with R=181km: distance between the ESA and the epicenter of which the agreement with the observed value of 2.38nT is rather well.