基于单频星载GPS数据的低轨卫星精密定轨

为满足搭载单频GPS接收机低轨卫星的精密定轨需求以及深化单频定轨研究,文中解决了单频星载GPS数据的周跳探测问题,并利用“海洋二号”（HY-2A）卫星及“资源三号”（ZY-3）卫星的单频星载GPS实测数据采用两种方法确定了二者的简化动力学轨道,并通过观测值残差分析、与双频精密轨道比较、激光测卫数据检核等方法对所得轨道精度进行评定。结果表明,在不考虑电离层延迟影响的情况下,HY-2A卫星定轨精度为2~3dm,ZY-3卫星为1m左右;而采用半和改正组合消除电离层延迟一阶项影响后,二者定轨精度均显著提高,HY-2A卫星三维精度提高至1dm左右,ZY-3卫星提高至1~2dm。文章的研究成果表明,搭载单频GPS接收机的低轨卫星也可获得厘米级的定轨精度。     

基于FY-3C掩星数据偶发E层的研究

利用FY-3C极轨卫星提供的2014年6月至2015年5月的GPS无线电掩星数据,统计分析了全球范围内抽样频率为50 Hz的C/A码SNR扰动情况,进而对偶发E层进行了研究.结果表明:偶发E层在夏季半球中纬地区的扰动强度远远大于冬季半球同一纬度地区的扰动强度,偶发E层在纬度40°附近扰动明显增强;在E层100 km高度附近,Es层在10:00 LT和22:00 LT达到峰值;Es层在夏季半球的出现率明显高于冬季半球;FY-3C卫星的掩星观测结果与COSMIC系统的观测结果较一致,可以利用FY-3C卫星的掩星数据研究电离层偶发E层等的变化.     

利用COSMIC掩星数据监测电离层的异常变化

分析了COSMIC掩星数据反演电子密度的方法,利用实例研究反演方法的特点,并采用ISR非相干散射雷达获取的电子密度数据进行验证,进而反演了长三角区域SHAO（IGS）站上空在日全食和太阳风暴期间的电子密度廓线图. 通过与平静日期间电离层电子密度进行比较,发现日全食及太阳风暴导致电离层发生的异常变化,从而提出COSMIC掩星数据反演电子密度在监测电离层变化时所具有的优势.     

月球电离层掩星探测研究

日本SELENE/KAGUYA探测任务提供了研究月球电离层的机会。采用无线电掩星探测技术和趋势外推算法,消除地球电离层和行星际等离子体的干扰影响,残余的信号相位信息的变化反映了月球电离层的信息,估算出月球周围附近近似对称分布的稀薄电离层中电子总含量(TEC)约为每立方米10-14个。     

月球卫星的精密定轨

发射月球卫星是探测月球的一个重要途径,而且为了找水,多采用极轨方式.由于作为中心天体的月球与地球有所不同,导致月球卫星的运动及其相应的轨道摄动解亦有明显差异.在这种背景下,本文选择采用一组适用于0≤e<1和0°≤i<180°的无奇点根数、充分利用数值法和分析法各自特点的定轨方法,给出月球卫星的精密定轨方案,并以测距的模拟资料进行了实算,计算结果与理论分析相符,有实用前景.     

基于固定模糊度PPP的低轨卫星实时精密定轨

传统动力学定轨法受制于动力学模型精度,传统几何定轨法精度受限,只能达到亚米级,而基于精密单点定位（PPP）模式的几何定轨法一般采用浮点解,定轨精度及可靠性较基于双差模式的相对定位较差。为提高PPP模式低轨定轨的定位性能,利用中国区域内外的IGS测站计算出当前所有卫星的宽巷和窄巷相位小数偏差产品,对经过中国大陆区域上空的国产低轨卫星海洋二号（HY-2）和资源三号 (ZY-3) 卫星进行固定模糊度PPP的定轨解算,与事后精密轨道结果进行比较,分析其外符合精度。结果表明：仅利用约10min弧段的HY-2和ZY-3卫星数据,切向与径向的定轨精度可达2cm左右,法向为5cm左右,较浮点解定轨精度大幅提升。基于固定模糊度PPP的定轨方法能够满足厘米级的实时精密定轨。     

QX-1 GNOS M掩星数据质量分析

QX-1 GNOS M是首台在轨实现北斗、GPS和Galileo三系统兼容的小型商业化全球导航卫星掩星探测仪. 2021年10月14日18:51 LT, 气象一号(QX-1)卫星发射并入轨. 自发射以来, 该卫星已收集大量观测数据. 基于QX-1 GNOS M的结构组成及性能, 统计分析了2022年8月17日全天的掩星事件及其全球分布情况. 通过将8月17日至9月3日的掩星数据与NCEP再分析模式对比, 评估QX-1掩星事件的探测穿透深度和折射率精度, 同时检验Galileo掩星数据的可靠性和一致性. 初步分析结果表明, QX-1 GNOS M在实现三系统兼容后, 掩星事件数量相较于仅使用GPS系统的情况增加了约1.5倍. 这一结果进一步证明, 不同全球导航卫星系统(GNSS)所提供的掩星数据在精度上具有一致性. 此外研究显示, 在多系统兼容的背景下, QX-1 GNOS M能够提供更为丰富和精确的气象数据.     

