Measurements of delays of the signals radiated by transmitters of navigational satellites allow us to obtain the total electron content (TEC). In addition, measurements of TEC allow solving problems such as development of local, regional, and global models of TEC and correction of ionospheric delay for increasing accuracy of positioning. Now, it is possible to set the task of calculation of critical frequency foF2 with the use of experimental values of TEC in a global scale. For this purpose it is necessary to know an equivalent slab thickness of the ionosphere which is a coefficient of proportionality between TEC and a maximum density of the ionosphere. The present paper is devoted to the analysis of investigation and utilization of this parameter. It is shown that (1) existing models of τ are not empirical and not always can provide an adequate accuracy of foF2 calculation, (2) experimental median (med) provides much larger accuracy of foF2 calculation than the empirical model and variations from day to day and allows filling gaps in the ionosonde data, and (3) it is possible to use a hyperbolic approximation and coefficient for development of a global model of .

1. Introduction

With the advent of navigation satellites, it became possible to measure the total electron content TEC of the ionosphere, which in turn is essential for the functioning of navigation, telecommunications, and other systems. For example, in paper [1], two categories of radio systems were described in terms of ionospheric dependence. The first category includes systems that rely on the critical frequency of the ionosphere foF2: MF and HF communication, HF broadcasting (“short-wave” listening), OTH radar surveillance, HFDF, and HF SIGINT. The operation of these systems requires knowing the behavior of this frequency. Systems of the second category, namely, satellite and navigation (e.g., GPS and GLONASS) communication, space-based radar and imaging, terrestrial radar surveillance, and tracking and others, are influenced by the ionosphere. If in the first case it is necessary to know the state of the ionosphere at the bottom side, namely, up to the maximum height of the layer F2, hmF2, in the second case it is necessary to know the density of electrons in the topside ionosphere (above hmF2) and higher up to the height of high-altitude satellites. A model profile of this part may be corrected by TEC. To solve this problem it is also necessary to know foF2; therefore, one of the most important applications of TEC measurements using navigation satellites is to estimate foF2. Traditionally, values of foF2 are measured by ground-based ionosondes, but their network is sparse. Measurement of TEC has several advantages: a large number of signal receivers, the possibility of continuous, global monitoring, and availability of initial data using network. These advantages allow scientists to study local, regional, and global features of ionospheric behavior independently, quickly, and simultaneously. Definition of morphological and perturbed features has allowed promoting essentially in understanding of spatial and temporary variations of the ionosphere. TEC is widely used for development of local, regional, and global models of TEC, for example, [24], which are necessary for correction of ionospheric delays to increase accuracy of positioning [5, 6]. With use of TEC ionospheric disturbance indexes were developed [7, 8]. In the present paper the basic attention is given to problem of estimation of foF2 using TEC measurements. To calculate foF2 using TEC it is necessary to know a proportionality coefficient between TEC and NmF2. This coefficient is an equivalent slab thickness τ of the ionosphere. This paper is devoted to joint use of the observational values of TEC(obs) and the equivalent slab thickness τ to estimate foF2. The opportunity of development of a global τ model is investigated.

2. Methodology and Used Data

