An approximation of the height of the O+ - H+ transition level for use in IRI

A set of data including over 1400 transition level (TL) encounters of the OGO-6 satellite is used to construct analytical presentations of this altitude upon dip-latitude and longitude. The formula represents the TL variation between ± 50° dip-latitude and 20 h and 04 h local time interval. The data were obtained in the high solar activity period in 1969 and 1970.     

Analytical expression of O---H ion transition surface for use in IRI

A global surface of O⁺---H⁺ transition level is constructed, based on published data from OGO-6, Intercosmos-2, Alouette-1, ISS-b, and TAIYO satellites. This surface covers ±60° dipole latitude, all longitudes, two levels of solar activity, summer and winter solstices, and 00 and 12 hours local time. The surface is used as input data to a mathematical model which calculates transition levels in 5-dimensional space: sunspot number (R), month (M), local time (LT), dipole latitude (DL), and longitude (LONG). This model is based on a generalized multivariable polynomial, using a system of linearly independent functions. Model transition levels are compared with averaged data from AE-E and AE-C, as well as rocket measurements from Vertical-6 and Vertical-10. The obtained analytical expression can be directly used in IRI.     

Parallel,adaptive, multi-object trajectory integrator for space simulation applications

Computer simulation is a very helpful approach for improving results from space born experiments. Initial-value problems (IVPs) can be applied for modeling dynamics of different objects – artificial Earth satellites, charged particles in magnetic and electric fields, charged or non-charged dust particles, space debris. An ordinary differential equations systems (ODESs) integrator based on applying different order embedded Runge–Kutta–Fehlberg methods is developed. These methods enable evaluation of the local error. Instead of step-size control based on local error evaluation, an optimal integration method is selected. Integration while meeting the required local error proceeds with constant-sized steps. This optimal scheme selection reduces the amount of calculation needed for solving the IVPs. In addition, for an implementation on a multi core processor and parallelization based on threads application, we describe how to solve multiple systems of IVPs efficiently in parallel.     

Modeling medium-scale TEC structures,observed by Belgian GPS receivers network

GALOCAD project “Development of a Galileo Local Component for the nowcasting and forecasting of atmospheric disturbances affecting the integrity of high precision Galileo applications” aims to perform a detailed study on ionospheric small- and medium-scale structures and to assess the influence of these structures on the reliability of Galileo precise positioning applications. GPS-derived TEC (total electron content) is obtained from the Belgium Dense Network (BDN), consisting of 67 permanent GPS stations. An empirical 3-D model is developed for studying these ionospheric structures. The model, named LLT model, described temporal variations of TEC in latitude/longitude frame (46°, 52°)N and (−1°, 11°)E. The spatial variations of TEC are modeled by Tchebishev base functions, while the temporal variations are described by a trigonometric basis. To fit the model to the data, the observed area is divided into bins with (1° × 1°) geographic scale and 6 min on time axis. LLT model is made flexible, with varying number of coefficients along each axis. This allows different degree of smoothing, which is the key element of the present approach. Model runs with higher number of coefficients, capturing in details medium-scale TEC structures are subtracted from results obtained with smaller number of coefficients; the latter represent the background ionosphere. The residual structures are localized and followed as they travel across the observed area. In this way, the size, velocity, and direction of the irregular structures are obtained.     

Comparison of calculation models for determination of the mesopause temperature using SATI images

The Spectral Airglow Temperature Imager is an instrument for ground-based spectroscopic measurements of the night-glow atmosphere emissions. This instrument was developed specially for gravity wave investigation. The measured airglow spectra are matched to synthetic spectra calculated in advance for determination of the temperature in the mesopause region where the radiation maximum of some О₂ emissions is situated. The synthetic spectra are transformed into a format which corresponds to the measured spectra in order to be matched. This transformation is based on the known values of the refractive index and the central wavelength of the interference filter used. A substantial part of the processing algorithms of the SATI images is connected with determination of these two filter parameters. The results of the original and newly-proposed algorithms for filter parameter calculation and their importance for the final results for temperature determination on the basis of the О₂ (864–868 nm) emission measurements are presented.     

DEMETER and DMSP satellite observations of the disturbed H/O ratio caused by Earth’s seismic activity in the Sumatra area during December 2004

