第23太阳活动周热层大气密度对太阳辐射指数F10.7响应的模拟研究

利用NCAR-TIEGCM计算了第23太阳活动周期间（1996—2008年）400km高度上的大气密度，并统计分析大气密度对太阳辐射指数FF_10.7的响应.结果表明，在第23太阳活动周内，大气密度的变化趋势与太阳辐射指数FF_10.7的变化趋势基本一致，但是大气密度在不同年份、不同月份对太阳辐射指数FF_10.7的响应存在差异.第23太阳活动周内太阳辐射极大值和极小值之比大于4，而大气密度的极大值与极小值之比则大于10.太阳辐射低年的年内大气密度变化不到2倍，而太阳辐射高年的年内大气密度变化可达2倍甚至3倍.大气密度与FF_10.7指数在北半球高纬的相关系数比南半球高纬的相关系数大.在低纬地区，太阳辐射高年大气密度与FF_10.7指数的相关系数比低年的大.不同纬度上，大气密度与太阳辐射指数FF_10.7的27天变化值之间的相关系数都大于其与81天变化值之间的相关系数.      

太阳活动与热层大气密度的相关性研究

为分析太阳活动对热层大气的影响,使用250km,400km,550km高度处热层大气密度与太阳F_10.7指数数据,研究了二者的周期变化及相关关系. 结果表明,热层大气密度的变化与太阳活动呈现相似的变化趋势;两者均具有显著的27天及11年周期变化特征,热层大气密度还存在7～11天及0.5年和1年的变化特征,且高度越高越明显;热层大气密度对太阳活动的最佳响应滞后为3天,无论何种地磁活动水平下,400km高度处相关性高于250km,550km处相关性最小,且太阳活动下降相期间高于上升相;250km,400km和550km高度处热层大气密度和太阳活动的统计结果分别为饱和、线性和放大关系;高度越高的热层大气密度对太阳活动响应越敏感.      

550～600 km高度上热层大气密度的中长期变化

选用了2005年8月20日至2006年7月28日高度550～600 km附近的热层大气密度探测数据,对表征太阳活动的F10.7值和表征地磁活动强度的Ap指数进行了相关特性的统计.分析结果表明,在无明显地磁扰动时热层大气密度日平均值的涨落呈现27日和准半年的周期性变化,但在地磁扰动期间这种变化的周期性会被削弱,且大气密度的周日变化幅度与F10.7值呈正相关关系.      

LSTM神经网络在太阳F10.7射电流量中期预报中的应用

F_10.7指数作为大气密度经验模型的重要输入参量，其预报精度直接影响航天器轨道预报精度.研究发现，太阳活动表现出长时间尺度上平均11年和中短时间尺度平均27天的周期性变化特征.依据这一观测事实，基于长短期记忆单元（Long Short-term Memory，LSTM）递归神经网络方法进行F_10.7指数未来27天的中期预报.利用一个连续长时段F_10.7数据作为训练数据，构建LSTM神经网络训练和预测模型，分别预测太阳活动高低年未来27天的F_10.7指数.结果表明，太阳活动高年的第27天F_10.7指数预报平均相对误差最优可达10%以内，低年最优可达2%以内.      

强磁暴、能量粒子暴与热层大气密度涨落之间的相关关系

利用1997-2007年由GOES8, GOES11和GOES12星载高能粒子探测器在地球同步轨道高度上所探测到的高能质子和高能电子通量探测数据以及高度560km左右星载大气密度探测器所得的热层大气密度探测数据, 统计分析了强地磁扰动、高能粒子通量跃变和热层大气密度涨落之间的相关关系, 初步获得强地磁扰动期间, 地球同步轨道(外辐射带外环)均出现了增幅大于三个数量级的高能质子通量(尤其是E>1MeV)强增强现象, 随后热 层大气密度强烈上涨, 表明三者之间是正相关关系. 在时间上地球同步轨道高能质子通量强增强现象先于日均Ap值(地磁活动程度)上涨约一天左右, 而热层大气密度强涨落现象又明显滞后于强地磁扰动事件.      

"神舟二号"大气密度探测器的探测结果Ⅱ.在太阳和地磁扰动期间高层大气密度的变化

