Histogram-based ionogram displays and their application to autoscaling

A simple method is described for displaying and auto scaling the basic ionogram parameters foF2 and h’F2 as well as some additional layer parameters from digital ionograms. The technique employed is based on forming frequency and height histograms in each ionogram. This technique has now been applied specifically to ionograms produced by the IPS5D ionosonde developed and operated by the Australian Space Weather Service (SWS). The SWS ionograms are archived in a cleaned format and readily available from the SWS internet site. However, the method is applicable to any ionosonde which produces ionograms in a digital format at a useful signal-to-noise level. The most novel feature of the technique for autoscaling is its simplicity and the avoidance of the mathematical imaging and line fitting techniques often used. The program arose from the necessity to display many days of ionogram output to allow the location of specific types of ionospheric event such as ionospheric storms, travelling ionospheric disturbances and repetitive ionospheric height changes for further investigation and measurement. Examples and applications of the method are given including the removal of sporadic E and spread F.     

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.     

Ground and satellite-based observations of ionospheric plasma bubbles and blobs at 5.65° latitude in the Brazilian sector

This investigation uses simultaneous observations from all-sky imager system and an ionosonde collocated at Araguatins (5.65° S, 48.07° W and dip-latitude of 4.17° S), a near-equatorial region in Brazil. These simultaneous observations were used to investigate the occurrence of plasma bubbles and blobs in the field of the imaging system and their association with atypical range Spread-F signature in ionograms. Also, in-situ observation of plasma density from Swarm satellites were used to support the ground-based observations. Using a few cases, a methodology will be established to identify in the plasma blobs (atypical ESF) in the ionograms when there is the simultaneous observation of plasma bubbles and blobs in the field of view of the ionosonde. For this purpose, simultaneous sequence of OI 630.0 nm nightglow images and ionograms are presented for different case studies; 1. when there is the absence of a plasma bubble or blob, 2. when there is only the occurrence of plasma bubbles and 3. when there is the occurrence of plasma bubbles and blobs, in order to compare traces in the ionogram in all these case studies. With these we can cover all kinds of signatures in the ionograms corresponding to no irregularities, plasma bubbles only and plasma bubbles-blobs. These OI 630.0 nm nightglow and ionograms recorded simultaneously make it possible to establish a novel methodology to recognize in ionograms cases when there is the occurrence of Spread-F signature associated with bubble-blob in the FOV of the ionosonde.     

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.     

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.     

Practical method for routine analysis of the valley parameters between E- and F-region of the ionosphere

Measurement of the virtual height and the frequency of the minimum of the F-region extraordinary trace in digital and analog ionograms can provide a world-wide survey of the main parameters for the valley between E- and F-region ionization. The two added quantities establish also the starting point for the true height analysis of the F-region ionization.     

频高图Es层回波的自动识别

提出了一种基于决策树算法的自动识别Es层回波的方法以提高Es识别效率.通过自适应二值化以及中值滤波算法对频高图进行预处理，并提取有效的Es回波区域.利用Es层在虚高轴分布的规律，通过图像投影的方法选择出特征，将这些特征以及时间信息作为决策树算法的输入；通过人工识别方式为每幅频高图标明是否存在Es回波的标签，与决策树算法的输出进行比较.对训练集进行训练、剪枝，获取最接近人工识别结果的决策树，使用测试集对获取的决策树进行测试验证.本文使用云南普洱站（22.7°N，101.5°E）记录的频高图数据作为训练集，分别使用云南普洱站（22.7°N，101.5°E）和四川乐山站（29.5°N，103.7°E）的频高图数据作为测试集对该方法进行测试验证，结果表明该方法对普洱站和乐山站两地Es二跳回波的识别均具有较高的准确率，分别达到84.2%和82.8%.      

Reliability of electron density profiles

Deviations from horizontal stratification in the F-region can cause significant errors in electron density profile calculations from ionograms. Such situations exist every day during sunrise and sunset. Angle of arrival measurements and studies of the variation of other F-region parameters indicate that gravity waves are frequently strong enough to produce effects of comparable magnitudes. Ray tracing model studies permit a first order estimate of the resulting errors which are largest for the peak parameters.     

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.     

The accuracy of data from ionosondes for the estimation of hmF2 and the identification of global change in the ionosphere

A long temporal series of simulated ionograms was generated with a superimposed secular variation of −14 km/century on the h_mF2 parameter. These ionograms were interpreted by the automatic scaling program Autoscala. By applying four different empirical formulas, four artificial series of h_mF2 were generated and then processed with the same methods used by other authors for real data sets. Data analysis of the simulated ionograms revealed the artificially imposed long-term trend. These results lead to the conclusion, that regardless of the empirical formula used, the accuracy of h_mF2 from ionosonde measurements would be adequate to observe a long-term trend of −14 km/century.     

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.     

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.     

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.     

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.     

