Effects of a “day-time” substorm on the ionosphere and radio propagation

Propagation mechanisms of lateral (non-great-circle) signals on a high-latitude HF radio path during magnetospheric substorms that occurred in the day-time have been considered. The path is equipped with oblique ionospheric sounding (OIS) from Murmansk to St. Petersburg. The OIS method gives the possibility to determine propagation modes, MOF (maximum observed frequency) values, signal delays, etc. Data of the CUTLASS radar, the IMAGE magnetometer system, the Finnish riometer chain, and the Tromso ionosonde were also used for the analysis. The main results are the following: (1) the lateral signal propagation takes place, as a rule, if the path midpoint is located near the irregularity region that moves sharply from high to low latitudes. The lateral signal propagation appearing during day-time is a new effect. (2) Formation of dense field-aligned irregularities during a substorm leads to decreasing F2MOF values on radio paths. These results can be useful for problems of radiolocation, HF communications and navigation.     

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.     

电离层Es对HF和VHF传播信号的影响

Es层是存在于电离层中的电子密度非常高的偶发E层（Sporadic E）,其电子密度可达常规E层的100倍.电离层Es能够反射原本穿透F层的VHF低频段（30~150MHz）无线电波,而且对HF高频段（10~30MHz）无线电波传播具有显著影响.运用垂测和斜测观测数据,研究HF频段Es层电波传播特征,得到了不同类型及不同高度Es层的衰减系数.根据f₀E_s的日变化规律,可得HF频段衰减系数的日变化规律,进而分析并得到Es层对短波传播的影响.不存在电离层Es时,通常无法通过电离层实现VHF超远距离通信.为了对VHF链路通过电离层Es的传播衰减进行定量分析,根据EBU多条链路的观测结果,拟合并建立了电离层Es衰减模型.将该模型、ITU模型和观测数据进行对比,发现本文建立的模型准确度更高.利用建立的模型,对电离层Es不同临频f₀E_s条件下接收信号场强和电压随传播距离的变化进行了计算,结果可为VHF链路设计及建立提供参考.      

Fading in the HF ionospheric channel and the role of irregularities

It is well known that the ionosphere affects radio wave propagation especially in the high frequency (HF) range. HF radio waves reflected by the ionosphere can reach considerable distances, often with changes in amplitude, phase, and frequency. The ionosphere is a dispersive in frequency and time, bi-refractive, absorbing medium, in which multipath propagation due to traveling irregularities is very frequent. The traveling irregularities undulate the reflecting ionospheric layer, introducing variations in signal amplitude (fading). In this multipath time variant channel fading is mainly considered, even though it is not the sole effect. Echo signals from a single reflection, as in ionospheric vertical sounding (VIS) techniques, are affected by a certain degree of variability even in quiet ionospheric conditions. In this work the behavior of the ionospheric channel is studied and characterized by observing the power variation of received echoes using the VIS technique. Multipath fading was analyzed quantifying the power variation of the signal echo due to irregularities on a temporal scale from 0.5 to 256 s. An experimental set-up derived from an ionosonde was implemented and the analysis was performed employing a special numerical algorithm operating off-line on the acquired time sequence of the signal. The gain-loss of the irregularity shapes are determined in some special cases.     

Impact of geomagnetic storm of September 7–8, 2017 on ionosphere and HF propagation: A multi-instrument study

The present study reveals the features of ionospheric parameters variations during the geomagnetic storm of September 7–8, 2017. In particular, parameters of vertical (foF2, foEs) and oblique ionospheric sounding (MOF, modes), absorption level, Total Electron Content (TEC) and particle fluxes at high altitudes were under analysis. The storm was characterized by two Dst-index mimima and can be considered as a sequence of two storms: first - with Dstmin?=??142?nT at 02 UT on September 8th and second - with Dstmin?=??122?nT and at 15 UT on September 8th. It was found that these two storms had different impacts on the ionosphere and HF propagation at mid- and high-latitudes of Northern Hemisphere. The signals of vertical and oblique ionospheric sounding were present in all ionograms before the first storm. Further, at the maximum of the first storm these signals were totally absorbed. Then, before the second storm and during its maximum the signals were detected again in the ionograms due to the low absorption. GOES satellite data showed the significant burst of electrons and protons only during the first storm and small particle fluxes - during the second storm. This feature was also confirmed with GPS data: TEC increased during the first storm and decreased during the second storm.     

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.     

