A new inversion algorithm for HF sky-wave backscatter ionograms

HF sky-wave backscatter sounding system is capable of measuring the large-scale, two-dimensional (2-D) distributions of ionospheric electron density. The leading edge (LE) of a backscatter ionogram (BSI) is widely used for ionospheric inversion since it is hardly affected by any factors other than ionospheric electron density. Traditional BSI inversion methods have failed to distinguish LEs associated with different ionospheric layers, and simply utilize the minimum group path of each operating frequency, which generally corresponds to the LE associated with the F₂ layer. Consequently, while the inversion results can provide accurate profiles of the F region below the F₂ peak, the diagnostics may not be so effective for other ionospheric layers. In order to resolve this issue, we present a new BSI inversion method using LEs associated with different layers, which can further improve the accuracy of electron density distribution, especially the profile of the ionospheric layers below the F₂ region. The efficiency of the algorithm is evaluated by computing the mean and the standard deviation of the differences between inverted parameter values and true values obtained from both vertical and oblique incidence sounding. Test results clearly manifest that the method we have developed outputs more accurate electron density profiles due to improvements to acquire the profiles of the layers below the F₂ region. Our study can further improve the current BSI inversion methods on the reconstruction of 2-D electron density distribution in a vertical plane aligned with the direction of sounding.     

返回散射电离图的前沿提取方法

提出了返回散射电离图前沿提取的方法. 其中, 一层电离图前沿和两层电离图下面一层前沿是在消除干扰和噪声后的电离图上采用边缘识别方法进行提取的, 两层电离图上面一层前沿则是在仅经过干扰抑制的电离图上采用基于高斯型加权系数的能量梯度算法进行提取的. 本文提出的前沿提取方法准确率较高, 并且易于实现, 已在工程中得到应用.      

Physical and model interpretation of HF radio propagation on the St. Petersburg–Longyearbyen (Svalbard) path

HF radio wave observations have been carried out with an oblique ionospheric sounding (OIS) method on the radio path from St. Petersburg to Longyearbyen (Svalbard), and experimental ionograms were obtained for December 2001. These ionograms have been analysed to investigate the impact of the main ionospheric trough (MIT) and magnetic disturbances on the signals on this path. The observations during weakly disturbed (K_р = 2) magnetic conditions on 14–15 December 2001 were compared with predictions from ray-tracing through a numerical model of the ionosphere. The ray-tracing computer program synthesizes the OIS ionograms by means of the “shooting method”. This method calculates trajectories of HF radio waves for different values of elevation angle and transmission frequency. There was a variety of calculated trajectories, from which we choose those which reach the receiver, and the selected paths provide a synthesis of the oblique ionograms. To simulate HF radio wave propagation, we apply a three-dimensional distribution of the electron density calculated with the mathematical model of the high-latitude ionosphere developed in the Polar Geophysical Institute (PGI). These numerical simulations permit us to interpret specific peculiarities of the OIS data such as abnormal propagation modes, increased delays of signals, enhanced MOF (maximum observed frequency) values etc. New results of the study are summarised as follows. (1) An unusual feature of the propagation along the path is the change of propagation mechanism during substorms on entering a path midpoint (or 1-hop reflection point) to the MIT. (2) Even weak substorms, having the distinguished intensities, lead to the appearance of different types of irregularities observed by the CUTLASS radar and therefore to the different propagation modes and F2MOF values. (3) The PGI model of the ionosphere was first used for ray-tracing at high latitudes. The model results are basically in a good qualitative agreement with experimental observations. This model provides the satisfactory agreement between the calculated and experimental F2MOF values while not correctly representing the fine structure of the experimental OIS ionograms at night. An agreement between the calculated and experimental data is better for day and evening hours than at night.     

电离层拟合方法及路径参数量化分析

针对多层电离层,采用反抛物层作为连接层来描述电离层之间区域的电子密度分布.以电子密度和其相对于高度的一阶偏导数连续为依据,拟合出连接层的临界频率、半厚度和底高等电离层参数.将电离层的准抛物模型和连接层的反抛物模型中计算出来的等离子频率代入Appleton-Hartree公式和Snell定理,结合射线的几何结构完成多层电离层中的射线追踪,同时计算群路径、相路径和传输距离这3种重要的路径参数.分析当给定频率时,路径参数与入射角之间的关系,群路径、相路径与传输距离之间的关系,频率对路径参数之间变化关系的影响;同时分析了给定入射角时,路径参数与射线工作频率之间关系和入射角对于路径参数之间变化关系的影响.     

