The confluence of recent instrumentation deployments in Africa with developments for the determination of plasmasphere electron content using Global Positioning System (GPS) receivers has provided new opportunities for investigations in that region. This investigation, using a selected chain of GPS stations, extends the method (SCORPION) previously applied to a chain of GPS stations in North America in order to separate the ionosphere and plasmasphere contributions to the total electron content (TEC) during a day (24 July) in 2011. The results span latitudes from the southern tip of Africa, across the Equator, to the southern Arabian Peninsula, providing a continuous latitudinal profile for both the ionosphere and plasmasphere during this day.
The peak diurnal vertical ionosphere electron content (IEC) increases from
about 14 TEC units (1 TEC unit
Comparisons to other GPS data analyses are also presented for TEC, indicating the influence of the PEC on the determination of latitudinal TEC variations and also on the absolute TEC levels, by inducing an overestimate of the receiver bias. The derived TEC latitudinal profiles, in comparison to global map profiles, tend to differ from the map results only about as much as the map results differ among themselves. A combination of ionosonde IEC and alternative GPS TEC measurements, which in principle permits a PEC determination through their difference, was compared to the composite and separate ionosphere and plasmasphere contributions derived solely by the SCORPION method for one station. Although there is considerably more scatter in the PEC values derived from the difference of the GPS TEC and ionosonde IEC measurements compared to the PEC values derived by the SCORPION method, the average overhead values for this day are comparable for the two methods, near 2 TEC units, at the South African site examined.
This initial investigation provides a basis for day-to-day TEC monitoring for Africa, with separate ionosphere and plasmasphere electron content determinations.
The SCORPION technique for calibration of Global Positioning System (GPS) total electron content (TEC) measurements was developed primarily to address the detrimental effects of the plasmasphere electron content (PEC) on the previous SCORE calibration technique (Lunt et al., 1999b), necessitating the development of a parametric representation for the plasmasphere as part of the calibration process (Mazzella et al., 2002, 2007). Previous studies included validations of the SCORPION technique for model data and corroborations against other methods for actual TEC measurements (Mazzella et al., 2007), and a more recent study demonstrated a method for resolving the remaining ambiguity between the GPS receiver bias and an omnipresent uniform PEC contribution (Mazzella, 2012). That method requires a contiguous latitudinal chain of GPS stations, permitting the alignment of vertical PEC between pairs of stations to determine the locally omnipresent uniform PEC contribution (the “plasmasphere baseline”).
Deployments of GPS receivers (and other ionospheric measurement systems) in
Africa have significantly increased since 2008 (Yizengaw et al., 2013),
permitting application of the same SCORPION analysis for a chain of stations
on that continent. As with the previous North America study (Mazzella,
2012), a strict longitudinal alignment of stations is not required, because
of the reliance on the local time dependence for the variations of the
ionosphere electron content (IEC) and PEC, but a reasonable bound is
required for the longitudinal differences between stations to minimize the
effects of intrinsic variations for the ionosphere and plasmasphere. An
initial bound of about 5
The sites used for the Africa chain study, from south to north, with their supporting networks, and their derived plasmasphere baseline values (in TEC units).
The GPS stations used for the Africa TEC profile for day 205, 2011,
with their associated ionosphere coverage regions (ovals), corresponding to
an elevation threshold of 35
An initial selection of GPS stations at strategic locations along the
prospective latitudinal chain was examined for general data quality,
including continuity and multipath. Replacement and supplementary stations
were included to establish a complete chain with overlapping coverage for
the ionosphere, using a 35
Ionosphere reference altitudes, used for the determination of the ionospheric penetration point (IPP) and the slant-to-vertical conversion for TEC, for each of the stations in the chain.
The final selection of GPS stations used for this study is displayed in Fig. 1 and listed in Table 1 with their network associations. The ovals
associated with each station indicate the regional coverage, for an
elevation threshold of 35
Prior to any extensive processing, an evaluation was conducted for two
parameters that have a significant role for the calibration process. These
parameters are the effective altitude for the (thin-layer) ionosphere, used
for the calculation of the ionospheric penetration point (IPP) coordinates,
and the slant factor altitude, used for the conversion from slant IEC to
equivalent vertical IEC. As noted by Mazzella (2009), these two altitude
quantities need not be equal. The evaluation of the two altitude quantities
was conducted using the Sheffield University Plasmasphere Ionosphere Model
(SUPIM; Bailey and Balan, 1996) for five of the stations in the chain
(HNUS, MFKG, TEZI, MOIU, and FRSN) for geophysical parameters typical of the
conditions during the period selected (F10.7
A notable feature in Fig. 2 is that the IPP altitude becomes greater than the slant factor altitude for TEZI and equatorward stations. This occurrence is attributed to the local TEC gradients and the limitations of the method for the altitude determinations, but contradicts the conjecture by Mazzella (2009) that the difference between the slant factor altitude and the IPP altitude can be interpreted as an effective thickness for the ionosphere.
The GPS receiver at Hermanus, South Africa (HNUS), was used for a survey of
the ionosphere and plasmasphere conditions for the entire period from 7
July to 4 August 2011, except for 20 July 2011, for which data were
unavailable. The most significant PEC occurred for 23 July 2011 (day 204),
but examination of other stations indicated active conditions for the
ionosphere, reducing both the accuracy of the calibrations and the
partitioning of TEC between the ionosphere and plasmasphere. Consequently,
the data were analyzed for 24 July 2011 (day 205), for which the PEC values
were only slightly lower and quiet ionospheric conditions (Ap
The initial processing of the RINEX (Receiver INdependent EXchange format) data files for each station was performed using the GPS Toolkit developed by the Applied Research Laboratory (ARL) of the University of Texas (Tolman et al., 2004), using data for two contiguous days to obtain complete satellite passes spanning a period of more than 24 h (Mazzella, 2012). Further examination and editing were conducted for each satellite pass to eliminate outliers for dispersive carrier phase and dispersive group delay, remove any residual phase discontinuities, and eliminate data segments affected by missing ephemeris data. The standard practice of aligning dispersive carrier phase to dispersive group delay was performed (e.g., Ma and Maruyama, 2003), without elevation-dependent weighting, followed by additional corrections for multipath to minimize bias discrepancies for different passes by the same satellite (Kee and Parkinson, 1994; Andreasen et al., 2002).
The diurnal TEC variations for the GPS stations, displaying, for each station, the magnetic latitude of the IPP in the bottom panel, the equivalent vertical ionospheric electron content (VIEC) in the middle panel, and the slant plasmaspheric electron content (SPEC) in the top panel, versus the magnetic local time at the IPP. Data segments are color-coded by GPS satellite.
All GPS stations were initially calibrated by SCORPION individually, but
preliminary inter-comparisons for adjacent stations indicated discrepancies
of up to 5 TEC units (1 TEC unit
The TEC determinations for all of the stations are presented in Fig. 3, displaying, from bottom to top panels, the magnetic latitude (MLAT) of the IPP, the equivalent vertical ionospheric electron content (VIEC), and the slant plasmaspheric electron content (SPEC) against the magnetic local time (MLT) at the IPP. As in the previous studies by Mazzella (2009, 2012), the boundary between the ionosphere and plasmasphere is considered to be at an altitude of 1000 km, for all latitudes, and is thus distinct from a boundary based on geophysical quantities. The peak VIEC increases significantly between TDOU and NURK, by more than a factor of 2, while the SPEC also increases somewhat. The latitudinal variations of SPEC (Fig. 3), both for each GPS station and across all stations, display characteristics that overlap with the previous North America study (Mazzella, 2012), but are distinctly different in the equatorial region. The smoothness of the PEC results is consistent with the large volume and diffuse density distribution of the plasmasphere.
The TEC comparisons for overlapping 1
The TEC comparisons for overlapping magnetic latitude bands, used to establish the plasmasphere baseline for each station except HNUS, are displayed in Fig. 4. Here, the vertical plasmaspheric electron content (VPEC) is displayed in the top section of each panel for each station. This quantity is not directly measured, but is determined by the parametric representation for the plasmasphere, as derived by SCORPION.
The correspondence for VIEC values generally reflects the usage of the referenced calibration method, but the comparisons between HNUS and SUTH are for independently calibrated stations. The plasmasphere baseline values, derived progressively beginning from HNUS, have been applied prior to generating this display, and eliminated any relative displacements of the paired VPEC profiles. These plasmasphere baseline values are presented in Table 1. The assignment of zero for the plasmasphere baseline for HNUS was confirmed by extrapolation from KERG, but is also consistent with the zero plasmasphere baseline determined for stations at similar magnetic latitudes for the North American chain (Mazzella, 2012).
Latitudinal vertical TEC profiles for the separate ionosphere (bottom panel) and plasmasphere (top panel), with the composite vertical TEC (middle panel), for the SCORPION results and the CODE, ESA, JPL, and UPC IONEX maps, for magnetic local times (clockwise from upper left) 06:00, 12:00, 18:00, and 24:00.
Receiver biases and error estimates, in TEC units, from SCORPION, IONEX maps, and the July 2011 CODE monthly report, indicating whether the coarse acquisition (C1) or precise code (P1) signal is used for the L1 frequency.
A further observation from Fig. 4 arises from the occurrence of vertical TEC (VTEC) profiles for southern and northern regions relative to an individual station in successive plots, indicating the transformation of the diurnal profile with latitude. The most significant change occurs for NURK, but UKAM and MOIU also display significant changes. These TEC changes with latitude are directly displayed in Fig. 5, for four magnetic local times (06:00, 12:00, 18:00, and 24:00, clockwise from upper left). For each local time, the VIEC appears in the bottom panel, the VPEC appears in the top panel, and their sum, the composite VTEC, appears in the center panel. Also displayed in these figures are IONEX TEC results from the Center for Orbit Determination in Europe (CODE) of the Astronomical Institute of the University of Bern, the European Space Agency (ESA), the Jet Propulsion Laboratory (JPL) of the California Institute of Technology, and the Research Group of Astronomy and Geomatics of the Technical University of Catalonia (UPC) analysis centers. The IONEX TEC results correspond to the composite VTEC, but are also displayed with the vertical ionosphere TEC for comparison.
In Fig. 5, the GPS data were selected over an interval of
Overlays of SCORPION vertical TEC (blue) and BC/ISR vertical TEC (ionosphere and plasmasphere) (red), for SCORPION VIEC (ionosphere only) in the left column, CEQV (composite equivalent vertical ionosphere and plasmasphere) in the center column, and CVTC (composite vertical ionosphere and plasmasphere) in the right column.
A review of the IONEX GPS station locations indicated an absence of coverage
for the African latitude band between about 10 and
15
Based on both the SCORPION and IONEX results, the equatorial ionosphere
anomaly is absent for this day over Africa, although the peak ionosphere
VTEC does occur as far south as 10
For the nighttime case (MLT
For a comparison of TEC results between SCORPION and a more conventional GPS
TEC determination method, the Boston College Institute for Scientific
Research (BC/ISR) calibration technique was chosen as an alternative
(Valladares et al., 2009). For the latter technique, the uncalibrated
dispersive, phase-aligned TEC measurements are corrected for the relative
GPS satellite biases reported by the University of Bern (
Four stations (HNUS, TEZI, MOIU, and FRSN) within the chain were selected for this comparison, and the quantities compared included the BC/ISR equivalent vertical TEC (VTEC) against the ionosphere equivalent vertical TEC (VIEC), the composite equivalent vertical TEC (CEQV), and the composite ionosphere and plasmasphere vertical TEC (CVTC) determined by SCORPION. Of the SCORPION quantities, the CEQV is closest in definition to the BC/ISR VTEC, for the conversion of the combined slant plasmasphere and ionosphere TEC to equivalent vertical, but differs through the biases and altitudes associated with the IPP and slant factor definitions. The CVTC from SCORPION is the closest quantity to true VTEC, being calculated as the combination of the VIEC and the SCORPION representation for the local vertical PEC (VPEC) at the IPP. (The VIEC, CVTC, and VPEC were previously displayed in Fig. 5, across latitudes.)
For the comparison of the BC/ISR VTEC against the SCORPION VIEC (Fig. 6, left column), the VIEC diurnal profiles for HNUS, TEZI, and FRSN are distinctly smoother than the BC/ISR VTEC profiles, with a lesser difference for MOIU. This difference is attributable to the plasmasphere, as noted by Mazzella (2009), and is further supported by the comparison of the BC/ISR VTEC against the SCORPION CEQV (Fig. 6, middle column), for which the variations at each local time are quite similar and the TEC differences are primarily for bias levels. A lesser variation for the SCORPION CVTC against the BC/ISR VTEC is evident (Fig. 6, right column), indicating the role of the IPP reference location for the plasmasphere vertical electron content determination in ameliorating the apparent gradients induced by converting between slant and vertical PEC using the same conversion as for the ionosphere.
An additional comparison for alternative measurements was provided by the
ionosonde and GPS receiver (GRHM) at Grahamstown, South Africa (McKinnell et
al., 2007). (Data from the ionosonde at Hermanus were unavailable for day
205, 2011.) The method utilized was previously described and implemented by
Belehaki et al. (2003), in which a topside extension of ionosonde
measurements to 1000 km altitude is used to determine the VIEC (Huang and
Reinisch, 2001) while the GPS receiver provides the combined IEC plus PEC.
The TEC data for GRHM were processed by SCORPION, and separately using the
BC/ISR method, and the latter results were combined with the ionosphere TEC
(ITEC) estimates from the ionosonde, after selecting GPS data for elevations
at or above 60
The derived relative satellite biases, resulting from the removal
of the receiver bias, separately for the receivers for which the coarse
acquisition (C1) signal is used
The receiver biases for each station were calculated as the average bias over all satellites, and are listed in Table 2, with corresponding values (if available) from tabulations in the IONEX files or reported by CODE. The ESA IONEX receiver bias for NURK is significantly different from the SCORPION, CODE IONEX, and CODE monthly report biases, and would produce negative TEC values if applied to the measurements for that station. Similar receiver bias values are reported in the ESA IONEX files for adjacent days, so this does not appear to be an isolated anomalous value, but may be indicative of the plasmasphere influence on the ESA calibration process. Such an overestimate of the biases for equatorial regions was noted by Mazzella et al. (2007), Anghel et al. (2009), and Carrano et al. (2009), for different calibration methods that did not incorporate determinations of the plasmasphere electron content.
The associated derived relative satellite biases, resulting from the removal
of the receiver bias, are displayed in Fig. 8, separately for the receivers
for which the coarse acquisition (C1) signal is used (Fig. 8a) instead of
the precise code (P1) signal (Fig. 8b). For comparison, the relative
satellite biases reported by CODE for July 2011 are also displayed. The
agreement among the stations and with CODE is generally within 3 TEC units
(with 1 TEC unit corresponding to 0.350 ns of differential delay),
although there are several discrepant relative bias differences exceeding
that bound, notably for NURK, a site for which the multipath was more
significant than most of the other sites. The PRN 1 CODE P1-P2 bias (23.0503
TEC units) for July 2011 is distinctly different from the CODE value
(30.6156 TEC units) reported for August 2011, possibly reflecting the
limited time for evaluation after the launch of PRN 1 on 16 July 2011
(
A quiet-day study for a latitudinal chain in Africa was conducted as an initial investigation for that region using the SCORPION method for partitioning the ionosphere and plasmasphere electron content contributions, both to improve the quality of the TEC calibrations and to evaluate more representative values for the VTEC. Standard procedures for converting STEC to equivalent VTEC utilize the combined slant IEC and PEC contributions along the line of sight. However, except for nearly vertical lines of sight, this PEC contribution is from a more distant region than the IPP, so the horizontal PEC gradient is compressed into a more steeply varying quantity by mapping this contribution to the IPP (Mazzella, 2009). This gradient enhancement appears as a “sawtooth” pattern of steeper, disjoint TEC gradients distinct from the general TEC trend across several sites (Mazzella, 2012). The absence of this sawtooth pattern in the composite VTEC latitudinal profiles (Fig. 5) is an indication of the quality of the current results. The smooth variation of the vertical PEC with latitude indicates the consistency of the regional plasmasphere determinations at each site, as enforced by both the referenced calibrations and the plasmasphere baseline matching process.
Although referenced calibrations were not required for the North America chain (Mazzella, 2012), and plasmasphere baseline adjustments were marginal for that investigation, both of those features of the SCORPION method had significant roles for the Africa chain, because of the nature of the ionosphere and plasmasphere at lower latitudes. Based on SUPIM studies, Mazzella (2009) indicated the possibility of different effective altitudes for the ionospheric penetration point and slant factor conversion at different local times, but similar evaluations for lower latitudes (also using SUPIM) indicate the prospect of varying effective altitudes with azimuth for a single local time, even with the exclusion of the plasmasphere contribution. In combination with the tasks of disentangling the bias effects and plasmasphere contributions from the ionosphere electron content, the possible variation for the slant factor conversion is an additional complication, although the apparent resolution by means of the referenced calibrations is somewhat encouraging. A caveat regarding the results is the possible progressive distortion of the TEC profiles along the chain of sites, but the individual diurnal and latitudinal profiles appear to indicate relatively restrictive limits regarding this possible occurrence.
A further favorable result in this study was the avoidance of negative TEC excursions for the equatorial region, as a consequence of the influence on the bias determinations by the plasmasphere when it is not explicitly treated within the calibration process (Mazzella et al., 2002, 2007; Anghel et al., 2009; Carrano et al., 2009; Mazzella, 2012). Consequently, the associated latitudinal trend for the receiver bias error (Mazzella, 2012) has also been significantly reduced. Accommodating latitudinally varying plasmasphere content for the accuracy of the calibration process provides some encouragement that the day-to-day temporal variation of the plasmasphere content will also have less influence on the derived receiver biases, so that a possibly spurious variability of receiver bias values can be eliminated.
The prospect for day-to-day TEC monitoring, with ongoing bias determinations and partitioning between the ionosphere and plasmasphere, requires addressing the detrimental effects of ionospheric activity on the calibration process. Such ionosphere activity invalidates the assumption of equal equivalent vertical IEC values for “common” IPP observations (with the operational definition of “common” being implemented by the parameters defining conjunctions). Use of relative (satellite) biases within the calibration process reduces the number of parameters to be determined by the calibration, allowing more restrictive conjunction limits for comparisons of equivalent vertical IEC values, especially for limiting universal time differences between conjunction pairs, thus reducing the time allowed for intrinsic ionospheric variations. For GPS stations that have previously been fully calibrated on a quiet day, as in this study, the set of relative satellite biases for the station may be more advantageous for use than a published globally derived set of relative satellite biases (such as those from CODE), because some effects ameliorating residual multipath will be advantageously incorporated. Using relative biases may also allow the activity occurrences within a calibration period to be excluded, while still allowing subsequent use of the TEC data for the active interval. However, further investigation and experiments are required for such an implementation, because of the uncertain effects on the partitioning between the ionosphere and plasmasphere.
Addressing the direct and indirect influences of the plasmasphere on TEC measurements for Africa is a significant development, for improving both the quality and the interpretation of the measurements. The influence of the plasmasphere content on the bias determinations and derived IEC latitudinal gradients has been noted previously (Mazzella et al., 2007; Mazzella, 2012), but the gradual replenishment and storm-related depletion of the plasmasphere produce circumstances for masquerading as seasonal or other medium-term effects, unless addressed by additional investigations, such as the superposed epoch analysis by Lunt et al. (1999c).
The current study included comparisons of SCORPION results against a more conventional single-station GPS TEC analysis method (BC/ISR), as well as against multi-station, global mapping TEC representations (IONEX) for four separate analysis centers (CODE, ESA, JPL, and UPC). The SCORPION comparisons to the BC/ISR results displayed the anticipated difference in VTEC levels between calibration methods with explicit or implicit treatments of the plasmasphere content, with the associated bias overestimate at lower latitudes when the plasmasphere content is not explicitly determined (Mazzella et al., 2007). The VTEC variations with latitude are also embellished when the plasmasphere content is not treated separately. The SCORPION CVTC latitudinal profiles, in comparison to the IONEX profiles, tend to differ from the IONEX results only about as much as the IONEX results differ among themselves, although the CODE and UPC results appear to match more closely with SCORPION, especially for magnetic local times between 12:00 and 16:00.
The ionosonde-with-GPS method (Belehaki et al., 2003) provides an alternative ground-based technique for separate ionosphere and plasmasphere electron content determinations, but requires both an accurate scale height determination from the ionosonde measurements and an accurate calibration for the GPS TEC in the presence of plasmasphere content. The large (5 TEC units) and irregular variations in the plasmasphere electron content derived by this method (Fig. 7b) indicate difficulties with either or both of those instrumental parameter determinations. The previous discussion concerning the GPS TEC determinations would indicate that the problem is not solely with that measurement. The discrepancies between the ionosonde TEC and GPS TEC noted by Fremouw et al. (2001) and Secan et al. (2005) for higher northern latitudes, where the plasmasphere contribution is insignificant, indicate considerable uncertainty for the ionosonde TEC determinations.
A remaining data source for PEC comparison data is the set of measurements from the GPS receiver aboard the JASON-2 satellite. (Comparisons of IEC, using the JASON-2 altimeter data, are considerably restricted for the chain of stations used in the current study, because of their inland locations.) A supplementary analysis similar to that by Yizengaw et al. (2008) and Lee et al. (2013) is being conducted specifically for Africa on the day of this case study to evaluate PEC values for comparison to the SCORPION results.
Data associated with Figs. 2–8 are provided as text tabulations in the Supplement.
Altitude for IPP location (HgtIPP):
Generate azimuth slant TEC (STEC) scans from a model, at 2 h
intervals, for a 35 Determine the median altitudes, using the cumulative STEC altitude
profile for the azimuth scans. For each 2 h interval, calculate an IPP altitude as an STEC-weighted
average of the median altitudes over all azimuths. Calculate the representative IPP altitude (HgtIPP) as the maximum of the
individual azimuthal STEC-weighted average IPP altitudes. (For calibrations,
this was considered to be the most appropriate choice, compared to a median
or average value.) Altitude for slant factor (AltIPP):
Generate meridian STEC scans from a model, for elevations from
0 to 90 Generate latitudinal true VTEC profiles from the model, for latitudes
between For elevations above a threshold (35 to solve for the slant factor altitudes (AltIPP) for every 2 h, for
Select the median value of all of the 2 h slant conversion altitudes
as the daily slant conversion altitude (AltIPP).
These processes are supplemented by various plot evaluations.
A. J. Mazzella Jr. is a developer for the SCORPION method. The other authors declare that they have no conflict of interest.
Part of E. Yizengaw's work was supported by AFOSR (FA9550-12-1-0437 and FA9550-15-1-0399) and NSF (AGS1450512) grants. Graham J. Bailey provided SUPIM. Development of SCORPION was conducted in
collaboration with G. Susan Rao and supported by the Air Force Research
Laboratory Space Vehicles Directorate under SBIR contracts FA8718-04-C-0009
and FA8718-05-C-0026 to NorthWest Research Associates. The GPS data were
obtained from the Crustal Dynamics Data Information System (CDDIS) at the
NASA Goddard Space Flight Center FTP site (