1998年9月24日地磁暴前宇宙线各向异性特性分析

分析了日本Nagoya 宇宙线闪烁体望远镜30°, 49°, 64° 倾角的东、西、南、北方向探测数据的变化特点, 运用小波分析方法定性地探讨了磁暴前后宇宙线南北、东西各向异性的变化特征. 研究发现, 当发 生大地磁暴时, 地面宇宙线强度的各向异性特征将发生非常大的变化, 这种变化一般在磁暴发生前10~20 h就开始出现. 当描述这种各向异性特征的各向异性指数的小波系数变化达到一定阈值时, 就可能有大地磁暴发生.      

质子地磁垂直截止刚度在2004年11月7-8日两个时刻处的数值模拟

利用单粒子轨道理论, 在T04和IGRF2000磁场模式建立的磁层模型基础上, 应用四阶龙格库塔方法, 模拟计算宇宙线带电质子在地球磁层中的运动以及沿天顶方向入射到达地球磁场内的某一特定位置, 得到了距地球表面450km高度处全球质子垂直截止刚度在2004年11月7-8日中两个时刻的计算值. 根据计算得到的684个不同位置处的截止刚度值, 分析了同一时刻地磁垂直截止刚度随磁纬和磁经的变化. 与此同时, 模拟计算了相同时刻下磁纬为30°, 磁经为0°, 45°, 90°, 135°, 180°, 225°, 270°, 315°处, 地磁有效垂直截止刚度随高度的变化情况. 结果表明, 在磁纬30°处, 其质子有效垂直截止刚度随距离地心高度的变化沿向阳处向磁尾处方向, 变化越来越缓慢.      

Normalized ionization yield function for various nuclei obtained with full Monte Carlo simulations

Several recent results important for production of ion pairs in the Earth atmosphere by various primary cosmic ray nuclei are presented. The direct ionization by various primary cosmic ray nuclei is explicitly obtained. The longitudinal profile of atmospheric cascades is sensitive to the energy and mass (charge) of the primary particle. In this study different cosmic ray nuclei are considered as primaries, namely Helium, Oxygen and Iron nuclei. The cosmic ray induced ionization is obtained on the basis of CORSIKA 6.52 code simulations using FLUKA 2006 and QGSJET II hadronic interaction models. The energy of the primary particles is normalized to GeV per nucleon. In addition, the ionization yield function Y is normalized as ion pair production per nucleon. The obtained ionization yield functions Y for various primaries are compared. The presented results and their application are discussed.     

A new dynamical atmospheric ionizing radiation (AIR) model for epidemiological studies.

A new Atmospheric Ionizing Radiation (AIR) model is currently being developed for use in radiation dose evaluation in epidemiological studies targeted to atmospheric flight personnel such as civilian airlines crewmembers. The model will allow computing values for biologically relevant parameters, e.g. dose equivalent and effective dose, for individual flights from 1945. Each flight is described by its actual three dimensional flight profile, i.e. geographic coordinates and altitudes varying with time. Solar modulated primary particles are filtered with a new analytical fully angular dependent geomagnetic cut off rigidity model, as a function of latitude, longitude, arrival direction, altitude and time. The particle transport results have been obtained with a technique based on the three-dimensional Monte Carlo transport code FLUKA, with a special procedure to deal with HZE particles. Particle fluxes are transformed into dose-related quantities and then integrated all along the flight path to obtain the overall flight dose. Preliminary validations of the particle transport technique using data from the AIR Project ER-2 flight campaign of measurements are encouraging. Future efforts will deal with modeling of the effects of the aircraft structure as well as inclusion of solar particle events.     

空间辐射场蒙特卡罗模拟计算方法研究

对利用蒙特卡罗方法对由银河宇宙射线引起的空间辐射场各成分进行计算的方法进行了调研,对计算模型的建立以及计算过程中通常使用的方差减小技术进行分析,给出了美国的Roesler等人利用FLUKA程序以及加拿大Anid等人利用MCNPX程序计算得到的由银河宇宙射线引起的空间辐射场各量值及其与实验结果的比较,验证了计算方法与计算模型的可靠性。对任意航线空间辐射场剂量分布预评估方法进行分析,给出了由银河宇宙射线引起的空间辐射场的基本特征。     

Distribution of secondary particles intensities over Earth’s surface: Effect of the geomagnetic field

We use the CORSIKA package (Heck et al., 1998) and AMS-01 flight data (Alcaraz et al., 2000) to evaluate the distribution of secondary particles in the Earth atmosphere. Distribution covers all longitudes and latitudes of STS-91 Space Shuttle flight trajectory to Mir Space Station. Moreover distribution covers all depth in the atmosphere in the evaluated area. We show distributions for e−, e+, μ+, μ−, gammas, hadrons and Cherenkov light from primary protons and helium component of cosmic rays flux. Our results compare favorably with other estimates made by different techniques.     

