由环电流的宇宙线效应验证二维对称环电流模型

本文利用不同纬度超中子堆记录的宇宙线强度, 对不同磁暴类型下环电流的宇宙线效应进行了分析。加强的西向环电流使宇宙线截止刚度下降, 在某些特大的磁暴中, 中低纬字宙线截止刚度下降可达1GV以上, 从而使宇宙线强度增加。利用二维对称环电流模式, 从宇宙线的角度计算了环电流的磁暴效应。与实际观测值的对比显示, 除非特大的磁暴过程, 利用该模式所得结果与实际观测符合得较好。      

南极中山站宇宙线观测

2014年由两个具有0.5m×0.5m闪烁体面积的探测器组成的宇宙线μ子望远镜在南极中山站建成,观测数据通过卫星链路传回并发布.通过对中山站μ子望远镜数据的气象效应分析,发现中山站的气压变化对宇宙线计数率影响显著并呈负相关.计算得到宇宙线高能质子在中山站的垂直截止刚度R_C为0.076GV,对比中国西藏羊八井和北京高能质子的垂直截止刚度,中山站非常有利于对太阳高能质子事件的观测.      

雷暴期间次级宇宙线粒子强度瞬时变化研究

雷暴期间大气电场强度变化及其伴随的宇宙线粒子增长的研究, 对于理解大气电场对宇宙线次级粒子的加速机制具有极其重要的意义. 2006年4月至 8月期间, 西藏羊八井宇宙线观测站记录到了20多次雷暴事件. 分析了雷暴期间, ARGO实验scaler模式下次级宇宙线计数与大气电场之间的相关性. 结果显示, 雷暴期间大气电场剧烈变化时, 多重数n=1, 2的次级宇宙线计数率有明显增长, 增幅在1%～9%之间, 然而n=3, n≥ 4的次级宇宙线计数率增长不明显, 甚至没有增长. 该结果为进一步研究雷暴期间大气电场对次级宇宙线的加速机制打下了基础.      

两个临近脉冲型太阳能量电子事件的观测与拟合

WIND飞船在2010年11月15日观测到两个临近的脉冲型太阳能量电子事件, 这两个事件在1AU处呈现出不同类型的通量-时间曲线. 事件一的通量表现出快速上升及快速下降的特性; 事件二则表现为缓慢上升, 更缓慢下降, 事件二的持续时间是事件一的5～17倍. 以往的解释认为事 件二中的电子在行星际受到了更强的散射. 本文引入等腰三角形的释放函数并 求解电报方程, 利用得到的解对1AU处的观测进行拟合. 根据最佳拟合结果 反推两事件在太阳上的释放函数和在行星际传播的平均自由程, 发现造成两事 件在1AU处能谱和通量-时间曲线形状不同的原因是太阳上电子 加速过程的不同而非行星际散射. 结合SOHO卫星的白光观测, 发现两事件可能 与一个CME (日冕物质抛射)相关, 并进一步推测了这两个太阳电子事件可能 的加速图像.      

太阳宇宙线在电离层D层中的电离

本文根据带电粒子对D层大气电离的理论,导出了太阳宇宙线在D层的电子产生率Q(h)的表达式,并计算了不同级别的太阳宇宙线事件、不同能谱参数下,Q(h)在极区随高度的分布。结果表明,不同级别、不同能谱的太阳宇宙线事件在极区产生的电离有显著的差别。同一级别,能谱指数γ越大,在较高的高度上电子产生率越大;能谱指数越小,在较低的高度上电子产生率越大。电子产生率的分布曲线出现明显的双峰,一个峰位于60公里左右,另一个峰位于85公里左右。前一个峰主要由太阳宇宙线中质子产生的,后一个峰主要是z≥2的重粒子成分产生的。本文所得结果明显好于Velinov等人的结果。      

低平流层温度对银河系宇宙线强度变化的响应

选用芬兰中子堆监测站从lop年7月至1989年3月的宇宙线强度月平均曲线数据,和在相对应期间我国境内海拉尔等4个高空气象观测站的30hPa和100hPa等位势高度上常规气温观测月平均资料,经统计分析,取得初步结果:低平流层气温与银河系宇宙线辐射强度之间呈现负相关的特性;与银河系宇宙线辐射强度峰区和谷区相对应会在低平流层中引起幅值为2-3℃的温度变化;不同纬度观测站所得观测数据表明,上述增温幅值随纬度降低而减少。      

宇宙线变化特征在强地磁暴预报中的应用

利用宇宙线中子探测数据定性分析了地面宇宙线多台站之间的相互联系以及大磁暴与宇宙线之间的响应关系. 以Irkutsk和Oulu宇宙线台站为例, 运用小波去噪技术提高数据的稳定性. 结果表明, 相同世界时条件下, 两站宇宙线通量相关性在事件发生时较高; 而相同地方时条件下, 相关性则在平静期较高. 进一步采用相同地方时条件对不同宇宙线台站的通量在平静期和扰动期的相对变化进行分析, 选取2004年7月强地磁暴典型事例进行直观分析, 发现大地磁暴前Irkutsk和Oulu台站的宇宙线相对通量发生明显差异, 可以尝试作为强地磁暴宇宙线先兆特征. 通过对2001年3月至2005年5月的强磁暴和中强磁暴进行统计, 得到与强地磁暴相关的适当宇宙线相对差异阈值. 将得到的阈值对2005年9月至2011年12月所有强磁暴及中强磁暴进行验证, 总成功率达到87.5%, 误报率为35.7%, 结果较好.      

