太阳E10.7指数的反演及预报方法

针对高层大气密度预报和轨道预报业务中对新型太阳紫外辐射指数E_10.7的需求,基于TIMED-SEE观测仪器提供的0.1～105 nm太阳辐射强度数据,开展了E_10.7指数反演和中期预报研究。 E_10.7指数是太阳光谱中波长为0.1～105 nm的辐射流量,单位与F_10.7指数相同（sfu,1 sfu=10–22 W·m–2·Hz–1）。 TIMED-SEE观测仪器提供的0.1～105 nm太阳辐射强度实测值具有高时间分辨率、延迟时间短和易获得的优势,利用最小二乘法拟合可反演出准实时的E_10.7指数,均方根误差为5.445 sfu。利用高阶自回归模型对E_10.7的中期预报效果尚佳,未来27天的预报值平均相对误差为7.83%。利用同样方法还开展了E_10.7指数81天中心滑动平均值未来27天预报试验,未来27天的预报值平均相对误差仅为3.63%。      

第23太阳活动周太阳极紫外辐射与磁暴相关性分析

统计第23个太阳活动周内中等及以上强度（Dst_min<-50nT）的磁暴事件,线性拟合分析磁暴主相DDst_min和达到DDst_min前一个表征太阳极紫外辐射强度的F_10.7之间的相关性.结果表明:随着太阳极紫外辐射增强,DDst_min<-50nT的磁暴出现的总数增多,在弱、中等和强太阳极紫外辐射条件下,其数量分别为56,84和85;随着太阳极紫外辐射增强,强磁暴（-200nT ≤ Dst_min<-100nT）和大磁暴（Dst_min<-200nT）发生的数量和相对发生率呈增长趋势,尤其是大磁暴数目（1,4,12）和相对发生率（1.79%,4.76%,14.12%）明显呈增长趋势;大磁暴（|Dst_min|）与太阳极紫外辐射（F_10.7）之间存在中度正相关关系,其相关系数为0.532,并且主要体现在大磁暴（|Dst_min|）与强太阳极紫外辐射（F_10.7）之间的中度正相关性,其相关系数为0.582.大磁暴与强太阳极紫外辐射之间的相关性可为空间天气预报提供参考依据.      

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%以内.      

一种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的计算误差.      

第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天变化值之间的相关系数.      

NRLMSISE-00大气模型与GRACE和CHAMP卫星大气密度数据的对比分析

利用GRACE(Gravity Recovery And Climate Experiment)和CHAMP(Challenging Mini-Satellite Payload)卫星2002-2008年的大气密度数据与NRLMSISE-00大气模型密度结果进行比较,分析了模型密度误差及其特点.结果显示,NRLMSISE-00大气模型计算的密度值普遍偏大,其相对误差随经纬度变化,在高纬度相对较小;相对误差随地方时变化,在02:00LT和15:00LT左右较大,10:00LT和20:00LT左右较小.通过模型密度相对误差与太阳F_10.7指数的对比分析发现,在太阳活动低年模型相对误差最大,而在太阳活动高年相对误差较小;将模型结果分别与GRACEA/B双星和CHAMP卫星的密度数据进行比较,发现对于轨道高度更高的GRACE卫星轨道,模型相对误差更大;在地磁平静期,相对误差与地磁ap指数(当前3h)相关性不强,但是在大磁暴发生时,误差急剧增大.      

基于活动区面积统计的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预测的准确性.      

基于SAMPEX卫星观测的南大西洋异常区高能质子动态分布特征

利用SAMPEX卫星1992年7月至2004年6月19～27MeV高能质子数据对南大西洋异常区的分布特征进行研究, 发现南大西洋异常区高能质子分布随高度及F_10.7的变化十分显著. 在540±25km高度上, 地磁较为平静时期南大西洋异常区高能质子微分通量随着F_10.7的增大而减小, 同时在F_10.7≥115sfu时减小趋势较为平缓. 对中等及以上磁暴进行统计分析发现, 磁暴期间南大西洋异常区高能质子微分通量和SYM-H指数的绝对值存在明显的反相关关系, 且地磁暴对南大西洋异常区高能质子微分通量存在明显的持续影响效应. 磁暴发生期间高能质子微分通量明显减少. 磁暴恢复相及其之后高能质子微分通量呈现较为显著的恢复过程.      

2007-2009谷年期间高热层大气密度的变化

在第23至第24太阳活动周的峰年之间,太阳活动谷年具有持续时间长,极低F_10.7太阳辐射通量（低至65）和超长期的零太阳黑子数记录等特点,因此是观测和研究在这种特殊背景下热层大气变化的极好机会.尤其是能充分理解和掌握在宁静环境下热层大气密度对弱太阳活动和小地磁扰动的响应特性.本文利用高度650 km以上星载大气密度探测器2007—2009年的连续探测数据进行分析,结果表明,在太阳辐射通量F_10.7极低值期间,较高热层大气密度对F_10.7的起伏具有更显著的响应变化.当F_10.7由70降至65时,日均大气密度会有4～5倍的显著降变,远大于通常大气模式中的降变值.同样在F_10.7极低值期间,较高热层大气密度对小地磁扰动也具有显著的响应增变,当日K_p指数之和由23增至30时,较高热层大气密度则会有80%～160%的强增变.     

太阳活动对电离层TEC变化影响分析ormalsize

