Observation of Terrestrial gamma-ray flashes in the RELEC space experiment on the <Emphasis Type="Italic">Vernov</Emphasis> satellite

The RELEС scientific payload of the Vernov satellite launched on July 8, 2014 includes the DRGE spectrometer of gamma-rays and electrons. This instrument comprises a set of scintillator phoswich-detectors, including four identical X-ray and gamma-ray detector with an energy range of 10 kev to 3 MeV with a total area of ~500 cm² directed to the atmosphere, as well as an electron spectrometer containing three mutually orthogonal detector units with a geometric factor of ~2 cm² sr. The aim of a space experiment with the DRGE instrument is the study of fast phenomena, in particular Terrestrial gamma-ray flashes (TGF) and magnetospheric electron precipitation. In this regard, the instrument provides the transmission of both monitoring data with a time resolution of 1 s, and data in the event-by-event mode, with a recording of the time of detection of each gamma quantum or electron to an accuracy of ~15 μs. This makes it possible to not only conduct a detailed analysis of the variability in the gamma-ray range, but also compare the time profiles with the results of measurements with other RELEC instruments (the detector of optical and ultraviolet flares, radio-frequency and low-frequency analyzers of electromagnetic field parameters), as well as with the data of ground-based facility for thunderstorm activity. This paper presents the first catalog of Terrestrial gamma-ray flashes. The criterion for selecting flashes required in order to detect no less than 5 hard quanta in 1 ms by at least two independent detectors. The TGFs included in the catalog have a typical duration of ~400 μs, during which 10–40 gamma-ray quanta were detected. The time profiles, spectral parameters, and geographic position, as well as a result of a comparison with the output data of other Vernov instruments, are presented for each of candidates. The candidate for Terrestrial gamma-ray flashes detected in the near-polar region over Antarctica is discussed.     

Estimate of the NO concentration in the auroral region based on emission intensities of 391.4, 557.7, and 630.0 nm

Using numerical modeling, the influence of the NO concentration on the intensity of 557.7 nm emission in aurora caused by electron precipitation has been studied. It has been shown that the O₂ ⁺ NO reaction, which reduces the contribution of the dissociative recombination of the O₂ ⁺ ion into the formation of the ¹S state of atomic oxygen, is the main channel of suppression of the intensity of the emission at 557.7 nm. A method of estimating the NO concentration in the aurora based on the data of photometric measurements of emissions at 391.4, 557.7, and 630.0 nm has been proposed. The method has been tested using the data of simultaneous rocket measurements of emissions at 391.4, 557.7, and 630.0 nm and the NO content in aurora. A good agreement of estimates of the NO concentrations performed by the method to the results of direct measurements has been obtained.     

Estimating the accuracy of the technique of reconstructing the rotational motion of a satellite based on the measurements of its angular velocity and the magnetic field of the Earth

The paper has studied the accuracy of the technique that allows the rotational motion of the Earth artificial satellites (AES) to be reconstructed based on the data of onboard measurements of angular velocity vectors and the strength of the Earth magnetic field (EMF). The technique is based on kinematic equations of the rotational motion of a rigid body. Both types of measurement data collected over some time interval have been processed jointly. The angular velocity measurements have been approximated using convenient formulas, which are substituted into the kinematic differential equations for the quaternion that specifies the transition from the body-fixed coordinate system of a satellite to the inertial coordinate system. Thus obtained equations represent a kinematic model of the rotational motion of a satellite. The solution of these equations, which approximate real motion, has been found by the least-square method from the condition of best fitting between the data of measurements of the EMF strength vector and its calculated values. The accuracy of the technique has been estimated by processing the data obtained from the board of the service module of the International Space Station (ISS). The reconstruction of station motion using the aforementioned technique has been compared with the telemetry data on the actual motion of the station. The technique has allowed us to reconstruct the station motion in the orbital orientation mode with a maximum error less than 0.6° and the turns with a maximal error of less than 1.2°.     

Dynamics of electron fluxes in the slot between radiation belts in November–December 2014 according to data of the <Emphasis Type="Italic">Vernov</Emphasis> satellite

The variations in the spatial structure and time in electron fluxes with E = 235–300 keV in the slot region (2 < L < 3) between the radiation belts in the period of November 1, 2014 through December 8, 2014 during weak and moderate geomagnetic disturbances (Kp < 4, Dst >–60 nT) are analyzed based on the data of the RELEC complex on board the Vernov satellite (the height and inclination of the orbit are from 640 to 830 km and 98.4°, respectively). Irregular increases in the fluxes of such electrons and formation of a local maximum at L ~ 2.2–3.0 were observed. It has been shown that the intensity of this maximum is inversely proportional to the L value and grows with an increase in the geomagnetic activity level. New features discovered for the first time in the dynamics of radiation belt electrons manifest in the variations in the local structure and dynamics of fluxes of subrelativistic electrons in the slot region.     

Hardware reconfiguration algorithm in multiprocessor systems of integrated modular avionics

Reconfiguration of multiprocessor systems makes it possible to improve their failure-resistance that is especially important for the integrated modular avionics systems. The algorithm considered in this paper allows minimizing the reservation and providing the better safety level and more effective flight completion or even its further execution in the case of failures of airborne equipment.     

Comparison of vertical E × B drift velocities and ground-based magnetometer observations of DELTA H in the low latitude under geomagnetically disturbed conditions