大气中子实时计算模式

为实时评估0~100km高度范围内的大气中子全球分布,对宇宙线在地磁场和大气中的传输过程进行了分析.利用蒙特卡罗方法工具包Geant4,预先计算不同能量的粒子在大气层中产生的次级粒子能谱分布,形成大气次级粒子数据库,并与相关模型进行对比,验证了该数据库的有效性和可靠性.以实测或预报的空间环境参数作为输入,计算同步轨道银河宇宙线和太阳质子事件能谱以及100km高度上的地磁垂直截止刚度,最终得到大气层顶上的粒子能谱.通过对大气次级粒子数据库的线性插值,实现1h分辨率的大气中子能谱和辐射剂量全球分布的实时计算.      

基于ACE飞船观测的银河宇宙线与太阳风变化的统计研究

基于ACE飞船的资料,通过时序迭加方法统计分析了最近两个太阳活动极小年时期（2007.0-2009.0和2016.5-2019.0年）的宇宙线计数与太阳风参数的关系.结果表明,宇宙线的计数受太阳风共转流相互作用区的强烈影响,宇宙线计数变化与快慢太阳风流界面的位置密切相关,例如流界面的穿越通常伴随着宇宙线计数的下降.分析表明,第一时段的具有“雪犁”效应的宇宙线计数下降对应于流界面附近的扩散系数急剧下降,而第二时段的非“雪犁”效应的计数下降可能是由穿越流界面后的扩散系数增大引起的.日球层电流片也与宇宙线计数变化存在一定的相关性,宇宙线粒子在日球层电流片附近存在一定程度的堆积.太阳风对宇宙线的作用机制表明,宇宙线的漂移和扩散效应决定了其在1AU附近的分布变化.      

磁暴期间环电流离子中性原子成像的模拟与观测比较研究

采用已经建立的环电流离子解析模型,结合Chamberlain地冕中性层模型,研究了2004年11月一次大磁暴期间的环电流区域中性原子(ENA)图像.结果表明,模拟的ENA图像与TC-2卫星搭载的中性原子成像仪(NUAUD)的观测图像在方位角或地方时分布、高度或纬度分布和能谱分布方面存在一定的差异.如果依据磁暴发展的不同阶段来选择环电流离子模型的方位角不对称因子和通量最大方向的方位角,增大地冕中性层在低高度区域的密度或者考虑氢(H)以外的其他中性成分,改进注入边界处的离子能谱分布函数,且考虑不同种类环电流离子的比例随磁暴发展可能发生的变化,该模型有望产生更符合观测的模拟ENA图像.     

Fluxes and nuclear abundances of cosmic rays inside the magnetosphere using a transmission function approach

At 1 AU and outside the Earth’s magnetosphere, the relative abundances to protons for He (He/p), C (C/p) and Fe (Fe/p) nuclei were calculated using the observation data of AMS-01 (for p and He) and HEAO-3 (for C and Fe) above 0.8 GeV/nucleon. In addition, the transmission function (TF) for the GCR propagation inside the magnetosphere was evaluated using the IGRF and T96 (introduced by Tsyganenko and Stern) models to obtain permitted and forbidden trajectories inside the magnetosphere. The TF allowed one to derive the primary He-nuclei fluxes in the same geomagnetic regions of AMS-01 observations. These fluxes were found in good agreement with the observation data. Furthermore inside the magnetosphere in addition to the flux of helium, it allowed one to obtain those of the primary p, C, and Fe nuclei and the relative abundances of He, C and Fe nuclei to protons from the same observation data of AMS-01 and HEAO-3 above ≈0.8 GeV/nucleon. Up to a geomagnetic latitude of ≈45.84°, the relative isotopic abundances were found to depend on the mass number I_isot and, on average, range from a factor ≈2.31 up to ≈3.35 larger than those outside the magnetosphere at 1 AU. Thus, the magnetospheric isotopic/nuclear relative abundances differ from those inside the solar cavity and those in the interstellar space. The usage of the TF approach can allow one to determine the nuclear abundances in the magnetosphere at any geomagnetic latitude and, thus, any orbit, provided that the CR spectra are determined at 1 AU.     

Effects of solar modulation on the cosmic ray positron fraction

We implemented a 2D Monte Carlo model to simulate the solar modulation of galactic cosmic rays. The model is based on the Parker’s transport equation which contains diffusion, convection, particle drift and energy loss. Following the evolution in time of the solar activity, we are able to modulate a local interstellar spectrum (LIS), that we assumed isotropic beyond the termination shock, down to the Earth position inside the heliosphere. In this work we focused our attention to the cosmic ray positron fraction at energy below ∼10 GeV, showing how the particle drift processes could explain different results for AMS-01 and PAMELA. We compare our modulated spectra with observations at Earth, and then make a prediction of the cosmic ray positron fraction for the AMS-02 experiment.     