地基行星大气掩星观测资料流程和数据整理

在中国火星探测萤火一号(YH-1)计划中, 包括了地基掩星观测反演火星大气的科研任务. 观测资料整理是反演流程的第一步. 本文描述了地基火星大气掩星观测处理软件系统的观测数据流程和观测资料整理模块,并详细介绍了观测资料整理模块的结构和功能, 其中包括时间系统转换、历表插值、坐标系变换、信号时延改正以及掩星平面建立. 利用行星数据系统公布的火星快车无线电科学数据和由SPICE得到的地球、火星历表以及火星快车的轨道数据, 结合本文的算法, 得到了一些实验结果.      

基于光学测角数据的风云四号同步轨道卫星精密定轨

针对风云四号同步卫星的精密定轨和精度评估需求,首先利用地面光学测角数据对FY-4A卫星进行精密定轨,定轨后方位角和高度角的残差rms分别为0.25＂和0.45＂。与基于测距数据的轨道相比,位置精度在有测角数据的弧段内小于50m。进一步联合测角数据和测距数据对FY-4A卫星进行联合定轨,定轨后轨道重叠精度优于15m。利用联合定轨结果评估了基于测距数据的实时轨道产品精度,可以明显发现轨道精度随着测距数据的积累而逐步提高。     

平流层爆发性增温期间大气变化特性的COSMIC掩星观测

利用高精度和高垂直分辨率的COSMIC掩星观测资料, 详细深入分析了2007年冬---2008年春平流层爆发性增温(SSW)期间10~60 km高度范围内大气的变化特性, 尤其是上平流层和低中间层大气的变化特性. 结果表明, 在SSW过程中, 温度场、风场和剩余环流都发生了明显的变化. 根据温度在主增温前和主增温盛期的变化特性, 在水平方向, 大约以55ºN为界, 在垂直方向, 大约以42 km为界, 可以将温度场在纬度-高度的分布分为4个区域: 高纬下层增温区, 增温幅度约高达25 K; 高纬上层降温区, 降温幅度约达30,K; 中纬下层降温区, 降温幅度约为几K; 中纬上层增温区, 增温也约为几K. SSW期间上下层大气纬向风场的变化规律基本相同. 在纬度方向以45ºN为界, 45ºN以北地区的西风减弱东风增强, 风场变化高达50 m/s; 45ºN以南地区西风增强东风减弱, 变化幅度比较小, 约10 m/s. 在2008年1月下旬到2月底, 大气温度和纬向风有明显的振荡现象, 周期约为12天. 剩余环流的环流圈在SSW期间会发生反转, 由此也表明, SSW期间大气中物质的输运方向也会发生改变.     

DORIS geodesy: A dynamic determination of geocentre location

Geoscience Australia contributed a multi-satellite, multi-year weekly time series to the International DORIS Service combined submission for the construction of International Terrestrial Reference Frame 2008 (ITRF2008). This contributing solution was extended to a study of the capability of DORIS to dynamically estimate the variation in the geocentre location. Two solutions, comprising different constraint configurations of the tracking network, were undertaken. The respective DORIS satellite orbit solutions (SPOT-2, SPOT-4, SPOT-5 and Envisat) were verified and validated by comparison with those produced at the Goddard Space Flight Center (GSFC), DORIS Analysis Centre, for computational consistency and standards. In addition, in the case of Envisat, the trajectories from the GA determined SLR and DORIS orbits were compared. The results for weekly dynamic geocentre estimates from the two constraint configurations were benchmarked against the geometric geocentre estimates from the IDS-2 combined solution. This established that DORIS is capable of determining the dynamic geocentre variation by estimating the degree one spherical harmonic coefficients of the Earth’s gravity potential. It was established that constrained configurations produced similar results for the geocentre location and consequently similar annual amplitudes. For the minimally constrained configuration Greenbelt–Kitab, the mean of the uncertainties of the geocentre location were 2.3, 2.3 and 7.6 mm and RMS of the mean uncertainties were 1.9, 1.2 and 3.5 mm for the X, Y and Z components, respectively. For GA_IDS-2_Datum constrained configuration, the mean of the uncertainties of the geocentre location were 1.7, 1.7 and 6.2 mm and RMS of the mean uncertainties were 0.9, 0.7 and 2.9 mm for the X, Y and Z components, respectively. The mean of the differences of the two DORIS dynamic geocentre solutions with respect to the IDS-2 combination were 1.6, 4.0 and 5.1 mm with an RMS of the mean 21.2, 14.0 and 31.5 mm for the Greenbelt–Kitab configuration and 4.1, 3.9 and 4.3 mm with an RMS 8.1, 9.0 and 28.6 mm for the GA_IDS-2_Datum constraint configuration. The annual amplitudes for each component were estimated to be 5.3, 10.8 and 11.0 mm for the Greenbelt–Kitab configuration and 5.3, 9.3 and 9.4 mm for the GA_IDS-2_Datum constraint configuration. The two DORIS determined dynamic geocentre solutions were compared to the SLR determined dynamic solution (which was determined from the same process of the GA contribution to the ITRF2008 ILRS combination) gave mean differences of 3.3, −4.7 and 2.5 mm with an RMS of 20.7, 17.5 and 28.0 mm for the X, Y and Z components, respectively for the Greenbelt–Kitab configuration and 1.1, −5.4 and 4.4 mm with an RMS of 9.7, 13.3 and 24.9 mm for the GA_IDS-2_Datum configuration. The larger variability is reflected in the respective amplitudes. As a comparison, the annual amplitudes of the SLR determined dynamic geocentre are 0.9, 1.0 and 6.8 mm in the X, Y and Z components. The results from this study indicate that there is potential to achieve precise dynamically determined geocentre from DORIS.     

