VLF / LF Amplitude Perturbations before Tuscany Earthquakes , 2013

The US Navy VLF/LF Transmitter’s NSY signal (45.9 kHz) transmitted from Niscemi, Sicily, Italy, and received at the Kiel Long Wave Monitor, Germany, was analyzed for the period of two months, May and June (EQ-month) of 2013. There were 12 earthquakes of magnitude greater than 4 that hit Italy in these two months, of which the earthquake of 21 June having magnitude of 5.2 and a shallow focal depth of 5 km was the major one. We studied the earthquake of 21 of June 2013, which struck Tuscany, Central Italy, (44.1713 ̊N and 10.2082 ̊E) at 10:33 UT, and also analyzed the effects of this earthquake on the sub-ionospheric VLF/LF signals. In addition, we also studied another earthquake, of magnitude 4.9, which hit the same place at 14:40 UT on 30 of June and had shallow focal depth of 10 km. We assessed the data using terminator time (TT) method and night time fluctuation method and found unusual changes in VLF/LF amplitudes/phases. Analysis of trend, night time dispersion, and night time fluctuation was also carried and several anomalies were detected. Most ionospheric perturbations in these parameters were found in the month of June, from few days to few weeks prior to the earthquakes. Moreover, we filtered the possible effects due to geomagnetic storms, auroras, and solar activities using parameters like Dst index, AE index, and Kp index for analyzing the geomagnetic effects, and Bz (sigma) index, sunspot numbers, and solar index F10.7 for analyzing the solar activities for the confirmation of anomalies as precursors.