The variability and predictability of the IRI B0, B1 parameters over Grahamstown,South Africa

The International Reference Ionosphere (IRI) parameters B₀ and B₁ provide a representation of the thickness and shape, respectively, of the F2 layer of the bottomside ionosphere. These parameters can be derived from electron density profiles that are determined from vertical incidence ionograms. This paper aims to illustrate the variability of these parameters for a single mid latitude station and demonstrate the ability of the Neural Network (NN) modeling technique for developing a predictive model for these parameters. Grahamstown, South Africa (33.3°S, 26.5°E) was chosen as the mid latitude station used in this study and the B₀ and B₁ parameters for an 11 year period were determined from electron density profiles recorded at that station with a University of Massachusetts Lowell Center for Atmospheric Research (UMLCAR) Digisonde. A preliminary single station NN model was then developed using the Grahamstown data from 1996 to 2005 as a training database, and input parameters known to affect the behaviour of the F2 layer, such as day number, hour, solar and magnetic indices. An analysis of the diurnal, seasonal and solar variations of these parameters was undertaken for the years 2000, 2005 and 2006 using hourly monthly median values. Comparisons between the values derived from measured data and those predicted using the two available IRI-2001 methods (IRI tables and Gulyaeva, T. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7(6), 39–48, 1987.) and the newly developed NN model are also shown in this paper. The preliminary NN model showed that it is feasible to use the NN technique to develop a prediction tool for the IRI thickness and shape parameters and first results from this model reveal that for the mid latitude location used in this study the NN model provides a more accurate prediction than the current IRI model options.     

Implementation of a new characteristic parameter into the IRI sub-peak electron density profile

The shape of electron density profile in the International Reference Ionosphere could be improved significantly if the height h_g and electron density N_g of the F region sub-peak inflexion point were entered in the set of the profile standard parameters. To study variations of these important parameters, the N(h) analysis of the statistically-summarized ionograms at the latitudes of 40–80°N of the Eastern hemisphere has been carried out for the two-hours intervals of local time, three seasons (winter, equinox and summer) and two levels of solar activities characterized by Covington indicesF_10.7 = 100 and 200. It is shown that the parameters of the inflexion point can be expressed in most cases via the peak parameters of the F2 layer ash_g= 0.8 hmF2 and N_g= 0.5 NmF2.     

Bottomside N(h) profiles over Vietnam and their comparison with IRI

Over 70 ionograms were selected from the ground-based ionosonde in Hanoi, Vietnam at times when the AE-E satellite passed above the station. N(h) profiles were calculated from ionograms and thus compared with the IRI model and local AE-E ion densities. For the purpose of the presentation, 15 profiles were selected to cover all seasons and local times. The comparison shows that the observed daytime N(h) profiles have significantly higher gradients below the F2 peak density while at night the observed and the IRI profiles almost coincide. The difference is more pronounced in summer than in equinox. Wintertime comparisons are quite limited and do not lead to reliable conclusions. N_mF2 values taken from the CCIR programme are compared with those observed. The comparison shows that the deviations are not so large.     

A statistical study of the F2 layer stratification at the northern equatorial ionization anomaly

Ionograms recorded at Puer station (PUR, 22.7°N, 101.05°E, Dip Latitude 12.9°N) in the Southwest of China from January 2015 to December 2016 were used to study characteristics of the F2 layer stratification at the northern equatorial ionization anomaly. Ionosonde observations show that the development of the F2 layer stratification is different under different conditions. Both the upward and downward movement of the F2 layer stratification could be observed. The F2 layer stratification could occur both at daytime and nighttime. The new cusp could originate from different positions on ionograms. Moreover, statistical results indicate that the F2 layer stratification occurred later in the winter than in other seasons at daytime, it occurred frequently in the local spring, and most of ionograms with the F2 layer stratification at post-midnight occurred in March and April. Our results also show that the F2 layer stratification has a correlation with solar activity.     

Controlling role of maximum usable frequencies in ionospheric informatics

Considering peak parameters of the ionospheric F region which might be observed/modelled/predicted a few relations connecting such data are pointed out. It is pointed out that the sub-peak semi-thickness of the profile depends on the peak height of the F2-layer. The dispersion range of monthly measurement of the MUFs (maximum usable frequency) inferred from M(3000)F2 is compared with that of the critical frequencies foF2 and found to be larger. The inverted Shimazaki formula yields straight-forward values of MUF from the F2 layer critical frequency and the peak height ZmF2. When comparing MUF-3000 values scaled routinely from ionograms with those obtained from foF2 and ZmF2 of the profile analysis, good agreement is obtained. It is felt that the parameter MUF-3000 is a valuable means for checking data obtained by measurements or by models or by predictions.