Tropospheric refractivity and zenith path delays from least-squares collocation of meteorological and GNSS data

Precise positioning requires an accurate a priori troposphere model to enhance the solution quality. Several empirical models are available, but they may not properly characterize the state of troposphere, especially in severe weather conditions. Another possible solution is to use regional troposphere models based on real-time or near-real time measurements. In this study, we present the total refractivity and zenith total delay (ZTD) models based on a numerical weather prediction (NWP) model, Global Navigation Satellite System (GNSS) data and ground-based meteorological observations. We reconstruct the total refractivity profiles over the western part of Switzerland and the total refractivity profiles as well as ZTDs over Poland using the least-squares collocation software COMEDIE (Collocation of Meteorological Data for Interpretation and Estimation of Tropospheric Pathdelays) developed at ETH Zürich. In these two case studies, profiles of the total refractivity and ZTDs are calculated from different data sets. For Switzerland, the data set with the best agreement with the reference radiosonde (RS) measurements is the combination of ground-based meteorological observations and GNSS ZTDs. Introducing the horizontal gradients does not improve the vertical interpolation, and results in slightly larger biases and standard deviations. For Poland, the data set based on meteorological parameters from the NWP Weather Research and Forecasting (WRF) model and from a combination of the NWP model and GNSS ZTDs shows the best agreement with the reference RS data. In terms of ZTD, the combined NWP-GNSS observations and GNSS-only data set exhibit the best accuracy with an average bias (from all stations) of 3.7 mm and average standard deviations of 17.0 mm w.r.t. the reference GNSS stations.     

A grid-based tropospheric product for China using a GNSS network

Tropospheric delay accounts for one source of error in global navigation satellite systems (GNSS). To better characterize the tropospheric delays in the temporal and spatial domain and facilitate the safety-critical use of GNSS across China, a method is proposed to generate a grid-based tropospheric product (GTP) using the GNSS network with an empirical tropospheric model, known as IGGtrop. The prototype system generates the GTPs in post-processing and real-time modes and is based on the undifferenced and uncombined precise point positioning (UU-PPP) technique. GTPs are constructed for a grid form (\(2.0{^{\circ }}\times 2.5{^{\circ }}\) latitude–longitude) over China with a time resolution of 5 min. The real-time GTP messages are encoded in a self-defined RTCM3 format and broadcast to users using NTRIP (networked transport of RTCM via internet protocol), which enables efficient and safe transmission to real-time users. Our approach for GTP generation consists of three sequential steps. In the first step, GNSS-derived zenith tropospheric delays (ZTDs) for a network of GNSS stations are estimated using UU-PPP. In the second step, vertical adjustments for the GNSS-derived ZTDs are applied to address the height differences between the GNSS stations and grid points. The ZTD height corrections are provided by the IGGtrop model. Finally, an inverse distance weighting method is used to interpolate the GNSS-derived ZTDs from the surrounding GNSS stations to the location of the grid point. A total of 210 global positioning system (GPS) stations from the crustal movement observation network of China are used to generate the GTPs in both post-processing and real-time modes. The accuracies of the GTPs are assessed against with ERA-Interim-derived ZTDs and the GPS-derived ZTDs at 12 test GPS stations, respectively. The results show that the post-processing and real-time GTPs can provide the ZTDs with accuracies of 1.4 and 1.8 cm, respectively. We also apply the GTPs in real-time kinematic GPS PPP, and the results show that the convergence time of the PPP solutions is shortened. These results confirm that the GTPs can act as an efficient information source to augment GNSS positioning over China.     

Comparative analysis of real-time precise point positioning zenith total delay estimates

