Ionospheric Stations Research Articles

Ionospheric precursor of a destructive earthquake that occurred on April 6, 2009 at L’Aquila (Italy)

On the basis of the 15-min data from a series of ground-based vertical ionospheric sounding stations, a study of variations of the foF2 critical frequency before the strong earthquake (M = 6.3) that occurred on April 6, 2009 at L’Aquila (Italy) was carried out. The earthquake epicenter was located 85 km north-eastward from Rome. Approximately 20 h prior to the earthquake, a well-pronounced statistically significant effect of foF2 increase relative to the average background for magnetically quiet days was observed for almost 1.5 h at the Rome ionospheric station. In this case, at control stations distanced from the earthquake epicenter, no statistically significant deviations of foF2 from the background values were detected during the same observations period. This fact provides grounds for consideration of the foF2 increase observed at Rome station as a possible ionospheric precursor of this earthquake.

Wave-like perturbations in the ionospheric F2-layer observed after the M s8.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.

Peculiarities of the nighttime winter foF2 increase over Irkutsk

An analysis of properties and peculiarities of the nighttime winter foF2 increases (NWI) in the East Siberia is made on data of ionospheric station Irkutsk in the periods 1958–1992 and 2002–2009 and the empirical model of the F2 layer critical frequency under the geomagnetic quiet conditions deduced from these data (model Q-F2). It is revealed, that the NWI is the stable regularity of the quiet ionosphere over Irkutsk. The amplitude of the NWI (the difference between maximum and minimum foF2 values at night hours) is the greatest in December–January and nearly the same at low and middle solar activity. It is a peculiarity of the quiet ionosphere in the East Siberia. Maximum in night foF2 under quiet geomagnetic conditions is observed mainly after midnight (02-04 LT) and is shifted to predawn hours as solar activity increases. At low solar activity the quiet ionosphere at ∼02–04 LT shows the following properties: (a) the fluctuations of foF2 and hmF2 are in the reverse correlation but this dependence is weak; (b) very strong fluctuations of foF2 (|δfoF2| > 30%) occur seldom (∼4% of events) and almost all of them are positive; an example of very strong fluctuations of foF2 up to 60% can be an extreme increase in the foF2 on 19.12.2008; (c) the very strong enhancements of foF2 in the NWI maximum can be observed at the low geomagnetic activity, they occur more often during substorms but very seldom during geomagnetic storms. Possible reasons of these properties of NWI are discussed.

Ionospheric transients observed at mid-latitudes prior to earthquake activity in Central Italy

Abstract. Ionograms from Rome (41.8N, 12.5E) and Sofia (42.4N, 23.2E) ionospheric stations during earthquake (EQ) activity with magnitude (M) between 5 and 6 in Central Italy are analyzed. It is found that several ionospheric disturbances occur in the intermediate E-F region before the EQ shock. In fact, besides sporadic E (Es) layer development (of type h) of short duration (transients), fmin increase, trace gaps near the critical frequencies, and E region trace disappearance are also observed within one to three hours before the EQ shock. Before the EQ shocks we find that the F2 region parameters are practically undisturbed. The only exception is the so-called fork trace that appears mostly near the critical frequency of the F2 region. Acoustic gravity waves (AGW) are suggested as one of the possible sources of transients observed in the ionosphere before the EQ shock.

Evaluation of the STORM model storm-time corrections for middle latitude

This paper presents results from the Storm-Time Ionospheric Correction Model (STORM) validation for selected Northern and Southern Hemisphere middle latitude locations. The created database incorporated 65 strong-to-severe geomagnetic storms, which occurred within the period 1995–2007. This validation included data from some ionospheric stations (e.g., Pruhonice, El Arenosillo) that were not considered in the development or previous validations of the model. Hourly values of the F2 layer critical frequency, foF2, measured for 5–7 days during the main and recovery phases of each selected storm were compared with the predicted IRI 2007 foF2 with the STORM model option activated. To perform a detailed comparison between observed values, medians and predicted foF2 values the correlation coefficient, the root-mean-square error (RMSE), and the percentage improvement were calculated. Results of the comparative analysis show that the STORM model captures more effectively the negative phases of the summer ionospheric storms, while electron density enhancement during winter storms and the changeover of the different storm phases is reproduced with less accuracy. The STORM model corrections are less efficient for lower-middle latitudes and severe geomagnetic storms.

