北斗掩星事件数量与分布的模拟研究

北斗卫星导航系统作为我国自主开发的卫星导航系统,在掩星探测领域有着广泛的应用前景.文章针对利用LEO星载掩星接收机进行北斗掩星探测的任务,在建立LEO卫星轨道模拟系统和北斗全星座模拟系统的基础上,通过仿真计算研究了不同LEO轨道参数条件下北斗掩星事件的特点.分析了对于单颗LEO卫星,北斗掩星事件的数量和分布随LEO轨道参数包括轨道近地点角距、升交点赤经、轨道高度和倾角而变化的规律.针对北斗导航系统由GEO、IGSO和MEO三种轨道卫星组成的特点,对不同类型北斗卫星的掩星事件进行了研究,并分别总结了三种轨道北斗卫星掩星事件的特点.研究结果对利用北斗导航系统进行掩星探测有参考作用.     

基于下一代四星转轮式编队系统精确和快速反演FSCF地球重力场

第一,由于重力卫星编队轨道的稳定性设计是建立下一代高精度和高空间分辨率地球重力场模型的关键,因此为保证下一代四星转轮式编队系统的稳定性,轨道根数的最优设计如下:(1)轨道半长轴a、轨道偏心率e、轨道倾角i和升交点赤经Ω保持不变;(2)每对卫星的近地点幅角ω和平近点角M分别相差180°;(3)初始近地点辐角ω设置于赤道处,初始平近点角M设计于极点处;(4)卫星编队系统椭圆轨道的半长轴和半短轴之比为2:1. 第二,基于下一代四星转轮式编队系统,利用星间速度插值法,通过相关系数(激光干涉测量系统的星间速度0.85、GPS接收机的轨道位置和轨道速度0.95、星载加速度计的非保守力0.90)、观测时间30天和采样间隔10 s,反演了120阶FSCF-1/2/3/4(Four-Satellite Cartwheel Formation)地球重力场,在120阶处累计大地水准面精度为1.162×10^-4 m,较目前GRACE地球重力场精度至少提高一个数量级. 第三,下一代四星转轮式编队系统具有低轨道高度、高精度测量、全张量观测、弱混频效应和强时变信号的优点.     

电离层GPS掩星观测改正TEC反演方法

电离层掩星观测中,当低轨卫星（LEO）轨道高度较低时,轨道以上的电离层电子总含量（TEC）对掩星反演的影响不能忽略.目前,一般采用指数函数等外推方法来处理该问题,对反演结果可能引起较大误差.为提高电离层掩星反演精度,本文研究利用LEO处于非掩星一侧GPS观测数据的改正TEC新反演方法.用三维射线追踪程序计算出电离层掩星观测模拟数据,分别应用改正TEC方法和外推方法进行反演,将反演结果与所用模式值进行比较.结果表明：对于轨道高度约800km的GPS/MET掩星模拟数据,外推方法和改正TEC方法反演结果都与模式值基本一致;对于轨道高度约400km的CHAMP掩星模拟数据,外推方法误差较大,改正TEC方法反演结果与模式值相符得较好.将改正TEC方法应用于GPS/MET实测数据的反演,取得了合理的结果.这些说明,改正TEC算法是一种有效的电离层掩星反演方法,尤其是对于轨道较低的LEO的电离层掩星观测反演特别有用.     

基于“张衡一号”卫星的GRO掩星载荷数据分析

"张衡一号"电磁卫星于2018年2月2日成功发射,该卫星搭载的全球导航卫星系统无线电掩星(GRO)接收机具备了全球定位系统(GPS)卫星和北斗导航系统(BD)卫星的掩星观测功能。对在轨测试期间数据进行分析,GRO掩星电离层反演结果合理、趋势正确,每天能够探测600个左右的电离层掩星事件。对连续2个多月共3万多个GRO掩星事件进行全球覆盖性分析,并与美国COSMIC掩星反演结果对比,GRO掩星与COSMIC掩星均能实现全球覆盖并且N_mF_2和H_mF_2的全球分布情况相近,随着纬度而变化明显。利用不同纬度的3个数字测高仪对时空一致的GRO掩星数据进行对比验证和分析,GRO掩星的N_mF_2相对数字测高仪的相对RM_SE优于9.41%,相关系数优于0.868 2。H_mF_2的相对RM_SE优于7.80%,相关系数优于0.706 6。     

低轨卫星轨道误差对中性层延迟量影响的模拟研究

