Neutron monitor asymptotic directions of viewing during the event of 13 December 2006

During the recent ground level enhancement of 13 December 2006, also known as GLE70, solar cosmic ray particles of energy bigger that ∼500 MeV/nucleon propagated inside the Earth’s magnetosphere and finally accessed low-altitude satellites and ground level neutron monitors. The magnitude and the characteristics of this event registered at different neutron monitor stations of the worldwide network can be interpreted adequately on the basis of an estimation of the solar particle trajectories in the near Earth interplanetary space. In this work, an extended representation of the Earth’s magnetic field was realized applying the Tsyganenko 1989 model. Using a numerical back-tracing technique the solar proton trajectories inside the magnetospheric field of the Earth were calculated for a variety of particles, initializing their travel at different locations, covering a wide range of energies. In this way, the asymptotic directions of viewing were calculated for a significant number of neutron monitor stations, providing crucial information on the Earth’s “magnetospheric optics” for primary solar cosmic rays, on the top of the atmosphere, during the big solar event of December 2006. The neutron monitor network has been treated, therefore, as a multidimensional tool that gives insights into the arrival directions of solar cosmic ray particles as well as their spatial and energy distributions during extreme solar events.     

Numerical modeling of propagation of breaking nonlinear acoustic-gravity waves from the lower to the upper atmosphere

Acoustic-gravity waves (AGWs) observed in the upper atmosphere may be generated near the Earth’s surface due to a variety of meteorological sources. Two-dimensional simulations of vertical propagation and breaking of nonlinear AGWs in the atmosphere are performed. Forcing near the Earth’s surface is used as the AGW source in the model. We use a numerical method based on finite-difference analogues of fundamental conservation laws for solving atmospheric hydrodynamic equations. This approach selects physically correct generalized solutions of the wave hydrodynamic equations. Numerical simulations are performed in a representative region of the Earth’s atmosphere up to altitude 500 km. Vertical profiles of temperature, density, molecular viscosity and heat conductivity were taken from the standard atmosphere model MSIS-90 for January. Calculations were made for different amplitudes and frequencies of lower boundary wave forcing. It is shown that after activating the tropospheric wave forcing, the initial pulse of AGWs may very quickly propagate to altitudes of 100 km and above and relatively slowly dissipate due to molecular viscosity and heat conduction. This may increase the role of transient nonstationary waves in effective energy transport and variations of atmospheric parameters and gas admixtures in a broad altitude range.     

December 2006 solar extreme events and their influence on the near-Earth space environment: “Universitetskiy-Tatiana” satellite observations

The Russian microsatellite “Universitetskiy-Tatiana” was launched on Jan. 20, 2005 and was both a scientific and educational mission. Its two main aims were declared as: (1) monitoring of the energetic particles dynamics in the near-Earth space environment after solar events and during quiet times, (2) educational activities based on experimental data obtained from the spacecraft. In this paper observations acquired during Dec. 5–16, 2006, known as “Solar Extreme Events 2006”, were analyzed. The “Universitetskiy-Tatiana” microsatellite orbit permits one to measure both solar energetic particle dynamics, variations of the boundary of solar particle penetration, as well as relativistic and sub-relativistic electrons of the Earth’s outer radiation belt during and after magnetic storms. Both relativistic electrons of the Earth’s outer radiation and solar energetic particles are an important source of radiation damage in near-Earth space. Therefore, the presented experimental results demonstrate the successful application of a small educational spacecraft both for scientific and educational programs.     

Solar neutron events as a tool to study particle acceleration at the Sun

The Sun provides unique opportunities to study particle acceleration mechanisms using data from detectors placed on the Earth’s surface and on board spacecrafts. Particles may gain high energies by several physical mechanisms. Differentiating between these possibilities is a fundamental problem of cosmic ray physics. Energetic neutrons provide us with information that keeps the signatures of the acceleration site. A summary of some representative solar neutron events observed on the Earth’s surface, including associated X and γ-ray observations from spacecrafts is presented. We discuss evidence of acceleration of particles by the Sun to energies up to several tens of GeV. In addition, a recent solar neutron event that occurred on September 7th 2005 and detected by several observatories at Earth is analyzed in detail.     

