Search

Large-scale ionospheric irregularities, which occurred on September 8, 2017, were detected by the ionosonde and the Global Navigation Satellite System (GNSS) receiver network in this letter. The Wuhan ionospheric sounding system (WISS), deployed at Puer (PUR, 22.7°N, 101.05°E), Leshan (LSH, 29.6°N, 103.7°E), and Zhangye (ZHY, 39.4°N, 100.13°E), detected these large-scale ionospheric irregularities. The 2-D maps of the rate of the GNSS total electron content (TEC) index (ROTI) further confirmed that the large-scale ionospheric irregularities can cover a broad area (from the equator to middle latitude, ~80°E to ~120°E, ~20°N to ~45°N). The ionospheric variation recorded by the WISS shows that the large-scale ionospheric irregularities were driven by the super eastward electric field, which was formed by the eastward prompt penetration electric field during a storm superimposed on the normal prereversal enhancement during the postsunset hours. Moreover, observations by the WISS show the higher the latitudes at which the ionospheric irregularities/spread F occurred, the shorter the duration time.

A novel land-based high-frequency (HF) geolocation system is proposed. The geolocating signals are transmitted from the stations with known geographic coordinates. Because of the ionospheric refraction, the signals propagate as skywave in the over-the-horizon (OTH) manner. Therefore, a wide spatial coverage can be achieved. The receiving station uses a small 2-D antenna array to estimate the angle of arrival (AOA) of each transmitter’s signal. An intersecting geolocation algorithm based on the geometric relations between the transmitters, and the receiver is described, which only relies on the coordinates of the transmission stations and the azimuth angles of the signals. A novel specific multichannel system structure is also described. It has a simple structure and low hardware cost, as well as good flexibility. The preliminary experimental results prove the feasibility of this geolocation method and its engineering application significance.

The generation and evolution of artificial plasma clouds is a complicated process that is strongly dependent on the background environment and release conditions. In this paper, based on a three-dimensional two-species fluid model, the evolution characteristics of artificial plasma clouds under various release conditions were analyzed numerically. In particular, the effect of ionospheric density gradient and ambient horizontal wind field was taken into account in our simulation. The results show that an asymmetric plasma cloud structure occurs in the vertical direction when a nonuniform ionosphere is assumed. The density, volume, and expansion velocity of the artificial plasma cloud vary with the release altitude, mass, and initial ionization rate. The initial release velocity can change the cloud's movement and overall distribution. With an initial velocity perpendicular to the magnetic field, an O+ density cavity and two bumps exist. When there is an initial velocity parallel to the magnetic field, the generated plasma cloud is bulb-shaped, and only one O+ density cavity and one density bump are created. Compared to the cesium case, barium clouds expand more rapidly. Moreover, Cs+ clouds have a higher density than Ba+ clouds, and the snowplow effect of Cs+ is also stronger.

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.

As a companion paper to Zhou RX et al. (2020) , this study describes application of the automatic detection and analysis module to identify all the tweek atmospherics detectible in the WHU ELF/VLF receiver data collected at Suizhou station during the period of 3 February through 29 February 2016. Detailed analysis of the identified low-latitude tweek events reveals that the occurrence rate varies considerably — from 800 to 6000 tweeks per day, and exhibits a strong diurnal and local time dependence, the peak occurring before local midnight. The diurnal variation of identified tweeks was similar to that of the lightning data obtained by the World-Wide Lightning Location Network (WWLLN).. Estimates of the propagation distance and ionospheric reflection height of tweek atmospherics suggest that the majority (~92%) of the low latitude tweeks originate from the lightning activity within a radius of 4000 km and that they are very likely to reflect from the lower ionospheric D-region at the height range of 75–85 km. At these lower ionospheric reflection altitudes, ~74% of the corresponding electron densities from the tweek spectral measurements are within 24.5–27.5 cm-3. The daily variation of estimated D-region electron densities in the considered period (February 2016) also exhibits a small overall increasing trend from early to later in the month.

As a dispersive wave mode produced by lightning strokes, tweek atmospherics provide important hints of lower ionospheric (i.e., D-region) electron density. Based on data accumulation from the WHU ELF/VLF receiver system, we develop an automatic detection module in terms of the maximum-entropy-spectral-estimation (MESE) method to identify unambiguous instances of low latitude tweeks. We justify the feasibility of our procedure through a detailed analysis of the data observed at the Suizhou Station (31.57°N, 113.32°E) on 17 February 2016. A total of 3961 tweeks were registered by visual inspection; the automatic detection method captured 4342 tweeks, of which 3361 were correct ones, producing a correctness percentage of 77.4% (= 3361/4342) and a false alarm rate of 22.6% (= 981/4342). A Short-Time Fourier Transformation (STFT) was also applied to trace the power spectral profiles of identified tweeks and to evaluate the tweek propagation distance. It is found that the fitting accuracy of the frequency–time curve and the relative difference of propagation distance between the two methods through the slope and through the intercept can be used to further improve the accuracy of automatic tweek identification. We suggest that our automatic tweek detection and analysis method therefore supplies a valuable means to investigate features of low latitude tweek atmospherics and associated ionospheric parameters comprehensively.

