Annals 48, 3, 2005+app1 491 ANNALS OF GEOPHYSICS, VOL. 48, N. 3, June 2005 Key words topside ionosphere – ionospheric mod- eling – ionospheric layer formulation – topside soundings 1. Introduction Different kinds of analytical functions are being used to model the electron density profile of the various ionospheric layers. Since the ex- perimental information of the topside electron density profile is very scarce, this part of the modeled profiles has been less tested against measured data. In the last few years a relatively large amount of topside electron density profiles has been made available to the scientific communi- ty through the Internet allowing to test models (Coïsson et al., 2002) and the different mathe- matical representations used. This work analyzes the ability of some ionos- pheric layers analytical formulation to reproduce the experimental profiles measured by some old satellites which carried topside sounders: the Russian Intercosmos-19 (IK19) and the Canadi- an-American International Satellite for Ionos- phere Studies 2 (ISIS2). Some analytical aspects (derivative, inflec- tion point and scale parameter) of common lay- er functions are investigated here in order to find better ways to model the lower part of the ionosphere topside. 2. Layer functions for the topside Some simple layer functions have been cho- sen to perform the comparisons: Epstein, Chap- man and other derived formulations. They are characterized by a scale parameter which deter- mines the thickness of the layer. Epstein layer function ( , ) ( / ) ( / ) exp exp h B h B h B Epstein 1 4 2 2 2 2 = +^ h (2.1) h being the height relative to the F2 maximum and B2 the scale parameter for the layer. Comparison of analytical functions used to describe topside electron density profiles with satellite data Carlo Fonda, Pierdavide Coïsson, Bruno Nava and Sandro M. Radicella The Abdus Salam International Centre for Theoretical Physics (ICTP), Trieste, Italy Abstract Electron density models of the ionosphere use different analytical formulations for the electron density vertical profile in the topside. The present paper compares some single-layer topside analytical descriptions (Chapman, Epstein, modified Epstein used in the NeQuick model) with experimental topside profiles obtained from meas- urements of IK19 and ISIS2 satellites. The limits of height range and shape for each formulation are described and analyzed and suggestions for the use of multiple layers solution to reproduce experimental results are given. Mailing address: Carlo Fonda, The Abdus Salam Inter- national Centre for Theoretical Physics (ICTP), Strada Co- stiera 11, 34014 Trieste, Italy; e-mail: cfonda@ictp.it 492 Carlo Fonda, Pierdavide Coïsson, Bruno Nava and Sandro M. Radicella «Modified Epstein» layer function ( , ) exp exp h H k k NeQuick 1 4 0 2= +^ h (2.2) . . k H h h H h 1 0 125 100 12 5 0 0 = + + b l (2.3) h being the height relative to the F2 maximum and H0 the scale parameter for the layer. It is used in the NeQuick model (Radicella and Leitinger, 2001). Chapman family of layer functions ( , , ) exp exph k k h k h Chapman 1= - - -a a bb l l; E (2.4) h being the height relative to the F2 maximum, k the scale parameter for the layer and α the shape parameter. This function is usually called according to the value of α: Chapman-alpha if α = 0.5 and Chapman-beta if α = 1. Examples of the described functions shapes are shown in fig. 1 in linear scale and in fig. 2 in logarithmic scale. All the profiles have been made to have an electron density of 0.001 NmF2 at 1000 km above hmF2. 3. Data from IK19 satellite The first set of data is the one from the IK19 satellite. This satellite was on an elliptical orbit with 74° inclination and it had perigee at 500 km and apogee at 980 km (Pulinets, 1989). The IK19 topside profiles cover height ranges form about 200 km up to about 750 km over the F2 peak. It is a data set which allows to analyze the lower part of the topside. The data available are from the high solar activity period March 1979-De- cember 1980. α-Chapman, β-Chapman, Epstein and modi- fied Epstein (NeQuick) layer functions were used to fit 7138 experimental electron density profiles derived from IK19 on board topside sounder. The fit was performed using least square methods, some examples are presented in figs. 3, 4 and 5. In addition the following parameters have also been evaluated: 1) height of the inflection point; 2) absolute difference between the integral of the fitting function and the integral of the experimen- tal profile, numerically obtained. 3.1. Analysis of integrals The four considered modeling functions have been analyzed, comparing for each profile the absolute difference ε between the integral of Fig. 1. Profile shapes for the different functions with intersection at 1000 km above hmF2: Chapman-alpha dashed, Chapman-beta dotted, Epstein dash-dotted, NeQuick solid gray. Linear scale. Fig. 2. Profile shapes for the different functions with intersection at 1000 km above hmF2: Chapman-alpha dashed, Chapman-beta dotted, Epstein dash-dotted, NeQuick solid gray. Logarithmic scale. 1 2 Fig. 5. IK19 experimental profile (solid black) and the fitted functions: Chapman-alpha dashed, Chapman-be- ta dotted, Epstein dash-dotted, NeQuick solid gray. Fig. 6. Distribution of ε for all IK19 profiles. Chapman-alpha dashed, Chapman-beta dotted, Epstein dash-dot- ted, NeQuick solid gray. 5 6 493 Comparison of analytical functions used to describe topside electron density profiles with satellite data the fitting function and the integral of the ex- perimental profile. This criterion has been used to evaluate how well the shape of these layer functions reproduce the behavior of experimen- tal profiles, in order to identify which function is better able to model this part of the topside. The obtained distribution, presented in fig. 6, shows that a Chapman function with the shape parameter α = 0.5 gives th e best result, as confirmed by simply looking at most of the profile plots. 3.2. Analysis of inflection points The height of the inflection point is used as an additional criterion to evaluate the capability of a function to reproduce the shape of the lowest part Fig. 3. IK19 experimental profile (solid black) and the fitted functions: Chapman-alpha dashed, Chapman-be- ta dotted, Epstein dash-dotted, NeQuick solid gray. Fig. 4. IK19 experimental profile (solid black) and the fitted functions: Chapman-alpha dashed, Chapman-be- ta dotted, Epstein dash-dotted, NeQuick solid gray. 3 4 Fig. 9. Distribution of difference between experimental and fitted inflection point heights of IK19 profiles. Chapman-alpha dashed, Chapman-beta dotted, Epstein dash-dotted. Fig. 10. ISIS2 experimental profile (solid black) and the fitted functions: Chapman-alpha dashed, Chapman- beta dotted, Epstein dash-dotted, NeQuick solid gray. 9 10 494 Carlo Fonda, Pierdavide Coïsson, Bruno Nava and Sandro M. Radicella of the topside (IK19 profiles usually have a max- imum height in the order of 500 km above the F2 peak). The distribution of these heights (fig. 7), compared to the positions of the inflection points as computed from the experimental profiles (fig. 8), again shows a better performance of the Chap- man-alpha with respect to the other functions (fig. 9). The inflection points cannot be analytically de- termined for the NeQuick topside function and they are not included in such analysis. 4. Data from ISIS2 satellite The ISIS2 satellite was in operation for more than a decade, it had a circular orbit at 1400 km with inclination 88° (Bilitza et al., 2002, 2003; Huang et al., 2002). The same method of profile fitting has been applied to the experimental data set derived from ISIS2 satel- lite, some examples are shown in figs. 10, 11 and 12. The results are consistent with those ob- Fig. 7. Distribution of inflection point heights of fitted IK19 profiles: Chapman-alpha dashed, Chapman-beta dotted, Epstein dash-dotted. Fig. 8. Distribution of inflection point heights of IK19 profiles. 7 8 495 Comparison of analytical functions used to describe topside electron density profiles with satellite data tained from IK19 profiles. The main difference being the average maximum height of the pro- files, higher for ISIS2, these results show better the behavior of model functions in a higher part of the topside (range of height: 400-1000 km above the F2 peak). The fits analysis also shows that the layer function used for the lower topside cannot be used «as it is» for this higher range of heights, but a new and different 2-layers topside formu- lation, at least, is needed in order to better re- produce the experimental shape. The three selected profiles shown here are plotted with a logarithmic vertical scale to clearly indicate the different behavior of the model functions and the experimental profiles for this higher range of heights. 5. Conclusions A few analytical aspects of two common- used topside layer functions, with some vari- ants were studied with a comparison among them and experimental profiles, using a least square fit approach. It appeared that all the con- sidered functions have difficulties in reproduc- ing experimental data in the required range of height. It is also evident that a different analyt- ical formulation is needed in order to reproduce the relevant part of the topside. This new for- mulation may possibly involve a Chapman lay- er function with α < 0.5 for heights lower than 400 km above the F2 maximum and should use a different functional approach for the higher part, up to 1000 km above hmF2. REFERENCES BILITZA, D., X. HUANG, B. REINISH and R.F. BENSON (2002): TOPIST - Automated processing of ISIS top- side ionograms, in Proceedings of 2002 Ionospheric Effects Symposium, 333-340. BILITZA, D., B. REINISH, R. BENSON, J. GREBOWSKY, N. PA- PISTASHVILI, X. HUANG, W. SCHAR and K. HILLS (2003): Online database of satellite sounder and insitu meas- urements covering two solar cycles, Adv. Space Res., 31 (3), 769-774. COÏSSON, P., S.M. RADICELLA and B. NAVA (2002): Compar- isons of experimental topside electron concentration profiles with IRI and NeQuick models, Ann. Geo- physics, 45 (1), 111-116. HUANG, X., B.W. REINISCH, D. BILITZA and R.F. Benson (2002): Electron density profiles of the topside iono- sphere, Ann. Geophysics, 45 (1), 125-130. PULINETS, S.A. (1989): Prospect of topside sounding, in WITS Handbook (SCOSTEP Publishing, Urbana, IL), vol. 2, ch. 3, 99-127. RADICELLA, S.M. and R. LEITINGER (2001): The evolution of the DGR approach to model electron density pro- files, Adv. Space Res., 27 (1), 35-40. 11 12 Fig. 11. ISIS2 experimental profile (solid black) and the fitted functions: Chapman-alpha dashed, Chapman- beta dotted, Epstein dash-dotted, NeQuick solid gray. Fig. 12. ISIS2 experimental profile (solid black) and the fitted functions: Chapman-alpha dashed, Chapman- beta dotted, Epstein dash-dotted, NeQuick solid gray.