Spectral line broadening by relativistic electrons in plasmas: Collision operator

In the present work we compute the real part of the impact collision operator for the electron broadening of ion lines in plasmas, taking into account relativistic effects in the dynamics of the perturbing electron. Specifically two relativistic effects are included: The modification of the trajectory due to non-Newtonian mechanics and the modification of the velocity distribution (Maxwell–Juttner). The results are compared to the non-relativistic case.     

A comparison of total reaction cross section models used in particle and heavy ion transport codes

Understanding the interactions and propagations of high energy protons and heavy ions are essential when trying to estimate the biological effects of Galactic Cosmic Rays (GCR) and Solar Particle Events (SPE) on personnel in space. To be able to calculate the shielding properties of different materials and radiation risks, particle and heavy ion transport codes are needed. In all particle and heavy ion transport codes, the probability function that a projectile particle will collide within a certain distance x in the matter depends on the total reaction cross sections, and the calculated partial fragmentation cross sections scale with the total reaction cross sections. It is therefore crucial that accurate total reaction cross section models are used in the transport calculations. In this paper, different models for calculating nucleon–nucleus and nucleus–nucleus total reaction cross sections are compared with each other and with measurements. The uncertainties in the calculations with the different models are discussed, as well as their overall performances with respect to the available experimental data. Finally, a new compilation of experimental data is briefly presented.     

On the Lagrange and Hill stability of motion in the three-body problem

In the three-body problem, we consider the Lagrange and Hill stability including the Lagrange stability for the manifold of symmetric motions that exists in the case where two of three bodies have equal masses. To analyze the stability, in addition to integrals of energy and angular momentum we use the Lagrange–Jacobi equality. We prove theorems on the Lagrange and Hill stability. The theorem on the Hill stability has effective application in the case where the mass of a body is much less than masses of two other bodies. In this case, as it is known, the model of the restricted three-body problem is usually applied.     

The extension of the INCL model for simulation of shielding in space

Radiation hazard for space missions is mainly due to cosmic ray protons, helium nuclei and light ions, whose energy spectrum is maximum around 1 GeV per nucleon but remains non-negligible for energies up to 15 GeV per nucleon. Nuclear reactions induced by high energy protons are often described by intranuclear cascade plus evaporation models. The attention is focused here on the Liège Intranuclear Cascade model (INCL), which has been shown to reproduce fairly well a great deal of experimental data for nucleon-induced reactions in the 200 MeV to 2 GeV range, when coupled with the ABLA evaporation-fission code. In order to extend the model to other conditions relevant for space radiation, three improvements of INCL are under development. They are reported on here. First, the reaction model has been extended to nucleon–nucleus reactions at incident energies up to 15 GeV, mainly by the inclusion of additional pion production channels in nucleon–nucleon collisions during the cascade. Second, a coalescence mechanism for the emission of light charged particles has been implemented recently. Finally, the model has been modified in order to accommodate light ions as projectiles. First results are shown and compared with illustrative experimental data. Implications for issues concerning radiation protection in space are discussed.     

Design and analysis of an extended mission of CE-2: From lunar orbit to Sun–Earth L2 region

Chang’E-2 (CE-2) has firstly successfully achieved the exploring mission from lunar orbit to Sun–Earth L2 region. In this paper, we discuss the design problem of transfer trajectory and at the same time analyze the visible segment of Tracking, Telemetry & Control (TT&C) system for this mission. Firstly, the four-body problem of Sun–Earth–Moon and Spacecraft can be decoupled in two different three-body problems (Sun–Earth + Moon Restricted Three-Body Problems (RTBPs) and Earth–Moon ephemeris model). Then, the transfer trajectory segments in different model are computed, respectively, and patched by Poincaré sections. The full-flight trajectory including transfer trajectory from lunar orbit to Sun–Earth L2 region and target Lissajous orbit is obtained by the differential correction method. Finally, the visibility of TT&C system at the key time is analyzed. Actual execution of CE-2 extended mission shows that the trajectory design of CE-2 mission is feasible.     

Observation of Mercury’s sodium tail using Fabry–Perot Interferometer

