磁层扰动对同步高度太阳质子截止刚度的影响

采用Williams—Mead磁场模式,考虑当行星际磁场存在南向分量时的磁层扰动情况,计算了一组同步卫星所在位置太阳质子截止刚度随地方时的变化。计算结果表明,当考虑行星际磁场存在5γ南向分量时,截止刚度值显著降低。午夜西方截止刚度极小值由原来的3.1兆电子伏降至1.9兆电子伏,中午西方截止刚度极大值由原来的22.1兆电子伏降至16.6兆电子伏。与同步高度太阳质子的观测结果更加接近。     

磁壳参数L与磁暴Dst指数和行星际条件的关系

用磁坐标L-/A来描述地球近地空间粒子特性和卫星位置是近年来空间物理探测研究和数据分析中的一个新趋势.利用T96磁场模型计算了L值,并比较了在地球磁层剧烈活动期间和不同行星际条件下,用偶极子模型,国际参照磁场(IGRF)模型,和T96磁场模型这三种方式计算得到的地球表面L-A磁坐标之间的区别.在地磁纬度大于30°时偶极子近似和IGRF磁场模型计算得到的L值差别开始增大.在地磁纬度大于50°时,用IGRF磁场模型和T96磁场模型计算得到的L值差别开始增大.由于T96磁场模型引入了行星际磁场南北和东西分量,计算的L值包含了行星际条件的影响,并具有了随地方时变化的特性.本工作对于辐射带粒子动态模式的建立,以及正确理解卫星磁坐标位置等具有重要意义.     

基于NOAA卫星数据的辐射带质子投掷角分布研究

利用NOAA-12卫星数据对空间环境平静时期太阳同步轨道处辐射带质子投掷角分布进行了研究. 根据投掷角分布的经验公式,计算出90°投掷角的质子方向强度和各向异性指数n. 质子投掷角分布按n的取值范围可分为三类,即90°峰值分布、平顶分布和蝴蝶形分布. 观测证实,对于辐射带质子,三种分布类型均存在并且具有明显的空间区域特征. 在内辐射带边缘地区90°峰值分布占主要优势;在外辐射带高L值区域,90°峰值分布明显减少,平顶分布和蝴蝶形分布逐渐增多. 针对90°峰值分布,研究了质子强度各向异性的区域分布特征,对于内辐射带区域,n值随L值的增大而增大,对于外辐射带,n值表现为逐渐下降的趋势. 为了研究质子投掷角分布对磁地方时的依赖关系,分析了能量为250～800keV的质子在两个不同磁地方时范围的投掷角分布规律. 结果显示,在内辐射带,质子强度的投掷角分布相对稳定,随磁地方时的变化并不显著;而在外辐射带的高L值区域,质子强度的投掷角分布随磁地方时变化明显,与磁地方时之间有明显的依赖关系.      

磁层-电离层电动耦合与中纬地磁指数的变化

本文探讨磁层一电离层耦合过程内中纬地磁指数的变化特点,并与极光电集流和赤道电集流(指数)变化相比较。相关分析和时序叠加分析均表明,高、中、低纬地磁指数变化可归结为磁层一电离层电动耦合的统一物理图象。有R事件的磁暴主相初期和无R事件的磁扰期内,赤道电集流和中纬地磁指数的变化各不相同。这再次证明,耦合分析中将磁层源扰动的直接穿透作用与经电离层内动力过程的效应二者加以区分和综合研究是很重要的。      

利用HRDI/UARS资料分析东亚区域中层大气纬向风气候特征

利用美国高层大气研究卫星(UARS)搭载的高分辨率多普勒测风仪(HRDI)获得的中层大气风场观测资料,对东亚区域中层大气纬向风的垂直分布与变化特征进行了分析研究.多年平均结果显示东亚区域中层大气纬向风具有显著的区域特征,与当前普遍使用的参考大气CIRA-86相比存在显著的不同.在冬季,东亚区域中间层西风急流中心位于25°-35°N之间的75 km高度,与CIRA-86相比,该中心纬度偏南5°,高度偏高10 km;在秋季,东亚区域低热层高度存在一个显著的从赤道到高纬度的东风带,而CIRA-86不存在.分析结果还表明,除了夏季中纬度地区,在东亚区域上空中高层大气各高度上均存在相当显著的区域尺度扰动结构.在热带,低热层高度纬向风无论冬夏,沿纬圈方向都表现出相当显著的不均匀性,夏季这种不均匀性进一步向下扩展到55 km高度.与上述热带扰动特征相比,中纬度地区夏季的纬向风在各个高度沿纬圈相当均匀,但是在冬季,中间层和低热层高度都存在沿纬圈方向显著的纬向风扰动结构.      

