Vertical shear at the equatorial F-region ionosphere during post-sunset hours

The characteristics of the equatorial F-region zonal plasma drift during post-sunset period have been investigated using the multi-frequency HF Doppler radar. The pattern of the zonal plasma drift is such that it starts with a westward drift during the pre-sunset hours, followed by an eastward drift shortly after the E-region sunset. The zonal plasma drift is characterized by the presence of a positive vertical shear around the post-sunset period and maximum shear is observed at the time of the peak of the pre-reversal enhancement in the vertical drift. The presence of vertical shear in the zonal drift is associated with the post-sunset velocity vortex existing at the equatorial F-region.     

Removing orbital debris with lasers

Orbital debris in low Earth orbit (LEO) are now sufficiently dense that the use of LEO space is threatened by runaway collision cascading. A problem predicted more than thirty years ago, the threat from debris larger than about 1 cm demands serious attention. A promising proposed solution uses a high power pulsed laser system on the Earth to make plasma jets on the objects, slowing them slightly, and causing them to re-enter and burn up in the atmosphere. In this paper, we reassess this approach in light of recent advances in low-cost, light-weight modular design for large mirrors, calculations of laser-induced orbit changes and in design of repetitive, multi-kilojoules lasers, that build on inertial fusion research. These advances now suggest that laser orbital debris removal (LODR) is the most cost-effective way to mitigate the debris problem. No other solutions have been proposed that address the whole problem of large and small debris. A LODR system will have multiple uses beyond debris removal. International cooperation will be essential for building and operating such a system.     

Realization of modes of satellite attitude motion with a small level of microaccelerations using electromechanical executive devices

Quasi-static microaccelerations are estimated for a satellite specially designed to perform space experiments in the field of microgravity. Three modes of attitude motion of the spacecraft are considered: passive gravitational orientation, orbital orientation, and semi-passive gravitational orientation. In these modes the lengthwise axis of the satellite is directed along the local vertical, while solar arrays lie in the orbit plane. The second and third modes are maintained using electromechanical executive devices: flywheel engines or gyrodynes. Estimations of residual microaccelerations are performed with the help of mathematical modeling of satellite’s attitude motion under the action of gravitational and aerodynamic moments, as well as the moment produced by the gyro system. It is demonstrated that all modes ensure rather low level of quasi-static microaccelerations on the satellite and provide for a fairly narrow region of variation for the vector of residual microacceleration. The semi-passive gravitational orientation ensures also a limited proper angular momentum of the gyro system.     

Activation of the tail open part during the magnetospheric storm

In each polar cap (PC) we mark out “old PC” observed during quiet time before the event under consideration, and “new PC” that emerges during the substorm framing the old one and expanding the PC total area. Old and new PCs are the areas for the magnetosphere old and new tail lobes, respectively. The new lobe variable magnetic flux Ψ₁ is usually assumed to be active, i.e. it provides the electromagnetic energy flux (Poynting flux) ɛ′ transport from solar wind (SW) into the magnetosphere. The old lobe magnetic flux Ψ₂ is supposed to be passive, i.e. it remains constant during the disturbance and does not participate in the transporting process which would mean the old PC electric field absolute screening from the convection electric field created by the magnetopause reconnection. In fact, screening is observed, but far from absolute. We suggest a model of screening and determine its quantitative characteristics in the selected superstorm. The coefficient of a screening is the β = Ψ₂/Ψ₀₂, where Ψ₀₂ = const is open magnetic flux through the old PC measured prior to the substorm, and Ψ₂ is variable magnetic flux through the same area measured during the substorm. We consider three various regimes of disturbance. In each, the coefficient β decreased during the loading phase and increased at the unloading phase, but the rates and amplitudes of variations exhibited a strong dependence on the regime. We interpreted decrease in β as a result of involving the old PC magnetic flux Ψ₂, which was considered to be constant earlier, in the Poynting flux ɛ′ transport process from solar wind into the magnetosphere. Transport process weakening at the subsequent unloading phase creates increase in β. Estimates showed that coefficient β during each regime and the computed Poynting flux ɛ′ varied manifolds. In general, unlike the existing substorm conception, the new scenario describes an unknown earlier of tail lobe activation process during a substorm growth phase that effectively increases the accumulated tail energy for the expansion and recovery phases.     

Long-term variations of fluxes of solar protons and helium isotopes

The fluxes of hydrogen and helium isotopes in the solar wind are reconstructed over a long time scale since the present time up to 600 million years back. Abundances of helium isotopes, obtained in the helium isotopic analysis made for 8 lunar soil samples, were used as initial data in the reconstruction procedure. Samples were taken off from various levels of the 1.6-m core of lunar soil delivered by the automatic Luna-24 station in 1976. The data on modern hydrogen and helium fluxes were used as well. The developed reconstruction procedure allowed one to select various solar wind components in a ??gross?? composition. Proton flux variations over the interval of 600 million years do not exceed a value of 40 %. Helium flux variations reach a value of 1.5?C2 relative to the average value. Most likely, this circumstance is caused by considerable variations of a number of coronal mass ejections (CME) enriched by helium. The arguments in favor of solar activity polycyclicity on a long time scale are discussed.     

