A new inversion algorithm for HF sky-wave backscatter ionograms

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 F₂ layer. Consequently, while the inversion results can provide accurate profiles of the F region below the F₂ 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 F₂ 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 F₂ 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.     

Application of neural networks to South African GPS TEC modelling

The propagation of radio signals in the Earth’s atmosphere is dominantly affected by the ionosphere due to its dispersive nature. Global Positioning System (GPS) data provides relevant information that leads to the derivation of total electron content (TEC) which can be considered as the ionosphere’s measure of ionisation. This paper presents part of a feasibility study for the development of a Neural Network (NN) based model for the prediction of South African GPS derived TEC. The South African GPS receiver network is operated and maintained by the Chief Directorate Surveys and Mapping (CDSM) in Cape Town, South Africa. Vertical total electron content (VTEC) was calculated for four GPS receiver stations using the Adjusted Spherical Harmonic (ASHA) model. Factors that influence TEC were then identified and used to derive input parameters for the NN. The well established factors used are seasonal variation, diurnal variation, solar activity and magnetic activity. Comparison of diurnal predicted TEC values from both the NN model and the International Reference Ionosphere (IRI-2001) with GPS TEC revealed that the IRI provides more accurate predictions than the NN model during the spring equinoxes. However, on average the NN model predicts GPS TEC more accurately than the IRI model over the GPS locations considered within South Africa.     

Electron density profile calculation technique for Autoscala ionogram analysis

An electron density profile model with free parameters is introduced. Initially the parameters are calculated on the basis of the ionospheric characteristics automatically obtained from the ionograms by Autoscala and considering the helio-geophysical conditions. The technique used to adjust the free parameters to the particular ionograms recorded is presented.     

Spatial distribution of TEC across India in 2005: Seasonal asymmetries and IRI prediction

Total electron content measured simultaneously at 10 locations over India during the low solar activity year 2005 is used to examine the temporal and spatial asymmetries and also to assess the predictability of the International Reference Ionosphere in respect of the observed asymmetrical distribution. The stations are distributed in latitude along 77°E and in longitude along 23°N forming a meridional and a zonal chain respectively. A longitudinal gradient positive towards east was observed in the daytime hours of equinox and summer. Equinoctial asymmetry was prevalent across India during this year. Within the crest and equator, winter anomaly has been observed. It is found that IRI 2012 (with Ne Quick option, URSI coefficients) is unable to fully capture the temporal variation and spatial gradients of the ionization density in the Indian sector during 2005. The amount of offset between the model and measurement varies with local time and location.     

The solar-terrestrial event of 23 February 1956

The solar flare of 23 February 1956 and the resulting geophysical disturbance ranks as one of the most remarkable solar-terrestrial events of the twentieth century. It sparked many papers and has seldom been equalled. Fifty years after the International Geophysical Year, it seems timely to review the observations of the event from today’s perspective, and to draw on the recollections of scientists who were active at the time.     

Influence of solar-geomagnetic disturbances on SABER measurements of 4.3 μm emission and the retrieval of kinetic temperature and carbon dioxide

Thermospheric infrared radiance at 4.3 μm is susceptible to the influence of solar-geomagnetic disturbances. Ionization processes followed by ion-neutral chemical reactions lead to vibrationally excited NO⁺ (i.e., NO⁺(v)) and subsequent 4.3 μm emission in the ionospheric E-region. Large enhancements of nighttime 4.3 μm emission were observed by the TIMED/SABER instrument during the April 2002 and October–November 2003 solar storms. Global measurements of infrared 4.3 μm emission provide an excellent proxy to observe the nighttime E-region response to auroral dosing and to conduct a detailed study of E-region ion-neutral chemistry and energy transfer mechanisms. Furthermore, we find that photoionization processes followed by ion-neutral reactions during quiescent, daytime conditions increase the NO⁺ concentration enough to introduce biases in the TIMED/SABER operational processing of kinetic temperature and CO₂ data, with the largest effect at summer solstice. In this paper, we discuss solar storm enhancements of 4.3 μm emission observed from SABER and assess the impact of NO⁺(v) 4.3 μm emission on quiescent, daytime retrievals of Tk/CO₂ from the SABER instrument.     

Observations of HF propagation on a path aligned along the mid-latitude trough

Observations of the direction of arrival and time of flight of HF signals propagating on a 1400 km path oriented along the mid-latitude trough are presented. At night, the signal commonly arrives from directions offset from the great circle bearing by up to 80° and these events have been categorised into five main types. Statistics indicating how often these categories of propagation were observed in the period August 2006 to September 2007 are presented. The physical mechanisms which result in the off great circle propagation are also discussed.     

Global empirical models of the density peak height and of the equivalent scale height for quiet conditions

Monthly average electron density profiles have been calculated from hourly electron density N(h) recorded in 26 digisonde stations distributed worldwide encompassing the time interval 1998–2006. The ionospheric electron density peak height of the F2 region, hmF2, and the effective scale height at the hmF2, Hm, deduced from average profiles have been analyzed to obtain the quiet-time behavior and have been analytically modeled by the spherical harmonic analysis (SH) technique using the modip latitude as the coordinate of the reference system. The coefficients of the SH models of hmF2 and Hm are bounded to the solar activity, and the temporal and seasonal variations are considered by Fourier expansion of the coefficients. The SH models provide a tool to predict hmF2 and Hm located anywhere in the range of latitudes between of 70°N and 70°S and at any time. The SH analytical model for hmF2 improves the fit to the observations by 10% in average compared to the IRI prediction, and it might improve the IRI prediction of hmF2 by more than 30% at high and low latitudes. The analytical model for Hm predicts the quiet behavior of the effective scale height with accuracy better than 15% in average which enables to obtain a good estimation of vertical profiles. These results could be useful to estimate information for the topside profile formulation.     

Interaction of oblique wave beam with ionospheric layer F2

Interaction of a powerful obliquely incident wave beam of decameter radio waves with the ionospheric F2 layer is analyzed. Much like the linear case, propagation through the natural anti-waveguide layer F2 splits the initial beam. Some part of its energy leaks through the ionospheric layer, the other part goes back along a downward trajectory. However, nonlinearity leads to further stratification of the ionospheric layer. A new feature, in comparison with the linear case, is appearing a narrow waveguide beneath the F2 layer maximum which traps a small part of the beam energy. We study the relationship between these parts of the wave field in a simplified model of parabolic F2 layer, with nonlinearity caused by thermal plasma expulsion from the high field intensity region. Analytical results are supplemented with numerical estimates of the effects.     

A simulation analysis of systematic errors in ion arrival angle and drift velocity from the satellite borne ion drift meter

As one payload of a Chinese seismic satellite program, an ion drift meter (IDM) will measure drift velocity of thermal ions at an altitude of 500 km. Previous works have shown that such instruments use biased grids to create nonuniform potential in the grid planes, which brings systematic errors to the inferred parameters. A commercial finite element analysis software is used to simulate this instrument in the exact size. The error sources from thermal velocity, nonuniform transparency of real grids and potential depression in the grid planes are explained. The simulation results show that the arrival angle and drift velocity will be underestimated in all the conditions and the maximal error will be about −0.87° and −121 m/s, respectively. Furthermore, the relative error of the inferred arrival angle and the drift velocity will be inversely correlated with the arrival angle because of the lensing effect of the potential depression. This simulation provides a quantificational method of evaluating and correcting the data during in situ operation.