A technique for routinely updating the ITU-R database using radio occultation electron density profiles

Well credited and widely used ionospheric models, such as the International Reference Ionosphere or NeQuick, describe the variation of the electron density with height by means of a piecewise profile tied to the F2-peak parameters: the electron density,NmF2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_m \mathrm{F2}$$\end{document}, and the height, hmF2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h_m \mathrm{F2}$$\end{document}. Accurate values of these parameters are crucial for retrieving reliable electron density estimations from those models. When direct measurements of these parameters are not available, the models compute the parameters using the so-called ITU-R database, which was established in the early 1960s. This paper presents a technique aimed at routinely updating the ITU-R database using radio occultation electron density profiles derived from GPS measurements gathered from low Earth orbit satellites. Before being used, these radio occultation profiles are validated by fitting to them an electron density model. A re-weighted Least Squares algorithm is used for down-weighting unreliable measurements (occasionally, entire profiles) and to retrieve NmF2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_m \mathrm{F2}$$\end{document} and hmF2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h_m \mathrm{F2}$$\end{document} values—together with their error estimates—from the profiles. These values are used to monthly update the database, which consists of two sets of ITU-R-like coefficients that could easily be implemented in the IRI or NeQuick models. The technique was tested with radio occultation electron density profiles that are delivered to the community by the COSMIC/FORMOSAT-3 mission team. Tests were performed for solstices and equinoxes seasons in high and low-solar activity conditions. The global mean error of the resulting maps—estimated by the Least Squares technique—is between 0.5×1010\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.5\times 10^{10}$$\end{document} and 3.6×1010\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3.6\times 10^{10}$$\end{document} elec/m−3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-3}$$\end{document} for the F2-peak electron density (which is equivalent to 7 % of the value of the estimated parameter) and from 2.0 to 5.6 km for the height (∼\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}2 %).