Relationship between the fluxes of relativistic electrons at geosynchronous orbit and the level of ULF activity on the Earth’s surface and in the solar wind during the 23rd solar activity cycle

We present the results of a cross-correlation analysis made on the basis of Spearman’s rank correlation method. The quantities to correlate are daily values of the fluence of energetic electrons at a geosynchronous orbit, intensities of ground and interplanetary ultra-low-frequency (ULF) oscillations in the Pc5 range, and parameters of the solar wind. The period under analysis is the 23rd cycle of solar activity, 1996–2006. Daily (from 6 h to 18 h of LT) magnetic data at two diametrically opposite observatories of the Intermagnet network are taken as ground-based measurements. The fluxes of electrons with energies higher than 2 MeV were measured by the geosynchronous GOES satellites. The data of magnetometers and plasma instruments installed on ACE and WIND spacecraft were used for analysis of the solar wind parameters and of the oscillations of the interplanetary magnetic field (IMF). Some results elucidating the role played by interplanetary ULF waves in the processes of generation of magneospheric oscillations and acceleration of energetic electrons are obtained. Among them are (i) high and stable correlation of ground ULF oscillations with waves in the solar wind; (ii) closer link of mean daily amplitudes of both interplanetary and ground oscillations with ‘tomorrow’ values of the solar wind velocity than with current values; and (iii) correlation of the intensity of ULF waves in the solar wind, normalized to the IMF magnitude, with fluxes of relativistic electrons in the magnetosphere.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

Computer modelling of a penetrator thermal sensor

The Philae lander is part of the Rosetta mission to investigate comet 67P/Churyumov-Gerasimenko. It will use a harpoon like device to anchor itself onto the surface. The anchor will perhaps reach depths of 1–2 m. In the anchor is a temperature sensor that will measure the boundary temperature as part of the MUPUS experiment. As the anchor attains thermal equilibrium with the comet ice it may be possible to extract the thermal properties of the surrounding ice, such as the thermal diffusivity, by using the temperature sensor data. The anchor is not an optimal shape for a thermal probe and application of analytical solutions to the heat equation is inappropriate. We prepare a numerical model to fit temperature sensor data and extract the thermal diffusivity. Penetrator probes mechanically compact the material immediately surrounding them as they enter the target. If the thermal properties, composition and dimensions of the penetrator are known, then the thermal properties of this pristine material may be recovered although this will be a challenging measurement. We report on investigations, using a numerical thermal model, to simulate a variety of scenarios that the anchor may encounter and how they will affect the measurement.     

Solar rhythms in the characteristics of the Arctic frontal zone in the North Atlantic

Long-term changes of the Arctic frontal zone characteristics near the south-eastern coasts of Greenland were considered, the NCEP/NCAR reanalysis data being used. It was found that in the cold half of the year the temperature gradients in the layer 1000–500 hPa in the region under study reveal strong ∼10-yr and ∼22-yr periodicities that seem to be related to solar activity cycles. The results obtained suggest the influence of solar activity and cosmic ray variations on the structure of the temperature field of the troposphere resulting in the changes of the temperature contrasts in the Arctic frontal zone that, in turn, may affect the intensity of cyclogenesis at middle latitudes. The detected effects seem to indicate an important part of frontal zones in the mechanism of solar activity and cosmic ray variation influence on the development of extratropical baric systems. It is suggested that the variations of the temperature gradients revealed in the Arctic frontal zone are due to the radiative forcing of cloudiness changes which may be associated with geomagnetic activity and cosmic ray variations.     

Studies of the geocenter motion using 16-years DORIS data

An accuracy of geocenter motion estimation is strongly dependent on the geodetic network size and stations distribution over the Earth’s surface. From this point of view DORIS system has an advantage, as its ground network of beacons consists of more than 50 sites, equally distributed over the Earth’s surface. Aiming to study variations of the geocenter movements, the results of DORIS data analysis for the time span 1993.0–2009.0 (inawd06.snx series), performed at the Analysis Centre of the Institute of astronomy of the Russian Academy of Sciences, have been used. DORIS data processing was made with GIPSY/OASIS II software, developed by Jet Propulsion Laboratory and modified for DORIS data processing by Institute Géographique National. Standard deviations of stations coordinates are estimated at the level 0.5–4.0 cm (internal consistency), depending on the number of satellites used in the solution. RMS of estimated components of the DORIS satellites orbits, compared with the solutions of other IDS analysis centres, do not exceed 1–2 cm. Weekly solutions for coordinates have been transformed from free network solutions (inawd06.snx series) to a well defined terrestrial reference frame ITRF2005 with the use of seven parameters of Helmert transformation, which were examined with a view to study variations of the geocenter movements (ina05wd01.geoc time series). In order to estimate linear trend, amplitudes, periods and phases of geocenter variation a method of linear regression was applied. The evaluated amplitudes of annual variations are of the order of 5–7 mm for X and Y components and 27–29 mm for Z component. Semi-annual amplitudes are also noticeable in all components (1–34 mm for X, Y and Z components). Secular trends in the DORIS geocenter coordinates are: −1.2, −0.1 and −0.3 mm/year for X, Y and Z directions respectively.