Fluctuations of cosmic rays and IMF in the vicinity of interplanetary shocks

Fluctuations of cosmic rays and interplanetary magnetic field upstream of interplanetary shocks are studied using data of ground-based polar neutron monitors as well as measurements of energetic particles and solar wind plasma parameters aboard the ACE spacecraft. It is shown that coherent cosmic ray fluctuations in the energy range from 10 keV to 1 GeV are often observed at the Earth’s orbit before the arrival of interplanetary shocks. This corresponds to an increase of solar wind turbulence level by more than the order of magnitude upstream of the shock. We suggest a scenario where the cosmic ray fluctuation spectrum is modulated by fast magnetosonic waves generated by flux of low-energy cosmic rays which are reflected and/or accelerated by an interplanetary shock.     

Universal scaling laws for fully-developed magnetic field turbulence near and far upstream of the Earth’s bow shock

We analyze the multifractal scaling of the modulus of the interplanetary magnetic field near and far upstream of the Earth’s bow shock, measured by Cluster and ACE, respectively, from 1 to 3 February 2002. The maximum order of the structure function is carefully estimated for each time series using two different techniques, to ensure the validity of our high-order statistics. The first technique consists of plotting the integrand of the pth order structure function, and the second technique is a quantitative method which relies on the power-law scaling of the extreme events. We compare the scaling exponents computed from the structure functions of magnetic field differences with the predictions obtained by the She–Lévêque model of intermittency in anisotropic magnetohydrodynamic turbulence. Our results show a good agreement between the model and the observations near and far upstream of the Earth’s bow shock, rendering support for the modelling of universal scaling laws based on the Kolmogorov phenomenology in the presence of sheet-like dissipative structures.     

Measurements of the radiation quality factor Q at aviation altitudes during solar minimum (2006–2008)

In radiation protection, the Q-factor has been defined to describe the biological effectiveness of the energy deposition or absorbed dose to humans in the mixed radiation fields at aviation altitudes. This particular radiation field is generated by the interactions of primary cosmic particles with the atoms of the constituents of the Earth’s atmosphere. Thus the intensity, characterized by the ambient dose equivalent rate H∗(10), depends on the flight altitude and the energy spectra of the particles, mainly protons and alpha particles, impinging on the atmosphere. These charged cosmic projectiles are deflected both by the interplanetary and the Earth’s magnetic field such that the corresponding energy spectra are modulated by these fields. The solar minimum is a time period of particular interest since the interplanetary magnetic field is weakest within the 11-year solar cycle and the dose rates at aviation altitudes reach their maximum due to the reduced shielding of galactic cosmic radiation. For this reason, the German Aerospace Center (DLR) performed repeated dosimetric on-board measurements in cooperation with several German airlines during the past solar minimum from March 2006 to August 2008. The Q-factors measured with a TEPC range from 1.98 at the equator to 2.60 in the polar region.     

Influence of Earth’s shadowing effects on space debris stability

Solar radiation pressure affects the evolution of high area-to-mass geostationary space debris. In this work, we extend the stability study of Valk et al. (2009) by considering the influence of Earth’s shadows on the short- and long-term time evolutions of space debris. To assess the orbits stability, we use the Mean Exponential Growth factor of Nearby Orbits (MEGNO), which is an efficient numerical tool to distinguish between regular and chaotic behaviors. To reliably compute long-term space debris motion, we resort to the Global Symplectic Integrator (GSI) of Libert et al. (2011) which consists in the symplectic integration of both Hamiltonian equations of motion and variational equations. We show how to efficiently compute the MEGNO indicator in a complete symplectic framework, and we also discuss the choice of a symplectic integrator, since propagators adapted to the structure of the Hamiltonian equations of motion are not necessarily suited for the associated variational equations. The performances of our method are illustrated and validated through the study of the Arnold diffusion problem. We then analyze the effects of Earth’s shadows, using the adapted conical and cylindrical Earth’s shadowing models introduced by Hubaux et al. (2012) as the smooth shadow function deriving from these models can be easily included into the variational equations. Our stability study shows that Earth’s shadows greatly affect the global behaviour of space debris orbits by increasing the size of chaotic regions around the geostationary altitude. We also emphasize the differences in the results given by conical or cylindrical Earth’s shadowing models. Finally, such results are compared with a non-symplectic integration scheme.     

