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.     

A comparative study of ionospheric spread-F and scintillation at low- and mid-latitudes in China during the 24th solar cycle

The spread-F echo of ionograms and scintillation of satellite signal propagation along the Earth-space path are two typical phenomena induced by ionospheric irregularities. In this study, we obtained spread-F data from HF (high frequency) digital ionosonde and scintillation index (S4) data from L-band and UHF receivers at low- and mid-latitudes in China during the 24th solar cycle. These four sites were located at Haikou (HK) (20°N, 110.34°E), Kunming (KM) (25.64°N, 103.72°E), Qingdao (QD) (36.24°N, 120.42°E), and Manzhouli (MZL) (49.56°N, 117.52°E). We used these data to investigate spread-F and scintillation occurrence percentages and variations with local time, season, latitude and solar activity. A comparative study of spread-F and scintillation occurrence rates has been made. The main conclusions are as follows: (a) FSF occurred mostly during post-midnight, while RSF and scintillation appeared mainly during pre-midnight at HK and KM; (b) FSF occurrence rates were larger at QD and MZL than expected; (c) the FSF occurrence percentages were anti-correlated with solar activity at HK and KM; meanwhile RSF and scintillation occurrence rates increased with the increase of solar activity at this two sites; (d) the highest FSF occurrence rates mostly appeared during the summer months, while RSF and scintillation occurred mostly in the equinoctial months at HK and KM; (e) the scintillation occurrence was usually associated with the appearance of RSF, probably due to a different physical mechanism comparing with FSF. Some of these results verified the conclusions of previous papers, whereas some show slight difference. These results are important in understanding ionospheric irregularities variations characteristic at low- and mid-latitudes in China.     

频高图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%.      

Spread-F occurrence during geomagnetic storms near the southern crest of the EIA in Argentina