A new Mars radiation environment model with visualization.

A new model for the radiation environment to be found on the planet Mars due to Galactic Cosmic Rays (OCR) has been developed at the NASA Langley Research Center. Solar modulated primary particles rescaled for Mars conditions are transported through the Martian atmosphere, with temporal properties modeled with variable timescales, down to the surface, with altitude and backscattering patterns taken into account. The Martian atmosphere has been modeled by using the Mars Global Reference Atmospheric Model--version 2001 (Mars-GRAM 2001). The altitude to compute the atmospheric thickness profile has been determined by using a model for the topography based on the data provided by the Mars Orbiter Laser Altimeter (MOLA) instrument on board the Mars Global Surveyor (MGS) spacecraft. The Mars surface composition has been modeled based on averages over the measurements obtained from orbiting spacecraft and at various landing sites, taking into account the possible volatile inventory (e.g., CO2 ice, H2O ice) along with its time variation throughout the Martian year. Particle transport has been performed with the HZETRN heavy ion code. The Mars Radiation Environment Model has been made available worldwide through the Space Ionizing Radiation Effects and Shielding Tools (SIREST) website, a project of NASA Langley Research Center.     

Be-7 nuclei produced by galactic cosmic rays and solar energetic particles in the earth''s atmosphere

Be-7 radioactive nuclei with a half-life of 53.3 days result from spallation reactions of galactic cosmic rays(GCR) and solar energetic particles (SEP) with N and O nuclei in the Earth's atmosphere. We calculate the average global production of Be-7 in the atmosphere by GCR and SEP The result indicates that an intense SEP event produces a large amount of Be-7 in the polar stratosphere and part of them could be transported to the surface at lower latitudes. The ground-level measurement of Be-7 in Japan exhibits the possibility of enhancement in the Be-7 radioactivity associated with the intense SEP event on July 14, 2000. In addition, the present experiment shows seasonal variations in the surface Be-7 concentration which peaks in spring and autumn. We discuss the possible air mass mixing between the stratosphere and troposphere to explain the measured seasonal variations. The surface concentration of Pb-210 nuclei indicates a similar trend to that of Be-7 and we suggest two possible explanations.     

Ionization effect of solar protons in the Earth atmosphere – Case study of the 20 January 2005 SEP event

The cosmic ray ground level enhancement on January 20, 2005 is among the largest recorded events in the history of cosmic ray measurements. The solar protons of MeV energies cause an excess of ionization in the atmosphere, specifically over polar caps following major solar disturbances. The ionization effect in the Earth atmosphere is obtained for various latitudes on the basis of solar proton energy spectra, reconstructed from GOES 11 measurements and subsequent full Monte Carlo simulation of cosmic ray induced atmospheric cascade. The estimation of ionization rates is based on a numerical model for cosmic ray induced ionization. The evolution of atmospheric cascade is performed with the CORSIKA 6.52 code using FLUKA 2006b and QGSJET II hadron interaction models. The atmospheric ion rate ionization is explicitly obtained for 40°N, 60°N and 80°N latitudes. The time evolution of obtained ion rates is presented. It is demonstrated that ionization effect is negative for 40°N and small for 60°N, because of accompanying Forbush decrease. The ionization effect is significant only in sub-polar and polar atmosphere during the major ground level enhancement of 20 January 2005.     

Polarization imaging of dust cloud particles: Improvements and applications of the PROGRA2 instrument

The imaging system of the -PROGRA2 instrument allows to obtain maps of polarization and brightness of levitating dust clouds with a theoretical resolution of 10 μm per pixel. The measurements are conducted in microgravity during parabolic flights and on the ground by air-draught. It is then possible to measure the contribution of individual particles (grains, aggregates and agglomerates.) The size distribution can be retrieved, as well as the variation of polarization for a given phase angle with size for particles larger than 10 microns. Two different kinds of particles are considered: compact grains and (aggregates and agglomerates of) fluffy particles. Opposite results are obtained for these two kinds of particles, concerning the dependence of polarization with size and color in the visible domain for gray materials. These results, coupled with such remote sensing observations in the solar system, can then help to better understand the physical properties of solid particles and their variation in cometary comae, as well as in the Earth's atmosphere.     

Recent results form measurements of the energy spectrum of cosmic-ray induced neutrons aboard an ER-2 airplane and on the ground.

