Diagnostics of equatorial and low latitude ionosphere by TEC mapping over Brazil

The total electron content (TEC) in the equatorial and low-latitude ionosphere over Brazil was monitored in two dimensions by using 2011 data from the ground-based global navigation satellite system (GNSS) receiver network operated by the Brazilian Institute for Geography and Statistics. It was possible to monitor the spatial and temporal variations in TEC over Brazil continuously during both day and night with a temporal interval of 10 min and a spatial resolution of about 400 km. The daytime equatorial ionization anomaly (EIA) and post-sunset plasma enhancement (PS-EIA) were monitored over an area corresponding to a longitudinal extension of 4000 km in South America. Considerable day-to-day variation was observed in EIA and PS-EIA. A large latitudinal and longitudinal gradient of TEC indicated a significant ionospheric range error in application of the GNSS positioning system. Large-scale plasma bubbles after sunset were also mapped over a wide range. Depletions with longitudinally separated by more than 800 km were observed. They were extended by more than 2000 km along the magnetic field lines and drifted eastward. It is expected that 2-dimensional TEC mapping can serve as a useful tool for diagnosing ionospheric weather, such as temporal and spatial variation in the equatorial plasma trough and crest, and particularly for monitoring the dynamics of plasma bubbles.     

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.     

Ionospheric path delay models for spaceborne GPS receivers flying in formation with large baselines

GPS relative navigation filters could benefit notably from an accurate modeling of the ionospheric delays, especially over large baselines (>100 km) where double difference delays can be higher than several carrier wavelengths. This paper analyzes the capability of ionospheric path delay models proposed for spaceborne GPS receivers in predicting both zero-difference and double difference ionospheric delays. We specifically refer to relatively simple ionospheric models, which are suitable for real-time filtering schemes. Specifically, two ionospheric delay models are evaluated, one assuming an isotropic electron density and the other considering the effect on the electron density of the Sun aspect angle. The prediction capability of these models is investigated by comparing predicted ionospheric delays with measured ones on real flight data from the Gravity Recovery and Climate Experiment mission, in which two satellites fly separated of more than 200 km. Results demonstrate that both models exhibit a correlation in the excess of 80% between predicted and measured double-difference ionospheric delays. Despite its higher simplicity, the isotropic model performs better than the model including the Sun effect, being able to predict double differenced delays with accuracy smaller than the carrier wavelength in most cases. The model is thus fit for supporting integer ambiguity fixing in real-time filters for relative navigation over large baselines. Concerning zero-difference ionospheric delays, results demonstrate that delays predicted by the isotropic model are highly correlated (around 90%) with those estimated using GPS measurements. However, the difference between predicted and measured delays has a root mean square error in the excess of 30 cm. Thus, the zero-difference ionospheric delays model is not likely to be an alternative to methods exploiting carrier-phase observables for cancelling out the ionosphere contribution in single-frequency absolute navigation filters.     

Global characteristics of ionospheric electric fields and disturbances during the first hours of magnetic storms

The ionospheric plasma density can be significantly disturbed during magnetic storms. In the conventional scenario of ionospheric storms, the negative storm phases with plasma density decreases are caused by neutral composition changes, and the positive storm phases with plasma density increases are often related to atmospheric gravity waves. However, recent studies show that the global redistribution of the ionospheric plasma is dominated primarily by electric fields during the first hours of magnetic storms. In this paper, we present the measurements of ionospheric disturbances by the DMSP satellites and GPS network during the magnetic storm on 6 April 2000. The DMSP measurements include the F region ion velocity and density at the altitude of ∼840 km, and the GPS receiver network provides total electron content (TEC) measurements. The storm-time ionospheric disturbances show the following characteristics. The plasma density is deeply depleted in a latitudinal range of ∼20° over the equatorial region in the evening sector, and the depletions represent plasma bubbles. The ionospheric plasma density at middle latitudes (20°–40° magnetic latitudes) is significantly increased. The dayside TEC is increased simultaneously over a large latitudinal range. An enhanced TEC band forms in the afternoon sector, goes through the cusp region, and enters the polar cap. All the observed ionospheric disturbances occur within 1–5 h from the storm sudden commencement. The observations suggest that penetration electric fields play a major role in the rapid generation of equatorial plasma bubbles and the simultaneous increases of the dayside TEC within the first 2 h during the storm main phase. The ionospheric disturbances at later times may be caused by the combination of penetration electric fields and neutral wind dynamo process.     

