基于垂直剖面函数式的全球对流层天顶延迟模型的建立

对流层延迟是无线电导航定位的主要误差源之一,其值对目标高程的变化敏感.在动态导航定位中,由于目标高程变化随机性强,延迟改正实时性需求高,已有的对流层延迟模型难以满足应用需求.本文利用2005到2006年ERA-Interim再分析气象资料积分方法计算的对流层天顶总延迟(ZTD)、天顶静力学延迟(ZHD)以及天顶非静力学延迟(ZWD)的垂直剖面研究了ZTD随高程变化的最佳拟合形式,并以此为基础建立了全球ZTD改正模型SHAO-H.该模型以大气中水汽的垂直分布规律为依据,将ZTD表示为高程的分段函数,进而再模制每段函数中各参数随时间的变化.精度评估显示:与积分ZTD相比,SHAO-H模型计算的ZTD在不同等压层上的平均bias大部分在±1 mm以内,随着高度的上升,平均RMS由39 mm减小至不足1 mm;与IGS (International GNSS Service)实测ZTD相比,SHAO-H模型的精度(bias为7.02 mm,RMS为38.50 mm)优于UNB3m模型(bias为14.67 mm, RMS为51.95 mm).SHAO-H模型具有精度稳定、计算简便等优点,适宜任意高度的用户使用.     

一种新的全球对流层天顶延迟模型GZTD

对流层延迟是GNSS导航定位主要误差源之一,主要受气象参数(如总气压、温度和水汽压等)的影响,具有变化随机性强的特点.本文利用 GGOS Atmosphere提供的2002-2009年全球天顶对流层延迟格网时间序列研究了全球对流层天顶延迟的时空变化特征.并以此为基础对全球天顶对流层延迟(Zenith Troposphere Delay, ZTD)进行建模,提出了一种基于球谐函数的全球非气象参数对流层天顶延迟改正模型--GZTD模型.实验对比结果表明考虑ZTD经纬向变化的GZTD模型内符合精度全球统计结果(bias:0.2 cm,RMS:3.7 cm)优于只考虑ZTD纬向变化的UNB3m (bias:3.4 cm,RMS:6.0 cm)、UNB4 (bias:4.7 cm,RMS:7.4 cm)、UNB3 (bias:4.0 cm,RMS:7.0 cm)和EGNOS (bias:4.5 cm,RMS:6.9 cm)等模型.使用全球385个IGS站进行外符合检验,统计结果表明GZTD模型(bias:-0.02 cm,RMS:4.24 cm)同样优于其它模型.GZTD模型具有改正效果良好、使用简单、所需参数少等优点.     

一种新的高精度全球对流层天顶延迟模型

对流层延迟是影响高精度卫星导航定位的关键因素,也是大气科学研究的重要数据.针对已有全球对流层延迟模型的模型方程未同时顾及高程、纬度和季节变化以及模型构建时仅使用单一格网点数据等问题,本文提出了一种对流层天顶延迟(ZTD)全球模型构建的新方法,即引入滑动窗口算法将全球剖分为大小一致的规则窗口,利用2008—2015年全球大地观测系统(GGOS)大气格网产品构建每个窗口同时顾及高程、纬度和季节因子的全球ZTD新模型(GGZTD模型).联合未参与建模的2016年全球GGOS格网产品和2016年全球316个IGS站精密ZTD产品,检验了GGZTD模型的精度和适用性.结果表明:以GGOS大气格网ZTD产品和IGS站ZTD产品为参考值,GGZTD模型在全球的精度分别为3.58 cm和3.62 cm,相对于UNB3m模型和目前标称精度最优的GPT2w模型计算的ZTD信息,GGZTD模型在全球表现出了最优的精度和稳定性,其精度相对于UNB3m模型具有显著的提升(精度提高了30%以上),相对于GPT2w模型仍具有一定的改善;在ZTD计算时GGZTD模型相对于GPT2w模型显著地减少了模型参数,尤其相对于GPT2w-1(减少了99%).GGZTD模型只需输入位置与时间和依赖相对较少的模型参数则能在全球获得高精度和稳定的ZTD信息,极大地提升了模型的计算效率.     

基于GNSS的再分析资料对流层延迟精度评估