为研究太阳活动对电离层TEC变化的影响,从整体到局部分析了2000—2016年的太阳黑子数、太阳射电流量F_10.7指数日均值与电离层TEC的关系,并重点分析了2017年9月6日太阳爆发X9.3级特大耀斑前后15天太阳活动与电离层TEC变化的相关性.结果表明:由2000—2016年的数据整体看来,太阳黑子数、太阳F_10.7指数、TEC两两之间具有很强的整体相关性,但局部相关性强弱不均;此次耀斑爆发前后太阳黑子数、太阳F_10.7指数和TEC具有很强的正相关特性,太阳活动对TEC的影响时延约为2天;太阳活动对全球电离层TEC的影响不同步,从高纬至低纬约有1天的延迟,且对低纬度的影响远大于中高纬度.太阳活动是影响电离层TEC变化的主要原因,但局部也可能存在其他重要影响因素.      

High accuracy satellite drag model (HASDM)

The dominant error source in force models used to predict low-perigee satellite trajectories is atmospheric drag. Errors in operational thermospheric density models cause significant errors in predicted satellite positions, since these models do not account for dynamic changes in atmospheric drag for orbit predictions. The Air Force Space Battlelab’s High Accuracy Satellite Drag Model (HASDM) estimates and predicts (out three days) a dynamically varying global density field. HASDM includes the Dynamic Calibration Atmosphere (DCA) algorithm that solves for the phases and amplitudes of the diurnal and semidiurnal variations of thermospheric density near real-time from the observed drag effects on a set of Low Earth Orbit (LEO) calibration satellites. The density correction is expressed as a function of latitude, local solar time and altitude. In HASDM, a time series prediction filter relates the extreme ultraviolet (EUV) energy index E_10.7 and the geomagnetic storm index a_p, to the DCA density correction parameters. The E_10.7 index is generated by the SOLAR2000 model, the first full spectrum model of solar irradiance. The estimated and predicted density fields will be used operationally to significantly improve the accuracy of predicted trajectories for all low-perigee satellites.     

太阳极紫外辐射对热层大气密度的影响

对2001-2021年SOHO卫星的极紫外辐射测量数据,以及CHAMP,GRACE-A和SWARM-C卫星资料推导出的高分辨率大气密度数据进行统计分析,发现大气密度与极紫外测量值的相关系数大于密度与F10.7指数的相关系数,证实极紫外辐射在不同地方时的影响程度存在显著差异,从而驱动大气密度的周日变化。利用三颗卫星的高度差异揭示极紫外辐射对大气密度的加热效应在350～500 km范围随着高度增加而减弱。统计得到极紫外辐射影响在地方时和纬度上的空间差异：对夏季半球的影响大于冬季半球;在白天,对中纬度地区的影响高于赤道和高纬度地区;在夜间,密度对辐射的斜率在夏季半球高纬度地区存在峰值,在冬季半球中纬度存在谷值,模型DTM2000和NRLMSISE00未能准确刻画。为了改进经验模型,提出基于球谐函数的拟合方法,优于主流模型周日效应采用的表达式,对周日效应建模和修正提供有益借鉴。利用昼夜间能量传输和热层大气经向环流机制探讨了统计结果的物理机制。     

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

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

E10.7指数在软X射线谱段的修正应用

表征EUV辐射通量的E10.7指数在越来越多的研究和应用中被用来代替传统的F10.7指数.X射线对地球D层和E层的电离起着重要作用,但由于D层观测数据的不足和E层电离源的多样性,难以被用来考虑X射线对电离层的影响.火星电离层下层的电离源几乎是单一的软X射线,这为研究X射线对电离层的作用提供了可能性.通过研究火星电离层下层的峰值电子浓度对E10.7的依赖关系,发现即便经过必要的修正,这种关系对不同的观测时段并不具备一致性.通过理论推导和数据分析,得到了一种特别用于描述太阳软X射线辐射通量的新指数,即Xs指数,用来替代E10.7指数.Xs指数在描述火星电离层下层对太阳辐射的依赖关系时,不同的观测时段有很好的一致性,表明Xs指数在表征太阳软X射线辐射强度方面比E10.7指数更加合适.      

Periodic Variations of Drag Coefficient for the ANDE Spherical Satellites During its Lifetime

A drag coefficient (C_D) inversion method is introduced to study the variations of the drag coefficient for orbital satellites with spherical geometry. Drag coefficients of the four micro satellites in the Atmospheric Neutral Density Experiment (ANDE) are compiled out with this new method. The Lomb-Scargle Periodgram (LSP) analysis of the four ANDE satellites' C_D series has shown that there are obvious 5, 7, 9, and 27 days' period in those data. Interesting results are found through comparing the LSP analysis with series of the daily solar radio flux at 10.7 cm (F_10.7 index), the Ap index, and the daily averaged solar wind speed at 1AU. All series in the same time interval have an obvious period of about 27 days, which has already been explained as the association with the 27 days' solar rotation. The oscillating periods less than 27 days are found in series of C_D, Ap and solar wind speed at 1AU, e.g., the 5, 7, 9 days period. However, these short periods disappeared in the time series of F_10.7 index. The same periodicities of 5, 7, 9 days in Ap and solar wind are presented at the same time interval during the declining phase of solar cycle 23. While in the ascending phase of solar cycle 24, these short oscillations are not so obvious as that in the declining phase of solar cycle 23. These results provide definite evidence that the C_D variations with period of 5, 7 and 9 days are produced by a combination of space weather effects caused by the solar wind and geomagnetic activity.      

Solar EUV and decimetric indices and thermospheric models

Data bases and limits of applicability of existing empirical thermospheric models are reviewed by using these models together with solar EUV irradiance data in studying the solar activity effect on composition, density and temperature. For two rather short aeronomy missions of the AEROS A and B satellites solar EUV indices as proposed by Schmidtke are used in comparison with the 10.7 cm solar flux F in determining the solar activity effect in in-situ composition measurements sampled by the same satellites at 250, 310 and 380 km altitude. No advantage of solar EUV indices over F could be determined.