基于星载GPS的HY-2卫星高精度精密定轨模拟研究(英文)

HY-2卫星是我国第一颗测高卫星,其径向定轨精度要求厘米量级,搭载了星载GPS接收机。目前HY-2还处于测试阶段,没有公布观测数据。为了确定基于星载GPS的HY-2精密定轨流程及其定轨精度,本文模拟了HY-2卫星星载GPS观测数据,结果表明HY-2星载GPS天线每个历元至少观测7颗GPS卫星。给出了基于星载GPS的精密定轨流程,分别采用简化动力学方法和动态几何法进行了精密定轨实验。对于相位1mm和3mm随机误差的相位观测数据,简化动力学法和动态几何法定轨都能够实现厘米量级的径向精密定轨,几何法定轨精度略低于简化动力定轨。地球重力场模型是影响HY-2卫星精密定轨的重要因素,本文对不同阶次的重力场模型EIGEN2、EGM96、TEG4和GEMT3进行了简化动力学定轨实验,高于50阶次的重力场模型都能够实现厘米级径向精密定轨,主要原因在于大量的高精度星载GPS观测数据和重力场模型精度的提高。     

Centimeter-level precise orbit determination for the HY-2A satellite using DORIS and SLR tracking data

The HY-2A satellite is the first ocean dynamic environment monitoring satellite of China. Centimeter-level radial accuracy is a fundamental requirement for its scientific research and applications. To achieve this goal, we designed the strategies of precise orbit determination (POD) in detail. To achieve the relative optimal orbit for HY-2A, we carried out POD using DORIS-only, SLR-only, and DORIS + SLR tracking data, respectively. POD tests demonstrated that the consistency level of DORIS-only and SLR-only orbits with respect to the CNES orbits were about 1.81 cm and 3.34 cm in radial direction in the dynamic sense, respectively. We designed 6 cases of different weight combinations for DORIS and SLR data, and found that the optimal relative weight group was 0.2 mm/s for DORIS and 15.0 cm for SLR, and RMS of orbit differences with respect to the CNES orbits in radial direction and three-dimensional (3D) were 1.37 cm and 5.87 cm, respectively. These tests indicated that the relative radial and 3D accuracies computed using DORIS + SLR data with the optimal relative weight set were obviously higher than those computed using DORIS-only and SLR-only data, and satisfied the requirement of designed precision. The POD for HY-2A will provide the invaluable experience for the following HY-2B, HY-2C, and HY-2D satellites.     

基于多颗低轨卫星的SLR台站坐标几何解算方法

卫星激光测距(SLR)技术作为卫星精密定轨手段和轨道检核重要方法,激光反射器已经成为重力卫星和测高卫星等低轨卫星的基本载荷.经典的SLR台站坐标是使用动力学方法计算的,本文根据多颗低轨卫星(LEO)多历元的激光观测数据,采用几何方法开展地面SLR测站坐标计算.通过组建低轨卫星群实现对全球激光站的动态观测,为了合理配置不同低轨卫星间观测值权重,削弱低轨卫星群可能存在的系统性偏差,提出采用方差分量估计组合的最小二乘法进行解算.实测结果显示,解算出SLR台站坐标框架解与SLRF2014差异平均值在25.1 mm,外符合精度达到1~2 cm.该方法避免了复杂动力学模型,SLR台站坐标的几何计算方法既可以作为激光测站框架解算手段之一,同时将LEO卫星群作为空间并址站实现不同技术地球参考框架间的融合.

     

基于多颗低轨卫星的SLR台站坐标几何解算方法

卫星激光测距（SLR）技术作为卫星精密定轨手段和轨道检核重要方法,激光反射器已经成为重力卫星和测高卫星等低轨卫星的基本载荷.经典的SLR台站坐标是使用动力学方法计算的,本文根据多颗低轨卫星（LEO）多历元的激光观测数据,采用几何方法开展地面SLR测站坐标计算.通过组建低轨卫星群实现对全球激光站的动态观测,为了合理配置不同低轨卫星间观测值权重,削弱低轨卫星群可能存在的系统性偏差,提出采用方差分量估计组合的最小二乘法进行解算.实测结果显示,解算出SLR台站坐标框架解与SLRF2014差异平均值在25.1 mm,外符合精度达到1~2 cm.该方法避免了复杂动力学模型,SLR台站坐标的几何计算方法既可以作为激光测站框架解算手段之一,同时将LEO卫星群作为空间并址站实现不同技术地球参考框架间的融合.     

Improvement of ERS-1 orbits using along-track accelerations from DORIS data on SPOT2