Variation in the ionospheric propagating factor M(3000)F2 at Ouagadougou,Burkina Faso

We use hourly monthly median values of propagation factor M(3000)F2 data observed at Ouagadougou Ionospheric Observatory (geographic12.4°N, 1.5°W; 5.9^o dip), Burkina Faso (West Africa) during the years Januar1987–December1988 (average F10.7 < 130 × 10⁻²² W/m²/Hz, representative of low solar flux conditions) and for January 1989–December1990 (average F10.7 ? 130 × 10⁻²² W/m²/Hz, representative of high solar epoch) for magnetically quiet conditions to describe local time, seasonal and solar cycle variations of equatorial ionospheric propagation factor M(3000)F2 in the African region. We show that that seasonal trend between solar maximum and solar minimum curves display simple patterns for all seasons and exhibits reasonable disparity with root mean square error (RMSE) of about 0.31, 0.29 and 0.26 for December solstice, June solstice and equinox, respectively. Variability Σ defined by the percentage ratio of the absolute standard deviation to the mean indicates significant dissimilarity for the two solar flux levels. Solar maximum day (10–14 LT) and night (22–02 LT) values show considerable variations than the solar minimum day and night values. We compare our observations with those of the IRI 2007 to validate the prediction capacity of the empirical model. We find that the IRI model tends to underestimate and overestimate the observed values of M(3000)F2, in particular, during June solstice season. There are large discrepancies, mainly during high solar flux equinox and December solstice between dawn and local midnight. On the other hand, IRI provides a slightly better predictions for M(3000)F2 between 0900 and 1500 LT during equinox low and high solar activity and equinox high sunspot number. Our data are of great importance in the area of short-wave telecommunication and ionospheric modeling.     

Comparative analysis of TEC disturbances over tropical cyclone zones in the North–West Pacific Ocean

In this paper, we study ionospheric total electron content (TEC) disturbances associated with tropical cyclones (TCs). The study relies on the statistical analysis of six cyclones of different intensity which occurred in the North–West Pacific Ocean in September–November 2005. We have used TEC data from the international network of two-frequency ground-based GPS receivers and NCEP/NCAR meteorological archive. TEC variations of different period ranges (02–20 and 20–60 min) are shown to be more intense during TC peaks under quiet geomagnetic conditions. The highest TEC variation amplitudes are registered when the wind speed in the cyclone and the TC area are maximum. The intensification of TEC disturbances is more pronounced when several cyclones occur simultaneously. We have revealed that the ionospheric response to TC can be observed only after the cyclone has reached typhoon intensity. The ionospheric response is more pronounced at low satellite elevation angles.     

The GFZ real-time GNSS precise positioning service system and its adaption for COMPASS

Motivated by the IGS real-time Pilot Project, GFZ has been developing its own real-time precise positioning service for various applications. An operational system at GFZ is now broadcasting real-time orbits, clocks, global ionospheric model, uncalibrated phase delays and regional atmospheric corrections for standard PPP, PPP with ambiguity fixing, single-frequency PPP and regional augmented PPP. To avoid developing various algorithms for different applications, we proposed a uniform algorithm and implemented it into our real-time software. In the new processing scheme, we employed un-differenced raw observations with atmospheric delays as parameters, which are properly constrained by real-time derived global ionospheric model or regional atmospheric corrections and by the empirical characteristics of the atmospheric delay variation in time and space. The positioning performance in terms of convergence time and ambiguity fixing depends mainly on the quality of the received atmospheric information and the spatial and temporal constraints. The un-differenced raw observation model can not only integrate PPP and NRTK into a seamless positioning service, but also syncretize these two techniques into a unique model and algorithm. Furthermore, it is suitable for both dual-frequency and sing-frequency receivers. Based on the real-time data streams from IGS, EUREF and SAPOS reference networks, we can provide services of global precise point positioning (PPP) with 5–10 cm accuracy, PPP with ambiguity-fixing of 2–5 cm accuracy, PPP using single-frequency receiver with accuracy of better than 50 cm and PPP with regional augmentation for instantaneous ambiguity resolution of 1–3 cm accuracy. We adapted the system for current COMPASS to provide PPP service. COMPASS observations from a regional network of nine stations are used for precise orbit determination and clock estimation in simulated real-time mode, the orbit and clock products are applied for real-time precise point positioning. The simulated real-time PPP service confirms that real-time positioning services of accuracy at dm-level and even cm-level is achievable with COMPASS only.     