选用了Sz—2大气密度探测器在2001年2—4月间的探测数据,进行日照和阴影区域热层大气密度对太阳和地磁活动程度的响应变化的探讨．结果表明,日照区大气密度峰值主要随F10.7值而变,在地磁扰动期间,阴影区大气密度对扰动的响应更明显,通常响应变化开始于高纬度地区,然后向低纬度地区推移．     

"神舟3号"运行高度上大气密度的变化

"神舟3号"(SZ-3)大气密度探测器搭载在SZ-3留轨舱上于2002年3月发射入轨,在轨运行期间获得了轨道舱运行高度范围(330-410km)内的大气密度数据.数据分析表明,无明显太阳和地磁扰动时,热层大气密度的主要变化之一是日照和阴影区域之间的涨落变化,最大涨落变化比约为3.0,变化比与太阳和地磁活动程度有关.在2002-04-17和2002-04-19的强地磁扰动时,全球热层大气密度上涨,同时在磁扰峰期探测获得30°N-40°N区域出现密度扰动异常现象.对强地磁扰动在运行轨道高度上大气密度最大涨幅约为60%左右,响应过程在时间上要比地磁扰动过程滞后6-7h,日照和阴影区域中大气密度的响应变化程度明显不同.在太阳活动程度发生变化时,热层大气密度会呈现出明显的正相关变化关系.     

一种NRLMSISE-00模型太阳辐射校正方法

根据空间天气的状态,调整大气模型的相关输入参数能够减小模型的计算误差.通过对比CHAMP卫星在轨大气密度探测数据与NRLMSISE-00模式的计算结果发现,通过调整F_10.7的输入,使轨道大气密度积分的模式计算结果与探测结果之间的误差达到最小,此时的F_10.7被称为理想F_10.7输入（F*）.进一步的分析发现,F*与太阳紫外辐射MgII指数存在很好的相关性,因此可以选择其他的太阳紫外辐射代理参数取代F_10.7,从而减小模型计算误差.本文采用神经网络技术,建立新的太阳紫外辐射代理参量F_euv与MgII,F_10.7等的对应模型,能够根据当日参数值计算F_euv.研究结果表明,新的代理参数能够有效减小NRLMSISE-00的计算误差.      

神舟4号大气成分探测的新结果

神舟4号(SZ-4)大气成分探测器搭载在SZ-4留轨舱上于2002年12月30日发射入轨，在在轨运行的前3个多月中，正值地球南半球处于夏半球季节，并发生了多次中低强度的地磁扰动事件，SZ-4大气成分探测器测得了轨道舱运行高度上(330-362km附近)大气成分的响应变化和异常现象新结果．探测数据表明，中低强度的地磁扰动事件也能引起热层大气中主要成分O和N2的数密度值增高的响应变化．同样在进入地磁扰动峰期后较高纬度处出现了N2的异常增变和O的异常降变以及大气密度异常扰动的现象，但此期间所出现异常现象的地域与SZ-3和SZ-2大气成分探测器探测结果相反，它仅位于南半球较高纬度地区．     

基于活动区面积统计的F10.7指数预报方法

太阳10.7cm射电流量（F_10.7）是反映太阳整体活动的重要指标,其主要源头是日面活动区.F_10.7指数与日面活动区具有显著的相关性,且不同面积的活动区与F_10.7并不遵循相同的线性关系.为进一步提高F_10.7预报的准确性,利用日面活动区面积与F_10.7的相关性,依据面积大小分类,提出F_10.7的预报公式并进行验证.采用2008-2018年SWPC （Space Weather Prediction Center）公布的活动区面积数据和CSWFC （Canadan Space Weather Forecast Center）公布的F_10.7实测数据计算预报公式系数,利用高年（2003年）和低年（1997年）的F_10.7预报验证其结果.研究结果表明,预报结果与实测值的相关系数分别为0.9318和0.9295,二者皆优于SWPC同时期的预报结果（相关系数分别为0.9186和0.8771）.本研究首次基于活动区的变化预测了F_10.7,提高了F_10.7预测的准确性.      

Consequences of the application of the streamer fluid model to the study of the sprite inception mechanism