In long-arc precise orbit determinations of altimetric satellites such as ERS-1, large errors may occur from mismodelling of aerodynamic drag and solar radiation pressure. Such surface forces for non-spherical satellites require accurate modelling of the effective area and particle-surface interactions, but the dominant source of error is neutral air density as derived from thermospheric models for aerodynamic drag. Several techniques can be employed to alleviate air-drag mismodelling but all require the solution of additional parameters from the tracking data. However, for ERS-1 the sparsity of laser range data limits the application of such empirical techniques. To overcome this, use can be made of the dense DORIS Doppler tracking for SPOT2 which is in a similar orbit to ERS-1. A recent investigation by CNES examined the use of drag scale factors from SPOT2 to constrain the ERS-1 orbit. An improvement to that methodology is to consider along-track mismodelling as observed by timing errors in the Doppler data for each pass of SPOT2. The along-track correction to the acceleration as derived from SPOT2 can then be applied to ERS-1 orbits, solving for a scale factor to absorb systematic errors - particularly that arising from the 50 km altitude difference. Results are presented of the associated improvement in ERS-1 orbits as derived from concurrent SPOT2 arcs. It will be seen that the procedure not only improves the laser range fit, but more importantly, leads to more precise radial positioning as evident in the altimeter and crossover residuals.     

APOD mission status and preliminary results

On September 20 th, 2015, twenty satellites were successfully deployed into a near-polar circular orbit at 520 km altitude by the Chinese CZ-6 test rocket, which was launched from the Tai Yuan Satellite Launch Center. Among these satellites, a set of 4 Cube Sats conform the atmospheric density detection and precise orbit determination(APOD) mission, which is projected for atmospheric density estimation from in-situ detection and precise orbit products. The APOD satellites are manufactured by China Spacesat Co. Ltd. and the payload instruments include an atmospheric density detector(ADD), a dual-frequency dualmode global navigation satellite system(GNSS) receiver(GPS and Beidou), a satellite laser ranging(SLR) reflector, and an S/Xband very long baseline interferometry(VLBI) beacon. In this paper, we compare the GNSS precise orbit products with colocated SLR observations, and the 3 D orbit accuracy shows better than 10 cm RMS. These results reveal the great potential of the onboard micro-electro-mechanical system(MEMS) GNSS receiver. After calibrating ADD density estimates with precise orbit products, the accuracy of our density products can reach about 10% with respect to the background density. Density estimates from APOD are of a great importance for scientific studies on upper atmosphere variations and useful for model data assimilation.     

海潮误差对GRACE时变重力场解算的影响研究

海潮误差是 GRACE 时变重力场反演中重要的误差源,目前发布的海潮模型中主要包含振幅较大的主潮波分量模型,在时变重力场反演中次潮波的影响也是不可忽略的,因此,GRACE 时变重力场反演中的海潮误差主要包括受限于海潮模型误差和次潮波影响.本文利用轨道模拟方法检测了短周期潮波的混频周期以及次潮波对ΔC20, ΔC30的时序特征,并进一步通过轨道模拟结果分析了海潮误差对时变重力场反演的影响,然后通过实测数据解算分析了海潮误差对当前 GRACE 时变重力场解算的影响,研究发现:(1) 利用轨道模拟能够有效地检测短周期潮波的混频周期;(2)时变重力场解算过程中,次潮波的影响大于海潮模型误差的影响;(3)海潮模型误差以及次潮波影响是当前 GRACE 没有达到基准精度的重要因素之一.     

Derivation Of Refractive Index And Temperature Gradients From Optical Scintillometry To Correct Atmospherically Induced Errors For Highly Precise Geodetic Measurements

Refraction effects of optical beams are generally caused by an inhomogeneous propagation medium and are a major source of systematic errors in the precise optical determination of angles and distances in the atmospheric surface layer. In this contribution a method for deriving vertical temperature and refractive index gradients from optical scintillation is presented. Knowledge of these gradients is required for the compensation of atmospherically induced errors for highly precise terrestrial geodetic measurements, like direct transfer and levelling. The advantage of the present optical method is, that temperature and refractive index gradients can be derived as line-averaged values over the propagation path, which is not possible by meteorological point measurements. Field observations have been carried out with a displaced-beam scintillometer over flat terrain and under different atmospheric conditions in order to verify this method. The experiments show, that this method allows to derive accurate correction values for precise terrestrial geodetic measurements.     

联合多代卫星测高和多源重力数据的局部大地水准面精化方法

