Wave-like perturbations in the ionospheric F2-layer observed after the Ms8.1 Samoa earthquake of September 29, 2009

A study of the critical frequency foF2 variations after the large earthquake (M_s = 8.1) which occurred on 29 September, 2009 in the region of Samoa Islands in the Pacific Ocean is carried out using data of the ionospheric station of Kwajalein. The epicenter of the earthquake was located at about 184 km southwest from Apia (the capital of West Samoa). It was found that wave-like perturbations of foF2 were observed for ∼3 h above the station (located approximately 3560 km northwest from the epicenter). The amplitude of the disturbance was as large as ∼20% of the average magnetic quiet day foF2 values. A comparison of the observed perturbations of foF2 with the ones detected at Stanford ionospheric station after the Alaska earthquake of 28 March 1964 (M_s = 8.4) showed a close similarity of the wave-like perturbations of foF2 in both cases.     

Surface latent heat flux anomalies quasi-synchronous with ionospheric disturbances before the 2007 Pu’er earthquake in China

One of various mechanisms of pre-earthquake lithosphere–atmosphere–ionosphere coupling as possible explanation of the seismo-ionospheric effects is connected with the release of latent heat. Abnormal variations of ionospheric electromagnetic parameters possibly related to the 2007 Ms 6.4 Pu’er earthquake in China were reported. This paper attempts to examine whether there were abnormal changes of surface latent heat flux (SLHF) linked with this pre-earthquake ionospheric disturbances. The spatio-temporal statistical analyzes of multi-years SLHF data from USA NCEP/NCAR Reanalysis Project reveal that local SLHF enhancements appeared 11, 10 and 7 days before the Pu’er earthquake, respectively. As contrasted to the formerly reported local ionospheric Ne enhancement 9 days before the shock observed by DEMETER satellite, it is discovered that the SLHF and Ne anomalies are quasi-synchronous and have good spatial correspondence with the epicentre and the local active faults. This is valuable for understanding the seismogenic coupling processes and for recognizing earthquake anomaly with multiple parameters from integrated Earth observation system.     

Statistical seismo-ionospheric precursors of M7.0+ earthquakes in Circum-Pacific seismic belt by GPS TEC measurements

The Circum-Pacific seismic belt is the region heavily affected by earthquakes in the world. The relationship between earthquake (e.g., the geographic location, occurrence time, magnitude, and focal depth) and ionospheric anomalies in the belt was investigated using 100 M7.0+ earthquakes during 2006–2015. The ground-based GPS measurements and global ionosphere map (GIM) data were used for the analyses of the ionospheric variations preceding the earthquakes. The results indicated that the occurrence rate of total electron content (TEC) anomalies was proportional to the magnitude and inversely proportional to the focal depth to a certain degree, and the occurrence frequency of anomalies had a rising trend with the days getting close to the main shock. The occurrence rate of TEC anomalies in the Southern hemisphere was larger than that in the Northern hemisphere. Besides, the spatial characteristics of TEC anomalies showed that the anomalies in low-middle latitudes did not coincide with the epicenter, sometimes the anomalies were also observed in the corresponding conjugated region. However, the TEC anomalies in the high latitude usually appeared around the epicenter and within the seismogenic zone while no TEC anomalies appeared in the conjugated area. These results may have potential applications to the earthquake prediction in the Circum-Pacific seismic belt.     

Statistical analysis of the ionospheric ion density recorded by DEMETER in the epicenter areas of earthquakes as well as in their magnetically conjugate point areas

