Modelling of ionospheric drifts in view of IRI

D/E region drift data, as worldwide obtained from 1957 to 1970 by methods D1 and D3, are statistically analyzed and a computer program describing the average variations in geomagnetic latitude and time is developed.     

Validating ionospheric models with measured electron density profiles

On behalf of an URSI Working Group 3 initiated study (VIM), three ionospheric models, IRI, PL/PRISM and FLIP, are compared with electron density profiles derived from ionograms Millstone Hill. Four months of data in 1989/90 were analyzed. For most of the time, N(h) profiles were available every 15 minutes providing a good statistical database for the evaluation of the ionospheric models in terms of diurnal and seasonal variations.     

Modelling studies of ionospheric convection in northern and southern polar regions

The realistic model of Quegan et al. has been used to investigate the convection paths of ionospheric plasma at 300 km altitude, for different polar cap radii and in both hemispheres. Taking the Northern magnetic dip pole to be at a co-latitude of 11° and the Southern magnetic dip pole at a co-latitude of 23°, these paths are presented in a Sun-Earth frame, with the position of the Earth's axis fixed as it is on 21 March, as polar plots centred on the magnetic pole. There are marked hemispheric differences between 13 and 23 L.T., particularly near the stagnation region at 18 to 21 L.T., but only minor differences between 00 and 12 L.T., when the radius of the polar cap exceeds 12°. For a smaller polar cap, the differences between the hemispheres are small at all local times. The time taken to perform a complete circuit is most dependent on the polar cap radius, and most variable - between 15 and 36 h - for convection paths starting near 60° latitude. The time that plasma convecting from noon to near midnight across the Northern polar cap spends within the 10° co-latitude circle increases from 6 h, for a polar cap radius of 10°, to 11.5 h at 17°. These results are compared and contrasted with other model calculation results and with some ground-based and satellite observations of plasma densities at high latitudes.     

Comparison of different N(h) profiles with mid-latitude ionospheric observations

After inversion, N(h) profiles obtained from ionograms that had been recorded during high solar activity at two mid-latitutde stations have been compared with those derived from IRI90 and DGR ionospheric models. A small data set has been selected such that both geomagnetically quiet and disturbed conditions are represented.     

Modelling ionospheric vertical drifts over Africa low latitudes using Empirical Orthogonal functions and comparison with climatological model

For the first time, empirical model of daytime vertical $\mathbf{E}\times \mathbf{B}$ drift based on Empirical Orthogonal functions (EOF) decomposition technique is presented. Day-to-day variability of $\mathbf{E}\times \mathbf{B}$ drift inferred from horizontal (H) geomagnetic field data around dip latitude for the period of 2008–2013 is used to both develop and validate the model. Results show that the EOF technique is promising with modelled values and data giving correlation coefficient values of at least 0.90 for geomagnetic conditions of both $K_p?3$ and $K_p>3$ within 2008–2013. Independent model validation shows that in situ $\mathbf{E}\times \mathbf{B}$ values from ion velocity meter (IVM) instrument on-board C/NOFS satellite are closer to model $\mathbf{E}\times \mathbf{B}$ estimates than the climatological Scherliess-Fejer (SF) model incorporated within the International Reference Ionosphere (IRI).     

Accuracy of IRI profiles of ionospheric density and temperatures derived from comparisons to Kharkov incoherent scatter radar measurements

The incoherent scatter radar (ISR) facility in Kharkov, Ukraine (49.6°N, 36.3°E) measures vertical profiles of electron density, electron and ion temperature, and ion composition of the ionospheric plasma up to 1100 km altitude. Acquired measurements constitute an accurate ionospheric reference dataset for validation of the variety of models and alternative measurement techniques. We describe preliminary results of comparing the Kharkov ISR profiles to the international reference ionosphere (IRI), an empirical model recognized for its reliable representation of the monthly-median climatology of the density and temperature profiles during quiet-time conditions, with certain extensions to the storm times. We limited our comparison to only quiet geomagnetic conditions during the autumnal equinoxes of 2007 and 2008. Overall, we observe good qualitative agreement between model and data both in time and with altitude. Magnitude-wise, the measured and modeled electron density and plasma temperatures profiles appear different. We discovered that representation accuracy improves significantly when IRI is driven by observed-averaged values of the solar activity index rather than their predictions. This result motivated us to study IRI performance throughout protracted solar minimum of the 24th cycle. The paper summarizes our observations and recommendations for optimal use of the IRI.     

