Tropospheric delay is one of the major factors affecting the accuracy of electromagnetic distance measurements. To provide wide-area real-time high precision zenith tropospheric delay (ZTD), the temporal and spatial variations of ZTD with altitude were analyzed on the bases of the latest meteorological reanalysis product (ERA-Interim) provided by the European Center for Medium-Range Weather Forecasts (ECMWF). An inverse scale height model at given locations taking latitude, longitude and day of year as inputs was then developed and used to convert real-time ZTD at GPS stations in Crustal Movement Observation Network of China (CMONOC) from station height to mean sea level (MSL). The real-time ZTD grid product (RtZTD) over China was then generated with a time interval of 5 min. Compared with ZTD estimated in post-processing mode, the bias and error RMS of ZTD at test GPS stations derived from RtZTD are 0.39 and 1.56 cm, which is significantly more accurate than commonly used empirical models. In addition, simulated real-time kinematic Precise Point Positioning (PPP) tests show that using RtZTD could accelerate the BDS-PPP convergence time by up to 32% and 65% in the horizontal and vertical components (set coordinate error thresholds to 0.4 m), respectively. For GPS-PPP, the convergence time using RtZTD can be accelerated by up to 29% in the vertical component (0.2 m).

GPS data during Typhoon Lekima at 700 stations in China were processed by the Precise Point Positioning (PPP) method. A refined regional T_m model was used to derive the precipitable water vapor (PWV) at these GPS stations. Spatio-temporal variations of PWV with the typhoon process were analyzed. As the typhoon approached, PWV at stations near the typhoon center increased sharply from about 50 mm to nearly 80 mm and then dropped back to about 40–50 mm as the typhoon left. Comparisons of GPS, radiosonde, the Global Data Assimilation System (GDAS) Global Forecast System (GFS) analysis products and ERA5 reanalysis products at four matched GPS-RS stations show overall overestimations of PWV from radiosonde, GFS and ERA5 compared with GPS in a statistical perspective. An empirical orthogonal functions (EOF) analysis of the PWV during the typhoon event revealed some different patterns of variability, with both the first EOF (~36.1% of variance) and second EOF (~30.3% of variance) showing distinctively large anomalies over the typhoon landing locations. The typhoon caused a large horizontal tropospheric gradient (HTG) with the magnitude reaching 5 mm and the direction pointing to the typhoon center when it made a landfall on mainland China. The magnitude and the consistency of the HTG direction decreased overall as the typhoon weakened.

The asymmetric delay has a considerable impact on Global Navigation Satellite Systems (GNSS) Positioning, Navigation and Timing (PNT) applications. In GNSS analyses, the impacts of the asymmetric delay are commonly compensated by using the classical methods with considering the north-south and east-west horizontal gradients. In this paper, we have initiatively proposed an extended method where the north-south and east-west horizontal gradients as well as the second-order horizontal gradients are included to better fit the asymmetric delay. The modeling accuracy of the extended method was evaluated at globally distributed 905 GNSS stations during 40 days in 2020. Significant performance of the extended method respect to the classical method was found, where the hydrostatic and wet modeling accuracy at 4° elevation angle was improved from 5.3 and 10.6 mm to 1.6 and 4.9 mm by 70% and 54%, respectively. The GNSS Precise Point Positioning (PPP) performance using the extended method was also validated at 107 Multi-GNSS Experiment (MGEX) stations. The superior performance on the coordinate repeatability and significant effectiveness on the coordinate and Zenith Total Delay (ZTD) estimations were also found, and the maximal vertical (U) coordinate and ZTD difference biases reached 8.6 and −4.5 mm. The extended method is therefore recommended to substitute the classical methods in the GNSS analyses, especially under severe atmospheric conditions.