对流层延迟对导航定位精度有着重要的影响,而再分析资料提供的高精度气象参数计算的对流层延迟可应用于定位过程中以提升定位精度.本文针对三种再分析资料计算的对流层延迟进行精度评估,并将其应用在精密单点定位中,分析其对定位精度的影响.首先,利用2020年全球范围内125个IGS(International GNSS Service)站的对流层天顶总延迟(Zenith Total Delay, ZTD)作为真值对三种再分析资料(ERA5、MERRA2、CRA40)计算的ZTD进行了精度评估,并分析其时空分布特性.研究结果表明：ERA5-ZTD的均方根误差(RMS)最小(12.1 mm),其次为CRA40-ZTD(15.8 mm)和MERRA2-ZTD(16.9 mm),整体上ERA5-ZTD的精度最高;据所选的IGS站点的比较结果发现赤道平均偏差(BIAS)呈现负值,在中高纬度地区CRA40的精度优于MERRA2,在低纬度地区则相反,而ERA5在各纬度平均精度均为最优;当考虑季节因素时,三者计算的ZTD-RMS在夏秋季较大,其中ERA5的RMS季节变化最稳定.之后还利用180个探空站点对三者计算...     

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

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

对流层天顶延迟模型研究进展及其在中国区域的精度分析

对流层延迟是空间大地测量技术的主要误差源之一,数据处理中对流层延迟的修正需要借助对流层延迟模型.本文首先从物理原理出发,梳理了对流层天顶延迟模型的研究历程和最新进展.按照时间顺序,对流层延迟模型先后经历了依赖实测气象资料的经典模型、不依赖实测参数的经验模型和以数值模型气象资料为基础的高分辨经验模型三个发展阶段.其次,本文利用中国区域内219个GPS测站2014—2015年两年实测的天顶延迟,对后两类经验模型中国际最新通用的代表模型UNB3m和GPT2w在中国境内的实际精度进行评估.精度评估结果显示:UNB3m模型在中国地区的平均Bias为-0.85 cm,平均RMSE为5.14 cm,其精度不随计算时间分辨率的变化而显著变化;模型参数的空间分辨率对GPT2w模型在中国地区的精度的影响不大,但GPT2w模型精度随计算时间分辨率的提高显著下降,2 h分辨率时GPT2w模型的平均RMSE分别为8.07 cm(1°参数文件)和7.97 cm(5°参数文件),1天分辨率时GPT2w模型的平均RMSE分别为3.49 cm(1°参数文件)和3.59 cm(5°参数文件);受水汽分布的影响,时间上,两个模型在冬季的精度相对最高,在夏季的精度相对最差,空间上,两个模型在高纬度和高海拔地区的精度相对较高.以上分析可为中国区域用户对流层延迟模型的选择提供参考.     

区域实时对流层网格模型构建及在精密单点定位中的应用

受制于对流层延迟建模方法和建模背景场精度及时空分辨率的影响,目前实时对流层延迟模型的精度和稳定性都有待进一步改善.本文利用甘肃及周围地基共计184个GNSS (Global Navigation Satellite System)站估算的ZTD (Zenith Troposphere Delay),构建了空间分辨率为0.25°×0.25°的甘肃地区实时ZTD网格模型.针对传统的高程归化模型及水平内插模型精度低的问题,本文提出了利用高斯指数函数模型将不同高程的GNSS/ZTD归化到统一的高度,再基于BP神经网络模型从网格顶点周围统一高度后的GNSS/ZTD中内插出网格顶点处的ZTD.为了验证甘肃ZTD网格模型的精度,选取2022年甘肃地区8个未参与建模的陆态网GNSS测站的数据进行了实验.统计结果显示：与事后PPP (Precise Point Positioning)处理GNSS估算的ZTD相比,甘肃ZTD网格模型与真值偏差的RMS优于1.52 cm.此外,将构建的实时ZTD格网模型用于约束PPP处理,对于PPP浮点解施加ZTD约束后U方向精度提升22.9%,U方向收敛时间缩短26.4...     

基于实时精密单点定位技术的暴雨短临预报

