Anthony J. Mannucci

A global mapping technique for GPS‐derived ionospheric total electron content measurements

A worldwide network of receivers tracking the transmissions of Global Positioning System (GPS) satellites represents a new source of ionospheric data that is globally distributed and continuously available. We describe a technique for retrieving the global distribution of vertical total electron content (TEC) from GPS-based measurements. The approach is based on interpolating TEC within triangular tiles that tessellate the ionosphere modeled as a thin spherical shell. The high spatial resolution of pixel-based methods, where widely separated regions can be retrieved independently of each other, is combined with the efficient retrieval of gradients characteristic of polynomial fitting. TEC predictions from climatological models are incorporated as simulated data to bridge significant gaps between measurements. Time sequences of global TEC maps are formed by incrementally updating the most recent retrieval with the newest data as it becomes available. This Kalman filtering approach smooths the maps in time, and provides time-resolved covariance information, useful for mapping the formal error of each global TEC retrieval. Preliminary comparisons with independent vertical TEC data, available from the TOPEX dual-frequency altimeter, suggest that the maps can accurately reproduce spatial and temporal ionospheric variations over latitudes ranging from equatorial to about ±65°.

Dayside global ionospheric response to the major interplanetary events of October 29--30, 2003 ''Halloween Storms''

We demonstrate extreme ionospheric response to the large interplanetary electric fields during the Halloween storms that occurred on October 29 and 30, 2003. Within a few (2-5) hours of the time when the enhanced interplanetary electric field impinged on the magnetopause, dayside total electron content increases of ∼40% and ∼250% are observed for the October 29 and 30 events, respectively. During the Oct 30 event, ∼900% increases in electron content above the CHAMP satellite (∼400 km altitude) were observed at mid-latitudes (±30 degrees geomagnetic). The geomagnetic storm-time phenomenon of prompt penetration electric fields is a possible contributing cause of these electron content increases, producing dayside ionospheric uplift combined with equatorial plasma diffusion along magnetic field lines to higher latitudes, creating a daytime super-fountain effect.

Monitoring of global ionospheric irregularities using the Worldwide GPS Network

A prototype system has been developed to monitor the instantaneous global distribution of ionospheric irregularities, using the worldwide network of Globa Positioning System (GPS) receivers. Case studies in this pape indicate that GPS receiver loss of lock of signal tracking may be associated with strong phase fluctuations. It is shown that a network-based GPS monitoring system will enable us to study the generation and evolution of ionospheric irregularities continuously around the globe under various solar and geophysical conditions, which is particularly suitable for studies of ionospheric storms, and for space weather research and applications.

Global ionosphere perturbations monitored by the Worldwide GPS Network

For the first time, measurements from the Global Positioning System (GPS) worldwide network are employed to study the global ionospheric total electron content (TEC) changes during a magnetic storm (November 26, 1994). These measurements are obtained from more than 60 world-wide GPS stations which continuously receive dual-frequency signals. Based on the delays of the signals, we have generated high resolution global ionospheric maps (GIM) of TEC at 15 minute intervals. Using a differential method comparing storm time maps with quiet time maps, we find that significant TEC increases (the positive effect) are the major feature in the winter hemisphere during this storm (the maximum percent change relative to quiet times is about 150%). During this particular storm, there is almost no negative phase. A traveling ionospheric disturbance (TID) event is identified that propagates from the northern subauroral region to lower latitudes (down to about 30°N) at a speed of ∼460 m/s. This TID is coincident with significant increases in the TEC around the noon sector. We also find that another strong TEC enhancement occurs in the pre-dawn sector in the northern hemispheric subauroral latitudes, in the beginning of the storm main phase. This enhancement then spreads into almost the entire nightside. The nighttime TEC increase in the subauroral region is also noted in the southern hemisphere, but is less significant. These preliminary results indicate that the differential mapping method, which is based on GPS network measurements, appears to be a powerful tool for studying the global pattern and evolution process of the entire ionospheric perturbation.

Subdaily Northern Hemisphere Ionospheric Maps Using an Extensive Network of GPS Receivers

