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.     

Comparing topside and bottomside-measured characteristics of the F2 layer peak

The ionospheric characteristics of the F2 layer peak have been measured with ionosondes from the ground or with satellites from space. The most common characteristics are the F2-peak density NmF2 and peak height hmF2. In addition to these two parameters this paper studies the F2-peak scale height. Comparing the median values of hmF2 and NmF2 obtained from topside and bottomside sounding shows good agreement in general. The Chapman scale height values for the F2 layer peak derived from topside profiles, H_m,top, are generally several times larger than H_m,bot derived from bottomside profiles.     

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.     

Effects of solar activity variations on the low latitude topside nighttime ionosphere

We have studied the topside nighttime ionosphere of the low latitude region using data obtained from DMSP F15, ROCSAT-1, KOMPSAT-1, and GUVI on the TIMED satellite for the period of 2000–2004, during which solar activity decreased from its maximum. As these satellites operated at different altitudes, we were able to discriminate altitude dependence of several key ionospheric parameters on the level of solar activity. For example, with intensifying solar activity, electron density was seen to increase more rapidly at higher altitudes than at lower altitudes, implying that the corresponding scale height also increased. The density increased without saturation at all observed altitudes when plotted against solar EUV flux instead of F10.7. The results of the present study, as compared with those of previous studies for lower altitudes, indicate that topside vertical scale height increases with altitude and that, when solar activity increases, topside vertical scale height increases more rapidly at higher altitudes than at lower altitudes. Temperature also increased more rapidly at higher altitudes than at lower altitudes as solar activity increased. In addition, the height of the F2 peak was seen to increase with increasing solar activity, along with the oxygen ion fraction measured above the F2 peak. These results confirm that the topside ionosphere rises and expands with increasing solar activity.     

The asymmetrical features in electron density during extreme solar minimum

The variations of plasma density in topside ionosphere during 23rd/24th solar cycle minimum attract more attentions in recently years. In this analysis, we use the data of electron density (Ne) from DEMETER (Detection of Electromagnetic Emissions Transmitted from Earthquake Regions) satellite at the altitude of 660–710 km to investigate the solstitial and equinoctial asymmetry under geomagnetic coordinate system at LT (local time) 1030 and 2230 during 2005–2010, especially in solar minimum years of 2008–2009. The results reveal that ΔNe (December–June) is always positive over Southern Hemisphere and negative over northern part whatever at LT 1030 or 2230, only at 0–10°N the winter anomaly occurs with ΔNe (December–June) > 0, and its amplitude becomes smaller with the declining of solar flux from 2005 to 2009. The ΔNe between September and March is completely negative during 2005–2008, but in 2009, it turns to be positive at latitudes of 20°S–40°N at LT 1030 and 10°S–20°N at LT 2230. Furthermore, the solstitial and equinoctial asymmetry index (AI) are calculated and studied respectively, which all depends on local time, latitude and longitude. The notable differences occur at higher latitudes in solar minimum year of 2009 with those in 2005–2008. The equinoctial AI at LT 2230 is quite consistent with the variational trend of solar flux with the lowest absolute AI occurring in 2009, the extreme solar minimum, but the solstitial AI exhibits abnormal enhancement during 2008 and 2009 with bigger AI than those in 2005–2007. Compared with the neutral compositions at 500 km altitude, it illustrates that [O/N₂] and [O] play some roles in daytime and nighttime asymmetry of Ne at topside ionosphere.     

Global validation of IRI TEC for high and medium solar activity conditions

The IRI model offers a choice of options for the computation of the electron density profile and electron content (TEC). Recently new options for the topside electron density profile have been developed, which have a strong impact on TEC. Therefore it is important to test massively the IRI and the new options with experimental data. A large number of permanent stations record dual frequency GPS data from which it is possible to obtain TEC values. Thirty-one worldwide distributed stations have been selected to investigate the capabilities of the IRI to reproduce experimental TEC. Data for years 2000 (high solar activity) and 2004 (medium solar activity) have been analyzed computing modeled values with the IRI-2001 and the IRI-2007-NeQuick topside options. It is found that IRI-2007-NeQuick option generally improves the estimate of the slant TEC, especially in the case of high latitudes stations during high solar activity.     

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.     