低纬(海南)电离层等离子体漂移对地磁活动的响应

利用2003-2016年期间子午工程海南站（19.5°N,109.1°E）数字测高仪观测到的电离层等离子体漂移数据,分析了高低两种太阳活动条件下纬向和垂直向漂移对近磁静、中等磁扰和强磁扰三种地磁活动水平的响应特性.结果表明:日间纬向漂移各季节均以西向为主,随地磁活动无明显变化,白天日出附近和夜间漂移在各季节均以东向为主,随地磁活动增强而减弱,减弱程度在分季最大,在夏季最小;日间垂直漂移在零值附近变化,且不受地磁活动和季节影响,日落附近漂移仅在分季受到地磁活动的抑制,午夜前垂直漂移在分季受到抑制,在冬季因强磁扰而反向,夏季无明显规律,子夜至日出后垂直漂移在各季节随地磁活动增强而减小.与赤道区Jicamarca相比,两地漂移对地磁活动的响应相近,但在幅度和相位上存在差异,这可能是两地区的地理位置、背景电场和风场结构等不同造成的.      

磁层－电离层电动耦合与中纬地磁指数的变化

本文探讨磁层一电离层耦合过程内中纬地磁指数的变化特点，并与极光电集流和赤道电集流（指数）变化相比较．相关分析和时序叠加分析均表明，高、中、低纬地磁指数变化可归结为磁层一电离层电动耦合的统一物理图象．有Ｒ事件的磁暴主相初期和无Ｒ事件的磁扰期内，赤道电集流和中纬地磁指数的变化各不相同．这再次证明，耦合分析中将磁层源扰动的直接穿透作用与经电离层内动力过程的效应二者加以区分和综合研究是很重要的．     

南大西洋异常区高能质子通量南北向差异

利用2012年8月NOAA N15,N16,N17,N18卫星中能质子和电子探测器（MEPED）>6.9MeV质子能档的数据,研究了800km高度处南大西洋异常区（SAA）质子分布的南北方向差异.数据显示,质子计数率沿经线随纬度的变化近似满足高斯分布,且向南运动的质子多于向北运动的质子.MEPED具有0°和90°两个探测方向,可对质子投掷角分布进行估算,从而消除在探测南向和北向质子时探测器方向角偏差所造成的影响.结果表明,在经度40°W、纬度13°S至23°S位置上,仍然存在南向质子比北向质子多20%~30%的现象.其原因可能与地磁异常使南向质子镜像点更低,更多质子散射进入大气层沉降有关.分析还发现,南大西洋异常区质子通量存在地方时变化.日侧比夜侧强约20%,这可能是由于地磁场日变化引起的.      

通量传输事件相关的动力学磁场湍流观测研究

THEMIS卫星观测到通量传输事件（FTE）的同时,也在磁层侧涡流区域观测到强磁场扰动现象.利用快速傅里叶变换分析磁场扰动频谱特征发现:大约在FTE的扰动频率（约0.1Hz）处,功率谱密度达到峰值;在质子回旋频率（约1Hz）至64Hz的频段内,功率谱密度随着频率的增大而减小,服从幂律分布P₀ f-α.因此,可以认为这些磁场扰动为低纬边界层中的动力学磁场湍流.研究结果表明,当低纬边界层（Low Latitude Boundary Layer,LLBL）中卫星相对磁层顶或FTE的位置越来越远时,功率谱密度与功率谱斜率α（幂律指数）降低,但FTE所在的方位角或低纬磁层顶的磁地方时对幂律指数α和功率谱密度没有显著影响.这些观测特征表明移动的FTE是磁场湍流的源.磁层顶上的大规模扰动（如FTE）和相关的磁场湍流从动力学尺度揭示了磁鞘与磁层的类黏滞相互作用.然而低纬边界层中FTE磁层侧涡流形成所需的黏滞性是否可由磁场湍流来提供还需要验证.      

基于嫦娥一号高能粒子数据的地球磁层屏蔽效应研究