The energy spectra of protons and helium measured with the ATIC experiment

The Advanced Thin Ionization Calorimeter (ATIC) balloon experiment is designed to investigate the composition and energy spectra of cosmic rays at the highest energies currently accessible by direct measurements, i.e., the region up to 100 TeV. The instrument consists of a silicon matrix for charge measurement, a graphite target (0.75 nuclear interaction length) to induce hadronic interactions, three layers of scintillator strip hodoscopes for triggering and trajectory reconstruction, and a Bismuth Germanate (BGO) crystal calorimeter (18 radiation lengths) to measure particle energies. ATIC has had two successful Long Duration Balloon (LDB) flights from McMurdo, Antarctica: one from 12/28/00 to 01/13/01 and the other from 12/29/02 to 01/18/03. We present the energy spectra of protons and helium extracted from the first flight, over the energy range from 100 GeV to 100 TeV, and compare them with the results from other experiments at both the lower and higher energies. ATIC-1 results do not indicate significant differences in spectral shape between protons and helium over the investigated energy range.     

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.     

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.     

Evaluation of solar proton spectra using balloon cosmic ray observations and Monte Carlo simulation results

We have developed a method to evaluate the spectrum of solar energetic protons at the top of the Earth’s atmosphere from the measurements of our balloon cosmic ray experiment. By using the Monte Carlo PLANETOCOSMICS code based on Geant4 we compute the interaction of solar protons [10 MeV–10 GeV] with the Earth’s atmosphere. We obtain the angular and energy distributions of secondary particles (p, e⁻, e⁺, photons, muons) at different atmospheric levels as a function of primary proton spectra. By comparing the calculated depth dependence of the particle flux with the data obtained by our balloon experiment we can deduce the parameters of the solar proton spectrum that best fit the observations. In this paper we discuss our solar proton spectrum estimation method, and present results of its application to selected solar proton events from 2001 to 2005.     

Statistical properties of SEP event flux declines

The interplanetary space is not a passive medium, which merely constitutes a scene for the propagation of previously accelerated energetic particles, but influences the distribution of particles by changing their energies as well due to interactions with magnetic field inhomogeneities. Such processes manifest themselves in the energy spectra of solar energetic particle (SEP) events. In this paper the fluxes of protons with energies of 4–60 MeV are investigated on the basis of two data sets. Both sets are homogeneous, obtained by the CPME instrument aboard the IMP 8 satellite between 1974 and 2001. The first includes all SEP events where the integral fluxes of >4 MeV protons exceeded 2 particle/cm² s sr. The other set consists of fluxes recorded in differential energy windows between 0.5 and 48 MeV. Important characteristics of SEP events include the rates of decrease of particle flux, which, as well as peak flux time, is an integral feature of the interplanetary medium within a considerable region, surrounding the observation point. The time intervals selected cover the decay phases of SEP events following flares, CMEs and interplanetary shocks of different origin. Only those parts of declines were selected, that could reasonably be described by exponential dependence, irrespective of the gradual/impulsive character of the events. It is shown that the average values of characteristic decay time, τ, and energy spectral index, γ, are all changing with the solar activity phase. Distributions of τ and γ values are obtained in SEPs with and without shocks and during different phases of events: just after peak flux and late after maximum.     

Resolving electrons from protons in ATIC

The Advanced Thin Ionization Calorimeter (ATIC) experiment is designed for high energy cosmic ray ion detection. The possibility to identify high energy primary cosmic ray electrons in the presence of the ‘background’ of cosmic ray protons has been studied by simulating nuclear-electromagnetic cascade showers using the FLUKA Monte Carlo simulation code. The ATIC design, consisting of a graphite target and an energy detection device, a totally active calorimeter built up of 2.5 cm × 2.5 cm × 25.0 cm BGO scintillator bars, gives sufficient information to distinguish electrons from protons. While identifying about 80% of electrons as such, only about 2 in 10,000 protons (@ 150 GeV) will mimic electrons. In September of 1999 ATIC was exposed to high-energy electron and proton beams at the CERN H2 beam line, and this data confirmed the electron detection capabilities of ATIC. From 2000-12-28 to 2001-01-13 ATIC was flown as a long duration balloon test flight from McMurdo, Antarctica, recording over 360 h of data and allowing electron separation to be confirmed in the flight data. In addition, ATIC electron detection capabilities can be checked by atmospheric gamma-ray observations.     

Peculiarities of the solar proton events of 19 October 1989 and 23 March 1991 according to the measurements onboard the Mir space station.

Flux and dose rate dynamics of solar cosmic rays were measured by the Lyulin dosimeter during the events 19 October 1989 and 23 March 1991. The maximum dose rate registered was 0.4, 0.12 and 0.01 cGy/hour, respectively. Based on the latitude distribution of particle flux a power law form for the energy spectra of solar protons in the anisotropic phase of the events on 19 October 1989 and 23 March 1991 was determined. It was obtained that after the development of geomagnetic storm protons with energies more than 1 GeV were registered.