Statistical analysis of spread F (SF), recorded at Puer (PUR, 22.7° N, 101.05° E, Dip Latitude 12.9° N) during 2015–2016, was carried out to reveal its characteristics at the northern equatorial ionization anomaly. In our study, SF was categorized into four types, frequency spread F (FSF), range spread F (RSF), mix spread F (MSF) and strong range spread F (SSF). The statistical results presented that FSF and MSF were dominant over Puer. Most types of SF appeared mostly in summer months, except the maximum occurrence of SSF in equinox months. Moreover, observations of SF events also showed that the solar activity and magnetic activity dependence of SF varied with seasons. Compared with observations at other region, the present results suggest that medium-scale traveling ionospheric disturbances (MSTIDs) may play a key role in generation of SF in low-latitude region.

Benefiting from the special design of BeiDou Navigation Satellite System (BDS) constellation, its Geostationary Earth Orbit (GEO) satellites observations can provide favorable conditions to investigate the ionospheric plasma bubble at low latitudes. This study proposes using BDS GEO Rate of Total Electron Content Index observations to analyze the equatorial plasma bubble characteristics. We use GEO satellites instead of ground-based GNSS stations as a reference and compare the onset time of plasma bubbles observed by different GEO satellites, i.e., the reference GEO, the west and east GEO satellites, during geomagnetically quiet days, and determine plasma bubbles that are generated locally or drifted from elsewhere. According to this strategy, it is found that there is a significant difference in the occurrence rates of plasma bubbles generated locally over two closely located stations, i.e., LALX (18.19 °N, 104.98 °E; geomagnetic latitude: 11.31 °N) and YONG (16.83 °N, 112.34 °E; geomagnetic latitude: 9.96 °N) during the whole year of 2014. Statistical results indicate that during March and September equinoxes, the occurrence rates of plasma bubbles generated locally at the eastern station YONG are 52% and 34%, while at the western station LALX, they are only 34% and 25%. Further analysis reveals a close relationship between the higher bubble generation rates and the active atmospheric inter-tropical convergence zone.

We present the observational evidence of $E$ - and $F$ -region field-aligned irregularities (FAIs), nighttime medium-scale traveling ionospheric disturbances (MSTIDs), and a sporadic $E$ ( $E_{\text {S}}$ )-layer using the Wuhan very-high-frequency (VHF) coherent scatter radar, Wuhan Global Navigation Satellite System (GNSS) network, and Wuhan ionosonde. We observed simultaneously E and F FAIs by VHF radar and the $E_{\text {S}}$ -layer and spread- $F$ by Wuhan ionosonde. We also observed MSTIDs using the Wuhan GNSS network in a southwestward direction of propagation and horizontal propagation velocity of less than 180 m/s, for a period of ~33 min. Simultaneous observations of an $E_{\text {S}}$ -layer and FAIs in the $E$ -region and FAIs in the $F$ -region suggested the existence of an electrodynamic coupling between the $E$ and $F$ regions in the mid- and low-latitude nighttime ionosphere. A polarized electric field associated with nighttime MSTIDs can generate the uplift movements of the $F$ -region’s electron density. This uplift effect consequently can excite the gradient drift instability (GDI) with an accompanying enhanced vertical gradient of the $F$ -layer’s electron density. Our results indicated that the $F$ -region MSTIDs excited by Perkins instability might be further modulated by the $E \times B$ effect and subsequently can evolve to spread-F like FAIs. Our observational investigation provided the first evidence of the full dynamics and links among different ionospheric disturbances in a mid- and low-latitude nighttime region of China.

We report a new observation of daytime ionospheric irregularities on December 27, 2009, at low latitudes in China using multistatic HF radar. The transmitter, with a horizontal polarization log-periodic antenna pointing westward, is located at Wanning (18.97° N, 110.50° E, dip angle 19°). Three separate receivers are in Sanya (18.20° N, 109.49° E), Qiongzhong (19.04° N, 109.77° E), and Tunchang (19.42° N, 110.13° E). Doppler spectra, frequency power spectra, and wavenumber power spectra from the three receivers were calculated. The main features of these observations are: 1) the irregularities observed at 12:10 (LT) had typical lengths on the scale of tens of meters and the altitude was in the range of 105–175 km, which is located in the E region and bottom of the F region; 2) the power distributions from the three receivers are consistent with the irregularities being field aligned; and 3) features consistent with a cascade of irregular spatial scales were observed. We suggest that the observed irregularity could be due to the intermediate layer modulated by electrodynamic instability.

By using the 2015–2016 data set of newly built Wuhan very high-frequency (VHF) radar, we investigate the occurrence of E-region field-aligned irregularities (FAIs) in mid-latitude China region for the first time. Wuhan E-region FAIs present the strongest occurrence of quasi-periodic (QP) echoes at postsunset period in the height range of 100–120 km and moderate occurrence of continuous echoes at postsunrise period in the height range of 90–100 km. Characteristics including temporal and altitudinal distribution, Doppler spectra and period of QP echoes are presented from a statistical point of view. We find that E-region FAIs observed by Wuhan VHF radar show strong linkage with the concurrent sporadic E layer observed by ionosonde. In general, Doppler spectra present type 2 feature with velocity in the range of ±50 m/s and width less than 50 m/s. Our results also indicate that mid-latitude E-region FAIs are closely associated with Es layer. Several mechanisms were proposed to interpret the mid-latitude E-region FAIs; however, further investigation is required to understand the generation with different ionospheric background and irregularity types.