Active stimulation of the ionospheric plasma

Laboratory experiments in which high power, pulsed electromagnetic waves interact with an inhomogeneous plasma indicate that the generated nonlinear plasma phenomena depend on peak incident power and not on pulse length. The electromagnetic waves can penetrate beyond the cutoff and produce large, enhanced electrostatic fields at the critical layer within 100 electron plasma periods. The enhanced electric field pressure can be comparable to the thermal pressure and can accelerate ions and electrons to velocities much greater than their thermal speed. Large density cavities (with δn/n ? 10%) can be created in a time shorter than the usual ion response time because of the accelerated ion dynamics. These laboratory results have been extended to create a new and generalized concept to actively stimulate space plasmas with high power pulses of short duration. A field experiment will be used for the stimulation of auroral ionospheric plasma. The ground-based system is modular, each module consisting of a 2 MW pulsed HF transmitter designed at UCLA and a crossed-dipole antenna element. Incoherent scatter radar and optical diagnostic methods are discussed.     

Modeling of ionospheric scintillation at low-latitude

The presence and movement of plasma density fluctuations in the F-region of the ionosphere are studied by monitoring phase and amplitude of radio waves propagating through the region. In this paper, we have used weak scattering theory and assumed the plasma density fluctuations to behave like phase changing diffraction screen. Appropriate relations for scintillation index S₄, and phase variance δ? are derived and computed for different parameters of the plasma density irregularities of the ionosphere. SROSS-C2 satellite in situ measurements of plasma density fluctuations, which provide direct information about the structure and morphology of irregularities that are responsible for scintillation of radio waves, were used first time to develop a scintillation model for low latitude. It is observed that the scintillation index S₄ and phase variance δ? depends on the strength of the plasma turbulence. Finally, the results obtained from modeling are compared and discussed with the available recent results.     

Magnetospheric and ionospheric manifestations of interplanetary foreshocks

The geometry of a typical interplanetary shock front in the vicinity of the Earth’s orbit predicts that the leading edge of the foreshock region comes into contact with the magnetosphere a few hours ahead of geomagnetic sudden impulses (SI). There is reason to believe that the interaction of the magnetosphere with the foreshock leads to magnetic and ionospheric disturbances, which can be detected by ground-based instruments. We searched for specific precursors of SIs in data from the Scandinavian riometer network and in the short period geomagnetic pulsation data from mid-latitude magnetometers. We found that SIs were preceded by the following three features: (1) an increase in riometric absorption, (2) excitation of Pcl magnetic pulsations and (3) a spectral broadening of the Pc3 magnetic pulsations. Our observations may be useful for the study of acceleration processes in the solar wind. These observations are also of potential forecasting interest.     

Lower-mesospheric inversion layers over brazilian equatorial region using TIMED/SABER temperature profiles

Lower-mesospheric inversion layers (MILs) were studied using the temperature profiles observed by TIMED/SABER over Cariri (7.5°S, 36.5°W), Brazil, in 2005. A total 175 MILs were identified with the maximum occurrence in April and October and the minimum in January and July. The lower MIL is located in a height region from 70 to 90 km, with the peak at around 83 ± 4 km with the temperature of 205 ± 5 K, and the thickness of 4–10 km. The results show large amplitudes of MILs during equinoxes and minimum in solstices, with a clear semiannual variation. A general feature of lower MIL in monthly mean profile was observed twice a year, one from February to May, and the other from August to October with a downward shift of the top level. These results suggest that formation and long persistence of MIL is an important factor to investigate propagation of atmospheric gravity waves in the mesosphere-lower thermosphere (MLT) region.     

Evaluation of COMPASS ionospheric model in GNSS positioning