GPS掩星监测大气方法中需要载有GPS接收机的卫星的精密定轨信息，而该卫星的轨道 精度对GPS掩星监测大气的效果进行分析是十分必要的. 基于地球大气模式CIRA1986，本文采用三维射线跟踪方法模拟了有代表性的5次完整的GPS掩星事件，得到了卫星的轨道精度 对GPS掩星测量中关键参数的影响的初步结果. 结果表明，该卫星的常规轨道误差对掩星测量的影响较小，完全可以采用差分定位方法来获取掩星测量中所需要的卫星定轨信息.     

GPS测量误差对大气掩星反演精度影响分析

在无线电GPS大气掩星观测中,GPS和LEO卫星的位置误差、速度误差、卫星钟的漂移、载波相位观测误差等会给其反演的大气参量带来一定的误差.为定量分析上述GPS测量误差对大气掩星反演精度的影响,本文采用模拟分析的方法,利用CHAMP实测轨道数据和大气与电离层经验模式,采用三维射线追踪方法模拟得到大气附加相位观测数据;然后根据各种GPS测量误差的不同性质,在模拟掩星观测大气附加相位数据中加入GPS测量误差并进行反演,将其反演结果与未加入误差时的反演结果进行比较,分析各测量误差所带来的掩星反演误差的特点;利用一天的151个分布全球的掩星事件的误差结果进行统计分析.结果表明:(1)±10 m的GPS径向位置误差引起的温度反演误差的平均误差在30 km的高度为(-+)0.8K,标准偏差为0.05K;±0.4 m的LEO径向位置误差引起的温度反演误差的平均误差在30 km的高度为(-+)1.2K,标准偏差为0.07K.(2)仅考虑卫星速度误差对掩星几何关系的影响,则±30 m/s的GPS速度误差引起的温度反演误差的平均误差最大为±0.02K,标准偏差为0.3K;±10 m/s的LEO速度误差引起的温度反演误差的平均误差最大为±0.4K,标准偏差为0.06K.(3)当卫星钟漂为1×10-9时,反演温度误差的平均误差在30 km的高度为263K,标准偏差为16K,当卫星钟漂为1×10-14时,反演温度误差的平均误差在30km的高度为0.019K,标准偏差为0.004K.(4)当L1和L2载波相位观测误差分别为1 mm和2 mm时,引起的温度反演的平均误差在30 km的高度为0.01K,标准偏差为0.4K.(5)利用较接近真实观测的三维大气和电离层背景下的模拟附加相位数据,同时加入国际最新指标的各种GPS测量误差进行反演,误差分析结果为:在30 km以下,折射率和密度的相对误差小于0.3%(1σ),压强的相对误差小于0.5%(1σ),温度误差小于1K(1σ).     

基于COSMIC资料分析电离层F层不规则体结构

根据电离层不规则体的产生会导致周围电子浓度发生起伏变化的原理,利用2007年COSMIC掩星系统的TEC数据,通过平滑滤波得到TEC的扰动值ΔTEC的变化,利用其研究F层不规则体的时空变化特征.统计结果表明:扰动较大的掩星事件主要发生在磁纬±20°之间和高纬地区,春季和秋季带状分布较为明显,不同经度地区较强扰动的掩星事件的分布也有不同特征;较强ΔTEC的掩星事件主要发生在地方时午夜前和午夜后两个时段,发生的高度主要在250~400km范围内.这些结果与已知的F层不规则体的时空分布特征较为一致,说明利用TEC的扰动量来分析电离层F层不规则体结构是可行的.     

Abel电离层掩星反演方法及误差分析

GNSS-LEO电离层无线电掩星技术是近年来发展的电离层探测新技术.为消除LEO轨道以上的电离层影响,改正TEC反演方法采用非掩星侧的观测数据进行电离层掩星反演.本文首次提出了一种新方法--基于历元差分的电离层反演方法;并将改正TEC与历元差分两种反演方法应用于模拟掩星观测数据反演,随后基于反演结果及误差分析得到一些有益的结论:历元差分反演精度较改正TEC反演精度均有所提高;不管是哪种方法,高轨(约800 km)反演结果要优于低轨(约500 km)的反演结果;随着剖面高度的降低,反演精度随之下降;反演误差主要集中在8至18时(当地时),主要分布在磁纬-30°至30°之间.     

MAGSAT地磁卫星及其应用