This study presents observations of nighttime spread F/ionospheric irregularities and spread Es at low and middle latitudes in the South East Asia longitude of China sectors during the recovery phase of the 7–9 September 2017 geomagnetic storm. In this study, multiple observations, including a chain of three ionosondes located about the longitude of 100°E, Swarm satellites, and Global Navigation Satellite System (GNSS) ROTI maps, were used to study the development process and evolution characteristics of the nighttime spread F/ionospheric irregularities at low and middle latitudes. Interestingly, spread F and intense spread Es were simultaneously observed by three ionosondes during the recovery phase. Moreover, associated ionospheric irregularities could be observed by Swarm satellites and ground-based GNSS ionospheric TEC. Nighttime spread F and spread Es at low and middle latitudes might be due to multiple off-vertical reflection echoes from the large-scale tilts in the bottom ionosphere. In addition, we found that the periods of the disturbance ionosphere are ~1 h at ZHY station, ~1.5 h at LSH station and ~1 h at PUR station, respectively. It suggested that the large-scale tilts in the bottom ionosphere might be produced by LSTIDs (Large scale Traveling Ionospheric Disturbances), which might be induced by the high-latitude energy inputs during the recovery phase of this storm. Furthermore, the associated ionospheric irregularities observed by satellites and ground-based GNSS receivers might be caused by the local electric field induced by LSTIDs.

This paper describes a novel sounding system for which the functions of the medium frequency (MF) radar and the ionosonde are integrated on the same hardware platform and antenna structure, namely the middle atmosphere-ionosphere (MAI) system. Unlike the common MF radar, MAI system adopts the pseudo-random (PRN) phase-coded modulation technology, which breaks the limitation of the traditional monopulse mode. Through the pulse compression, only a small peak power is needed to achieve the signal-to-noise ratio (SNR) requirement. The excellent anti-jamming performance is also very suitable for the ionospheric sounding. One transmitting and six receiving modes are adopted for the MF sounding. While neglecting the structure of the T/R switches, the coupling interference between the transmitter and the receiver may also be avoided. Moreover, by employing a miniaturized antenna array composed of progressive-wave antennas for the MF receiving and ionospheric sounding, the MAI system takes account of the requirements of the inversion algorithms of MF radar and the large bandwidth need for the ionospheric sounding concurrently. Such an antenna structure can also greatly simplify the system structure and minimize the difficulty of deployment. The experiments verified the availability of the system scheme and its engineering application significance. Through further analysis of the sounding data, the wind field of the mesosphere, the electron density of D layer and electron density profile from layers E to F were obtained at the identical location. The capability of MAI system can play an important role in studying the interaction and coupling mechanism between the mesosphere and ionosphere.

In this paper, a new dedicated multi-channel ionospheric oblique backscatter sounder is described. A miniaturized L-shaped antenna array was employed for the receiving of the oblique backscatter echoes in the present system. Firstly, two typical adaptive beamforming algorithms were introduced to improve the anti-jamming ability. Then, simulations were carried out to verify the beamforming performance in azimuth and elevation simultaneously. Furthermore, the experimental results by the present sounding system were used to test the performance of the adaptive beamformers. Results show that the radio frequency interference and the interference of the vertical echoes can be effectively suppressed by the adaptive beamformers. In particular, the use value of the beamforming in the receiving of the ionospheric oblique backscatter sounding is described in detail through the analysis of the signal sequences in several typical frequencies. And after the constant false-alarm rate processing, the oblique backscatter ionograms processed by the adaptive beamforers have clearer and continuous leading edge compared with the original ionograms. As a result, the adaptive beamformers of great significance to improve the detection ability of the ionospheric oblique backscatter sounding.

Spread F is a commonly observed phenomenon on ionograms caused by plasma irregularities or wave-like structures in the ionosphere. In general, Spread F could induce fluctuations of amplitude and phase of radio waves which travel through the ionosphere. Therefore, investigation of Spread F could be used to not only reveal ionospheric electrodynamics process, but also have a significant engineering application. Due to a large amount of ionograms recorded by ionosondes, it is a challenge work to manually identify ionograms with Spread F to study characteristics of Spread F. Thus, much work has been devoted to automatic identification of Spread F. In the present study, a machine learning method related to decision tree was adopted to automatically identify Spread F from ionograms. First, ionograms were processed by image method and projection techniques to provide input parameters for decision tree. The output of the proposed decision tree is whether Spread F is present or not on ionograms. Then, a set of ionograms was used to construct a decision tree. At last, a set of ionograms was adopted to validate the performance of this decision tree. In this study, ionograms recorded at Puer station (Geographic latitude and longitude: 22.7°N, 101.5°E; Geomagnetic latitude: 12.8°N) in the Yunnan province were used. Results indicate that the decision tree performed well in automatic identification of Spread F on ionograms. The accuracy of automatic identification in a set of ionograms with Spread F was reached up to 89%. It inspires us to continually improve the performance of automatic identification of Spread F in the future work.