The continuous evolution of global navigation satellite systems (GNSS) meteorology has led to an increased use of associated observations for operational modern low-latency numerical weather prediction (NWP) models, which assimilate GNSS-derived zenith total delay (ZTD) estimates. The development of NWP models with faster assimilation cycles, e.g., 1-h assimilation cycle in the rapid update cycle NWP model, has increased the interest of the meteorological community toward sub-hour ZTD estimates. The suitability of real-time ZTD estimates obtained from three different precise point positioning software packages has been assessed by comparing them with the state-of-the-art IGS final troposphere product as well as collocated radiosonde (RS) observations. The ZTD estimates obtained by BNC2.7 show a mean bias of 0.21 cm, and those obtained by the G-Nut/Tefnut software library show a mean bias of 1.09 cm to the IGS final troposphere product. In comparison with the RS-based ZTD, the BNC2.7 solutions show mean biases between 1 and 2 cm, whereas the G-Nut/Tefnut solutions show mean biases between 2 and 3 cm with the RS-based ZTD, and the ambiguity float and ambiguity fixed solutions obtained by PPP-Wizard have mean biases between 6 and 7 cm with the references. The large biases in the time series from PPP-Wizard are due to the fact that this software has been developed for kinematic applications and hence does not apply receiver antenna eccentricity and phase center offset (PCO) corrections on the observations. Application of the eccentricity and PCO corrections to the a priori coordinates has resulted in a 66 % reduction of bias in the PPP-Wizard solutions. The biases are found to be stable over the whole period of the comparison, which are criteria (rather than the magnitude of the bias) for the suitability of ZTD estimates for use in NWP nowcasting. A millimeter-level impact on the ZTD estimates has also been observed in relation to ambiguity resolution. As a result of a comparison with the established user requirements for NWP nowcasting, it was found that both the G-Nut/Tefnut solutions and one of the BNC2.7 solutions meet the threshold requirements, whereas one of the BNC2.7 solution and both the PPP-Wizard solutions currently exceed this threshold.     

Assessment of ZTD derived from ECMWF/NCEP data with GPS ZTD over China

The accuracy and feasibility of computing the zenith tropospheric delays (ZTDs) from data of the European Center for Medium-Range Weather Forecasts (ECMWF) and the United States National Centers for Environmental Prediction (NCEP) are studied. The ZTDs are calculated from ECMWF/NCEP pressure-level data by integration and from the surface data with the Saastamoinen model method and then compared with the solutions measured from 28 global positioning system (GPS) stations of the Crustal Movement Observation Network of China (CMONOC) for 1 year. The results are as follows: (1) the error of the integration method is 1–3 cm less than that of the Saastamoinen model method. The agreement between the ECMWF ZTD and GPS ZTD is better than that between NCEP ZTD and GPS ZTD; (2) the bias and root mean square difference (RMSD), especially the latter, have a seasonal variation, and the RMSD decreases with increasing altitude while the variation with latitude is not obvious; and (3) when using the full horizontal resolution of 0.5° × 0.5° of the ECMWF meteorological data in place of a reduced 2.5° × 2.5° grid, the mean RMSD between GPS and ECMWF ZTD decreases by 4.5 mm. These results illuminated the accuracy and feasibility of computing the tropospheric delays and establishing the ZTD prediction model over China for navigation and positioning with ECMWF and NCEP data.     

Estimation and evaluation of real-time precipitable water vapor from GLONASS and GPS

The revitalized Russian GLONASS system provides new potential for real-time retrieval of zenith tropospheric delays (ZTD) and precipitable water vapor (PWV) in order to support time-critical meteorological applications such as nowcasting or severe weather event monitoring. In this study, we develop a method of real-time ZTD/PWV retrieval based on GLONASS and/or GPS observations. The performance of ZTD and PWV derived from GLONASS data using real-time precise point positioning (PPP) technique is carefully investigated and evaluated. The potential of combining GLONASS and GPS data for ZTD/PWV retrieving is assessed as well. The GLONASS and GPS observations of about half a year for 80 globally distributed stations from the IGS (International GNSS Service) network are processed. The results show that the real-time GLONASS ZTD series agree quite well with the GPS ZTD series in general: the RMS of ZTD differences is about 8 mm (about 1.2 mm in PWV). Furthermore, for an inter-technique validation, the real-time ZTD estimated from GLONASS-only, GPS-only, and the GPS/GLONASS combined solutions are compared with those derived from very long baseline interferometry (VLBI) at colocated GNSS/VLBI stations. The comparison shows that GLONASS can contribute to real-time meteorological applications, with almost the same accuracy as GPS. More accurate and reliable water vapor values, about 1.5–2.3 mm in PWV, can be achieved when GLONASS observations are combined with the GPS ones in the real-time PPP data processing. The comparison with radiosonde data further confirms the performance of GLONASS-derived real-time PWV and the benefit of adding GLONASS to stand-alone GPS processing.     