美国宇航局（NASA）1979年10月30日在西海岸成功地发射了一颗地磁卫星MAGSAT,卫星的近地点和远地点高度分别为352km和561km,轨道倾角为96.76°,运行周期为93.726 min。卫星上载有磁通门式向量磁强计,每秒取样16次,Cs蒸汽标量磁强计,每秒取样8次。为了避免卫星磁场的影响,磁强计装置在6m横杆的末端。卫星姿态由两个测星相机和一个太阳传感器进行测量,姿态测量精度为20″。这样在300-     

ZH-1卫星得到的中国及邻区顶部电离层背景特征

利用ZH-1(CSES)卫星LAP载荷原位电子密度数据对中国及邻区(0°-54°N,70°-140°E)的顶部电离层背景分布及随季节变化进行了详细分析,研究结果显示:(1)研究区赤道异常的纬度延伸范围、随经度分布形态及它们的季节变化,具有与其他研究结果一致的规律性.(2)中纬度区,白天电子密度存在一个低值带,夜间电子密...     

Global sounding of sporadic E layers by the GPS/MET radio occultation experiment

The GPS radio occultation technique is sensitive for layered structures with horizontal scales of around 100 km and with vertical scales of a few 100 m or more at the Earth's limb. These structures cause strong fluctuations of the GPS L1 and L2 phase paths which have been measured by a GPS receiver onboard of Microlab-1 satellite in 730 km orbit during the GPS/Meteorology experiment (GPS/MET of UCAR, Boulder). By means of GPS/MET radio occultation data, profiles of electron density fluctuations are derived for the mesosphere/lower thermosphere region with a height resolution of around 1 km. Data analysis of 1900 radio occultation events in June/July 1995, 1540 events in October 1995, and 2690 events in February 1997 confirms seasonal dependence of sporadic E layers. The meridian slices of average sporadic E activity show a dominance of plasma irregularities in the summer hemisphere. The irregularities mainly occur at heights 90–110 km. Auroral and equatorial sporadic E, electron density depletions, and multiple ionization layers are also present in the high resolution GPS/MET data. The multiple layers often have a distance of around 5–10 km in height, and appear up to a height of 140 km (upper height limit for 50 Hz sampling rate of GPS receiver). For February and June, the GPS/MET observations are compared to ground-based observations of the Asia/Australia ionosonde chain.     

Radio occultation data analysis by the radioholographic method

The radioholographic method is briefly described and tested by using data of 4 radio occultation events observed by the GPS/MET experiment on 9 February 1997. The central point of the radioholographic method (Pavelyev, 1998) is the generation of a radiohologram along the LEO satellite trajectory which allows the calculation of angular spectra of the received GPS radio wave field at the LEO satellite. These spectra are promising in view of detection, analysis and reduction of multipath/diffraction effects, study of atmospheric irregularities and estimation of bending angle error. Initial analysis of angular spectra calculated by the multiple signal classification (MUSIC) method gives evidence that considerable multibeam propagation occurs at ray perigee heights below 20 km and at heights around 80–120 km for the 4 GPS/MET occultation events. Temperature profiles obtained by our analysis (radioholographic method, Abel inversion) are compared with those of the traditional retrieval by the UCAR GPS/MET team (bending angle from slope of phase front, Abel inversion). In 3 of 4 cases we found good agreement (standard deviation σ_T∼1.5°K between both retrievals at heights 0–30 km).     

A methodology for using optimal MSIS parameters retrieved from SSULI data to compute satellite drag on LEO objects

A key application to be derived from Space Weather research will be to forecast atmospheric drag on low Earth orbit (LEO) satellites with significantly better accuracy than is attainable today. The recently launched STP P91-1 ARGOS mission will serve as a testbed for the use of future operational ultraviolet remote-sensing data to achieve such an improvement. This paper describes the associated methodology, which uses discrete inverse theory in conjunction with the data to derive correction factors in near real time for the MSISE-90 empirical thermospheric model. To simulate the application of this technique to orbit prediction, we use the Jacchia-71 operational model to generate an evolving “ground truth” upper atmospheric state over a 48 h time period. This permits a state-of-the-art Satellite Tool Kit orbit propagator to synthesize a corresponding “ground truth” orbit on a standard LEO test object at 350 km altitude. Our tests show that, for orbit prediction, the “data-enhanced” MSIS density specification can provide significant improvement over the uncorrected MSIS specification. However, for orbit prediction, the results are sensitive to the strategy selected for applying the correction factors. We contrast our results for orbit prediction with those of Marcos et al. (1998. Astrodynamics. Vol. 97(1). AAS, San Diego, pp. 501–513) for precision orbit determination. An important result in the context of Space Weather is that the Jacchia and MSIS models can show significant point-to-point disagreement, which has major implications for operational specification of thermospheric drag.     