INTRODUCTION
Earthquakes still exist as a threatening disaster.Being able to predict earthquakes will certainly help prevent substantial loss of life and property.Perhaps, Very Low Frequency (VLF/LF) signals can be estab-Natural Science lished as a tool to predict earthquake.VLF/LF band (3 -30 kHz) falls under radio spectrum and is widely known for its application in navigation and communication.Due to low attenuation and long wavelength, VLF/LF signals-transmitted from different VLF/LF transmitters established for research or other application purposes-can travel long distance by undergoing consecutive reflections between lower ionosphere and surface of the Earth.This system of ionosphere and earth that guides the VLF/LF wave is called Earth-ionosphere waveguide.Many research works present that perturbations in VLF/LF amplitudes or phases have been observed prior to earthquake day [1][2][3][4].However, the results are often controversial for there are many other factors that might influence the VLF/LF waves.For instance, the ionospheric layer is highly influenced by geomagnetic winds, solar flares, cosmic rays and gamma-rays from the space, coronal mass ejection, so on and so forth, which can also bring about variations in VLF/LF wave propagation.
Anomalous shift in Terminator times, during sunrise and sunset, have been reported few days prior to earthquake day majorly in the works of [3,[5][6][7][8][9][10][11].Several studies show shifts in D-layer formation and D-layer disappearance time prior to seismic activities [12].Extensive use of the methods of trend, dispersion and night time fluctuation have been done to analyze the VLF/LF signals, where the anomaly is characterized by the simultaneous decrease in the trend and the increases in dispersion and night time fluctuation [13][14][15][16][17]. Thus, changes observed in terminator times, trend, dispersion and night time fluctuation have been considered as precursor to earthquakes [2, 5, 13-15, 18, 19].The analysis of the Central Italy Earthquake of 2016 was performed in our previous work for the same Transmitter-Receiver Great Circle Path (TRGCP) where significant anomalies in terminator times and night time data were found [14].This work for Tuscany earthquake is an extension of previous work for earthquakes in Italy.
We carried out analysis of VLF/LF signals for the month of May and June (EQ month) 2013 to see if any perturbations occurred in those signals before the Tuscany earthquakes of 21 st of June of magnitude 5.2 M and 30 th of June of magnitude 4.9 M. The earthquakes of magnitude <5.5 are considered incapable of perturbing the ionosphere.A long study had observed changes in VLF/LF signals only due to earthquakes that had magnitude at least greater than 5.5 M [17].However, in our Transmitter-Receiver Great Circle Path (TRGCP), we observed significant perturbations in VLF/LF amplitude despite having EQ of magnitude <5.5 M which makes our work an important one.Reference [16] had carried out similar research in finding the relation between earthquake and ionospheric perturbations using very long propagation path of length 11,400 km.The study was done on the earthquakes in the Indonesia region and anomalies due to an earthquake of magnitude 5.8, which had a depth of 53 km and epicenter ~45 km away from the propagation path, were observed.Considerable deviation in terminator times and changes in the trend and night time fluctuation were reported.In contrast to [16], our propagation path was short with path distance ~1937 km; the propagation path is said to be short when the propagation path is <2500 -3000 km [20].The earthquake on 21 st June 2013 (hereafter EQ-1) had shallow depth of 5 km while that of 30 th June 2013 (hereafter EQ-2) was 10 km and the propagation path almost crossed over the epicenters.Thus, we assume shallow depth of the earthquake to be the reason that the seismic activities could cause disturbances in the lower ionosphere.Similarly, these disturbances were received by the signal before they got attenuated due to the short propagation path length and closeness of epicenters to the propagation path; and hence, even the earthquakes of magnitude 5.2 and 4.9 could have caused considerable perturbations.

TRANSMITTER AND RECEIVER
The VLF/LF recording system of Kiel Long Wave monitor consists of a receiving antenna, a receiver antenna (amplifier), and a computer.The setup of Kiel Long Wave monitor consists of a Thin-Client computer along with identical ferrite loop antennae.The antenna runs at a constant location in a constant direction and nulls at 50 and 230 degrees.The receivers are broadband of 10 -90 kHz running 90 Fast Fourier Transforms (FFT) per second.Each datum used in this case was a one minute spectrum that was an average from approximately 5000 raw spectra.The Kiel Long wave Monitor records the VLF/LF spectrum between 10 and 96 kHz by means of two receivers running at 6 km of distance: the H-field probe (54˚23'46"N, 10˚03'13"E) and the E-field (54˚22'24"N, 10˚08'15"E).Figure 1 gives the outline of the procedures carried out in Kiel Long Wave Monitor.A total of 12 earthquakes greater than magnitude 4 in Natural Science

METHODS OF ANALYSIS
Broadly, two method of analysis used are as follows: 1) Terminator time (TT) method; 2) Night time Fluctuation (NF) method.Natural Science The terminator time method takes into account two parameters: Sunrise terminator time (SRT) and Sunset terminator time (SST).We observed one minimum for each morning and evening around sunrise and sunset.We will simply use SRT for sunrise terminator time and SST for sunset terminator time henceforth.SRT is defined as the first minimum which occurs after the sharp weakening of signals around sunrise and SST is defined as the last weakening of the signal before its sharp rise around sunset (marked by arrows in Figure 2) [12].Moreover, we have marked T A to be time from which the signal starts to weaken sharply towards SRT (marked by arrows and dashed vertical blue lines) and T B to be time where signal becomes more or less constant after a period of sharp increase from SST (marked by arrows and dashed vertical green lines).The parameter SRT − T A (SRT minus T A ) is D-layer Preparation time (DLPT), time taken to lower D-layer boundary in the early morning and T B − SST (T B minus SST) is D-layer disappearance time (DLDT), time taken to raise D layer boundary again in the evening [11,12].There are two other minima in between SRT and SST which are fluctuating rapidly and randomly throughout the month, due to solar effects.
There are two methods of nighttime fluctuation method.In method by [4], we calculate the difference between maximum and minimum fluctuations of signal amplitude during nighttime.In our case, the night time taken is 21:15 UT to 2:00 UT.We have plotted this difference for 19 days centered on 21 st June, EQ-1.In method by [13], we have computed various parameters namely trend, nighttime fluctuation and dispersion based on local nighttime.Here, the procedure to calculate ±2σ lines and to calculate trend and night time Figure 2. Variation of amplitude of VLF/LF signal on 16 th of June.The first minimum around sunrise is sunrise terminator time (SRT) and the second minimum around sunset is sunset terminator time (SST).D-layer preparation time (DLPT) is the time period between the blue and red dashed lines and D-layer disappearance time (DLDT) is the time period between the pink and green dashed lines.T A and T B are the points where blue and green lines intersect the graph respectively.Natural Science fluctuation is escaped as it has been explained thoroughly in our previous work [14].However, we have mentioned the formulas that we used and discussed about new a parameter being used in this paper.
Considering initial terminator times as T i and final terminator time as T f , trend, night time dispersion, and night time fluctuations can be mathematically calculated as: ( ) Nighttime Fluctuation NF ; 0 where, = deviation of instantaneous amplitude from monthly average (at the same time). .Night time dispersion considers deviation from night time average rather than deviation from daily average.We couldn't calculate the daily averages because we had poor daytime data (owing to transmitter shutdowns).So, we simply used the night time data we had in order to find night time averages and eventually calculated the deviations from those night time averages.This new parameter was then called night time dispersion, denoted as dispersion*.This parameter would measure how the night time data is fluctuating from the average of the particular night.We expect this parameter, dispersion*, to decrease simultaneously with anomalies in trend and night time fluctuation because significant decrease in amplitude has been found to be signature of seismic activities [2,17].
The terminator times and nighttime amplitudes are supposed to be more or less constant on normal days, days unaffected by solar and geomagnetic activities, lightning, CG Flashes, LAIC and so on.The data points must lie closely around the line of best fit.Standard deviation, represented by σ, is a measure of how far the data points are deviated over whole set of data, implying a state where all data points vary between ±σ lines.However, if data crosses ±σ lines, it indicates that the ionosphere is perturbed to some extent by various events like solar and geomagnetic activities.On crossing ±2σ lines, it indicates a significant perturbation caused in ionosphere.It is believed that the deposition of ground seismic energy or radioactive decay of radon and other isotopes into the ionosphere can cause excessive ionization which may enhance the conductivities of these regions much before the seismic event [21].The heating by radon decay may also produce slow expansion of atmosphere which can perturb the ionosphere.Moreover, fracturing of tectonic plate emits ions and electron to affect the conductivity of the ionosphere.However, real physical explanation is least understood due to less progress in lithosphere-ionosphere coupling.This way, the ionosphere is perturbed strongly from few days to few weeks before the earthquake day causing the data of various parameters to shift crossing ±2σ lines.