To obtain ionospheric conditions, in particular, values of foF2, empirical models are widely used, for example, IRI and NeQuick. These models define an average condition of the ionosphere; however for functioning of modern systems it is necessary to know instantaneous values and variations of foF2 from day to day. In the present paper, the instantaneous values of TEC measured by means of navigating satellites are used to estimate these values and variations. The methodology of foF2 estimation, used in this paper, consists of joint use of observational values TEC(obs) and an equivalent slab thickness τ of the ionosphere. The following options are considered: (1) τ(IRI) of the IRI model [9], (2) of the Neustrelitz models, (3) the model of Muslim and coauthors [10], and (4) median (med) of observational values τ(obs). τ(IRI) is traditionally used [1113]. Besides, the IRI model is the international standard to which results of other models are usually compared. The model = TEC(NGM)/NmF2(NGM) is actually used in [14] on the basis of Neustrelitz models for TEC and NmF2 [15, 16], but without sufficient validation. This model has been tested in [17]. It was obtained that use of τ(NGM) did not improve results of use of τ(IRI) in most cases. Authors of [18] complained that researchers develop ionospheric models but not a model of τ. Meanwhile, authors of [17] only have shown the values for the 6 stations in different regions of the globe. Moreover, values of τ(obs) were calculated from ionospheric electron content IEC values, that is, in the range of heights from the bottom of the ionosphere to the height of 500 km while TEC is measured for the satellite altitude (about 20,000 km). The latest step has been done in [10], where a model of average values of τ was developed by expansion in Fourier series according to the TEC of 21 stations. For TEC, monthly averages of the global CODE map were taken, and for foF2, the monthly medians. To test the model, data of 13 stations were used in such a way as to get results for various latitude zones (low, mid, and equatorial). The results were obtained for quiet and disturbed conditions by comparison with the results of the IRI model that takes into account STORM factor. The assumptions made in developing the model are as follows: (1) a linear dependence of all parameters (TEC, foF2, τ) on solar activity, (2) absence of longitudinal dependence of these parameters at the same local time LT, (3) transition from a geographic to a geomagnetic coordinates does not affect the dependence of parameter variations on LT, and (4) invariance of τ in quiet and disturbed conditions. Results were obtained for 5 magnetic storms of varying intensity in the period 2000–2014. They are described in detail for several stations and individual disturbances. The general conclusions are as follows: (1) the new model enables improved compliance with measurements compared with the model IRI-STORM in the mid- and low latitudes and (2) the model gives the deterioration in equatorial latitudes in quiet and disturbed conditions. However, as could be seen from Table of paper [10], deterioration takes place in quiet conditions even for the midlatitude Chilton station and for the low-latitude Del’Ebre station. Probably, it is connected with that the assumptions, taken as a principle of model, are not absolutely strict. In particular, it is possible to show that τ(med) depends on solar activity in the nonlinear way. The authors themselves point out that the presence of longitudinal dependence can be the cause of the deterioration. Illustration of violation of the assumption (2) is given in Figure 1 for stations in different zones for March 2015. Figures are given for τ(med) and τ(IRI) in (a) midlatitude zone, (b) low latitudes, and (c) an equatorial region. Latitudes of couples of stations are very close. Couple Juliusruh-Novosibirsk belongs to midlatitudes and couple Nicosia-Kokubunji is located in low latitudes. Couple Sao Luis-Fortaleza is in the equatorial zone. Authors of [10] do not apply their model to high latitudes and auroral zones; however, as in the papers [19, 20] the possibility of using the IRI model in these areas was shown; results are given also for couple Tromsö-Amderma. The parameter, for which the plot is under construction, is moved in heading together with units of measurements. On an axis , local time LT is postponed.

It is seen as a good agreement between the values of τ(IRI) of two stations; however, large differences are between τ(med). Violation of condition (4) is illustrated in Figure 2 where values of τ(med) and τ(IRI), close to values in quiet conditions, and values τ(obs) in disturbed days are compared. Results are given for the etalon Juliusruh station and two moderate disturbances in July 2004 with a minimum of Dst = −197 nT and in December 2006 with a minimum of Dst = −147 nT.

Unlike the model [10] it is offered to investigate two possibilities: (1) hyperbolic approximation of τ(obs) and (2) use of a median τ(med) in comparison with τ(IRI). Hyperbolic approximation [21] was chosen on the basis of visual similarity and simplicity of analytical representation. The illustration is given in Figure 3. In case of the ESA map dispersion is stronger.

Comparison is carried out by means of several types of deviations of the calculated values of foF2 from observational foF2(obs). The first deviation concerns values of foF2(IRI) for the initial IRI model IRI. The second deviation concerns value of the critical frequency calculated with use of TEC(obs) and τ(IRI). Values of τ(IRI) are obtained from a relation τ(IRI) = TEC(IRI)/NmF2(IRI). The use of experimental values TEC(obs) provides calculated values NmF2(calc) = TEC(obs)/τ(IRI) and foF2(τIRI) = (NmF2(calc)). The third deviation concerns value of the critical frequency calculated with use of TEC(obs) and τ(med). The following expressions are used: τ(med) = med(TEC(obs)/NmF2(obs)), NmF2(calc) = TEC(obs)/τ(med), and foF2(τmed) = (NmF2(calc)). Corresponding deviations are calculated as follows. The magnitude = represents the first deviation. The deviation = is the second deviation, that is, a difference between foF2(obs) and values calculated with use of τ(IRI) and TEC(obs). The third deviation is written as = and is defined as a difference between foF2(obs) and values calculated with joint use of TEC(obs) and τ(med). Monthly averages were calculated. To compare these deviations an efficiency coefficients of use of instantaneous TEC(obs) for calculation of instantaneous values of foF2 are introduced. Coefficient τIRI = is an efficiency coefficient of joint use of TEC(obs) and τ(IRI). Coefficient eff = is an efficiency coefficient of joint use of TEC(obs) and τ(med). Thus, coefficients τIRI and eff characterize how many times the compliance between the calculated and experimental values of foF2 using τ(IRI) and τ(med) as compared to the initial model is increased.