Ionospheric total electron content (TEC) data derived from dual-frequency Global Positioning System (GPS) signals from 30 globally distributed network sites are fit to a simple ionospheric shell model, yielding a map of the ionosphere in the northern hemisphere every 12 hours during the January 1–15, 1993 period, as well as values for the satellite and receiver instrumental biases. Root-mean-square (RMS) residuals of 2–3 TEC units are observed over the 20°–80° latitude band. Various systematic errors affecting the TEC estimates are discussed. The capability of using these global maps to produce ionospheric calibrations for sites at which no GPS data are available is also investigated.

Ionospheric total electron content perturbations monitored by the GPS global network during two northern hemisphere winter storms

The global evolution of two major ionospheric storms, occurring on November 4, 1993, and November 26, 1994, respectively, is studied using measurements of total electron content (TEC) obtained from a worldwide network of ground-based GPS receivers. The time-dependent features of ionospheric storms are identified using TEC difference maps based on the percent change of TEC during storm time relative to quiet time. The onset of each ionospheric storm is indicated by the appearance of auroral/subauroral TEC enhancements which occur within 1 hour of the beginning of the geomagnetic storm main phase. Significant TEC enhancements (> 100%) are observed in the winter northern hemisphere. The rate at which TEC enhancements appear is found to correlate with gradients in the Dst index. The large scale ionospheric structures identified during the storms are (1) nightside auroral/subauroral enhancements which surround the auroral oval, (2) dayside (around noon) high-latitude and middle-latitude enhancements associated with traveling ionospheric disturbances, and (3) conjugate latitudinal enhancements. For the November 1993 storm, a short positive phase (about 15 hours) is followed by a long negative phase (∼60 hours). In the November 1994 storm, we have identified the clear signature of a traveling ionospheric disturbance (TID) which propagated at a speed of ∼460 m/s from ∼60° N to ∼40° N. The motion of this disturbance appears to conserve angular momentum.

Detecting ionospheric TEC perturbations caused by natural hazards using a global network of GPS receivers: The Tohoku case study

Recent advances in GPS data processing have demonstrated that ground-based GPS receivers are capable of detecting ionospheric TEC perturbations caused by surface-generated Rayleigh, acoustic and gravity waves. There have been a number of publications discussing TEC perturbations immediately following the M 9.0 Tohoku earthquake in Japan on March 11, 2011. Most investigators have focused on the ionospheric responses up to a few hours following the earthquake and tsunami. In our research, in addition to March 11, 2011 we investigate global ionospheric TEC perturbations a day before and after the event. We also compare indices of geomagnetic activity on all three days with perturbations in TEC, revealing strong geomagnetic storm conditions that are also apparent in processed GEONET TEC observations. In addition to the traveling ionospheric disturbances (TIDs) produced by the earthquake and tsunami, we also detect “regular” TIDs across Japan about 5 hours following the Tohoku event, concluding these are likely due to geomagnetic activity. The variety of observed TEC perturbations are consistent with tsunami-generated gravity waves, auroral activity, regular TIDs and equatorial fluctuations induced by increased geomagnetic activity. We demonstrate our capabilities to monitor TEC fluctuations using JPL’s real-time Global Assimilative Ionospheric Model (GAIM) system. We show that a real-time global TEC monitoring network is able to detect the acoustic and gravity waves generated by the earthquake and tsunami. With additional real-time stations deployed, this new capability has the potential to provide real-time monitoring of TEC perturbations that could potentially serve as a plug-in to enhance existing early warning systems.

A comparative study of ionospheric total electron content measurements using global ionospheric maps of GPS, TOPEX radar, and the Bent model