Introduction
Well credited and widely used ionospheric models, such as the International Reference Ionosphere (IRI) (Bilitza 2001) or the NeQuick (Nava et al. 2008), describe the variation of the electron density with height by means of a piecewise profile tied to the F2-peak parameters: the electron density N m F2, and the height h m F2.Accurate values of these parameters are crucial for retrieving reliable electron density estimations from those models (Bilitza et al. 2012).The best option is to measure these parameters with an ionosonde or a incoherent scatter radar but, when this possibility is not available, IRI and NeQuick rely on the so-called ITU-R (Radio Communication Sector of the International Telecommunications Union) maps (ITU-R 1997).These maps are based on a pioneering work published in the early 1960s by William Jones and Roger Gallet (Jones andGallet 1962, 1965;Jones and Obitts 1970), who developed an ad-hoc mapping technique (particularly well adapted to the computers of their time) to compute monthly median values of the critical frequency f 0 F2, and propagation factor M 3000 F2, of the F2-peak (see Sect. 2.1).IRI or NeQuick models compute the values for the N m F2 and h m F2 driving parameters (see Sect. 2.2) using the ITU-R maps.
The ITU-R technique for computing the f 0 F2 and M 3000 F2 parameters depends on two sets of numerical coefficients (hereafter called the ITU-R database), one for low and the other for high solar activity.These coefficients were adjusted to data which were registered in the late 1950s on a network of ionosondes whose geographic distribution was rather inhomogeneous.For a solar activity level in between (outside) those defined for the database, the coefficients must be linearly interpolated (extrapolated).The ITU-R technique also requires the use of the 'modip' latitude (Rawer 1984), which depends on the geomagnetic inclination.Then to be consistent with the data used to build the maps, the geomagnetic field model of epoch 1960 should be used for computing the modip.
The present contribution is a work-in-progress report on a technique that we have developed to routinely update the ITU-R database using radio occultation electron density profiles (RO profiles) from the COSMIC/Formosat 3 mission (see Sect. 4.1).This development was guided by the interest in establishing a robust technique that could be automatically implemented in a computer and run autonomously and routinely to keep the database permanently updated.The technique takes advantage from the following facts: (i) the data used for updating the database are representative of the actual ionospheric conditions and related to the present geomagnetic field); (ii) at monthly intervals the data coverage is quite homogeneous both in space and in time; (iii) it avoids the linear interpolation/extrapolation between low and high solar activity periods, and (iv) at mid-latitudes the N m F2 and h m F2 values retrieved from RO profiles have proven to be in good agreement with ionosonde observations (Jakowski et al. 2004;Angling 2008).The main disadvantage is that RO profiles become less accurate at low latitudes.This effect is particularly evident in the equatorial anomaly region, where the electron density distribution presents strong spatial gradients (Schreiner et al. 2007;Yue et al. 201).
To mitigate the impact of this loss of accuracy, we have developed a special processing strategy consisting on fitting every single RO profile to our electron density model (i.e.: the La Plata Ionospheric Model, LPIM; see Sect.2.2 and Brunini et al. 2011Brunini et al. , 2013b)), and carefully checking the statistical significance of the residuals obtained with the adjustment.This procedure eliminated a great amount of single electron density data marked as outliers and/or entire profiles that were considered as unreliable.
In Brunini et al. (2013b), we explored the possibility of using a spherical harmonic expansion with time-dependent coefficients-instead of the Jones and Gallet formalism-to directly map h m F2-instead of M 3000 F2.This way of modeling proved to accurately reproduce the h m F2 values derived from the ITU-R maps, with the benefits of: (i) being better suited than the Jones and Gallet formalism for implementing a data assimilation scheme into the IRI or the NeQuick models; and (ii) avoid the intermediate step of transforming M 3000 F2 into h m F2.Based on that experience, for this work we decided to develop the technique to directly derive maps of N m F2 and h m F2 (instead of f 0 F2 and M 3000 F2) from the RO profiles.
The paper is organized as follows: Sect. 2 presents the background of this research, i.e.: the Jones and Gallet mapping technique, the relations that link M 3000 F2 and f 0 F2 to h m F2 and N m F2, and the ionospheric model used to assess the quality of the RO profiles.Section 3 describes the technique developed to obtain the updated databases for computing global maps of h m F2 and N m F2.This is done in a two-step procedure: (i) h m F2 and N m F2 values are retrieved from every RO profile and statically validated; (ii) the retrieved values are used to compute an ITU-R like database.Section 4 presents the dataset used in this research, the results obtained from these data after applying the technique previously developed and a short discussion about the differences between the h m F2 and N m F2 maps obtained with our technique and with the ITU-R database.The last section summarizes the research carried out and the conclusions extracted from the analysis.
Before closing this introduction, we would like to credit the efforts that have been made over the years to develop predictive models of the F2-peak parameters at global, regional or single site scale.These models are based on measurements provided by ground-based ionosondes in most cases, by incoherent scatter radars and space-based ionosondes in a fewer cases, and by GPS receivers onboard of low-orbiter satellites in the most recent studies.Among others, it is worth mentioning: Bilitza et al. (1979), Fox and McNamara (1988), Rush et al. (1989), Altinay et al. (1997), Wintoft and Cander (1999), Oyeyemi et al. (2005), Gulyaeva et al. (2008), Liu et al. (2008), Zhang et al. (2009), and Hoque andJakowski (2011, 2012).
We also wish to point out that the development of the our technique has been guided by the intention of producing a database that can be used by the IRI and NeQuick models, without having to introduce significant changes in their source codes.

