Epochwise prediction of GPS single differenced ionospheric delays of formation flying spacecraft

Both single and dual frequency GPS relative navigation filters may benefit from proper predictions of single differenced ionospheric delays. In this article, the single differenced ionospheric delays of GPS observations are predicted for the GRACE formation during the switch manoeuvre.Two prediction methods are considered. The first is based on a Taylor expansion to first order of a mapping function that maps slant total electron content measurements to vertical total electron content estimates. The second method fits a shape profile through undifferenced ionospheric data available. It then raytraces through this profile to estimate the difference in total electron content along the path of the GPS signals.Continuously changing ionospheric conditions hamper the assessment of the quality of the predictions. Comparison of both methods shows that the raytracing method performs better. The difference of predictions and measurements generally shows a smaller RMS than the measurements alone. However, both methods suffer from a number of systematically unpredicted observations, which arise from small unaccounted differential variations in electron densities along the path of the GPS signals. These prediction methods perform better when spacecraft separation is small. Baselines considered here range from tens of kilometres down to several hundred metres. When smallest spacecraft separation occurs (0.4 km), the single differenced ionospheric delays exhibit RMS values of 0.0089 m. The first method shows a difference between measurements and predictions with an RMS of 0.0081 m. For the second method the difference RMS is found to be 0.0067 m.     

Medium-scale traveling ionospheric disturbances (MSTID) modeling using a dense German GPS network

The MSTIDs are wave-like perturbations of the ionospheric plasma, which cause the most common ionospheric disturbances in mid-latitude regions. Generally the MSTIDs have velocities of several hundred meters per second and wavelengths of several hundred kilometers. The wave-like effect of the MSTID is one of the main obstacles for accurate interpolation of ionospheric corrections in a medium-scale reference GPS network. In this paper we show a new method of detecting and modeling MSTIDs using dense German GPS network. The between-epoch single difference ionospheric delays from a medium scale dense GPS network are used to estimate the parameter of the MSTID e.g. amplitude, wavelength and velocity. The efficiency of the approach is tested with data from about 320 GPS stations in and near Germany. A MSTID wave moving from east to west across Germany was observed at September 27 in 2009. Its wavelength is about 302 km, with a period of ∼7 min and velocity of about 700 m/s.     

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.     

Ionospheric differential error determination using ray tracing for a short baseline

Since the United States government discontinued Selective Availability (SA) on 1 May 2000, ionospheric effects have been responsible for the largest errors in GPS systems. The standard Differential GPS (DGPS) method is incapable of completely eliminating the ionospheric error. This paper describes a new approach to determine the differential ionospheric error between geographically distributed receiver stations. The ray paths of GPS signals were simulated using a modified Jones 3D ray tracing programme that includes the effect of the geomagnetic field. A Nelder–Mead optimisation algorithm was embedded in the program to precisely determine the satellite-to-station path. A realistic ionospheric model is essential for accurate ray tracing results and for estimates of differential error that are accurate on sub-centimetre scales. Here, the ionospheric model used in the ray tracing programme was developed by fitting realistic ionosphere profiles with a number of exponential functions. Results were compared to the theoretical approach. Results show that the differential delay is about 1–5 cm at low elevation angles for a short baseline of 10 km, as reported in other literature. This delay is often neglected in DGPS application. The differential delay also shows a pattern similar to that predicted by the Klobuchar model. The method proposed here can be used to improve future GPS applications.     

TEC derived from some GPS stations in Nigeria and comparison with the IRI and NeQuick models

Total electron content (TEC) measured simultaneously using Global Positioning System (GPS) ionospheric monitors installed at some locations in Nigeria during the year 2011 (Rz = 55.7) was used to study the diurnal, seasonal, and annual TEC variations. The TEC exhibits daytime maximum, seasonal variation and semiannual variations. Measured TEC were compared with those predicted by the improved versions of the International Reference Ionosphere (IRI) and NeQuick models. The models followed the diurnal and seasonal variation patterns of the observed values of TEC. However, IRI model produced better estimates of TEC than NeQuick at all locations.     

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.     

