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.     

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.     

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.     

High-resolution, coupled thermosphere–ionosphere models for space weather applications

Current first-principles global models of the coupled thermosphere–ionosphere (T–I) system use grids that are too coarse to simulate the mesoscale and small-scale structures that occur in this complex system. These small-scale and mesoscale structures have a great effect on global-scale neutral and plasma distributions and have important consequences for daily space weather. In this paper, we present a new first-principles, high-resolution, T–I nested grid (TING) polar cap model that incorporates multiple nesting levels and two-way interaction. The TING model simulation of the electron densities and temperatures demonstrates the importance of high spatial resolution. It is found that both the mid-latitude electron density trough and its associated dawn electron temperature peak are more pronounced and structured in the nested grid than in the coarse grid. Using the TING model to simulate ionospheric F₂ region electron density variations with geomagnetic activity and universal time (UT) is also discussed.     

Physical implications of discrepancy between summer and winter PC indices observed in the course of magnetospheric substorms

The PC index based on a statistically justified relationship between the polar cap magnetic activity and the interplanetary electric field E_KL has been derived as a value standardized for the E_KL intensity regardless of season, UT and hemisphere. As a result, the summer and winter PC indices are consistent with one another under ordinary conditions. Discrepancies between the summer and winter PC indices arising in the course of magnetospheric substorms are analyzed in this paper. It is argued that the channel of enhanced conductivity, formed in the auroral oval owing to intense auroral particle precipitation, strongly improves the conditions for closure of the Region 1 field-aligned currents in the winter dark polar region but only trivially affects the conditions of the Region 1 FAC closure in the summer sunlit ionosphere. Since the coefficients describing the relationship between E_KL and the polar cap magnetic activity were derived for statistically justified (i.e., mean) conditions, their application to such abnormal situation, as intense field-aligned currents in the winter dark polar region, leads to overestimation of the winter PC index. The summer and winter PC indices level off as soon as the intense auroral particle precipitation terminates and the auroral ionosphere in the winter and summer polar caps returns to the ordinary (statistically justified) state.     

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.     

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.     

Estimation of ground range on the sweep frequency backscatter leading edge

The high frequency management system with backscatter radar supplies the real time ionosphere channel conditions to high frequency users, which leads to the demand for the ground range between the radar location and the scatters on the distant ground. The ionosphere electron density profile is usually inversed to obtain the ground range. An inversion algorithm, with which the ground range on the leading edge of the backscatter ionograms can be obtained without electron density, is presented in this paper. The ray path geometry of the backscatter sounding and the change in the group path on the leading edge with operating frequency are used to derive the ground range. Synthesized backscatter ionogram and experimental backscatter ionograms are processed to validate the algorithm. The results indicate that the algorithm is usable for high frequency management system.     

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.     

On high-altitude electric fields in the auroral downward current region and their coupling to ionospheric electric fields

The downward field-aligned current region plays an active role in magnetosphere–ionosphere coupling processes associated with aurora. A quasi-static electric field structure with a downward parallel electric field forms at altitudes between 800 km and 5000 km, accelerating ionospheric electrons upward, away from the auroral ionosphere. Other phenomena including energetic ion conics, electron solitary waves, low-frequency wave activity, and plasma density cavities occur in this region, which also acts as a source region for VLF saucers. Results are presented from high-altitude Cluster observations with particular emphasis on the characteristics and dynamics of quasi-static electric field structures. These, extending up to altitudes of at least 4–5 Earth radii, appear commonly as monopolar or bipolar electric fields. The former occur at sharp boundaries, such as the polar cap boundary whereas the bipolar fields occur at softer boundaries within the plasma sheet. The temporal evolution of quasi-static electric field structures, as captured by the pearls-on-a-string configuration of the Cluster spacecraft, indicates that the formation of electric field structures and of ionospheric plasma density cavities are closely coupled processes. A related feature of the downward current is a broadening of the current sheet with time, possibly related to the depletion process. Preliminary studies of the coupling of electric fields in the downward current region, show that small-scale structures are typically decoupled from the ionosphere, similar to what has been found for the upward current region. However, exceptions are also found where small-scale electric fields couple perfectly between the ionosphere and Cluster altitudes. Recent FAST results indicate that the degree of coupling differs between sheet-like and curved structures, and that it is typically partial. The electric field coupling further depends on the current–voltage relationship, which is highly non-linear in the downward current region, and still unrevealed, as to its specific form.     