Terminator Time (TT) Method
This is the first method of data analysis in which we plotted the graphs of sunrise terminator time Natural Science (SRT) and sunset terminator time (SST) against days for the month of May and June 2013.The each graph possess three colored lines, the upper green line with small triangles which represents +2σ line, the middle red line with small circles indicating the trend line (the line of best fit), and lastly the lower blue line with small triangles which symbolizes −2σ line.The values that cross the green line (+2σ line) or blue line (−2σ line) are considered "anomalous values".
In the month of May, the month prior to earthquake month, we found SRT well within ±2σ lines (the green and blue lines) as shown in the Figure 3(a).Similarly, it can be seen in Figure 3(b) that SST shows shift on 3 rd of May 2013 as it crosses −2σ line (the blue line).D-layer preparation time (DLPT)does not show abnormality but D-layer disappearance time (DLDT) was found to cross 2σ line (lower blue line) on 3 rd of May as presented in Figure 4.This anomaly in DLDT on 3 rd May couldn't be the precursor to the earthquake of June 21 st and 30 th which were still 48 -57 days away.However, the shift in 3 rd of May could be the precursor to the earthquake of 4 th of May of magnitude 4.1 (and even due to the effects of other earthquakes like that of 8 May having magnitude 4.4).It can't be ignored that the earthquakes given in Table 1 have very shallow depth varying from 4 km to 10 km; and Tuscany is an area that falls almost within the propagation path of the signals i.e., the VLF/LF NSY signals from Niscemi nearly passes over Tuscany on the way to Kiel, Germany, where the signal are received.Thus, the shallow depth could be the reason that the weak earthquakes of magnitude just as great as 4 could cause disturbances and were easily detected.And the solar activities were also relatively high on 3 rd May as shown by the increase in solar  F10.7 index.The effects of solar and geomagnetic activities will be dealt again in the section 5 in more detailed way.Thus, it is also likely that the anomaly was caused due to the effect of earthquakes (occurring in the month of May) as well as that of solar activities.
In the month of June (EQ-month), greater perturbations in the VLF/LF amplitudes were observed.The sunrise terminator time (SRT) shifted toward the night as it crossed −2σ line (lower blue line) on 23 rd of June, 7 days prior toEQ-2, as seen in Figure 5(a).The SRT greatly shift towards the night on 7 th , 8 th and 18 th June (days before EQ-1) though does not cross the 2σ lines.The sunset terminator time SST crossed the ±2σ lines on 1 st , 7 th , 9 th , 11 th , 13 th , 17 th , 18 th , 21 st , 22 nd , 23 rd , 24 th and 27 th of June, as presented in Figure 5(b).However, the entire anomalous shifts might not have been the effect of just EQ-1 and EQ-2.The shift in SST on 1 st , 7 th , 9 th , 11 th , and 13 th might have occurred due to earthquake of 16 June 2013 which had magnitude of 4.1 and epicenter ~98 km away from the propagation path.It had shallow depth as well i.e., <10 km.And the anomalies on 17 th and 18 th June might have been due to the seismic activities that might have disturbed the lower ionosphere before the tremor was felt on 21 st June (5.2M, EQ-1).Similarly, the shift of SST on 21 st and 22 nd could be precursor to the earthquake of 23 rd of June (4.7 M and 4.0M) and that of 23 rd , 24 th , and 27 th might be precursor to earthquake of 30 th June (4.7 M and 4.9 M, EQ-2).These results are consistent with our claim that changes in terminator times are seen before earthquakes.Focusing on the main earthquakes under study, we find that anomalies were seen from 3 days (18 th June) to 4 days (17 th June) prior to EQ-1 and 7 days (23 rd June) to 3 days (27 th June) prior to EQ-2.However, these aren't filtered data.Natural Science However, the conclusion section provides a table where the summary of probable precursors after casting out the days that were influenced by solar or geomagnetic activities.
It was unusual not to get any anomalies before the earthquake of 21 st May of magnitude 4.5 (greater than that of 4 th and 5 th May or 16 th June).It was also surprising to find such a great shift on 27 th of June.On evaluating the geomagnetic and solar activities, no significant changes in the parameters were found on 27 th of June except for Bz (sigma), which seemed relatively higher.Thus, the anomaly might have resulted from combined effect of the earthquakes-of magnitude 4.7 and 4.9 that struck Tuscany on 30 th June-and solar activities.Most of our results agreed that anomalous shift in terminator times could be considered possible precursors of the earthquakes.From the Figure 6(a), we can observe that D-layer preparation time (DLPT) shows anomaly 5 days before EQ-2 as it just crosses the +2σ lines (upper green lines in graph) on 25 th of June while D-layer disappearance time (DLDT) shows anomaly 11 days before of EQ-1 and 2 days before EQ-2 for it just touches the ±2σ on 10 th and 28 th of June 2013.However, we didn't observe any anomalous fluctuations in value of DLPT before EQ-1.
On the basis of these results for the month of June, it is clear that perturbations were seen in terminator time before the Tuscany EQ-1 and EQ-2.Thus, perturbations in the terminator times are likely to be Natural Science observed before the earthquake which can help in predicting it.Similarly, DLPT and DLDT also show perturbations before earthquake and can be considered as a precursor to the earthquake.