Four-peak longitudinal distribution of the equatorial plasma bubbles observed in the topside ionosphere: Possible troposphere tide influence

In this paper we consider an idea of the troposphere tide influence on the character of the longitudinal variations in the distribution of the equatorial plasma bubbles (EPBs) observed in the topside ionosphere. For this purpose, the obtained EPB longitudinal patterns were compared with the thermosphere and ionosphere characteristics having the prominent “wave-like” longitudinal structures with wave number 4, which are uniquely associated with the influence of the troposphere DE3 tides. The characteristics of the equatorial mass density anomaly (EMA), equatorial ionization anomaly (EIA), zonal wind and pre-reversal E?×?B drift enhancement (PRE) were used for comparison. The equinox seasons during high solar activity were under consideration. It was obtained that the longitudinal patterns of the EMA and zonal wind show the surprising similarity with the EPB distributions (R???0.8, R???0.72). On the other hand, the resemblance with the ionosphere characteristics (EIA, PRE) is rather faint (R???0.37, R???0.12). It was shown that the thermosphere zonal winds are the most possible transfer mediator of the troposphere DE3 tide influence. The most successful moment for the transfer of the troposphere DE3 tide energy takes place in the beginning of the EPB production, namely, during the seed perturbation development.     

利用COSMIC数据研究顶部电离层O+扩散流的统计特征

从等离子体运动方程出发, 利用COSMIC星座的掩星数据, 借助相关经验模式, 计算了太阳活动低年顶部电离层O+场向扩散速度和扩散通量, 并分析了其全球分布和日变化特征. 结果表明, 白天等离子体扩散速度的方向随高度增加由向下(极向)逐步变为向上(赤道向), 方向转变的高度一般在hmF2+80 km以下. 在白天较高高度, 南北磁纬10o ~20 o存在着向上方向的最大扩散速度和扩散通量; 而在夜间, 南北磁纬30o~40 o存在向下方向的最大扩散速度和扩散通量. 在分点, 南北半球的扩散通量和扩散速度大致对称; 而在至点, 扩散通量存在着明显的南北半球不对称现象. 另外, 不同纬度的扩散速度有着不同的日变化特征.      

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.     

Analysis of a topside ionospheric model using GPS and ionosonde observables

Ground-based vertical incidence soundings are well suited to model the bottom-side ionosphere but are not so good for dependably modelling the topside ionosphere. This study aims to combine vertical incidence sounding and dual-frequency GPS measurements to reconstruct the topside profile. The reconstruction technique relays on the use of the so-called vary-Chap approach that use an α-Chapman function with a continuously varying scale height.     

An analysis of the topside ionosphere parameters based on the long-duration Irkutsk incoherent scatter radar measurements

The topside ionosphere parameters are studied based on the long-duration Irkutsk incoherent scatter radar (52.9N, 103.3E) measurements conducted in September 2005, June and December 2007. As a topside ionosphere parameter we chose the vertical scale height (VSH) related to the gradient of the electron density logarithm above the peak height. For morphological studies we used median electron density profiles. Besides the median behavior we also studied VSH disturbances (deviations from median values) during the magnetic storm of September 11th 2005. We compared the Irkutsk incoherent scatter radar data with the Millstone Hill and Arecibo incoherent scatter radar observations, the IRI-2007 prediction (using the two topside options) and VSH derived from the Irkutsk DPS-4 Digisonde bottomside measurements.     

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.     

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.     

Response of the mid-latitude D-region ionosphere to the total solar eclipse of 22 July 2009 studied using VLF signals in South Korean peninsula

