由大气密度的季变化所造成的中、低层大气电离率的季变化

本文对CIRA 1972 模式下的中、低层大气密度季变化所造成的电离率的季变化作了理论计算, 利用宇宙线电离源函数对全球中、低层大气的电离率作了数值计算.计算结果表明, 在70km以下的大气电离率的季变化幅度为百分之几十.在电离峰值高度(约15km)以下幅度为约百分之十, 高度越高变幅越大, 到70km处达35%.在20km以上电离率的季变化趋势与大气密度的季变化趋势大致相同, 夏季比冬季大, 20km以下变化稍为复杂.      

Ion production and ionization effect in the atmosphere during the Bastille day GLE 59 due to high energy SEPs

The influence of high energy particles, specifically cosmic rays, on atmospheric physics and chemistry is highly discussed. In most of the proposed models the role of ionization in the atmosphere due to cosmic rays is not negligible. Moreover, effect(s) on minor constituents and aerosols are recently observed, specifically over the polar regions during strong solar particle events. According to the recent findings for such effects it is necessary an essential increase of ion production, specifically during the winter period. The galactic cosmic rays are the main source of ionization in the Earth’s stratosphere and troposphere. Occasionally, the atmospheric ionization is significantly enhanced during strong solar energetic particles events, specifically over the polar caps. During the solar cycle 23 several strong ground level enhancements were observed. One of the strongest was the Bastille day event occurred on 14 July 2000. Using a full Monte Carlo 3-D model, we compute the atmospheric ionization, considering explicitly the contribution of cosmic rays with galactic and solar origin, focusing on high energy particles. The model is based on atmospheric cascade simulation with the PLANETOCOSMICS code. The ion production rate is computed as a function of the altitude above the sea level. The ion production rate is computed on a step ranging from 10 to 30?min throughout the event, considering explicitly the spectral and angular characteristics of the high energy part of solar protons as well as their time evolution. The corresponding event averaged ionization effect relative to the average due to galactic cosmic rays is computed in lower stratosphere and upper troposphere at various altitudes, namely 20?km, 15?km, 12?km and 8?km above the sea level in a sub-polar and polar regions. The 24^h and the weekly ionization effects are also computed in the troposphere and low stratosphere. Several applications are discussed.     

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.     

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.     

Model for induced ionization by galactic cosmic rays in the Earth atmosphere and ionosphere

The ionization profiles produced by galactic cosmic rays in the Earth atmosphere and ionosphere are obtained on the basis of Monte Carlo simulations. Cascade processes in the atmosphere are simulated using CORSIKA 6.52 code with FLUKA 2006 and QGSJET II hadronic interaction subroutines. Proton induced showers are considered using a realistic atmospheric model (US Standard Atmosphere). The energy deposit from different components is taken into account, namely electromagnetic, hadron and muon components. The curvature of the atmosphere is considered in the computer code. On the basis of the computed ionization yield function the ion pair production rate in the atmosphere is obtained for different conditions and locations. The model is applicable to the entire atmosphere, from ground level to upper atmosphere. Several applications of the obtained results are discussed. The Monte Carlo simulation model considers nuclear interactions below the altitude of 35 km. It is compared with analytical–numerical electron production rate model. The latter model which takes into account the electromagnetic interactions above altitudes of 35 km has two main regions of application: above 50 km (thin target model) and between 35 and 50 km (intermediate target model). A good agreement between the CORSIKA results and analytical–numerical model results is found above altitude of 35 km.     

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

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

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.     

The atmospheric cosmic- and solar energetic particle radiation environment at aircraft altitudes

Galactic cosmic rays interact with the solar wind, the earth's magnetic field and its atmosphere to produce hadron, lepton and photon fields at aircraft altitudes. In addition to cosmic rays, energetic particles generated by solar activity bombard the earth from time to time. These particles, while less energetic than cosmic rays, also produce radiation fields at aircraft altitudes which have qualitatively the same properties as atmospheric cosmic rays. We have used a code based on transport theory to calculate atmospheric cosmic-ray quantities and compared them with experimental data. Agreement with these data is seen to be good. We have then used this code to calculate equivalent doses to aircraft crews. We have also used the code to calculate radiation doses from several large solar energetic particle events which took place in 1989, including the very large event that occurred on September 29^th and 30^th of that year. The spectra incident on the atmosphere were determined assuming diffusive shock theory.     