A comparative study of decision tree,random forest,and convolutional neural network for spread-F identification

Ionospheric spread-F (SF) is a commonly observed phenomenon of electron density perturbation in the F-layer. The ionospheric irregularities structure has an adverse effect on the propagation of electromagnetic waves in the ionosphere. The automatic identification of ionospheric spread-F and statistical study of the formation of spread-F are of great significance to the study of the physical mechanism of ionospheric inhomogeneity and for prediction of ionospheric irregularities. In this paper, we describe and implement three automatic identification methods of spread-F based on machine learning: decision tree, random forest, and convolutional neural network (CNN). The performance of these automatic identification methods was verified using a large set of test data. Results show that the accuracy of all three methods on identifying ionograms with spread-F exceeded 90%. After comparing the results of the three methods, we found that the decision tree method was the simplest and with the structure easiest to be understood, and it required the shortest interpretation time. In terms of the identification results, the random forest method provided better results than the decision tree method, and the CNN method was the best at accurately identifying ionograms with spread-F.     

Use of varying shell heights derived from ionosonde data in calculating vertical total electron content (TEC) using GPS – New method

The dispersive nature of the ionosphere makes it possible to measure its total electron content (TEC). Thus Global Positioning System, which uses dual-frequency radio signals, is an ideal system to measure TEC. When data from an ionosonde situated in polar region was observed, the height of an approximated thin shell of electrons (shell height) used in GPS studies was seen not to be fixed but rather changing with time. Here we introduce a new method in which we included the varying shell heights derived from the ionosonde to map the slant total electron content from GPS to obtain a more precise vertical total electron content of the ionosphere contrary to some previous methods which used fixed shell heights. In this paper we also compared the ionosonde derived TEC with the GPS derived vertical TEC (vTEC) values. These GPS vTEC values were obtained from GPS slant TEC (sTEC) measurements using both fixed shell height and varying shell heights (from ionosonde measurements). For the polar regions, the varying shell height approach produced better results than the fixed shell height and compared to exponential function, Chapman function seems to be a better function to model the topside ionosphere.     

一种数字电离层测高仪系统中O波与X波的分离方法

电离层测高仪系统中,O波与X波的分离非常重要.本文提出一种新的O波与X波分离的实现方法,根据电离层回波的极化特性,通过在接收电路上采用数字方法合成圆极化波的方式,实现了对O波与X波的有效分离.与现有DPS-4测高仪系统采用模拟域电信号合成圆极化波的方法相比,本文方法通过在数字下变频处理过程中引入±90°相移,消除了在模拟域电信号合成方法中相移器的带宽限制和非线性问题.另外,该方法在实现发射和接收信号的极化状态转换时不需要连续切换多个模拟开关,从而提高了整个系统的稳定性.      

太阳活动区R9077引起的强烈吸收事件

利用南极中山站的观测数据分析了太阳活动区R9077所引起的强烈吸收事件，其中2000年7月14日的太阳质子事件引起了持续3天多强烈的极盖吸收，同时，激烈的太阳活动使磁层处于极度扰动状态，磁层高能粒子沉降使许多持续时间较短的吸收峰叠加在极盖吸收背景之上，最突出的是7月1日的吸收增强事件，其最大值达26dB。这是自1997年2月中山站安装成像式宇宙噪声接收机以来观测到的最强的吸收，另一个较突出的吸收峰发生在14日1753UT前后，本文还讨论了产生这些吸收的原由。     

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.     

South polar permanent CO2 ice cap presentation in the Global Mars Multiscale Model

The atmospheric influence caused by the Martian permanent south CO₂ ice cap is examined to improve the Global Mars Multiscale Model (GM3) to see if it can significantly improve the representation of south polar meteorology. However, the seasonal carbon dioxide ice in the polar regions is presented in the surface ice simulation by the Global Mars Multiscale Model but the model does not produce a permanent south CO₂ ice cap, and the physics code must modify to capture the realistic physical such as ice process detail; probably makes a bias in terms of total CO₂ ice and meteorological processes in the model aside from ice formation. The permanent south CO₂ ice cap in the model can significantly improve the representation of south polar meteorology for example in predicted surface temperatures, surface pressures, horizontal and zonal winds over the south cap and possible initiation of dust storms at south polar region during the southern summer period.     