We present the results of a streamer-fluid model used to investigate the electrodynamical coupling between the troposphere and upper atmosphere due to the penetration of lightning electric fields into the mesosphere and the lower ionosphere, generating sprites. The model solves the continuity equation for electrons and ions coupled to Poisson equation. The dominant physical response of the atmosphere is the formation of a screening-ionization wave. The wave shields the atmosphere above it from the action of the lightning field and, together with the conductivity reduction below it due to attachment, the wave amplifies the total field below it, allowing for the penetration of intense electric fields in the mesosphere as it propagates downwards into regions of higher density that compress the wave. This is the key physical mechanism for sprite inception. We evaluated the effects of the thundercloud charge geometry, lightning current waveshape, atmospheric conductivity, via different electron density profiles, and the effect of ionization, attachment and electron mobility coefficients in the electrical breakdown process, related to halo production, and sprite streamer initiation. The results showed that electrons with higher mobility are more efficient in shielding the lightning electric field before breakdown, causing delay, and they contribute to the formation of the streamer seed after breakdown, anticipating the sprite streamer inception. Similarly, a higher effective ionization rate, produced by modifications in the attachment and ionization coefficients, anticipates sprite inception. The simulations with 6 different electron density profiles, and therefore conductivities, spanning 4 orders of magnitude, showed that the altitude of breakdown and sprite initiation, as well as their time delays from the lightning discharge are directly related to atmospheric conductivity: higher conductivities produce halo and sprite inception at lower altitudes with longer delays and may hinder sprite formation. We document that variations of 30 times in the lightning current leads to sprite initiation altitudes in the range 66.0–73.5 km, with delays between 1.550 and 34.500 ms, while variations of 4 orders of magnitude in the conductivity profile lead to initiation altitudes 61.0–70.6 km, with delays in the range 3.825–9.825 ms. Consequently, we suggest that lightning characteristics dominate over atmospheric parameters in determining sprites’ initiation altitude and delay. The simulation of a −CG, with a constant current of 30 kA, did not produce a sprite seed, confirming an asymmetry in the response of the atmosphere to positive and negative lightning. This is due to the free electron drift direction that is away from the screening ionization wave, preventing the formation of the streamer seed for the great majority of −CGs. The same does not apply to halos, which depend on the occurrence of breakdown and can be produced by discharges of both polarities.     

落塔微重力水平激光干涉测量方法中的干扰因素分析

落塔是获得微重力环境的重要设施,落塔微重力水平的测量对微重力科学实验的研究至关重要. 激光干涉是测量落塔微重力水平的一种新方法,这种方法的基本原理是让一个参考落体在落舱中自由下落,落舱由于受到服外空气阻力的作用将与自由落体运动略有差异,1 11. 肫内的参考落体则更接近理想的自由落体运动,这就使得落舱与参考落体之间存在着加速度差,这种加速度差便反映了落服的微重力水平,其所导致的相对运动则可通过激光干涉的方法测量出来. 本文对落塔微重力水平的激光干涉测量方法中将会遇到的一些主要干扰因素进行了分析,计算结果表明,这些干扰网素所造成的总误差约为 1.2 x 10-7g,低于微重力水平的预测值 10-4~lO-6 g,因此该方法是一种比较可行的测量方法.      

Possible physical properties of the flow in the vicinity of the heliopause

