A topside investigation over a mid-latitude digisonde station in Cyprus

We examine the systematic differences between topside electron density measurements recorded by different techniques over the low-middle latitude operating European station in Nicosia, Cyprus (geographical coordinates: 35.14^oN, 33.2^oE), (magnetic coordinates 31.86^oN, 111.83 ^oE). These techniques include space-based in-situ data by Langmuir probes on board.European Space Agency (ESA) Swarm satellites, radio occultation measurements on board low Earth orbit (LEO) satellites from the COSMIC/FORMOSAT-3 mission and ground-based extrapolated topside electron density profiles from manually scaled ionograms. The measurements are also compared with International Reference Ionosphere Model (IRI-2016) topside estimations and IRI-corrected NeQuick topside formulation (method proposed by Pezzopane and Pignalberi (2019)). The comparison of Swarm and COSMIC observations with digisonde and IRI estimations verifies that in the majority of cases digisonde underestimates while IRI overestimates Swarm observations but in general, IRI provides a better topside representation than the digisonde. For COSMIC and digisonde profiles matched at the F layer peak the digisonde systematically underestimates topside COSMIC electron density values and the relative difference between COSMIC and digisonde increases with altitude (above hmF2), while IRI overestimates the topside COSMIC electron density but after a certain altitude (~150 km above hmF2) this overestimation starts to decrease with altitude. The IRI-corrected NeQuick underestimates the majority of topside COSMIC electron density profiles and relative difference is lower up to approximately 100 km (above the hmF2) and then it increases. The overall performance of IRI-corrected NeQuick improves with respect to IRI and digisonde.     

Sporadic sodium layers and the average vertical distribution of atmospheric sodium

We have used atmospheric sodium profiles obtained at São José dos Campos over the period 1993–2004 to study the influence of sporadic layers on the average sodium distribution. All the profiles for the time period 19–22 h local standard time were individually classified as either normal or sporadic. Profiles were considered sporadic if they contained sharp increases above the background profile by more than a factor of 2, occurring on the topside of the layer, with a width of less than 5 km. For each night’s data, the averages of the normal and sporadic profiles were stored separately. A comparison between the long-term means of these two average profiles, for those days on which both normal and sporadic profiles were observed, showed no significant difference between them. This indicates that, on average, there is no additional sodium in the sporadic profiles. This analysis suggests that at least part of the topside of the normal sodium layer originates in a source that also gives rise to sporadic layers, which would occur when the source itself becomes layered. The most likely source is sodium-containing ions, which would be layered by the tide and/or gravity wave driven wind-shear mechanism generally accepted as being responsible for some types of ionospheric sporadic-E.     

Mapping seasonal trends of electron temperature in the topside ionosphere based on DEMETER data

The diurnal, seasonal and latitudinal variations of the electron temperature in the Earth‘s topside ionosphere during relatively low solar activity period of 2005 – 2008 are investigated. In order to examine seasonal variations and morphology of the topside ionospheric plasma temperature, CNES micro-satellite DEMETER ISL data are used. Presented study is oriented on the dataset gathered in 2005 and 2008. Within conducted analysis, global maps of electron temperature for months of equinoxes and solstices have been developed. Furthermore, simultaneous studies on two-dimensional time series based on DEMETER measurements and predictions obtained with the IRI-2012 model supply examination of the topside ionosphere during recent deep solar minimum. Comparison with the IRI-2012 model reveals discrepancies between data and prediction, that are especially prominent during the periods of very low solar activity.     

Adaptation of the bottomside electron density profile of the IRI model to data of topside radiosounding

A method is proposed for reconstructing the electron density profiles N(h) of the IRI model from ionograms of topside satellite sounding of the ionosphere. An ionograms feature is the presence of traces of signal reflection from the Earth's surface. The profile reconstruction is carried out in two stages. At the first stage, the N(h) –profile is calculated from the lower boundary of the ionosphere to the satellite height (total profile) by the method presented in this paper using the ionogram. In this case, the monotonic profile of the topside ionosphere is calculated by the classical method. The profile of the inner ionosphere is represented by analytical functions, the parameters of which are calculated by optimization methods using traces of signal reflection, both from the topside ionosphere and from the Earth. At the second stage, the profile calculated from the ionogram is used to obtain the key parameters: the height of the maximum hmF2 of the F2 layer, the critical frequency foF2, the values of B0 and B1, which determine the profile shape in the F region in the IRI model. The input of key parameters, time of observation, and coordinates of sounding into the IRI model allows obtaining the IRI-profile corrected to real experimental conditions. The results of using the data of the ISIS-2 satellite show that the profiles calculated from the ionograms and the IRI profiles corrected from them are close to each other in the inner ionosphere and can differ significantly in the topside ionosphere. This indicates the possibility of obtaining a profile in the inner ionosphere close to the real distribution, which can significantly expand the information database useful for the IRTAM (IRI Realmax Assimilative Modeling) model. The calculated profiles can be used independently for local ionospheric research.     