提出了一种将实时精密单点定位(Precise Point Positioning, PPP)技术用于暴雨短临预报的新方法.该方法首先基于GPS连续运行参考站网(Continuously Operating Reference Stations, CORS)实时估计的精密卫星钟差完成PPP解算,再以实时获取的对流层延迟(Zenith Tropospheric Delay, ZTD)及其增量变化为依据进行暴雨短临预报.研究结果表明:一般雷暴天气来临之前的2~6 h,ZTD增量表现为先后突破±5 mm/5 min,且后续记录到的实际降水量大小与ZTD维持在高水平阶段的时间长短有较好的对应关系;就热带气旋而言,在强风作用下,ZTD增量变化表现的异常活跃和复杂,规律性较弱,但对短临预报强降雨仍有一定的指示作用.     

两种精化的对流层延迟改正模型

对流层延迟是全球导航卫星系统(Global Navigation Satellite System,GNSS)导航定位中的重要误差源,其量值主要受气象条件影响.采用传统对流层建模思路,利用GPT2模型来提供相对准确的气温、气压和相对湿度,然后利用Saastamoinen模型来计算天顶对流层延迟,由此构建了GPT2+Saas模型;采用新的对流层建模思路,直接针对天顶对流层延迟的时空特性建模,构建了与GPT2+Saas模型相匹配的GZTDS格网模型.以GGOS Atmosphere格网数据为参考,GPT2+Saas模型(Bias:0.2cm;RMS:4.2cm)和GZTDS模型(Bias:0.2cm;RMS:3.7cm)较UNB3m模型精度分别提升34%和43%.以IGS(International GNSS Service)数据为参考,GPT2+Saas(Bias:0.5cm;RMS:4.7cm)和GZTDS(Bias:-0.3cm;RMS:3.8cm)相对UNB3m模型精度分别提升10%和27%.针对GPT2+Saas模型在少数测站出现精度异常的情况进行了研究,探讨了可能的原因.在两种不同思路构建的精化对流层模型中,GZTDS模型不仅表现出更高的精度,而且在时间稳定性和地理稳定性上也表现出优越性.     

基于改进的BP神经网络构建区域精密对流层延迟模型

利用神经网络算法挖掘海量数据的规律已成为科技发展的一种趋势,本文针对卫星信号的天顶对流层延迟进行建模.对流层延迟是影响卫星定位精度的重要因素之一,建立精密区域对流层模型对高精度定位有着重要的意义.对区域测站对流层延迟数据的分析,考虑到实时建模中传统BP（Back Propagation）神经网络计算量大,易出现"过拟合"现象、不稳定等因素,通过改进的BP神经网络建立了区域精密对流层模型.详细介绍了新模型的建立过程,并与常用的对流层区域实时模型进行了对比.还讨论了建模测站数目对预报精度的影响.相比现有的其他对流层延迟模型,基于改进的BP神经网络构建的区域精密对流层延迟模型无论在拟合和预报方面都有较好的精度,且随着测站数目的增加模型精度趋于平稳.改进的模型参数较少,可以进行实时的区域精密对流层延迟改正;需要播发的信息量小,适用于连续运行参考站系统（Continuously Operating Reference Stations,CORS）的应用.研究表明：改进的BP神经网络模型能够更好的充分利用大规模历史数据描述卫星信号对流层延迟的空间分布情况,适用于实时大区域精密对流层建模.基于日本地区2005年近1000多个测站的NCAR（National Center Atmospheric Research）对流层数据进行区域对流层延迟建模,结果表明改进的BP神经网络模型在拟合和预报精度上都有较大提升,RMSE（Root Mean Square Error）分别为：7.83 mm和8.52 mm,而四参数模型拟合、预报RMSE分别18.03 mm和16.60 mm.     

Signature of atmospheric oscillations in GPS-measured tropospheric delay

Ground-based GPS finds potential applications in many atmospheric studies such as spatial distribution of columnar water vapor as well as the tidal oscillations in the atmosphere. The zenith tropospheric delay (ZTD) derived from GPS data at two Indian IGS stations are used to establish its potential for studying the atmospheric tidal, intra-seasonal and planetary oscillations. The major tidal oscillations observed in ZTD data are diurnal, semi-diurnal and their harmonics. Prominent intra-seasonal oscillations observed in ZTD are reported for the first time in this context. These intra-seasonal oscillations are Madden–Julian Oscillation (30–70 days, 60–90 days, 100–120 days) and planetary waves (like 27, 16 and 5–10 days periodicities). Quantification of these periodicities will provide a useful handle to improve the empirical models employed in the estimation of tropospheric delay.     

