不同强度磁暴期间热层大气密度响应特征分析

利用CHAMP卫星数据，对2002-2008年12个不同强度磁暴事件期间的热层大气密度变化特征进行分析，并研究对应磁暴期间大气模式NRLMSISE-00分布特征.结果表明，大磁暴期间日侧大气密度峰值从高纬到低纬的时间延迟为2h，中小磁暴期间的延迟时间为3～4h；春秋季暴时大气密度分布基本呈南北对称分布，而夏冬季大气密度的分布是夏半球大于冬半球，春秋季暴时大气密度大于夏冬季；NRLMSISE-00大气模式得到的热层大气密度很好的体现了半球分布以及季节分布的特征，但模式模拟结果偏小；Dst指数峰值比ap指数峰值更能反应大气密度的变化情况.      

地磁暴前地面宇宙线强度变化Morlet小波特征分析

用Morlet小波变换对Oulu台站和Apatity台站(主要是Apatity台站)1998-2002年间宇宙线静日和地磁暴前的地面宇宙线强度变化特征进行分析，得到：在宇宙线静日期间普遍存在准24h周期变化特征，并且在当地时间0200，1400左右分别出现最小值和最大值；对12个例子的分析可以看到地面宇宙线强度在地磁暴之前l-2天均出现了不同于宇宙线静日期间变化特征，或有小尺度周期出现，或周期变化完全消失，或有异常24h周期变化，这种变化特征在一定程度上可作为地磁预报的先兆特征之一。     

神经网络方法在太阳质子事件短期预报中的应用

在大量统计结果的基础上,深入研究了太阳质子事件预报机理．总结了质子事件爆发与太阳活动区面积、位置、McIntosh结构、磁结构以及前两天活动区爆发耀斑事件数目之间的关系．然后,在神经网络的基础上建立了太阳质子事件短期预报模型,并对2000年以后12个未参加训练的样本进行测试,结果对事件预报的准确率为83％．此外,我们还利用该模型对2002年1-4月发生的几次质子事件进行了预报试验,结果发现,这期间发生的6次事件都被预报．其中3次质子事件系统预报提前了3天,两次事件预报提前了2天,一次事件提前1天预报．     

太阳10.7 cm射电辐射流量预报方法初探

引进“相似周”方法,通过两种方式对第23周太阳10.7cm辐射流量（F10.7）月均值进行预报和预报,由“相似周”方法得到的第23周太阳黑子数月均平滑预测值来预测F10.7月均值和直接采用“相似周”方法对F10.7月均平滑值进行预测,通过对两种预报试验结果分析,得到以下结论。（1）两种预测结果与实际都比较吻合,都得到双峰结构。（2）直接通过相似周预测的F10.7的月均值结果较间接通过太阳黑子数的预测结果所推断的F10.7预测结果更接近实际观测结果。（3）使用“相似周”,预报方法,可以给出F10.7比较精细的剖面结构,这是其他普通预报方法很难做到的。     

基于实时观测数据的大气密度模式修正

针对国际大气密度模式NRLMSISE-00, 以中国神舟飞船探测数据为基础, 提出一种基于实时大气密度观测数据的模式修正方法. 通过计算分析模式计算结果与探测数据的误差分布特征, 针对地磁相对平静期(Ap≤ 30)模式计算的误差特点, 建立了一种平均误差修正方法, 即认为在相对平静期, 在相同纬度和地方时, 模式误差基本相同, 某一时刻模式预测误差可以近似用与其相同纬度和地方时的平均误差来替代, 从而对模式预测结果进行修正. 以神舟4号探测数据为基础, 通过对模式预测结果采用两种方式进行修正, 可以看到模式误差得到了一定的改善. 采用误差库累积准实时修正, 修正后的误差由原来的20 %降至6 %; 采用误差库5天滑动预报修正后, 模式提前1, 2, 3天的预测误差由原来的20 %分别降至7.8 %, 9.4 %和10.5%.      

利用相似周方法预测第24活动周开始时间

利用相似周方法对第24活动周的开始时间与第23活动周下降相后期的太阳黑子数进行了预报.根据第23周已经出现的特征参量和下降相的形态特征,选取9,10,11,15,17和20等六个太阳活动周作为第23周下降相的相似周,对第24周开始时间进行预报.预报结果显示,第24活动周的开始时间为2007年5±1月,黑子数平滑月均极小值为7.1±2.6,第23太阳活动周长度为11.1年.与其他研究者的预报结果相比较,本文给出的结果与文献[11]和[12]及MSFC的结果比较一致.通过对相似周方法在下降相预测太阳活动周结束时间的研究讨论,及对第23周上升阶段的太阳黑子数和F10.7平滑月均值预报结果的评估,可以看出,相似周预报方法在太阳活动周长期预报中是很有应用价值的.      

大气密度模型球谐修正研究

大气模型修正是提高模型精度的一种重要方法.利用CHAMP卫星高精度加速仪反演的密度数据,采用球谐函数的形式对NRLMSISE-00模型进行修正.为了消除轨道高度变化对密度修正结果的影响,将密度数据同化到同一高度处,计算修正之后的密度误差,进而对未来三天的密度进行预报.结果表明,经球谐修正后,修正误差和预报误差均有显著降低.在太阳活动高年,修正误差可降至10%左右,提前1~3天预报精度分别提高31.34%,21.39%和13.75%;太阳低年时修正误差可降至14%左右,提前1~3天预报精度分别提高55.03%,47.79%和43.60%.      

太阳黑子数及Ap指数周期变化特征的小波分析

应用Morlet小波变换方法从多个变化尺度上对1932—2000年的太阳黑子数及Ap指数的变化特征进行分析．(1)太阳黑子数存在准11年、准32年的周期变化特征及Ap指数存在准32年、准11年、准6个月、准27天和准13．9天的周期变化特征；(2)太阳黑子数及Ap指数有着相似的准11年周期变化，但Ap指数极值的出现要比太阳黑子数极值出现滞后1—2年；(3)Ap指数准27天的周期变化在太阳黑子活动高、低年不同，在太阳活动低年，Ap指数有着较稳定的准27天周期变化，但在太阳活动高年，27天的周期变化几乎消失，这种周期变化的消失和出现时间可在Morlet小波变换图中体现出来。     

中长期轨道预报中大气阻力系数补偿算法的研究

利用地磁平静期的CHAMP卫星精密星历对多个弧段的大气阻力系数进行解算,找到了适用于轨道预报的最优系数,分析了最优阻力系数与地磁指数、迹向残差的关系,利用线性回归分析建立了大气阻力系数补偿算法。针对中长期轨道预报的需求,该算法能够降低定轨弧段较短条件下阻力系数的解算误差。将该补偿算法应用于不同时期CHAMP卫星和天宫一号的轨道预报,验证了该算法的正确性和普遍适用性,结果表明预报精度能够提高20%以上。     

太阳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%。