As important products of GNSS navigation message, ionospheric delay model parameters are broadcasted for single-frequency users to improve their positioning accuracy. GPS provides daily Klobuchar ionospheric model parameters based on geomagnetic reference frame, while the regional satellite navigation system of China’s COMPASS broadcasts an eight-parameter ionospheric model, COMPASS Ionospheric Model(CIM), which was generated by processing data from continuous monitoring stations, with updating the parameters every 2 h. To evaluate its performance, CIM predictions are compared to ionospheric delay measurements, along with GPS positioning accuracy comparisons. Real observed data analysis indicates that CIM provides higher correction precision in middle-latitude regions, but relatively lower correction precision for low-latitude regions where the ionosphere has much higher variability. CIM errors for some users show a common bias for in-coming COMPASS signals from different satellites, and hence ionospheric model errors are somehow translated into the receivers’ clock error estimation. In addition, the CIM from the China regional monitoring network are further evaluated for global ionospheric corrections. Results show that in the Northern Hemisphere areas including Asia, Europe and North America, the three-dimensional positioning accuracy using the CIM for ionospheric delay corrections is improved by 7.8%–35.3% when compared to GPS single-frequency positioning ionospheric delay corrections using the Klobuchar model. However, the positioning accuracy in the Southern Hemisphere is degraded due apparently to the lack of monitoring stations there.     

The calculation of ionospheric absorption with modern computers

New outcomes are proposed for ionospheric absorption starting from the Appleton–Hartree formula, in its complete form. The range of applicability is discussed for the approximate formulae, which are usually employed in the calculation of non-deviative absorption coefficient. These results were achieved by performing a more refined approximation that is valid under quasi-longitudinal (QL) propagation conditions. The more refined QL approximation and the usually employed non-deviative absorption are compared with that derived from a complete formulation. Their expressions, nothing complicated, can usefully be implemented in a software program running on modern computers. Moreover, the importance of considering Booker’s rule is highlighted. A radio link of ground range D = 1000 km was also simulated using ray tracing for a sample daytime ionosphere. Finally, some estimations of the integrated absorption for the radio link considered are provided for different frequencies.     

Estimated errors in a global gravity wave climatology from GPS radio occultation temperature profiles

In a previous paper by Schmidt et al. (2008), from CHAllenging Minisatellite Payload (CHAMP) Global Positioning System (GPS) radio occultation data, a comparison was made between a Gaussian filter applied to the “complete” temperature profile and to its “separate” tropospheric and stratospheric height intervals, for gravity wave analyses. It was found that the separate filtering method considerably reduces a wave activity artificial enhancement near the tropopause, presumably due to the isolation process of the wave component. We now propose a simple approach to estimate the uncertainty in the calculation of the mean specific wave potential energy content, due exclusively to the filtering process of vertical temperature profiles, independently of the experimental origin of the data. The approach is developed through a statistical simulation, built up from the superposition of synthetic wave perturbations. These are adjusted by a recent gravity wave (GW) climatology and temperature profiles from reanalyses. A systematic overestimation of the mean specific wave potential energy content is detected and its variability with latitude, altitude, season and averaging height interval is highlighted.     

Computerized ionospheric tomography with the IRI model

Computerized ionospheric tomography (CIT) is a method to estimate ionospheric electron density distribution by using the global positioning system (GPS) signals recorded by the GPS receivers. Ionospheric electron density is a function of latitude, longitude, height and time. A general approach in CIT is to represent the ionosphere as a linear combination of basis functions. In this study, the model of the ionosphere is obtained from the IRI in latitude and height only. The goal is to determine the best representing basis function from the set of Squeezed Legendre polynomials, truncated Legendre polynomials, Haar Wavelets and singular value decomposition (SVD). The reconstruction algorithms used in this study can be listed as total least squares (TLS), regularized least squares, algebraic reconstruction technique (ART) and a hybrid algorithm where the reconstruction from the TLS algorithm is used as the initial estimate for the ART. The error performance of the reconstruction algorithms are compared with respect to the electron density generated by the IRI-2001 model. In the investigated scenario, the measurements are obtained from the IRI-2001 as the line integral of the electron density profiles, imitating the total electron content estimated from GPS measurements. It has been observed that the minimum error between the reconstructed and model ionospheres depends on both the reconstruction algorithm and the basis functions where the best results have been obtained for the basis functions from the model itself through SVD.     

Comparing ionospheric models with mid-latitude ionosonde observations

The purpose of this research work is to validate the ionospheric models (IRI and CHIU) to assess its suitability and usefulness as an operational tool. The ionospheric model is a computer model designed to predict the state of the global ionosphere for 24 h. The scope was limited to conduct comparisons between the predicted F₂ layer critical frequencies (f₀F₂) against observed ionosonde data. The ionospheric prediction model (IPM) was designed to predict by using monthly median sunspot number, while the observation data are taken from two digital ionospheric sounding stations (Okinawa, 26.28N, 127.8E and Wakkanai, 45.38N, 141.66E) which lies within the mid-latitude region of the globe. Analysis of the f₀F₂ data from stations for year (2001) with high solar activity and year (2004) with low solar activity, four months (March, June, September and December) chosen based primarily on data availability. From results it seen that the ratio between monthly median predicted and observed f₀F₂ values for each model used in this research work and for the chosen months was nonlinear with local time, so the empirical formula for applying correction factors were determined, these formula can be used to correct the error occurred in predicted f₀F₂ value.     