月球绕地球运行轨道约有1/4位于地球磁层内,因此,地球磁层是否会为月球轨道附近高能粒子提供足够的磁场屏蔽对于探索月球活动具有重要影响.嫦娥一号是中国首颗绕月人造卫星,其绕月飞行的工作轨道距离月球表面200 km.通过对嫦娥一号高能粒子探测器(HPD)的探测数据进行分析,比较了当月球位于地球磁层内外6个不同能道(能量范围4～400 MeV)时质子通量的变化,发现当月球位于地球磁层内时,这些能道的质子通量并没有发生显著减少,结果表明地球磁层不能为月球轨道附近高能粒子提供显著的磁屏蔽.     

The local time dependence of the anisotropic solar cosmic ray flux.

The distribution of the solar cosmic radiation flux over the earth is not uniform, but the result of complex phenomena involving the interplanetary magnetic field, the geomagnetic field and latitude and longitude of locations on the earth. The latitude effect relates to the geomagnetic shield; the longitude effect relates to local time. For anisotropic solar cosmic ray events the maximum particle flux is always along the interplanetary magnetic field direction, sometimes called the Archimedean spiral path from the sun to the earth. During anisotropic solar cosmic ray event, the locations on the earth viewing "sunward" into the interplanetary magnetic field direction will observe the largest flux (when adjustments are made for the magnetic latitude effect). To relate this phenomena to aircraft routes, for anisotropic solar cosmic ray events that occur during "normal quiescent" conditions, the maximum solar cosmic ray flux (and corresponding solar particle radiation dose) will be observed in the dawn quadrant, ideally at about 06 hours local time.     

The Tien-Shan mountain cosmic ray station of the Ionosphere Institute of Kazakhstan Republic

The mountain cosmic ray (CR) station of the Ionosphere Institute of Kazakhstan Republic (CR station Alma-Ata B, 43.1N latitude, 76.6E longitude, geomagnetic rigidity cutoff 6.69 GV) is a center for an experimental study of the non-stationary processes caused by cosmic rays of different origin in the interplanetary and near-Earth space, so as their influence on the state of the Earth’s magnetosphere and the upper atmosphere layers. This paper summarizes efforts performed over the years by the Almaty CR group till its present status.     

On the possibility to check the magnetosphere’s model by CR: The strong geomagnetic storm in November 2003

We report on a study of cosmic ray cutoff rigidity variations during the strong geomagnetic storm of 18–24 November 2003. We employed the most recent Tsyganenko magnetospheric model to represent the very strong disturbed Magnetosphere. We used this magnetic field for the cosmic ray trajectory calculations to determine the geomagnetic cutoff rigidity throughout this period of severe geomagnetic disturbance. We determine the cutoff rigidity changes during this period by two methods, by trajectory calculations and by the spectrographic global survey method. The values of the change of cutoff rigidities obtained by two different methods are correlated with the D_st and interplanetary magnetic field and plasma parameters and result in correlation coefficients in the range 0.63–0.84 for the various cosmic ray stations. The result of this study indicates that the most significant contributions to the cutoff rigidity changes are due to D_st variation although the influence of solar wind density and B_z and B_y components of IMF variations is significant.     

World grid of cosmic ray vertical cut-off rigidity for the last decade

Cosmic ray cut-off rigidity tables and maps over the world concerning the epochs 2010, 2015 and the current one 2020 have been constructed. These maps display the effective cut-off rigidity in every five degrees in latitude and in longitude at the altitude of 20 km above the surface of the international reference ellipsoid. The values of the geomagnetic cut-off rigidity were calculated in every 5° in latitude and in every 15° in longitude applying the well-known method of particle trajectory calculations resulted from the theory of the particle motion in the Earth's magnetic field. The applied software employed the 12th Generation of the International Geomagnetic Reference Field (IGRF 12) and trajectories were calculated at 0.01 GV intervals in order to determine the vertical cut-off rigidity for each location. Beyond the use of the calculated cut-off rigidity values as a basic reference of charged particle access to different geographical locations during quiet and/or more intense geomagnetic periods, these results can be used for a long- term forecasting of the geomagnetic conditions variations.     

Preliminary validation of computational procedures for a new atmospheric ionizing radiation (AIR) model.