Crews of future high-altitude commercial aircraft may be significantly exposed to atmospheric cosmic radiation from galactic cosmic rays (GCR). To help determine such exposures, the Atmospheric Ionizing Radiation Project, an international collaboration of 15 laboratories, made simultaneous radiation measurements with 14 instruments on a NASA ER-2 high-altitude aircraft. The primary instrument was a sensitive extended-energy multisphere neutron spectrometer, which was also used to make measurements on the ground. Its detector responses were calculated for neutrons and charged hadrons at energies up to 100 GeV using the radiation transport code MCNPX. We have now recalculated the detector responses including the effects of the airplane structure. We are also using new FLUKA calculations of GCR-induced hadron spectra in the atmosphere to correct for spectrometer counts produced by charged hadrons. Neutron spectra are unfolded from the corrected measured count rates using the MAXED code. Results for the measured cosmic-ray neutron spectrum (thermal to >10 GeV), total neutron fluence rate, and neutron dose equivalent and effective dose rates, and their dependence on altitude and geomagnetic cutoff generally agree well with results from recent calculations of GCR-induced neutron spectra.     

Trapped antiprotons produced by cosmic rays in the Earth's magnetosphere.

The existence of significant fluxes of antiparticles in the Earth magnetosphere has been predicted on theoretical considerations in this article. These antiparticles (positrons or antiprotons) at several hundred kilometers of altitudes, we believe are not of direct extraterrestrial origin, but are the natural products of nuclear reactions of the high energy primary cosmic rays (CR) and trapped protons (TP) confined in the terrestrial radiation belt, with the constituents of terrestrial atmosphere. Extraterrestrial positrons and antiprotons born in nuclear reactions of the same CR particles passing through only 5-7 g/cm2 of interstellar matter, exhibit lower fluxes compared to the antiprotons born at hundreds of g/cm2 in the atmosphere, which when confined in the magnetic field of the Earth (in any other planet), get accumulated. We present the results of the computations of the antiproton fluxes at 10 MeV to several GeV energies due to CR particle interactions with the matter in the interstellar space, and also with the residual atmosphere at altitudes of approximately 1000 km over the Earth's surface. The estimates show that the magnetospheric antiproton fluxes are greater by two orders of magnitude compared to the extraterrestrial fluxes measured at energies <1-2 GeV.     

Preliminary evidence for drizzle in the middle cloud layer of Venus

Vertical trajectory data from the VEGA 1 and 2 balloon flights in the atmosphere of Venus is re-analyzed. A previously employed helium leak rate profile invoked to entirely account for the decrease in equilibrium float altitude of VEGA 2 is questioned and deemed ad hoc. As an alternative, it is proposed that both VEGA 1 and 2 experienced in-flight mass increases due to the deposition of cloud particles onto their envelopes, as well as losing helium at a reduced rate consistent with the pre-flight prediction. Particle deposition rates are estimated and found to be compatible with this alternative scenario. Possible evidence for drizzle is also presented. Preliminary experiments to derive aerosol deposition rate on a flat plate and the maximum feasible liquid mass that may be accumulated on a near-spherical envelope are briefly described. Further experimental work is recommended to constrain the deposition efficiency values involved and the maximum feasible drizzle fluxes that could have been encountered by both VEGA 1 and 2.     

Inner radiation belt source of helium and heavy hydrogen isotopes

Nuclear interactions between inner zone protons and atoms in the upper atmosphere provide the main source of energetic H and He isotopes nuclei in the radiation belt. This paper reports on the specified calculations of these isotope intensities using various inner zone proton intensity models (AP-8 and SAMPEX/PET PSB97), the atmosphere drift-averaged composition and density model MSIS-90, and cross-sections of the interaction processes from the GNASH nuclear model code. To calculate drift-averaged densities and energy losses of secondaries, the particles were tracked in the geomagnetic field (modelled through IGRF-95) by integrating numerically the equation of the motion. The calculations take into account the kinematics of nuclear interactions along the whole trajectory of trapped proton. The comparison with new data obtained from the experiments on board RESURS-04 and MITA satellites and with data from SAMPEX and CRRES satellites taken during different phases of solar activity shows that the upper atmosphere is a sufficient source for inner zone helium and heavy hydrogen isotopes. The calculation results are energy spectra and angular distributions of light nuclear isotopes in the inner radiation belt that may be used to develop helium inner radiation belt model and to evaluate their contribution to SEU (single event upset) rates.     

Time variations of cosmic-ray helium isotopes with BESS-Polar I