Total electron content measurements in ionospheric physics

With the advent of modern global networks of dual-frequency Global Positioning System (GPS), total electron content (TEC) measurements along slant paths connecting GPS receivers and satellites at 22,000 km have become the largest data set available to ionospheric scientists. The TEC can be calculated from the time and phase delay in the GPS signal using the GPS Toolkit, but an unknown bias will remain. In addition, UHF/VHF radio beacons on board low-Earth-orbiting satellites can also be used to measure the electron content. However, the TEC measurements are obtained by integrating TEC differences between slant paths, but also contain biases. It is often necessary to use data assimilative algorithms like the Ionospheric Data Assimilation Three-Dimensional (IDA3D), and to treat both GPS- and LEO-beacon TEC measurements as relative data in order to conduct ionospheric studies.     

A different view of the ionospheric winter anomaly

The Winter Anomaly (WA) is an ionospheric phenomenon, particularly related to the F2-layer, that is classically defined as the situation in which the ionization level during winter is higher than during summer for a certain location. This situation is anomalous because it contradicts what would be naturally expected with summer’s ionization levels higher than winter ones. This phenomenon has been a matter of study since the early decades of the XX century.This contribution tackles the study of the WA based on long time series (up to two solar cycles) of Total Electron Content (TEC) measurements over a globally distributed network of stations. The work done relies on a threefold strategy including: a classical approach based on the comparison of the winter and summer maximum TEC values, that confirmed the results previously documented; an intermediate approach with maximum TEC values modelled as linear functions of the solar radiation level, that allowed to identify stations where the WA is likely to be observed (and at what solar radiation level) and where is not; a final original approach where the maximum TEC values are modelled using Chapman’s function for the seasonal variations, linear regressions for the solar activity dependence and the introduction of site and month dependent equivalent ionization coefficients for the remaining effects. The main conclusions are that the coefficients for winter months seem to depend mainly on geomagnetic latitude, to increase towards high latitudes and to be asymmetric between hemisphere. In accordance to these findings, the occurrences of WA effects would be regulated by proper combinations of these three effects.     

Impact of heliogeophysical disturbances on ionospheric HF channels

The article presents the results of the observation of a strong magnetic storm and two X-ray flares during the summer solstice in 2015, and their impact on the HF signals characteristics in ionospheric oblique sounding. It was found that the negative phase of the magnetic storm led to a strong degradation of the ionospheric channel, ultimately causing a long blackout on paths adjacent to subauroral latitudes. On mid-latitude paths, the decrease in 1FMOF reached ～50% relative to the average values for the quiet ionosphere. It is shown that the propagation conditions via the sporadic Es layer during the magnetic storm on a subauroral path are substantially better than those for F-mode propagation via the upper ionosphere. The delay of the sharp decrease in 1FMOF during the main phase of the magnetic storm allowed us to determine the propagation velocity of the negative phase disturbances (～100 m/s) from subauroral to mid-latitude ionosphere along two paths: Lovozero – Yoshkar-Ola and Cyprus – Nizhny Novgorod. It is shown that both the LOF and the signal/noise ratio averaged over the frequency band corresponding to the propagation mode via the sporadic Es layer correlate well with the auroral AE index. Using an over-the-horizon chirp radar with a bistatic configuration on the Cyprus – Rostov-on-Don path, we located small-scale scattering irregularities responsible for abnormal signals in the region of the equatorial boundary of the auroral oval.     

Main ionospheric trough studied from intercosmos-Bulgaria-1300

The main ionospheric trough is a phenomenon in the mid and high latitude ionosphere, characterized by an abrupt decrease of the electron and ion density and increase of the electron temperature. We here examine the behaviour of the trough for different geomagnetic conditions based on data from the Intercosmos-Bulgaria-1300 satellite. The dependence of the polar trough wall on the boundary of the soft electron precipitation from the plasma layer is shown, and a possible universal time dependence is examined.