Global ionospheric mapping (GIM) is a new and emerging technique for determining global ionospheric TEC (total electron content) based on measurements from a worldwide network of Global Positioning System (GPS) receivers. In this study, GIM accuracy in specifying TEC is investigated by comparison with direct ionospheric measurements from the TOPEX altimeter. A climatological model (Bent model) is also used to compare with the TOPEX altimeter data. We find that the GIM technique has much better agreement with TOPEX in TEC measurements, compared with the predictions of the climatological model. The difference between GIM and TOPEX in TEC measurements is very small (less than 1.5 TEC units (TECU)) within a 1500-km range from a reference GPS station. The RMS gradually increases with increasing distance from the station, while the Bent model shows a constant large RMS, unrelated to any station location. Within a 1000-km distance of a GPS site (elevation angle > 25°), GIM has a good correlation (R > 0.93) to TOPEX with respect to TEC measurements. The slope of the linear fitting line to the data set from two TOPEX cycles is 44.5° (near the ideal 45°). In the northern hemispheric regions, ionospheric specification by GIM appears to be accurate to within 3-10 TECU up to 2000+ km away from nearest GPS station (corresponding to ∼1° elevation angle cutoff). Beyond 2000 km, GIM accuracy, on average, is reduced to the Bent model levels. In the equatorial region, the Bent model predictions are systematically lower (∼5.0 TECU) than TOPEX values and often show a saturation at large TEC values. During ionospheric disturbed periods, GIM sometimes shows differences from TOPEX values due to transient variations of the ionosphere. Such problems may be improved by the continuous addition of new GPS stations in data-sparse regions. Thus, over a GPS stations measurement realm (up to 2000 km in radius), GIM can produce generally accurate TEC values. Through a spatial and temporal extrapolation of GPS-derived TEC measurements, the GIM technique provides a powerful tool for monitoring global ionospheric features in near real time.

Space Weather Observations by GNSS Radio Occultation: From FORMOSAT-3/COSMIC to FORMOSAT-7/COSMIC-2

The joint Taiwan-United States FORMOSAT-3/COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate) mission, hereafter called COSMIC, is the first satellite constellation dedicated to remotely sense Earths atmosphere and ionosphere using a technique called Global Positioning System (GPS) radio occultation (RO). The occultations yield abundant information about neutral atmospheric temperature and moisture as well as space weather estimates of slant total electron content, electron density profiles, and an amplitude scintillation index, S4. With the success of COSMIC, the United States and Taiwan are moving forward with a follow-on RO mission named FORMOSAT-7/COSMIC-2 (COSMIC-2), which will ultimately place 12 satellites in orbit with two launches in 2016 and 2019. COSMIC-2 satellites will carry an advanced Global Navigation Satellite System (GNSS) RO receiver that will track both GPS and Russian Global Navigation Satellite System signals, with capability for eventually tracking other GNSS signals from the Chinese BeiDou and European Galileo system, as well as secondary space weather payloads to measure low-latitude plasma drifts and scintillation at multiple frequencies. COSMIC-2 will provide 4–6 times (10–15X in the low latitudes) the number of atmospheric and ionospheric observations that were tracked with COSMIC and will also improve the quality of the observations. In this article we focus on COSMIC/COSMIC-2 measurements of key ionospheric parameters.

CEDAR Electrodynamics Thermosphere Ionosphere (ETI) Challenge for systematic assessment of ionosphere/thermosphere models: NmF2, hmF2, and vertical drift using ground‐based observations

[1] Objective quantification of model performance based on metrics helps us evaluate the current state of space physics modeling capability, address differences among various modeling approaches, and track model improvements over time. The Coupling, Energetics, and Dynamics of Atmospheric Regions (CEDAR) Electrodynamics Thermosphere Ionosphere (ETI) Challenge was initiated in 2009 to assess accuracy of various ionosphere/thermosphere models in reproducing ionosphere and thermosphere parameters. A total of nine events and five physical parameters were selected to compare between model outputs and observations. The nine events included two strong and one moderate geomagnetic storm events from GEM Challenge events and three moderate storms and three quiet periods from the first half of the International Polar Year (IPY) campaign, which lasted for 2 years, from March 2007 to March 2009. The five physical parameters selected were NmF2 and hmF2 from ISRs and LEO satellites such as CHAMP and COSMIC, vertical drifts at Jicamarca, and electron and neutral densities along the track of the CHAMP satellite. For this study, four different metrics and up to 10 models were used. In this paper, we focus on preliminary results of the study using ground-based measurements, which include NmF2 and hmF2 from Incoherent Scatter Radars (ISRs), and vertical drifts at Jicamarca. The results show that the model performance strongly depends on the type of metrics used, and thus no model is ranked top for all used metrics. The analysis further indicates that performance of the model also varies with latitude and geomagnetic activity level.