Deterministic pion and muon transport in Earth’s atmosphere

An accurate understanding of the physical interactions and transport of space radiation is important for safe and efficient space operations. Secondary particles produced by primary particle interactions with intervening materials are an important contribution to radiation risk. Pions are copiously produced in the nuclear interactions typical of space radiations and can therefore be an important contribution to radiation exposure. Charged pions decay almost exclusively to muons. As a consequence, muons must also be considered in space radiation exposure studies. In this work, the NASA space radiation transport code HZETRN has been extended to include the transport of charged pions and muons. The relevant transport equation, solution method, and implemented cross sections are reviewed. Muon production in the Earth’s upper atmosphere is then investigated, and comparisons with recent balloon flight measurements of differential muon flux are presented. Muon production from the updated version of HZETRN is found to match the experimental data well.     

Assessment of the consequences of the Fengyun-1C breakup in low Earth orbit

On 11 January 2007, the People’s Republic of China conducted a successful anti-satellite test against one of their defunct polar-orbiting weather satellites. The target satellite, called Fengyun-1C, had a mass of 880 kg and was orbiting at an altitude of about 863 km when the collision occurred. Struck by a direct-ascent interceptor at a speed of 9.36 km/s, the satellite disintegrated, spreading the cataloged fragments between 200 and 4000 km, with the highest concentration near the breakup height. By the end of April 2008, 2377 pieces of debris, including the original payload remnant, had officially been cataloged by the US Space Surveillance Network. Of these, nearly 1% had reentered the Earth’s atmosphere. This deliberate act is the largest debris-generating event on record, and its consequences will adversely affect circumterrestrial space for many years.     

The IGY and beyond: A brief history of ground-based cosmic-ray detectors

The state of art of ground-based cosmic-ray research from its discovery to present is reviewed. After discovery of cosmic rays by Hess in 1912, the nature of the primary and secondary radiation was established from recordings by a variety of instruments, sensitive to various components of cosmic rays and operated at different latitudes, longitudes and altitudes, including instruments carried by balloons. The IGY formalized international co-operation and coordinated study of cosmic rays, which is vital for meaningful interpretation of cosmic-ray data. Data collected at different geographic locations require an effective cutoff rigidity as a data ordering parameter. This parameter is obtained from tracing trajectories of primary cosmic rays in the Earth’s magnetic field. After 50 years the world’s neutron monitor network remains still the backbone for studying intensity variations of primary cosmic rays in the rigidity ranges between 1 and 15 GV, associated with transport and with transient events. Also the penetrating muon and neutrino components of secondary cosmic rays have a long history of recording and fundamental problem investigations. Valuable data about composition and spectrum of primary cosmic rays in ever increasing high-energy regions have been obtained during the years of investigations with various configurations and types of extensive air shower detectors. The culture of personal involvement of the physicist in carrying out experiments and data acquisition characterized the continued vitality of cosmic-ray investigations ranging from its atmospheric, geomagnetic and heliospheric transport through to its solar and astrophysical origins.     

A review of global satellite-derived snow products

Snow cover over the Northern Hemisphere plays a crucial role in the Earth’s hydrology and surface energy balance, and modulates feedbacks that control variations of global climate. While many of these variations are associated with exchanges of energy and mass between the land surface and the atmosphere, other expected changes are likely to propagate downstream and affect oceanic processes in coastal zones. For example, a large component of the freshwater flux into the Arctic Ocean comes from snow melt. The timing and magnitude of this flux affects biological and thermodynamic processes in the Arctic Ocean, and potentially across the globe through their impact on North Atlantic Deep Water formation.     