Observations and predictions of secondary neutrons on Space Shuttle and aircraft.

The Cosmic Radiation Effects and Activation Monitor has flown on six Shuttle flights between September 1991 and February 1995 covering the full range of inclinations as well as altitudes between 220 and 570 km, while a version has flown at supersonic altitudes on Concorde between 1988 and 1992 and at subsonic altitudes on a SAS Boeing 767 between May and August 1993. The Shuttle flights have included passive packages in addition to the active cosmic ray monitor which comprises an array of pin diodes. These are positioned at a number of locations to investigate the influence of shielding and local materials. Use of both metal activation foils and scintillator crystals enables neutron fluences to be inferred from the induced radioactivity which is observed on return to Earth. Supporting radiation transport calculations are performed to predict secondary neutron spectra and the energy deposition due to nuclear reactions in silicon pin diodes and the induced radioactivity in the various scintillator crystals. The wide variety of orbital and atmospheric locations enables investigation of the influence of shielding on cosmic ray, trapped proton and solar flare proton spectra.     

大气中子实时计算模式

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

Applications and usage of the real-time Neutron Monitor Database

A high-time resolution Neutron Monitor Database (NMDB) has started to be realized in the frame of the Seventh Framework Programme of the European Commission. This database will include cosmic ray data from at least 18 neutron monitors distributed around the world and operated in real-time. The implementation of the NMDB will provide the opportunity for several research applications most of which will be realized in real-time mode. An important one will be the establishment of an Alert signal when dangerous solar cosmic ray particles are heading to the Earth, resulting into ground level enhancements effects registered by neutron monitors. Furthermore, on the basis of these events analysis, the mapping of all ground level enhancement features in near real-time mode will provide an overall picture of these phenomena and will be used as an input for the calculation of the ionization of the atmosphere. The latter will be useful together with other contributions to radiation dose calculations within the atmosphere at several altitudes and will reveal the absorbed doses during flights. Moreover, special algorithms for anisotropy and pitch angle distribution of solar cosmic rays, which have been developed over the years, will also be set online offering the advantage to give information about the conditions of the interplanetary space. All of the applications will serve the needs of the modern world which relies at space environment and will use the extensive network of neutron monitors as a multi-directional spectrographic detector. On top of which, the decreases of the cosmic ray intensity – known as Forbush decreases – will also be analyzed and a number of important parameters such as galactic cosmic ray anisotropy will be made available to the users of NMDB. A part of the NMDB project is also dedicated to the creation of a public outreach website with the scope to inform about cosmic rays and their possible effects on humans, technological systems and space-terrestrial environment. Therefore, NMDB will also stand as an informative gate on space research through neutron monitor’s data usage.     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.     

Interannual and seasonal variations of diurnal tide,gravity wave,ozone, and water vapor as observed by MLS during 1991–1994

The diurnal tide in the mesosphere and lower thermosphere (MLT) shows large seasonal and interannual variations. Despite recent modeling investigations, the underlying physical mechanisms for causing these variations remain unclear. This paper provides further observational constraints to tide-sensitive variables (H₂O, O₃, and gravity wave variances) used by the models, which are obtained simultaneously by upper atmosphere research satellite microwave limb sounder at altitudes below the MLT region. The strong quasi biannual oscillation and semiannual oscillation variations in these measurements reveal good correlations between the diurnal tide with other tide-sensitive variables, which should be taken into account for further modeling studies.     

非相干雷达探测低电离层的结果分析