The use of ionospheric tomography and elevation masks to reduce the overall error in single-frequency GPS timing applications

Signals from Global Positioning System (GPS) satellites at the horizon or at low elevations are often excluded from a GPS solution because they experience considerable ionospheric delays and multipath effects. Their exclusion can degrade the overall satellite geometry for the calculations, resulting in greater errors; an effect known as the Dilution of Precision (DOP). In contrast, signals from high elevation satellites experience less ionospheric delays and multipath effects. The aim is to find a balance in the choice of elevation mask, to reduce the propagation delays and multipath whilst maintaining good satellite geometry, and to use tomography to correct for the ionosphere and thus improve single-frequency GPS timing accuracy. GPS data, collected from a global network of dual-frequency GPS receivers, have been used to produce four GPS timing solutions, each with a different ionospheric compensation technique. One solution uses a 4D tomographic algorithm, Multi-Instrument Data Analysis System (MIDAS), to compensate for the ionospheric delay. Maps of ionospheric electron density are produced and used to correct the single-frequency pseudorange observations. This method is compared to a dual-frequency solution and two other single-frequency solutions: one does not include any ionospheric compensation and the other uses the broadcast Klobuchar model. Data from the solar maximum year 2002 and October 2003 have been investigated to display results when the ionospheric delays are large and variable. The study focuses on Europe and results are produced for the chosen test site, VILL (Villafranca, Spain). The effects of excluding all of the GPS satellites below various elevation masks, ranging from 5° to 40°, on timing solutions for fixed (static) and mobile (moving) situations are presented. The greatest timing accuracies when using the fixed GPS receiver technique are obtained by using a 40° mask, rather than a 5° mask. The mobile GPS timing solutions are most accurate when satellites at lower elevations continue to be included: using a mask between 10° and 20°. MIDAS offers the most accurate and least variable single-frequency timing solution and accuracies to within 10 ns are achieved for fixed GPS receiver situations. Future improvements are anticipated by combining both GPS and Galileo data towards computing a timing solution.     

Effects of physical correlations on long-distance GPS positioning and zenith tropospheric delay estimates

The global positioning system (GPS) has become an essential tool for the high precision navigation and positioning. The quality of GPS positioning results mainly depends on the model’s formulations regarding GPS observations, including both a functional model, which describes the mathematical relationships between the GPS measurements and unknown parameters, and a stochastic model, which reflects the physical properties of the measurements. Over the past two decades, the functional models for GPS measurements have been investigated in considerable detail. However, the stochastic models of GPS observation data are simplified, assuming that all the GPS measurements have the same variance and are statistically independent. Such assumptions are unrealistic. Although a few studies of GPS stochastic models were performed, they are restricted to short baselines and short time session lengths. In this paper, the stochastic modeling for GPS long-baseline and zenith tropospheric delay (ZTD) estimates with a 24-h session is investigated using the residual-based and standard stochastic models. Results show that using the different stochastic modelling methods, the total differences can reach as much as 3–6 mm in the baseline component, especially in the height component, and 10 mm in the ZTD estimation. Any misspecification in the stochastic models will result in unreliable GPS baseline and ZTD estimations. Using the residual-based stochastic model, not only the precision of GPS baseline and ZTD estimation is obviously improved, but also the baseline and ZTD estimations are closer to the reference value.     

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.     

Evaluation of COMPASS ionospheric model in GNSS positioning