Solar proton enhancements in different energy channels and coronal mass ejections during the last solar cycle

The main properties of 11622 coronal mass ejections (CMEs) observed by the Solar and Heliospheric Observatory (SOHO) mission’s Large Angle and Spectrometric Coronagraph (LASCO-C2) from January 1996 through December 2006 are considered. Moreover, the extended database of solar proton enhancements (SPEs) with proton flux >0.1 pfu at energy >10 MeV measured at the Earth’s orbit is also studied. A comparison of these databases gives new results concerning the sources and acceleration mechanisms of solar energetic particles. Specifically, coronal mass ejections with width >180° (wide) and linear speed >800 km/s (fast) seem they have the best correlation with solar proton enhancements. The study of some specific solar parameters, such as soft X-ray flares, sunspot numbers, solar flare index etc. has showed that the soft X-ray flares with importance >M5 may provide a reasonable proxy index for the SPE production rate. From this work, it is outlined that the good relation of the fast and wide coronal mass ejections to proton enhancements seems to lead to a similar conclusion. In spite of the fact that in the case of CMEs the statistics cover only the last solar cycle, while the measurements of SXR flares are extended over three solar cycles, it is obvious for the studied period that the coronal mass ejections can also provide a good index for the solar proton production.     

A lunar cargo mission design strategy using variable low thrust

A design technique for a near optimal, Earth–Moon transfer trajectory using continuous variable low thrust is proposed. For the Earth–Moon transfer trajectory, analytical and numerical methods are combined to formulate the trajectory optimization problem. The basic concept of the proposed technique is to utilize analytically optimized solutions when the spacecraft is flying near a central body where the transfer trajectories are nearly circular shaped, and to use a numerical optimization method to match the spacecraft’s states to establish a final near optimal trajectory. The plasma thruster is considered as the main propulsion system which is currently being developed for crewed/cargo missions for interplanetary flight. The gravitational effects of the 3rd body and geopotential effects are included during the trajectory optimization process. With the proposed design technique, Earth–Moon transfer trajectory is successfully designed with the plasma thruster having a thrust direction sequence of “fixed-varied-fixed” and a thrust acceleration sequence of “constant-variable-constant”. As this strategy has the characteristics of a lesser computational load, little sensitivity to initial conditions, and obtaining solutions quickly, this method can be utilized in the initial scoping studies for mission design and analysis. Additionally, derived near optimal trajectory solution can be used as for initial trajectory solution for further detailed optimization problem. The demonstrated results will give various insights into future lunar cargo trajectories using plasma thrusters with continuous variable low thrust, establishing approximate costs as well as trajectory characteristics.     

Plasma flow variations and energetic protons upstream of the earth’s bow shock: A statistical study

Foreshock is a special region located upstream of the Earth’s bow shock characterized by the presence of various plasma waves and fluctuations caused by the interaction of the solar wind plasma with particles reflected from the bow shock or escaping from the magnetosphere. On the other hand, foreshock fluctuations may modify the bow shock structure and, being carried through the magnetosheath, influence the magnetopause. During the years 1995–2000, the INTERBALL-1 satellite made over 10,000 hours of plasma and energetic particles measurements in the solar wind upstream of the Earth’s bow shock. We have sorted intervals according to the level of solar wind ion flux fluctuations and/or according to the flux of back-streaming energetic protons. An analysis of connection between a level of ion flux fluctuations and fluxes of high-energy protons and their relation to the IMF orientation is presented.     

The ground GNSS tomography – unconstrained approach

The GNSS signal along its more than 20,000 km line of sight is bended, attenuated and delayed. These effects are results of the Earth’s atmosphere, and Sun radiation. Amongst all mentioned effects the one considered in this paper is a signal phase delay in troposphere.     

