Unusual Change in Critical Frequency of F 2 Layer during and Prior to Earthquakes

The unusualness in the critical frequency of different layers of earth’s ionosphere is commensurate to be associated with seismic events. We present study of critical frequency of F2 layer (denoted as f0F2) during some major earthquakes in South American region. We use the semi-empirical Barbier’s theorem of air-glow and define a parameter using the critical frequency and virtual height of F2 layer and named it as “F Parameter”. To investigate the variation of this parameter, we consider five large earthquakes in the junction of Nazca plate and South American plate having magnitude greater than M > 6.5 and study the temporal variation of F parameter during these earthquakes. The F Parameter is measured using the ionograms as recorded from the Ionosonde in Jicamarca Radio Observatory (lat. 11.95 ̊S, long 76.87 ̊W) in Chile which lies within the earthquake preparation zones of these five earthquakes. We examine the F Parameter within a span of ±15 days during earthquakes and observed significant change in the evaluated F Parameter in 12 to 3 days prior to the earthquakes. The increment is over +3σ from the normal variation. We also observe significant changes during aftershock events. The solar geomagnetic indices were found to be low which ensures that these anomalies in F Parameter are due to seismic events.


Introduction
This understanding of seismicity is a vast, complicated, anisotropic and multi parametric problem.Predictions of earthquakes from perceived precursory phenomena with accurate specification of time, location and magnitude of a future earthquake with sufficient precision is therefore not very easy.It is well recognized that electromagnetic wave propagation technique through earth-ionosphere waveguide could be an important tool to predict occurrences of seismic hazards.
The mechanism of seismicity can create a significant thermal, mechanical and electrical perturbation in ionospheric layers through the so-called Lithosphere-Atmosphere-Ionosphere Coupling (LAIC).Thus any physical or chemical changes resulting from such a coupling can be used as a precursory tool for seismic hazards.
Ionosphere is divided in mainly three distinct layers, namely, D, E and F. F is also divided in two sub layers, namely, F 1 and F 2 .The behavior of F 2 layer is described in terms of the extent of its departure from that of a hypothetical Chapman layer.The departure can be termed as an "anomaly".Geographical anomaly, Diurnal anomaly, Seasonal anomaly, winter anomaly are some of such anomalies observed usually in the experimental measurement of Critical frequency of F 2 layer.The Critical frequency is the highest frequency which could be reflected from ionosphere.Above this frequency, the radio wave penetrates into the upper ionosphere.The critical frequency of F 2 layer is denoted by f 0 F 2 .When a radio wave is reflected from a perfect reflector, it is likely to be reflected from a point.But in reality the signal is continuously bent or refracted as it travels through ionosphere, since the reflection and penetration both occurs simultaneously till a height is reached where the residual ray is totally reflected.However, though the path is bent gradually, one can conceive that the incident and reflected rays, when extrapolated, meet at a point.This height could be called the virtual reflection height.The virtual height of an ionized atmospheric layer is measured by the time interval between the transmission of a radio signal and the receipt of the return of its echo.
The existence of a phenomenon called "airglow" was probably discovered before 1800.Yntema [1] was the first person to photometrically establish the phenomenon of airglow which he termed as Earthlight.Following a suggestion of Otto Struve, Elvey [2] introduced the name airglow for the first time.From the basic physics of airglow, related chemical kinetics and excitation mechanism and above all, the ionospheric physics and chemistry, it is obvious that there exists a direct relationship, either of complex or of straightforward nature between ionospheric parameters and airglow emissions.The well-known Barbier formula establishes a direct relationship between the airglow intensity and the critical frequency of the F 2 layer.The variable part of this formula deals with the critical frequency and the virtual reflection height of F layer.The airglow intensity variation depends on the charge densities and type of ions which recombines to emit the radiation.As a whole, the critical frequency is associated with the charge density profile of the ionospheric layers.So can be used as a tool to study the variabilities of the charge density profiles.LAIC mechanism predicts anomalous changes of charge density profile before seismic events.Ionospheric disturbance associated with seismic activities has been largely studied since the Great Alaska Earthquake in March 27, 1964 [3]- [30].Correlations of seismicity with anomaly in F-layer were achieved using different methodologies from spectral analysis [31], satellite observation [32], equatorial ionization anomaly [33], etc.The first publications deal with ionospheric characteristics variations as seismic precursors were Antselevich [34] study of the variations of f 0 E parameter before the Tashkent earthquake 1966.The peak electron density in the F 2 layer appears to be one of the most sensitive parameters connected to seismic activity.Spatial and Temporal variation of electron concentration based on topside ionosonde data during seismic events was done by Pulinets et al. [35].Several studies were carried out regarding the f 0 F 2 variation by Gaivoronskaya & Zelenova [36]; Dupuev & Zelenova [37]; Chuo et al. [38]; Pulinets et al. [39] [40]; Liu et al. [41]; Silina et al. [42]; Pulinets & Legen'ka [43]; Rios et al. [44]; Pulinets & Boyarchuk [16]; Hobara & Parrot [45]; Liperovskaya et al. [46].
In this paper, we consider five different earthquakes from South American region near Peru and Chile.We compute a variable parameter from Barbier theorem using the (f 0 F 2 ) and virtual height profile and study this parameter for a duration of ±15 days around those five earthquakes.The plan of the paper is as following: in the next Section, we explain the methodology and data analysis; in Section 3, we present our results; and in Section 4, we draw our conclusion.