基于GPS非差技术近实时获取对流层延迟

本文提出一种近实时获取高精度对流层延迟(Zenith Tropospheric Delay,ZTD)的方法。该方法基于全球GPS参考站网络的非差观测值数据和IGU预报产品(卫星轨道和ERP),使用GPS非差技术和卡尔曼滤波近实时估计各参考站的高精度ZTD。将该方法应用于香港GPS参考站网得到的ZTD与基于事后技术计算的ZTD进行比较后发现,两者的平均偏差均优于5mm,RMS均优于6mm。     

基于GNSS基准站网的对流层延迟建模

针对在卫星导航定位中,通常采用对流层模型进行,对流层延迟误差修正的现状,该文研究了一种基于GNSS基准站网的对流层延迟建模方法,并基于此方法利用日本地区GENET参考网约737个站5a的GNSS-ZTD序列建立了区域对流层模型ZTDM-JPN,并将ZTDM-JPN模型应用于GPS及北斗定位实验,分析了其在GPS及北斗定位中的实际应用性能。通过与国际上常用的对流层模型EGNOS、UNB3m作比较,结果表明,ZTDM-JPN模型的模拟精度较相同条件下的EGNOS与UNB3m模型分别提升约26%和21%,从而验证了该建模方法的可行性与优越性。     

Characterization of diurnal cycles in ZTD from a decade of global GPS observations

The diurnal cycle of the tropospheric zenith total delay (ZTD) is one of the most obvious signals for the various physical processes relating to climate change on a short time scale. However, the observation of such ZTD oscillations on a global scale with traditional techniques (e.g. radiosondes) is restricted due to limitations in spatial and temporal resolution. Nowadays, the International GNSS Service (IGS) provides an important data source for investigating the diurnal and semidiurnal cycles of ZTD and related climatic signals. In this paper, 10 years of ZTD data from 1997 to 2007 with a 2-hour temporal resolution are derived from global positioning system (GPS) observations taken at 151 globally distributed IGS reference stations. These time series are used to investigate diurnal and semidiurnal oscillations. Significant diurnal and semidiurnal oscillations of ZTD are found for all GPS stations used in this study. The diurnal cycles (24 hours period) have amplitudes between 0.2 and 10.9 mm with an uncertainty of about 0.5 mm and the semidiurnal cycles (12 h period) have amplitudes between 0.1 and 4.3 mm with an uncertainty of about 0.2 mm. The larger amplitudes of the diurnal and semidiurnal ZTD cycles are observed in the low-latitude equatorial areas. The peak times of the diurnal cycles spread over the whole day, while the peak value of the semidiurnal cycles occurs typically about local noon. These GPS-derived diurnal and semidiurnal ZTD signals are similar with the surface pressure tides derived from surface synoptic pressure observations, indicating that atmospheric tides are the main driver of the diurnal and semidiurnal ZTD variations.     

Variational data assimilation for deriving global climate analyses from GNSS radio occultation data

