利用OH夜气辉反演原子氧数密度时输入参数的不确定度对反演结果的影响

借助OH夜气辉辐射的光化学模式,由OH夜气辉辐射反演中间层-低热层区域的原子氧数密度时,输入参数的不确定性将导致反演得到的原子氧数密度具有不确定性.以在sudden death猝灭模式下通过OH(8-3)振动带体辐射率反演原子氧数密度为例,分别研究了大气参数和OH气辉辐射率的不确定度引起的反演不确定度、化学反应速率常数的不确定度引起的反演不确定度,以及所有输入参数的不确定度共同引起的反演不确定度,找出其不确定度对反演结果影响最大的参数.结果表明,三种反演不确定度均随着高度的升高而增大,温度和体辐射率的不确定度对第一种反演不确定度的贡献最大,反应速率常数b(8)和A(8-3)、的不确定度对第二种反演不确定度的贡献最大.     

冕洞相关地磁Ap指数中短期预报方法研究

地磁Ap指数是描述全球地磁活动水平的重要指数, 过去许多参考大气模式中都用Ap指数来表述地磁活动状态, 大气模式的运行需要输入地磁Ap指数, 因此, 地磁Ap指数的预报一直是空间环境预报中一个非常重要的内容. 针对太阳活动低年冕洞引起的地磁扰动具有明显27天重现的特性, 利用修正的自回归方法, 对地磁Ap指数进行了提前27天的预报; 采用从SOHO/EIT观测资料发展出来的描述冕洞特性的Pch因子, 进行了提前三天的地磁Ap指数预报. 结果显示, 将统计方法与物理分析相结合, 进行地磁Ap指数的中短期数值预报, 可以得到较好的预报效果.      

Advanced ensemble modeling method for space object state prediction accounting for uncertainty in atmospheric density

For objects in the low Earth orbit region, uncertainty in atmospheric density estimation is an important source of orbit prediction error, which is critical for space traffic management activities such as the satellite conjunction analysis. This paper investigates the evolution of orbit error distribution in the presence of atmospheric density uncertainties, which are modeled using probabilistic machine learning techniques. The recently proposed “HASDM-ML,” “CHAMP-ML,” and “MSIS-UQ” machine learning models for density estimation (Licata and Mehta, 2022b; Licata et al., 2022b) are used in this work. The investigation is convoluted because of the spatial and temporal correlation of the atmospheric density values. We develop several Monte Carlo methods, each capturing a different spatiotemporal density correlation, to study the effects of density uncertainty on orbit uncertainty propagation. However, Monte Carlo analysis is computationally expensive, so a faster method based on the Kalman filtering technique for orbit uncertainty propagation is also explored. It is difficult to translate the uncertainty in atmospheric density to the uncertainty in orbital states under a standard extended Kalman filter or unscented Kalman filter framework. This work uses the so-called “consider covariance sigma point (CCSP)” filter that can account for the density uncertainties during orbit propagation. As a test-bed for validation purposes, a comparison between CCSP and Monte Carlo methods of orbit uncertainty propagation is carried out. Finally, using the HASDM-ML, CHAMP-ML, and MSIS-UQ density models, we propose an ensemble approach for orbit uncertainty quantification for four different space weather conditions.     

Adaptive sliding mode control of spacecraft attitude-orbit dynamics on SE(3)

In this paper, to solve the problem of parameters uncertainty in spacecraft tracking control, an adaptive controller based on sliding mode is proposed for the relative spacecraft attitude-orbit dynamics on the Lie group SE(3). The dynamic equations of relative attitude orbit error for two spacecraft are established in the framework of Lie group SE(3). Considering the uncertainty of spacecraft parameters, a formal decomposition of known and unknown parameters, the state variables and control variables is firstly made in the original system. An online estimator is designed to evaluate the unknown parameters. A sliding mode controller is developed to actuate the spacecraft to track the target spacecraft. Then a Lyapunov function of tracking error and parameters estimated error is designed to prove the stability of the closed-loop system. Finally, the simulation results and analysis are presented to verify the effectiveness and feasibility of the proposed method.     

Orbital error propagation considering atmospheric density uncertainty

Due to the influence of various errors, the orbital uncertainty propagation of artificial celestial objects while orbit prediction is required, especially in some applications such as conjunction analysis. In the orbital error propagation of artificial celestial objects in low Earth orbits (LEOs), atmospheric density uncertainty is one of the important factors that require special attention. In this paper, on the basis of considering the uncertainties of position and velocity, the atmospheric density uncertainty is also taken into account to further investigate the orbital error propagation of artificial celestial objects in LEOs. Artificial intelligence algorithms are introduced, the MC Dropout neural network and the heteroscedastic loss function are used to realize the correction of the empirical atmospheric density model, as well as to provide the quantification of model uncertainty and input uncertainty for the corrected atmospheric densities. It is shown that the neural network we built achieves good results in atmospheric density correction, and the uncertainty quantization obtained from the neural network is also reasonable. Moreover, using the Gaussian mixture model - unscented transform (GMM-UT) method, the atmospheric density uncertainty is taken into account in the orbital uncertainty propagation, by adding a sampled random term to the corrected atmospheric density when calculating atmospheric density. The feasibility of the GMM-UT method considering atmospheric density uncertainty is proved by the further comparison of abundant sampling points and GMM-UT results (with and without considering atmospheric density uncertainty).     

