Kinematic Zenith Tropospheric Delay Estimation with GNSS PPP in Mountainous Areas

The use of global navigation satellite systems (GNSS) precise point positioning (PPP) to estimate zenith tropospheric delay (ZTD) profiles in kinematic vehicular mode in mountainous areas is investigated. Car-mounted multi-constellation GNSS receivers are employed. The Natural Resources Canada Canadian Spatial Reference System PPP (CSRS-PPP) online service that currently processes dual-frequency global positioning system (GPS) and Global’naya Navigatsionnaya Sputnikovaya Sistema (GLONASS) measurements and is now capable of GPS integer ambiguity resolution is used. An offline version that can process the above and Galileo measurements simultaneously, including Galileo integer ambiguity resolution is also tested to evaluate the advantage of three constellations. A multi-day static data set observed under open sky is first tested to determine performance under ideal conditions. Two long road profile tests conducted in kinematic mode are then analyzed to assess the capability of the approach. The challenges of ZTD kinematic profiling are numerous, namely shorter data sets, signal shading due to topography and forests of conifers along roads, and frequent losses of phase lock requiring numerous but not always successful integer ambiguity re-initialization. ZTD profiles are therefore often only available with float ambiguities, reducing system observability. Occasional total interruption of measurement availability results in profile discontinuities. CSRS-PPP outputs separately the zenith hydrostatic or dry delay (ZHD) and water vapour content or zenith wet delay (ZWD). The two delays are analyzed separately, with emphasis on the more unpredictable and highly variable ZWD, especially in mountainous areas. The estimated delays are compared with the Vienna Mapping Function 1 (VMF1), which proves to be highly effective to model the large-scale profile variations in the Canadian Rockies, the main contribution of GNSS PPP being the estimation of higher frequency ZWD components. Of the many conclusions drawn from the field experiments, it is estimated that kinematic profiles are generally determined with accuracy of 10 to 20 mm, depending on the signal harshness of the environment.