本文研究了基于泊松小波径向基函数融合多代卫星测高及多源重力数据精化大地水准面模型的方法.分别以沿轨垂线偏差和大地水准面高高差作为卫星测高观测量,研究了使用不同类型测高数据对于大地水准面建模精度的影响.针对全球潮汐模型在浅水区域及部分开阔海域精度较低的问题,引入局部潮汐模型研究了不同潮汐模型对于大地水准面的影响.数值分析表明:相比于使用沿轨垂线偏差作为测高观测量,基于沿轨大地水准面高高差解算得到的大地水准面模型的精度更高,特别是在海域区域,其精度提高了2.3 cm.由于使用沿轨大地水准面高高差作为测高观测量削弱了潮汐模型长波误差的影响,采用不同潮汐模型对大地水准面解算的影响较小.总体而言,船载重力及测高观测数据在海洋重力场的确定中呈现互补性关系,联合两类重力场观测量可以提高局部重力场的建模精度.     

卫星重力梯度测量与地球引力场的精度研究

本文根据地球引力位的球谐函数展开式,利用重力梯度张量各分量导出了位系数模型的精度估计公式.从三方面进行了研究:假定卫星重力梯度仪测量精度,探讨用重力梯度数据确定地球重力场模型的精度;求出位系数模型和大气阻力引起的重力梯度卫星的轨道误差;最后,反求轨道误差和位系数误差对重力梯度测量值的影响.数值计算表明,与地面技术和常规卫星方法相比,卫星梯度测量可使重力场模型的精度至少提高3—5倍;利用重力梯度张量全分量求得的重力值精度比单用径向分量Vrr的结果提高40％以上;若仅顾及位系数模型和大气阻力误差,则轨道误差对梯度测量值的影响△Vi3（i=3,2,1）至少可分别在1/4和1/3弧圈内达到△Vi3≤σ（仪器精度）.     

卫星重力梯度测量与地球引力场的精度研究

本文根据地球引力位的球谐函数展开式,利用重力梯度张量各分量导出了位系数模型的精度估计公式.从三方面进行了研究:假定卫星重力梯度仪测量精度,探讨用重力梯度数据确定地球重力场模型的精度;求出位系数模型和大气阻力引起的重力梯度卫星的轨道误差;最后,反求轨道误差和位系数误差对重力梯度测量值的影响.数值计算表明,与地面技术和常规卫星方法相比,卫星梯度测量可使重力场模型的精度至少提高3-5倍;利用重力梯度张量全分量求得的重力值精度比单用径向分量V_rr的结果提高40％以上;若仅顾及位系数模型和大气阻力误差,则轨道误差对梯度测量值的影响△V_i3（i=3,2,1）至少可分别在1/4和1/3弧圈内达到△V_i3≤σ（仪器精度）.     

Regional ocean tide loading modelling around the Iberian Peninsula

We developed a new 1/12° resolution oceanic tide model in the complex region that surrounds the Iberian Peninsula. The model, named IBER01, allows us to obtain more accurate tidal loading computations for precise geodetic and gravimetric observations in this area. The modelling follows the scheme of data assimilation (coastal tide gauge, bottom pressure sensors and TOPEX/Poseidon altimetry) into a hydrodynamical model, which is based on two-dimensional barotropic depth averaged shallow-water equations. Detailed bathymetry data and quadratic bottom friction with a specific drag coefficient for the region have been considered. Improved ocean load maps for the Iberian Peninsula are obtained for eight harmonic constituents (Q1, P1, O1, K1, N2, M2, S2 and K2), after computing the load effect (Newtonian attraction and elastic contribution) using IBER01 and six present-day global oceanic tide models for comparison. The results achieved verify the quality of the new model. Our ocean loading computations reduce considerably the discrepancies between the theoretical Earth tide parameters and those from observations at the level of 0.3%.     

固定非差整数模糊度的PPP快速精密定位定轨

从GPS基本观测模型出发,给出并推导了分离相位小数偏差求解非差整数模糊度的精密单点定位数学模型和算法.利用少量IGS跟踪站组成服务端观测网计算未检校的相位小数偏差改正信息,用于改正用户端接收机的相位观测值,实现了固定非差整数模糊度的快速精密单点定位与定轨.试验结果表明: 利用30 min的地面动态或静态观测数据进行精密单点定位,非差模糊度固定成整数后,其定位结果较PPP浮点解有明显改善,水平方向提高了近一个数量级,可达到1 cm甚至毫米级的精度;高程方向与对流层延迟解算精度也改善了20%~60%.与浮点解相比,固定解能显著改善PPP的定轨精度,仅用15 min的短弧段观测数据,切向与法向的定轨精度可达到1 cm左右;径向方向为3~5 cm左右,较浮点解定轨精度改善了50%~70%.因此,固定非差整数模糊度后的PPP能够满足毫米至厘米级的快速精密定位和定轨的要求,这在GPS(准)实时应用与服务中具有很好的应用前景.     