基于卫星加速度数据反演的热层大气密度与NRLMSISE-00模式结果的比较研究

利用CHAMP/STAR加速度数据反演的热层大气密度与NRLMSISE-00模式反演的热层大气密度进行比较, 结果表明, 热层大气密度在春秋季期间高于冬夏季, 并且太阳活动高年比低年更加显著; 日照面和阴影区大气密度的比值在低纬地区由太阳活动高年的4下降到太阳活动低年的2左右, 中纬地区大约由3变化到1.5, 高纬地区变化较小; NRLMSISE-00模式能够较好地模拟热层大气密度的变化趋势, 但是磁暴期间模式精度较差. 统计结果表明, 模式整体比反演结果偏高, 2002-2008年相对偏差分别为16.512%, 20.004%, 18.915%, 18.245%, 25.161%, 33.261%和41.980%; NRLMSISE-00模式在高纬地区的相对偏差为27.337%, 高于中低纬地区的24.047%; 模式在中等太阳活动水平相对偏差较为稳定, 基本在15%左右.      

压力传感器动态特性参数不确定度评定

压力传感器动态特性参数的不确定度是表征其动态测量性能的重要指标。提出了一种压力传感器动态特性参数的不确定度评定方法。首先,使用激波管动态校准系统产生阶跃压力信号激励压力传感器,得到传感器的输出信号;其次,采用基于经验模态分解（EMD）的传感器输出信号预处理方法,减小动态校准过程中噪声的影响;然后,根据传感器的输入输出信号,采用自适应最小二乘法建立压力传感器的数学模型,进而得到其时频域动态特性参数;最后,针对重复校准实验得到的动态特性参数序列的小样本特点,采用自助法计算参数的扩展不确定度和相对不确定度。采用激波管系统对压力传感器进行多次重复动态校准实验,计算时频域动态特性参数的不确定度,并与现有方法进行对比。实验结果表明:本文方法可以弥补贝塞尔法在处理小样本量数据中的不足,且与蒙特卡罗法的不确定度评定结果相对误差小于10%,说明本文方法可以有效地评定压力传感器动态特性参数的不确定度。分析时频域动态特性参数的相对不确定度得到传感器的工作频带和超调量受噪声的影响较大,为动态校准实验条件的改善提供了重要依据。      

Modelling bottom and topside electron density and TEC with profile data from topside ionograms

The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.     

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

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

Thermospheric storm effects

A brief review is given of our current understanding of the atmospheric perturbations in the thermosphere and exosphere that are related to geomagnetic disturbances and of current efforts to represent these in empirical models of the upper atmosphere. A particular model, based on ESRO4 mass spectrometer observations of neutral composition and density, is presented in detail. This model gives the effects on the principal constituents of the upper atmosphere as a function of the geomagnetic coordinates and the K_p geomagnetic index. It is a modification of an earlier model, the most important difference being the inclusion of the variation with magnetic local time.     

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)相关性不强,但是在大磁暴发生时,误差急剧增大.      

太阳活动突发过程对临近空间大气影响的模拟

空间天气对地球及近地空间具有重要影响,大的空间天气事件对中上层大气动力学和成分具有不同的影响。利用全大气耦合模式WACCM,针对太阳耀斑、太阳质子、地磁暴三类事件,以太阳活动平静期2015年5月10-14日的GEOS-5数据为模式背景场,通过F_10.7、离子产生率、Kp及Ap指数设置,分别模拟三类事件对临近空间大气温度、密度和臭氧的影响。结果表明耀斑事件在三类事件中对临近空间大气温度和密度的影响最为显著。平流层大气温度增加是由耀斑辐射增强引起平流层臭氧吸收紫外辐射发生的光化学反应所致,耀斑事件引起平流层和低热层温度增加约为2～3 K,低热层大气相对密度增加在6%以内;太阳质子事件及磁暴事件主要影响低热层,但太阳质子事件和磁暴事件对低热层温度扰动不大于1 K。     

基于IMM-UKF的航天器相对位姿确定算法