We observed sodium emission from Mercury’s atmosphere using a Fabry–Perot Interferometer at Haleakala Observatory on June 14, 2006. The Fabry–Perot Interferometer was used as a wavelength-tunable filter. The spectra of the surface reflection were subtracted from the observed spectra because sodium emission is contaminated by the surface reflection of Mercury. The image obtained in our observation shows the sodium exosphere extended to the anti-solar direction. The lifetime of a sodium atom was estimated to be 1.6 × 10⁴ to 1.9 × 10⁵ s with an error by a factor of 3–4.     

An empirical model for the variation of rain drop size distribution with rain rate at a few locations in southern India

A Joss–Waldvogel impact type disdrometer was installed at four different locations in the Indian peninsula during various periods from 2001 till date. The data are analysed to study the nature of rain drop size distribution (DSD) in this region. Out of the three well known distributions that describe DSD, namely, the Marshall–Palmer, Gamma and Lognormal, it has been found that Lognormal distribution fits the DSD in this region better than the other ones. Lognormal distributions for different rain rates were then derived by fitting the lognormal function to the data using a curve fitting software. Then the variation of fit parameters with rain rate was evaluated. Incorporating these variations, into the Lognormal distribution, an empirical equation that describes the DSD in this region for different rain rates was derived. Then this equation was tested with sample data from each of these stations. The data used for validation were not used for fitting lognormal equation to derive the fit parameters. The correlation between the DSD measured and derived using the empirical model was found to be quite good (0.9) except in some cases where the coefficient dropped to 0.75. The empirical model can be updated when more data are available.     

Plasma electron temperature variability in lunar surface potential and in electric field under average solar wind conditions

The Moon is immersed in plasma environment. The most interesting challenge of the lunar plasma– field environment is that it is alternatively dominated by the extended but variable outer atmosphere of the Earth – the magnetosphere – and by the extended but highly variable solar atmosphere – the solar wind. Understanding the plasma environment and its interaction with the lunar surface will be beneficial to both manned and robotic surface exploration activities and to scientific investigations. Presented is a preliminary map of variations of lunar surface electric potential over the day side and night side using probe equations and a discussion on dust dynamics in this E-field structure using the data from Electron Reflectometer in Lunar Prospector spacecraft during 1998–1999. On the day side, potential is around 5 V and on the night side it reaches up to −82 V. On the night side region, only highly energetic electrons can overcome this large negative potential. The variation at electron temperature (T_e) strongly reflects in the surface potential. The potential reaches to a value of −82 V for T_e = 58 eV. Surface charging causes the electrostatic transport of charged dust grains. Dust grain size of 0.1 μm shows a levitation height of 4.92 m on lunar day side, 748 m on terminator region and 3.7 km on the night side. The radius of maximum sized grain to be lofted, R_max, peaks at the terminator region (R_max = 0.83 μm). At the transition region dust levitation is almost absent. This region is most suited for exploration activities as the region is free from hazards caused by lunar dust.     

Dependence of the 27-day variation of cosmic rays on the global magnetic field of the Sun

We show that the higher range of the heliolongitudinal asymmetry of the solar wind speed in the positive polarity period (A > 0) than in the negative polarity period (A < 0) is one of the important reasons of the larger amplitudes of the 27-day variation of the galactic cosmic ray (GCR) intensity in the period of 1995–1997 (A > 0) than in 1985–1987 (A < 0). Subsequently, different ranges of the heliolongitudinal asymmetry of the solar wind speed jointly with equally important corresponding drift effect are general causes of the polarity dependence of the amplitudes of the 27-day variation of the GCR intensity. At the same time, we show that the polarity dependence is feeble for the last unusual minimum epoch of solar activity 2007–2009 (A < 0); the amplitude of the 27-day variation of the GCR intensity shows only a tendency of the polarity dependence. We present a three dimensional (3-D) model of the 27-day variation of GCR based on the Parker’s transport equation. In the 3-D model is implemented a longitudinal variation of the solar wind speed reproducing in situ measurements and corresponding divergence-free interplanetary magnetic field (IMF) derived from the Maxwell’s equations. We show that results of the proposed 3-D modeling of the 27-day variation of GCR intensity for different polarities of the solar magnetic cycle are in good agreement with the neutron monitors experimental data. To reach a compatibility of the theoretical modeling with observations for the last minimum epoch of solar activity 2007–2009 (A < 0) a parallel diffusion coefficient was increased by ∼40%.     

Effects of sector structure of the interplanetary magnetic field on the upper mesosphere–lower thermosphere dynamics

