Solution of the flyby problem for large space debris at sun-synchronous orbits

the paper considers the flyby problem related to large space debris (LSD) objects at low earth orbits. The data on the overall dimensions of known last and upper stages of launch vehicles makes it possible to single out five compact groups of such objects from the NORAD catalog in the 500–2000 km altitude interval. The orbits of objects of each group have approximately the same inclinations. The features of the mutual distribution of the orbital planes of LSD objects in the group are shown in a portrait of the evolution of deviations of the right ascension of ascending nodes (RAAN). In the case of the first three groups (inclinations of 71°, 74°, and 81°), the straight lines of relative RAAN deviations of object orbits barely intersect each other. The fourth (83°) and fifth (97°–100°) LSD groups include a considerable number of objects whose orbits are described by straight lines (diagonals), which intersect other lines many times. The use of diagonals makes it possible to significantly reduce the temporal and total characteristic velocity expenditures required for object flybys, but it complicates determination of the flyby sequence. Diagonal solutions can be obtained using elements of graph theory. A solution to the flyby problem is presented for the case of group 5, formed of LSD objects at sun-synchronous orbits.     

Families of periodic solutions in three-dimensional restricted three-body problem

Planar orbits of three-dimensional restricted circular three-body problem are considered as a special case of three-dimensional orbits, and the second-order monodromy matrices M (in coordinate z and velocity v _z ) are calculated for them. Semi-trace s of matrix M determines vertical stability of an orbit. If |s| ≤ 1, then transformation of the subspace (z, v _z ) in the neighborhood of solution for the period is reduced to deformation and a rotation through angle φ, cosφ = s. If the angle ? can be rationally expressed through 2π,φ = 2π·p/q, where p and q are integer, then a planar orbit generates the families of three-dimensional periodic solutions that have a period larger by a factor of q (second kind Poincareé periodic solutions). Directions of continuation in the subspace (z, v _z ) are determined by matrix M. If |s| < 1, we have two new families, while only one exists at resonances 1: 1 (s = 1) and 2: 1 (s = ?1). In the course of motion along the family of three-dimensional periodic solutions, a transition is possible from one family of planar solutions to another one, sometimes previously unknown family of planar solutions.     

Forecasting the velocity of quasi-stationary solar wind and the intensity of geomagnetic disturbances produced by it

A brief review is given of contemporary approaches to solving the problem of medium-term forecast of the velocity of quasi-stationary solar wind (SW) and of the intensity of geomagnetic disturbances caused by it. At the present time, two promising models of calculating the velocity of quasi-stationary SW at the Earth’s orbit are realized. One model is the semi-empirical model of Wang-Sheeley-Arge (WSA) which allows one to calculate the dependence V(t) of SW velocity at the Earth’s orbit using measured values of the photospheric magnetic field. This model is based on calculation of the local divergence f _S of magnetic field lines. The second model is semi-empirical model by Eselevich-Fainshtein-Rudenko (EFR). It is based on calculation in a potential approximation of the area of foot points on the solar surface of open magnetic tubes (sources of fast quasistationary SW). The new Bd-technology is used in these calculations, allowing one to calculate instantaneous distributions of the magnetic field above the entire visible surface of the Sun. Using predicted V(t) profiles, one can in EFR model calculate also the intensity of geomagnetic disturbances caused by quasi-stationary SW. This intensity is expressed through the K _p index. In this paper the EFR model is discussed in detail. Some examples of epignosis and real forecast of V(t) and K _p (t) are discussed. A comparison of the results of applying these two models for the SW velocity forecasting is presented.     

Large-scale irregularities of the winter polar topside ionosphere according to data from <Emphasis Type="Italic">Swarm</Emphasis> satellites