As important products of GNSS navigation message, ionospheric delay model parameters are broadcasted for single-frequency users to improve their positioning accuracy. GPS provides daily Klobuchar ionospheric model parameters based on geomagnetic reference frame, while the regional satellite navigation system of China’s COMPASS broadcasts an eight-parameter ionospheric model, COMPASS Ionospheric Model(CIM), which was generated by processing data from continuous monitoring stations, with updating the parameters every 2 h. To evaluate its performance, CIM predictions are compared to ionospheric delay measurements, along with GPS positioning accuracy comparisons. Real observed data analysis indicates that CIM provides higher correction precision in middle-latitude regions, but relatively lower correction precision for low-latitude regions where the ionosphere has much higher variability. CIM errors for some users show a common bias for in-coming COMPASS signals from different satellites, and hence ionospheric model errors are somehow translated into the receivers’ clock error estimation. In addition, the CIM from the China regional monitoring network are further evaluated for global ionospheric corrections. Results show that in the Northern Hemisphere areas including Asia, Europe and North America, the three-dimensional positioning accuracy using the CIM for ionospheric delay corrections is improved by 7.8%–35.3% when compared to GPS single-frequency positioning ionospheric delay corrections using the Klobuchar model. However, the positioning accuracy in the Southern Hemisphere is degraded due apparently to the lack of monitoring stations there.     

Assessment of precision in ionospheric electron density profiles retrieved by GPS radio occultations

The Constellation Observing System for Meteorology Ionosphere and Climate (COSMIC) is a six satellite radio occultation mission that was launched in April 2006. The close proximity of these satellites during some months after launch provides a unique opportunity to evaluate the precision of Global Positioning System (GPS) radio occultation (RO) retrievals of ionospheric electron density from nearly collocated and simultaneous observations. RO data from 30 consecutive days during July and August 2006 are divided into ten groups in terms of daytime or nighttime and latitude. In all cases, the best precision values (about 1%) are found at the F peak height and they slightly degrade upwards. For all daytime groups, it is seen that electron density profiles above about 120 km height exhibit a substantial improvement in precision. Nighttime groups are rather diverse: in particular, the precision becomes better than 10% above different levels between 120 and 200 km height. Our overall results show that up to 100–200 km (depending on each group), the uncertainty associated with the precision is in the order of the measured electron density values. Even worse, the retrieved values tend sometimes to be negative. Although we cannot rely directly on electron density values at these altitudes, the shape of the profiles could be indicative of some ionospheric features (e.g. waves and sporadic E layers). Above 200 km, the profiles of precision are qualitatively quite independent from daytime or latitude. From all the nearly collocated pairs studied, only 49 exhibited a difference between line of sight angles of both RO at the F peak height larger than 10°. After analyzing them we find no clear indications of a significant representativeness error in electron density profiles due to the spherical assumption above 120 km height. Differences in precision between setting and rising GPS RO may be attributed to the modification of the processing algorithms applied to rising cases during the initial period of the COSMIC mission.     

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.     

Plasmaspheric electron content derived from GPS TEC and FORMOSAT-3/COSMIC measurements: Solar minimum condition

The plasmaspheric electron content (PEC) was estimated by comparison of GPS TEC observations and FORMOSAT-3/COSMIC radio occultation measurements at the extended solar minimum of cycle 23/24. Results are retrieved for different seasons (equinoxes and solstices) of the year 2009. COSMIC-derived electron density profiles were integrated up to the height of 700 km in order to retrieve estimates of ionospheric electron content (IEC). Global maps of monthly median values of COSMIC IEC were constructed by use of spherical harmonics expansion. The comparison between two independent measurements was performed by analysis of the global distribution of electron content estimates, as well as by selection specific points corresponded to mid-latitudes of Northern America, Europe, Asia and the Southern Hemisphere. The analysis found that both kinds of observations show rather similar diurnal behavior during all seasons, certainly with GPS TEC estimates larger than corresponded COSMIC IEC values. It was shown that during daytime both GPS TEC and COSMIC IEC values were generally lower at winter than in summer solstice practically over all specific points. The estimates of PEC (h > 700 km) were obtained as a difference between GPS TEC and COSMIC IEC values. Results of comparative study revealed that for mid-latitudinal points PEC estimates varied weakly with the time of a day and reached the value of several TECU for the condition of solar minimum. Percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 50–60%) during night-time and lesser values (25–45%) during day-time.     

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).     