Adjustable box-wing model for solar radiation pressure impacting GPS satellites

One of the major uncertainty sources affecting Global Positioning System (GPS) satellite orbits is the direct solar radiation pressure. In this paper a new model for the solar radiation pressure on GPS satellites is presented that is based on a box-wing satellite model, and assumes nominal attitude. The box-wing model is based on the physical interaction between solar radiation and satellite surfaces, and can be adjusted to fit the GPS tracking data.     

卫星与测站钟差支持条件下的GEO卫星精密定轨

在基于伪距的GEO卫星精密定轨中, GEO卫星的静地特性导致定轨解算无法对星地组合钟差进行有效估计, 需要独立的时间同步支持. 本文讨论了卫星和测站钟差支持条件下的GEO卫星定轨原理, 利用仿真数据系统地分析了中国区域网跟踪条件下GEO卫星的定轨精度, 从定性和定量角度分析了钟差二次项、星地时间同步精度、站间时间同步精度及系统差等因素对定轨精度的影响.     

GPS single-frequency orbit determination for low Earth orbiting satellites

The determination of high-precision orbits for Low Earth Orbiting (LEO) satellites (e.g., CHAMP, GRACE, MetOp-A) is based on dual-frequency tracking data from on-board GPS (Global Positioning System) receivers. The two frequencies allow it to eliminate the first order ionosphere effects. Data screening and precise orbit determination (POD) procedures are optimized under the assumption of the availability of two frequencies.     

The International DORIS Service (IDS): Toward maturity

DORIS is one of the four space-geodetic techniques participating in the Global Geodetic Observing System (GGOS), particularly to maintain and disseminate the Terrestrial Reference Frame as determined by International Earth rotation and Reference frame Service (IERS). A few years ago, under the umbrella of the International Association of Geodesy, a DORIS International Service (IDS) was created in order to foster international cooperation and to provide new scientific products. This paper addresses the organizational aspects of the IDS and presents some recent DORIS scientific results. It is for the first time that, in preparation of the ITRF2008, seven Analysis Centers (AC’s) contributed to derive long-term time series of DORIS stations positions. These solutions were then combined into a homogeneous time series IDS-2 for which a precision of less than 10 mm was obtained. Orbit comparisons between the various AC’s showed an excellent agreement in the radial component, both for the SPOT satellites (e.g. 0.5–2.1 cm RMS for SPOT-2) and Envisat (0.9–2.1 cm RMS), using different software packages, models, corrections and analysis strategies. There is now a wide international participation within IDS that should lead to future improvements in DORIS analysis strategies and DORIS-derived geodetic products.     

2016HO3探测任务星载光学观测量建模及定轨

以中国首次小行星探测任务为背景,根据星载相机获取的光学影像构造三种观测量,分别为小行星相对于航天器的高度角和方位角、赤经赤纬以及探测器与行星/小行星之间的夹角,分析了其在探测器定轨中的作用。仿真定轨结果表明,观测时长为100 h,探测器三轴位置误差小于50 km,满足工程上对巡航段的轨道精度要求,但x和y方向的位置和速度分量具有较强自相关性。此外还发现,使用单一观测数据类型比联合观测量的定轨精度低3~4个量级,第三类观测量相对于其他两类观测量在定轨精度方面具有显著作用,这表明在2016HO3探测中,利用太阳系大天体的位置信息有助于约束探测器轨道,提高探测器的定轨精度。