Recently, Jiang et al., (2016) reported daytime spread F at Puer station (22.7oN, 101.05oE, Dip Lat 12.9oN). As a continual work, daytime spread F at middle latitudes was investigated in this study through ionograms recorded by an ionosonde installed at Zhangye (ZHY, 39.4oN, 100.0oE, Dip Lat 29.6oN) in the Northwest of China. The ionosonde at ZHY station recorded daytime spread F three times on 6 January, 2017. The first case occurred at approximately 06:45 LT and disappeared at about 08:35 LT. The second and third cases lasted from 09:50 LT to 10:30 LT and from 11:10 LT to 11:30 LT, respectively. In addition, Swarm satellite observations were used to reveal the possible mechanism in this study. Results show that daytime spread F observed in this study might be attributed to Traveling Ionospheric Disturbances (TIDs)/atmospheric gravity waves induced by a geomagnetic storm. However, the physical processes might be different for these cases. Ionospheric instabilities produced by local F region electric field might be contributed to the first case. The second and third cases might be attributed to off-vertical reflections from wave-like ionospheric structure generated by TIDs/atmospheric gravity waves.

In this paper, we investigate the ionospheric heating by oblique incidence of powerful high-frequency (HF) radio waves using three-dimensional numerical simulations. The ionospheric electron density and temperature perturbations are examined by incorporating the ionospheric electron transport equations and ray-tracing algorithm. The energy distribution of oblique incidence heating waves in the ionosphere is calculated by the three-dimensional ray-tracing algorithm. The calculation takes into consideration the electric field of heating waves in the caustic region by the plane wave spectral integral method. The simulation results show that the ionospheric electron density and temperature can be disturbed by oblique incidence of powerful radio waves, especially in the caustic region of heating waves. The oblique ionospheric heating with wave incidence parallel and perpendicular to the geomagnetic field in the mid-latitude ionosphere is explored by simulations, results of which indicate that the ionospheric modulation is more effective when the heating wave propagates along the magnetic field line. Ionospheric density and temperature striations in the caustic region due to thermal self-focusing instability are demonstrated, as well as the time evolution of the corresponding fluctuation spectra. Keywords. Ionosphere (active experiments)

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

We investigate the electric field penetration of the lithosphere–atmosphere–ionosphere coupling (LAIC) problem to study abnormal seismo-ionospheric disturbance. By directly solving the LAIC electric field penetration model at the high-latitude region, we find that the additional current induced at the ground surface flows into the ionosphere completely and further generates an abnormal ionospheric electric field. Therefore, we reasonably suggest that the electric field penetration of LAIC at middle- and low-latitude regions can be solved from the perspective of the ionospheric electric field model. The current from the downward atmosphere is treated as the source term. The simulation results demonstrate the following principal findings: (a) for the high-latitude region, the horizontal electric field in the ionosphere does not change with height and the vertical electric field can be neglected; (b) for the middle- and low-latitude regions, the intensity of the total horizontal electric field increases with the latitude and the vertical electric field is more obvious at low latitudes; and (c) the penetration height of the LAIC electric field in the ionosphere is lower at low latitudes than at high latitudes. We also find that according to the diurnal change of the ionospheric conductivity, the most efficient time for electric field penetration is between 00:00 and 04:00 local time.

This study mainly investigates equatorial and low-latitude ionospheric response to a great geomagnetic storm occurred on the 17 March 2015. We found that there were some interesting ionospheric phenomena, e.g. short-term ionospheric positive effect, daytime spread F, and morning Equatorial Ionization Anomaly (EIA) in the topside ionosphere, emerged at equatorial and low-latitude region along the longitude of about 100oE. Ground-based ionosondes and in situ satellite (Swarm) were utilized to study the possible mechanisms for these ionospheric phenomena. We found that vertical downward transport of plasma or neutral induced by traveling ionospheric disturbances (TIDs) or traveling atmospheric disturbances (TADs) might make a contribution to the short-term ionospheric positive effect at the main stage of this great storm. Additionally, results suggested that the occurrence of daytime spread F at low latitudes might be due to the diffusion of equatorial ionospheric irregularities in the topside ionosphere along geomagnetic field lines. Moreover, observational evidence shows that TIDs also might be the main driver for morning EIA-like feature recorded by Swarm B satellite in the topside ionosphere. These ionospheric phenomena mentioned above could make us to gain a better understanding of ionospheric storm effects at equatorial and low-latitude region.

Electromagnetic extremely low frequency (ELF) waves play an important role in modulating the Earth’s radiation belt electron dynamics. High-frequency (HF) modulated heating of the ionosphere acts as a viable means to generate artificial ELF waves. The artificial ELF waves can reside in two different plasma regions in geo-space by propagating in the ionosphere and penetrating into the magnetosphere. As a consequence, the entire trajectory of ELF wave propagation should be considered to carefully analyze the wave radiation properties resulting from modulated ionospheric heating. We adopt a model of full wave solution to evaluate the Poynting vector of the ELF radiation field in the ionosphere, which can reflect the propagation characteristics of the radiated ELF waves along the background magnetic field and provide the initial condition of waves for ray tracing in the magnetosphere. The results indicate that the induced ELF wave energy forms a collimated beam and the center of the ELF radiation shifts obviously with respect to the ambient magnetic field with the radiation power inversely proportional to the wave frequency. The intensity of ELF wave radiation also shows a weak correlation with the size of the radiation source or its geographical location. Furthermore, the combination of ELF propagation in the ionosphere and magnetosphere is proposed on basis of the characteristics of the ELF radiation field from the upper ionospheric boundary and ray tracing simulations are implemented to reasonably calculate magnetospheric ray paths of ELF waves induced by modulated ionospheric heating.