Equatorial predictions from a new neural network based global foF2 model

A new neural network (NN) based global empirical model for the foF2 parameter, which represents the peak ionospheric electron density, has been developed using extended temporal and spatial geophysical relevant inputs. It has been proposed that this new model be considered as a suitable replacement for the International Union of Radio Science (URSI) and International Radio Consultative Committee (CCIR) model options currently used within the International Reference Ionosphere (IRI) model for the purpose of F2 peak electron density predictions. The most recent version of the model has incorporated data from 135 global ionospheric stations including a number of equatorial stations. This paper concentrates on the ability of this new model to predict foF2 for the equatorial sector, an area that has been identified as problematic within the current IRI peak prediction setup. The improvement in the predictions of the foF2 parameter by the new model as compared to the URSI and CCIR model options of the IRI is demonstrated and the requirement for additional foF2 data from the equatorial zone for the purpose of global modeling of foF2 is highlighted in this paper.

Satellite traces, range spread-F occurrence, and gravity wave propagation at the southern anomaly crest

Abstract. Range spread-F (RSF) and occurrence of "satellite" traces prior to RSF onset were studied at the southern peak of the ionospheric equatorial anomaly (EA). Ionograms recorded in September 2007 at the new ionospheric station of Tucumán, Argentina (26.9° S, 294.6° E, dip latitude 15.5° S), by the Advanced Ionospheric Sounder (AIS) developed at the Istituto Nazionale di Geofisica e Vulcanologia (INGV), were considered. Satellite traces (STs) are confirmed to be a necessary precursor to the appearance of an RSF trace on the ionograms. Moreover, an analysis of isoheight contours of electron density seems to suggest a relationship between RSF occurrence and gravity wave (GW) propagation.

F2 layer response to geomagnetic disturbances in Eastern Asia under the low solar activity

In this paper, the peculiarities of ionospheric response to geomagnetic disturbances observed at the decay and minimum of solar activity (SA) in the period 2004–2007 are investigated with respect to different geomagnetic conditions. Data from ionospheric stations and results of total electron content (TEC) measurements made at the network of GPS ground-based receivers located within the latitude–longitude sector (20–70°N, 90–160°Е) are used in this study. Three groups of anomalous ionospheric response to geomagnetic disturbances have been observed during low solar activity. At daytime, the large-scale traveling ionospheric disturbances (LSTIDs) could generally be related to the main phase of magnetic storm. Quasi-two-days wavelike disturbances (WLDs) have been also observed in the main phase independent of the geomagnetic storm intensity. Sharp electron density oscillations of short duration (OSD) occurred in the response to the onset of both main and recovery phases of the magnetic storm in the daytime at middle latitudes. A numerical model for ionosphere–plasmasphere coupling was used to interpret the occurrence of LS TIDs. Results showed that the LSTIDs might be associated with the unexpected lifting of F2 layer to the region with the lower recombination rate by reinforced meridional winds that produces the increase of the electron density in the F2 layer maximum.

Ionospheric precursors for crustal earthquakes in Italy