Results of a statistical variation of total ion density observed in the vicinity of epicenters as well as around magnetically conjugated points of earthquakes are presented in this paper. Two data sets are used: the ion density measured by DEMETER during about 6.5?years and the list of strong earthquakes (M_W?≥?4.8) occurring globally during this period (14,764 earthquakes in total). First of all, ionospheric perturbations with 23–120?s observation time corresponding to spatial scales of 160–840?km are automatically detected by a software (64,287 anomalies in total). Second, it is checked if a perturbation could be associated either with the epicenter of an earthquake or with its magnetically conjugated point (distance?<?1500?km and time?<?15?days before the earthquake). The index Kp?<?3 is also considered in order to reduce the effect of the geomagnetic activity on the ionosphere during this period. The results show that it is possible to detect variations of the ionospheric parameters above the epicenter areas as well as above their conjugated points. About one third of the earthquakes are detected with ionospheric influence on both sides of the Earth. There is a trend showing that the perturbation length increases as the magnitude of the detected EQs but it is more obvious for large magnitude. The probability that a perturbation appears is higher on the day of the earthquake and then gradually decreases when the time before the earthquake increases. The spatial distribution of perturbations shows that the probability of perturbations appearing southeast of the epicenter before an earthquake is a little bit higher and that there is an obvious trend because perturbations appear west of the conjugated point of an earthquake.     

Study of the 2013 Lushan M7.0 earthquake coseismic ionospheric disturbances

On April 20, 2013, an earthquake of M7.0 occurred in Lushan, Sichuan province, China. This paper investigates the coseismic ionospheric anomalies using GPS (Global Positioning System) data from 23 reference stations in Sichuan province that are a part of the Crustal Movement Observation Network of China (CMONOC). The recorded results show that a clear ionospheric anomaly occurred within 15 min after the earthquake near the epicenter, and the occurrence time of the anomalies recorded by various stations is related to the distance from the epicenter. The maximum anomaly is 0.25 TECu, with a 2 min duration and the distance of the recording station to the epicenter is 83 km. Acoustic waves generated by the crustal vertical movement during the earthquake propagate up to the height of the ionosphere lead to the ionospheric anomaly, and the propagation speed of the acoustic wave is calculated as 0.72 ± 0.04 km/s based on the propagation time and propagation distance, consistent with the average speed of sound waves within a 0–450 km atmospheric height.     

Magnetospheric disturbances associated with the 13 December 2006 solar flare and their ionospheric effects over North-East Asia

We present an observational study of magnetospheric and ionospheric disturbances during the December 2006 intense magnetic storm associated with the 4В/Х3.4 class solar flare. To perform the study we utilize the ground data from North–East Asian ionospheric and magnetic observatories (60–72°N, 88–152°E) and in situ measurements from LANL, GOES, Geotail and ACE satellites. The comparative analysis of ionospheric, magnetospheric and heliospheric disturbances shows that the interaction of the magnetosphere with heavily compressed solar wind and interplanetary magnetic field caused the initial phase of the magnetic storm. It was accompanied by the intense sporadic E and F2 layers and the total black-out in the nocturnal subauroral ionosphere. During the storm main phase, LANL-97A, LANL 1994_084, LANL 1989-046 and GOES_11 satellites registered a compression of the dayside magnetosphere up to their orbits. In the morning–noon sector the compression was accompanied by an absence of reflections from ionosphere over subauroral ionospheric station Zhigansk (66.8°N, 123.3°E), and a drastic decrease in the F2 layer critical frequency (foF2) up to 54% of the quite one over subauroral Yakutsk station (62°N, 129.7°E). At the end of the main phase, these stations registered a sharp foF2 increase in the afternoon sector. At Yakutsk the peak foF2 was 1.9 time higher than the undisturbed one. The mentioned ionospheric disturbances occurred simultaneously with changes in the temperature, density and temperature anisotropy of particles at geosynchronous orbit, registered by the LANL-97A satellite nearby the meridian of ionospheric and magnetic measurements. The whole complex of disturbances may be caused by radial displacement of the main magnetospheric domains (magnetopause, cusp/cleft, plasma sheet) with respect to the observation points, caused by changes in the solar wind dynamic pressure, the field of magnetospheric convection, and rotation of the Earth.     

Ionospheric disturbances on 8 September, 2010: Was it connected with the incoming moderate Chongqing earthquake?