The Balloon-borne Experiment with a Superconducting Spectrometer (BESS) is configured with a solenoidal superconducting magnet and a suite of precision particle detectors, including time-of-flight hodoscopes based on plastic scintillators, a silica-aerogel Cherenkov detector, and a high resolution tracking system with a central jet-type drift chamber. The charges of incident particles are determined from energy losses in the scintillators. Their magnetic rigidities (momentum/charge) are measured by reconstructing each particle trajectory in the magnetic field, and their velocities are obtained by using the time-of-flight system. Together, these measurements can accurately identify helium isotopes among the incoming cosmic-ray helium nuclei up to energies in the GeV per nucleon region. The BESS-Polar I instrument flew for 8.5 days over Antarctica from December 13th to December 21st, 2004. Its long-duration flight and large geometric acceptance allow the time variations of isotopic fluxes to be studied for the first time. The time variations of helium isotope fluxes are presented here for rigidities from 1.2 to 2.5 GV and results are compared to previously reported proton data and neutron monitor data.     

A method for operative calculation of charged particle penetration into the LEO

We develop a simple method for calculating the effective vertical cutoff rigidity of charged particles, taking into account the Kp-index and the local time, on the basis of generalization of the results of extensive trajectory calculations for trial particles moving in the geomagnetic field. The vertical cutoff rigidities, calculated by the Tsyganenko-89 model, are presented as an International Geomagnetic Reference Field (IGRF) model calculated and thereafter corrected in accordance with the geomagnetic disturbance and local time conditions. The fits from the proposed method agree with the results of cutoff rigidity measurements carried out by satellites. The method is intended for applications using cutoff calculations, such as evaluating particle penetration of spatial boundaries, calculating magnetospheric transmissions for low-orbital spacecrafts flights and interpreting the results of orbital experiments.     

Neutron spectra in the atmosphere from interactions of primary cosmic rays

Spectra of neutrons from interactions of primary cosmic rays in the earth's atmosphere are calculated with the Monte Carlo model fluka for various depths down to sea level. We discuss the environmental models describing the primary cosmic ray spectrum and details of the calculations. Neutron energy spectra are presented for different depths in the atmosphere and for different geographical locations. By comparing results of calculations to measurements on neutron spectra it is shown that fluka may serve as an important tool for the estimation of the radiation environment in the atmosphere.     

The relationship between plasma effects and cosmic radiation studied with TriTel-LMP measurements during the ESEO mission

The main point of the paper is to use the simultaneous measurements of the energetic particle flux by TriTel and those of electron density by a Langmuir probe to study the question of to what extent solar electromagnetic and corpuscular radiation (galactic cosmic rays, particle precipitation from the radiation belts) are responsible for the ionization of the atmosphere. The electron density measured by the Langmuir probe is the sum of the ionization produced by the solar electromagnetic radiation and that due to the corpuscular radiation. The ionization produced by the solar electromagnetic radiation may be computed. The flux of energetic particles in an energy range may be determined by taking the difference between the threshold energy of the TriTel telescopes and the energy corresponding to the local cut-off rigidity. As the ESEO satellite will have a quasi-polar and circular orbit, the cut-off rigidity will change from low to high latitudes, thus enabling the assignment of different energy bands for the telescopes. Thus, it will be possible to determine which energy bands of particle produce ionization at different latitudes.     

Computation of radiation environment during ground level enhancements 65, 69 and 70 at equatorial region and flight altitudes

The radiation environment in the troposphere of the Earth is governed by cosmic rays of galactic and solar origin. During major solar energetic particles events the radiation environment changes dramatically. As a results the risk of biological effects due to exposure to ionizing radiation of aircrew increases. Here we present a numerical model for computation of absorbed dose in air due to cosmic rays of galactic and solar origin. It is applied for computation of radiation environment at flight altitude in the equatorial region during several major ground level enhancements, namely GLE65 on 28 October 2003, GLE69 on 20 January 2005 and GLE70 on 13 December 2006. The model is based on a full Monte Carlo simulation of cosmic ray induced atmospheric cascade. The cascade simulation is carried out with CORSIKA 6.990 code with corresponding hadron generators FLUKA 2011 and QGSJET II. The contribution of different cascade components, namely electromagnetic, hadron and muon is explicitly obtained. The spectra of arriving solar energetic particles are calculated from ground level measurements with neutron monitors and satellite data from GOES. The obtained results are discussed.     

The physics of the FLUKA code: Recent developments

FLUKA is a Monte-Carlo code able to simulate interaction and transport of hadrons, heavy ions and electromagnetic particles from few keV (or thermal neutron) to cosmic ray energies in whichever material. The highest priority in the design and development of the code has always been the implementation and improvement of sound and modern physical models. A summary of the FLUKA physical models is given, while recent developments are described in detail: among the others, extensions of the intermediate energy hadronic interaction generator, refinements in photon cross sections and interaction models, analytical on-line evolution of radio-activation and remnant dose. In particular, new developments in the nucleus–nucleus interaction models are discussed. Comparisons with experimental data and examples of applications of relevance for space radiation are also provided.