Abstract. Crustal earthquakes with magnitude 6.0>M≥5.5 observed in Italy for the period 1979–2009 including the last one at L'Aquila on 6 April 2009 were considered to check if the earlier obtained relationships for ionospheric precursors for strong Japanese earthquakes are valid for the Italian moderate earthquakes. The ionospheric precursors are based on the observed variations of the sporadic E-layer parameters (h'Es, fbEs) and foF2 at the ionospheric station Rome. Empirical dependencies for the seismo-ionospheric disturbances relating the earthquake magnitude and the epicenter distance are obtained and they have been shown to be similar to those obtained earlier for Japanese earthquakes. The dependences indicate the process of spreading the disturbance from the epicenter towards periphery during the earthquake preparation process. Large lead times for the precursor occurrence (up to 34 days for M=5.8–5.9) tells about a prolong preparation period. A possibility of using the obtained relationships for the earthquakes prediction is discussed.

Detection of large scale TIDs associated with the dayside cusp using SuperDARN data

Variations in the dayside ionosphere parameters as a result of a large-scale acoustic gravity wave (LS AGW) were studied for the 17 February 1998 substorm using the super dual auroral radar network (SuperDARN) measurements. This event was characterised by a sharp rise in the AE index with a maximum of ∼900 nT. The source of the disturbance responsible for the LS AGW appears to have been located within the plasma convection throat and in the dayside cusp region. The location of the source was obtained from studies of a number of datasets including high-latitude convection maps, data from 4 DMSP satellites and networks of ground-based magnetometers. The propagation of the LS AGWs caused quasi-periodic variations in the skip distance (with an amplitude up to 220–260 km) of the ground backscatter measured by up to 6 SuperDARN radars, including Goose Bay and Kapuskasing, resulting in two large-scale travelling ionospheric disturbances (LS TIDs). The LS TIDs had wave periods of 1.5 and 2 h, a velocity of ∼400 m/s for both, and wavelengths of 2200 and 2900 km, respectively. These quasi-periodic variations were also present in the peak electron density and height of the F2 layer measured by the Goose Bay ionosonde. The numerical simulation of the inverse problem show good agreement between Goose Bay radar and Goose Bay ionosonde measurements. But these LS TIDs would be difficult to deduce from the ground based ionospheric station data alone, because h mF2 variations were 10–40 km only and f OF2 variations between 10% and 20%. The results demonstrate how important SuperDARN radars can be, and that this is a more powerful technique than routine ground-based sounding for studies of weak quasi-periodic variations in the dayside subauroral ionosphere related to LS AGW.

Development of an Ionospheric Numerical Assimilation Nowcast and Forecast System Based on Gauss‐Markov Kalman Filter–an Observation System Simulation Experiment taking Example for China and Its Surrounding Area

AbstractIn this paper, we constructed an ionospheric data assimilation system based on Gauss‐Markov Kalman filter and gave some test results. We chose some ionosphere stations (including meridional project stations, China lithosphere deformation GPS network, part of IGS stations) in China and its surrounding areas as the observational system to do the simulation experiment. The International Reference Model (IRI) is chosen to be the background model, while NeQuick model output is taken to be the observations. Our assimilation results show that it can get good estimation of ionosphere electron density by ingesting the observed slant TEC data into the model by Kalman filter. It illustrates that our assimilation algorithm is feasible and the selected parameters are reasonable. We also obtained reasonable short‐time forecast results by Gauss‐Markov assumption.

Comparison of foF2 values in the daytime and after sunset

The ratio of the values of the F2, layer critical frequencies 2 h after the sunset and at 1400 LT has been considered. It is assumed that this ratio is a good indicator of the intensity of wind processes in the thermosphere. For 18 ionospheric stations, the dependences of this ratio on the Fwl, solar activity index have been constructed for the period before 1980. It has been found that the standard deviation (S), from this dependence has been increased after 1980 at all stations. An increase in S, is interpreted as an indication of the existence of long-term trends in th e wind regime after approximately 1980.

Electric fields and large-scale undulations in the evening sector of the diffuse auroral zone