Earthquake prediction stimulates the searches for a correlation between seismic activity and ionospheric anomalies. Contrary to common focuses on strong earthquakes, we report the ionospheric disturbances, 2 days before a moderate Ms = 4.7 Chongqing earthquake (29.4°N, 105.5°E, depth = 7.0 km, occurred at 21:21 LT, 10 September, 2010) with the data of ground-based ionosondes and IGS receivers. The data covering the period under the quiet geomagnetic conditions and a geomagnetic storm was analyzed with upper and lower bounds. It is found that there were significant enhancements of foF2 and total electron content (TEC) on the afternoon of 8 September, 2010, with a limited area close to the epicentre, which was different from the feature of ionospheric perturbations triggered by the geomagnetic storm on 15 September. Taking into account the heliogeomagnetical condition, we conclude that the observed ionospheric enhancements were very likely associated with the forthcoming moderate Chongqing earthquake, which implies that the relationship between the amplitudes of ionospheric disturbances and earthquakes is very complicated.     

An annual variation analysis of the ionospheric spatial gradient over a regional area for GNSS applications

An ionospheric spatial gradient represents the ionosphere delay difference between different locations, and its variation over a specific area is important for implementing differential GNSS systems. An estimation method for the ionospheric spatial gradient over a small regional area is proposed. A plate map model is implemented for the direct estimation of the gradients. Nine years of GPS data were processed to figure out the annual variation of the mean gradient at the mid-geomagnetic latitude of 30° N. Gradients along the north–south direction have a mean of 0.65 mm/km and follow solar-cycle variations.     

Ground-based and satellite DC-ULF electric field anomalies around Wenchuan M8.0 earthquake

The ground-based and satellite DC-ULF electric field data were analyzed around Wenchuan M8.0 earthquake on May 12, 2008 in China. The results show that ground electric field anomalies occurred at 3 stations located to the north and south of the epicenter with the amplitude of 3–100 mV/km. The change shapes and their amplitude of ground electric field anomalies are different largely due to their individual underground layer conductivity, water level and so on. The analysis of long time series illustrates that the abnormal geoelectric field started since March 2008. Onboard the DEMETER satellite, the ULF waveforms of electric field were collected and processed by wavelet transform method. The disturbances in the ionosphere were about 3–5 mV/m at a frequency band lower than 0.5 Hz. When the ground and space electric field anomalies were compared, their occurrence time and spatial distribution points are consistent with each other, including the long time anomalies from March 2008 and the short term ones 1–2 days before the Wenchuan earthquake. Finally, the coupling mechanism was discussed.     

地磁场与电离层异常现象及其与地震的关系

利用中国地磁台网与电离层台站资料,总结了大地震前出现的地磁低点位移、地磁日变异常及电离层f0F2(F2层临界频率)异常现象.对比研究了1997年11月8日玛尼7.5级与2001年11月14日昆仑山口西8.1级地震前磁场与电离层异常分布及特征.结果显示,两次巨大地震前磁场与电离层短临异常时空分布特征有较好的一致性,震中周围出现日变异常、拉萨台出现电离层f0F2明显异常;震前约1个月出现地磁低点位移,其突变分界线通过震中地区.      

Ionospheric anomalies of local earthquakes detected by GPS TEC measurements using data from Tashkent and Kitab stations

This paper reports the ionospheric anomalies observed during strong local earthquakes (M?5.0) which occurred mostly in and around Uzbekistan in seismically active zones, during years 2006 to 2009 within approximately 1000 km distance from the observing GPS stations located in Tashkent and Kitab, Uzbekistan. The solar and geomagnetic conditions were quiet during occurrence of the selected strong earthquakes. We produce Total Electron Content (TEC) time series over both sites and apply them to detect anomalous TEC signals preceding or accompanying the local earthquakes. The results show anomalous increase or decrease of TEC before or during the earthquakes. In general the anomalies occurred 1–7 days before the earthquakes as ionospheric electromagnetic precursors. To identify the anomalous values of TEC we calculated differential TEC (dTEC). dTEC is obtained by subtracting monthly averaged diurnal vTEC from the values of observed vTEC at each epoch. This procedure removes normal diurnal variations of vTEC. The present results are in good agreement with the previous observations on ionospheric earthquake precursors reported by various researchers.     