Multi-beam sounding ionograms in the polar cap region: Absorption induced by proton precipitations

The simulation of the multi-beam ionograms in the polar cap region, assessing absorption effect is performed. It is reasonable to distinguish among four different mechanisms responsible for absorption: regular absorption due to solar UV illumination, absorption associated with energetic particles precipitation, absorption connected with X-rays flare and additional absorption in Auroral oval area. In this paper the absorption attributed to proton precipitations is envisaged. The computational model of the high-latitude ionosphere with irregularities oriented to application for the high frequency wave propagation problem was elaborated (Zaalov et al., 2005). A number of the quasi-vertical ionograms in the polar cap region were simulated on the basis of this model. A well-known algorithm (Sauer and Wilkinson, 2008) is applied for the absorption effects calculation. The simulated high-latitude ionograms with the absorption effect and the measured ionograms exhibit quite a good resemblance. This paper illustrates the importance of the understanding and taking into account the absorption effect in the presence of the various structural features in the polar ionosphere (in particular, patches of enhanced electron density) in interpreting ionosonde data.     

Numerical study of traveling ionosphere disturbances with vertical incidence data

High frequency ionosphere vertical sounding, as an important and representative application for detecting the ionosphere and studying the characteristics of radio propagation, can be utilized to monitor the ionosphere continuously variation and to acquire the ionosphere asymmetrical features of diverse scale above the ionosphere vertical sounding stations. This is a first article on real time application of numerical methods to obtain the parameters of traveling ionosphere disturbances (TIDs) using vertical incident ionograms. In this paper, the distribution of ionosphere electron density with TIDs is constructed using a background ionosphere model superimposed a perturbation theory model. The background ionosphere electron density is modelled by the inversion of vertical incident ionograms which are observed before the appearance of the disturbance. Based on the fourth order Adams-Bashforth-Moulton (the so-called ABM) predictor corrector method, instead of Runge-Kutta method, the fast digital ray tracing method is established. According to process of the disturbed trace simulation and parameters inversion, the characteristic parameters of ionosphere disturbance at different detection time can be obtained in real time. The numerical analysis of TIDs is then captured completely.     

时延-角度特性研究

本文在球不对称局域准抛物电离层模型下,设计了一套射线追踪算法的计算机程序。研究了时延随角度变化特性。指出了最小时延的传统求解方法的适应范围,讨论了时延-角度特性曲线与返回散射电离图的关系及其能量分布特点。      

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.     

A realistic simulation framework to evaluate ionospheric tomography

Observation of the 3-dimensional (3-D) electron density of the ionosphere is useful to study large-scale physical processes in space weather events. Ionospheric data assimilation and ionospheric tomography are methods that can create an image of the 3-D electron density distribution. While multiple techniques have been developed over the past 30 years, there are relatively few studies that show the accuracy of the algorithms. This paper outlines a novel simulation approach to test the quality of an ionospheric tomographic inversion. The approach uses observations from incoherent scatter radar (ISR) scans and extrapolates them spatially to create a realistic ionospheric representation. A set of total electron content (TEC) measurements can then be simulated using real geometries from satellites and ground receivers. This data set, for which the ‘truth’ ionosphere is known, is used as input for a tomographic inversion algorithm to estimate the spatial distribution of electron density. The reconstructed ionospheric maps are compared with the truth ionosphere to calculate the difference between the images and the truth.To demonstrate the effectiveness of this simulation framework, an inversion algorithm called MIDAS (Multi-Instrument Data Analysis Software) is evaluated for three geographic regions with differing receiver networks. The results show the importance of the distribution and density of GPS receivers and the use of a realistic prior conditioning of the vertical electron density profile. This paper demonstrates that when these requirements are met, MIDAS can reliably estimate the ionospheric electron density. When the region under study is well covered by GPS receivers, as in mainland Europe or North America, the errors in vertical total electron content (vTEC) are smaller than 1 TECu (2–4%) . In regions with fewer and more sparsely distributed receivers, the errors can be as high as 20–40%. This is caused by poor data coverage and poor spatial resolution of the reconstruction, which has an important effect on the calibration process of the algorithm.     

基于一维模型的高频电波加热电离层F层数值研究