The synchronous observations of strong electric fields and large-scale undulations observed on December 12, 2004, in the evening sector of the diffuse auroral zone 0900-1000 UT (~1700-1800 MLT) have been analyzed. The appearance of strong northward electric field at ~0900 UT was almost simultaneously registered at Tixie Bay ionospheric station (71.6° N, 128.9° E, L =, 5.6) and on the DMSP F15 satellite. At 0910-1000 UT, the all-sky TV camera at Tixie Bay and the DMSP satellites (F13, F14, and F15) registered eight undulations propagating westward at a velocity of 0.7—0.8 km/s. The undulation parameters registered during the TV observations agree with the satellite measurements. The distinctive feature of the analyzed event consists in that an intense electric field and undulations were localized within the diffuse zone in the region of increased precipitation of keV electrons. A comparison of the ground-based and satellite measurements made it possible to draw the conclusion on the necessary conditions for formation of diffuse undulations.

A qualitative study of the seismo-ionospheric precursors prior to the 6 April 2009 earthquake in L'Aquila, Italy

Abstract. In this paper we use the Cross Correlation analysis method in conjunction with the Empirical Mode Decomposition to analyze foF2 signals collected from Rome, Athens and San Vito ionospheric stations, in order to verify the existence of seismo-ionospheric precursors prior to M=6.3 L'Aquila earthquake in Italy. The adaptive nature of EMD allows for removing the geophysical noise from the foF2 signals, and then to calculate the correlation coefficient between them. According to the cross correlation coefficient theory, we expect the stations which located inside the earthquake preparation area, as evaluated using Dobrovolsky equation, to capture the ionospheric disturbances generated by the seismic event. On the other hand the stations outside of this area are expected to remain unaffected. The results of our study are in accordance with the theoretical model, evidencing ionospheric modification prior to L'Aquila earthquake in a certain area around the epicenter. However, it was found that the selection of stations at the limits of the theoretically estimated earthquake preparation area is not the best choice when the cross correlation method is applied, since the modification of the ionosphere over these stations may not be enough for the ionospheric precursors to appear. Our experimental results also show that when a seismic event constitutes the main shock after a series of pre-seismic activity, precursors may appear as early as 22 days prior to the event.

Variability of foF2 in the African equatorial ionosphere

This paper presents the impact of diurnal, seasonal and solar activity effects on the variability of ionospheric foF2 in the African equatorial latitude. Three African ionospheric stations; Dakar (14.8°N, 17.4°W, dip: 11.4°N), Ouagadougou (12.4°N, 1.5°W, dip: 2.8°N) and Djibouti (11.5°N, 42.8°E, dip: 7.2°N) were considered for the investigation. The overall aim is to provide African inputs that will be of assistance at improving existing forecasting models. The diurnal analysis revealed that the ionospheric critical frequency (foF2) is more susceptible to variability during the night-time than the day-time, with two peaks in the range; 18–38% during post-sunset hours and 35–55% during post-midnight hours. The seasonal and solar activity analyses showed a post-sunset September Equinox maximum and June Solstice maximum of foF2 variability in all the stations for all seasons. At all the stations, foF2 variability was high for low solar activity year. Overall, we concluded that equatorial foF2 variability increases with decreasing solar activity during night-time.

10.1007/s11478-008-1011-x

The results of observations of the solar eclipse ionospheric effects on March 29, 2006, are presented. The observations were conducted using the partial reflection method near Nizhni Novgorod and the vertical sounding method at the automatic ionospheric station near Murmansk. It has been obtained that the electron density at altitudes of 77 and 91 km decreases by a factor of more than 4; in this case the response of the ionosphere at an altitude of 91 km lags behind the eclipse maximum phase on the Earth by approximately 20 min. It has been established that the eclipse in the E and F1 regions of the polar ionosphere causes a change in the electron density by 15–20%. The delay time of this effect varies from 12 to 24 min depending on the altitude. It has been registered that the reflection virtual altitude at altitudes of the ionospheric F region increases in Murmansk and Nizhni Novgorod.

A dynamic system to forecast ionospheric storm disturbances based on solar wind conditions