A new computational procedure to determine particle fluxes in the Earth's atmosphere is presented. The primary cosmic ray spectrum has been modeled through an analysis of simultaneous proton and helium measurements made on high altitude balloon flights and spacecraft. An improved global fit to the data was achieved through applying a unique technique utilizing the Fokker-Plank equation with a non-linear rigidity-dependent diffusion coefficient. The propagation of primary particles through the Earth's atmosphere is calculated with a three-dimensional Monte Carlo transport program called FLUKA. Primary protons and helium nuclei (alphas) are generated within the rigidity range of 0.5 GV-20 TV uniform in cos2 theta. For a given location, primaries above the effective cutoff rigidity are transported through the atmosphere. Alpha particles are initially transported with a separate package called HEAVY to simulate fragmentation. This package interfaces with FLUKA to provide interaction starting points for each nucleon originating from a helium nucleus. Results from this calculation are presented and compared to measurements.     

Variations in cosmic ray cutoff rigidities during the great geomagnetic storm of November 2004

A very strong interplanetary and magnetospheric disturbance observed on 7–13 November 2004 can be regarded as one of the strongest events during the entire period of space observations. In this paper we report on the studies of cosmic ray cutoff rigidity variations during 7–13 November 2004 showing how storm conditions can affect the direct cosmic ray access to the inner magnetosphere. Effective cutoff rigidities have been calculated for selected points on the ground by tracing trajectories of cosmic ray particles through the magnetospheric magnetic field of the “storm-oriented” Tsyganenko 2003 model. Cutoff rigidity variations have also been determined by the spectrographic global survey method on the basis of experimental data of the neutron monitor network. Relations between the calculated and experimental cutoff rigidities and the geomagnetic Dst-index and interplanetary parameters have been investigated. Correlation coefficients between the cutoff rigidities obtained by the trajectory tracing method and the spectrographic global survey method have been found to be in the limits 0.76–0.89 for all stations except the low-latitude station Tokyo (0.35). The most pronounced correlation has been revealed between the cutoff rigidities that exhibited a very large variation of ∼1–1.5 GV during the magnetic storm and the Dst index.     

Cosmic ray anisotropies in the outer heliosphere

The observation of the directional distribution of energetic and cosmic ray particles has been done with the Voyager spacecraft over a long period. Since 2002, when the first flux enhancements of charged particles associated with the approach of Voyager 1 to the solar wind termination shock were observed, these anisotropy measurements have become of special interest. They play an important role to understand the magnetic field and shock structure and the basics of the modulation of cosmic ray and anomalous particles at and beyond the termination shock. They also serve as motivation to study the spatial behavior of galactic and anomalous cosmic ray anisotropies with numerical modulation models in order to illustrate how the radial anisotropy, at different energies, change from upstream to downstream of the termination shock. Observations made by Voyager 1 indicate that the termination shock is a complicated region than previously thought, hence the effects of the latitude dependence of the termination shock’s compression ratio and injection efficiency on the radial anisotropies of galactic and anomalous protons will be illustrated. We find that the magnitude and direction of the radial anisotropy strongly depends on the position in the heliosphere and the energy of particles. The effect of the TS on the radial anisotropy is to abruptly increase its value in the heliosheath especially in the A > 0 cycle for galactic protons and in both polarity cycles for anomalous protons. Furthermore, the global effect of the latitude dependence of the shock’s compression ratio is to increase the radial anisotropy for galactic protons throughout the heliosphere, while when combined with the latitude dependence of the injection efficiency this increase depends on modulation factors for anomalous protons and can even alter the direction of the radial anisotropy.     

Cosmic ray exposure factors for Shuttle altitudes derived from calculated cut-off rigidities.

The allowed cosmic radiation flux accessible to an earth-orbiting spacecraft is a complex function of the satellite position and the geomagnetic cutoff characteristics at each zenith and azimuth angle at each position. We have determined cosmic ray exposure factors for the galactic cosmic ray spectrum for typical shuttle altitudes and inclinations up to 50 degrees. We have utilized d world grid of trajectory-derived cutoff rigidity calculations at 400 km altitude to determine geomagnetic transmission functions that permit a simple and direct calculation of the allowed cosmic ray spectrum to a 400 km satellite orbit. If the interplanetary cosmic ray spectrum is multiplied by the orbit-averaged geomagnetic transmission function the result is the allowed cosmic ray spectrum at the spacecraft.     

大气中子实时计算模式

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