A comprehensive global navigation satellite system (GNSS) based radio occultation (RO) data set is available for meteorology and climate applications since the start of GNSS RO measurements aboard the CHAllenging Mini-satellite Payload (CHAMP) satellite in February 2001. Global coverage, all-weather capability, long-term stability and accuracy not only makes this innovative use of GNSS signals a valuable supplement to the data set assimilated into numerical weather prediction (NWP) systems but also an excellent candidate for global climate monitoring. We present a 3D variational data assimilation (3D-Var) scheme developed to derive consistent global analysis fields of temperature, specific humidity, and surface pressure from GNSS RO data. The system is based on the assimilation of RO data within 6 h time windows into European Centre for Medium-Range Weather Forecasts (ECMWF) short-term (24 h, 30 h) forecasts, to derive climatologic monthly mean fields. July 2003 was used as a test-bed for assessing the system’s performance. The results show good agreement with climatologies derived from RO data only and recent NWP impact studies. These findings are encouraging for future developments to apply the approach for longer term climatologic analyses, validation of other data sets, and atmospheric variability studies.     

利用中国区域MERRA-2资料计算ZTD和ZWD的精度分析

对流层延迟是影响高精度定位与导航的主要误差之一,也是全球导航卫星系统（global navigation satellite system,GNSS）水汽探测的关键参数。美国航空航天局发布了最新一代的大气再分析资料（MERRA-2资料）,其可用于计算高时空分辨率的对流层延迟产品,但是目前尚无文献对利用MERRA-2资料计算天顶对流层延迟（zenith tropospheric delay,ZTD）和天顶湿延迟（zenith wet delay,ZWD）的精度进行分析。因此,联合2015年中国陆态网214个GNSS站ZTD产品和分布于中国区域的87个探空站资料,对利用MERRA-2资料在中国区域计算ZTD/ZWD的精度进行评估。结果表明：（1）以陆态网ZTD为参考值,利用MERRA-2资料积分计算ZTD的年均偏差和均方根误差（root mean square error,RMSE）分别为0.32 cm和1.21 cm,且偏差和RMSE均表现出一定的季节变化,总体上呈现为夏季精度低、冬季精度高;在空间分布上,偏差随纬度和高程的变化趋势并不明显,但RMSE随纬度和高程的增加总体上呈现递减的趋...     

Troposphere delays from space geodetic techniques, water vapor radiometers, and numerical weather models over a series of continuous VLBI campaigns

Continuous, very long baseline interferometry (VLBI) campaigns over 2 weeks have been carried out repeatedly, i.e., CONT02 in October 2002, CONT05 in September 2005, CONT08 in August 2008, and CONT11 in September 2011, to demonstrate the highest accuracy the current VLBI was capable at that time. In this study, we have compared zenith total delays (ZTD) and troposphere gradients as consistently estimated from the observations of VLBI, Global Navigation Satellite Systems (GNSS), and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) at VLBI sites participating in the CONT campaigns. We analyzed the CONT campaigns using the state-of-the-art software following common processing strategies as closely as possible. In parallel, ZTD and gradients were derived from numerical weather models, i.e., from the global European Centre for Medium-Range Weather Forecasts (ECMWF) analysis fields, the High Resolution Limited Area Model (European sites), the Japan Meteorological Agency-Operational Meso-Analysis Field (MANAL, over Japan), and the Cloud Resolving Storm Simulator (Tsukuba, Japan). Finally, zenith wet delays were estimated from the observations of water vapor radiometers (WVR) at sites where the WVR observables are available during the CONT sessions. The best ZTD agreement, interpreted as the smallest standard deviation, was found between GNSS and VLBI techniques to be about 5–6 mm at most of the co-located sites and CONT campaigns. We did not detect any significant improvement in the ZTD agreement between various techniques over time, except for DORIS and MANAL. On the other hand, the agreement and thus the accuracy of the troposphere parameters mainly depend on the amount of humidity in the atmosphere.     

GGOS对流层延迟产品精度分析及在PPP中的应用