Reconstruction of topside density profile by using the topside sounder model profiler and digisonde data

To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (H_s), O⁺–H⁺ transition height (H_T), and their ratio Rt = H_s/H_T, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.     

Storm time behavior of topside scale height inferred from the ionosphere–plasmasphere model driven by the F2 layer peak and GPS-TEC observations

The International Reference Ionosphere model extended to the plasmasphere, IRI-Plas, presents global electron density profiles and total electron content, TECiri, up to the altitude of the GPS satellites (20,000 km). The model code is modified by input of GPS-derived total electron content, TECgps, so that the topside scale height, Hsc, is obtained minimizing in one step the difference between TECiri and TECgps observation. The topside basis scale height, Hsc, presents the distance in km above the peak height at which the peak plasma density, NmF2, decays by a factor of e (∼2.718). The ionosonde derived F2 layer peak density and height and GPS-derived TECgps data are used with IRI-Plas code during the main phase of more than 100 space weather storms for a period of 1999–2006. Data of seven stations are used for the analysis, and data from five other stations served as testing database. It is found that the topside basis scale height is growing (depressing) when the peak electron density (critical frequency foF2) and electron content are decreased (increased) compared to the median value, and vice versa. Relative variability of the scale height, rHsc, and the instantaneous Hsc are inferred analytically in a function of the instantaneous foF2, median fmF2 and median Hmsc avoiding a reference to geomagnetic indices. Results of validation suggest reliability of proposed algorithm for implementation in an operational mode.     

Examining the use of the NeQuick bottomside and topside parameterizations at high latitudes

An examination of the high latitude performance of the bottomside and topside F-layer parameterizations of the NeQuick electron density model is presented using measurements from high latitude ionosonde and Incoherent Scatter Radar (ISR) facilities.For the bottomside, we present a comparison between modeled and measured B2Bot thickness parameter. In this comparison, it is seen that the use of the NeQuick parameterization at high latitudes results in significantly underestimated bottomside thicknesses, regularly exceeding 50%. We show that these errors can be attributed to two main issues in the NeQuick parameterization:(1) through the relationship relating foF2 and M3000F2 to the maximum derivative of F2 electron density, which is used to derive the bottomside thickness, and (2) through a fundamental inability of a constant thickness parameter, semi-Epstein shape function to fit the curvature of the high latitude F-region electron density profile.For the topside, a comparison is undertaken between the NeQuick topside thickness parameterization, using measured and CCIR-modeled ionospheric parameters, and that derived from fitting the NeQuick topside function to Incoherent Scatter Radar-measured topside electron density profiles. Through this comparison, we show that using CCIR-derived foF2 and M3000F2, used in both the NeQuick and IRI, results in significantly underestimated topside thickness during summer periods, overestimated thickness during winter periods, and an overall tendency to underestimate diurnal, seasonal, and solar cycle variability. These issues see no improvement through the use of measured foF2 and M(3000)F2 values. Such measured parameters result in a tendency for the parametrization to produce a declining trend in topside thickness with increasing solar activity, to produce damped seasonal variations, and to produce significantly overestimated topside thickness during winter periods.     

Validate the IRI2007 model by the COSMIC slant TEC data during the extremely solar minimum of 2008

During 2008, the solar activity is extremely low. The satellite observations show that the ionospheric height and electron density is much lower than the predictions by the international reference ionosphere (IRI) model. In this paper, we compared the slant total electron content (TEC) observed by the COSMIC satellites during 2008 with the IRI model results. It is found that the IRI model with IRI2001 and IRI2001 Cor. topside options will always overestimate the electron density in both lower and higher altitudes. But the rest two topside options (NeQuick, and TTS) tend to overestimate the electron density in the F layer and underestimate it in the topside altitudes. The switch altitude between overestimation and underestimation and the latitude-local time distribution of the model deviation depend on the topside option. The current investigation might be useful for the model improvement as well as data assimilation work based on the IRI model and the LEO TEC data.     

Evaluating the IRI topside model for the South African region: An overview of the modelling techniques