An interface between the fully ionized hydrogen plasma of the solar wind (SW) and the partially ionized hydrogen gas flow of the local interstellar medium (LISM) is formed as a region where there is a strong interaction between these two flows. The interface is bounded by the solar wind termination shock (TS) and the LISM bow shock (BS) and is separated on two regions by the heliopause (HP) separating the solar wind and charged component of the LISM (plasma component below). The BS is formed due to the deceleration of the supersonic LISM flow relative to the solar system. Regions of the interface between the TS and HP and between the HP and BS were in literature named as the inner and outer heliosheaths, respectively. An investigation of the structure and physical properties of the heliosheath is at present especially interested due to the fact that Voyager-1 and Voyager-2 have crossed the TS in December 2004 (Burlaga, L.F., Ness, N.F., Acuna, M.Y., et al. Crossing the termination shock into the the heliosheath. Magnetic fields. Science 309, 2027–2029, 2005; Fisk, L.A. Journey into the unknown beyond. Science 309, 2016–2017, 2005; Decker, R.B., Krimigis, S.M., Roelof, E.C., et al. Voyager 1 in the foreshock, termination shock and heliosheath. Science 309, 2020–2024, 2005; Stone, E.C., Cummings, A.C., McDonald, F.B., et al. Voyager 1 explores the termination shock region and the heliosheath beyond. Science 309, 2017–2020, 2005) and in September 2007 (Jokipii, J.R. A shock for Voyager 2. Nature 454, 38–39, 2008; Gurnett, D.A., Kurth, W.S. Intense plasma waves at and near the solar wind termination shock. Nature 454, 78–80, 2008. doi: 10.1038/nature07023; Wang, L., Lin, R.P., Larson, D.E., Luhmann, J.G. Domination of heliosheath pressure by shock-accelerated pickup ions from observations of neutral atoms. Nature 454, 81–83, 2008. doi: 10.1038/nature07068.14; Burlaga, L.F., Ness, N.F., Acuna, M.H., et al. Magnetic fields at the solar wind termination shock. Nature 454, 75–77, 2008. doi: 10.1038/nature07029; Richardson, J.D., Kasper, J.C., Wang, C., et al. Cool heliosheath plasma and deceleration of the upstream solar wind at the termination shock. Nature 454, 63–66, 2008. doi: 10.1038/nature07024; Stone, E.C., Cummings, A.C., McDonald, F.B., et al. An asymmetric solar wind termination shock. Nature 454, 71–74, 2008. doi: 10.1038/nature07022; Decker, R.B., Krimigis, S.M., Roelof, E.C., et al. Mediation of the solar wind termination shock by non-thermal ions. Nature 454, 67–70, 2008. doi: 10.1038/nature 07030), respectively, and entered to the inner heliosheath.     

Testing the improvement of performance of the IRI model in the estimation of TEC over the mid-latitude American regions

This paper mainly discusses the improvement of performance of the International Reference Ionosphere (IRI) model in estimating the variation of the Vertical Total Electron Content (VTEC) over the mid latitude American regions during the relatively low (2008–2010) and relatively high (2012) solar activity years. This has been conducted employing the VTEC values obtained from the dual frequency ground based Global Positioning System (GPS) receivers located at Mineral Area Community College, MACC (37.85°N, 269.52°W) and Mississippi County Airport, MAIR (36.85°N, 270.64°W), and the latest versions of the IRI online model (IRI 2007, IRI 2012 and IRI 2016). The study mainly focuses to compare the trend of variability of the monthly and seasonal modeled VTEC values (IRI 2007 VTEC, IRI 2012 VTEC and IRI 2016 VTEC) with the corresponding measured VTEC values (GPS VTEC). The overall results show that the IRI VTEC values (almost in all versions of the model) are generally smaller than the GPS VTEC except after about 15:00 UT (09:00 LT) in the December solstice when the Sun shifts to the high solar activity. On the contrary, overestimations of the VTEC values by the model are observed in traversing from the low solar activity (2008) to high solar activity (2012) phase, especially after about 15:00 UT (09::00 LT) with the IRI 2016 version showing the highest. In general, the IRI 2007 and IRI 2012 versions show similar monthly and seasonal underestimations or overestimations showing that the two versions have almost similar performance. The IRI 2016 version is generally better in capturing both the diurnal and arithmetic mean GPS VTEC values with some exceptional months and seasons as compared to those of the IRI 2007 and IRI 2012 versions.     

Features of the occurrence of the additional stratification on the bottom-side F region over the equatorial location of Trivandrum

The general features of occurrence of an additional layer on the bottom side of F region, referred to as F0.5 layer in the pre noon period, over the magnetic equatorial location of Trivandrum (8.5° N; 77° E; dip lat of 0.5° N) in India during the period from 2004 to 2007 are presented using ionosonde observations. The F0.5 layer has a June (northern summer) solsticial maximum probability of occurrence with secondary maxima during December (northern winter) solstice. The seasonal as well as the day-to-day variability in the occurrence of F0.5 layer as mentioned in this paper seems to be a result of the variations in the amplitude and phases of the tides and gravity waves, and inventory of the metallic ions of meteoric origin. This study brings out an important manifestation of morning time F layer base region dynamics.     