Study of the ionosphere during a magnetospheric storm using data on the ground and in space

For the magnetospheric storm of May 14–16, 1997 geophysical data of satellites DMSP and IMP-8 are compared with data of radio propagation on the high-latitude HF radio path of Heiss Island – St. Petersburg and data from European ionosondes. Peculiarities of variations of the operational frequencies range MOF–LOF (maximum and lowest observed frequencies) on the path were considered. The range has been determined by the method of oblique ionospheric sounding (OIS). The latter is more informative for observations during a magnetic storm compared to the vertical sounding method. Nevertheless, an analysis of variations of the critical frequency of the ionospheric F2 layer from the chain of European ionosondes was carried out. For interpretation of results, data of magnetospheric parameters, AE-indexes and riometer data were used. The variations of both frequency range on the path and critical frequencies of the F2 layer through the ionosondes chain during the disturbed period had certain regularities of behaviour. These regularities are being explained from the physical point of view. The analysis of the satellite DMSP data has showed that a magnetospheric disturbance causes displacement equatorward of precipitation and some growth of its width and energy.     

Anthropogenic trigger of substorms and energetic particles precipitations

The high-frequency (HF) emission in near-Earth space from various powerful transmitters (radio communications, radars, broadcasting, universal time and navigation stations, etc.) form an integral part of the modern world that it cannot do without. In particular, special-purpose research facilities equipped with powerful HF transmitters are used successfully for plasma experiments and local modification of the ionosphere. In this work, we are using the results of a complex space-ground experiment to show that exposure of the subauroral region to HF emission can not only cause local changes in the ionosphere, but can also trigger processes in the magnetosphere–ionosphere system that result in intensive substorm activity (precipitations of high-energy particles, aurorae, significant variations in the ionospheric parameters and, as a consequence, in radio propagation conditions).     

Studying characteristics of traveling ionospheric disturbances using U-shaped traces on vertical incidence ionograms

This paper presents a case study when due to the descending additional U-shaped trace on vertical incidence ionograms, increased critical frequency stabilizes. This corresponds to an ionospheric disturbance that moves toward the ionosonde and then stays overhead.Within a 2D model, traveling ionospheric disturbances (TIDs) are superimposed on the inverted background ionosphere. So ray tracing is used to obtain propagation paths through non-stratified ionosphere thus synthesizing the disturbed ionogram traces. Investigated are changes in the cusp shape caused by varying TID parameters. A cusp-fitting method to determine the TID amplitude, spatial scale, and horizontal drift velocity are shown.     

Lower ionosphere response to external forcing: A brief review

There are two ways of external forcing of the lower ionosphere, the region below an altitude of about 100 km: (1) From above, which is directly or indirectly of solar origin. (2) From below, which is directly or indirectly of atmospheric origin. The external forcing of solar origin consists of two general factors – solar ionizing radiation variability and space weather. The solar ionization variability consist mainly from the 11-year solar cycle, the 27-day solar rotation and solar flares, strong flares being very important phenomenon in the daytime lower ionosphere due to the enormous increase of the solar X-ray flux resulting in temporal terminating of MF and partly LF and HF radio wave propagation due to heavy absorption of radio waves. Monitoring of the sudden ionospheric disturbances (SIDs – effects of solar flares in the lower ionosphere) served in the past as an important tool of monitoring the solar activity and its impacts on the ionosphere. Space weather effects on the lower ionosphere consist of many different but often inter-related phenomena, which govern the lower ionosphere variability at high latitudes, particularly at night. The most important space weather phenomenon for the lower ionosphere is strong geomagnetic storms, which affect substantially both the high- and mid-latitude lower ionosphere. As for forcing from below, it is caused mainly by waves in the neutral atmosphere, i.e. planetary, tidal, gravity and infrasonic waves. The most important and most studied waves are planetary and gravity waves. Another channel of the troposphere coupling to the lower ionosphere is through lightning-related processes leading to sprites, blue jets etc. and their ionospheric counterparts. These phenomena occur on very short time scales. The external forcing of the lower ionosphere has observationally been studied using predominantly ground-based methods exploiting in various ways the radio wave propagation, and by sporadic rocket soundings. All the above phenomena are briefly mentioned and some of them are treated in more detail.     