Modelling the fine-structure of the geoid: Methods,data requirements and some results

The requirements for precise geoid models on local and regional scales have increased in recent years, primarily due to the ongoing developments in height determination by GPS on land, but also due to oceanographic requirements in using satellite altimetry for recovering dynamic sea-surface topography. Suitable methods for geoid computations from gravity data include Stokes integration, FFT methods, and least-squares collocation. Especially the FFT methods are efficient in handling large amounts of gravity data, and new variants of the methods taking earth curvature rigorously into account provide attractive methods for obtaining continental-scale, high-resolution geoid models. The accuracy of such models may be from 2–5 cm locally, to 50–100 cm on regional scales, depending on gravity data coverage, long wave-length gravity field errors, and datum problems. When approaching the cm-level geoid basic geoid definition questions (geoid or quasigeoid?) become very significant, especially in rugged areas. In the paper the geoid modelling methods and problems are reviewed, and some investigations on local data requirements for cm-level geoid prediction are presented. Some actual results are presented from Scandinavia, where a recent regional high-resolution geoid model yields apparent accuracies of 2–10 cm over GPS baselines of 50 to 2000 km.     

Contributions of Satellite Laser Ranging to Past and Future Radar Altimetry Missions

Satellite laser ranging (SLR) has proven avery efficient method for contributingto the tracking of altimetric satellites anddetermining accurately their orbitalthough hampered by the non-worldwide coverageand the meteorologicalconditions. Indeed, in some cases it is the onlymethod available to determinethe satellite orbit (e.g., the orbits of the ERS-1and Geosat-Follow-On missions).Moreover, any operational and non-weather dependenttechniques, like GPS,DORIS, PRARE, can exhibit systematic errors inpositioning and orbitography. Acomparison with SLR results allows to evidence sucherrors and vice versa. Fordoing that, two different approaches for determiningprecise orbits can beconsidered: one based on global orbit determination,the other on a short-arctechnique used to locally improve a global orbitdetermined by another trackingtechniques, such as DORIS or GPS. We can thusvalidate a global orbit andachieve orbit quality control to a level of2 to 3 centimeters at present and expectto achieve a level of 1 to 2 centimeters inthe near future. Errors induced bystation coordinates or by the gravity field(geographically correlated errors, forexample) can be estimated from SLR tracking data.Colocation experiments withdifferent techniques in the same geodetic siteplay also a key role to ensure preciserelationships between the geodetic referenceframes linked to each technique. Inparticular, the role of the SLR technique is tostrengthen the vertical component(including velocity) of the positioning, whichis crucial for altimetry missions.The role of SLR data in the modelling of the firstterms of the gravity field has finally to be emphasized,which is of primary importance in orbitography,whatever the tracking technique used.Another application of SLR technology is thesatellite altimeter calibration. Examples of past calibrationand future experiments are given, including theaccuracy we can expect from the Jason-1 and EnviSatspace oceanography missions.     

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.     

潮汐效应对沿海多条带时序InSAR地表形变监测的影响

潮汐效应对沿海大范围、高精度的合成孔径雷达干涉测量(Interferometric Synthetic Aperture Radar, InSAR)技术地表形变监测的影响不可忽略. 潮汐效应中的固体潮(Solid Earth Tide, SET)位移和海洋潮汐负荷(Ocean Tide Loading, OTL)位移会在时间和空间域上对沿海大范围InSAR地表形变监测产生分米级的误差, 而且不同条带中地表潮汐形变时空差异性会对多条带InSAR影像拼接产生较大误差. 本文针对沿海地区多条带InSAR形变时序, 重点分析了潮汐效应在不同条带InSAR地表沉降监测中的影响, 并采用多种潮汐位移估计方法对多条带地表潮汐形变进行改正. 结果表明, 研究区域不同条带时序InSAR中地表潮汐形变具有较大的时空差异性, 其造成沿海地区相邻条带时序InSAR地表沉降速率差异能达到1~2 cm·a^-1; 通过固体潮模型、海潮负荷模型或GPS参考站网海潮位移改正, 能够消除不同条带中地表潮汐形变时空差异性对时序InSAR形变结果拼接产生的空间高阶非线性误差, 弥补了传统拼接方法在拟合潮汐位移偏差的不足. 在大气延迟误差改正的基础上, 时序InSAR形变残差的标准差(Standard Deviation, STD)由潮汐改正前的2.3 cm减少至0.75 cm; 本文研究揭示了沿海地表潮汐形变改正对多条带时序InSAR形变结果拼接的重要性, 可大幅提高沿海地区广域时序InSAR地表形变监测的精度.

     