Nighttime Fluctuation (NF) Method
In Figure 7, we plot the absolute difference between maximum amplitude signal and minimum amplitude signal of particular night over 19 days centered to EQ-1.We computed data based on each night's data as in [4].We noted that the peak occurs anomalously 2 days before the event implying that the fluctuation becomes very high.In Figure 7(a), the dashed line represents the average nighttime amplitude fluctuation and the short dotted line represents the +3σ line (where σ is the standard deviation), computed on each night.The red vertical dashed line shows the day of earthquake (EQ-1, 5.2 M).Clearly, the graphs in Figure 7(a) shows that the fluctuations peak crossed +3σ line on 19th of June, two days prior to EQ-1.
Interestingly, the day 19th June was international quiet day (see Section 6) i.e. geomagnetic variations were quite low, which strengthens the fact that anomalies were there as a precursor.This is consistent with previous findings by the authors [4].However no such fluctuation peak was observed before the EQ-2 implying less perturbation during nighttime.

Trend, Night Time Fluctuation, and Night Time Dispersion
In night time fluctuation method, we estimated the values of trend and night time fluctuation to show Natural Science that both of these parameters show unusual behavior simultaneously prior to the earthquake as described in [13,[15][16][17] and we expect the night time fluctuation to cross the 2σ lines (i.e. to show anomaly) with simultaneous decrease in trend and increase in dispersion.However, in case of the dispersion, instead of taking the daily averages we only took the night time averages and determined the value for night time dispersion (dispersion*).Therefore, in our case, we expect the dispersion* to decrease simultaneously on the days when trend and night time fluctuation cross 2σ lines.These parameters have been explained in Section 3. Figure 8 shows non-normalized trend, night time dispersion, and night time fluctuations, while Figure 9 shows the normalized trend, night time fluctuations and night time dispersion for the month of May and June.For obtaining values for normalized trend, we subtract the average of month from the trend of corresponding days and then, divide it by their standard deviation.In the same way, we compute normalised night time dispersion and normalised night time fluctuation.In normalized graphs, it automatically brings in sigma lines, making analysis easier.
The anomalous values, the values that cross 2σ line (horizontal dashed line), are denoted by red filled columns in corresponding graphs in Figure 8 and Figure 9.The trend value declined and crossed the 2σ line Natural Science  In brief, the night time fluctuation crosses the 2σ lines on 31 st of May, 17 th of June and 25 th of June and simultaneously the trend decreases abruptly even crossing 2σ line.Not only this, the night time dispersion also decreases crossing 2σ lines except on 25 th of June where it only decreases and therefore, the results were as per our expectation.Since EQ-I was not strong in magnitude; it is very unlikely that it perturbed the ionosphere 20 days before, on 31 st May.The Kp index was relatively high on 31 st May which explains the geomagnetic influences as the reason for the perturbation.Similarly, other factors like CG flashes, which have not been analyzed in this paper, could have also led to this unexpected variation.However, the anomaly seen on 17 th of June, which occurred4 days before EQ-1and anomaly seen on 25 th of June, which occurred 5 days prior to EQ-2, could be the precursor of those earthquakes.[2,3,15,19].We determined the values of d(A) for night time, in the same way as mentioned in section 3 (Methods of Analysis), for the period of 15 th -25 th June (Figure 10).The solid horizontal line (y = 0 line) is the average night time amplitude and the values of the night time fluctuations in different times of a day were plotted.While analyzing, we were more concerned with d(A) < 0 because the average night time amplitude were found to decrease around the days of earthquake [14,17].The night time fluctuation decreased 2 -4 days before EQ-1, starting from 17 June through 19 June, as shown by the red dotted eclipses in Figure 10.Thus, changes in night time fluctuations also seems likely to indicate the upcoming seismic activity.

GEOMAGNETIC AND SOLAR ACTIVITIES
In order to verify that the anomalies that we observed were result of seismo-ionospheric event (EQ) and not due to the geomagnetic storms, auroral activities or other geomagnetic activities we plotted the graph of Disturbance Storm-Time (Dst) index, AuroralElectrojet (AE) index and Kp index for the months May and June.We extracted the required data from http://omniweb.gsfc.nasa.gov/form/dx1.htmlandconfirmed them from http://wdc.kugi.kyoto-u.ac.jp/.We then evaluated if the days of anomalies (the blue dashed lines) coincides with the days having high geomagnetic activities with help of Figure 11   no significant fall in Dst index on the day of anomaly (Figure 11).So, it lends support to our assumption that the anomalous variation on 3 rd May was due to various earthquakes occurring in the month of May i.e., 4 th and 8 th of May (Table 1).We also saw that on 31 st May, when the anomalies on trend, night time dispersion and night time fluctuations were seen, value for AE and Dst index was normal but value of Kp index was relatively high.This clarifies that the perturbations on 31 st of May was due to the effect of geomagnetic storms.
Similarly, no rise in AE and Kp index and no fall in Dst index were seen on anomalous days of the month of June as shown in Figure 12, except for 24 th of June where AE index and Kp index have greater value as compared to other anomalous days.Thus, most of the disturbances in VLF signals that caused led to anomalies in the parameters we used for analysis were due the seismic activities and their effect on the ionosphere before the earthquake.
The solar activities like solar flares, coronal mass ejection, solar cycles, etc., can also affect the ionosphere and cause unwanted disturbances in VLF/LF amplitude.Hence, to check whether the observed anomalies were due to different solar activities, we plotted the graphs of Bz (sigma), sunspot numbers, and solar index 10.7 (it provides climatology of solar activity over different solar cycles).
Blue dashes lines passing through the anomalous days have been given in Figure 13 and Figure 14.The terms in blue at the top and bottom of the graph show the corresponding parameters that were found anomalous.We found that Bz was normal for all anomalous day in May but was relatively higher on 31 st Natural Science May, the day in which trend, night time dispersion, and night time fluctuation were unexpectedly anomalous, (Figure 13(a)).The solar index F10.7 was relatively higher on 3 rd May, the day in which anomaly in SST and DLDT was seen (Figure 13(b)).Thus, we doubt that the anomaly on 3 rd might also be affected by solar activities along with seismo-ionospheric effects.
For the month of June, we found that Bz sigma was relatively high on 1 st and 27 th of June among the anomalous days of June as shown in Figure 14(a).The solar index F10.7 and sunspot number was relatively high on 17 th , 18 th , 21 st , 22 nd , and 23 rd of June as presented in Figure 14(b) and Figure 14(c).
Among the large number of anomalies we noted, few anomalies seem to suffer effects of geomagnetic or solar activities.So in these days, the perturbations might have occurred as a combined effect of seismic effects and the geomagnetic and/or solar effects.All the other anomalies (on 3 rd , 7 th , 9 th , 10 th , 11 th , 13 th , 25 th , and 28 th ) were caused by some other factor than geomagnetic or solar activities.This can be confirmed from the Figures 11-14, as the values of all the indices that we used to measure geomagnetic and solar activities were normal in many anomalous days of the month of June.Thus, the perturbations on the rest of the anomalous days can be said to be caused by the various seismicactivities, which influenced the sub-ionospheric layer, before the earthquake was felt.Thus, it can be concluded that various shifts in terminator times and changes in night time fluctuation, trend, night time dispersion, and average night time amplitudes that occur prior to the earthquake are the possible precursors to the earthquake.

10 INTERNATIONAL QUIETEST DAYS (Q1-Q10) AND 5 INTERNATIONAL MOST DISTURBED DAYS (D1-D5)
International disturbed days are the days in which maximum geomagnetic variations takes place.Natural Science Based on the table above, only 7 th of June, among the anomalous days, was international disturbed day.However, among the days after 16 th June, which we are more concerned with, none of the days are international disturbed days.This data also helps to support our assertion that the anomalies on other days were no due to geomagnetic activities.We referred to the site http://www.ga.gov.au/oracle/geomag/iqd_form.jsp for getting the information about international quiet and disturbed days and confirmed from http://wdc.kugi.kyoto-u.ac.jp/qddays/.Following Table 2 shows the corresponding days for the month of May and June, 2013.

CONCLUSIONS
On analyzing the VLF/LF amplitude data of May and June 2013, many unusual perturbations were observed especially in the month of June, the month in which the Tuscany earthquake took place.Following conclusions were drawn from this study.We have taken the considered the anomalous days within 2 weeks' time while making this table."Clear" precursors in Table 3 indicate the anomalies that were filtered out after considering the effects of solar and geomagnetic effects."Unclear" precursors mentioned in table must have one of the following possibilities 1) Combined effect of seismic activities before the earthquake tremor and the particular high geomagnetic or solar activities.2) Complete effect of high geomagnetic or solar activities.3) Complete effect of seismic events.Most of the "unclear" precursors in our case might be combined effect.It is not necessary that high geomagnetic or solar activity always affects the ionosphere severely.However, to be on safe side, we removed all the possible effects by geomagnetic or solar activities to confirm the anomalies to be clear precursors.Moreover, one can notice most of the anomalies are seen in terminator times and most of clear precursors are seen in nighttime fluctuation method.The anomalies before 7 days of earthquake have been considered the precursor to other earthquakes of the month as mentioned in Section 4.1.
Thus, VLF/LF signals were perturbed few days about 2 weeks before the day of earthquake.Hence, these aforementioned anomalies can be considered possible precursors to the earthquake.We assume that the shortness of propagation path, nearness of epicenter to propagation path, and shallow focal depth to be the reasons that the perturbations in VLF signal were easily received and eventually recorded by the Kiel Long Wave Monitor even when the earthquake was not strong enough.

Figure 1 .
Figure 1.Block diagram of kiel long wave monitor where ferrite antennae receives the VLF/LF signals, amplifies, and processes them to obtain ASCII file finally containing minute average data and second average data.It uploads the data into web via SSH.
night time average (night of particular date).Here, integration represents the area under the curve dA over time (dt) which is further divided by the interval between the chosen terminator times (f i T T − ).The parameter night time dispersion is the variation from the night time average of amplitudes ( ) night A

Figure 3 .
Figure 3. (a) Daily variation in sunrise terminator time (SRT) for the month of May (b) Daily variation in sunset terminator time (SST) for the month of May.The upper green line signifies +2σ line, the middle red line represents trend line (line of best fit), and the lower blue line symbolizes the −2σ line.SST shifted on 3 rd of May and crossed −2σ line (lower blue line).

Figure 4 .
Figure 4. (a) Daily variation in D-layer preparation time (DLPT) for the month of May (b) Daily variation in D-layer disappearance time (DLDT) for the month.The upper green line signifies +2σ line, the middle red line represents trend line (line of best fit), and the lower blue line symbolizes the −2σ line.On 3 rd of May, DLDT just crossed 2σ line (upper green line).

Figure 5 .
Figure 5. (a) Daily variation in sunrise terminator times (SRT) for the month of June (b) Daily variations in sunset terminator times (SST) for the month of June.The upper green line signifies +2σ line, the middle red line represents trend line (line of best fit), and the lower blue line symbolizes the −2σ line.The red dashed lines marks the earthquake days of 21 st June and 30 th June.

Figure 6 .
Figure 6.(a) Daily variation in D-layer preparation time (DLPT) for the month of June (b) Daily variation in D-layer disappearance time (DLDT) for the month of June.The upper green line signifies +2σ line, the middle red line represents trend line (line of best fit), and the lower blue line symbolizes the −2σ line.The red dashed lines marks the earthquake days of 21 st June and 30 th June.

Figure 7 .
Figure 7. (a) The amplitude of the nighttime VLF fluctuations (closed squares) signals are plotted as a function of day number.The horizontal dashed line represents the average nighttime fluctuations while the horizontal short dotted line represent the +3σ line.The vertical red dashed line indicates EQ-1.The fluctuation 2 days before the event crossed the +3σ line; (b) The nighttime VLF amplitude fluctuations for 19 days centered on earthquake of 21 st June 2013.The peak appears 2 days before the event.

Figure 8 .
Figure 8. Graphical representations of non-normalized trend, night time fluctuation (NF), and night time dispersion.The horizontal black dashed lines represent the ±2σ lines.The red filled columns represent anomalous values (values crossing the 2σ lines) while the green filled column in dispersion* (8(c)) indicates unusual value relative to the trend and night time fluctuation on 25 th June.The earthquakes on 21 st June and 30 th June are marked by vertical red dashed lines.

Figure 9 .
Figure 9. Graphical representation of normalized trend, night time fluctuation (NF), and dispersion.The horizontal black dashed lines represent the ±2σ lines.The red filled columns represent anomalous values (crossing the 2σ lines) while the green filled column in normalized night time fluctuation, Figure 10(b), indicates unusual value relative to the normalized dispersion and night time fluctuation on 25 th June.The earthquakes on 21 st June and 30 th June are marked by vertical red dashed lines.

Figure 10 .
Figure 10.Fluctuations in night time amplitude during the period of 15 -22 June 2013.The red dotted eclipses shows the condition of d(A) < 0 prior to the earthquake day.The earthquake day has been marked in red letter on the x-axis.The night time fluctuation unusually decreased 2 -4 days prior to EQ-1 of 21 st June.i.e., on 17, 18 and 19 June.
Figure 11.Variations in the (a) AE index, (b) Dst index, and (c) Kp index for the month of May.Blue dashed lines pass through anomalous days of May, 2013.On 3 rd May, SST and DLDT showed anomalies.

Figure 12 .
Figure 12.Variations in the (a) AE index, (b) Dst index, and (c) Kp index for the month of June.Blue dashed lines pass through the anomalous days of June, 2013.The terms in blue represent the parameters that were anomalous on the corresponding days.TND and NFM stand for trend and night time fluctuation method respectively.

Figure 13 .
Figure 13.Variations in the (a) Bz (sigma), (b) solar index F10.7,and (c) sunspot number for the month of May.Blue dashed lines pass through anomalous days of May, 2013.On 3 rd May SST and DLDT show anomalies.

Figure 14 .
Figure 14.Variations in the (a) Bz (sigma), (b) solar index F10.7,and (c) sunspot number for the month of June.Blue dashed lines pass through the anomalous days of June, 2013.The terms in blue represent the parameters that were anomalous on the corresponding days.

Table
Table showing earthquakes of magnitude greater than 4 in Richter scale in the month of May and June.

Table 2 .
Table showing international disturbed and quiet days for the month of May and June.