In the present paper, plasma probe data taken from DEMETER and DMSP-F15 satellites were used to study the ion density and temperature disturbances in the morning topside ionosphere, caused by seismic activity at low latitudes. French DEMETER (Detection of Electro-Magnetic Emissions Transmitted from Earthquake Regions) micro-satellite mission had been especially designed to provide global scale observations in the topside ionosphere over seismically active regions. Onboard the DEMETER satellite, the thermal plasma instrument called “Instrument Analyser de Plasma” (IAP) provides ion mass and densities, ion temperature, three component ion drift and ion density irregularities measurements. As a part of “Defense Meteorological Satellite Program”, DMSP-F15 satellite is on orbit operation since 1999. It provides ionospheric plasma diagnostics by means of the “Special Sensor-Ion, Electron and Scintillations” (SSIES-2) instrument. We examined few examples of possible seismic effects in the equatorial ionosphere, probably associated with seismic activity during December month in the area of Sumatra Island, including main shock of giant Sumatra event. It is found that the localized topside ionospheric disturbances appear close to the epicenters of certain earthquakes in the Sumatra region. In two cases, ion H⁺/O⁺ ratio rises more than one hour before the main shock, due to the O⁺ density decrease at the winter side of the geomagnetic equator, with longitudinally closest location to the epicenter of the earthquakes. These anomalous depletions in O⁺ density do exist in all cases of SSIES-2 data. Particularly for Sumatra main event, more than one hour after the main shock, we observe large-scale depletion in O⁺ density northward of the geomagnetic equator at winter side hemisphere. Associated with O⁺ depletion, ion temperature latitudinal profile around the geomagnetic equator shows enhanced asymmetry with minimum at the summer side and maximum in positive Ti deviation from mean value at the winter side. This disturbance lasted for more than three hours, later in time observed at the same place by IAP/DEMETER.     

Comparison of subpeak electron density profiles deduced from ionograms with the International Reference Ionosphere (IRI)

A comparison is made between the subpeak electron density profiles, obtained at selected local hours by vertical ionospheric sounding at the ionospheric station at Sofia (42.6°N; 23.3°E) and the IRI profiles for spring, summer, winter and two levels of solar activity (R = 10 and 100). It is demonstrated that the ionospheric profiles above Sofia are in rather good agreement with the values computed with IRI.     

Reconstruction of topside density profile by using the topside sounder model profiler and digisonde data

To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (H_s), O⁺–H⁺ transition height (H_T), and their ratio Rt = H_s/H_T, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.     

Comparison of topside ionosphere scale height modeled by the Topside Sounder Model and incoherent scatter radar ionospheric model

The topside ionosphere scale height extracted from two empirical models are compared in the paper. The Topside Sounder Model (TSM) provides directly the scale height (H_T), while the incoherent scatter radar ionospheric model (ISRIM) provides electron density profiles and its scale height (H_R) is determined by the lowest gradient in the topside part of the profile. H_T and H_R are presented for 7 ISR locations along with their dependences on season, local time, solar flux F10.7, and geomagnetic index ap. Comparison reveals that H_T values are systematically lower than respective H_R values as the average offset for all 7 stations is 55 km. For the midlatitude stations Arecibo, Shigaraki, and Millstone Hill this difference is reduced to 43 km. The range of variations of H_R is much larger than that of H_T, as the H_T range overlaps the lower part of the H_R range. Dependences on ap, DoY and LT are much stronger in the ISRIM than in TSM. This results in much larger values of H_R at higher ap. Diurnal amplitude of H_R is much larger than that of H_T, with large maximum of H_R at night. The present comparison yields the conclusion that the ISR measurements provide steeper topside Ne profiles than that provided by the topside sounders.     

Seasonal dependence of the ion density irregularities at the magnetic equator from OGO-6 retarding potential analyzer data

The morphology of the nighttime F-region irregularities is a subject of several satellite in-situ and ground-based optical and radar measurements. In this paper, we use the RMS variations of the total ion current as measured by retarding potential analyzer (RPA) aboard the OGO-6 satellite. More than 600 orbits were used for a statistical examination of the occurrence frequency of the irregularities at night (between 19 and 06 h LT) over the magnetic equator for spring and summer conditions in 1970. We discuss a possible statistical relationship between these OGO-6 irregularity data and thunderstorm activity data from OSO-5 for the same period.