In this paper sources of observational data are IONEX files from ftp://cddis.gsfc.nasa.gov/pub/gps/products/ionex/ for TEC available with step 2 hour and SPIDR database http://spidr.ngdc.noaa.gov/spidr/index.jsp for foF2. The JPL map is used in most cases.

3. Distinctions between Critical Frequencies foF2 Estimated with Use of TEC(obs) and Various Options of τ

It is obvious that differences of critical frequencies foF2, estimated with use of TEC(obs), are defined by distinction of values τ. Values of τ(med) for the JPL map and τ(IRI) are presented in Figure 4 for several stations with long arrays of observations in various regions of the globe: Juliusruh (mid latitudes), Thule (a polar zone), and Ascension Island (an equatorial region) as the monthly averages for December and July of the corresponding year. Dots of various colours specify years. On an axis , time UT is postponed.

The essential difference between τ(IRI) and τ(med) is obvious in some cases. This essential difference leads to corresponding difference between derived values of foF2. If values of τ(IRI) and τ(med) are close, deviations of foF2 will be small. Such cases are characteristic for the Thule station. Deviations of calculated foF2 from experimental values (Figure 5(a)) and efficiency coefficients eff (Figure 5(b)) are shown in Figure 5 for stations Juliusruh, Thule, and Ascension Island. On left graphs black dots concern values of the initial IRI model (option of IRI2001 to distinguish from other options corr and NeQuick of the IRI2012 model). On Figure 5(b) this case corresponds to the curve . Dark blue circles represent results for a case of use TEC(obs) and τ(IRI). Red dots correspond to a case of use of TEC(obs) and τ(med). If the ratio is equal to 1, the use of the experimental value of TEC and the equivalent slab thickness provides the same results as the model itself without using TEC(obs). If the ratio is higher than 1 then the use of TEC gives better results than the model. If the ratio is less than 1 the use of TEC deteriorates results compared with the model. Coefficients eff concern with a case of use of TEC(obs) and τ(med). Coefficients τIRI describe results for a case of use of TEC(obs) and τ(IRI).

In all cases the least deviations are provided at the use of TEC(obs) and τ(med), and coefficients eff always are higher than 1. It testifies that the use of TEC(obs) provides better results than the initial model. Situation with use of TEC(obs) together with τ(IRI) is ambiguous. In most cases use of TEC(obs) and τ(IRI) improves results of the initial model (curves shown by dark blue circles in comparison with black curves) but at every station there are periods when the coefficients τIRI are less than 1 and the use of TEC(obs) deteriorates results of the initial model. These results lead to the following conclusions: (1) in the equatorial zone deviations are larger than in the mid- and high-latitude regions, (2) use of TEC(obs) together with τ(IRI) does not always improve correspondence between the calculated and experimental values of foF2 in comparison with the initial IRI model, (3) use of τ(med) leads to more accurate values of foF2, and (4) the coefficient eff is higher than 1.

4. About the Possibility of Developing a Global Model of τ(med)

To develop a model of τ(med) on a global scale there are possible several approaches: (1) the construction of surface function such as kriging from values of τ(med) at several points, (2) two-parameter model based on hyperbolic approximation [21] and the use of the coefficient , and (3) the model NGM. Regarding the first option, in paper [22], some doubts were expressed. Since the construction of the model using values of τ(obs) is not possible due to the high variability of values (in particular, presunrise peak at some latitudes), we make an attempt to use a hyperbolic approximation and coefficient . Results are illustrated for the case of March 2015 for the largest number of stations. For the hyperbolic dependence, coefficients and of were modeled. The results are shown for the two regions (2 and 4) and two wider zones (Lat1 and Lat2). Zone 2 (15°E < λ < 40°E) includes 8 stations in the European and Russian regions. Zone 4 (110°E < λ < 170°E) includes 9 east stations. Area Lat1 includes 13 stations, mostly over American continent of northern and southern hemispheres. Area Lat2 includes 20 stations from European, Siberian, and south-eastern regions. Behavior of the coefficients and for these regions is shown in Figure 6. The axis shows a latitude of stations.