Data and Methodologies
In our entire analysis, we gathered the Ionosonde data from Jicamarca station from http://digisonde.com.In Figure 1  We found the all the chosen earthquakes are near the junction of two tectonic plates, the Nazca plate to the West and the South American plate to the East.
The South American Plate is in motion, moving westward away from the Mid-Atlantic Ridge.The eastward-moving and denser Nazca Plate is sub-ducting under the western edge of the South American Plate along the Pacific coast of the continent at a rate of 77 mm per year.Subduction zones such as the South America arc are geologically complex and generate numerous earthquakes from a variety of tectonic processes that cause deformation of the western edge of South America.The Ionosonde observatory near this place is the Jicamarca Radio Observatory (JRO) (latitude 11.95˚ South, longitude 76.87˚ West).The altitude of JRO is about 520 meters above the sea level.The location of the Jicamarca observatory (blue square) and the 5 earthquakes (>6.5) in the South American region (red circle) are shown in Figure 2.
We also calculate the radius of earthquake preparation zone for these chosen earthquakes using Dobrovolsky formula [47] 0.43 10 km where ρ is the radius of the earthquake preparation zone and M is the Richter magnitude of the earthquake.We found that the distance of Jicamarca Radio Observatory from the epicenter of 5 earthquake lies within the preparation zone.
The detail information of the earthquakes is given in Table 1.In Figure 3 The airglow intensity variation can be can mathematically represented as, [ ] ( ) Barbier [49] [50] was the first to establish a semi empirical formula for OI 6300 Å airglow emission which is given by, ( ) ( ) where, H is the scale height in terms of oxygen and was assumed by Barbier himself to be equal to 80 km.We only consider the variable part of this semi empirical formula and modified this equation to We name the variable quantity ( ) ( ) calculate this for all the five earthquakes to check possible correlations.

Result
We compute the F value for all the earthquakes using the f 0 F 2 from the ionograms.Figure 4 shows the variation of f 0 F 2 for a duration of 15 days from 18 September to 2 October 2001.The black curves are the actual values of f 0 F 2 recorded as an interval of 30 minutes.So for each single day, there are 48 data points.The red curves are the average variation of f 0 F 2 for the same time period  There is an enhancement of f 0 F 2 , four days before the earthquake.The peak after the earthquake day is due to the major aftershocks after the main quake.
for which there were no significant seismic event.Therefore the red curve is a basic calibration of the regular variation of f 0 F 2 in a seismically quiet condition.
An earthquake of Richter scale magnitude (M = 7.5) occurred on 25 September, 2005.It is clear from the above Figure that the value of f 0 F 2 increases unusually four days before the earthquake.There is a similar enhancement of f 0 F 2 just after the earthquake.After the main shock there were a series of aftershocks up to 3 to 4 days with average magnitude more than 5.5.The origin of the second peak is possibly due to these aftershocks.We computed F parameter for all the earthquakes under consideration using the formula mentioned above.Figure 5 shows variation of F parameter as a function of time in days for a span of 21 (10 days before and after) days around the earthquake.The thick curve is the average value of F parameter.We calculate the standard deviation from the mean value and plotted the ±σ, ±2σ and ±3σ with the thick dashed, dotted-dashed and dotted curves respectively with the average value.
Figure 5 shows an enhancement of the F parameter six days before the earthquake.The zero of the X-axis represents the day of the earthquake.It is clear that the F parameter increases anomalously with an order of more than 5σ from the average value.The variation of F parameter does not follow the actual f 0 F 2 variation as presented in Figure 4.There is no secondary maximum after the earthquake day.So the entire effect is pre-seismic.
In Figure 6, we present variation of F parameter for the rest of the four earthquakes in a single grid.The for graphs represents the earthquake as   Figure 6 shows a similar enhancement of F parameter for all the cases.For (a), (b) and (d), the F parameter suffers from an enhancement on 3 to 9 days before the earthquake.For the case of 28 October, 2011 (d), the maximum is quite sharp and reaches more than 7σ level.For (a) 1 April, 2014 and (b) 15 August, 2007, the enhancement is not so sharp but crosses 3σ and 2σ level respectively.For the case of (c) 23 June, 2001, the behavior of F parameter is still anomalous but rather different.The F value becomes minimum on the day of the earthquake.Before the earthquake the F value is quite higher than that for the earthquake day but the values just cross the σ level.The value has a pre-seismic maxima on 11 days prior to the maxima but have significant secondary maxima on 2 to 5 days after the earthquake day.There are two possible reasons behind this post-seismic maxima.First, there are a series of aftershocks occurred for this earthquake.Secondly, there was another strong earthquake with Magnitude M = 7.6 occurring on 7 July, 2001.The secondary peak can be due to pre-seismic effects of the second quake as the second main shock occurs within the next 14 days.So we believe the post-seismic shocks are due to the combined effects of these two factors.To check the solar geomagnetic condition during the earthquake and its associated days we plot the geomagnetic k p index for ±7 days around the earthquake day for all the 5 earthquakes.Figure 7 shows variation of kp index for a duration of 15 days.
It is clear from Figure 7 that during all the earthquakes and their surrounding days, the value k p index ranges between 0 to 4.29 which implies geomagnetically quiet condition (kp < 5).So for the signal has no perturbation due to solar geomagnetic activities and the anomaly in the signal is due to the seismic events.