A global empirical model for estimating zenith tropospheric delay

Tropospheric delay acts as a systematic error source in the Global Navigation Satellite Systems (GNSS) positioning. Empirical models UNB3, UNB3m, UNB4 and EGNOS have been developed for use in Satellite-Based Augmentation Systems (SBAS). Model performance, however, is limited due to the low spatial resolution of the look-up tables for meteorological parameters. A new design has been established in this study for improving performance of the tropospheric delay model by more effectively eliminating the error produced by tropospheric delay. The spatiotemporal characteristics of the Zenith Tropospheric Delay (ZTD) were analyzed with findings that ZTD exhibits different annual variations at different locations and decreases exponentially with height increasing. Spherical harmonics are utilized based on the findings to fit the annual mean and amplitude of the ZTD on a global scale and the exponential function is utilized for height corrections, yielding the ZTrop model. On a global scale, the ZTrop features an average deviation of -1.0 cm and Root Mean Square (RMS) of 4.7 cm compared with the International GNSS Service (IGS) ZTD products, an average deviation of 0.0 cm and RMS of 4.5 cm compared with the Global Geodetic Observing System (GGOS) ZTD data, and an average deviation of -1.3 cm and RMS of 5.2 cm compared with the ZTD data from the Constellation Observing System of Meteorology, Ionosphere, and Climate (COSMIC). The RMS of the ZTrop model is 14.5% smaller than that of UNB3, 6.0% smaller than that of UNB3m, 16% smaller than that of UNB4, 14.5% smaller than that of EGNOS and equivalent to the sophisticated GPT2+Saas model in comparison with the IGS ZTD products. The ZTrop, UNB3m and GPT2+Saas models are finally evaluated in GPS-based Precise Point Positioning (PPP), as the models act to aid in obtaining PPP position error less than 1.5 cm in north and east components and relative large error (>5 cm) in up component with respect to the random walk approach.     

Assessment of GPS data for meteorological applications over Africa: Study of error sources and analysis of positioning accuracy

The aim of this study is to assess the availability and quality of data from the International GNSS Service (IGS) Global Positioning System (GPS) network in Africa, especially for retrieving zenith tropospheric delay (ZTD), from which precipitable water vapour (PWV) can be derived, in view of application to the African Monsoon Multidisciplinary Analysis (AMMA) project. Three major error sources for the GPS data analysis evaluating PWV in Africa are the accuracy of the satellite orbits, the correction for the radio delay induced by the ionosphere and the vertical site displacements due to ocean loading. The first part of this study examines these error sources and the validity of GPS data for meteorological applications in Africa in dedicated analyses spanning the year 2001. These analyses were performed using the IGS precise orbits. Weak degradation of baseline precision with increasing baseline lengths suggests that the average orbital error is not limiting the GPS analysis in Africa. The impact of the ionosphere has been evaluated during a maximum of solar activity in 2001. The loss of L₂ data has actually been observed. It amounts to 2% on average for 2001, with maxima of 8% during magnetic storm events. A slight decrease in formal accuracy of ZTD seems to be related to the loss of L₂ data at the end of the day. This indicates that scintillation effects are present in the GPS observations but however are not a major limitation. The impact of ocean loading is found to be significant on ZTD estimates (up to ±2 mm in equivalent PWV). The use of a proper ocean loading model eliminates this effect.The second aspect of this study concerns the IGS analysis quality for the African stations. The accuracy has been assessed through position dispersion between individual solutions and the most recent version of the IGS combined solution IGb00, and residuals from the transformation of the IGS combined solution into the International Terrestrial Reference Frame 2005. The positioning performance of the IGS analysis is consistent with an accuracy in ZTD of ±6 mm (±1 mm in PWV), as requested for meteorological applications such as planned in AMMA.     

4D tropospheric tomography using GPS slant wet delays