The Jones and Gallet mapping technique
The mapping technique developed by Jones and Gallet (1962) provides values for f 0 F2 and M 3000 F2 anywhere in the world at any time and for any given level of solar activity.Jones and Gallet preferred to map M 3000 F2 instead of h m F2 because its applicability in radio-propagation studies.M 3000 F2 is defined by Eq. (1) where f 0 F2 is the critical frequency of the F2 layer and maximum usable frequency (MUF) is the maximum frequency at which an electromagnetic signal emitted from the ground can reach another point on the ground 3000 km away from the transmitter.
M 3000 F2 = MUF/ f 0 F2. (1) The sets { f 0 F2, M 3000 F2} and {N m F2, h m F2} are related through a system of equations that will be presented in the Sect.2.2.Jones and Gallet (1962) represented the daily variations of the parameters with Fourier series expansions as: (2) where is the parameter ( f 0 F2 or M 3000 F2) to be mapped, t is the Universal Time and J is the maximum number of harmonics for mapping the diurnal variation: J = 6 for f 0 F2 and J = 4 for M 3000 F2.The spatial variation of these parameters is taken into account through the dependence of the Fourier coefficients on the geographic coordinates: where K = 75 for f 0 F2 and K = 49 for M 3000 F2, U 1,0 , . . ., U 2J,K are the numerical coefficients of the expansion (988 for f 0 F2 and 441 for M 3000 F2), and G 0 , . . ., G K are special functions whose explicit form depends on the k index, for example where φ and λ are the geographic latitude and longitude and μ is the so-called modip latitude (Rawer 1984), defined by: where I is the magnetic inclination at an altitude of 350 km above the Earth's surface.
The numerical coefficients U 1,0 , . . ., U 2J,K of Eq. (3) were estimated by means of a Least Squares adjustment to a set of observations collected from 1954 to 1958 at a network of around 150 ionosondes with an uneven global distribution.This procedure led to establish two sets of coefficients for each parameter to be mapped for every month of the year (one set for low and the other for high solar activity).The solar activity level enters into the model with the so-called R 12 parameter which is the 12-month running mean value of the monthly mean sunspot number.For a given month and a R 12 value between R 12 = 0 (low) and R 12 = 100 (high solar activity), the coefficients are linearly interpolated from the tabulated values.For R 12 > 100 the coefficients are extrapolated.