Average values were calculated. They are 250.62 km and 4757.36 m−2 for zone 2, 280.21 km and 4386.01 m−2 for zone 4, 282.92 km and 5581.81 m−2 for Lat1 area, and 257.63 km and 4276.64 m−2 for Lat2 zone. Results are shown in Table 1. The table contains the following data. Column 1 indicates the station name and its location (corresponding region). The second column shows the coefficients of the hyperbolic dependence of τ(obs) according to the respective station. The third column specifies the conditions to which the two rows of values belong. The top line indicates average for all days of the month, at the bottom, the average for disturbed days (from 16 to 21 March). Symbol means full month and concerns magnitudes averaged on disturbed days. The fourth column shows the results for the initial IRI model and the fifth one the absolute difference between the experimental values foF2(obs) and values calculated using TEC(obs) and τ(med).

Column 6 contains deviations of frequencies calculated using the coefficients and of the hyperbolic approximation for a given station. The rest columns give results using the coefficients of the areas referred to in the column header. All of these values should be compared with values of the IRI model in bold.

It can be seen that all values are larger in disturbed days and the largest differences correspond to the initial IRI model. The table illustrates opportunity to use the coefficients of one region for the calculation of foF2 of other zones. This demonstrates a global character of the τ(med) model.

One important issue is the dependence of the coefficients on solar activity. Figure 7 shows the coefficients and for different years, arranged in descending order of solar activity.

The use of coefficients can be another method of constructing the global τ(med) model. Definite advantage of this model might be in the fact what in its denominator is the value of τ(IRI), which has a global character, and a small change of in areas with close latitudes. Distinctive feature is a development of a model for each hour. The degree of proximity of values is best illustrated in the circular charts. An example of some charts is shown in Figure 8 for region 4, for UT = 0, 6, 12, and 18 in March 2015. The red line shows the value of the coefficient ; green points are the average values and the blue marks mean a circle with radius . Besides dots values of latitudes of used stations are given.

The model is the average value (mean). The algorithm of its use is reduced to what the new value is computed and this new value is used together with the experimental value of the TEC to calculate foF2. To test the algorithm, (mean) averages were calculated for 7 stations and used to calculate foF2 for the 8th station. Results in the form of deviations of calculated foF2 values from foF2(obs) are presented for the four experimental stations: Pruhonice, Gorkovskaya, Tunguska, and Ramey in Figure 9.

Black dots show differences between values of foF2(IRI) of the IRI model and observational values foF2(obs) (symbol IRI). Blue dots indicate deviations in the case of joint use of TEC(obs) and τ(IRI) (symbol τ(IRI)). Green triangles represent the results for the case of joint use of TEC(obs) and τ(med) (symbol τ(med)). Red diamond shows deviations for joint use of TEC(obs) and a new coefficient (symbol new). Zones, which included stations Pruhonice, Gorkovskaya, and Ramey, consist of 8 stations. The coefficient was calculated by the values of of seven stations. For the Tunguska station, area includes 6 stations, and the coefficient was calculated according to 5 stations.

Except disturbed days (17–20 March) average foF2(IRI) corresponds to instantaneous values well. The use of TEC(obs) with τ(IRI) improves the results of the initial model (previous curve) except case of the Ramey Station. The use of the TEC(obs) with a median τ(med) gives the better coincidence with observational values foF2(obs). Naturally, the value of for the model is slightly larger than for τ(med), but less than for the other options. This testifies to a possibility of development of the global model of τ(med).

5. Additional Comments Concerning Use of TEC(obs)