利用COSMIC低轨卫星对GPS信号的顶部TEC观测资料研究等离子体层电子含量的变化特征

本文尝试利用COSMIC低轨卫星对GPS信号的顶部TEC观测资料研究等离子体层电子含量(简称PEC)的变化规律.首先介绍从低轨卫星对GPS的顶部TEC观测资料提取等离子体层垂直电子含量的方法,然后利用该方法提取2008年全年的PEC数据,进而研究了2008年这一太阳活动低年PEC随地磁纬度(MLAT)、磁地方时(MLT)以及不同季节的变化规律.此外,还利用提取的120°E和300°E经度链上的数据对比研究了PEC的经度变化情形.研究结果表明:(1)PEC主要集中分布在磁赤道±45°之间的一个绕地球的环带状区域中;(2)PEC表现出以下的昼夜变化规律特征:白天时段之值高于夜间,约在12—16MLT之间达到最高峰值,而最小PEC值出现在日出前大约4—5MLT左右的时段;(3)相比其他季节月份而言,PEC在北半球夏季月份(5—8月)具有最小值;(4)PEC存在明显的经度变化,不同经度链上的PEC存在不同的季节变化特征.     

近地空间环境的GNSS无线电掩星探测技术

从GPS/MET计划开始,基于GNSS的无线电掩星技术已成为一种强大的近地空间环境探测手段.截至到目前,已经有20多颗发射的低轨道卫星带GPS掩星接收机,其中COSMIC是首个专门用于掩星探测的卫星星座.这些掩星数据被广泛应用于气象预报、气候与全球变化研究、及空间天气监测和电离层研究.由于COSMIC的成功,相关合作单位目前正积极推动COSMIC-2计划,该计划将总共有12颗卫星,于2016年与2019年各发射6颗.COSMIC-2将携带一个高级的GNSS掩星接收机,它将接受GPS与GLONASS信号,并具备接受其他可获得信号源的能力(如中国北斗定位信号),其每日观测的掩星数量将是COSMIC的4~6倍.同时COSMIC-2还将携带两个空间天气载荷,加强空间天气的监测能力.本文以COSMIC与COSMIC-2计划为主线,对掩星的发展历史、技术要点进行了简单介绍,并简要综述了COSMIC取得的部分科学成果,同时对未来包括技术发展和众多的掩星观测进行了展望.     

Properties of the sunspot latitudinal distribution skewness

The properties of the sunspot latitudinal distributions related to skewness have been studied based on the data of the extended Greenwich catalog for 1874–2011. The results of the performed analysis indicate that a significant skewness is present in most annual latitudinal distributions of the sunspot index. In this case, the distribution skewness increases near the 11-year cycle maximum phase. An increase in the sunspot group number is also accompanied by an increase in skewness. In particular, when the sunspot index is large, the number of groups located below midlatitudes is mostly larger than the number of groups above these latitudes and this imbalance increases with increasing total sunspot activity level. In medium and large 11-year cycles, the average distribution skewness for a cycle is always positive and its value is related to the cycle amplitude. This results agree with the theoretical models of the 11-year cycle, where the specific features of the low-latitude meridional circulation are related to the sunspot activity level.     

地基GNSS VTEC约束的大气掩星电离层误差修正方法

介绍了一种考虑电离层沿GNSS掩星射线路径分布不对称且兼顾一阶和二阶项的大气掩星电离层误差修正新方法,它综合地基GNSS VTEC的水平变化信息和电离层模式的垂直变化信息,在沿"入射线"与"出射线"双边局部球对称假设下,估算GNSS掩星射线路径上的电子密度;进而计算一阶和二阶项弯曲角电离层误差廓线.采用太阳活动低年的2008年7月15日和太阳活动较高年的2013年7月15日两天的MetOp-A和GRACE掩星观测资料和IGS的GNSS VTEC数据,计算了GNSS大气掩星弯曲角一阶和二阶项电离层误差廓线.对比分析表明：新方法经验模型和理论模型的二阶项弯曲角电离层残差,以及经验模型与实测数据的一阶和二阶项弯曲角电离层误差均具有良好的一致性,因此该方法可用于L2信号数据质量较差或L2信号中断的掩星事件,同时修正大气掩星一阶和二阶电离层误差,从而提高弯曲角精度和掩星观测资料的利用率.