Ionospheric precursors observed at low latitudes around the time of koyna earthquake

On December 11, 1967 at 05:21 LT, an immense earthquake of magnitude 6.7 struck Koyna, the Indian province of Maharashtra. Its epicenter was located at geographic latitude 17.37°N and longitude 73.75°E with depth of about 3 km. Ground based measurements show variation in the critical frequency of ionospheric F2 layer (foF2) before and after the shock. In the present study the behavior of F2-region of ionosphere has been examined over the equatorial and low latitudinal region ionosphere during the month of December 1967 around the time of Koyna earthquake. For this purpose, the ionospheric data collected with the help of ground-based ionosondes installed at Hyderabad (located close to the earthquake epicenter) Ahmedabad, Trichirapulli, Kodaikanal and Trivendrum have been utilized. The upper and lower bound of Interquartile range (IRQ) are constructed to monitor the variations in foF2 other than day-to-day and diurnal pattern for finding the seismo-ionospheric precursors. Some anomalous electron density variations are observed between post midnight hours to local pre-noon hours at each station. These anomalies are strongly time dependent and appeared a couple of days before the main shock. The period considered in this study comes under the quiet geomagnetic conditions. Hence, the observed anomalies (which are more than the usual day-to-day variability) over all stations are likely to be associated with this imminent earthquake. The possible mechanism to explain these anomalies is the effect of seismogenic electric field generated just above the surface of earth within the earthquake preparation zone well before the earthquake due to emission of radioactive particles and then propagated upward, which perturbs the F-region ionosphere.     

Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

Difficulties in the study of cosmic radio noise absorption at 30 MHz using riometer at low latitude station,Kolhapur (Lat-16.8°N,Long-74.25°E)

A dual dipole antenna has been installed at low latitude station Kolhapur (Geographic 16.8°N, 74.25°E), Maharashtra, India for the study of cosmic radio noise absorption using Solid State Riometer (which operates at 30 MHz) during pre phase of 24th solar maxima. The aim for this type of study over Kolhapur was to know the response of lower (D region) ionosphere over low latitude by cosmic radio noise absorption using riometer technique during quite period as well as sudden ionospheric disturbances (SID). The observations are being taken for 3 years. Two different sites (∼40 km away from each other) were used for the installation of riometer equipment assuming minimum local noise. It is found that solar noise to cosmic radio noise hence resulting in signal saturation. The night time signal is relatively free of interference but sometimes local noise is responsible for spike-like signatures. Hence it is concluded that Kolhapur (a low latitude station) is not suitable for the study of cosmic radio noise absorption on 30 MHz with riometer and dual dipole antenna. Proper choice for operating frequency of riometer and antenna gain is suggested for low latitude use of this technique for ionospheric deviative and nondeviative absorption studies.     

Study of mid-latitude ionospheric convection during quiet and disturbed periods using the SuperDARN Hokkaido radar

Westward ionospheric convective flows around midnight are frequently observed at mid-latitudes. They can be generated by so-called disturbance dynamo mechanisms working mainly in the mid-latitudes. To understand the influence of disturbance dynamo effects in the mid-latitudes, we studied the latitudinal distribution of westward flows in association with several kinds of geomagnetic disturbances using the SuperDARN Hokkaido radar. This radar creates high temporal resolution (1 s to 2 min), two-dimensional observations measuring the line-of-sight velocities of ionospheric plasma irregularities, which can be regarded as line-of-sight velocities of ionospheric convection in the mid-latitude region from 40° to 50°. This region could not be monitored using preexisting SuperDARN radars. In this study, we used ionospheric echo data obtained by the SuperDARN Hokkaido radar over 5 years (December 2006 to November 2011). We identified westward flows around midnight at about 40° to 55° geomagnetic latitude. Additionally, the data showed that the westward flow around midnight intensified under high geomagnetic activity (high Kp). This suggests that the disturbance dynamo could affect the mid-latitude ionospheric convection. We performed Superposed Epoch Analysis (SEA) to study the influences from the geomagnetic disturbances on mid-latitude ionospheric convection. We found no obvious influence during major storms (minimum Dst below −60 nT). SEA was also used to study the temporal and latitudinal dependence on the influences from substorms. From analysis of 36 events of AL-defined substorms, we saw that the influence of substorms lasted from 5 to 20 h after the onset between 44° and 53° geomagnetic latitude. The westward flow at mid-latitude grew to a maximum at 12 h after the geomagnetic substorm onset. This is consistent with the results of past numerical simulation studies of the disturbance dynamo effects.     