A software tool named ionoScaler (ionogram Scaler) was developed to carry out manual and automatic scaling of vertical incidence ionograms. Additionally, ionosonde Total Electron Content (ITEC) can be estimated by this tool through reconstructing the bottomside and topside electron density profiles of ionograms. Ionograms recorded at Puer (22.7°N, 101.05°E) and Wuhan (30.5°N, 114.3°E) were adopted to test features of the ionoScaler. Results show that the software tool performed well for scaling of ionograms and estimation of ITEC. However, bad quality ionograms would reduce the accuracy of automatic scaled (autoscaled) values. The ionoScaler presented in this paper will greatly expedite the interpretation of ionograms, and it is also an alternative method for estimation of TEC values. Moreover, this tool will be useful for research in the ionosphere and space weather.

We report observations of ionospheric disturbances in response to a North Korean underground nuclear explosion (UNE) on 3 September 2017. By using data from IGS (International GNSS Service) stations and Swarm satellites, geomagnetic conjugate ionospheric disturbances were observed. The observational evidence showed that UNE-generated ionospheric disturbances propagated radially from the UNE epicenter with a velocity of ∼280 m s−1. We propose that the ionospheric disturbances are results of electrodynamic process caused by LAIC (lithosphere–atmosphere–ionosphere coupling) electric field penetration. The LAIC electric field can also be mapped to the conjugate hemispheres along the magnetic field line and consequently cause ionospheric disturbances in conjugate regions. The UNE-generated LAIC electric field penetration plays an important role in the ionospheric disturbances in the region of the nuclear test site nearby and the corresponding geomagnetic conjugate points.

We report a statistical study of medium-scale traveling ionospheric disturbances (MSTIDs) during 2014–2017 in the low-latitude China region using the Hong Kong Continuously Operating Reference Stations network with a baseline length of 10–15 km. Polynomial fitting and the multichannel maximum entropy method are utilized to derive MSTID parameters from global positioning system total electron content data. We find that MSTIDs can be sorted into three types in this region. One type is daytime MSTIDs, which are mainly distributed during 1000–1700 local time (LT) in spring, autumn and winter. The second type is nighttime MSTIDs, which mainly occur from 2200 to 0300 LT in spring and summer. The third type is morning MSTIDs, mainly occurring from 0500 to 0700 LT in spring and autumn. No clear difference in the MSTID parameters is distinguishable between the cyclone period and non-cyclone period, which may be due to the distance restriction for cyclone identification and different propagation distances of the primary gravity waves (GWs) and secondary GWs.

In this paper, complete complementary code (CCC) sequences are applied to a High Frequency (HF) ionospheric sounding network. Ionosondes distributed at multiple locations use the mutually orthogonal CCC sequences to conduct vertical soundings synchronously. At the same time, thanks to the omnidirectional antennas, every station can receive the oblique echoes transmitted from the others. Due to the orthogonality between the code sequences, both vertical and oblique ionograms can be simultaneously obtained and completely separated. Through this method, the sounding efficiency can be enhanced, and the inversion difficulty can be reduced. Further, by using the data assimilation method, vertical and oblique sounding results can be combined to obtain a wide range of regional ionospheric characteristics. To verify the performance of this kind of sounding network, validation experiments are implemented to demonstrate that vertical and oblique ionograms can be obtained independently at the same time and inverted separately and that the maps of foF2 parameters obtained by using the data assimilation method provide more details than single vertical or oblique sounding.

Motivated by the observations and linear theory analysis of ionospheric irregularities at Mars, we performed numerical simulations of the nonlinear evolution of the electromagnetic gradient drift instability in the Martian ionospheric dynamo region. The seeding source of ionospheric irregularities is perturbation zonal neutral wind. We found that the perturbation electric fields induced by the gradient drift instability can convect lower density plasma into higher density plasma at higher altitudes. Then, the associated perturbation magnetic field and electric field can cause the velocity shear of the plasma, which induced the Kelvin–Helmholtz instability at higher altitude. The Kelvin–Helmholtz instabilities furthermore lead to smaller-scale irregularities in plasma density, magnetic field, and electric field in the Martial ionosphere. Key points: (1) Nonlinear simulation of small-scale ionospheric irregularities at Mars was present. (2) Ionospheric irregularities at Mars can be seeded by the perturbation neutral winds. (3) Model results are comparable to linear theory analysis and satellite observations.

A local maximum value of ionospheric total electron content (TEC) in midlatitude regions, named midlatitude ionospheric peak in this study, has been frequently observed during both the daytime and nighttime. We have carried out a modeling study of midlatitude ionospheric peaks in TEC using the Thermosphere Ionosphere Electrodynamic General Circulation Model (TIEGCM) to investigate possible mechanisms driving these peaks. Comparisons between model results and Global Position System (GPS) TEC observations in the solar minimum year of 2009 show that TIEGCM can reproduce the nighttime and daytime midlatitude ionospheric peaks. Model results indicate that midlatitude ionospheric peaks in the daytime and nighttime are likely caused by different mechanisms. Nighttime midlatitude ionospheric peaks are attributed to the downwards-moving plasma ambipolar diffusive flux from the plasmasphere. However, the daytime upwards-moving plasma ambipolar diffusive flux and poleward meridional winds, which can cause ionospheric depletion, might play dominated roles in producing midlatitude ionospheric peaks in the early morning in the winter season.