基于电离层一维仿真加热模型,详细介绍了模型中电子的动量方程、连续性方程、能量方程和各类参量表达式,利用对角矩阵追赶法数值求解电离层F层加热过程,分析了不同时次电子数密度和电子温度的变化,讨论了不同频率和不同功率电波加热的情形.结果表明:当高频电波加热高电离层时,电子温度迅速增加,并很快趋近于饱和状态;电子密度的变化较为迟缓,但在加热过程中其变化幅度却越来越大;电子密度变化量在沿磁场方向上形成空洞和上下稠团两峰一谷构型;频率越大、功率越高的电波加热时,电子密度的变化也越大,但存在一适值频率的电波对电子温度的影响最小.      

N(h) profiles from ground and “Interkosmos-19” satellite data

During the years of high solar activity (1979–1980) simultaneous soundings of the ionosphere were made by the Bulgarian ground stations at Sofia and Michurin and by the satellite Interkosmos 19. The N(h) profiles for the ground stations were computed in association with IZMIRAN, Moscow and those for the satellite by SIB-IZMIRAN.During some of the satellite transit orbits the electron densities at the maximum of the F2 layer coincide precisely with those obtained from ground ionograms. In other cases differences were observed, perhaps because of the differences in time; the ionospheric stations at Sofia and Michurin worked on the ordinary 15-minute regime.     

New possibilities for mesospheric electricity diagnostics

A possible cause of large variations in the electron collision frequency could be the effect of strong external electric fields of atmospheric origin. This provides a new opportunity to take measurements of electric fields in the lower ionosphere using remote sensing instruments employing radio wave techniques. It has been proposed the technique for making estimates of strong mesospheric electric field intensities on the lower edge of the ionosphere using MF radar data on the effective electron collision frequency, and the data has been presented. The technique described permits a real-time derivation from MF radar data of changes in mesospheric electric field intensities, and estimates of electric current densities. Our results give proof that the source of strong mesospheric electricity is very likely to be a current source.     

Progress in ionospheric informatics based on electron-density profile analysis of ionograms

Evolution of techniques for electron-density profile analysis of ionograms attests to a close coupling with ionospheric modelling. Advanced achievements in both fields are used in the International Reference Ionosphere, IRI, the selection being made by criteria of internal consistency. This paper presents the results of comparisons of modern methods of calculating electron-density profiles from ground-based vertical-incidence ionograms. It is shown that the shapes of the resulting electron-density profiles depend upon the peak height and electron density of the E and F-layers. Formulae for the F2-layer sub-peak electron density, the E-F valley and their relation to the peak height of the ionosphere are presented as functions of solar-zenith angle. Implementation of such relationships allows refinement and considerable simplification of representations in the IRI, depending on the solar and geophysical controlling parameters.     

Daytime VHF amplitude scintillations recorded at an Indian low-latitude station,Waltair (17.7°N, 83.3°E) during 1997–2003

In this research, it is presented the daytime amplitude scintillations recorded at VHF frequency (244 MHz) at an Indian low-latitude station, Waltair (17.7°N, 83.3°E) during seven continuous years (1997–2003). Contrary to the nighttime scintillation seasonal trends, the occurrence of daytime scintillations maximizes during summer followed by winter and the equinox seasons. The fade depths, scintillation indices and the patch durations of daytime scintillations are meager when compared with their nighttime counterparts. A co-located digital high frequency (HF) ionosonde radar confirms the presence of sporadic (Es) layers when daytime scintillations are observed. The presence of daytime scintillations is evident when the critical frequency of the Es-layer (foEs) is ≥4 MHz and Es-layers are characterized by a highly diffuse range spread Es echoes as can be seen on ionograms. It is surmised that the gradient drift instability (GDI) seems to be the possible mechanism for the generation of these daytime scintillations. It is quite likely that the spread Es-F-layer coupling is done through polarization electric fields (Ep) that develop inside the destabilized patches of sporadic E layers, which are mapped up to the F region along the field lines as to initiate the daytime scintillations through the GDI mechanism. Further, the presence of additional stratification of ionosphere F-layer, popularly known as the F3-layer, is observed on ionograms once the Es-layers and daytime scintillations are ceased.     

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.     

Application of small-size antennas for estimation of angles of arrival of HF signals scattered by ionospheric irregularities

