CC BY
72
Proceedings of the MSTU, Vol. 15, No. 2, 2012 pp.471-476
UDC 550.388.2
Modeling of the ionosphere response on the earthquake preparation
M.I. Karpov, O.V. Zolotov, A.A. Namgaladze
Polytechnic Faculty of MSTU, Physics Department
Abstract. Seismo-ionosphere coupling processes have been investigated considering the GPS observed anomalous ionospheric Total Electron Content (TEC) variations before strong earthquakes as their precursors. The numerical simulations' results of the TEC response on the vertical electric currents flowing between the Earth and ionosphere during the earthquake (EQ) preparation time have been performed. Model experiments have been carried out using the Upper Atmosphere Model. The following currents' parameters were varied in: (i) direction (to or from the ionosphere); (ii) latitudinal zone of the sources' (EQ epicenters) location; (iii) currents' configuration: (1) grid nodes with "straight" currents were surrounded by "border" grid points with currents of opposite direction ("return" currents); (2) the "return" currents were spread out over the globe; (3) without "return" currents. Numerical simulations have shown that electric currents with density of 4Т0-8 A/m2 over the area of about 200 km in longitude and 2500 km in latitude produce both positive and negative TEC disturbances with magnitude up to 35 % in agreement with GPS TEC observations before EQs.
Аннотация. Исследованы проявления сейсмо-ионосферных связей - аномальных возмущений полного электронного содержания (ПЭС) ионосферы перед землетрясениями как их предвестников. Представлены результаты численного моделирования полного электронного содержания (ПЭС) в зависимости от вертикальных электрических токов между Землей и ионосферой перед землетрясениями. Модельные эксперименты проводились с помощью модели верхней атмосферы UAM. Варьировались следующие параметры токов: направление, широтное расположение источников (эпицентров
землетрясений), конфигурации токов: узлы с "прямыми" токами были окружены пограничными узлами с токами противоположного направления ("возвратными" токами); "возвратные" токи были распределены по всей земной поверхности; без "возвратных" токов. Расчеты показали, что токи плотностью 4Т0-8 A/м2 над областью около 200 км по долготе и 2500 км по широте вызывают положительные и отрицательные возмущения ПЭС магнитудой до 35 % и более, что соответствует GPS наблюдениям перед землетрясениями.
Key words: Total Electron Content, seismogenic electric current, earthquake precursors, ionosphere, UAM, modeling
Ключевые слова: полное электронное содержание, сейсмический электрический ток, предвестники землетрясений, ионосфера,
UAM, моделирование
1. Introduction
Seismogenic electric fields generated near the tectonic faults are important elements of the global electric circuit (e.g., Harrison et al., 2010; Pulinets, Boyarchuk, 2004, and references therein; Freund et al., 2009; Pulinets, Ouzounov, 2011; Sorokin et al., 2005; 2006; 2007). Radioactive radon emanation (Pulinets, Boyarchuk, 2004) and "positive holes" mechanism (Freund et al., 2009) both are acting as air ionization sources over the faults. According to Freund (2009) stressed rocks could generate electric current with density of 0.51.25 pA/m2. Sorokin et al. (2005; 2006; 2007) calculated that such currents within an area about 200 km in radius are able to create ionosphere electric field disturbances of several mV/m. Such electric fields in turn may disturb the ionosphere F-2 electron density and total electron content (TEC) up to 50 % and more (Namgaladze et al., 2007; 2009a; 2009b).
Many other studies, observations and experiments demonstrate that vertical electric currents with such densities exist in nature. For example, Krider and Musser (1982) used data of time variations in thunderstorm electric fields, both aloft and at the ground, provided by a large field mill network. They estimated that the average Maxwell current density under a small Florida thunderstorm ranged from 1 to 4 nA/m2 for a short period of time. The air conductivity and the vertical electric field were measured over thunderstorms from U -2 airplane by Blakeslee, Christian (1989). It was shown that the current with density ranging from 0 to 20 nA/m2 flows upward from thunderstorm areas. Le Mouel et al. (2010) observed transient electric potential variations in a standing poplar tree equipped with electrodes along its height. It was obtained that current with typical value of 2.5 pA flows through the average cross section of the trunk of 0.5 m2. It corresponds to current density of 5 pA/m2.
471
Karpov M.I. et al. Modeling of the ionosphere response...
Fig. 1. Magnetic maps of model TEC deviations (%) relative to the quiet conditions for "straight" currents located at 5° lat. (from top to bottom): "return" currents spread out over the globe except points with "straight" currents; without "return" currents; the same as at top but currents flow from the ionosphere to the Earth; "return" currents are located at the grid nodes nearest to the "straight" currents. Stars - grid points with "straight" currents. Diamonds - the magnetically conjugated points. Orange circle - the subsolar point. Black curves -
terminator lines
-25-20-15-10 -5 0 5 10 15 20 25
11111______I-1
о
60' -
45' -
я 30' -
— 15 -
о O' -
о -15
С ftfl -30' -
я -45'
s -60 -
60'
45' -
30' -
15‘
O' -
-15’
-30' -
-45
-60' -
60'
45' -
30 -
15
‘i ' I TEC diff. (%) maps
u:i)T 04UT 06UT 08UT 10UT 12UT 14UT
18UT :uci 221 I
j —. c 1 > 1 1 1 1 ■.. 1 dtMrM TTT' i *4 / T' ' 1 i■ ' t f 4-r-r-r-r-r-r-f r^ "1 ' ' ' 1 1 ' t t ' 1 ' 1 ' ' ' ' ' ' 1 1 ' ' 1 "1 ' ' ' 1 ' ' ’ ' ' ' 1 i ' 1 ’ \mmi \ * 1 ' ' I ' ' '
i — c ' 1 ' 1 1 1 *■1 '1 1 * TtT i | .4 мм/ У f * : : 1 ’ - ’ r^ ' 1 1 1 1 1 1 Г •им 1 1 if : y. U X ' ' 1 ’ 1 ' 1 '
I ■■■'■■ i I ■•■■■■ 1 • 1 3 •МИ». ! J / ■ ■■•!■• 1 ***"
- . —*** J...I ...1.. .1.. f ; V • 1 I i.M.nrr.rpt \ ««««« \
O'
-15'
-30'
-45'
-60
3 60'
45'
я 30'
— 15
о 0'
о -15
с зо -30
я -45
2 -60
-40^20’ 0“ 20"
Magnetic longitude
-40^20” O' 20' -40^20' 0" 20” -4СИ-20' 0” 20” -4CF-20' O' 20'
Magnetic longitude Magnetic longitude
Fig. 2. Magnetic maps of model TEC deviations (%) relative to the quiet conditions for "straight" currents located at 15° lat.
472
Magnetic latitude Magnetic latitude Magnetic latitude Magnetic latitude
Proceedings of the MSTU, Vol. 15, No. 2, 2012 pp.471-476
-35 -25 -15 -5 5 15 25 35
.......................1.......................
TEC dilT. (%) maps
-40-20 0 20
Magnetic longitude
-40^20 0 20 -40-20 O' 20'
Magnetic longitude
-40-20' O' 20' -40-20' 0 20
Magnetic longitude
Fig. 3. Magnetic maps of model TEC deviations (%) relative to the quiet conditions for "straight" currents located at 30° lat.
Fig. 4. Magnetic maps of model TEC deviations (%) relative to the quiet conditions for "straight" currents located at 45° lat.
473
Karpov M.I. et al. Modeling of the ionosphere response...
Previous numerical simulations (Namgaladze, 2010; Namgaladze et al., 2011; Namgaladze, Zolotov, 2011) showed that point electric current sources larger than 10-6A/m2 generated too strong TEC disturbances. Similar point sources but with densities less than 10"9-10"8 A/m2 generated TEC disturbances not exceeding 10-20 %. Those values were smaller than ones usually taken into consideration in pre-earthquake TEC modifications studies. The vertical electric current with density of about 2T0"8 A/m2 flowing between the fault and the ionosphere set at the area of about ~200 km x -4000 km may create electric fields generating the TEC increases up to ~50 % at night-time as observed before the Haiti Jan. 12, 2010 earthquake. Appearance of the terminator and subsolar point caused depression down to full destruction of electric potential generated by external electric current. TEC disturbances were also modified: "escaped" from terminator to dark-side of the ionosphere with following full destruction, but with time lag relative to the electric potential. Next night the electric potential and corresponding TEC disturbances restored.
The principle goal of this paper is to present the results of numerical simulations of ionosphere electric fields and corresponding TEC disturbances depending on various parameters of electric currents between the Earth and ionosphere: (1) latitudinal position, (2) general spatial orientation of the vertical electric currents' distribution and (3) direction (to or from the ionosphere).
2. Simulation methods
Model experiments were carried out using the Upper Atmosphere Model - global numerical model of the Earth's upper atmosphere (thermosphere, ionosphere, plasmasphere and inner magnetosphere). The model calculates main physical parameters of upper atmosphere such as densities, temperatures and velocities of the neutral (O, O2, N2, H) and charged (O2+, NO+, O+, H+ and electrons) species by numerical integration of the continuity, momentum and heat balance equations as well as the equation of the electric field potential (Namgaladze et al., 1988; 1991; 1998a; 1998b).
To build model difference maps of the TEC we first performed a regular calculation without any additional electric current sources to use those results as quiet background values. Then, an external electric current ("straight" current) flowing between the lower atmosphere and the ionosphere presumably of seismogenic origin was used as a model input for the calculations of the ionospheric field and the corresponding TEC variations. These additional sources of the electric current with the magnitude of 4-10-8 A/m2 were switched on at the lower boundary (80 km) in the UAM electric potential equation, which was solved numerically jointly with all other UAM equations for neutral and ionized gases.
Currents sources were switched on at five numerical grid nodes which correspond to the size of the area about 200 km along the parallel and about 3000 km along the meridian depending on latitude. The latitudinal position (5°, 15°, 30° and 45° mag. lat.) and direction of the currents (to or from the ionosphere) were varied in. The grid nodes with vertical electric currents ("straight" currents) were surrounded by boundary grid points with currents of the opposite direction ("return" currents) so that the total current in the global electric circuit didn't change. In other investigated configuration "return" currents were spread out over the globe except the points with "straight" current. The configuration without "return" currents were studied as well.
3. Discussion and summary
According to UAM simulations (Fig. 1 -4) sources of the additional electric currents flowing between the Earth and ionosphere presumably of seismogenic origin with current density 4T0-8 A/m2 over the area of about 200 km in latitude and 3000 km in longitude produce both positive and negative TEC disturbances. The size and geographical position of the disturbed area didn't change for a few hours. Approaching of the sunlit ionosphere is followed by disappearance of the disturbances.
The results of the numerical computations are in good agreement with the observations of TEC before strong earthquakes at middle and low latitudes in the spatial scales and amplitude characteristics (Liu et al., 2000; 2006; 2011; Pulinets, Boyarchuk, 2004; Pulinets, Ouzhounov, 2011; Zakharenkova, 2006; 2008; Zolotov et al., 2012).
As Figs. 1-4 show changing the direction of the currents causes mirror reverse of the positive and negative disturbance relatively magnetic meridian crossing the epicenter area and magnetically conjugated area except the case for currents at 45° lat.
Various configurations of the currents have different influence on the TEC disturbances. The configuration without "return" currents creates TEC disturbances very similar to those produced by configuration with "return" currents spread out over the whole globe. The magnitude of the TEC disturbances produced by "return" currents located in boundary grid nodes is lower than magnitudes of the disturbances created by other configurations.
The most powerful magnitudes (35 % and larger) are produced by currents located at 30° and 45° mag. lat. The shift of the vertical electric current from 30° magnetic latitude towards the magnetic equator leads to less
474
Proceedings of the MSTU, Vol. 15, No. 2, 2012 pp.471-476
strong disturbances by magnitude. Sources for latitudes less than 30° generate pronounced modification of the Appleton anomaly.
It should be noted that similar numerical calculations were performed by (Kuo et al., 2011). The difference is that their currents flow the over area 200 km x 30 km that is much smaller than in our case. "Rock currents" were switched on near ground surface (not at the lower boarder of the ionosphere). Then they used the empirical model of the atmosphere's conductivity modified from some considerations. They simulated only TEC response at the epicenter region and did not consider magnetically conjugated ones. The simulations were local (not for the globe as in our case) with "single" point sources. According to their simulations using SAMI2 model the current density of 0.2-10 pA/m2 in the earthquake fault zone can cause TEC variations of up to 2-25 % in the daytime ionosphere and the density of 0.01-1 pA/m2 can lead to nighttime TEC variations of 1-30 %. Such huge currents are required due to a too small area of the current generation in their simulation. Beside of this, our and their model results do not contradict each other.
References
Blakeslee R.J., Christian H.J. Electrical measurements over thunderstorms. J. Geophys. Res., v.94, N D11, p. 13,135-13,140, 1989.
Freund F.T., Kulahci I.G., Cyr G., Ling J., Winnick M., Tregloan-Reed J., Freund M.M. Air ionization at rock surfaces and pre-earthquake signals. Journal of Atmospheric and Solar-Terrestrial Physics, v.71, p.1824-1834, 2009.
Harrison R.G., Aplin K.L., Rycroft M.J. Atmospheric electricity coupling between earthquake regions and the ionosphere. J. Atmos Solar-Terr Phys., v.72, N 5-6, p.376-381, 2010.
Krider E.P., Musser J.A. Maxwell currents under thunderstorms. J. Geophys. Res., v.87, N C13, p. 11,17111,176, 1982.
Kuo C.L., Huba J.D., Joyce G., Lee L.C. Ionosphere plasma bubbles and density variations induced by pre-earthquake rock currents and associated surface charges. J. Geophys. Res., v.116, A10317, doi:10.1029/2011JA016628, 2011.
Le Mouel J.-L., Gibert D., Poirier J.-P. On transient electric potential variations in a standing tree and atmospheric electricity. ComptesRendus Geoscience, v.342, Issue 2, p.95-99, 2010.
Liu J.Y., Chen Y.I., Pulinets S.A., Tsai Y.B., Chuo Y.J. Seismo-ionospheric signatures prior to M>6.0 Taiwan earthquakes. Geophys. Res. Letters, v.27, N 19, p.3113-3116, doi:10.1029/2000GL011395, 2000.
Liu J.Y., Chen Y.I., Chuo Y.J., Chen C.S. A statistical investigation of preearthquake ionospheric anomaly.
J. Geophys. Res., v.111, p.A05304, doi:10.1029/2005JA011333, 2006.
Liu J.Y., Le H., Chen Y.I., Chen C.H., Liu L., Wan W., Su Y.Z., Sun Y.Y., Lin C.H., Chen M.Q. Observations and simulations of seismoionospheric GPS total electron content anomalies before the 12 January 2010 M7 Haiti earthquake. J. Geophys. Res., v.116(A04302), doi:10.1029/2010JA015704, 2011. Namgaladze A.A., Korenkov Yu.N., Klimenko V.V., Karpov I.V., Bessarab F.S., Surotkin V.A., Glushchenko T.A., Naumova N.M. Global model of the thermosphere-ionosphere-protonosphere system. Pure and Applied Geophysics, v.127, N 2/3, p.219-254, doi:10.1007/BF00879812, 1988. Namgaladze A.A., Korenkov Yu.N., Klimenko V.V., Karpov I.V., Surotkin V.A., Naumova N.M. Numerical modeling of the thermosphere-ionosphere-protonosphere system. Journal of Atmospheric and Terrestrial Physics, v.53, N 11/12, p.1113-1124, doi:10.1016/0021-9169(91)90060-K, 1991.
Namgaladze A.A., Martynenko O.V., Volkov M.A., Namgaladze A.N., Yurik R.Yu. High-latitude version of the global numeric model of the Earth's upper atmosphere. Proc. of the MSTU, v.1, N 2, p.23-84, 1998a. URL: http://goo.gl/8x9f2
Namgaladze A.A., Martynenko O.V., Namgaladze A.N. Global model of the upper atmosphere with variable latitudinal integration step. Int. J. of Geomagnetism and Aeronomy, v.1, N 1, p.53-58, 1998b.
Namgaladze A.A., Klimenko M.V., Klimenko V.V., Zakharenkova I.E. Physical mechanism and mathematical modeling of earthquake ionospheric precursors registered in total electron content. Geomagnetism and Aeronomy, v.49, N 2, p.252-262, doi:10.1134/S0016793209020169, 2009a. Namgaladze A.A., Shagimuratov I.I., Zakharenkova I.E., Zolotov O.V., Martynenko O.V. Possible mechanism of the TEC enhancements observed before earthquakes. XXIV IUGG General Assembly, Perugia, Italy, 02-13 July, 2007.
Namgaladze A.A., Zolotov O.V., Zakharenkova I.E., Shagimuratov I.I., Martynenko O.V. Ionospheric total electron content variations observed before earthquakes: Possible physical mechanism and modeling. Proc. of the MSTU, v.12, N 2, p.308-315, 2009b. ArXivID: 0905.3313 URL: http://goo.gl/A8cLx Namgaladze A.A. Physical model of earthquake ionospheric precursors. Abstract NH24A-03 presented at 2010 Fall Meeting, AGU, San Francisco, Calif., 2010.
475
Karpov M.I. et al. Modeling of the ionosphere response...
Namgaladze A.A., Zolotov O.V. Ionospheric effects from different seismogenic electric field sources. XXXth URSI General Assembly and Scientific Symposium, Istanbul, Turkey. IEEE, August. doi:10.1109/URSIGASS.2011.6051040, 2011.
Namgaladze A.A., Zolotov O.V., Prokhorov B.E. The TEC signatures as strong seismic event precursors.
XXXth URSI General Assembly and Scientific Symposium. Istanbul, Turkey. IEEE, August. doi:10.1109/URSIGASS.2011.6051048, 2011.
Pulinets S.A., Boyarchuk K. Ionospheric precursors of earthquakes. Berlin, Springer, 315 p., 2004.
Pulinets S., Ouzounov D. Lithosphere-atmosphere-ionosphere coupling (LAIC) model - an unified concept for earthquake precursors validation. J. Asian Earth Sci., v.41, N 4-5, p.371-382, doi:10.1016/j.jseaes.2010.03.005, 2011.
Sorokin V.M., Chmyrev V.M., Yaschenko A.K. Theoretical model of DC electric field formation in the ionosphere stimulated by seismic activity. J. Atmos. Solar-Terr. Phys., v.67, p.1259-1268, 2005.
Sorokin V.M., Yaschenko A.K., Hayakawa M. Formation mechanism of the lower-ionosphere disturbances by the atmosphere electric current over a seismic region. J. Atmos. Solar-Terr. Phys., v.68, p.1260-1268, 2006.
Sorokin V.M., Yaschenko A.K., Hayakawa M. A perturbation of DC electric field caused by light ion adhesion to aerosols during the growth in seismic-related atmospheric radioactivity. Nat. Haz. Earth System Sci., v.7, p.155-163, 2007.
Zakharenkova I.E., Krankowski A., Shagimuratov I.I. Modification of the low-latitude ionosphere before the 26 December 2004 Indonesian earthquake. Nat. Haz. Earth Syst. Sci., v.6 (5), p.817-823, doi:10.5194/nhess-6-817-2006, 2006.
Zakharenkova I.E., Shagimuratov I.I., Tepenitzina N.Yu., Krankowski A. Anomalous modification of the ionospheric total electron content prior to the 26 September 2005 Peru earthquake. J. of Atmos. and Solar-Terr. Phys., v.70, N 15, p.1919-1928, doi:10.1016/j.jastp.2008.06.003, 2008.
Zolotov O.V., Namgaladze A.A., Zakharenkova I.E., Martynenko O.V., Shagimuratov I.I. Physical interpretation and mathematical simulation of ionospheric precursors of earthquakes at midlatitudes. Geomagnetism andAeronomy, v.52, N 3, p.390-397, doi: 10.1134/S0016793212030152, 2012.
476