Onboard orbit determination using GPS observations based on the unscented Kalman filter

Spaceborne GPS receivers are used for real-time navigation by most low Earth orbit (LEO) satellites. In general, the position and velocity accuracy of GPS navigation solutions without a dynamic filter are 25 m (1σ) and 0.5 m/s (1σ), respectively. However, GPS navigation solutions, which consist of position, velocity, and GPS receiver clock bias, have many abnormal excursions from the normal error range for space operation. These excursions lessen the accuracy of attitude control and onboard time synchronization. In this research, a new onboard orbit determination algorithm designed with the unscented Kalman filter (UKF) was developed to improve the performance. Because the UKF is able to obtain the posterior mean and covariance accurately by using the second-order Taylor series expansion through the sampled sigma points that are propagated by using the true nonlinear system, its performance can be better than that of the extended Kalman filter (EKF), which uses the linearized state transition matrix to predict the covariance. The dynamic models for orbit propagation applied perturbations due to the 40 × 40 geo-potential, the gravity of the Sun and Moon, solar radiation pressure, and atmospheric drag. The 7(8)th-order Runge–Kutta numerical integration was applied for orbit propagation. Two types of observations, navigation solutions and C/A code pseudorange, can be used at the user’s discretion. The performances of the onboard orbit determination were verified using real GPS data of the CHAMP and KOMPSAT-2 satellites. The results of the orbit determination were compared with the precision orbit ephemeris (POE) of the CHAMP and KOMPSAT-2 satellites.     

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.     

A new method for improving the performance of an ionospheric model developed by multi-instrument measurements based on artificial neural network

There are remarkable ionospheric discrepancies between space-borne (COSMIC) measurements and ground-based (ionosonde) observations, the discrepancies could decrease the accuracies of the ionospheric model developed by multi-source data seriously. To reduce the discrepancies between two observational systems, the peak frequency (foF2) and peak height (hmF2) derived from the COSMIC and ionosonde data are used to develop the ionospheric models by an artificial neural network (ANN) method, respectively. The averaged root-mean-square errors (RMSEs) of COSPF (COSMIC peak frequency model), COSPH (COSMIC peak height model), IONOPF (Ionosonde peak frequency model) and IONOPH (Ionosonde peak height model) are 0.58 MHz, 19.59 km, 0.92 MHz and 23.40 km, respectively. The results indicate that the discrepancies between these models are dependent on universal time, geographic latitude and seasons. The peak frequencies measured by COSMIC are generally larger than ionosonde’s observations in the nighttime or middle-latitudes with the amplitude of lower than 25%, while the averaged peak height derived from COSMIC is smaller than ionosonde’s data in the polar regions. The differences between ANN-based maps and references show that the discrepancies between two ionospheric detecting techniques are proportional to the intensity of solar radiation. Besides, a new method based on the ANN technique is proposed to reduce the discrepancies for improving ionospheric models developed by multiple measurements, the results indicate that the RMSEs of ANN models optimized by the method are 14–25% lower than the models without the application of the method. Furthermore, the ionospheric model built by the multiple measurements with the application of the method is more powerful in capturing the ionospheric dynamic physics features, such as equatorial ionization, Weddell Sea, mid-latitude summer nighttime and winter anomalies. In conclusion, the new method is significant in improving the accuracy and physical characteristics of an ionospheric model based on multi-source observations.     

Performance evaluation of IRI-2007 at equatorial latitudes and its Matlab version for GNSS applications

International Reference Ionosphere (IRI) model is the widely used empirical model for ionospheric predictions, especially TEC which is an important parameter for radio navigation and communication. The Fortran based IRI-2007 does not support real-time interactive visualization and debugging. Therefore, the source code is converted into Matlab and is validated for the purposes of this study. This facilitates easy representation of results and for near real-time implementation of IRI in the applications including spacecraft launching, now casting, pseudolite based navigation systems etc. In addition, the vertical delay results over the equatorial region derived from IRI and GPS data of three IGS stations namely Libreville (Garbon, Africa), Brasilia (Brazil, South America) and Hyderabad (India, Asia) are compared. As the IRI model does not account for plasmasphere TEC, the vertical delays are underestimated compared to vertical delays of GPS signals. Therefore, the model should be modified accordingly for precise TEC estimation.     