对流层延迟是卫星导航定位的主要误差源,GNSS广域增强需要高精度的对流层延迟产品进行误差修正。对流层延迟可通过GNSS进行实时估计,也可通过融合多源数据的数值气象预报模型获取。IGS发布的全球对流层天顶延迟产品由GNSS解算,其精度可达4mm,时间分辨率为5min,但其分布不均匀,在广袤的海洋区域无数据覆盖。GGOS Atmosphere基于ECMWF 40年再分析资料,可提供1979年以来时间分辨率为6h、空间分辨率为2.5°×2°的全球天顶对流层总延迟格网数据。本文通过2015年全球IGS测站的ZTD资料对GGOS的ZTD产品进行了评估,研究了GGOS Atmosphere对流层延迟产品与IGS发布ZTD资料之间的系统差,通过线性拟合估计出每个测站GGOS-ZTD与IGSZTD系统差系数(包括比例误差a和固定误差b),然后对比例误差a、固定误差b进行球谐展开,建立了两种ZTD数据源之间的系统差模型。选取IGS测站和陆态网测站,对附加系统偏差改正后的GGOSZTD产品对PPP的收敛速度的影响进行研究。本文研究结果表明:GGOS-ZTD与IGS-ZTD间存在系统偏差,其bias平均为-0.54cm;两者之间较差的RMS平均为1.31cm,说明GGOS-ZTD产品足以满足广大GNSS导航定位用户对对流层延迟改正的需要。将改正了系统差后的GGOS-ZTD产品用于ALBH、DEAR、ISPA测站、PALM测站、ADIS测站、YNMH测站、WUHN测站进行PPP试验,发现可明显提高定位收敛速度,尤其是在U方向上,收敛速度分别提高10.58%、31.68%、15.96%、43.89%、51.46%、14.69%、18.40%。     

An enhanced strategy for GNSS data processing of massive networks

Although the computational burden of global navigation satellite systems (GNSS) data processing is nowadays already a big challenge, especially for huge networks, integrated processing of denser networks with data of multi-GNSS and multi-frequency is desired in the expectation of more accurate and reliable products. Based on the concept of carrier range, in this study, the precise point positioning with integer ambiguity resolution is engaged to obtain the integer ambiguities for converting carrier phases to carrier ranges. With such carrier ranges and pseudo-ranges, rigorous integrated processing is realized computational efficiently for the orbit and clock estimation using massive networks. The strategy is validated in terms of computational efficiency and product quality using data of the IGS network with about 460 stations. The experimental validation shows that the computation time of the new strategy increases gradually with the number of stations. It takes about 14 min for precise orbit and clock determination with 460 stations, while the current strategy needs about 82 min. The overlapping orbit RMS is reduced from 27.6 mm with 100 stations to 24.8 mm using the proposed strategy, and the RMS could be further reduced to 23.2 mm by including all 460 stations. Therefore, the new strategy could be applied to massive networks of multi-GNSS and multi-frequency receivers and possibly to achieve GNSS data products of higher quality.     

京津冀地区GNSS对流层延迟空间插值研究

对流层延迟差异影响合成孔径雷达干涉测量技术(InSAR)形变测量精度;水汽的变化影响天气变化．对流层延迟与水汽具有较好的对应,因此有必要开展全球导航卫星系统(GNSS)对流层延迟的插值研究．以京津冀地区为例,针对GNSS对流层延迟,开展对流层延迟的空间插值研究．首先开展了GNSS对流层延迟与水汽的比较分析,两者存在显著正相关特性,相关性超过91.7％,论证了对流层延迟取代水汽的可行性．然后利用反距离权重法对京津冀地区2016年9月至2017年8月的12组GNSS测站对流层延迟进行空间插值,通过提取插值点对流层延迟与GNSS站点对流层延迟比较验证空间插值精度．全年数据平均偏差最大为1.12 cm,均方根误差最大为0.89 cm;未发生降水过程平均偏差最大为1.25 cm,均方根误差最大为0.82 cm;发生降水过程平均偏差最大为1.08 cm,均方根误差最大为1.38 cm．京津冀平原区域的GNSS对流层延迟空间插值结果精度满足气象等应用要求,可为气象预报和InSAR大气校正提供参考.      

Using external tropospheric corrections to improve GNSS positioning of hot-air balloon