F-region ionospheric parameters observed in the equatorial and low latitude regions during medium solar activity in the Brazilian sector and comparison with the IRI-2007 model results

The ionospheric sounding observations using the Canadian Advanced Digital Ionosondes (CADIs) operational at Palmas (PAL; 10.2°S, 48.2°W; dip latitude 6.6°S; a near-equatorial station), and São José dos Campos (SJC, 23.2°S, 45.9°W; dip latitude 17.6°S; a low-latitude station located under the southern crest of the equatorial ionospheric anomaly), Brazil, are analyzed during the different seasons viz., winter (June and July 2003), spring (September and October 2003), summer (December 2003 and January 2004), and fall (March and April 2004). The period used has medium solar activity (sunspot number between 77.4 and 39.3). The seasonal mean variations (using only geomagnetically quiet days) of the ionospheric parameters foF2 (critical frequency of the F-region), hpF2 (virtual height at 0.834 foF2; considered to be close to hmF2 (peak height of the F-region)), and h’F (minimum virtual height of the F-region) are calculated and compared between PAL and SJC. The prominent differences between PAL and SJC are as follows: h’F variations show strong post-sunset enhancement at PAL during the seasons of spring, summer, and fall; hpF2 variations show pre-sunrise uplifting of the F-layer at both stations during all the seasons and the hpF2 values during the daytime are lower at SJC compared with PAL during all the seasons; the foF2 variations show mid-day bite-out at PAL during all the seasons and SJC shows strong equatorial ionospheric anomaly during summer and fall seasons. Also, the seasonal variations of the ionospheric parameters foF2 and hpF2 (with ±1 standard deviation) observed at PAL and SJC are compared with the IRI-2007 model results of foF2 and hmF2. In addition, variations of the foF2 and hpF2 observed at SJC are compared with the IRI-2001 model results of foF2 and hmF2. It should be pointed out that the ionospheric parameter hpF2 is much easier to obtain using computer program developed at UNIVAP compared with hmF2 (using POLAN program). During the daytime due to underlying ionization hpF2 estimated is higher (approximately 50 km) than the true peak height hmF2. During the nighttime hpF2 is fairly close to hmF2. The comparison between the foF2 variations observed at PAL and SJC with the IRI-2007 model results shows a fairly good agreement during all the seasons. However, the comparison between the hpF2 variations observed at PAL and SJC with the hmF2 variations with the IRI-2007 model results shows: (1) a fairly good agreement during the nighttime in all the seasons; (2) the model results do not show the pre-sunrise uplifting of the F-layer at PAL and SJC in any season; (3) the model results do not show the post-sunset uplifting of the F-layer at PAL; (4) considering that, in general, hpF2 is higher than hmF2 during the daytime by about 50 km, the model results are in good agreement at PAL and SJC during all the seasons except summer at SJC, when large discrepancies in the observed hpF2 and modeled hmF2 are observed. Also, it has been observed that, in general, hmF2 values for SJC calculated using IRI-2001 are higher than IRI-2007 during the daytime in winter, summer, and fall. However, hmF2 values for SJC calculated using IRI-2001, are lower than IRI-2007 during the nighttime in spring.     

Validation of the STORM model in IRI-2001 at a high latitude station

The International Reference Ionosphere IRI-2001 model contains geomagnetic activity dependence based on an empirical storm time ionospheric correction (STORM model). An extensive validation of the STORM model for the middle latitude region has been performed. In this paper the ability of the STORM model to predict foF2 values at high latitudes is analyzed. For this, ionosonde data obtained at Base Gral. San Martin (68.1°S, 293°E) are compared with those obtained by the IRI-2001 model with or without storm correction during four geomagnetic storms that occurred in 2000 (Rz₁₂ = 117) and 2001 (Rz₁₂ = 111). The results show that predicted values with the STORM model follow the behaviour of foF2 experimental data better than without the STORM model. The relative deviation between measured and predicted foF2 reaches values of up to 24% and 43% with and without the STORM model in IRI-2001, during the main phase of the storms. In order to explain increases of electron density that occurred prior to the storm onset and also decreases of electron density observed during the first part of the recovery of the storm, possible physical mechanisms are discussed.     