A modification of the Doppler Interferometry Technique is suggested to enable estimating angles of arrival of comparatively broadband HF signals scattered by random irregularities of the ionospheric plasma with the use of small-size weakly directional antennas. The technique is based on the measurements of cross-spectra phases of the probe radiation recorded at least in three spatially separated points. The developed algorithm has been used to investigate the angular and frequency-time characteristics of HF signals propagating at frequencies above the maximum usable one (MUF) for the direct radio path Moscow-Kharkiv. The received signal spectra show presence of three families of spatial components attributed, respectively, to scattering by plasma irregularities near the middle point of the radio path, ground backscatter signals and scattering of the sounding signals by the intense plasma turbulence associated with auroral activations. It has been shown that the regions responsible for the formation of the third family components are located well inside the auroral oval. The drift velocity and direction of the auroral ionosphere plasma have been determined. The obtained estimates are consistent with the classical conception of the ionospheric plasma convection at high latitudes and do not contradict the results of investigations of the auroral ionosphere dynamics using the SuperDARN network.     

Validation of FORMOSAT-3/COSMIC radio occultation electron density profiles by incoherent scatter radar data

Low Earth Orbiting satellites carrying a dual frequency GPS receiver onboard offer a unique opportunity to remote sensing of the global ionosphere on a continuous basis. No other profiling technique unifies profiling through the entire F2-layer with global coverage. The FORMOSAT-3/COSMIC data can make a positive impact on the global ionosphere study providing essential information about the height electron density distribution and particularly over regions that are not accessible with ground-based measuring instruments such as ionosondes and GPS dual frequency receivers. Therefore, it is important to verify occultation profiles with other techniques and to obtain experience in the reliability of their derivation. In the given study we present results of comparison of the electron density profiles derived from radio occultation measurements on-board FS-3/COSMIC and from the Kharkov incoherent scatter radar sounding.     

Method of the HF wave absorption evaluation based on GIRO network data

High frequency (HF) communication is strongly dependent on the state of the ionosphere, which specifies the mode structure of the radio wave propagating in ionosphere. Another core factor defining the strength of the HF signal at the receiving site is the ionospheric absorption. Accurate modelling the effect of absorption is an essential part of many studies of the HF propagation in the ionosphere.This paper proposes a method for estimating the absorption. The method is based on analysis of vertical sounding ionograms. The main idea of the approach is to compare the main parameters retrieved from measured and simulated ionograms. The combination of Global Ionospheric Radio Observatory (GIRO, http://giro.uml.edu) data and ionograms modelling allows for developing the empiric absorption model available at near real-time. The ionogram simulation taking into account absorption utilizes the NIM-RT (North Ionospheric Model and Ray Tracing) software. As a result, the proposed technique provides more reliable and accurate evaluation of minimum frequency at which echoes are observed in vertical incidence ionosonde soundings. The values of these frequencies should be used in the following simulation to optimize parameters in the empirical formulae for defining absorption HF wave in ionosphere.The ultimate objective of this work is the designing the method, which allows the simulating of HF radio channel accounting for regular absorption due to UV radiation of the Sun. Eventually it could be considered as some kind of the HF propagation forecasting.     

Ionosphere tomography based on spaceborne SAR

Two models of ionosphere tomography based on spaceborne SAR (Synthetic Aperture Radar) are proposed. For HF-SAR the signal with sweeping frequency lower than the characteristic frequency of ionosphere will be scatted during the ionosphere propagation and completely reflected at a corresponding height. The ionospheric electron density isolines looked as series of random surfaces can be reconstructed from the HF-SAR echoes by using the inverse scattering technique for layered rough surfaces and the method of moment (MoM). The numerical simulation show that due to the MoM can provide a full wave solution, the ionosphere tomography with high resolution can be obtained as long as enough sampling data of HF-SAR echoes are used. For VHF/UHF/P/L-band SAR the TEC (Total Electron Content) can be obtained from the SAR echoes scattered by some strong point targets (such as the calibrators, etc.) appeared in the SAR imaged ground region, and the ionosphere tomography can be performed by CT technique.     

昆明上空电离层D区电子密度季节变化的首次观测分析

通过分析2008年8月至2009年7月昆明站(25.6°N, 103.8°E) 中频(MF)雷达观测数据, 研究了太阳活动低年电离层D区电子密度的季节变化特性,发现D区电子密度主要呈现半年变化特征, 即在春秋季电子密度较大, 而在夏冬季则较小, 这与国际参考电离层(IRI)预测的年变化趋势不一致, 但与昆明站电离层测高仪的最低回波频率f_min的观测结果相符. 同时比较了D区电子密度半年变化与纬向风半年变化的关系, 发现二者之间保持了非常一致的变化趋势并对这种一致性的内在原因进行了分析.