Comparison of model and experimental ionospheric parameters at high latitudes

Modern use and study of the auroral region needs to attract a wider class of models for describing conditions of radio wave propagation in the ionosphere. In this paper the possibilities of the International Reference Ionosphere model, well-proven and widespread in the mid-latitudes, are investigated in the high latitude zone. Model and measured values of the critical frequency foF2 for two mid-latitude stations (Juliusruh and Goosebay) and four high-latitude ones (Loparskaya, Sodankyla, Sondrestrom, Thule) are compared. Deviations of medians, variations from day to day and solar activity trends seemed to be similar for both areas. This similarity is irrespective of the RZ12 index. Special attention is paid to the TEC parameter and its determination using 6 versions of models, a new version of the model IRI2010 (IRI-Plas) among them. It is shown that the IRI-Plas model significantly improves the definition of TEC in contrast to the versions of IRI2007 and the new model NeQuick. The use of the median of the experimental equivalent slab thickness, together with the current values of the TEC, increases by a factor of two the agreement between calculated and measured foF2 values as compared with the variations from day to day. This allows foF2 to be defined in near-real time.     

Longitudinal behaviors of the IRI-B parameters of the equatorial electron density profiles retrieved from FORMOSAT-3/COSMIC radio occultation measurements

The electron density profiles in the bottomside F₂-layer ionosphere are described by the thickness parameter B0 and the shape parameter B1 in the International Reference Ionosphere (IRI) model. We collected the ionospheric electron density (N_e) profiles from the FORMOSAT-3/COSMIC (F3/C) radio occultation measurements from DoY (day number of year) 194, 2006 to DoY 293, 2008 to investigate the daytime behaviors of IRI-B parameters (B0 and B1) in the equatorial regions. Our fittings confirm that the IRI bottomside profile function can well describe the averaged profiles in the bottomside ionosphere. Analysis of the equatorial electron density profile datasets provides unprecedented detail of the behaviors of B0 and B1 parameters in equatorial regions at low solar activity. The longitudinal averaged B1 has values comparable with IRI-2007 while it shows little seasonal variation. In contrast, the observed B0 presents semiannual variation with maxima in solstice months and minima in equinox months, which is not reproduced by IRI-2007. Moreover, there are complicated longitudinal variations of B0 with patterns varying with seasons. Peaks are distinct in the wave-like longitudinal structure of B0 in equinox months. An outstanding feature is that a stable peak appears around 100°E in four seasons. The significant longitudinal variation of B0 provides challenges for further improving the presentations of the bottomside ionosphere in IRI.     

Feasibility of developing an ionospheric E-region electron density storm model using TIMED/SABER measurements

We present a new technique for improving ionospheric models of nighttime E-region electron densities under geomagnetic storm conditions using TIMED/SABER measurements of broadband 4.3 μm limb radiance. The response of E-region electron densities to geomagnetic activity is characterized by SABER-derived NO⁺(v) 4.3 μm Volume Emission Rates (VER). A storm-time E-region electron density correction factor is defined as the ratio of storm-enhanced NO⁺(v) VER to a quiet-time climatological average NO⁺(v) VER, which will be fit to a geomagnetic activity index in a future work. The purpose of this paper is to demonstrate the feasibility of our technique in two ways. One, we compare storm-to-quiet ratios of SABER-derived NO⁺(v) VER with storm-to-quiet ratios of electron densities measured by Incoherent Scatter Radar. Two, we demonstrate that NO⁺(v) VER can be parameterized by widely available geomagnetic activity indices. The storm-time correction derived from NO⁺(v) VER is applicable at high-latitudes.     

Near real-time assimilation in IRI of auroral peak E-region density and equatorward boundary

The paper describes the method and initial results of assimilating the auroral peak E-region density (NmE) and the auroral equatorward boundary (EB) into the International Reference Ionosphere (IRI). The NmE and EB are obtained using a FUV based auroral model or FUV measurements in near real-time. Initial results show that the auroral NmE is often significantly larger than the NmE due to the solar EUV. This indicates the importance of including the contribution of precipitating electrons in IRI. The global equatorial boundary helps to improve the specification of the sub-auroral ionosphere trough in IRI. An IDL software package has been developed to interactively display the IRI parameters with assimilated NmE and EB. It can serve as an operational tool for space weather monitoring.     

Time and spatial variations of the foF2(night)/foF2(day) values

The ratio of daytime and nighttime values of the foF2 critical frequency is analyzed on the basis of the data of 35 ionospheric stations of the eastern hemisphere and 7 stations of the western hemisphere. The analysis confirms the results of Danilov [Danilov, A.D. Variations in the ratio of nighttime and daytime critical frequencies of the F2 layer. Geomagn. Aeronomy, 47 (6), 551–561 (in Russian) 2007; Danilov, A.D. Time and spatial variations in the ratio of nighttime and daytime critical frequencies of the F2 layer, J. Atmos. Sol-Terr. Phys. 70, 1201–1212, 2008.] that after about 1980, a systematic behavior of this ratio with time is observed: an increase with time (a positive trend) or a decrease with time (a negative trend). The sign of this trend is shown to be governed by the signs of the magnetic declination D and magnetic inclination I at the given ionospheric station. This fact makes it possible to assume that the above indicated trend is caused by long-term variations in the zonal component V_ny of the horizontal wind in the thermosphere, the latter component contributing into the vertical drift velocity W.     

Mid-latitude amplitude scintillation of GPS signals and GPS performance slips