Theoretical aspects of ionosphere-magnetosphere-solar wind coupling

We initially review magnetospheric plasma distributions emphasising the explanative power of the open model of the magnetosphere. Next we turn attention to the role of MHD wave phenomena both as transients and as standing structures in the coupled flow system. Finally the physics of the magnetopause flux transfer events is discussed. The logic of the existing interpretation suggests a miniature twin vortex convection system should be created in part of the polar cap ionosphere magnetically connecting to the magnetopause flux tube.     

短波射线追踪技术中的电离层混合建模方法

在国际参考电离层模型和多层准抛物模型的基础上,提出了一种混合应用两种模型进行电离层建模的新方法.利用射线追踪技术,分别对混合模型和传统国际参考电离层模型下短波射线在电离层中的轨迹进行了仿真,得到了电波群路径.通过与实测电波群路径的对比,结果表明:对中国中纬度地区在电离层混合模型下的射线追踪精度优于传统国际参考电离层模型下的射线追踪精度,同时混合建模方法降低了多层准抛物模型对输入条件的要求,扩展了多层准抛物模型在射线追踪技术中的应用范围.     

Magnetospheric electric field variations caused by storm-time shock fronts

On January 20, 2005 there was an X 7.1 solar flare at 0636 UT with an accompanied halo coronal mass ejection (CME). The resultant interplanetary shock impacted earth ∼36 h later. Near earth, the Advanced Composition Explorer (ACE) spacecraft observed two impulses with a staircase structure in density and pressure. The estimated earth-arrival times of these impulses were 1713 UT and 1845 UT on January 21, 2005. Three MINIature Spectrometer (MINIS) balloons were aloft on January 21st; one in the northern polar stratosphere and two in the southern polar stratosphere. MeV relativistic electron precipitation (REP) observed by all three balloons is coincident (<3 min) with the impulse arrivals and magnetospheric compression observed by both GOES 10 and 12. Balloon electric field data from the southern hemisphere show no signs of the impulse electric field directly reaching the ionosphere. Enhancement of the balloon-observed convection electric field by as much as 40 mV/m in less than 20 min during this time period is consistent with typical substorm growth. Precipitation-induced ionospheric conductivity enhancements are suggested to be (a) the result of both shock arrival and substorm activity and (b) the cause of rapid (<6 min) decreases in the observed electric field (by as much as 40 mV/m). There is poor agreement between peak cross polar cap potential in the northern hemisphere calculated from Super Dual Auroral Radar Network (SuperDARN) echoes and horizontal electric field at the MINIS balloon locations in the southern hemisphere. Possible reasons for this poor agreement include (a) a true lack of north–south conjugacy between measurement sites, (b) an invalid comparison between global (SuperDARN radar) and local (MINIS balloon) measurements and/or (c) radar absorption resulting from precipitation-induced D-region ionosphere density enhancements.     

南极GPS定位中的电离层延迟分析

基于南极地区国际GNSS服务组织(IGS,International GNSS Service)跟踪站的全球定位系统(GPS,Global Position System)双频实测数据,分析了南极地区电离层延迟的变化情况及其二阶项延迟对南极GPS定位结果的影响.结果表明:南极地区的总电子含量(TEC,Total Electron Content)日间波动频繁,其日间TEC最大值变化较中纬度地区剧烈;在南极地区夏季,电离层二阶项延迟对GPS定位结果的影响可达cm级.同时,由于欧洲定轨中心(CODE,The European Center for Orbit Determination)提供的全球电离层模型(GIM,Global Ionosphere Maps)在南极区域应用的局限性,通过选取南极地区6个IGS跟踪站作为基准站建立了区域电离层TEC模型(RIM,Regional Ionosphere Model).经实测数据计算证明,对于南极地区,RIM的定位精度在一定程度上优于全球电离层模型GIM.