In this paper, we analyze VLF signals received at Busan to study the the D-region changes linked with the solar eclipse event of 22 July 2009 for very short (∼390 km) transmitter–receiver great circle path (TRGCP) during local noon time 00:36–03:13 UT (09:36–12:13 KST). The eclipse crossed south of Busan with a maximum obscuration of ∼84%. Observations clearly show a reduction of ∼6.2 dB in the VLF signal strength at the time of maximum solar obscuration (84% at 01:53 UT) as compared to those observed on the control days. Estimated values of change in Wait ionospheric parameters: reflection height (h′) in km and inverse scale height parameter (β) in km⁻¹ from Long Wave Propagation Capability (LWPC) model during the maximum eclipse phase as compared to unperturbed ionosphere are 7 km and 0.055 km⁻¹, respectively. Moreover, the D-region electron density estimated from model computation shows 95% depletion in electron density at the height of ∼71 km. The reflection height is found to increase by ∼7 km in the D-region during the eclipse as compared to those on the control days, implying a depletion in the Lyman-α flux by a factor of ∼7. The present observations are discussed in the light of current understanding on the solar eclipse induced D-region dynamics.     

Latitudinal variation of the topside electron temperature at different levels of solar activity

A database of electron temperature (T_e) measurements comprising of most of the available satellite measurements in the topside ionosphere is used for studying the solar activity variations of the electron temperature T_e at different latitudes, altitudes, local times and seasons. The T_e data are grouped into three levels of solar activity (low, medium, high) at four altitude ranges, for day and night, and for equinox and solstices. We find that in general T_e changes with solar activity are small and comparable in magnitude with seasonal changes but much smaller than the changes with altitude, latitude, and from day to night. In all cases, except at low altitude during daytime, T_e increases with increasing solar activity. But this increase is not linear as assumed in most empirical T_e models but requires at least a parabolic approximation. At 550 km during daytime negative as well as positive correlation is found with solar activity. Our global data base allows to quantify the latitude range and seasonal conditions for which these correlations occur. A negative correlation with solar activity is found in the invdip latitude range from 20 to 55 degrees during equinox and from 20 degrees onward during winter. In the low latitude (20 to −20 degrees invdip) F-region there is almost no change with solar activity during solstice and a positive correlation during equinox. A positive correlation is also observed during summer from 30 degrees onward.     

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.     

Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.     

Estimation of the Love and Shida numbers: h2, l2 using SLR data for the low satellites

In this paper we present results for the global elastic parameters: Love number h₂ and Shida number l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two low satellites STELLA (H = 800 km) and STARLETTE (H = 810 km) observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analyzed. The analysis was done separately for the two satellites. We do a sequential analysis and study the stability and convergence of the estimates as a function of length of the data set used.     

Variation of F-region critical frequency (foF2) over equatorial and low-latitude region of the Indian zone during 19th and 20th solar cycle

The effect of solar cycle and seasons on the daytime and nighttime F-layer ionization has been investigated over the equatorial and low-latitude region during 19th (1954–1964) and 20th (1965–1976) solar cycle. The F-layer critical frequency (foF2) data observed from the three Indian Ionosonde stations has been used for the present study. The dependence of foF2 on solar cycle has been examined by performing regression analysis between the foF2 values and R₁₂ (twelve month running average sunspot number). The result shows that the magnitude of the cycle, seasons and the location of station has considerable effects on foF2. There is a significant nonlinear relationship between the foF2 values and R₁₂ during 19th solar cycle as compared to 20th solar cycle. Further, the nighttime saturation effect is prominently seen during the 19th solar cycle and summer season. It is also observed that the most profound saturation effect appears at the equatorial ionization anomaly crest region. Seasonally, it is seen that all the stations exhibits semiannual anomaly. The phenomenon of winter anomaly decays as we move higher along the latitude and is prominently seen during the intense solar activity.