Determination of Total Electron Content at Equatorial Region—Thailand Using Radio Occultation Technique ()
1. Introduction
Ionosphere is the upper layer of the Earth atmosphere which has a great impact on signal from telecommunication system. It consists of free electrons and ions that are mostly concentrated in the F-region and produced mainly by solar radiation [1] . It causes delay of electromagnetic signals passing through it being a dispersive medium [2] . Ionosphere influences GPS signals mostly in the equatorial region within latitude of 20˚N and 20˚S which results in error measurements and hence, the need to analyze the variation of the total electron content in the ionosphere. Radio occultation data processing system (see Figure 1(a)) in essence was developed for the purpose of numerical weather prediction and climate benchmarking [3] . Right from the pioneer GPS/MET mission, low earth orbit (LEO) based GNSS radio occultation (RO) technique has been a powerful technique in ionosphere monitoring. Subsequently, many LEO satellites were launched with radio occultation payload, which included: CHAMP, GRACE, SACC/D, COSMIC, C/NOFS, Metop-A/B, and TerraSAR-X/TanDEM-X [4] [5] .
The fundamental problem addressed in this paper although very difficult is the ionospheric electron density and TEC fluctuations monitoring. Determining the amount of free electrons present within the ionospere over Thailand was carried out and then integrated to obtain total electron content. It shows to be more violent during the summer for the period of four years. This study is vital to space applications related organizations such as aviation, military and space weather, for numerical weather prediction and security purposes. The paper outline is as follows: Section 2 presents the data description and methodology. Section 3 analyzes the results and finally Section 4 is the conclusion.
2. Data Description and Methodology
2.1. Data Description
The Global Positioning System Ionospheric Radio Occultations are obtained from the COSMIC Data Analysis and Archive Centre (CDAAC) website. The ionospheric level 12 operational data products are made available within 24 hours by UCAR/CDAAC. The data products output used in this paper are ionospheric profiles of electron density (ionprf) whose accuracy is generally about 104 - 105 cm−3 and in NetCDF format. Observation shows that, cycle slips may affect some of these profiles [6] . UCAR/CDAAC is the only primary data processing centre particularly dedicated for COSMIC mission. This data is made available to the public both ftp and http server through CDAAC. One of the most important radio occultation products is the electron density profile for ionosphere. It is achieved from the calculation of excess phase file for each occultation with the aid of Abel inversion.
2.2. Methodology
COSMIC satellite is provided with four antennas each with which ionospheric electron density measurements is taken. Two antennas are used one for rising and one for setting occultation during this period. The GPS satellite receiver onboard LEO can track up to 13 GPS satellites through the two antennas and measure the difference GPS phase data on L1 and L2 every second. The bend angle, α (of the satellite or radio signals) obtained through
(a) (b)
Figure 1. Schematic diagram illustrating radio occultation, geometric of ray, bending angle α and impact parameter a [10] .
Abel inversion method is used to compute the refractivity. The radio occultation electron density profile retrieval approach is based on a few assumptions and approximations such as straight line signal propagation, coplanar of GNSS ray, spherical symmetry of electron density, circular LEO orbit and first order estimation of orbital electron density. Even though radio occultation technique is characterized with high precision, provides global profiles and maps of the atmospheric boundary layer, no satellite to satellite bias, self calibration, independent of processing centre, the spherical symmetry assumptions used in Abel inversion is known to be the most significant source of error [7] [8] . TEC calculation involves the integration of electron density along the straight line (signal path) between the LEO receiver and GPS transmitter (as demonstrated in Equation (11)). Measuring the phase delay of radio waves at L1 and L2 frequencies from GPS as they are occulted by the earth atmosphere is illustrated in Figure 1(a) and note that the sample rate per second is 1 Hz. The simulated STEC (Slant TEC) is then inverted into electron density profile along the tangent points by the same Abel inversion software (Bernese v5.0 software package) used for CDAAC electron density profiles retrieval [9] .
Ionospheric Bending and Refractivity
The refractivity of the radio signals passing through the ionosphere depends majorly on the amount of electron density and ions, time of the day and carrier frequency. It is expressed as shown in Equation (1) in which the refractive index is proportional to the electron density and inversely proportional to the square of the carrier frequency:
(1)
where the index of refraction is denoted with n2, N is the electron density (el cm−3) and f is the signal frequency (Hz).
As mentioned earlier; the constraint problem when constructing a 3-D refractivity field from 1-D observational data is the assumption; spherical symmetry of refractivity around the extended region of the perigees of the sounding rays. The Snell’s law is obtained when the 3-D ray equations are integrated under this assumption. Therefore, the bending angle α and excess phase s are given as follows:
(2)
(3)
where a is defined as the impact parameter (), ρ is the distance from the ray tangent point to the centre of sphericity of the refractivity field, represent the fractional radius (r is radius) and LGL denotes the distance between the GPS and LEO satellites. Figure 1(b) shows the bending angle α and impact parameter a that are used in Equations (2) and (3) as well as the ray geometry.
Ionospheric refractivity is negligible at GPS orbit altitudes in the first term of Equation (2) while the second term (LEO) is not. The bending angle does not depend on the positions of the GPS and LEO satellite, excess phase and ray separation from a straight line connecting the satellites. Calculations based on Figure 1(b) were done by numerical integration of Equations (2) and (3) where L1 frequency (f = f1) and L2 (f = f2).
2.3. Theory of Abel Inversion and Techniques
Abel Inversions
In this section under the assumption of spherical symmetry, the formulations for the construction of electron density profile through the use of bending angle data and excess phase (TEC) data was conducted.
1) Abel Inversion Obtained from Bending Angle Data
To determine the bending angle from excess phase, the relationship between Doppler shift of the carrier frequency, and the projections of satellite velocities at the GPS and LEO positions is as follows:
(4)
For and represent 2-D projections of GPS and LEO satellite 3-D velocities respectively on the occultation plane, c is defined as the velocity of light in a vacuum, angles β and f are defined in Figure 1(b) while nLEO and nGPS denotes the indices of refraction. It is observed that Equation (4) has insufficient information to solve for fGPS and fLEO as it contain only Doppler data. Based on the assumption made, Snell’s law is incorporated as shown below to complement Equation (4):
(5)
Solving both Equations (4) and (5) iteratively to obtain the impact parameter, a and the bending angle, α with the aid of Equation (6):
(6)
Solving Equations (4) and (5) causes an error due to the assumption of the refractivity at GPS and LEO positions to be unity.
2) Abel Inversion Obtained from Total Electron Content Data
Total electron content is a measure of propagation delay time of the radio signal transmitted from the satellite to the receiver. The total electron content (TEC), T along the signal path can be related to electron density N, excess phase S and index of refraction n as:
(7)
S is measured in metres (m) and T in (TEC (0.1 TECu).
If bending is neglected, TEC can be obtained from as shown in Equation (8):
(8)
It is advantageous using or to estimate TEC as it automatically eliminate the orbit and clock errors from the difference while the disadvantage is that when anti-spoofing is activated L2 noise is introduced which may degrade the TEC inversion results [10] [11] . Just like the inversions through bending angles using the assumption of spherical symmetry, it is thus adopted in electron density for TEC inversions. Then, TEC in Equation (8) is related to electron density as shown in Equation (9):
(9)
r0 is defined as the impact distance of the straight line connecting the GPS and LEO satellites. Hence, TEC data calibration is done along the ray path as demonstrated in Figure 2 and expressed as in Equation (10):
(10)
The calculation of r0 (the straight line impact distances) is carried out for all observational data with respect to the point where there is maximum impact distance. Then, the interpolation of uncalibrated TEC as a function of impact distance with the aids of cubic splines. Calibration of TEC is the next stage, i.e.. Now, the calibrated TEC Equation (10) becomes:
Figure 2. Total electron content and bending angle calibration geometry diagrammatic illustration.
(11)
Equation (11) can be further expressed in terms of as:
(12)
3. Results and Discussions
3.1. Variability of the Ionospheric Monthly TEC and Electron Density
Monthly TEC and electron density variation analysis are highly recommended and of great importance in the investigation of electron density and TEC since ionospheric irregularities are natural phenomena and contribute to a large extent ionospheric communication outages of radio signals. The only way to minimize its effect is to monitor the electron density and total electron content (TEC) though a very difficult task in real-time as it is constantly changing over time and space. The intensity and effects of the sun activity varies from month to month and even the hours of solar activity are not the same. The study was carried out during the increasing solar activity according to 24th solar cycle prediction by National Aeronautics and Space Administration as shown in Figure 3.
COSMIC radio occultation technique was used to retrieve electron density profile (TEC) data in Thailand from 2010 to 2013 to extract the monthly mean TEC and mean electron density as shown in Figure 4 and Figure 5 respectively.
There were periodical characteristic in the inter-monthly variation of the TEC and electron density in Thailand. It can be observed in 2010 that the mean peak TEC value was 140 (TECu) in April while the months of May, September and October were high above 100 (TECu). 2011 showed a higher mean peak TEC value of 260 (TEC) in the month of December while the Months of April, May, September and October was a bit higher than in 2010. It can be deduced from monthly TEC variation analysis in 2012 and 2013 that, the mean peak TEC value further increased to 275 (TECu) and 375 (TECu) respectively in the month of April and May. Although there was data loss during the months of February, September, October, November and December in 2013 which could be due to system failure (hardware or software problem). A critical look and comparison of Figure 4 and Figure 5 shows electron density is proportional to total electron content.
Figure 3. National aeronautics and space administration 24th solar cycle prediction [12] .
Figure 4. Variability of the monthly total electron content.
Figure 5. Variability of monthly electron density.
Generally, from 2010 to 2013, there was consistence and drastic increment in the values of TEC and electron density. Take for instance, the mean peak TEC and mean Peak electron density values increased from 140 (TECu) to 375 (TECu) and 3.8 × 105 cm3 to 8.2 × 105 cm3 respectively. It means that the monthly mean TEC and electron density characteristics correlate with 24th solar cycle prediction by National Aeronautics and Space Administration in which the solar activity was on the increase within the period of this investigation [13] . Except in 2011, it was also clear that the months that made up summer season were observed to be very high.
3.2. Seasonal Variability of Ionospheric TEC and Electron Density Monitoring
Since seasonal variation plays a vital role in electron density and TEC variation, it is thus pertinent in this research to carry out seasonal ionospheric analysis to investigate the increasing solar activity and TEC behaviour during this period. In this paper, there are three distinct seasons which are winter, summer and rainy in Thailand. Winter comprises November, December, January and February. The summer is made up of March, April, and May while rainy season consists of June, July, August September and October. For each season, the mean TEC and electron density values were determined and characteristics displayed.
TEC and electron density values changed dynamically from one season to the other and the observation made in the course of this research was that the ionospheric mean TEC and electron density was at peak during the summer period and the least ionospheric mean TEC and electron density was experienced during the winter season throughout the period of this investigation as displayed in Figure 6 and Figure 7 respectively. These results are in conformity with theoretical background that the impacts of solar activity on the earth’s atmosphere must be taken into account [13] . In summer solar activities was significantly high while during the winter the energy released by the sun was very low. Therefore, there was high concentration of free electrons interfering with the radio signal in summer than in winter. Other paper that gave similar results includes [14] .
Based on the results obtained, the seasonal mean TEC and electron density increased from 2010 to 2013 except for winter and rainy seasons that declined in 2013 which may be due to data loss, cycle slips or multipath. Following the trend generally, there was increment in the TEC and electron density characteristics and that means it complied with NASA solar cycle prediction curve.
Figure 6. Seasonal variation of ionospheric total electron content.
3.3. Annual Ionospheric TEC and Electron Density Monitoring
Solar activity differs annually as clearly predicted by NASA curve and therefore required to be studied in order to monitor and determine why there are annual variations in electron density and TEC. Annual TEC and electron density variations from year 2010 to 2013 are shown in Figure 8. The trend shows the mean annual TEC and electron density kept steadily increasing. It means that characteristic of electron density and TEC variations were influenced by solar activity monthly, seasonally and annually. The monitoring of electron density and TEC characteristics as seen helps space application industries in their planning and gives awareness of the kind of the problem they are likely to face at a particular time and enable them proffer solution to it. Therefore, the signal degradation and interference was on the increase from 2010 to 2013.
Figure 7. Seasonal variation of ionospheric electron density.
Figure 8. Annual electron density and total electron content variation monitoring.
4. Conclusions
This research work has investigated the dynamic behaviours of ionospheric electron density and TEC within Thailand geographical region. Ionospheric electron density and TEC variations were monitored monthly, through the seasons and on annual basis from 2010 to 2013. This research was conducted to monitor the instability of solar activity and also to guide against serious adverse effects of electron density concentration and TEC from time to time. The results showed an increasing trend of ionospheric electron density and TEC values except in 2013 when the values declined for the rainy and winter seasons. It could be as a result of data loss, cycle slips and possibly multipath within the period. It was obviously established that electron density was proportional to total electron.
Based on the results and observations, it is therefore concluded that telecommunications, geodetic and navigation systems will experience more error during summer (strong solar activity). Thus, users of GPS applications such as security agencies, military and aviation industries in Thailand may take the advantage of this clue to avoid the period of atmospheric disturbance and error.
Acknowledgements
The authors express their sincere thanks and gratitude to University Corporation for Atmospheric Research (UCAR) Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) Program Office for their support and permission to freely use data from their website. Also profound gratitude goes to the management of African Regional Centre for Space Science and Technology Education in English, Ile-Ife, Nigeria for the great opportunity given to be a part of 2014 GNSS Space Application Training in Beihang University, China.