Comparison of Antarctic riometer radio wave absorption and THEMIS mission energetic electron fluxes

Simultaneous observations of in situ plasma properties in the tail of the Earth’s magnetosphere and of ground based instruments, lying on the same geomagnetic field lines, have recently proved to yield significant new results. In most cases magnetosphere ionosphere interactions during the night-time northern hemisphere conditions are studied. Here, observations of energetic electrons in the tail of the Earth’s magnetosphere made by the THEMIS mission satellites are compared with auroral radio wave absorption determined by riometers in the Antarctic for sunlit conditions. Days for which satellites and riometers are connected by the same geomagnetic field line are selected using a geomagnetic field model. The six days analysed show clear associations between fluxes and absorptions in some cases. However, these do not necessarily correspond to conjugacy intervals. Hours of positive associations are 1.65 times those for negative associations, all hours and days considered (1.42–3.6 on five days and 0.58 on the other day). These computations are assumed appropriate since the footprints of the satellites used approximately follow corrected geomagnetic parallels for all six days studied. The use of a finer parameterization of geomagnetic models to determine conjugacy may be needed.     

A demonstration of the potential of Cryosat-2 to contribute to mesoscale observation

Cryosat-2 was designed for its primary scientific objectives, i.e. for cryosphere science. As far as oceanography is concerned, various mission design choices make it less accurate than missions designed to comply with ocean surface topography requirements such as Jason-2 or ENVISAT. Cryosat-2-specific errors are equivalent to more than 50% of the sea surface height variability over 40% of the oceans. Cryosat-2’s sampling pattern is also suboptimal for mesoscale observation because the satellite tracks from any consecutive period of 2 to 20 days (e.g. the most recent and most valuable data for near real time mesoscale observation) are aggregated in 500 km wide bands which are interleaved with 500 km wide observation gaps.     

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.     

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.     

Integral method for calculating a turbofan engine thrust reverser jet interacting with an external stream

An integral method is proposed for calculating a jet propagating from the turbofan engine thrust reverser and interacting with a stream formed as a result of the after-landing aircraft run. The calculation results for the PS-90 engine are presented. The calculation data obtained show that the mathematical model developed adequately describes a qualitative pattern of main parameter variation in the sector jet propagating in a stream.     

Investigation of gas dynamics numerical schemes based on the solution of one-dimensional test problems

The application of the optimized dispersion relation preserving scheme (DRP-scheme) in combination with the explicit optimized two-layer Runge-Kutta scheme is presented to solve a system of one-dimensional and quasi-one-dimensional Euler equations using as an example the solution of four test problems, namely, discontinuity disintegration in a tube (Sod’s problem); transfer of the lowamplitude Gaussian pulse; acoustic wave propagation through the transonic nozzle; acoustic wave-shock interaction. Also given are the comparison of the calculation results using different schemes: DRP, CABARET, CE-SE and the standard Lax-Wendroff schemes as well as the solutions obtained with the use of software packages.     

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.     

High-speed streams of the solar wind near the earths’s orbit and their sources on the sun according to stereoscopic observations in the minimum of the 23rd cycle

Based on the results of plasma and magnetic measurements at three different points of the heliosphere and telescopic observations of the Sun from these points we study simultaneously high-speed streams (HSS) of the solar wind (SW) near the Earths’s orbit and coronal holes (CH) that have generated them. The data from spacecraft STEREO-A, STEREO-B, ACE, and SOHO are used together with ground-based observations from March 2007 to May 2008. In this period there existed HSS whose sources represented CH of various polarity, geometry, and location relative to the heliographic and heliomagnetic equators. Dependence of SW parameters on mutual positions of spacecraft with respect to CH and heliospheric current sheet, and also on heliolatitude and geometry of the CH is revealed. A difference of more than 5° in locations of spacecraft with respect to the heliospheric current sheet in November 2007 allowed us to discover a heliolatitude velocity gradient of the SW streams between the STEREO-A and STEREO-B spacecraft. On the average this gradient at that time was equal to 20 km/s per degree. Substantial variations in SW streams associated with variations of the HSS SW sources during a few hours or days were also observed. This variability makes it difficult to use the data of spacecraft STEREO-B for sufficiently accurate prediction of SW properties in the near-Earth space by the method of simple advanced ti me shift due to heliolongitude difference between a spacecraft and the Earth even in solar activity minimum.