位于波多黎各的Arecibo非相干雷达可以获得低电离层电子和离子密度, 利用此非相干雷达数据对中纬度低电离层的运动特征进行研究. 得到了电子密度随时间和高度的变化 情况, 结果显著呈现出周日变化特征, 并分析了电子密度随高度的变化规律. 进一步对数据进行频谱分析, 深入研究低电离层电子密度的周日变化效应. 得到电子密度的高度剖面, 发现从F层底部到E层有明显的等离子体沉降. 低电离层的层结构特征及电子密度变化表明, 在该区域还存在不同程度的等离子体扰动, 由此对低电离层的作用因素 进行分析, 认为大气潮汐或声重波可能对低电离层产生扰动, 即低电离层与大气存在一定程度的耦合作用.      

Preliminary results of fair-weather atmospheric electric field in the proximity of Main Boundary Thrust,Northern Pakistan

The atmospheric electric Potential Gradient (PG) variation of fair-weather days at Muzaffarabad (MZF, Northern Pakistan) station is presented for the period of January 2015–October 2017. The present investigations focus on the diurnal variation of atmospheric PG on the seasonal and annual time scale. The fair-weather seasonal and annual average PG variation has revealed two peaks, i.e. a primary maximum peak and a secondary maximum peak. The average maximum PG value found is ～410?Vm^?1. The results are then compared with the well-known standard oceanic Carnegie curve. The diurnal curve of PG is found to deviate from the Carnegie curve. The seasonal PG variation at MZF shows lower values during the summer (monsoon) and autumn (post-monsoon) as compared to the winter and spring (pre-monsoon) which could be due to local aerosol concentration. Global comparison of PG (%) values of the annual mean at different longitude, with MZF observatory, shows local morning and evening peaks depicting local sunrise and sunset effects on the PG. The overall outcomes will certainly contribute to further investigate the Global Electric Circuit (GEC).     

Long-term (50 years) measurements of cosmic ray fluxes in the atmosphere

Since the middle of 1957 till present time the group of researchers of P.N. Lebedev Physical Institute of the Russian Academy of Sciences has carried out the regular balloon borne measurements of charged particle fluxes in the atmosphere. The measurements are performed at polar (northern and southern) and middle latitudes and cover the interval of heights from the ground level up to 30–35 km. Standard detectors of particles (gas-discharged counters) have been used. More than 80,000 measurements of cosmic ray fluxes in the atmosphere have been performed to the present time. In the data analysis the geomagnetic field and the Earth’s atmosphere are used as cosmic ray spectrometers.     

日全食期间中间层大气臭氧光化过程试验

利用包括活O_x、HO_x、和NO_x等成份在内的一维时变光化模式,研究了日全食时上中层大气臭氧和O₂(1△g)的变化特性,并通过与观测资料的比较进行了光化模式的试验.结果表明,日全食时太阳辐照变化对1.27μm气晖影响的高度比曙暮期间明显偏低,利用1.27μm气晖光度计测量大气臭氧的高度范围约为52-67km.1997年3月9日漠河日全食时测量的大气臭氧含量变化幅度为50％.这与考虑O_x、HO_x、和NO_x等成份的光化模式理论模拟结果基本一致.它表明在52-67km高度范围不存在臭氧亏损现象.      

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.     

大气中层顶区域波相互作用的一个观测个例

利用SOUSY VHF雷达的观测数据分析了极区中层顶83.4-91.2km范围内大气风场波动的非线性相互作用。大气风场的谱在不同高度上均有明显的潮汐分量峰值，纬向风分量中35h波、半日潮和8.9h惯必重力波构成共振相互作用对，经向风分量中33h波、半日潮和19h惯性重力波构成共振相互作用对。双谱分析表明，这些共振对在许多高度上都发生耦合，35h或33h波振幅的极小值与半日潮的极大值出现的高度几乎相同，呈现出明显的非线性相互作用在空间上不是局域的，而是存在于中层顶区域的几乎所有高度上，这种相互作用不仅导致半日潮振幅随时间的变化，也使半日潮的振幅随空间变化。35h和33h波动可能是在其他时段或其他位置通过行星波与周日潮相互作用产生的，然后传播到观测点并与半日潮发生相互作用。