In the present work, we analyzed the daytime vertical E × B drift velocities obtained from Jicamarca Unattended Long-term Ionosphere Atmosphere (JULIA) radar and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca, and Piura in Peru for 22 geomagnetically disturbed events in which either SC has occurred or Dst_max < ?50 nT during the period 2006–2011. The ΔH component of geomagnetic field is measured as the differences in the magnitudes of horizontal H component between magnetometer placed directly on the magnetic equator and one displaced 6–9° away. It will provide a direct measure of the daytime electrojet current, due to the eastward electric field. This will in turn gives the magnitude of vertical E × B drift velocity in the F region. A positive correlation exists between peak values of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of the events. It was observed that 45% of the events have daytime vertical E × B drift velocity peak in the magnitude range 10–20 m/s and 20–30 m/s and 20% have peak ΔH in the magnitude range 50–60 nT and 80–90 nT. It was observed that the time of occurrence of the peak value of both the vertical E × B drift velocity and the ΔH have a maximum (40%) probability in the same time range 11:00–13:00 LT. We also investigated the correlation between E × B drift velocity and D_st index and the correlation between delta H and D_st index. A strong positive correlation is found between E × B drift and D_st index as well as between delta H and D_st Index. Three different techniques of data analysis – linear, polynomial (order 2), and polynomial (order 3) regression analysis were considered. The regression parameters in all the three cases were calculated using the Least Square Method (LSM), using the daytime vertical E × B drift velocity and ΔH. A formula was developed which indicates the relationship between daytime vertical E × B drift velocity and ΔH, for the disturbed periods. The E × B drift velocity was then evaluated using the formulae thus found for the three regression analysis and validated for the ‘disturbed periods’ of 3 selected events. The E × B drift velocities estimated by the three regression analysis have a fairly good agreement with JULIA radar observed values under different seasons and solar activity conditions. Root Mean Square (RMS) errors calculated for each case suggest that polynomial (order 3) regression analysis provides a better agreement with the observations from among the three.     

Case study of inclined sporadic E layers in the Earth’s ionosphere observed by CHAMP/GPS radio occultations: Coupling between the tilted plasma layers and internal waves

We have used the radio occultation (RO) satellite data CHAMP/GPS (Challenging Minisatellite Payload/Global Positioning System) for studying the ionosphere of the Earth. A method for deriving the parameters of ionospheric structures is based upon an analysis of the RO signal variations in the phase path and intensity. This method allows one to estimate the spatial displacement of a plasma layer with respect to the ray perigee, and to determine the layer inclination and height correction values. In this paper, we focus on the case study of inclined sporadic E (E_s) layers in the high-latitude ionosphere based on available CHAMP RO data. Assuming that the internal gravity waves (IGWs) with the phase-fronts parallel to the ionization layer surfaces are responsible for the tilt angles of sporadic plasma layers, we have developed a new technique for determining the parameters of IGWs linked with the inclined E_s structures. A small-scale internal wave may be modulating initially horizontal E_s layer in height and causing a direction of the plasma density gradient to be rotated and aligned with that of the wave propagation vector k. The results of determination of the intrinsic wave frequency and period, vertical and horizontal wavelengths, intrinsic vertical and horizontal phase speeds, and other characteristics of IGWs under study are presented and discussed.     

Relative equilibria in quasi-homogeneous planar three body problems

In this paper we find the families of relative equilibria for the three body problem in the plane, when the interaction between the bodies is given by a quasi-homogeneous potential. The number of the relative equilibria depends on the values of the masses and on the size of the system, measured by the moment of inertia.     

Comparison of GPS derived TEC with the TEC predicted by IRI 2012 model in the southern Equatorial Ionization Anomaly crest within the Eastern Africa region

We have compared the TEC obtained from the IRI-2012 model with the GPS derived TEC data recorded within southern crest of the EIA in the Eastern Africa region using the monthly means of the 5 international quiet days for equinoxes and solstices months for the period of 2012 – 2013. GPS-derived TEC data have been obtained from the Africa array and IGS network of ground based dual-frequency GPS receivers from four stations (Kigali (1.95°S, 30.09°E; Geom. Lat. 11.63°S), Malindi (2.99°S, 40.19°E; Geom. Lat. 12.42°S), Mbarara (0.60°S, 30.74°E; Geom. Lat. 10.22°S) and Nairobi (1.22°S, 36.89°E; Geom. Lat. 10.69°S)) located within the EIA crest in this region. All the three options for topside Ne of IRI-2012 model and ABT-2009 for bottomside thickness have been used to compute the IRI TEC. Also URSI coefficients were considered in this study. These results are compared with the TEC estimated from GPS measurements. Correlation Coefficients between the two sets of data, the Root-Mean Square Errors (RMSE) of the IRI-TEC from the GPS-TEC, and the percentage RMSE of the IRI-TEC from the GPS-TEC have been computed. Our general results show that IRI-2012 model with all three options overestimates the GPS-TEC for all seasons and at all stations, and IRI-2001 overestimates GPS-TEC more compared with other options. IRI-Neq and IRI-01-corr are closely matching in most of the time. The observation also shows that, GPS TEC are underestimated by TEC from IRI model during noon hours, especially during equinoctial months. Further, GPS-TEC values and IRI-TEC values using all the three topside Ne options show very good correlation (above 0.8). On the other hand, the TEC using IRI-Neq and IRI-01- corr had smaller deviations from the GPS-TEC compared to the IRI-2001.     

Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.