Stability of the seasonal variations in diurnal and semidiurnal components of mid-latitude F2 layer parameters

We report work utilizing 15-min resolution ionospheric data obtained with DPS-4 digisonde in 2003–2011 to study the seasonal variations in amplitudes and phases of the most powerful spectral components of the F2 layer critical frequency (foF2) and peak height (hmF2) fluctuations over Irkutsk (52.5°N, 104.0°E). We show that fluctuations of both parameters contain quasi-harmonic components with periods of T_n = 24/n h (n = 1–7). The number of distinct spectral peaks varies from 3 in summer to 7 in winter. Amplitude and phase characteristics of the diurnal (n = 1) and semidiurnal (n = 2) components is studied using the data sets extracted from the original data sets with band-pass filter. It has been found that the amplitudes of diurnal/semidiurnal foF2 and diurnal hmF2 components are maximum in winter and minimum in summer. Amplitudes of the diurnal components vary gradually; those of the foF2 semidiurnal one, abruptly, thus forming a narrow winter maximum in November–January. The phase (local time of maximum) of the diurnal foF2 component increases gradually by 4–6 h from winter to summer. The phase of the semidiurnal foF2 component is nearly stable in winter/summer and sharply decreases (increases) by 2–3 h near the spring (autumn) equinox. The phase of the diurnal component of hmF2 (local time of minimum) varies slightly between 1130 and 1300 LT; that of the semidiurnal one decreases (increases) by 4–6 h from January to March (from September to November). The results obtained show that the main features of seasonal variations in the diurnal and semidiurnal components of the mid-latitude F2 layer parameters recur consistently during the solar activity growth and decline phases.     

Variations of ionospheric total electron content before three strong earthquakes in the Qinghai-Tibet region

We investigate the ionospheric total electron content (TEC) anomalies occurred in the Qinghai-Tibet region before three large earthquakes (M > 7.0). The temporal and spatial TEC variations were used to detect the ionospheric possible precursors of these earthquakes. We identified two TEC enhancements in the afternoon local time 9 days and 2–3 days before each earthquake, between which a TEC decrement occurred 3–6 days before earthquakes. These anomalies happened in the area of about 30° in latitude and the maximum is localized equatorward from the epicenters. These TEC anomalies can be found in all three earthquakes regardless the geomagnetic conditions. The features of these anomalies have something in common and may have differences from those caused by geomagnetic storms. Our results suggest that these ionospheric TEC perturbations may be precursors of the large earthquakes.     

GPS derived TEC and foF2 variability at an equatorial station and the performance of IRI-model

The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.     

Modeling medium-scale TEC structures,observed by Belgian GPS receivers network

GALOCAD project “Development of a Galileo Local Component for the nowcasting and forecasting of atmospheric disturbances affecting the integrity of high precision Galileo applications” aims to perform a detailed study on ionospheric small- and medium-scale structures and to assess the influence of these structures on the reliability of Galileo precise positioning applications. GPS-derived TEC (total electron content) is obtained from the Belgium Dense Network (BDN), consisting of 67 permanent GPS stations. An empirical 3-D model is developed for studying these ionospheric structures. The model, named LLT model, described temporal variations of TEC in latitude/longitude frame (46°, 52°)N and (−1°, 11°)E. The spatial variations of TEC are modeled by Tchebishev base functions, while the temporal variations are described by a trigonometric basis. To fit the model to the data, the observed area is divided into bins with (1° × 1°) geographic scale and 6 min on time axis. LLT model is made flexible, with varying number of coefficients along each axis. This allows different degree of smoothing, which is the key element of the present approach. Model runs with higher number of coefficients, capturing in details medium-scale TEC structures are subtracted from results obtained with smaller number of coefficients; the latter represent the background ionosphere. The residual structures are localized and followed as they travel across the observed area. In this way, the size, velocity, and direction of the irregular structures are obtained.     

