Motion of deformable planets and stability of their stationary rotations

A mechanical system consisting from N deformable spheres interacting according to the law of gravity is considered as a model of planetary system. Deformations of the viscoelastic spheres are described according to the model of the theory of elasticity of small deformations, the Kelvin-Voigt model of viscous forces, and occur under the action of gravitational fields and fields of centrifugal forces. Approximate equations describing motions of the centers of mass of the spheres and their rotations relative to the centers of mass are constructed by the method of separation of motions on the basis of solving quasistatic problems of the theory of viscoelasticity with allowance made for smallness of sphere deformations. Using the first integral of conservation of the angular momentum of the system relative to its center of mass, the expression for the changed potential energy is obtained with the use of the Routh method. An investigation of stationary rotations is carried out, and it is shown that all of them are unstable, if the number of planets is more than two.     

Radiation environment on the Mir orbital station during solar minimum.

The Mir station has been in a 51.65 degrees inclination orbit since March 1986. In March 1995, the first US astronaut flew on the Mir-18 mission and returned on the Space Shuttle in July 1995. Since then three additional US astronauts have stayed on orbit for up to 6 months. Since the return of the first US astronaut, both the Spektr and Priroda modules have docked with Mir station, altering the mass shielding distribution. Radiation measurements, including the direct comparison of US and Russian absorbed dose rates in the Base Block of the Mir station, were made during the Mir-18 and -19 missions. There is a significant variation of dose rates across the core module; the six locations sampled showed a variation of a factor of nearly two. A tissue equivalent proportional counter (TEPC) measured a total absorbed dose rate of 300 microGy/day, roughly equally divided between the rate due to trapped protons from the South Atlantic Anomaly (SAA) and galactic cosmic radiation (GCR). This dose rate is about a factor of two lower than the rate measured by the thinly shielded (0.5 g cm-2 of Al) operational ion chamber (R-16), and about 3/2 of the rate of the more heavily shielded (3.5 g cm-2 of Al) ion chamber. This is due to the differences in the mass shielding properties at the location of these detectors. A comparison of integral linear energy transfer (LET) spectra measured by TEPC and plastic nuclear track detectors (PNTDs) deployed side by side are in remarkable agreement in the LET region of 15-1000 keV/micrometer, where the PNTDs are fully efficient. The average quality factor, using the ICRP-26 definition, was 2.6, which is higher than normally used. There is excellent agreement between the measured GCR dose rate and model calculations, but this is not true for trapped protons. The measured Mir-18 crew skin dose equivalent rate was 1133 microSv/day. Using the skin dose rate and anatomical models, we have estimated the blood-forming organ (BFO) dose rate and the maximum stay time in orbit for International Space Station crew members.     

Corotating Interaction Regions in the Outer Heliosphere

We discuss the structure and evolution of CIRs and their successors in the outer heliosphere. These structures undergo significant evolution as they are convected to greater heliocentric distances. A progression of different types of structure are observed at increasing distance from the Sun. Similar structures are observed at similar heliocentric distance at different portions of the solar cycle. CIRs and their successors are associated with many important physical processes in the outer heliosphere. We discuss the relationship between these structures and recurrent phenomena such as cosmic ray variations, and review some of the associated theoretical models on the role of corotating structures and global merged interaction regions (GMIRs) in global cosmic ray modulation. We also discuss some outstanding questions related to the origin of non-dispersive quasi-periodic particle enhancements associated with CIRs and their successors in the outer heliosphere. This revised version was published online in August 2006 with corrections to the Cover Date.     

Modeling the pre-earthquake electrostatic effect on the F region ionosphere

This paper presents the results of modeling the ionospheric effect of the seismogenic electrostatic field (SEF) seen at the earth’s surface as a perturbation of the vertical atmospheric electrostatic field in the earthquake preparation zone. The SEF distribution at ionospheric altitudes is obtained as an analytical solution of the continuity equation for the electric current density. It is shown that at night, the horizontally large scale SEF can efficiently penetrate into the ionosphere and produce noticeable changes in the horizontal distribution of the F region electron density. The results suggest that the seismogenic electrostatic field could be a possible source for the ionospheric variations observed over Taiwan before the strong Chi Chi earthquake of September 21, 1999.     

Using altimetry and seafloor pressure data to estimate vertical deformation offshore: Vanuatu case study