In Introduction, a series of advantages of the use of TEC for calculation of foF2 was marked; however there is also a series of shortages. It is possible to name among them the following: dependence of TEC on type of receivers [23], a large number of methods, and TEC values which can deviate up to 2 and more times [24]. In papers [17, 20, 25] of authors of this article, the large attention was given to influence of differences of TEC on values of foF2. On an example of global maps JPL, CODE, UPC, and ESA it was shown that, despite the large differences of TEC, deviations |ΔfoF2| were less than deviations for the initial IRI model practically for all maps. It testifies to good calibration properties of the equivalent slab thickness τ, reducing influence of differences of TEC. In most cases, the least values of |ΔfoF2| are provided by the JPL map. Very close results were obtained for the CODE map. The reasonable compromise is the utilization of values IGS which are robust averages from values of 4 maps [26]. For them values of |ΔfoF2| lie between values for JPL and CODE maps. The use of τ(med) can play an important role for an estimation of foF2 during disturbances and for filling of gaps of the data. Examples for these cases are given in Figures 1012 for disturbances of April 2001 according to stations Juliusruh, Thule, and Loparsk. Results are presented in the form of daily dependences of TEC and foF2 for several days. Experimental values of TEC are given together with medians; results for foF2 are given for 4 options: (1) experimental instantaneous values (symbol obs), (2) values of the initial IRI model (symbol IRI), (3) values calculated with use of TEC(obs) and τ(med) and designated by symbol rec, and (4) monthly medians of foF2 (symbol med).

Figure 10 shows two disturbances: positive and negative when values of foF2 essentially differ from medians. The proximity of values for the IRI model to medians testifies to high quality of the model which is median one and cannot reflect instantaneous values. Values of the third option are closest to foF2(obs) and reflect disturbed conditions. Figure 11 concerns the auroral Thule station.

Perturbations also have positive and negative character. Both have a global character. In this case data of foF2(obs) were available only in the first period and here values for the third option are closest to foF2(obs). However in this case values of the initial model essentially differ from medians. In the second period, observational data of foF2 were not available, but strong negative perturbation in behavior of TEC is seen. It is possible to assume that the critical frequency could experience strong negative perturbation during these periods according to the values reconstructed with use of TEC(obs) and τ(med). For confirmation of such character of behavior, it is possible to use data of other high-latitude stations. Data of the Tromso was available for three days but, in these days, some values were absent; for the Sodankula station all the days had the incomplete information.

The Loparsk station was chosen having the fullest data though often the F2 layer was shielded by sporadic layers Es. Results of calculations are presented in Figure 12. One can see the positive and negative disturbances. In the first period, the IRI model and case of joint use of TEC(obs) and τ(med) showed values exceeding the median. Nevertheless, there is a better match with foF2(obs) than the median. In the second period only TEC(obs) together with τ(med) has provided negative deviations from the median and a better match with foF2(obs). Average coefficients eff and τIRI were equal to 1.2 and 1.1 in the first period and 3.5 and 2.9 in the second period.

Thus, the use of τ(med) and TEC(obs) allows defining character of perturbation and to calculate values of foF2 more accurately.

6. Conclusions

With the advent of the navigation satellites providing measurement of the total electron content of the ionosphere, this parameter has got the large practical value. In the presented paper, opportunity of use of observational values TEC(obs) to estimate disturbed conditions is investigated, and, as an example, stations with long-term arrays of observation are chosen in all regions of globe from high latitudes to equatorial area. The key role belongs to an equivalent slab thickness of the ionosphere. Results are presented as comparison of deviations of the calculated values of the critical frequency foF2 from experimental magnitudes foF2(obs) for 3 options: (1) initial model, (2) joint use of TEC(obs) and τ(IRI), and (3) joint use of TEC(obs) and τ(med). The main conclusions are as follows.(1)It is shown that the model IRI provides good conformity of calculated foF2 with experimental medians; however values of τ(IRI) strongly differ from an experimental median τ(med) both on magnitude and on character of daily variations.(2)Joint use of TEC(obs) and τ(IRI) not always improves results of the initial model.(3)Joint use of TEC(obs) and τ(med) always improves results of the initial model. Thus, because of calibration properties of τ(med) (good match to TEC(obs)) results for all global maps, TEC values which can differ 1.5–2 times are better than for the initial model. The map JPL provides the best conformity with foF2(obs) in most cases. Close results can correspond to the map CODE.(4)Experimental values of foF2 are often absent during disturbances. In this case, TEC is a unique source of the information on character of disturbances, and joint use of TEC(obs) and τ(med) allows estimating foF2 and filling gaps of vertical sounding data.(5)The positive fact is attempt of construction of empirical model of τ (NGM, Muslim et al. [10]) which has shown that for any model there are regions and conditions of solar activity at which the model provides the better results than the initial model. Two approaches, offered in the given paper, specify opportunity of construction of model τ, improving results of the initial model on a global scale.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.


The authors thank the scientists who provided data of SPIDR, global maps of TEC, and operation and modification of the IRI model. The work of Olga Maltseva was supported by Southern Federal University Grant no. 213.01-11/2014-22.