This paper presents a new method to reconstruct the vertical electron density profile based on vertical Total Electron Content (TEC) using the simulated annealing algorithm. The present technique used the Quasi-parabolic segments (QPS) to model the bottomside ionosphere. The initial parameters of the ionosphere model were determined from both International Reference Ionosphere (IRI) (Bilitza et al., 2014) and vertical TEC (vTEC). Then, the simulated annealing algorithm was used to search the best-fit parameters of the ionosphere model by comparing with the GPS-TEC. The performance and robust of this technique were verified by ionosonde data. The critical frequency (foF2) and peak height (hmF2) of the F2 layer obtained from ionograms recorded at different locations and on different days were compared with those calculated by the proposed method. The analysis of results shows that the present method is inspiring for obtaining foF2 from vTEC. However, the accuracy of hmF2 needs to be improved in the future work.

The ionograms, which were recorded by the ionosonde located at Pu’er station (PUR, 22.7°N, 101.05°E, Dip Latitude 12.9°N) in the Southwest of China in the year of 2016, were used to study the ionospheric behavior of the ordinary critical frequency of the F2 layer (foF2) in the region of the northern equatorial ionization anomaly. To verify the performance of the International Reference Ionosphere (IRI) over the Southwest of China, a comparative study of the observed foF2 and the latest version of the International Reference Ionosphere (IRI-2016) was carried out. We found that the foF2 in equinox months is greater than summer and winter. Moreover, a higher frequency of the observed bite-out of foF2 in January and April than other months and the IRI-2016 cannot represent the bite-out of foF2 in diurnal variations. Compared to the observations at Pu’er Station, the IRI-2016 underestimated foF2 for most time of the year. The IRI with the International Radio Consultative Committee (CCIR) option overestimated foF2 is higher than that with the International Union of Radio Science (URSI) option. Furthermore, the normalized root mean square error of foF2 from the IRI-2016 with the CCIR option is less than that with the URSI.

Powerful high-frequency (HF) radio waves can be used to efficiently modify the upper-ionospheric plasmas of the F region. The pressure gradient induced by modulated electron heating at ultralow-frequency (ULF) drives a local oscillating diamagnetic ring current source perpendicular to the ambient magnetic field, which can act as an antenna radiating ULF waves. In this paper, utilizing the HF heating model and the model of ULF wave generation and propagation, we investigate the effects of both the background ionospheric profiles at different latitudes in the daytime and nighttime ionosphere and the modulation frequency on the process of the HF modulated heating and the subsequent generation and propagation of artificial ULF waves. Firstly, based on a relation among the radiation efficiency of the ring current source, the size of the spatial distribution of the modulated electron temperature and the wavelength of ULF waves, we discuss the possibility of the effects of the background ionospheric parameters and the modulation frequency. Then the numerical simulations with both models are performed to demonstrate the prediction. Six different background parameters are used in the simulation, and they are from the International Reference Ionosphere (IRI-2012) model and the neutral atmosphere model (NRLMSISE-00), including the High Frequency Active Auroral Research Program (HAARP; 62.39° N, 145.15° W), Wuhan (30.52° N, 114.32° E) and Jicamarca (11.95° S, 76.87° W) at 02:00 and 14:00 LT. A modulation frequency sweep is also used in the simulation. Finally, by analyzing the numerical results, we come to the following conclusions: in the nighttime ionosphere, the size of the spatial distribution of the modulated electron temperature and the ground magnitude of the magnetic field of ULF wave are larger, while the propagation loss due to Joule heating is smaller compared to the daytime ionosphere; the amplitude of the electron temperature oscillation decreases with latitude in the daytime ionosphere, while it increases with latitude in the nighttime ionosphere; both the electron temperature oscillation amplitude and the ground ULF wave magnitude decreases as the modulation frequency increases; when the electron temperature oscillation is fixed as input, the radiation efficiency of the ring current source is higher in the nighttime ionosphere than in the daytime ionosphere.

'Many observations in the ionospheric heating experiment, by a powerful high frequency electromagnetic wave with ordinary polarization launched from a ground-based facility, is attributed to parametric instability (PI). In this paper, the general dispersion relation and the threshold of the PI excitation in the heating experiment are derived by considering the inhomogeneous spatial distribution of pump wave field. It is shown that the threshold of PI is influenced by the effective electron and ion collision frequencies and the pump wave frequency. Both collision and Landau damping should be considered in the PI calculation. The derived threshold expression has been used to calculate the required threshold for excitation of PI for several ionospheric conditions during heating experiments conducted employing EISCAT high frequency transmitter in Tromso, Norway, on 2nd October 1998, 8th November 2001, 19th October 2012 and 7th July 2014. The results indicate that the calculated threshold is in good agreement with the experimental observations. © 2018 Hefei Institutes of Physical Science, Chinese Academy of Sciences and IOP Publishing.

In this letter, a novel ionospheric oblique-incidence sounding network consisting of the ionospheric oblique backscatter sounder and the parasitic oblique-incidence sounder is presented. With the new configuration, the oblique ionogram and the backscatter ionogram can be provided at the same time. Due to the usage of a novel waveform combining the pseudorandom phase coding and the chirp modulation, about 130-dB signal processing gain can be usually obtained to improve the signal-to-noise ratio for the low-power transmission. Furthermore, through choosing different types of pseudorandom sequences for each transmitter, many oblique ionograms that depict the information of different propagation paths can be simultaneously achieved at a single receiver site. In addition, the digital beamforming technique is implemented on the receiving antenna array to simultaneously receive the echoes produced by each transmitter, to achieve a good space isolation to enhance the echo discrimination, and to suppress the strong interferers appearing in the high-frequency band. Details of the system configuration, the ambiguity function of the sounding waveform, and the signal processing algorithm are described. The initial experimental results show that multiple oblique ionograms can be simultaneously obtained at a single receiver site with this sounding network.