Reducing variability in short term orbital lifetime prediction

Within the last year three major re-entries occurred. The satellites UARS, ROSAT and Phobos-Grunt entered Earth’s atmosphere with fragments reaching the surface. Due to a number of uncertainties in propagating an object’s trajectory the exact place and time of a satellite’s re-entry is hard to determine. Major influences when predicting the re-entry time are the changing precision of the available orbital data, the satellite’s ballistic coefficient, the activity of the sun which influences the Earth’s atmosphere and the underlying quality of the atmospheric model. In this paper a method is presented which can reduce the variability in short-term orbital lifetime prediction induced by fluctuating orbital data accuracies. A re-entry campaign is used as a reference for this purpose. For a window of a few weeks before the re-entry the position data of a synthetic object is disturbed considering different degrees of orbital data errors. As a result different predictions will exist for the generated position data of a given day. Using a regression algorithm on the available data an average position is obtained, which is then used for the orbital lifetime prediction. The effect of this measure is a more consistent prediction of the orbital lifetime. The paper concludes with the comparison of the generated re-entry windows in various test cases for the original and the averaged data.     

Model for the VLF/LF radio signal anomalies formation associated with earthquakes

The influence of quasi-static electric field of seismic origin on the characteristics of the internal gravity waves (IGWs) in the Earth’s ionosphere is considered. The electric field in the ionosphere arises due to the injection of charged aerosols into the atmosphere, formation of an EMF in the near Earth atmosphere and perturbation of the conductive electric current in the global electric circuit. Amplification of the electric current in seismic zone is accompanied by the formation of perturbation of the lower ionosphere that affects the amplitude and phase of VLF/LF signals. The action of the electric field on the IGWs is connected with the appearance of the Ampere’s force in the ionosphere. In the spectral range of these waves the latter acts on the neutral component of the ionosphere plasma. As the result of this interaction the ionosphere starts to support the discrete spectrum of oscillations. Periods of their maximums increase as numbers of natural sequence. The existence of such peculiarities of the waves in the ionosphere is confirmed by observations.     

Model of energetic populations at Ganymede,implications for an orbiter

A model is developed to study the energetic particle populations in Ganymede’s magnetosphere. The main objective is to estimate to what extent the moon could protect an orbiter from radiations. Using Liouville’s theorem, the phase space density of particles coming from Jupiter’s magnetosphere is calculated at any point of Ganymede’s environment. Up to energies of ∼50–100 keV for ions and ∼10–20 MeV for electrons, Ganymede’s magnetic field appears to be able to form distinctive populations as loss-cones over the polar caps and radiation belts. At larger energies, these features are blurred by Larmor radius effects; the moon absorption simply creates a quasi-isotropic layer of ∼500 km thickness where the flux is reduced by ∼40–50%. The predictions are compared to Galileo measurements. In particular, we demonstrate the importance of the moon sweeping in reducing the flux over the polar caps. Interestingly, this can be accounted for by assuming that the particles bouncing between Jupiter and Ganymede are ideally scattered in pitch angle and permanently re-fill the loss-cone, which increases the precipitation on Ganymede’s polar cap. In overall, it is estimated that the radiation dose received by an orbiter of Ganymede will be reduced by more than 50–60% compared to the expected dose at Jupiter/Ganymede distance. This should have a positive impact on the design of a future orbiter of Ganymede.     

磁场重联的二维粒子模拟研究

利用二维全粒子模拟方法研究了无碰撞等离子体中的磁场重联过程，得到了不同区域的离子和电子速度分布．计算结果表明，电子和离子在扩散区中的不同动力学特性产生的Hall电流使磁场的y分量By呈现四极形分布．离子和电子的速度分布偏离了初态时的Maxwell分布，呈现非局域的多重分布．同时由于磁场重联而产生的电场使电子在X点附近得到加速和加热，因而在电子的能谱分布中形成-高能尾。