An analysis of the electron density measurements (Ne) along the flyby trajectories over the high-latitude region of the Northern Hemisphere under winter conditions in 2014 and 2016 has shown that the main large-scale structure observed by Swarm satellites is the tongue of ionization (TOI). At the maximum of the solar cycle (F_10.7 = 160), the average value of Ne in the TOI region at an altitude of 500 km was 8 × 10⁴ cm^–3. Two years later, at F_10.7 = 100, Ne ~ 5 × 10⁴ cm^–3 and Ne ~2.5 × 10⁴ cm^–3 were observed at altitudes of 470 and 530 km, respectively. During the dominance of the azimuthal component of the interplanetary magnetic field, the TOI has been observed mainly on the dawn or dusk side depending on the sign of B _y . Simultaneous observations of the convective plasma drift velocity in the polar cap show the transpolar flow drift to the dawn (B_y < 0) or dusk side (B _y < 0). Observations and numerical simulation of the Ne distribution have confirmed the significant role of the electric field of the magnetospheric convection in the generation of large-scale irregularities in the polar ionosphere.     

Fluxes of Quasi-trapped Electrons with Energies >0.08 MeV in the Near-Earth Space on Drift Shells <Emphasis Type="Italic">L</Emphasis> < 2

We study the characteristics of fluxes of electrons with energy >80 keV in the near-Earth space regions corresponding to the drift shells L = 1.7, 1.4, and 1.1 observed during the entire period of the GRIF experiment onboard the Spectr module of the Mir orbital station from October 1995 to June 1997. The obtained geographic maps of the distribution of electron fluxes at the height of the station flight (400 km) and, also, the estimates of the spectra indicate that the South-Atlantic Anomaly provides for a mechanism of stable replenishment for shells with L < 1.5. The mechanism of stable replenishment of shells with L < 1.5 may be due to the scattering, in the residual atmosphere, of electrons from the inner radiation belt precipitating into the region of the South-Atlantic Anomaly.     

Statistical investigation of Heliospheric conditions resulting in magnetic storms: 2

This work is a continuation of investigation [1] of the behavior of the solar wind’s and interplanetary magnetic field’s parameters near the onset of geomagnetic storms for various types of solar wind streams. The data of the OMNI base for the 1976–2000 period are used in the analysis. The types of solar wind streams were determined, and the times of beginning (onsets) of magnetic storms were distributed in solar wind types as follows: CIR (121 storms), Sheath (22 storms), MC (113 storms), and “uncertain type” (367 storms). The growth of variations (hourly standard deviations) of the density and IMF magnitude was observed 5–10 hours before the onset only in the Sheath. For the CIR-, Sheath-and MC-induced storms the dependence between the minimum of the IMF B _z-component and the minimum of the D _st -index, as well as the dependence between the electric field E _y of solar wind and the minimum of the D _st -index are steeper than those for the “uncertain” solar wind type. The steepest D _st vs. B _z dependence is observed in the Sheath, and the steepest D _st vs. E _y dependence is observed in the MC.     

Short increases of the fluxes of electrons with energies > 0.08 MeV at low-latitude (<Emphasis Type="Italic">L</Emphasis> < 2) regions of near-Earth space

We present the characteristics of short (duration less than 1 min) increases of the counting rate of electrons with energies >0.08 MeV observed in low-latitude (L < 2.0) regions of near-Earth space in the course of the GRIF experiment on the Spektr module of the Mir orbital station. The measurements were carried out using a set of instruments including X-ray and gamma-ray spectrometers, as well as detectors of electrons, protons, and nuclei with large and small geometrical factors, which allowed one to detect the fluxes of charged particles both in the region of the Earth’s radiation belts and in regions close to the geomagnetic equator. As a result of more than 1.5 years of observation, it is demonstrated that short increases in the intensity of electrons of subrelativistic energies are detected not only in the regions of the near-Earth space known as “precipitation zones” (1.7 < L < 2.5), but in high-latitude regions (up to the geomagnetic equator, L < 1.1) as well. Two types of increases of the electron counting rate are found: either fairly regular increases repeating on successive orbits or increases local in time. The latter type of increases can be caused by a short enhancement of electron flux on a given drift shell. The results of our measurements have shown that the duration of the detected increases in intensity can be rather short, as little as 20–30 s. Therefore, in the case of large amplitudes, such increases of the counting rate of electrons can imitate astrophysical events of the type of cosmic gamma-ray bursts in the detectors of hard X-ray and gamma radiation.     

Reconstruction of spacecraft rotational motion using a Kalman filter