High accurate global navigation satellite systems (GNSS) require to correct a signal delay caused by the troposphere. The delay can be estimated along with other unknowns or introduced from external models. We assess the impact of the recently developed augmentation tropospheric model on real-time kinematic precise point positioning (PPP). The model is based on numerical weather forecast and thus reflects the actual state of weather conditions. Using the G-Nut/Geb software, we processed GNSS and meteorological data collected during the experiment using a hot-air balloon flying up to an altitude of 2000 m. We studied the impacts of random walk noise setting of zenith total delay (ZTD) on estimated parameters and the mutual correlations, the use of external tropospheric corrections, the use of data from a single or dual GNSS constellation and the use of Kalman filter and backward smoothing processing methods. We observed a significant negative correlation of the estimated rover height and ZTD which depends on constraining ZTD estimates. Such correlation caused a degraded performance of both parameters when estimated simultaneously, in particular for a single GNSS constellation. The impact of ZTD constraining reached up to 50-cm differences in the rover height. Introducing external tropospheric corrections improved the PPP solution regarding: (1) shortened convergence, (2) better overall robustness, particularly, in case of degraded satellite geometry, (3) less adjusted parameters with lower correlations. The numerical weather model-driven PPP resulted in 9–12- and 5–6-cm uncertainties in the rover altitude using the Kalman filter and the backward smoothing, respectively. Compared to standard PPP, it indicates better performance by a factor of 1–2 depending on the availability of GNSS constellations, the troposphere constraining and the processing strategy.     

Using a regional numerical weather prediction model for GNSS positioning over Brazil

The global navigation satellite system (GNSS) can provide centimeter positioning accuracy at low costs. However, in order to obtain the desired high accuracy, it is necessary to use high-quality atmospheric models. We focus on the troposphere, which is an important topic of research in Brazil where the tropospheric characteristics are unique, both spatially and temporally. There are dry regions, which lie mainly in the central part of the country. However, the most interesting area for the investigation of tropospheric models is the wet region which is located in the Amazon forest. This region substantially affects the variability of humidity over other regions of Brazil. It provides a large quantity of water vapor through the humidity convergence zone, especially for the southeast region. The interconnection and large fluxes of water vapor can generate serious deficiencies in tropospheric modeling. The CPTEC/INPE (Center for Weather Forecasting and Climate Studies/Brazilian Institute for Space Research) has been providing since July 2012 a numerical weather prediction (NWP) model for South America, known as Eta. It has yield excellent results in weather prediction but has not been used in GNSS positioning. This NWP model was evaluated in precise point positioning (PPP) and network-based positioning. Concerning PPP, the best positioning results were obtained for the station SAGA, located in Amazon region. Using the NWP model, the 3D RMS are less than 10 cm for all 24 h of data, whereas the values reach approximately 60 cm for the Hopfield model. For network-based positioning, the best results were obtained mainly when the tropospheric characteristics are critical, in which case an improvement of up to 7.2 % was obtained in 3D RMS using NWP models.     

Monitoring precipitable water vapor in real-time using global navigation satellite systems

This paper addresses real-time monitoring of the precipitable water vapor (PWV) from GNSS measurements and presents some results obtained from 6-month long GNSS PWV experiments using international and domestic GNSS networks. In the real-time GNSS PWV monitoring system a server/client structure is employed to facilitate formation of PWV networks and single-differenced GNSS measurements are utilized to mitigate errors in GNSS satellites’ orbits and clocks. An issue relating to baseline length between the server and clients is discussed in detail and as a result the PWV monitor is configured to perform in two modes depending on the baseline length. The server estimates sequentially the zenith wet delay of the individual stations, which is then converted into the PWV of the stations. We evaluate system performance by comparing the real-time PWV solution with reference solutions including meteorological measurements obtained with radiosondes and deferred-time precision GNSS PWV solutions. Results showed that the standard deviation of difference between the real-time PWV and the reference solutions ranged from 2.1 to 3.4 mm in PWV for a 6-month long comparison, which was improved to 1.4 to 2.9 mm by reducing comparison period to 20 days in winter.     

基于MGEX站多系统GNSS反演大气可降水量精度评估