APOD卫星大气密度数据处理与标校

APOD卫星是我国首颗以热层大气密度探测与精密定轨为科学目标的微纳卫星,搭载大气密度探测器、双频GNSS接收机等载荷,于2015年9月20日发射入轨,2015年10月27日进入轨道高度460 km、轨道倾角97.4°、降交点地方时6:20的工作轨道,各项载荷随即展开例行观测.本文给出了APOD卫星大气密度探测器的基本原理和数据处理流程,采用基于双行根数（TLE）反演获取的密度数据,对2015年12月至2016年12月的就位探测数据进行了标校,并与经验密度模式进行了比较.结果表明,反演密度与APOD卫星就位探测数据的线性相关性达到0.943,采用线性拟合与二次函数拟合的残差水平基本相当.两种不同方法标校密度相对于NRLMSIS00模式日均值误差的均值和标准偏差为10.1%、18.2%和5.1%、17.1%,二次函数标校略优于线性标校;相对于JB2008模式日均值误差的均值和标准偏差为0.6%、14.9%和3.9%、16.9%,线性标校略优于二次函数标校.总体而言,APOD卫星大气密度就位探测数据与常用经验模式精度基本一致,可为开展大气密度变化规律及应用研究提供数据基础.

     

Progress in technology for the 2005 height determination of Qomolangma Feng (Mt. Everest)

Qomolangma Feng——Mt. Everest (QF in shorter) is located in the east part of the China-Nepal boundary. QF is the highest mountain in China, also in the world. It lies in the collision zone along the boundary of Eurasian and Indian plates. It is one of focuses in geo- sciences study[1―7]. The northern slope area of the QF (QF area in shorter) is in Tibet, China, and the altitude of snow line in the QF area is around 6000 m. Spring in the QF area is March, April and May every year, a…     

The Integrated Global Geodetic Observing System (IGGOS) viewed from the perspective of history

Towards the end of the 19th century, geodetic observation techniques allowed it to create geodetic networks of continental size. The insight that big networks can only be set up through international collaboration led to the establishment of an international collaboration called “Central European Arc Measurement”, the predecessor of the International Association of Geodesy (IAG), in 1864. The scope of IAG activities was extended already in the 19th century to include gravity.At the same time, astrometric observations could be made with an accuracy of a few tenths of an arcsecond. The accuracy stayed roughly on this level, till the space age opened the door for milliarcsecond (mas) astrometry. Astrometric observations allowed it at the end of the 19th century to prove the existence of polar motion. The insight that polar motion is almost unpredictable led to the establishment of the International Latitude Service (ILS) in 1899.The IAG and the ILS were the tools (a) to establish and maintain the terrestrial and the celestial reference systems, including the transformation parameters between the two systems, and (b) to determine the Earth's gravity field.Satellite-geodetic techniques and astrometric radio-interferometric techniques revolutionized geodesy in the second half of the 20th century. Satellite Laser Ranging (SLR) and methods based on the interferometric exploitation of microwave signals (stemming from Quasars and/or from satellites) allow it to realize the celestial reference frame with (sub-)mas accuracy, the global terrestrial reference frame with (sub-)cm accuracy, and to monitor the transformation between the systems with a high time resolution and (sub-)mas accuracy. This development led to the replacement of the ILS through the IERS, the International Earth Rotation Service in 1989.In the pre-space era, the Earth's gravity field could “only” be established by terrestrial methods. The determination of the Earth's gravitational field was revolutionized twice in the space era, first by observing geodetic satellites with optical, Laser, and Doppler techniques, secondly by implementing a continuous tracking with spaceborne GPS receivers in connection with satellite gradiometry. The sequence of the satellite gravity missions CHAMP, GRACE, and GOCE allow it to name the first decade of the 21st century the “decade of gravity field determination”.The techniques to establish and monitor the geometric and gravimetric reference frames are about to reach a mature state and will be the prevailing geodetic tools of the following decades. It is our duty to work in the spirit of our forefathers by creating similarly stable organizations within IAG with the declared goal to produce the geometric and gravimetric reference frames (including their time evolution) with the best available techniques and to make accurate and consistent products available to wider Earth sciences community as a basis for meaningful research in global change. IGGOS, the Integrated Global Geodetic Observing System, is IAG's attempt to achieve these goals. It is based on the well-functioning and well-established network of IAG services.