Optimization of planning the experiments to be carried out onboard orbital stations

The technique and algorithms for optimization of planning the program of experiments carried out onboard an orbiting spacecraft are described taking into account the execution of service operations. A general approach to optimization of planning the experiments is used, developed for investigations onboard the Salyut and Mir space stations, and on the International Space Station (ISS). The approach is based on formalization of the problem in the form of an integer linear programming problem. In this approach, the spacecraft orbit is considered to be known, and the optimization of the planning of experiments is reduced to composing the optimum sequence of zones for the performance of experiments. The list of experiments, service operations, and tasks to be solved during the planning interval are assumed to be specified.     

Study of singular surfaces of systems of unsupported gyrodynes by the parameter continuation method

A system of two-degree-of-freedom force gyroscopes (gyrodynes) is considered to be used for spacecraft attitude control. Possible values of its total angular momentum form some finite region P in the frame of reference rigidly connected with the spacecraft. Near the boundary of this region and singular surfaces located inside it the control of the angular momentum is complicated or impossible. The program angular momentum of the gyrodynes, realizing the law of variation of the spacecraft orientation, should lie inside P and outside some singular surfaces, and due to this fact the boundary and internal singular surfaces should be studied. This work is dedicated to the numerical construction of region P and its internal singular surfaces by the method of parameter continuation. Using the results by E.N. Tokar’ we formalize sufficient conditions which in some cases allow one to determine the type of the singular surface. As an example, a system of six gyrodynes is considered, for which the regions of variations of the intrinsic angular momentum and singular surfaces are constructed. The possibilities of the system are demonstrated for the case when one gyrodyne fails.     

Estimation of dynamic characteristics of the <Emphasis Type="Italic">International Space Station</Emphasis> from measurements of microaccelerations

We describe the results of determining the mass of the International Space Station using the data of MAMS accelerometer taken during correction of the station orbit on August 20, 2004. The correction was made by approach and attitude control engines (ACE) of the Progress transporting spacecraft. The engines were preliminary calibrated in an autonomous flight using the onboard device for measuring apparent velocity increment. The method of calibration is described and its results are presented. The error in station mass determination is about 1%. The same data of MAMS and similar data obtained during the orbit correction on August 26, 2004 were used for the analysis of high-frequency vibrations of the station mainframe caused by operation of the ACE of Progress. Natural frequencies of the ACE are determined. They lie in the frequency band 0.024–0.11 Hz. ACE operation is demonstrated to result in a substantial increase of microaccelerations onboard the station in the frequency range 0–1 Hz. The frequencies are indicated at which disturbances increase by more than an order of magnitude. The study described was carried out as a part of the Tensor technological experiment.     

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.     

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.     

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.     

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.     

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.     

Disturbance of the upper atmosphere after crash launch of the <Emphasis Type="Italic">Proton</Emphasis> carrier rocket

The data of inclined sounding obtained on July 5, 1999 after explosion of the Proton carrier launcher in the upper atmosphere are discussed. A comparative analysis is performed of the detected disturbances with disturbances recorded during explosion of the Soyuz rocket in the middle atmosphere and with disturbances typical for the standard mode of the rocket flight.     

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.     

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.     

Sсheme of rendezvous mission to lunar orbital station by spacecraft launched from Earth

In recent years, great experience has been accumulated in manned flight astronautics for rendezvous in near-Earth orbit. During flights of Apollo spacecraft with crews that landed on the surface of the Moon, the problem of docking a landing module launched from the Moon’s surface with the Apollo spacecraft’s command module in a circumlunar orbit was successfully solved. A return to the Moon declared by leading space agencies requires a scheme for rendezvous of a spacecraft launched from an earth-based cosmodromee with a lunar orbital station. This paper considers some ballistic schemes making it possible to solve this problem with minimum fuel expenditures.     

Study of microaccelerations onboard the <Emphasis Type="Italic">International Space Station</Emphasis> with the DAKON-M convection sensor

The results of experiments with the DAKON-M convection sensor onboard the Russian orbital segment of the International Space Station are described. A comparison of the sensor measurements with the results of calculation of the quasistatic microacceleration component at the point of installation is made. For this comparison we have used three measurement intervals of the experiments in 2009, during which spacecraft were docked with the station, undocked from it, and actuation of jet engines of the attitude control system took place. When calculating microacceleration, we use the measurement data of the low-frequency MAMS accelerometer, installed on the American segment, and the telemetry data on the ISS rotational motion. This information allowed one to convert the MAMS measurements to the point of installation of the DAKON-M convection sensor. A comparison of sensor measurements with calculated microaccelerations showed sufficiently accurate coincidence between the calculated and measured data.     

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.     

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.     

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.     

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.     

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.     

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.     

Study and comparison of the parameters of five hot flow anomalies at a bow shock front

Five hot flow anomalies (HFA) recorded by the Tail Probe of the INTERBALL satellite in 1996 are analyzed in present work. For the five chosen events the authors determined the characteristics of current sheets whose interaction with the bow shock front led to formation of an HFA, as well as the directions of external electric fields and the directions of motion of these HFAs over a shock front. The analysis of plasma convection in an HFA body is carried out; the average velocities of plasma motion in the HFA are determined in a coordinate system linked with the normal to a current layer and with the normal to the bow shock. According to the character of plasma convection in an HFA body, these five events may be divided into two types, which also differ in the direction of the motion over the front of the bow shock. In the first-type HFAs, the convection of plasma has a component directed from the intermediate region confirming its identification as a source of energy for the formation of an HFA. In the second-type HFAs, plasma motion from the intermediate region in leading and trailing parts is less expressed. This fact, as well as the great variation of peculiar velocities in the body of anomalies, allowed the assumption that second-type anomalies are nonstationary. Evidence is presented that the anomalies considered in the paper are bordered with shocks formed in solar wind passing a large-scale, decelerated body of heated plasma.