Conclusion
The LAIC mechanism relates properties of apparently distant and distinct components of Earth system science ranging from lithosphere, lower ionospheric D-layer to upper F 2 layer.Starting from the tropospheric thermal anomalies, lower ionospheric electron density variation to upper ionospheric critical frequency modulation, LAIC takes into account a wide range of geochemical and geo-physical phenomena which could be affected simultaneously.In this paper, we focused our study on the ionospheric F 2 layer where we observe the critical frequency variation for a span of three weeks during some strong seismic events in South America region.We compute a parameter (F) which contains the critical frequency of F 2 layer (f 0 F 2 ) and the virtual reflection height (h') from the historic Barbier's airglow equation and study the behavior of this parameter during those earthquake days.We observe significant increase of this parameter on three to nine days prior to those seismic events.We observe the effects of the aftershocks in the direct observation of both f 0 F 2 and F parameter.We also presented the geomagnetic kp indices for all the five earthquakes and found low geomagnetic condition during all the earthquakes, suggesting our assumption of LAIC mechanism.As yet, we have no clear idea of how and why exactly the lithospheric changes percolate into ionospheric changes.However, our study proves that such changes do occur.Study of physical mechanisms behind the LAIC mechanism and implementing acquired knowledge for future earthquake prediction, regular monitoring of such parameters is absolutely essential to achieve our goal.At the same time, we need to increase number of receiving stations so as to compare the anomalies from different points in order to locate the epicenter with accuracy.
, we present schematic diagram of radio wave propagation technique to indicate the reflection height and critical frequency.

Figure 1 .
Figure 1.Schematic diagram of radio wave propagation technique through earth-ionosphere waveguide depicting the concept of critical frequency and virtual reflection height.

Figure 2 .
Figure 2. The location of Jicamarca Ionosonde observatory (blue square) and epicenters of five earthquakes (red circles) under consideration.
, we present a typical ionogram as obtained from the Jicamarca station.The ionogram is a six-dimensional display, with sounding frequency as the abscissa, virtual reflection height (simple conversion of time delay to range assuming propagation at 3 × 10 8 m/sec) as the ordinate, signal amplitude as the dot size, and echo status (Polarization, Doppler shift, and angle of arrival) mapped into 12 available distinct colors.The wave polarizations are shown as two different color groups.Firstly, the green scale, "neutral" colors showing extraordinary polarization and the red scale, "demanding attention" colors showing ordinary polarization.The angle of arrival is shown by different colors (using the "warm" scale for South and the "cold" scale for North) and the Doppler shift is indicated by the color shades.The left side of Figure shows a table of ionospheric characteristics scaled automatically by the ARTIST software [48].

Figure 4 .
Figure 4.The variation of f 0 F 2 as a function of day number around the earthquake day on 25 September, 2005.The black curves are the actual f 0 F 2 variations where the red curves are the average value of f 0 F 2 when there was no earthquake.There was strong earthquake on 25 September, 2005.There is an enhancement of f 0 F 2 , four days before the earthquake.The peak after the earthquake day is due to the major aftershocks after the main quake.
(a) 01/04/2014; (b) 15/08/2007; (c) 23/06/2001 and (d) 28/10/2011.We follow the same convention for calculating the σ and put the different σ level envelop with the same line style we use in Figure 5.The zero represents the earthquake day marked as E.

Figure 5 .
Figure 5.The variation of F parameter as a function of day number for ±10 days around the earthquake day on 25 September, 2005.The thick solid curve is the average value of the F parameter.The dashed, dotted-dashed and dotted curves are for ±σ, ±2σ and ±3σ level respectively added with the average curve.E is the day of the earthquake.F parameter shows an enhancement before the earthquake.