大功率无线电波与低电离层的相互作用

地面入射的大功率无线电波能加热电离层等离子体，引起电离层电子温度和密度的扰动，实现电离层的地面人工变态．本文中，着重考虑电波和电离层相互作用过程中的自吸收，构造一个自治的相互作用模型，在一定功率和频率的加热电波作用下，利用该模型计算了白天低电离层电子温度和由温度的变化而引起的电子密度的变化．计算结果表明，在白天低电离层，电波的自吸收在90km以下比较显著，而最大温度变化在70km高度上，大约增加了2倍．在α复合的假设下，电子密度变化幅度随高度的增加而减少，在70km处，大约增加55％、120km处则为4％左右．     

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.     

Impact of heliogeophysical disturbances on ionospheric HF channels

The article presents the results of the observation of a strong magnetic storm and two X-ray flares during the summer solstice in 2015, and their impact on the HF signals characteristics in ionospheric oblique sounding. It was found that the negative phase of the magnetic storm led to a strong degradation of the ionospheric channel, ultimately causing a long blackout on paths adjacent to subauroral latitudes. On mid-latitude paths, the decrease in 1FMOF reached ～50% relative to the average values for the quiet ionosphere. It is shown that the propagation conditions via the sporadic Es layer during the magnetic storm on a subauroral path are substantially better than those for F-mode propagation via the upper ionosphere. The delay of the sharp decrease in 1FMOF during the main phase of the magnetic storm allowed us to determine the propagation velocity of the negative phase disturbances (～100 m/s) from subauroral to mid-latitude ionosphere along two paths: Lovozero – Yoshkar-Ola and Cyprus – Nizhny Novgorod. It is shown that both the LOF and the signal/noise ratio averaged over the frequency band corresponding to the propagation mode via the sporadic Es layer correlate well with the auroral AE index. Using an over-the-horizon chirp radar with a bistatic configuration on the Cyprus – Rostov-on-Don path, we located small-scale scattering irregularities responsible for abnormal signals in the region of the equatorial boundary of the auroral oval.     

Supertyphoon Hagibis action in the ionosphere on 6–13 October 2019: Results from multi-frequency multiple path sounding at oblique incidence