The relation between (
N m F2 and f 0 F2 parameters are related through a simple relation (Bilitza 2002): where The relation between h m F2 and M 3000 F2 has been an active subject of study.Shimazaki (1955) proposed an empirical relationship (based on the strong anti-correlation observed between the mentioned parameters) with the form: where h m F2 is in km.This equation evolved to a more complex one after the works of Wright and Mcduffie (1960), Bradley andDudeney (1973), Eyfrig (1973), Dudeney (1975), Bilitza et al. (1979), Obrou et al. (2003), and Rawer and Eyfrig (2004).In this paper we used the equation established by Bilitza et al. (1979): with the correction factor given by: CF R being the sunspot number, ϑ the geomagnetic latitude and f 0 E the critical frequency of the ionospheric E layer.

The LPIM electron density vertical profile
The model used to represent the electron density vertical profile used in this work is the one implemented in the La Plata Ionospheric Model (LPIM) (Brunini et al. 2013a).It uses three α-Chapman layers (corresponding to the three characteristic ionospheric layers) to compute the electron density N as a function of the height h: where the sub-index i = 1, 2, 3, identifies the layers E, F1 and the bottom-side of the F2, and N mi , h mi and H i are the maximum electron density, the height of that maximum, and the scale height of the corresponding layer (i.e., N m3 ≡ N m F2, h m3 ≡ h m F2, and The electron density in the top-side is represented with a vary-Chapman function (Reinisch and Huang 2001): with ζ going from the lower to the upper integration limit.H 4 (h) is the height-varying scale height given by: where h T is the transition height where the dominant ion specie changes from O + to H + , H T is the scale height at the transition height and p is the shape parameter of the topside profile.All of them are evaluated according to Meza et al. (2008).In summary, the electron density computed with the LPIM profile at any given height can be expressed as: with N i (h) given by Eq.( 10).
The shape of the LPIM profile is determined by the F2 peak parameters (N m3 ≡ N m F2, h m3 ≡ h m F2, H 3 ≡ H F2) because once they are fixed, LPIM computes all the other parameters according to the ITU-R recommendations: where χ is the solar zenith angle, A 1 , A 2 , n 1 and n 2 are functions that depend on the geographic latitude and the solar activity level.The explicit forms of them are contained in the ITU-R (1997) recommendations.For a given location (latitude, longitude and height) and time, the LPIM electron density prediction is given by a function that depends on the three parameters, N m F2, h m F2 and H F2; symbolically: In Brunini et al. (2013a) we described the LPIM model together with a data assimilation technique that is the basis for the data analysis procedure described in the next section.In that paper we also assessed the model performances by comparing its predictions to a dataset independent of the one assimilated into the model.After assimilating groundbased GNSS measurements and RO profiles we compared the LPIM-derived vertical total electron content (vTEC) estimation to Jason1 (J1) mission determinations.The comparison was performed for a complete year of low solar activity and showed a quite good agreement between model and measurements.99 percent of the differences-LPIM minus J1 vTEC-did not exceed the range of [−0.5, −3.4] total electron content units (TECu).The mean difference showed a systematic bias of −1.The process of estimating the N m F2 and h m F2 parameters from the RO profiles consisted on adjusting (by Least Squares) the LPIM profile to every RO profile available for the time period of interest.For doing this, we used the linearized form of Eq. ( 16) to obtain the following equation of observation: where N obs is the electron density value from the RO profile; N pri is the corresponding electron density value computed with a priori values N m F2 pri , h m F2 pri and are the correction to the a priori that would be estimated with the Least Squares adjustment; and v is the deviation between the measured and a posteriori computed electron density value.
The a priori values for N m F2 and h m F2 (i.e.: N m F2 and h m F2 pri ) were computed using the Jones and Gallet (1962) mapping technique to obtain f 0 F2 pri and M 3000 F2 pri , and then, the Eqs.( 6), ( 8) and ( 9) to convert f 0 F2 pri and M 3000 F2 pri into N m F2 pri and h m F2 pri .The a priori value for H F2 (i.e.: H F2 pri ) was computed from f 0 F2 pri and M 3000 F2 pri , using the empirical relation adopted by LPIM (Brunini et al. 2013a), which are in accordance with the ITU-R (1997) recommendations: where f 0 F2 is measured in MHz and H F2in km.
The corrections N m F2, h m F2 and H F2 were estimated along with their standard deviations by means of a re-weighted Least Squares procedure.All the weights were initially set to the same value and then were iteratively modified according to the difference between the RO electron density determination and value predicted by the model with previous solution.The weighting criterion was the 'bisquare' function described by Huber (1981).The iteration was stopped when the changes in the solution became negligible.

Building the database
To overcome the computational limitations of their time, Jones and Gallet (1962) constructed the special basis of mathematical functions like the one presented in Eq. ( 4).This basis fulfilled the condition of orthogonality with respect to the geographic configuration of the data sample, which avoided the inversion of the normal matrix.Based on the encouraging results obtained in (Brunini et al. 2013b), we changed the Jones and Gallet's formulation by a spherical harmonic expansion series with time-dependent coefficients, and directly mapped the N m F2 and h m F2 instead of f 0 F2 and M 3000 F2 parameters (primary parameters of the IRI and NeQuick models).This way of modeling implies that Eqs.
(3) shall be replaced by: where P lm are the associated Legendre's polynomials, μ is the modip latitude (see Eq. 5) and L is the maximum degree of the expansion.For mapping the daily variation of N m F2 and h m F2, we kept the Fourier series expansions with maximum order J given by Eq. ( 2).
The ITU-R method to compute monthly median values of f 0 F2 uses two sets of coefficients (U 1,0 , . . ., U 2J,K ), one for low and the other for high solar activity.The proposed approach requires one set of (J + 1) • (L + 1) 2 coefficients (A c 0,0,0 , . . ., A s J,L ,L , B c 0,0,0 , . . ., B s J,L ,L ) (Eqs. ( 2), ( 19)) for computing monthly mean values of N m F2 for the solar activity level corresponding to the adjusted period.The same applies to the computation of h m F2.
Summarizing the procedure, the N m F2 and h m F2 (and their standard deviations) estimated in the first step (see Sect. 3.1) are adjusted on a monthly basis by a weighted least squares technique to the functions determined by Eqs.
(2) and ( 19), then obtaining a new set of coefficients the map the F2 peak parameters.

The used dataset
The database used in this research is composed of electron density determinations provided by the Constellation Observing System for Meteorology/Formosa Satellite 3 mission (COSMIC/FORMOSAT-3) (Anthes 2008).This mission is a collaborative project between the National Space Organization (NSPO) of Taiwan and the University Corporation for Atmospheric Research (UCAR) of the USA and consists in a constellation of six microsatellites that were launched on April 2006 at an altitude of 512 km of altitude with an inclination of 72 • (Cheng et al. 2006).During the first 17 months after the launch the satellites were gradually moved into their final orbits at ∼800 km, with a separation of 30 • in longitude between the nodes of neighboring orbital planes.The primary payloads of the satellites are GPS radio-occultation receivers, which allow computing profiles of atmospheric refractivity (Kuo et al. 2004).These profiles can be used later to estimate temperature profiles in the stratosphere, temperature and water vapor profiles in the troposphere and electron density profiles (RO profiles) in the ionosphere (Yunck et al. 2000).
The database used in this work was downloaded from the Data Analysis and Archival Center (CDAAC) site (http:// cdaac-www.cosmic.ucar.edu/cdaac/products.html), at the UCAR.It includes RO profiles from one continuous month for each solstice and equinox for 2007 (low solar activity) and 2011 (moderate to high solar activity).Between 24×10 3 and 64 × 10 3 profiles, comprising 9.4 × 10 6 to 20.8 × 10 6 single electron density values were effectively processed for each season.This database provides a very good coverage either in local or in universal time, latitude, longitude and heightfrom the lower limit of the ionosphere up to the altitude of the satellites (∼800 km).

Results
Table 1 contains the number of profiles, N_PR, and single electron density determinations, N_ED, available (AVAIL column) in the CDAAC site and the number that remained after the quality filter step and were processed (PROC) to obtain the final results.Approximately 20 % of the measurements (including 35 × 10 3 complete profiles) received negligible weights and for this were eliminated.We have  devoted considerable effort to understand the reasons why these data were discarded.About half of the cases coincided with data, and even complete profiles, that indisputably appeared unrealistic.For the second half, it was difficult to ascertain whether the misalignment between data and models should be attributed to data, the model or both.Nevertheless we have verified that the discarded determinations did not exert significant influence on the final results.
The first step of our analysis consisted in the evaluation of the quality of the adjustment described in Sect.3.1, i.e. the agreement between the LPIM fitted profiles and the RO electron density determinations.Figure 1a shows a map with RMS of the adjustments for the 18-20 UT interval, the September equinox of 2007.The interval was selected because during these hours the Ionospheric Equatorial Anomaly passes over the South Atlantic geomagnetic anomaly and represents a challenging scenario for modeling.The 2007 September equinox is the period with the lowest solar activity level analyzed in this work.
Figure 1b shows the corresponding results for the December of 2011 (the period of highest solar activity analyzed in this work).On average, the RMS of the differences presented in these maps are at about 3 % of the corresponding N m F2 values.Similar statistics were obtained for the other equinoxes and solstices of 2007 and 2011.
The next step in our analysis was the assessment of the propagated errors of the estimated N m F2 parameters.Figure 2a shows the N m F2 estimations obtained with our technique for the 18-20 UT interval least squares adjusting data from 31 days around the September equinox, 2007.Figure 2b shows the corresponding formal errors (standard deviation) obtained in the adjustment.Figure 2c, d shows equivalent maps for the December solstice, 2011.Figure 3 shows the corresponding figures for h m F2.The mean errors computed from these maps are lower than 1 % of the estimated values of the parameters.Similar values were obtained for the other UT intervals and the other equinoxes and solstices of 2007 and 2011.Before proceeding to compute the coefficients of the spherical harmonics expansion given by the Eq. ( 19) we needed to fix the values for J and L (the maximum degrees of the expansion) in order to get the best possible results from the available dataset.For this reason, we applied an analysis of the variance (ANOVA) technique (Weisberg 1980) to evaluate the statistical significance of the terms added to the expansion series, when J and L were increased.For the expansion of N m F2 the ANOVA led us to set J = 12 and L = 9, implying the estimation of (J + 1)•(L + 1) 2 = 1200 coefficients (around 20 % more coefficients than the 988 determined by Jones and Gallet (1962)).For the expansion of h m F2, the ANOVA tests gave J = 6 and L = 8, implying the estimation of 567 coefficients (around 30 % more than the 441 determined by Jones and Gallet 1962).These maximum degrees permits a temporal/spatial resolution of approximately 1 h /20 • for N m F2 and 3 h /22 • for h m F2.
Considering the different periods analyzed in this work, the unit-weights standard deviation (Fan 2010), resulting from the Lest Squares adjustments, ranged from ±0.5 × 10 10 to ±1.2 × 10 10 elec/m −3 for N m F2, and from ±1.7 to ±3.4 km for h m F2.The χ 2 statistical test applied to the obtained unit-weight standard deviations did not reject the working hypothesis within a level of significance of 97.5 %.

Discussion
Figure 4a shows a global map of N m F2 monthly median values computed with the Jones and Gallet (1962) technique for 19 UT, September equinox of 2007. Figure 4b shows the corresponding map computed with our technique and Fig. 4c shows the map for the estimated errors.Figure 4d-f shows the equivalent maps for December solstice of 2011. Figure 5 shows the equivalent results for h m F2.
The analysis of Figs. 4 and 5 (and other similar figures representative of the other time intervals, seasons and years considered in this work) shows that the largest errors in the estimation of both N m F2 and h m F2 are found in the region of the Equatorial Anomaly.This effect is a direct consequence of the difficulty to model the electron density distribution in that region and, up to some extent, to the limitations of the RO technique in regions of steep horizontal gradients.The geographic behavior of the errors is also modulated by the data coverage.For example, a clear amplification of the errors over the modip equator in the African continent is produced by the scarcity of RO profiles in that region.
The global mean errors ranged from 0.5 × 10 10 to 3.6 × 10 10 elec/m −3 for N m F2 (approximately 7 % of the value of the estimated parameter) and from 2.0 to 5.6 km for h m F2 (∼2 %).
For N m F2, the largest differences between the Jones and Gallet (1962) and our maps are observed over the Central/South American, Atlantic and African sectors, in a band approximately within ±45 • of modip latitude, when the Equatorial Anomaly is developed over those sectors.The differences range from approximately ±80 × 10 10 elec/m −3 for high solar activity to approximately +20 × 10 10 /−50 × 10 10 elec/m −3 for low solar activity.For h m F2, the largest differences are around ±80 km for both, high and low solar activity conditions.These largest differences appear in different parts of the globe, depending on the time interval, season and year, but a positive difference (our h m F2 values are larger than the Jones and Gallet values 1962) tends to occur over the Equatorial Anomaly close to the solar terminator sector.

Summary and conclusion
This work presents a technique that we have developed to routinely update the ITU-R F2 peak dataset, using RO electron density profiles derived from GPS receivers flying on low Earth orbiter satellites.The technique is used to compute two sets of coefficients for computing monthly mean values of N m F2 and h m F2, for any given latitude, longitude and UT, and for the current solar activity level.Along with the development we have slightly modified the mathematical formulation recommended by the ITU-R so that IRI and NeQuick models could implement our dataset without the need of significant changes on the source codes.
A critical point of our work was the quality evaluation of the RO profiles (i.e. the acceptance of a profile as a valid one) that are the primary input of our technique.After quality filtering the input data, every RO profile available for the analyzed period was fitted to the LPIM profile using a Least Squares adjustment and an iteratively re-weighting scheme.The previous step provided N m F2 and h m F2 values accompanied with estimations of their errors (standard deviations) for every RO profile.These set of values were used to compute, in a monthly basis and with a Least Squares method, two sets of numerical coefficients for the global mapping the N m F2 and h m F2.The geographical variation of the parameters was modeled with a spherical harmonic expansion series and the diurnal variation with Fourier expansion series.The maximum degrees of these expansions were established through an ANOVA test, which led to maximum degrees equal to 12 (geographical variation) and 9 (diurnal variation) for N m F2, and 6 and 8 for h m F2.
We tested our technique using the RO profiles produced by the COSMIC/FormoSat 3 team, for the 4 months (solstices and equinoxes) of one low and one high solar activity year (each month at a time).The RMS of the differences between RO and modeled values did not exceed 3 % of the N m F2 value of the corresponding profiles We have also analyzed error propagation into the estimated parameters and found errors that did not exceed 1 %.The unit-weight standard deviation resulting from the Least Squares adjustments ranged from ±0.5 × 10 10 to ±1.2 × 10 10 elec/m −3 for N m F2, and from ±1.7 to ±3.4 km for h m F2.A χ 2 test was applied to these values and the result is that we cannot reject the proposed model (spherical harmonic and Fourier expansion series, Sect.3.2) with a significance level of 97.5 %.Furthermore, the error propagation law was conducted to estimate the errors of the maps, which ranged from 0.5 × 10 10 to 3.6 × 10 10 elec/m −3 for N m F2 (approximately 7 % of the mapped parameter) and from 2.0 to 5.6 km for h m F2 (approximately 2 % of the mapped parameter).
Regarding the differences of N m F2 and h m F2 between the values computed with the Jones and Gallet (1962) and our technique, we saw that the largest errors for both parameters are found in region of the Equatorial Anomaly.This fact could be related to difficulties to model the electron density distribution in that region and, up to some extent, to the limitations of the RO technique to accurate values in regions of steep horizontal gradients.The largest differences for the N m F2 ranged from approximately ±80 × 10 10 elec/m −3 for high solar activity to +20 × 10 10 /−50 × 10 10 elec/m −3 for low solar activity.The corresponding differences for the h m F2 maps are around ± 80 km for both, high and low solar activity conditions.

Fig. 1
Fig. 1 Global representation of the RMS difference between the observed and a posteriori computed electron density for the profiles within the 18-20 UT interval of a the September equinox of the year 2007; and b the December solstice of the year 2011

Fig. 2
Fig. 2 Global representation of the N m F2 estimated value and error for the profiles within the 18-20 UT interval: a N m F2 for the 2007 September equinox; b N m F2 error for the same period; c N m F2 for the 2011 December solstice; and d N m F2 error for the same period

Fig. 4
Fig. 4 Global map of a monthly median values of N m F2 computed with the Jones and Gallet technique, at 19 UT, for the September equinox of the year 2007; b monthly mean values of N m F2 derived with our technique from the profiles between the 18-20 UT of the September equinox of the year 2007; c error of our map; d monthly median value

Fig. 5
Fig. 5 Global map of: a monthly median values of h m F2 computed with the Jones and Gallet technique, at 19 UT, for the September equinox of the year 2007; b monthly mean value of h m F2 derived with our technique form the profiles within the 18-20 UT interval of the September equinox of the year 2007; c error of our map; d monthly median val- Codrescu et al. 2001;Hernandez-Pajares 2003;Delay and Doherty 2004;Brunini et al. 2005;Azpilicueta and Brunini 2008)u et al. 2001;Hernandez-Pajares 2003;Delay and Doherty 2004;Brunini et al. 2005;Azpilicueta and Brunini 2008).

Table 1
General description of the used dataset