The representation of the topside ionosphere (the region above the F2 peak) is critical because of the limited experimental data available. Over the years, a wide range of models have been developed in an effort to represent the behaviour and the shape of the electron density (Ne) profile of the topside ionosphere. Various studies have been centred around calculating the vertical scale height (VSH) and have included (a) obtaining VSH from Global Positioning System (GPS) derived total electron content (TEC), (b) calculating the VSH from ground-based ionosonde measurements, (c) using topside sounder vertical Ne profiles to obtain the VSH. One or a combination of the topside profilers (Chapman function, exponential function, sech-squared (Epstein) function, and/or parabolic function) is then used to reconstruct the topside Ne profile. The different approaches and the modelling techniques are discussed with a view to identifying the most adequate approach to apply to the South African region’s topside modelling efforts. The IRI-2001 topside model is evaluated based on how well it reproduces measured topside profiles over the South African region. This study is a first step in the process of developing a South African topside ionosphere model.     

Comparing TOPEX TEC measurements with IRI predictions

TEC values obtained from TOPEX satellite were compared with the International Reference Ionosphere (IRI) 2001 model estimates. The present work also shows results of the IRI model with the option of a new topside electron density distribution (NeQuick model). TOPEX TEC measurements, which include years of high and middle to low solar activity (2000 and 2004), were analyzed by binning the region covered by the satellite (±66°) every five degrees of modip. In general, there is good agreement between IRI predictions and Topex measurements. Cases with large disagreements are observed at low and high latitudes during high solar activity. Comparing the model predictions using the default IRI2001 model and the NeQuick topside option show that the default IRI 2001 version represents the observed data in a more realistic way, but appears to be less reliable at high and low latitudes in some cases.     

On the relative abundance of helium ions in the topside ionosphere

Based on data from satellite INTERCOSMOS-BULGARIA-1300, the latitudinal distribution of oxygen and helium ions in the topside ionosphere is discussed for night-time equinox at high solar activity. A comparison with the corresponding IRI-79 distribution is made. The vertical IRI ion composition profile is checked with measurements made with VERTICAL-10 rocket. Some recommendations are made in order to improve the IRI-modelling of the ion composition in the topside ionosphere.     

Evaluation of the IRI-2007 model options for the topside electron density

The International Reference Ionosphere (IRI) 2007 provides two new options for the topside electron density profile: (a) a correction of the IRI-2001 model, and (b) the NeQuick topside formula. We use the large volume of Alouette 1, 2 and ISIS 1, 2 topside sounder data to evaluate these two new options with special emphasis on the uppermost topside where IRI-2001 showed the largest discrepancies. We will also study the accurate representation of profiles in the equatorial anomaly region where the profile function has to accommodate two latitudinal maxima (crests) at lower altitudes but only a single maximum (at the equator) higher up. In addition to IRI-2001 and the two new IRI-2007 options we also include the Intercosmos-based topside model of Triskova, Truhlik, and Smilauer [Triskova, L., Truhlik, V., Smilauer, J. An empirical topside electron density model for calculation of absolute ion densities in IRI. Adv. Space Res. 37 (5), 928–934, 2006] (TTS model) in our analysis. We find that overall IRI-2007-NeQ gives the best results but IRI-2007-corrected provides a more realistic representation of the altitudinal–latitudinal structure in the equatorial anomaly region. The applicability of the TTS model is limited by the fact that it is not normalized to the F2 peak density and height.     

Modelling bottom and topside electron density and TEC with profile data from topside ionograms

The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.     

基于Volterra级数的功放行为建模与简化

Volterra级数行为模型因其高复杂度而只适用于弱非线性系统,对于具有强非线性和长记忆效应的功率放大器,Volterra级数模型将包含大量的待辨识参数.对射频功率放大器的行为建模进行了研究,提出了功率放大器Volterra级数行为模型的简化方法,通过引入邻对角核控制因子和动态偏离阶数控制因子,合理移除了对信号输出影响较少甚至无关的辨识参数,从而在保证模型准确度的前提下有效地降低了模型的复杂度和参数提取时间.采用了一个10 W的射频功率放大器实例用于模型提取和验证,模型与仿真实验的数据结果对比表明,简化的Volterra行为模型能够很好地描述功放的非线性特性及记忆效应,验证了简化行为模型的有效性和灵活性.     

On the use of NeQuick topside option in IRI-2007