Quasi-static microaccelerations of four satellites of the Foton series (nos. 11, 12, M-2, M-3) were monitored as follows. First, according to measurements of onboard sensors obtained in a certain time interval, spacecraft rotational motion was reconstructed in this interval. Then, along the found motion, microacceleration at a given onboard point was calculated according to the known formula as a function of time. The motion was reconstructed by the least squares method using the solutions to the equations of satellite rotational motion. The time intervals in which these equations make reconstruction possible were from one to five orbital revolutions. This length is increased with the modulus of the satellite angular velocity. To get an idea on microaccelerations and satellite motion during an entire flight, the motion was reconstructed in several tens of such intervals. This paper proposes a method for motion reconstruction suitable for an interval of arbitrary length. The method is based on the Kalman filter. We preliminary describe a new version of the method for reconstructing uncontrolled satellite rotational motion from magnetic measurements using the least squares method, which is essentially used to construct the Kalman filter. The results of comparison of both methods are presented using the data obtained on a flight of the Foton M-3.     

The use of the Earth’s gravitational potential for interplanetary flights

An efficient scheme of the use of the Earth’s gravity in interplanetary flights is suggested, which opens up new opportunities for exploration of the solar system. The scheme of the gravitational maneuver allows one to considerably reduce the spacecraft mass consumption for a flight and the time of flight. An algorithm of the gravitational maneuver is suggested that takes into account the restriction on the altitude of a planet flyby. Estimates are made of transport capabilities for delivery of a spacecraft to the orbits of Jupiter, Saturn, and Uranus. The spacecraft is based on a middle-class carrier launcher of the Soyuz type and includes chemical and electric plasma jet engines of the SPD-140 type, which use solar energy.     

Solar activity index for long-term ionospheric forecasts

Based on the comparison of solar activity indices (annual average values of the relative number of sunspots Rz₁₂ and solar radio emission flux at a wavelength of 10.7 cm F₁₂) with the ionospheric index of solar activity IG₁₂ for 1954–2013, we have found that the index F₁₂ is a more accurate (than Rz₁₂) indicator of solar activity for the long-term forecast of foF2 (the critical frequency of the F2-layer). This advantage of the F₁₂ index becomes especially significant after 2000 if the specific features of extreme ultraviolet radiation of the Sun are additionally taken into account in the minima of solar cycles, using an appropriate correction to F₁₂. Qualitative arguments are given in favor of the use of F₁₂ for the long-term forecast of both foF2 and other ionospheric parameters.     

Observations of electromagnetic emissions inside the Earth’s plasmasphere from the <Emphasis Type="Italic">Interball-1</Emphasis> satellite

The spectrum analyzer AKR-X onboard the Interball-1 satellite at the beginning (August–October 1995) and at the end (August–October 2000) of satellite operation in perigees of its orbital motion recorded and analyzed electromagnetic emissions of the inner regions of the Earth’s plasmasphere in the frequency band 100–1500 kHz at distances of 1.1–1.8 R _E. The observations have shown that the electromagnetic modes (the Z and LO modes escaping the magnetosphere) which are formed at the altitudes 600–4000 km are associated with the subauroral nonthermal continuum and with the recently discovered kilometric continuum. There are noticeable differences in the spectral character of these emissions during the minimum (1996) and maximum (2000) solar activity, when, as a rule, the LO mode escaping the plasmaphere and the continua are not present.     

Observations of corotating solar wind structures at radio sounding by signals of the <Emphasis Type="Italic">Rosetta</Emphasis> and <Emphasis Type="Italic">Mars Express</Emphasis> spacecraft

In the implementation of the space projects Rosetta and Mars Express, a large-scale series of experiments has been carried out on radio sounding circumsolar plasma by decimeter (S-band) and centimeter (X-band) signals of the Rosetta comet probe (from October 3 to October 31, 2010) and the Mars Express satellite of Mars (from December 25, 2010 to March 27, 2011). It was found that in the phase of ingress the spacecraft behind the Sun, the intensity of the frequency fluctuations increases in accordance with a power function whose argument is the solar offset distance of radio ray path, and when the spacecraft is removed from the Sun (the egress phase), frequency fluctuations are reduced. Periodic strong increases in the fluctuation level, exceeding by a factor of 3–12 the background values of this value determined by the regular radial dependences, are imposed on the regular dependences. It was found that increasing the fluctuations of radio waves alternates with the periodicity m × T or n × T, where m = 1/2, n = 1, аnd T is the synodic period of the Sun’s rotation (T ≈ 27 days). It was shown that the corotating structures associated with the interaction regions of different speed fluxes are formed in the area of solar wind acceleration and at distances of 6–20 solar radii already have a quasi-stationary character.     