Measuring ground deformation underwater is essential for understanding Earth processes at many scales. One important example is subduction zones, which can generate devastating earthquakes and tsunamis, and where the most important deformation signal related to plate locking is usually offshore. We present an improved method for making offshore vertical deformation measurements, that involve combining tide gauge and altimetry data. We present data from two offshore sites located on either side of the plate interface at the New Hebrides subduction zone, where the Australian plate subducts beneath the North Fiji basin. These two sites have been equipped with pressure gauges since 1999, to extend an on-land GPS network across the plate interface. The pressure series measured at both sites show that Wusi Bank, located on the over-riding plate, subsides by 11 ± 4 mm/yr with respect to Sabine Bank, which is located on the down-going plate. By combining water depths derived from the on-bottom pressure data with sea surface heights derived from altimetry data, we determine variations of seafloor heights in a global reference frame. Using altimetry data from TOPEX/Poseidon, Jason-1, Jason-2 and Envisat missions, we find that the vertical motion at Sabine Bank is close to zero and that Wusi Bank subsides by at least 3 mm/yr and probably at most 11 mm/yr.     

On the nature of coronal loops above the quiet sun network

The structure and dynamics of a box in a stellar corona can be modeled employing a 3D MHD model for different levels of magnetic activity. Depending on the magnetic flux through the surface the nature of the resulting coronal structures can be quite different. We investigate a model of an active region for two sunspots surrounded by magnetic field patches comparable in magnetic flux to the sunspots. The model results in emission from the model corona being concentrated in loop structures. In Gudiksen and Nordlund (2005) the loops seen in EUV and X-ray emission outline the magnetic field, following the general paradigm. However, in our model, where the magnetic field is far from a force-free state, the loops seen in X-ray emission do not follow the magnetic field lines. This result is of interest especially for loops as found in areas where the magnetic field emerging from active regions interacts with the surrounding network.     

Model of the slot region of earth''s electron radiation belt depending on the heliospheric parameters

Relativistic electrons in the slot region of Earth's electron radiation belt are studied using CORONAS-I satellite data obtained in March–May 1994. The strong shifts of the slot latitudinal location (from L≈3 to L≈2) were found. These shifts are associated with the Earth crossing sector boundaries formed by sector corotating structures of the solar wind and interplanetary magnetic field (IMF). The quantitative analysis of the relationship between properties of the slot region and values of the solar wind and IMF parameters was undertaken. The empirical model of the slot region dynamics was developed by means of Artificial Neural Network (ANN).     

Computational modeling of interference between helicopter main and tail rotors based on nonlinear blade vortex model

The Mi-2 helicopter was taken as an example for modeling the main and tail rotors operation with account for the interference at low speed sliding flight. Aerodynamic characteristics were calculated. The form of nonlinear vortical wake was obtained.     

Study of the mesosphere using wide-field twilight polarization measurements: Early results beyond the polar circle

This paper discusses the results of early measurements of temperature and dust in the mesosphere on the basis of wide-field twilight sky polarimetry, which began in 2015 in Apatity (North of Russia, 67.6° N, 33.4° E) using the original entire-sky camera. These measurements have been performed for the first time beyond the Polar Circle in the winter and early spring period. The general polarization properties of the twilight sky and the procedure for identifying single scattering are described. The key results of the study include the Boltzmann temperature values at altitudes higher than 70 km and the conclusion on a weak effect of dust on scattering properties of the mesosphere during this period.     

Study of the mode of angular velocity damping for a spacecraft at non-standard situation

Non-standard situation on a spacecraft (Earth’s satellite) is considered, when there are no measurements of the spacecraft’s angular velocity component relative to one of its body axes. Angular velocity measurements are used in controlling spacecraft’s attitude motion by means of flywheels. The arising problem is to study the operation of standard control algorithms in the absence of some necessary measurements. In this work this problem is solved for the algorithm ensuring the damping of spacecraft’s angular velocity. Such a damping is shown to be possible not for all initial conditions of motion. In the general case one of two possible final modes is realized, each described by stable steady-state solutions of the equations of motion. In one of them, the spacecraft’s angular velocity component relative to the axis, for which the measurements are absent, is nonzero. The estimates of the regions of attraction are obtained for these steady-state solutions by numerical calculations. A simple technique is suggested that allows one to eliminate the initial conditions of the angular velocity damping mode from the attraction region of an undesirable solution. Several realizations of this mode that have taken place are reconstructed. This reconstruction was carried out using approximations of telemetry values of the angular velocity components and the total angular momentum of flywheels, obtained at the non-standard situation, by solutions of the equations of spacecraft’s rotational motion.