The latest version of IRI includes various options for the computation of the topside electron density profile. One of the possible choices is based on NeQuick model. Its inclusion in IRI has been made transferring all the formulations used in NeQuick model. In details, an Epstein layer function is used to describe the electron density profile and the topside shape is controlled by an empirical parameter, connected to the NeQuick F2 bottomside thickness parameter, B2bot. It is computed also in this IRI topside option in order to maintain self-consistency with its original formulation. This paper analyses the possibility of using the IRI bottomside parameters for this option and its impact on the profile and TEC. The case of experimental peak values given as input is also analysed.     

A simplified D-region ion chemistry scheme and its possible use for IRI lower ionosphere modelling

Ion composition of the D region is principally characterized by the existence of two distinct regions of predominant molecular ions and predominant cluster ions, separated from each other by a rather sharp ‘transition height’, which is proposed to be included in the IRI as an additional parameter, supplementing the electron density models. It is possible to predict the position of this ‘transition height’ at a given place and time with the aid of a simplified ion chemistry scheme which is shown to be satisfactorily compatible with experimental ion composition data available in the literature. Our suggested method of this prediction makes use of the (IRI or experimental) electron density profile at the location and season in question, together with an effective clustering rate coeeficient calculated from corresponding temperature and density profiles taken from a suitable reference model of the neutral atmosphere.     

Neural network prediction of the topside electron content over the Euro-African sector derived from Swarm-A measurements

This study presents the first prediction results of a neural network model for the vertical total electron content of the topside ionosphere based on Swarm-A measurements. The model was trained on 5 years of Swarm-A data over the Euro-African sector spanning the period 1 January 2014 to 31 December 2018. The Swarm-A data was combined with solar and geomagnetic indices to train the NN model. The Swarm-A data of 1 January to 30 September 2019 was used to test the performance of the neural network. The data was divided into two main categories: most quiet and most disturbed days of each month. Each category was subdivided into two sub-categories according to the Swarm-A trajectory i.e. whether it was ascending or descending in order to accommodate the change in local time when the satellite traverses the poles. Four pairs of neural network models were implemented, the first of each pair having one hidden layer, and the second of each pair having two hidden layers, for the following cases: 1) quiet day-ascending, 2) quiet day-descending, 3) disturbed day-ascending, and 4) disturbed day-descending. The topside vertical total electron content predicted by the neural network models compared well with the measurements by Swarm-A. The model that performed best was the one hidden layer model in the case of quiet days for descending trajectories, with RMSE = 1.20 TECU, R = 0.76. The worst performance occurred during the disturbed descending trajectories where the one hidden layer model had the worst RMSE = 2.12 TECU, (R = 0.54), and the two hidden layer model had the worst correlation coefficient R = 0.47 (RMSE = 1.57).In all cases, the neural network models performed better than the IRI2016 model in predicting the topside total electron content. The NN models presented here is the first such attempt at comparing NN models for the topside VTEC based on Swarm-A measurements.     

Towards a new reference model of hmF2 for IRI

A numerical model of the peak height of the F2 layer, hmF2_top, is derived from the topside sounding database of 90,000 electron density profiles for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The model of regular hmF2 variations is produced in terms of local time, season, geomagnetic latitude, geodetic longitude and solar radio flux. No geomagnetic activity trends were discernible in the topside sounding data. The nighttime peak of hmF2_top evident for mid-latitudes disappears near the geomagnetic equator where a maximum of hmF2_top occurs at sunset hours when it can exceed 500 km at solar maximum. The hmF2 given by the IRI exceeds hmF2_top at the low solar activities. The hmF2_top, obtained by extrapolation of the first derivative of the topside profile to zero shows saturation similar to foF2 the greater the solar activity. The proposed model differs from hmF2 given by IRI based on M(3000)F2 to hmF2 conversion by empirical relationships in terms of foF2, foE and R12 with these quantities mapped globally by the ITU-R (former CCIR) from ground-based ionosonde data. The differences can be attributed to the different techniques of the peak height derivation, different epochs and different global distribution of the source data as well as the different mathematical functions involved in the maps and the model presentation.     

赤道地区和中纬度地区扩展F的特点和比较

通过数值模拟研究了赤道地区和中纬度地区扩展F的触发机制间的联系,发现相同的触发机制在不同纬度的条件下所演化出的结构形态不相同,因此观测上会有所不同;研究了西向电场在赤道地区的作用,发现在适合的电离层参数组合下会产生观测到的电离层F区顶部扰动出现率高于底部扰动出现率的状况;在数值模拟的基础上提出了一个中纬度地区电波散射的模型,用以解释观测到的中纬度地区电波散射的出现率高于线性理论预期的现象.