Automatic scaling ionogram can get the parameters of ionogram which are vital to ionosphere detecting. In this paper, a new method is proposed to scale F2 layer trace automatically from oblique ionogram based on morphological operator and inversion technique. This method is verified through the comparison of actual detecting data with statistical analysis. The results show that the proposed automatic scaling method has high acceptable rate and is suitable for scaling oblique ionogram with different high angle wave states. It is fast and precise to fit O-mode echoes in F2 layer without the influence from F1 layer. This method could be applied in real-time ionospheric oblique sounding research with high reliability and versatility.

The F2 layer critical frequency of the ionosphere (foF2) is one of the most significant parameters for studying the ionosphere. To investigate the large-scale characteristics of the ionosphere over particular regions, modeling foF2 is an effective method. In this paper, we use both the Kriging (KG) and neural network (NN) methods to reconstruct foF2 maps over North China. The neural network is trained by the genetic algorithm (GA) to avoid the ‘local minimum’ phenomenon in most NN applications. We then carry out a comparison between foF2 provided by both the KG and NN methods with vertical model operation of ionosonde data including Beijing, Qingdao, Suzhou, and Changchun. All of the foF2 data used in the comparison are obtained from the oblique and vertical mode operation of ionosonde from the China Ground-based Seismo-ionospheric Monitoring Network. To allow for a possible seasonal and diurnal variation, data obtained from summer, winter, and equinox months are applied in the present comparison. In addition, we make a comparison during a magnetic storm period. The results of our comparisons demonstrate that both the KG and NN methods are appropriate tools for modeling foF2 maps. However, when the data set is spare, the performance of the NN method is better than the KG method. On the other hand, the KG method is more robust than the NN method during a magnetic storm.

In this paper, a novel oblique-incidence ionosonde (Wuhan Ionospheric Oblique-Incidence Sounding System) and its new application in the localization of the ionospheric irregularities are presented. Due to the usage of the binary-phase-coded waveform, a large signal processing gain, a high Doppler and range resolution, and a large unambiguous detection range can be achieved in this ionosonde. This ionosonde also adopts the peripheral component interconnect extensions for instruments (PXI) bus technology and is designed as a small-sized PXI-based system. Furthermore, a high-performance oven-controlled crystal oscillator that is disciplined by the Global Positioning System is used to achieve a good time and frequency synchronization. With multichannel digital receiver and multiple receiving sites, this ionosonde can be applied in the localization of the ionospheric irregularities. The details of the system configuration, the ambiguity function of the sounding waveforms, the signal processing algorithm, and the time and frequency synchronization method are described. The experimental results show that the virtual height along with the ground position of the ionospheric field-aligned irregularities can be preliminarily localized with this ionosonde.

The solar eclipse effects on the ionosphere are very complex. Except for the ionization decay due to the decrease of the photochemical process, the couplings of matter and energy between the ionosphere and the regions above and below will introduce much more disturbances. Five ionosondes in the Northeast Asia were used to record the midlatitude ionospheric responses to the solar eclipse of 20 May 2012. The latitude dependence of the eclipse lag was studied first. The foF2 response to the eclipse became slower with increased latitude. The response of the ionosphere at the different latitudes with the same eclipse obscuration differed from each other greatly. The plasma flux from the protonsphere was possibly produced by the rapid temperature drop in the lunar shadow to make up the ionization loss. The greater downward plasma flux was generated at higher latitude with larger dip angle and delayed the ionospheric response later. The waves in the foEs and the plasma frequency at the fixed height in the F layer are studied by the time period analytic method. The gravity waves of 43–51 min center period during and after the solar eclipse were found over Jeju and I-Cheon. The northward group velocity component of the gravity waves was estimated as ~108.7 m/s. The vertical group velocities between 100 and 150 km height over the two stations were calculated as ~5 and ~4.3 m/s upward respectively, indicating that the eclipse-induced gravity waves propagated from below the ionosphere.

Many observational results have shown electromagnetic abnormality in the ionosphere before large earthquakes. The theoretical simulation can help us to understand the internal mechanism of these anomalous electromagnetic signals resulted from seismic regions. In this paper, the horizontal and vertical components of electric and magnetic field at the topside ionosphere are simulated by using the full wave method that is based on an improved transfer matrix method in the lossy anisotropic horizontally stratified ionosphere. Taken account into two earthquakes with electric field perturbations recorded by the DEMETER satellite, the numerical results reveal that the propagation and penetration of ULF (ultra-low-frequency) electromagnetic waves into the ionosphere is related to the spatial distribution of electron and ion densities at different time and locations, in which the ion density has less effect than electron density on the field intensity. Compared with different frequency signals, the minimum values of electric and magnetic field excited by earthquakes can be detected by satellite in current detection capability have also been calculated, and the lower frequency wave can be detected easier.