Ionospheric scintillations at Guilin detected by GPS ground-based and radio occultation observations

The occurrence of ionospheric scintillations with S4 ? 0.2 was studied using GPS measurements at Guilin, China (25.29°N, 110.33°E; geomagnetic: 15.04°N, 181.98°E), a station located near the northern crest of the equatorial anomaly. The results are presented for data collected from January 2009 to March 2010. The results show that nighttime amplitude scintillations only took place in February and March of the considered years, while daytime amplitude scintillations occurred in August and December of 2009. Nighttime amplitude scintillations, observed in the south of Guilin, always occurred with phase scintillations, TEC (Total Electron Content) depletions, and ROT (Rate Of change of TEC) fluctuations. However, TEC depletions and ROT fluctuations were weak during daytime amplitude scintillations, and daytime amplitude scintillations always took place simultaneously for most of the GPS satellites which appeared over Guilin in different azimuth directions. Ground-based GPS scintillation/TEC observations recorded at Guilin and signal-to-noise-ratio (SNR) measurements obtained from GPS-COSMIC radio occultation indicate that nighttime and daytime scintillations are very likely caused by ionospheric F region irregularities and sporadic E, respectively. Moreover, strong daytime amplitude scintillations may be associated with the plasma density enhancements in ionospheric E region caused by the Perseid and Geminid meteor shower activities.     

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.     

The COMPLEIK subroutine of the IONORT-ISP system for calculating the non-deviative absorption: A comparison with the ICEPAC formula

The present paper proposes to discuss the ionospheric absorption, assuming a quasi-flat layered ionospheric medium, with small horizontal gradients. A recent complex eikonal model (Settimi et al., 2013b) is applied, useful to calculate the absorption due to the ionospheric D-layer, which can be approximately characterized by a linearized analytical profile of complex refractive index, covering a short range of heights between h₁ = 50 km and h₂ = 90 km. Moreover, Settimi et al. (2013c) have already compared the complex eikonal model for the D-layer with the analytical Chapman’s profile of ionospheric electron density; the corresponding absorption coefficient is more accurate than Rawer’s theory (1976) in the range of middle critical frequencies. Finally, in this paper, the simple complex eikonal equations, in quasi-longitudinal (QL) approximation, for calculating the non-deviative absorption coefficient due to the propagation across the D-layer are encoded into a so called COMPLEIK (COMPLex EIKonal) subroutine of the IONORT (IONOspheric Ray-Tracing) program ( Azzarone et al., 2012). The IONORT program, which simulates the three-dimensional (3-D) ray-tracing for high frequencies (HF) waves in the ionosphere, runs on the assimilative ISP (IRI-SIRMUP-P) discrete model over the Mediterranean area ( Pezzopane et al., 2011). As main outcome of the paper, the simple COMPLEIK algorithm is compared to the more elaborate semi-empirical ICEPAC formula (Stewart, undated), which refers to various phenomenological parameters such as the critical frequency of E-layer. COMPLEIK is reliable just like the ICEPAC, with the advantage of being implemented more directly. Indeed, the complex eikonal model depends just on some parameters of the electron density profile, which are numerically calculable, such as the maximum height.     

Fading in the HF ionospheric channel and the role of irregularities

It is well known that the ionosphere affects radio wave propagation especially in the high frequency (HF) range. HF radio waves reflected by the ionosphere can reach considerable distances, often with changes in amplitude, phase, and frequency. The ionosphere is a dispersive in frequency and time, bi-refractive, absorbing medium, in which multipath propagation due to traveling irregularities is very frequent. The traveling irregularities undulate the reflecting ionospheric layer, introducing variations in signal amplitude (fading). In this multipath time variant channel fading is mainly considered, even though it is not the sole effect. Echo signals from a single reflection, as in ionospheric vertical sounding (VIS) techniques, are affected by a certain degree of variability even in quiet ionospheric conditions. In this work the behavior of the ionospheric channel is studied and characterized by observing the power variation of received echoes using the VIS technique. Multipath fading was analyzed quantifying the power variation of the signal echo due to irregularities on a temporal scale from 0.5 to 256 s. An experimental set-up derived from an ionosonde was implemented and the analysis was performed employing a special numerical algorithm operating off-line on the acquired time sequence of the signal. The gain-loss of the irregularity shapes are determined in some special cases.     

Improved predictions of equatorial TEC for real time applications