This work presents, for the first time, the analysis of the occurrence of ionospheric irregularities during geomagnetic storms at Tucumán, Argentina, a low latitude station in the Southern American longitudinal sector (26.9°S, 294.6°E; magnetic latitude 15.5°S) near the southern crest of the equatorial ionization anomaly (EIA). Three geomagnetic storms occurred on May 27, 2017 (a month of low occurrence rates of spread-F), October 12, 2016 (a month of transition from low to high occurrence rates of spread-F) and November 7, 2017 (a month of high occurrence rates of spread-F) are analyzed using Global Positioning System (GPS) receivers and ionosondes. The rate of change of total electron content (TEC) Index (ROTI), GPS Ionospheric L-band scintillation, the virtual height of the F-layer bottom side (h'F) and the critical frequency of the F2 layer (foF2) are considered. Furthermore, each ionogram is manually examined for the presence of spread-F signatures.The results show that, for the three events studied, geomagnetic activity creates favorable conditions for the initiation of ionospheric irregularities, manifested by ionogram spread-F and TEC fluctuation. Post-midnight irregularities may have occurred due to the presence of eastward disturbance dynamo electric fields (DDEF). For the May storm, an eastward over-shielding prompt penetration electric field, (PPEF) is also acting. A possibility is that the PPEF is added to the DDEF and produces the uplifting of the F region that helps trigger the irregularities. Finally, during October and November, strong GPS L band scintillation is observed associated with strong range spread-F (SSF), that is, irregularities extending from the bottom-side to the topside of the F region.     

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.     

In situ studies of electron density during equatorial spread-F

A Langmuir probe operated in fixed bias mode was launched onboard a RH-560 rocket from the Sriharikota Range (SHAR, Lat.13° 42'N, Geog. Long.80° 14'E, dip 10°) India on October 1, 1980 at 21h03 IST, to study the electron density profile and the electron density irregularities in the equatorial spread-F. The payload was designed to study medium and large scale irregularities. A highly variable and structured electron density profile was obtained. This was the first rocket launch in the Indian zone during spread-F condition.     

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.     

Morphological studies on ionospheric VHF scintillations over an Indian low latitude station during a solar cycle period (2001–2010)

The amplitude scintillations data recorded at 244 MHz from the geostationary satellite, FLEETSAT (73°E) at a low latitude station, Waltair (17.7°N, 83.3°E) during the ten year period of high to low solar activity from 2001 to 2010 is considered to study the occurrence characteristics of the VHF scintillations. A close association between the intense scintillations on VHF signals during pre-midnight hours, associated with range type of spread-F on ionograms and a relatively weak and slow fading scintillations during post-midnight hours associated with frequency type of spread-F is observed during the relatively high sunspot years from 2001 to 2004, whereas during the low sunspot years from 2005 to 2010 the scintillation activity as well as spread-F activity are found to be minimum. During both the high and low sunspot years, it is observed that the maximum scintillation activity occurs during equinoctial months followed by winter with the minimum occurrence during summer months. The annual mean percentage occurrence of scintillations is found to be clearly associated with the variations in the annual mean sunspot number. The nocturnal variations in the occurrence of scintillations show the onset of scintillation activity starts from 19:00 h LT with maximum of occurrence around 21:00 h LT. A clear semiannual variation in the occurrence of scintillations is observed during pre-midnight hours with two peaks in equinoctial months of March/April and September/October. The number of scintillation patches observed is found to be more during pre-midnight hours compared to those during post-midnight hours. The most probable scintillation patch duration lies around 30 min. Further, it is also found that the number of scintillation patches with durations of 60 min and more decreases with the increase in the patch duration. It is also observed in general that the scintillation activity is inhibited during geomagnetic disturbed days.     

Intercosmos-19 observations of an additional topside ionization layer: the F3 layer

Spatial properties of an additional ionization layer in the topside ionosphere were investigated using Intercosmos-19 satellite ionograms. The data under analysis were choosen for equinoctial conditions of the high solar activity period (1979–1981). The F3 layer was detected in a narrow longitude sectors (about 60°) between the equatorial anomaly crests. Its intensity has a maximum just above the equator and decreases poleward within ±10° dip. A nighttime F3 layer was observed as well as the daytime events.     

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.     

同步卫星讯号显示的电离层闪烁特性

本文利用1983年5—8月,1984年5—12月在武昌(114.4°E,30.6°N)对日本ETS-Ⅱ卫星(130.0°E)发出的136.1124MHz讯号的观测资料进行了统计分析。结果表明:(1)武昌电离层闪烁不但有日变化,而且有季变化。每年5—7月为闪烁最大活动期,在这些月份的夜间常出现法拉弟旋转角类波扰动伴随有强闪烁现象。武昌电离层闪烁是属于中纬闪烁型;(2)闪烁指数与法拉弟旋转角起伏密切相关,它们出现率之间的相关系数为0.8以上;夜间闪烁与扩展F层,白天闪烁与突发E层出现率之间的相关系数分别为0.6和0.55。      

Seismoionospheric fountain-effect as analogue of active space experiment

Space investigation includes modelling and initiation of natural processes in ionospheric-magnetospheric plasma by active experimental methods with artificial influence as a specific part. For the first time we present here the possibility of triggering the unique equatorial ionosphere phenomenon known as fountain-effect by lithospheric processes before an earthquake. Anomalous influence of earthquake preparation processes with epicenter close to magnetic equator on latitudinal distribution of f_oF - F-region penetration frequencies and spread-F is discussed and considered by us as an analogue of an active experiment. Using Alouette-1 data it is shown that approximately one day before the earthquake f_oF2 dependence on magnetic latitude looks like a double-crest curve with minimum close to epicenter. Simul taneously oblique and spread traces appeared at topside ionograms. The effect is evident during evening hours. Such kind of f_oF2-dependence is unusual and is not observed during the preceding days, nor some days later. So, a possible reason for the observations might be connected with generation of an anomalous electric field near the epicenter during earthquake preparation stage, which initiates a phenomenon similar to the natural fountain-effect.     

海口站与Huancayo站扩展F的差异

本文统计分析了我国海口站太阳黑子高年和低年的扩展F的资料,得出了该站扩展F的出现率随地方时、季节和太阳活动周期的某些变化规律并与Huancayo站扩展F的资料进行了比较,指出了两站扩展F的差异。      

Nighttime ionosphere–thermosphere coupling observed during an intense geomagnetic storm

The electrodynamics of the ionosphere in the tropical region presents various scientific aspects, which remain subject of intensive investigations and debates by the scientific community. During the year 2002, in a joint project between the Universidade do Vale do Paraíba (UNIVAP) and Universidade Luterana do Brasil (ULBRA), a chain of three Canadian Advanced Digital Ionosondes (CADIs) was established nearly along the geomagnetic meridian direction, for tropical ionospheric studies, such as, changes and response due to geomagnetic disturbances and thermosphere–ionosphere coupling and the generation and dynamics of ionospheric irregularities, in the Brazilian sector. The locations of the three ionosondes stations are São José dos Campos (23.2°S, 45.9°W, dip latitude 17.6°S – under the southern crest of equatorial ionospheric anomaly), Palmas (10.2°S, 48.2°W, dip latitude 5.5°S – near the magnetic equator) and Manaus (2.9°S, 60.0°W, dip latitude 6.4°N – between the geographic and geomagnetic dip equators). It should be pointed out that Palmas and Manaus are located on the opposite sides of the magnetic equator but both are south of the geographic equator. The three CADIs work in time-synchronized mode and obtain ionograms every 5 min. This configuration of the ionospheric sounding stations allowed us to study the F-region dynamics during geomagnetically disturbed period in the meridional direction. Just after the installation and testing of the three CADIs, on September 05, 2002 a coronal mass ejection (CME) left the Sun and about 2 days after the CME left the Sun, it reached the Earth’s magnetosphere and complex and multi step events took place during the period September 07–09. In the study we note that the equatorial stations located north (Manaus, dip latitude 6.4°N) and south (Palmas, dip latitude 5.5°S) of the dip equator presented significant F-layer height asymmetries during the storm main phase. In addition, the low-latitude station SJC (dip latitude 17.6°S) presented decrease in the F-layer densities (negative phase), whereas Palmas presented increase in the F-layer densities (positive phase) during the main phase. This was followed by positive phase at both the stations. During the first night of the recovery phase a strong formation and evolution of large-scale ionospheric irregularities (equatorial spread-F (ESF)) was observed, but on the second night of the recovery phase, there was strong and almost simultaneous sporadic E (Es) formation at all three stations. During the presence of Es, spread-F formation is not observed, indicating the suppression of spread-F, possibly by sporadic E.     

Predicting the probability of occurrence of spread-F over Brazil using neural networks

The probability of occurrence of spread-F can be modeled and predicted using neural networks (NNs). This paper presents a feasibility study into the development of a NN based model for the prediction of the probability of occurrence of spread-F over selected equatorial stations within the Brazilian sector. The input space included the day number (seasonal variation), hour (diurnal variation), sunspot number (measure of the solar activity), magnetic index (measure of the magnetic activity) and magnetic position. Twelve years of spread-F data from Brazil (covering the period 1978–1989) measured at the equatorial site Fortaleza (3.9°S, 38.45°W) and low latitude site Cachoeira Paulista (22.6°S, 45.0°W) are used in the development of an input space and NN architecture for the model. Spread-F data that is believed to be related to plasma bubble developments (range spread-F) was used in the development of the model. The model results show the probability of spread-F occurrence as a function of local time, season and latitude. Results from the Brazilian Sector NN (BSNN) based model are presented in this paper, as well as a comparative analysis with a Brazilian model developed for the same purpose.     

IRI-2007对扩展F的预测与观测比较研究

通过对赤道地区两台站的扩展F 实际观测值与IRI-2007 模型预测值的比较, 研究了该模型预测结果的地域局限性和预测的准确度. 该模型建立在美洲扇区赤道附近台站资料基础上, 本文选取非洲扇区大体同纬度两个台站资料进行对比研究. 统计结果表明, 无论在发生概率日均值的年变化趋势上, 还是在扩展F 发生率年均值随地方时变化以及随太阳活动的统计中, 都存在着明显的差异; 而对非洲两站实测值的分析和IRI-2007 模型对南美巴西地区的预测又分别与当前已有的研究结果相符. 分析结果表明, (1) 建立在巴西台站数据基础上的IRI-2007 模型的统计建模极具参考价值; (2) 赤道和低纬度电离层扩展F 现象有强烈的经度效应, 在该模型建模方法基础上, 结合不同扇区实际资料的统计, 能改进该模型全球应用的普适性; (3) IRI-2007模型能正确反映巴西地区扩展F 的发生概率, 用该模型和其他扇区不同条件下的实测值对比, 提供了一个细致研究扩展F经度效应的基础, 有利于确认导致扩展F发生的基本条件和多种影响因素的具体作用.      

Investigation of Satellite Trace (ST) and Multi-reflected Echo (MRE) ionogram signatures and its possible correlation to nighttime spread F development from Cyprus over the solar mini-max (2009–2016)

Multi-reflected echoes (MREs) and satellite traces (STs) are referred in literature as ionogram signatures of Travelling Ionospheric Disturbances (TIDs) which is a phenomenon that apparently drives spread F development mainly at nighttime mid-latitude ionosphere. A long-term statistical study has been undertaken to investigate the morphological aspect of these signatures over the lower midlatitude European station of Nicosia, Cyprus (35.19°N, 33.38°E geographic; magnetic dip. 29.38°N) by inspecting all ionograms recorded by the DPS-4D digisonde in the interval 2009–2016. The results underline the systematic manifestation of these TID signatures over Cyprus with a possible (although not quite clear) solar activity dependence and a distinctive seasonal and diurnal occurrence rate with a seasonal peak of STs during summer and of MREs during January to April. Based on the experimental results of the present study, the seasonal occurrence rate of MREs and STs is found to increase by 75% and 56% during high solar activity periods. Satellite traces are well known ionogram signatures of TIDs and mostly correlated to the nighttime spread F formation. The occurrence of mid-latitude spread Fs over European longitude sector normally increases during summer. The occurrences of TIDs are also prominent at this interval of the year over nighttime mid-latitude ionosphere. The presence of MREs as an ionogram signature of TIDs over mid-latitude ionosphere is unique in nature.     

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.     

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.     

Role of neutral winds in generating irregularities in equatorial F-region

Electron density and neutral wind velocity measurements were carried out by rocketborne probes from rocket ranges in India. The experiments were carried out at the time of the onset of spread-F at sunset hours. The results show that a neutral wind velocity in north-south direction greater than 100 m/sec is required to trigger spread-F. It is suggested that spread-F is generated by the interaction of neutral gas with ionospheric plasma.