A method of calculation of vertical cutoff rigidity in the geomagnetic field

On the basis of generalization of the results of extensive trajectory calculations for trial charged particles moving in the geomagnetic field the method of calculation of effective vertical cutoff rigidity, taking into account the values of K _p -index and local time, is developed. The IGRF and Tsyganenko-89 models are used for the geomagnetic field. A comparison of the results of model simulations with the experimental data on penetration of charged particles into near-Earth space is made, and penetration functions for typical spacecraft orbits are calculated.     

Heating and acceleration of charged particles during magnetic dipolarizations

In this paper, we analyzed the thermal and energy characteristics of the plasma components observed during the magnetic dipolarizations in the near tail by the Cluster satellites. It was previously found that the first dipolarization the ratio of proton and electron temperatures (T _p/T _e) was ~6–7. At the time of the observation of the first dipolarization front T _p/T _e decreases by up to ~3–4. The minimum value T _p/T _e (~2.0) is observed behind the front during the turbulent dipolarization phase. Decreases in T _p/T _e observed at this time are associated with an increase in T _e, whereas the proton temperature either decreases or remains unchanged. Decreases of the value T _p/T _e during the magnetic dipolarizations coincide with increase in wave activity in the wide frequency band up to electron gyrofrequency f _ce. High-frequency modes can resonantly interact with electrons causing their heating. The acceleration of ions with different masses up to energies of several hundred kiloelectron-volts is also observed during dipolarizations. In this case, the index of the energy spectrum decreases (a fraction of energetic ions increases) during the enhancement of low-frequency electromagnetic fluctuations at frequencies that correspond to the gyrofrequency of this ion component. Thus, we can conclude that the processes of the interaction between waves and particles play an important role in increasing the energy of plasma particles during magnetic dipolarizations.     

Peculiarities of long-wave radio bursts from solar flares preceding strong geomagnetic storms

Radio bursts in the frequency range of 100–1500 kHz, recorded in 1997–2000 on the INTERBALL-1 satellite during the solar flares preceding the strong geomagnetic storms with D _st < ?100 nT, are analyzed in this paper. The observed long-wave III-type radio bursts of solar origin at frequencies of 1460 and 780 kHz were characterized by large values of the flux S _f = 10^?15 ?10^?17 W/m² Hz and duration longer than 10 min. The rapid frequency drift of a modulated radio burst continued up to a frequency of 250 kHz, which testified that the exciting agent (a beam of energetic electrons) propagated from the Sun to the Earth. All such flares were characterized by the appearance of halo coronal mass ejections, observed by the LASCO/SOHO, and by the presence of a southward B_z-component of the IMF, measured on the ACE and WIND spacecraft. In addition, shortly after radio bursts, the INTERBALL-1 satellite has recorded the fluxes of energetic electrons with E > 40 keV.     

Diurnal and longitudinal variations in the earth’s ionosphere in the period of solstice in conditions of a deep minimum of solar activity

The results of studies of longitudinal and LT variations in parameters of the ionosphere–plasmasphere system, obtained using the Global Self-Consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP), assimilation ionospheric model IRI Real-Time Assimilation Mapping (IRTAM), and satellite and ground-based observational data are presented in the paper. The study of the main morphological features of longitudinal and LT variations in the critical frequency of the ionospheric F2 layer (foF2) and total electron content (TEC) depending on latitude in the winter solstice during a solar-activity minimum (December 22, 2009) is carried out. It is shown that the variations in foF2 and TEC, on the whole, are identical, and so mutually substitutable, while creating empirical models of these parameters in quiet geomagnetic conditions. The longitudinal and LT variations in both foF2 and TEC are within an order of magnitude everywhere except for the equator anomaly region, where LT variation is larger by an order of magnitude than longitudinal variation. According to the results of the study, in the American longitudinal sector at all latitudes of the Southern (summer) Hemisphere, maxima of foF2 and TEC are formed. The near-equatorial and high-latitudinal maxima are separated out from these. The estimate of the contribution into the longitudinal variation in foF2 and TEC for various local time sectors and at various latitudes has been obtained for the first time. In the Southern (summer) Hemisphere, longitudinal variation in foF2 and TEC is formed in the nighttime.     