A nonlinear background removal method for seismo-ionospheric anomaly analysis under a complex solar activity scenario: A case study of the M9.0 Tohoku earthquake

A precise determination of ionospheric total electron content (TEC) anomaly variations that are likely associated with large earthquakes as observed by global positioning system (GPS) requires the elimination of the ionospheric effect from irregular solar electromagnetic radiation. In particular, revealing the seismo-ionospheric anomalies when earthquakes occurred during periods of high solar activity is of utmost importance. To overcome this constraint, a multiresolution time series processing technique based on wavelet transform applicable to global ionosphere map (GIM) TEC data was used to remove the nonlinear effect from solar radiation for the earthquake that struck Tohoku, Japan, on 11 March, 2011. As a result, it was found that the extracted TEC have a good correlation with the measured solar extreme ultraviolet flux in 26–34 nm (EUV_26–34) and the 10.7 cm solar radio flux (F_10.7). After removing the influence of solar radiation origin in GIM TEC, the analysis results show that the TEC around the forthcoming epicenter and its conjugate were significantly enhanced in the afternoon period of 8 March 2011, 3 days before the earthquake. The spatial distributions of the TEC anomalous and extreme enhancements indicate that the earthquake preparation process had brought with a TEC anomaly area of size approximately 1650 and 5700 km in the latitudinal and longitudinal directions, respectively.     

Ionospheric disturbances under low solar activity conditions

The paper is focused on ionospheric response to occasional magnetic disturbances above selected ionospheric stations located at middle latitudes of the Northern and Southern Hemisphere under extremely low solar activity conditions of 2007–2009. We analyzed changes in the F2 layer critical frequency foF2 and the F2 layer peak height hmF2 against 27-days running mean obtained for different longitudinal sectors of both hemispheres for the initial, main and recovery phases of selected magnetic disturbances. Our analysis showed that the effects on the middle latitude ionosphere of weak-to-moderate CIR-related magnetic storms, which mostly occur around solar minimum period, could be comparable with the effects of strong magnetic storms. In general, both positive and negative deviations of foF2 and hmF2 have been observed independent on season and location. However positive effects on foF2 prevailed and were more significant. Observations of stormy ionosphere also showed large departures from the climatology within storm recovery phase, which are comparable with those usually observed during the storm main phase. The IRI STORM model gave no reliable corrections of foF2 for analyzed events.     

Response of the ionosphere to super geomagnetic storms: Observations and modeling

The relative importance of the main drivers of positive ionospheric storms at low-mid latitudes is studied using observations and modeling for the first time. In response to a rare super double geomagnetic storm during 07–11 November 2004, the low-mid latitude (17°–48°N geomag. lat.) ionosphere produced positive ionospheric storms in peak electron density (NmF2) in Japan longitudes (≈125°–145°E) on the day of main phase (MP1) onset (06:30 LT) and negative ionospheric storms in American longitudes (≈65°–120°W) on the following day of MP1 onset (13:00–16:00 LT). The relative effects of the main drivers of the positive ionospheric storms (penetrating daytime eastward electric field, and direct and indirect effects of equatorward neutral wind) are studied using the Sheffield University Plasmasphere Ionosphere Model (SUPIM). The model results show that the penetrating daytime (morning–noon) eastward electric field shifts the equatorial ionisation anomaly crests in NmF2 and TEC (total electron content) to higher than normal latitudes and reduces their values at latitudes at and within the anomaly crests while the direct effects of the equatorward wind (that reduce poleward plasma flow and raise the ionosphere to high altitudes of reduced chemical loss) combined with daytime production of ionisation increase NmF2 and TEC at latitudes poleward of the equatorial region; the later effects can be major causes of positive ionospheric storms at mid latitudes. The downwelling (indirect) effect of the wind increases NmF2 and TEC at low latitudes while its upwelling (indirect) effect reduces NmF2 and TEC at mid latitudes. The net effect of all main drivers is positive ionospheric storms at low-mid latitudes in Japan longitude, which qualitatively agrees with the observations.