2020年6月北斗卫星导航系统(BDS)完成全面组网,为分析其解算水汽信息的精度,选用15个MGEX (Multi-GNSS Experiment)测站2021年10月至11月的观测数据进行水汽反演. 利用GAMIT软件分别解算BDS、GPS、Galileo和GLONASS的观测数据,将得到的对流层天顶延迟(ZTD)与国际GNSS服务(IGS)发布的结果进行对比,并将解算的大气可降水量(PWV)分别与探空数据、ERA5数据计算得到的PWV对比. 实验结果表明:截止高度角设置为5°时,4个卫星系统估计的ZTD均方根 (RMS)均小于13 mm,GPS-PWV、BDS-PWV、Galileo-PWV、GLONASS-PWV与无线电探空可降水量(RS-PWV)相比,RMS平均值分别为2.25 mm、2.46 mm、2.52 mm和2.84 mm,RMS均小于3 mm;与ERA5-PWV相比,RMS平均值分别为1.63 mm、1.86 mm、1.76 mm和1.99 mm,RMS均小于2 mm. GPS探测水汽的精度最高,BDS探测水汽的精度低于GPS和Galileo,高于GLONASS,均满足气象学应用需求.      

一种优化的基于神经网络的经验ZTD模型

目前,经验对流层天顶延迟（ZTD）模型已经有了飞速的发展,因为它们在使用时无需任何测量的实时地面气象数据,这给GNSS用户提供了极大方便。神经网络技术在实测参数型的ZTD建模中已经取得了一定的成果。与此同时,国内虽然有学者构建了神经外网络的经验ZTD模型,其最大的缺点是忽略了ZTD时间变化且只能单独预报ZTD。本文针对这些缺点构建了优化的神经网络经验ZTD模型。试验结果表明,本文提出的神经网络模型可以分别预报天顶干延迟ZHD和天顶湿延迟ZWD,且具有良好的精度：ZHD的Bias和RMSE分别为-3.7和19.8 mm;ZWD的Bias和RMSE分别为-0.6和34.2 mm。本文的神经网络模型预报的ZHD和ZWD的精度均与目前世界著名的GPT2w格网模型相当。另外,与GPT2w模型相比较,神经网络模型最大的优点就是无需庞大的预存格网数据作为输入,在使用时仅需要知道一个训练好的神经网络即可,该特点为GNSS用户提供了极大的方便。     

Climate applications of a global, 2-hourly atmospheric precipitable water dataset derived from IGS tropospheric products

A global, 2-hourly atmospheric precipitable water (PW) dataset is produced from ground-based GPS measurements of zenith tropospheric delay (ZTD) using the International Global Navigation Satellite Systems (GNSS) Service (IGS) tropospheric products (~80–370 stations, 1997–2006) and US SuomiNet product (169 stations, 2003–2006). The climate applications of the GPS PW dataset are highlighted in this study. Firstly, the GPS PW dataset is used as a reference to validate radiosonde and atmospheric reanalysis data. Three types of systematic errors in global radiosonde PW data are quantified based on comparisons with the GPS PW data, including measurement biases for each of the fourteen radiosonde types along with their characteristics, long-term temporal inhomogeneity and diurnal sampling errors of once and twice daily radiosonde data. The comparisons between the GPS PW data and three reanalysis products, namely the NCEP-NCAR (NNR), ECMWF 40-year (ERA-40) and Japanese reanalyses (JRA), show that the elevation difference between the reanalysis grid box and the GPS station is the primary cause of the PW difference. Secondly, the PW diurnal variations are documented using the 2-hourly GPS PW dataset. The PW diurnal cycle has an annual-mean, peak-to-peak amplitude of 0.66, 0.53 and 1.11 mm for the globe, Northern Hemisphere, and Southern Hemisphere, respectively, with the time of the peak ranging from noon to late evening depending on the season and region. Preliminary analyses suggest that the PW diurnal cycle in Europe is poorly represented in the NNR and JRA products. Several recommendations are made for future improvements of IGS products for climate applications.