Degradation of transionospheric radio signals and operation failures during ionospheric disturbances constitute a crucial factor of space weather influence on radio engineering satellite systems performance. We found that during the main phase of strong magnetic storms in 2000–2003 when the auroral oval expands into mid-latitudes, its southern boundary develops a region with intense small-scale electron density irregularities. Such irregularities may cause strong amplitude scintillations of GPS signals at both GPS operating frequencies. The another consequence of it was significant random GPS signal phase fluctuations, breaking-down of signal tracking, and sharp increasing of GPS positioning errors as a result.     

New field of application of the IRI modeling – Determination of ionosphere transfer characteristic for radio astronomical signals

We suggest a new field of application of IRI modeling – determination of ionosphere transfer characteristic (ITC) for radio astronomical signals (RAS). VHF and HF RAS are widely used for observations of the Sun and pulsars. It is necessary to take into account possible distortions of RAS in the Earth ionosphere. However, in contrast to modern navigation systems (GPS, GLONASS, GALILEO), where very accurate reconstruction of ionosphere parameters is a built-in function, in present-day radio astronomy a retrieve of ITC has not been appropriately worked out yet. It collides with increasing requirements to accuracy of the analysis of RAS amplitude profile and to the angular and polarizing resolution of radio telescopes of new generation. We have developed a method and software for calculation of the ionosphere measure of rotation (RM) and the measure of dispersion (DM). We used the ionosphere model IRI-2001, magnetic-field model IGRF-10 and values of ionosphere total electron content as deduced from GPS measurements. The obtained values of the ionosphere DM and RM were recalculated into characteristics of phase delay, Faraday amplitude modulation and polarization changes. We made calculations for different levels of geomagnetic activity and for different angular position of radio sources as well.     

The high-latitude ionosphere structure on 22 March, 1979 magnetic storm from multi-satellite and ground-based observation

Variations in the high-latitude ionosphere structure during March 22, 1979 geomagnetic storm are examined. Electron density Ne and temperature Te from the Cosmos-900 satellite, NmF2, Ne and He⁺ from the ISS-b satellite, precipitation of soft electrons from the Intercosmos-19 satellite, and the global picture of the auroral electron precipitation from the DMSP, TIROS and P78 satellites are used. These multi-satellite databases allow us to investigate the storm-time variations in the locations of the following ionospheric structures: the day-time cusp, the equatorial boundary of the diffuse auroral precipitation (DPB), the main ionospheric trough (MIT), the day-time trough, the ring ionospheric trough (RIT) and the light ions trough (LIT). The variations in NmF2, Ne, He⁺ and Te in the high-latitude ionosphere for the different local time sectors are analyzed also. The features of the high-latitude ionospheric response to a strong magnetic storm are described.     

Anthropogenic trigger of substorms and energetic particles precipitations

The high-frequency (HF) emission in near-Earth space from various powerful transmitters (radio communications, radars, broadcasting, universal time and navigation stations, etc.) form an integral part of the modern world that it cannot do without. In particular, special-purpose research facilities equipped with powerful HF transmitters are used successfully for plasma experiments and local modification of the ionosphere. In this work, we are using the results of a complex space-ground experiment to show that exposure of the subauroral region to HF emission can not only cause local changes in the ionosphere, but can also trigger processes in the magnetosphere–ionosphere system that result in intensive substorm activity (precipitations of high-energy particles, aurorae, significant variations in the ionospheric parameters and, as a consequence, in radio propagation conditions).     

First results of operational ionospheric dynamics prediction for the Brazilian Space Weather program

It is shown the development and preliminary results of operational ionosphere dynamics prediction system for the Brazilian Space Weather program. The system is based on the Sheffield University Plasmasphere–Ionosphere Model (SUPIM), a physics-based model computer code describing the distribution of ionization within the Earth mid to equatorial latitude ionosphere and plasmasphere, during geomagnetically quiet periods. The model outputs are given in a 2-dimensional plane aligned with Earth magnetic field lines, with fixed magnetic longitude coordinate. The code was adapted to provide the output in geographical coordinates. It was made referring to the Earth’s magnetic field as an eccentric dipole, using the approximation based on International Geomagnetic Reference Field (IGRF-11). During the system operation, several simulation runs are performed at different longitudes. The original code would not be able to run all simulations serially in reasonable time. So, a parallel version for the code was developed for enhancing the performance. After preliminary tests, it was frequently observed code instability, when negative ion temperatures or concentrations prevented the code from continuing its processing. After a detailed analysis, it was verified that most of these problems occurred due to concentration estimation of simulation points located at high altitudes, typically over 4000 km of altitude. In order to force convergence, an artificial exponential decay for ion–neutral collisional frequency was used above mentioned altitudes. This approach shown no significant difference from original code output, but improved substantially the code stability. In order to make operational system even more stable, the initial altitude and initial ion concentration values used on exponential decay equation are changed when convergence is not achieved, within pre-defined values. When all code runs end, the longitude of every point is then compared with its original reference station longitude, and differences are compensated by changing the simulation point time slot, in a temporal adjustment optimization. Then, an approximate neighbor searching technique was developed to obtain the ion concentration values in a regularly spaced grid, using inverse distance weighting (IDW) interpolation. A 3D grid containing ion and electron concentrations is generated for every hour of simulated day. Its spatial resolution is 1° of latitude per 1° of longitude per 10 km of altitude. The vertical total electron content (VTEC) is calculated from the grid, and plotted in a geographic map. An important feature that was implemented in the system is the capacity of combining observational data and simulation outputs to obtain more appropriate initial conditions to the ionosphere prediction. Newtonian relaxation method was used for this data assimilation process, where ionosonde data from four different locations in South America was used to improve the system accuracy. The whole process runs every day and predicts the VTEC values for South America region with almost 24 h ahead.