The purpose of the present paper is to describe the observations of the variations in the parameters of HF radio waves propagating through the ionosphere when the action of the super typhoon Hagibis on 6–13 October 2019 occurred. The observations have been made with the Harbin Engineering University (the People's Republic of China) multi-frequency multiple path radio system involving the software-defined technology. The action of the super typhoon has been shown to be accompanied by enhanced atmospheric wave activity acting to generate wave processes with periods of 10 to 120 min. Coupling in the atmosphere–upper-atmosphere–ionosphere system has been confirmed to be carried out with atmospheric gravity waves. The ionosphere underwent the greatest impact on those days when the supertyphoon had maximum energy, on 8, 10, and especially 9 October 2019, and when it was found to be in an ～2,500–3,000-km distance range from the propagation path midpoints. Under the action of wave processes, the height of the reflection region was observed to oscillate within the ±(30–90 km) limits. The amplitude of the quasi-periodic variations in the ionospheric F-region electron density was estimated to be 10–12% for periods of ～20 min, and 30–60% for periods of ～60–120 min. The joint action of the dusk terminator and the supertyphoon has been confirmed to enhance wave activity in the ionosphere. Similar effects for the dawn terminator have not been detected.     

电离层中小尺度扰动的多频多普勒快速探测与分析

电离层反射回波的高频多普勒频移测量是研究电离层扰动的重要方法之一.利用基于单频信号相位变化率测量的实时探测方法,获取连续高精度多频多普勒频高图,以实现中小尺度电离层扰动的快速探测研究.该方法被应用于敏捷式高性能电离层无线电多频探测系统样机平台.经验证在该平台上获取的多普勒频移分辨率可达0.039Hz,频高图探测周期最短小于1min.基于这种快速探测方法和平台在武汉观测站进行了较长时间的观测和数据采集,获得高精度多普勒频移并反演出电子浓度等值面法向运动速度,得到电离层反射寻常波的多普勒频高图和垂直扰动速度等信息,进而推演出电离层扰动随时间和空间的实时变化特征.对这些时域信号进行频谱分析,初步结果显示这些扰动主要是由极区活动激发的中国中部地区冬季出现率较高的中尺度TID.此外,对三种常用的电离层扰动反演分析方法进行了对比研究,结果显示电离层扰动的变化趋势基本一致,说明观测数据和研究方法可靠稳定,为多频多点电离层扰动的传播特性观测研究提供了基础.      

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.     

A method for automatic detection of equatorial spread-F in ionograms

A method is presented for automatic detection of spread-F. The method is based on an image recognition technique and is applied to ionograms recorded at the ionospheric station of Tucumán (26.9°S, 294.6°E). The performance achieved is statistically evaluated and demonstrated with significant examples. The proposed method improves Autoscala's ability to reject ionograms with insufficient information, including those featuring Spread-F. Automatic identification of cases of spread-F is of additional interest in Space Weather applications, when it helps detect degraded radio propagation conditions.The present data analysis is a retrospective study but forms the basis for real-time application as an extension of Autoscala’s capabilities.     

大功率高频无线电波对赤道E区双流不稳定性的影响--PIC静电粒子模拟研究

地面入射的大功率高频无线电波(泵波)和电离层等离子体之间的参数相互作用，能够引起静电波的激发，在一定条件下，产生不稳定性．本文用PIC静电粒子模拟方法，研究泵波与赤道电离层E区等离子体的相互作用．研究结果表明，泵波能够控制双流不稳定性的发生，在不同条件下，泵波对双流不稳定性起着稳定与不稳作用，模拟结果定性地与理论研究结果相符合，这为我们对不规则体产生的地面人工控制提供了依据．     

Updating of the IRI-1990 with respect to substorms

The method of updating the IRI-1990 model for the effects of large scale travelling ionospheric disturbances (LSTIDs) of the F-region have been established as a result of the statistical analysis of VI sounding data of European- Asian sector for 380 substorms. The effect of a LSTID occurring during a substorm is modelled as an additive correction to the quiet IRI foF2 and hmF2 values which depends on the time since the onset of the substorm and on the maximum AE-index value reached in it. In case of a sequence of substorms the additive corrections ΔfoF2 and ΔhmF2 are taking into account. The method for night-time ionosphere is described in detail, the listing of the FORTRAN-77 computer program is available on the floppy disk. The method is workable as correction to any ionospheric model.     

The COMPLEIK subroutine of the IONORT-ISP system for calculating the non-deviative absorption: A comparison with the ICEPAC formula

The present paper proposes to discuss the ionospheric absorption, assuming a quasi-flat layered ionospheric medium, with small horizontal gradients. A recent complex eikonal model (Settimi et al., 2013b) is applied, useful to calculate the absorption due to the ionospheric D-layer, which can be approximately characterized by a linearized analytical profile of complex refractive index, covering a short range of heights between h₁ = 50 km and h₂ = 90 km. Moreover, Settimi et al. (2013c) have already compared the complex eikonal model for the D-layer with the analytical Chapman’s profile of ionospheric electron density; the corresponding absorption coefficient is more accurate than Rawer’s theory (1976) in the range of middle critical frequencies. Finally, in this paper, the simple complex eikonal equations, in quasi-longitudinal (QL) approximation, for calculating the non-deviative absorption coefficient due to the propagation across the D-layer are encoded into a so called COMPLEIK (COMPLex EIKonal) subroutine of the IONORT (IONOspheric Ray-Tracing) program ( Azzarone et al., 2012). The IONORT program, which simulates the three-dimensional (3-D) ray-tracing for high frequencies (HF) waves in the ionosphere, runs on the assimilative ISP (IRI-SIRMUP-P) discrete model over the Mediterranean area ( Pezzopane et al., 2011). As main outcome of the paper, the simple COMPLEIK algorithm is compared to the more elaborate semi-empirical ICEPAC formula (Stewart, undated), which refers to various phenomenological parameters such as the critical frequency of E-layer. COMPLEIK is reliable just like the ICEPAC, with the advantage of being implemented more directly. Indeed, the complex eikonal model depends just on some parameters of the electron density profile, which are numerically calculable, such as the maximum height.