In this paper we study the influence of the interplanetary magnetic field (IMF) polarity changes caused by the Earth passing through the IMF sector boundary on the dynamic processes taking place in neutral atmosphere within the altitude interval of the upper mesosphere–lower thermosphere (83–101 km). The analysis has revealed the influence of the IMF sector structure on dynamics of the upper mesosphere–lower thermosphere. There has been a significant seasonal variation of the wind reaction to the IMF polarity changes observed. The influence of the IMF polarity changes on neutral atmosphere dynamics within the altitude range of 83–101 km is most pronounced in the zonal component of neural wind when the IMF polarity changes from negative to positive in all the seasons except for spring and when IMF polarity changes from positive to negative – in spring only.     

基于源项辨识的飞机座舱污染浓度动态预测

随着大型民机飞行时间的延长,座舱空气污染事故发生概率也随之增大,快速准确的污染浓度预测对保证乘客生命安全具有重要意义.座舱各污染浓度的动态预测和污染源项强度辨识是实现座舱空气质量实时预测的关键技术.污染源项散发强度辨识,如采用最小二乘算法,参数估计是静态的,一般延迟较大;如采用单模卡尔曼滤波算法,虽能实现动态辨识,但不能同时兼顾稳态和过渡过程(突发污染)的参数估计性能,导致误差较大.为解决上述难题,本文提出基于源项辨识的飞机座舱污染浓度动态预测方法,同时完成污染源散发强度动态辨识和污染浓度状态实时预测.该算法由2个滤波器组成,分别用于匹配系统的稳态和突发过渡过程特征,提高浓度方程参数估计和状态预测性能,保证飞机座舱空气质量态势预测的快速性和准确性.仿真结果证实了该算法的有效性.     

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.     

Decadal variability of tropical Pacific temperature in relation to solar cycles

We use the 8-year long satellite temperature data (2002–2010) from Atmospheric InfraRed Sounder (AIRS) and Atmospheric Microwave Sounding Unit (AMSU) on the Aqua satellite to identify temperature trends in the troposphere and low stratosphere over the Niño 3.4 region of the Tropical Pacific Ocean in the most recent 11-year solar cycle. Employing more extended sea surface temperature (SST) data for five solar cycles (1950–2009) in this region we show that the satellite trends reflect a typical decrease of the sea surface temperature (SST) in the Niño 3.4 region in the declining phase of the solar cycle. The magnitude of the SST decrease depends on the solar cycle and ranges between 0.07 K/yr and 0.27 K/yr for the last five solar cycles.     

Extension of INCL4 between 2 and 15 GeV

The intranuclear cascade model INCL4 has been shown to be very successful for describing, without adjustable parameters, a whole set of data for p-induced reactions in the 40 MeV–2 GeV energy range. In view of its possible application to cosmic ray interactions, the INCL4 code has been extended to the 2–15 GeV energy range, so covering a large part of the spectrum of the incident energy of the cosmic rays.     

Relative amplitude of the total electron content variations depending on geomagnetic activity

We developed a method of estimation of a relative amplitude dI/I of the total electron content (TEC) variations in the ionosphere as deduced from the data of the global GPS receivers network. To obtain statistically significant results we picked out three latitudinal belts provided in the Internet by the maximum number of GPS sites. They are high-latitudinal belt (50–80°N, 200–300°E; 59 sites), mid latitude belt (20–50°N, 200–300°E; 817 sites), and equatorial belt (±20°N, 0–360°E; 76 sites). The results of the analysis of the diurnal and latitudinal dependencies of dI/I and dI/I distribution probability for 52 days with different levels of geomagnetic activity are presented. It was found that on average the relative amplitude of the TEC variations varies within the range 0–10% proportionally to the value of the K_p geomagnetic index. In quiet conditions the relative amplitude dI/I of the TEC variations at night significantly exceeds the daytime relative amplitude. At high levels of magnetic field disturbances, the geomagnetic control of the amplitude of TEC variations at high and middle latitudes is much more significant than the regular diurnal variations. At the equatorial belt, on average, the amplitude of TEC variations in quiet and disturbed periods almost does not differ. The obtained results may be useful for development of the theory of ionospheric irregularities.     

基于相干度优化的极化顺轨干涉SAR慢小目标CFAR检测