For the reliable performance of technologically advanced radio communications systems under geomagnetically disturbed conditions, the forecast and modelling of the ionospheric response during storms is a high priority. The ionospheric storm forecasting models that are currently in operation have shown a high degree of reliability during quiet conditions, but they have proved inadequate during storm events. To improve their prediction accuracy, we have to take advantage of the deeper understanding in ionospheric storm dynamics that is currently available, indicating a correlation between the Interplanetary Magnetic Field (IMF) disturbances and the qualitative signature of ionospheric storm disturbances at middle latitude stations. In this paper we analyse observations of the foF2 critical frequency parameter from one mid-latitude European ionospheric station (Chilton) in conjunction with observations of IMF parameters (total magnitude, Bt and Bz-IMF component) from the ACE spacecraft mission for eight storm events. The determination of the time delay in the ionospheric response to the interplanetary medium disturbances leads to significant results concerning the forecast of the ionospheric storms onset and their development during the first 24 h. In this way the real-time ACE observations of the solar wind parameters may be used in the development of a real-time dynamic ionospheric storm model with adequate accuracy.

The new ionospheric station of Tucumán: first results

An Advanced Ionospheric Sounder, built at the Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy, was installed at Tucuman, Argentina, particularly interesting for its location, near the southern peak of the ionospher- ic equatorial anomaly. The aim of this installation is to collect a large number of continuous data useful both to study the dynamics of the equatorial ionospheric plasma and to develop reliable regional ionospheric prediction models. Moreover this ionosonde will contribute to the ionospheric database and real time knowledge of South- ern Hemisphere ionospheric conditions for space weather applications. The ionosonde is completely program- mable and two PCs support the data acquisition, control, storage and on-line processing. In this work the first results, in terms of ionograms and autoscaled characteristics, are presented and briefly discussed.

Inversion of the F2 layer critical frequency with the solar zenith angle

Normalization of the F2 layer critical frequency is put forward with a model related to the solar zenith angle at the time of observation and its local noon value. The method implies a physically based replacement of variable representing one of the key controlling parameters of the modern empirical iono� spheric models. The inversion of critical frequency is analyzed with the hourly values of foF2 of seven iono� spheric stations for 2000-2006 (half the solar activity cycle). The normalized critical frequency fnF2 reveals an improved correlation between the data of different stations and improved interseasonal correlation of the data of a particular station as compared with the initial foF2 values. A correlation radius in summer is two times greater for normalized critical frequency than for observations. The analytical model of expansion of the diurnal variation of critical frequency with four hourly values at 00, 06, 12 and 18 h of local time exhibits improved statistical features of normalized critical frequency as compared with observed critical frequency.

Comparative analysis of disturbances in the midlatitude ionospheric F 2 layer during strong and extreme magnetic storms in March 2001 according to the data from magnetically conjugate ground ionospheric stations

Ionospheric disturbances at heights of the F2 layer maximum during the strong magnetic storm (the minimum value of the Dst index was −149 nT) and the magnetic superstorm (the minimum value of the Dst index was −387 nT) have been compared based on the data from two pairs of magnetically conjugate midlatitude ground stations for ionospheric vertical sounding. The storms began on March 19, 2001, and March 31, 2001, respectively. It has been obtained that almost only negative ionospheric disturbances were observed in the Northern and Southern hemispheres in both cases. The maximum effect in changes in the layer critical frequency (foF2) in both hemispheres has a time delay relative to the moment of the maximum disturbance in the Dst index on the order of 3–4 h for the strong storm and about 1 h for the superstorm. The disturbed variations in the foF2 critical frequency in different hemispheres correlate well with each other in the plane of one magnetic meridian, but the correlation substantially weakens at different magnetic longitudes. An assumption is made that the revealed features of the behavior of the disturbed midlatitude ionospheric F2 layer are caused by the complex character of the thermospheric response to the energy release in the auroral zone during the considered magnetic storms.