Tomographic techniques are successfully applied to obtain 4D images of the tropospheric refractivity in a local dense network of global positioning system (GPS) receivers. We show here how GPS data are processed to obtain the tropospheric slant wet delays and discuss the validity of the processing. These slant wet delays are the observables in the tomographic processing. We then discuss the inverse problem in 4D tropospheric tomography making extensive use of simulations to test the system and define the resolution and the impact of noise. Finally, we use data from the Kilauea network in Hawaii for February 1, 1997, and a local 4 × 4 × 40 voxel grid on a region of 400 km² and 15 km in height to produce the corresponding 4D wet refractivity fields, which are then validated using forecast analysis from the European Center for Medium Range Weather Forecast (ECMWF). We conclude that tomographic techniques can be used to monitor the troposphere in time and space.     

Inversion, error analysis, and validation of GPS/MET occultation data

The global positioning system meteorology (GPS/MET) experiment was the first practical demonstration of global navigation satellite system (GNSS)-based active limb sounding employing the radio occultation technique. This method measures, as principal observable and with millimetric accuracy, the excess phase path (relative to propagation in vacuum) of GNSS-transmitted radio waves caused by refraction during passage through the Earth’s neutral atmosphere and ionosphere in limb geometry. It shows great potential utility for weather and climate system studies in providing an unique combination of global coverage, high vertical resolution and accuracy, long-term stability, and all-weather capability. We first describe our GPS/MET data processing scheme from excess phases via bending angles to the neutral atmospheric parameters refractivity, density, pressure and temperature. Special emphasis is given to ionospheric correction methodology and the inversion of bending angles to refractivities, where we introduce a matrix inversion technique (instead of the usual integral inversion). The matrix technique is shown to lead to identical results as integral inversion but is more directly extendable to inversion by optimal estimation. The quality of GPS/MET-derived profiles is analyzed with an error estimation analysis employing a Monte Carlo technique. We consider statistical errors together with systematic errors due to upper-boundary initialization of the retrieval by a priori bending angles. Perfect initialization and properly smoothed statistical errors allow for better than 1 K temperature retrieval accuracy up to the stratopause. No initialization and statistical errors yield better than 1 K accuracy up to 30 km but less than 3 K accuracy above 40 km. Given imperfect initialization, biases ≫ 2 K propagate down to below 30 km height in unfavorable realistic cases. Furthermore, results of a statistical validation of GPS/MET profiles through comparison with atmospheric analyses of the European Centre for Medium-range Weather Forecasts (ECMWF) are presented. The comparisons indicate the high utility of the occultation data in that very good agreement of upper troposphere/lower stratosphere temperature (better than 1.5 K rms, ≪ 0.5 K bias) is found for a region (Europe+USA) where the ECMWF analyses are known to be good, but poorer agreement for a region (Southern Pacific) where the analyses are known to be degraded.     

Preliminary assessment of the navigation and positioning performance of BeiDou regional navigation satellite system

BeiDou regional navigation satellite system(BDS)also called BeiDou-2 has been in full operation since December 27,2012.It consists of 14 satellites,including 5 satellites in Geostationary Orbit(GEO),5 satellites in Inclined Geosynchronous Orbit(IGSO),and 4 satellites in Medium Earth Orbit(MEO).In this paper,its basic navigation and positioning performance are evaluated preliminarily by the real data collected in Beijing,including satellite visibility,Position Dilution of Precision(PDOP)value,the precision of code and carrier phase measurements,the accuracy of single point positioning and differential positioning and ambiguity resolution(AR)performance,which are also compared with those of GPS.It is shown that the precision of BDS code and carrier phase measurements are about 33 cm and 2 mm,respectively,which are comparable to those of GPS,and the accuracy of BDS single point positioning has satisfied the design requirement.The real-time kinematic positioning is also feasible by BDS alone in the opening condition,since its fixed rate and reliability of single-epoch dual-frequency AR is comparable to those of GPS.The accuracy of BDS carrier phase differential positioning is better than 1 cm for a very short baseline of 4.2 m and 3 cm for a short baseline of 8.2 km,which is on the same level with that of GPS.For the combined BDS and GPS,the fixed rate and reliability of single-epoch AR and the positioning accuracy are improved significantly.The accuracy of BDS/GPS carrier phase differential positioning is about 35 and 20%better than that of GPS for two short baseline tests in this study.The accuracy of BDS code differential positioning is better than 2.5 m.However it is worse than that of GPS,which may result from large code multipath errors of BDS GEO satellite measurements.