For decades, high-frequency (HF) radar has played an important role in sensing the Earth’s environment. Advances in radar technology are providing opportunities to significantly improve the performance of HF radar, and to introduce more applications. This paper presents a low-power, small-size, and multifunctional HF radar developed by the Ionospheric Laboratory of Wuhan University, referred to as the Wuhan Ionospheric Oblique Backscattering Sounding System (WIOBSS). Progress in the development of this radar is described in detail, including the basic principles of operation, the system configuration, the sounding waveforms, and the signal and data processing methods. Furthermore, its various remote sensing applications are briefly reviewed to show the good performance of this radar. Finally, some suggested solutions are given for further improvement of its performance.

In this letter, a new high-frequency (HF) skywave multistatic radar that is used for ionospheric sensing and aircraft detection is presented. With multiple receiving sites, this radar can detect the aircraft in much larger surveillance with high probability and can obtain the characteristics of many ionospheric paths simultaneously. In addition, this radar adopts a binary-phase-coded continuous waveform to achieve a high Doppler resolution, a large signal processing gain, and a low probability of interception. Details of the system configuration, the ambiguity function of the sounding waveform, and the signal processing algorithm are described. The experimental results show that backscatter ionograms and oblique ionograms can be simultaneously obtained and that the aircraft with different flying conditions can be clearly detected with an HF skywave multistatic experimental radar.

Oblique backscatter sounding is a powerful tool for detecting and monitoring the ionosphere continuously at a remote distance. High-frequency (HF) backscatter ionograms provide the amplitudes of backscatter signals with respect to group path or time delay against operating frequency. Application of inversion algorithm to a backscatter ionogram can extract useful information regarding the ionospheric electron density along the propagation paths. The present study proposes a new inversion algorithm on basis of simulated annealing method to acquire the leading edge of sweep-frequency ionogram, which is subsequently validated by ionospheric vertical sounding data. Quantitative comparisons between the vertical sounding measurements and the inversion results obtained from oblique backscatter sounding indicate that the new algorithm enables us to overcome the instability issue that traditional inversion algorithm faces and output reliable information of ionospheric inversion with satisfactory efficiency, thus providing a robust alternative for ionospheric detection based on oblique backscatter ionograms especially when the ionosphere is calm with slow changes.

Vertical incidence ionograms indicate ionospheric characteristics over the ionosonde station, from which electron density profiles can be derived. This paper describes a method for the automatic calculation of electron density profiles from vertical incidence ionograms. First, the method calculates the initial parameters of the quasi-parabolic segments (QPS) model by using the International Reference Ionosphere (IRI) model, the Nequick2 model, image processing techniques, and the Empirical Orthogonal Function (EOF). Once the initial parameters have been calculated the method then adjusts those to obtain the electron density profiles and synthesized traces that match the recorded ionograms. The algorithm then selects the best-fit synthesized trace and corresponding parameters as output from the candidate ones. Furthermore, the corresponding electron density profile of the recorded ionogram is calculated by using the best-fit parameters of the QPS model. To further test the feasibility of the proposed method, we apply it to some ionograms that were recorded at Wuhan during different seasons. As a result, our results demonstrate that the proposed method is feasible for the automatic calculation of electron density profiles.

[1] Recent advances in computing technologies have renewed interest in the intelligent systems for automatic interpretation of ionograms, images obtained by remote sensing of the ionospheric plasma. The ionogram “autoscaling” techniques based on the template matching method, previously rendered unrealistic for their computing complexity, have now become feasible. This work presents an automatic scaling technique for extracting the main features of the F1 and F2 layers of the ionosphere, such as critical frequency and virtual height, from vertical incidence ionograms that do not distinguish O/X polarization of the echoes. The proposed technique uses the quasi-parabolic segments (QPS) to model the electron density profile shapes that are used to synthesize a pool of candidate traces. Moreover, the empirical orthogonal functions and image technique are applied to reduce the size of the candidate traces so that the auto-scaling algorithm can run in realistic time. With the template matching algorithm, the technique will provide the initial parameters of the QPS model for the F1 and F2 layers, which are then fine-tuned to obtain the better fitting parameters. In order to evaluate the performance of this technique, a large data set of ionograms recorded in Wuhan at daytime and nighttime in winter, summer, and equinoctial months, are analyzed and investigated. The automatic scaling results are compared with manually scaling results. Our results indicate that the proposed technique described in this paper is reliable and efficient and will facilitate the statistical study of temporal and spatial ionospheric characteristics over Wuhan.

The Wuhan Ionospheric Oblique Backscattering Sounding System (WIOBSS) is a monostatic high-frequency sky-wave radar used for ionospheric remote sensing. The sweep frequency backscatter ionogram (SFBI) recorded by the WIOBSS contains both backscatter echo scattered by distant terrestrial surface and vertical incidence (VI) echo reflected by local ionosphere over the sounding station. The approach for SFBI inversion introduced in this letter requires input of leading edge and peak height derived from local VI echo. The final output of this SFBI inversion approach is the 2-D electron density distribution in a vertical plane aligned in the direction of sounding. In addition, the time and geographic variation of foF2 can be obtained from a series of SFBI inversion results. Two experiments have been utilized to validate the SFBI inversion approach, and the SFBI inversion results are found to be very close to the ionosonde data. The fast-converged feature of the approach makes it possible for applications in real time. The success of the SFBI inversion approach enables us to use the WIOBSS as a vehicle-mounted system to obtain ionosphere electron density profile over a large geographic area.