It is important to use models developed specifically for the equatorial ionospheric estimation for real-time applications, particularly in Satellite Navigation. This work demonstrates a methodology for improved predictions of VTEC in real time using the model developed for the equatorial ionosphere by the authors. This work has been done using TEC data of the low solar activity period of 2005 obtained using dual frequency GPS receivers installed under the GAGAN project of ISRO. For the purpose, the model is first used in conjunction with Kriging technique. Improvement in accuracy is observed when compared with the estimations from the model alone using the measurements as true reference. Further improvement is obtained by Bayesian combination of these estimates with independent Neural Network based predictions. Statistical performance of improvement is provided. An improvement of ∼1 m in confidence level of estimation of VTEC is obtained.     

Modeling the ionospheric gradient by using VSS base functions and GPS data over Iran

The spatial and temporal variations of ionosphere play an important role in positioning and navigation by the space geodetic techniques. Therefore, the ionospheric gradient should be calculated, analyzed, and applied in both space and time. Spatial gradients of the ionosphere have remarkable delay on the propagation of electromagnetic waves. This study intends to propose a new method for simultaneous modeling of the spatial gradients of ionosphere and VTECs in the local scale for Iran. Vector Spherical Slepian (VSS) base functions are used for the development of this method.Five VSS models with the maximal degrees of L = 30, 35, 40, 45 and 50 are taken into account. For implementing the VSS method, 24 permanent GPS stations from the Iranian Permanent GPS Network (IPGN) have been used. The unknown coefficients are estimated with the observations of these stations with least squares technique. Four other stations are used for evaluating the accuracy of the models. Repeatability of baselines is the measure that is used for this purpose. Based on the results obtained, L = 40 is the optimum degree for the VSS model with this input data over Iran.The baselines’ repeatability showed that ionospheric gradients have more influence on the north–south component. Moreover, the spatial gradient is negligible in the east–west and up-down component when a short baseline is processed. In other words, this kind of ionospheric modeling has significant application for baseline, which is longer than 1000 km. In the study, proposed method has improved the long baselines' solution by more than 12%, 18% and 10% in east–west, north–south and up-down components, respectively.     

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.     

Demonstration of the use of the Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS) measurements to validate GPS ionospheric imaging

There is a lack of independent ionospheric data that can be used to validate GPS imaging results at mid latitudes over severe storm times. Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS), a global network of dual-frequency ground to satellite observations, provides this missing data and here is employed as verification to show the accuracy of the ionospheric GPS images in terms of the total electron content (TEC). In this paper, the large-scale ionospheric structures that appeared during the strong geomagnetic storm of 20 November 2003 are reconstructed with a GPS tomographic algorithm, known as MIDAS, and validated with DORIS TEC measurements. The main trough shown in an extreme equatorward position in the ionospheric imaging over mainland Europe is confirmed by DORIS satellite measurements. Throughout the disturbed day, the variations of relative slant TECs between DORIS data and MIDAS results agree quite well, with the average of the mean differences about 2 TECu. We conclude that as a valuable supplement to GPS data, DORIS ionospheric measurements can be used to analyse TEC variations with a relatively high resolution, ∼10 s in time and tens of kilometres in space. This will be very helpful for identification of some highly dynamic structures in the ionosphere found at mid-latitudes, such as the main trough, TID (Travelling Ionospheric Disturbances) and SED (Storm Enhanced Density), and could be used as a valuable auxiliary data source in ionospheric imaging.     

Variability of the ionospheric electron density at fixed heights and validation of IRI-2007 profile’s prediction at Ilorin

A study on the variability of the equatorial ionospheric electron density was carried out at fixed heights below the F2 peak using one month data for each of high and low solar activity periods. The data used for this study were obtained from ionograms recorded at Ilorin, Nigeria, and the study covers height range from 100 km to the peak of the F2 layer for the daytime hours and height range from 200 km to the peak of the F2 layer for the nighttime hours. The results showed that the deviation of the electron density variation from simple Chapman variation begins from an altitude of about 200 km for the two months investigated. Daytime minimum variability of between 2.7% and 9.0% was observed at the height range of about 160 and 200 km during low solar activity (January 2006) and between 3.7% and 7.8% at the height range of 210 and 260 km during high solar activity (January 2002). The nighttime maximum variability was observed at the height range of 210 and 240 km at low solar activity and at the height range of 200 and 240 km at high solar activity. A validation of IRI-2007 model electron density profile’s prediction was also carried out. The results showed that B0 option gives a better prediction around the noontime.