为改善干涉合成孔径雷达（SAR）系统对慢动小目标的检测性能,研究了全极化顺轨干涉SAR（AT-POLINSAR）实现慢动目标恒虚警（CFAR）检测的方法。首先,通过对单基线AT-POLINSAR的系统设计,明确了其在现有技术条件下的可实现性,并对其信号形式与极化干涉回波进行了建模分析。然后,针对AT-POLINSAR 6维极化干涉矢量提出了以背景杂波平均相干度为优化准则的极化降维新方法,构建了一种统计分布类型与单极化干涉数据相同的次优极化标量干涉回波,从而使目前单极化顺轨干涉SAR（AT-INSAR）慢动目标CFAR检测方法可直接扩展至全极化情形。最后,通过检测实验对次优极化与单极化的慢动目标检测性能进行了对比分析。结果表明,次优极化方法能充分利用全极化信息提高INSAR对慢小目标的检测概率。      

The fast debris evolution model

The ‘particles-in-a-box’ (PIB) model introduced by Talent [Talent, D.L. Analytic model for orbital debris environmental management. J. Spacecraft Rocket, 29 (4), 508–513, 1992.] removed the need for computer-intensive Monte Carlo simulation to predict the gross characteristics of an evolving debris environment. The PIB model was described using a differential equation that allows the stability of the low Earth orbit (LEO) environment to be tested by a straightforward analysis of the equation’s coefficients. As part of an ongoing research effort to investigate more efficient approaches to evolutionary modelling and to develop a suite of educational tools, a new PIB model has been developed. The model, entitled Fast Debris Evolution (FADE), employs a first-order differential equation to describe the rate at which new objects ?10 cm are added and removed from the environment. Whilst Talent [Talent, D.L. Analytic model for orbital debris environmental management. J. Spacecraft Rocket, 29 (4), 508–513, 1992.] based the collision theory for the PIB approach on collisions between gas particles and adopted specific values for the parameters of the model from a number of references, the form and coefficients of the FADE model equations can be inferred from the outputs of future projections produced by high-fidelity models, such as the DAMAGE model.     

Sunspot numbers,interplanetary magnetic field,and cosmic ray intensity at earth: Nexus for the twentieth century

The pivotal role played by the interplanetary magnetic field (B) in modulating galactic cosmic ray (GCR) intensity in the heliosphere is described. We show that the inverse correlation observed by Forbush (1958) between GCRs and sunspot numbers (SSNs) is reflected in high correlation between SSNs and B (cc = 0.94). The SSN data are available since 1700 and the derived B data since 1835. The paleo-cosmic ray data are available for several millennia in the form of ¹⁰Be radionuclide sequestered in polar ice. The data of the ion chambers (ICs) at the Cheltenham–Fredericksburg–Yakutsk (CFY) sites are combined to create a data string for 1937–1988. In turn, these data are used to extend the measurements of the low energy GCR ions (>0.1 GeV) at balloon altitudes at high latitudes in Russia to 1937. These data are then correlated to B and the fit parameters are used to extend the low energy ion data to 1900, creating the instrumental era GCR time series for the twentieth century. The derived GCR time series is compared to ¹⁰Be measured at two sites in Greenland, namely Dye 3 and NGRIP for 1900–2000 to check the internal consistency of datasets for the long-term trend. We find that the annual mean rate (%) for 1965 at NGRIP is an outlier. We replace it with the mean of 1964 and 1965 rates and construct a new re-normalized time series at NGIP, improving the agreement with the derived instrumental era GCR time series for the twentieth century as well. This should encourage its use by heliophysics community for varied applications.     

Occurrence of equatorial plasma bubbles over Kolhapur

This paper reports the nightglow observations of OI 630.0 nm emissions, made by using all sky imager operating at low latitude station Kolhapur (16.8°N, 74.2°E and dip lat. 10.6°N) during high sunspot number years of 24th solar cycle. The images are analyzed to study the nocturnal, seasonal and solar activity dependence occurrence of plasma bubbles. We observed EPBs in images regularly during a limited period 19:30 to 02:30 LT and reach maximum probability of occurrence at 22:30 LT. The observation pattern of EPBs shows nearly no occurrence during the month of May and it maximizes during the period October–April. The equinox and solstice seasonal variations in the occurrence of plasma bubbles show nearly equal and large differences, respectively, between years of 2010–11 and 2011–12.