在航天器相对导航过程中,相对距离测量信息容易受到干扰,测量误差有较大的不确定性,通常基于单一模型的滤波算法无法对噪声进行辨识,很难获得精确的导航结果。针对应用Clohessy-Wiltshire(C-W)方程受到圆轨道假设的限制问题,研究了建立在惯性坐标系下的近距离相对运动方程(Lawden方程),建立了基于这两个方程的模型集。根据导航系统测量敏感器的特点,设计基于Rodrigues参数及无迹卡尔曼滤波(UKF)的交互式多模型(IMM)视觉相对位姿动态估计算法(IMM-UKF),在保证计算效率的前提下,确保相对轨道姿态确定的稳定性和精确性。数值仿真验证了算法的有效性和先进性。     

多学科设计优化中的不确定性建模

不确定性因素对复杂工程系统设计的影响贯穿整个设计过程.对多学科设计优化(MDO,Multidisciplinary Design Optimization)中不确定性的来源进行了总结和归类;建立了包含2个耦合学科的MDO系统分析输出结果的不确定性与设计变量输入数据误差和模型误差之间关系的数学模型,用于计算各种不确定性对系统性能的影响.算例分析结果和实验结果吻合,表明该数学模型有效可靠.经过扩展可用于估算包含多个耦合学科的MDO系统分析输出结果的不确定性,为鲁棒多学科设计优化方法的构建提供理论基础.      

2005年8月24日强磁暴事件对高层大气密度的扰动

对2005年8月24日发生的突发型强磁暴(Kp峰值达到9)事件,利用星载大气密度探测器在轨实时的连续探测数据进行了处理和分析.结果表明,此次强磁暴事件期间,引起560 km高度附近大气密度剧烈扰动,并存在着两种响应过程.一种是跟随地磁扰动程度变化的全球性大气密度涨落变化,响应时间滞后6h左右, 最大涨落变化比为2.5;另一种为磁暴峰期出现在高纬地区的大气密度突发性跃增,增变比高达5.5.后者存在着区域上的不对称性及时间上的突发性和增幅的差异.此次强磁暴峰期还同时出现了南北半球高纬地区的大气密度跃增双峰.同时还表明这种增变峰可能存在着由高纬向低纬地区迅速推移的现象,在中纬地区推移速度可达15°/h(纬度)左右.      

利用人工神经网络提前1h预报电离层TEC

提出了一种利用人工神经网络提前1h预报电离层TEC的简便方法. 考虑到实际工程应用要求, 没有使用其他空间天气参数, 而是选择电离层TEC观测数据本身作为输入参数. 输入参数为当前时刻TEC、一阶差分、相对差分和时间, 输出参数为预报时刻TEC. 利用文中介绍的GPS/TEC处理方法解算厦门站2004年电离层TEC观测数据, 对预报方法进行评估, 全年平均相对误差为9.3744%, 预报结果与观测值相关性达到了0.96678. 结果表明, 利用人工神经网络方法提前1h预报电离层TEC有很好的应用前景.      

基于轨道预报的热层模式评估

在太阳活动高低年的地磁平静/扰动环境下,利用不同热层大气模式J77,DTM78,MSIS00,JB2008和CHAMP加速度计反演密度,分析有无先验信息条件下的轨道预报误差.结果表明无先验信息的精密轨道预报中,热层模式的性能可能被弹道系数等参数偏差干扰,此时预报误差不能作为模式性能的评价标准.先验信息对轨道预报精度提升非常明显,尤其是地磁扰动期先进热层模式性能得以展现,轨道预报误差为无先验信息情况下的10%～25%.目前热层模式的主要缺陷存在于地磁扰动期.各模式之间的差异是:JB2008模式可以通过线性和单一频率周期项补偿,而J77及DTM78等模式还存在更多频率的误差.本文对不同情况下精密轨道预报的研究结果可为空间碎片碰撞预警等工程实践提供参考.      

Storm-time atmospheric density modeling using neural networks and its application in orbit propagation

Upper atmospheric densities during geomagnetic storms are usually poorly estimated due to a lack of clear understanding of coupling mechanisms between the thermosphere and magnetosphere. Consequently, the orbit determination and propagation for low-Earth-orbit objects during geomagnetic storms have large uncertainties. Artificial neural networks are often used to identify nonlinear systems in the absence of rigorous theory. In the present study, an attempt has been made to model the storm-time atmospheric density using neural networks. Considering the debate over the representative of geomagnetic storm effect, i.e. the geomagnetic indices ap and Dst, three neural network models (NNM) are developed with ap, Dst and a combination of ap and Dst respectively. The density data used for training the NNMs are derived from the measurements of the satellites CHAMP and GRACE. The NNMs are evaluated by looking at: (a) the mean residuals and the standard deviations with respect to the density data that are not used in training process, and (b) the accuracy of reconstructing the orbits of selected objects during storms employing each model. This empirical modeling technique and the comparisons with the models NRLMSIS-00 and Jacchia-Bowman 2008 reveal (1) the capability of neural networks to model the relationship between solar and geomagnetic activities, and density variations; and (2) the merits and demerits of ap and Dst when it comes to characterizing density variations during storms.