Method for the prediction of the effective dose equivalent to the crew of the International Space Station

This paper describes a methodology for assessing the pre-mission exposure of space crew aboard the International Space Station (ISS) in terms of an effective dose equivalent. In this approach, the PHITS Monte Carlo code was used to assess the particle transport of galactic cosmic radiation (GCR) and trapped radiation for solar maximum and minimum conditions through an aluminum shield thickness. From these predicted spectra, and using fluence-to-dose conversion factors, a scaling ratio of the effective dose equivalent rate to the ICRU ambient dose equivalent rate at a 10 mm depth was determined. Only contributions from secondary neutrons, protons, and alpha particles were considered in this analysis.     

电离层电场对磁层活动响应时延的模式研究

本文用极盖边界上电离层驱动电位φ₀随时间变化的不同模式计算了场向电流J_∥2及电离层对流电场E的演化过程。计算表明,当φ₀(t)先升后降有极大值时,J_∥2(t)和E(t)也表现出类似的趋势。但它们的极值滞后于φ_0max出现的时刻,即在一段时间内,φ₀虽已开始下降,J_∥2及高纬E却继续增大。一般说,E先于J_∥2达到极值,但相差甚小。φ₀变化形式不同时,滞后时间亦不同。当φ₀陡升缓降时,E、J_∥2的极值相对于φ_0max的时延可超过一小时。这与持续时间较长的磁暴期间所观测到的电离层场强响应时延量级是一致的。时延大小还受电离层电导率的制约。低纬电离层场强的响应与高纬不同,其升降趋势与φ₀同步。      

Influence of the Ionosphere on the Observation of the Mid-Latitude Pi2 Pulsations at the Global Scale

Results of our investigation showed that occurrence frequency of Pi2 over a 24 hour period undergoes seasonal variations in time coincidence with foF2. In the winter months, at sunrise and sunset (when foF2gradients are the largest) the observation probability of these oscillations is minimal. At periods of summer solstice when the F2-layer persists almost round the clock, no effect of Pi2 pulsation attenuation is observed at sunrise and sunset. The pulsation amplitudes behave in a similar manner. Results of this study suggest the conclusion that the propagation of signal from the Pi2 sourse into the mid-latitudes, and also the parameters of these pulsations are essentially affected by electron density in the ionospheric F2-layer.      

Effects of the position of the solar wind termination shock and the heliopause on the heliospheric modulation of cosmic rays

The effects of changing the position of the solar wind termination shock and the position of the heliopause, and therefore the extent of the heliosheath, on the modulation of cosmic ray protons are illustrated. An improved numerical model with diffusive termination shock acceleration, a heliosheath and drifts is used. The modulation is computed in the equatorial plane and at 35 heliolatitude using recently derived diffusion coefficients applicable to a number of cosmic ray species during both magnetic polarity cycles of the Sun. It was found that qualitatively the modulation results for the different heliopause positions are similar although they differ quantitatively, e.g., clearly different radial gradients are predicted for the regions beyond the termination shock compared to inside the shock. The difference between the modulation for the two solar polarity cycles are less significant at a heliolatitude of 35° than in the equatorial plane. We found that moving the termination shock from 90 to 100 AU, with the heliopause fixed at 120 AU, caused only quantitative differences so that the exact position of the TS in the outer heliosphere seems not crucially important to global modulation. Moving the heliopause outwards, to represent the modulation in the tail region of the heliosphere, causes overall decreases in the cosmic ray intensities but not linearly as a function of energy, e.g., at 1 GeV the effect is insignificant. We conclude from this modelling that the modulation of protons in the heliospheric nose and tail regions are qualitatively similar although, clear quantitative and interesting differences occur.     

Comprehensive investigation of the basic parameters of the upper atmosphere at the time of the flight of the geophysical rocket “vertical-6”

Ion temperature and total ion concentration measured on 25th October 1977 during the flight of the geophysical rocket “Vertical-6” are analyzed. The solar EUV fluxes determined in five wave-length bands with a photoelectron analyzer are also given. The observed anomalous variation of ion temperature between 700 and 900 km and the measured ion concentration can be explained, if the charge exchange reactions H⁺ ? O⁺ and diffusion are taken into account.