Study of notches in the Earth’s plasmasphere based on data of the <Emphasis Type="Italic">MAGION-5</Emphasis> satellite

Depleted narrow (localized in longitude) regions (field tubes) in the plasmasphere, recently discovered in He⁺ radiation measurements on the IMAGE spacecraft, were first directly observed by the Magion-5 satellite. The low-density regions (notches) occupy <～ 10–30° in longitude and extend from L ～ 2–3 to the plasmasphere boundary in neighboring plasmasphere regions with larger densities. The Magion-5 data give evidence that in the low-density regions temperature is enhanced as compared to the neighboring denser plasmasphere regions. Formation of notches in the plasmasphere is, apparently, associated with AE intensification during weak magnetic storms, while strong magnetic storms usually result in the overall reduction of plasmasphere dimensions. However, even a strong magnetic storm on April 6–7, 2000 (max K _p = 9-and min D _st ～ ?290 nT), but accompanied by an isolated AE impulse, resulted in a density decrease only in the longitudinally limited post-midnight sector of the plasmasphere.     

Turbulent fluctuations of plasma and magnetic field parameters in the magnetosheath and the low-latitude boundary layer formation: Multisatellite observations on March 2, 1996

The results of simultaneous analysis of plasma and magnetic field characteristics measured on the INTERBALL/Tail Probe, WIND and Geotail satellites on March 2, 1996, are presented. During these observations the INTERBALL/Tail Probe crossed the low-latitude boundary layer, and the WIND and Geotail satellites measured the solar wind’s and magnetosheath’s parameters, respectively. The plasma and magnetic field characteristics in these regions have been compared. The data of the Corall, Electron, and MIF instruments on the INTERBALL/Tail Probe satellite are analyzed. Fluctuations of the magnetic field components and plasma velocity in the solar wind and magnetosheath, measured onboard the WIND and Geotail satellites, are compared. The causes resulting in appearance of plasma jet flows in the low-latitude boundary layer are analyzed. The amplitude of magnetic field fluctuations in the magnetosheath for a studied magnetosphere boundary crossing is shown to exceed the magnetic field value below the magnetopause near the cusp. The possibility of local violation of pressure balance on the magnetopause is discussed, as well as penetration of magnetosheath plasma into the magnetosphere, as a result of magnetic field and plasma flux fluctuations in the magnetosheath.     

Sporadic structures and small-scale irregularity in the nighttime polar ionosphere in the period of high solar activity according to the data of radio occultation measurements on satellite-to-satellite paths

Results of the analysis of 327 sessions of radio occultation on satellite-to-satellite paths are presented. The data are taken in the nighttime polar ionosphere in the regions with latitudes of 67°–88°, and in the period of high solar activity from October 26, 2003 to November 9, 2003. Typical ionospheric changes in the amplitude and phase of decimeter radio waves on paths GPS satellites-CHAMP satellite are presented. It is demonstrated that these data make it possible to determine characteristics of the sporadic E _s structures in the lower ionosphere at heights of 75–120 km. Histograms of distribution of the lower and upper boundaries, thickness, and intensity of the E _s structures are presented. Dispersion and spectra of amplitude fluctuations of decimeter radio waves caused by small-scale irregularity of the ionospheric plasma are analyzed. The relation of the polar E _s structures and intensity of small-scale plasma irregularity to various manifestations of solar activity is discussed. The efficiency of monitoring the ionospheric disturbances caused by shock waves of the solar wind by the radio occultation method on satellite-to-satellite paths is demonstrated.