Response of the ionospheric F-region in the Brazilian sector during the super geomagnetic storm in April 2000 observed by GPS

The response of the ionospheric F-region in the equatorial and low latitude regions in the Brazilian sector during the super geomagnetic storm on 06–07 April 2000 has been studied in the present investigation. The geomagnetic storm reached a minimum D_st of −288 nT at 0100 UT on 07 April. In this paper, we present vertical total electron content (VTEC) and phase fluctuations (in TECU/min) from GPS observations obtained at Imperatriz (5.5°S, 47.5°W; IMPZ), Brasília (15.9°S, 47.9°W; BRAZ), Presidente Prudente (22.12°S, 51.4°W; UEPP), and Porto Alegre (30.1°S, 51.1°W; POAL) during the period 05–08 April. Also, several GPS-based TEC maps are presented from the global GPS network, showing widespread and drastic TEC changes during the different phases of the geomagnetic storm. In addition, ion density measurements on-board the satellite Defense Meteorological Satellite Program (DMSP) F15 orbiting at an altitude of 840 km and the first Republic of China satellite (ROCSAT-1) orbiting at an altitude of 600 km are presented. The observations indicate that one of the orbits of the DMSP satellite is fairly close to the 4 GPS stations and both the DMSP F15 ion-density plots and the phase fluctuations from GPS observations show no ionospheric irregularities in the Brazilian sector before 2358 UT on the night of 06–07 April 2000. During the fast decrease of D_st on 06 April, there is a prompt penetration of electric field of magnetospheric origin resulting in decrease of VTEC at IMPZ, an equatorial station and large increase in VTEC at POAL, a low latitude station. This resulted in strong phase fluctuations on the night of 06–07 April, up to POAL. During the daytime on 07 April during the recovery phase, the VTEC observations show positive ionospheric storm at all the GPS stations, from IMPZ to POAL, and the effect increasing from IMPZ to POAL. This is possibly linked to the equatorward directed meridional wind. During the daytime on 08 April (the recovery phase continues), the VTEC observations show very small negative ionospheric storm at IMPZ but the positive ionospheric storm effect is observed from BRAZ to POAL possibly linked to enhancement of the equatorial ionospheric anomaly.     

An alternative Shell inversion technique - Analysis and validation based on COSMIC and ionosonde data

Multi-channel Global Positioning System (GPS) carrier phase signals, received by the six low Earth orbiting (LEO) satellites from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) program, were used to undertake active limb sounding of the Earth’s atmosphere and ionosphere via radio occultation. In the ionospheric radio occultation (IRO) data processing, the standard Shell inversion technique (SIT), transformed from the traditional Abel inversion technique (AIT), is widely used, and can retrieve good electron density profiles. In this paper, an alternative SIT method is proposed. The comparison between different inversion techniques will be discussed, taking advantage of the availability of COSMIC datasets. Moreover, the occultation results obtained from the SIT and alternative SIT at 500 km and 800 km, are compared with ionosonde measurements. The electron densities from the alternative SIT show excellent consistency to those from the SIT, with strong correlations over 0.996 and 0.999 at altitudes of 500 km and 800 km, respectively, and the peak electron densities (N_mF₂) from the alternative SIT are equivalent to the SIT, with 0.839 vs. 0.844, and 0.907 vs. 0.909 correlation coefficients when comparing to those by the ionosondes. These results show that: (1) the N_mF₂ and h_mF₂ retrieved from the SIT and alternative SIT are highly consistent, and in a good agreement with those measured by ionosondes, (2) no matter which inversion technique is used, the occultation results at the higher orbits (∼800 km) are better than those at the lower orbits (∼500 km).