The tropospheric delay is a significant error source in Global Navigation Satellite System (GNSS) positioning and navigation. It is usually projected into zenith direction by using a mapping function. It is particularly important to establish a model that can provide stable and accurate Zenith Tropospheric Delay (ZTD). Because of the regional accuracy difference and poor stability of the traditional ZTD models, this paper proposed two methods to refine the Hopfield and Saastamoinen ZTD models. One is by adding annual and semi-annual periodic terms and the other is based on Back-Propagation Artificial Neutral Network (BP-ANN). Using 5-year data from 2011 to 2015 collected at 67 GNSS reference stations in China and its surrounding regions, the four refined models were constructed. The tropospheric products at these GNSS stations were derived from the site-wise Vienna Mapping Function 1 (VMP1). The spatial analysis, temporal analysis, and residual distribution analysis for all the six models were conducted using the data from 2016 to 2017. The results show that the refined models can effectively improve the accuracy compared with the traditional models. For the Hopfield model, the improvement for the Root Mean Square Error (RMSE) and bias reached 24.5/49.7 and 34.0/52.8 mm, respectively. These values became 8.8/26.7 and 14.7/28.8 mm when the Saastamoinen model was refined using the two methods. This exploration is conducive to GNSS navigation and positioning and GNSS meteorology by providing more accurate tropospheric prior information.

The delays of radio signals transmitted by global navigation satellite system (GNSS) satellites and induced by neutral atmosphere, which are usually represented by zenith tropospheric delay (ZTD), are required as critical information both for GNSS positioning and navigation and GNSS meteorology. Establishing a stable and reliable ZTD model is one of the interests in GNSS research. In this study, we proposed a regional ZTD model that makes full use of the ZTD calculated from regional GNSS data and the corresponding ZTD estimated by global pressure and temperature 3 (GPT3) model, adopting the artificial neutral network (ANN) to construct the correlation between ZTD derived from GPT3 and GNSS observations. The experiments in Hong Kong using Satellite Positioning Reference Station Network (SatRet) were conducted and three statistical values, i.e., bias, root mean square error (RMSE), and compound relative error (CRE) were adopted for our comparisons. Numerical results showed that the proposed model outperformed the parameter ZTD model (Saastamoinen model) and the empirical ZTD model (GPT3 model), with an approximately 56%/52% and 52%/37% RMSE improvement in the internal and external accuracy verification, respectively. Moreover, the proposed method effectively improved the systematic deviation of GPT3 model and achieved better ZTD estimation in both rainy and rainless conditions.

The ground-based Global Navigation Satellite System (GNSS) water vapor tomography is increasingly important in GNSS meteorology. As the multi-GNSS and more ground-based GNSS sites can be incorporated into the regional water vapor tomographic model, determining the tomographic window and sampling rate is crucial for the modelling of the water vapor tomography. These two factors affect not only the number of available signal rays from the satellites, but also the number of tomographic voxels crossed by the signal rays. This study uses Hong Kong as the research area to explore the impact of 12 schemes with different tomographic window and sampling rate on the three water vapor tomography methods, including Least squares, Kalman filtering, and Multiplicative Algebraic Reconstruction Technique (MART). Numerical results show that the tomographic results with the three methods get better as the width of the tomographic window decreases and the sampling rate increases in these 12 schemes, and it is found that the Least squares method is most affected by the two factors, followed by Kalman filtering and MART methods. It is recommended to set a tomographic window width of 10 min and a sampling rate of 300 s in a GNSS water vapor tomographic experiment with dense GNSS site like Hong Kong.

The global navigation satellite system (GNSS)-based attitude determination system has attracted more and more attention with the advantages of having simplified algorithms, a low price and errors that do not accumulate over time. However, GNSS signals may have poor quality or lose lock in some epochs with the influence of signal fading and the multipath effect. When the direct attitude determination method is applied, the primary baseline may not be available (ambiguity is not fixed), leading to the inability of attitude determination. With the gradual popularization of low-cost receivers, making full use of spatial redundancy information of multiple antennas brings new ideas to the GNSS-based attitude determination method. In this paper, an attitude angle conversion algorithm, selecting an arbitrary baseline as the primary baseline, is derived. A multi-antenna attitude determination method based on primary baseline switching is proposed, which is performed on a self-designed embedded software and hardware platform. The proposed method can increase the valid epoch proportion and attitude information. In the land vehicle test, reference results output from a high-accuracy integrated navigation system were used to evaluate the accuracy and reliability. The results indicate that the proposed method is correct and feasible. The valid epoch proportion is increased by 16.2%, which can effectively improve the availability of attitude determination. The RMS of the heading, pitch and roll angles are 0.52°, 1.25° and 1.16°.

Due to the great influence of multipath effect, noise, clock and error on pseudorange, the carrier phase double difference equation is widely used in high-precision indoor pseudolite positioning. The initial position is determined mostly by the known point initialization (KPI) method, and then the ambiguities can be fixed with the LAMBDA method. In this paper, a new method without using the KPI to achieve high-precision indoor pseudolite positioning is proposed. The initial coordinates can be quickly obtained to meet the accuracy requirement of the indoor LAMBDA method. The detailed processes of the method follows: Aiming at the low-cost single-frequency pseudolite system, the static differential pseudolite system (DPL) method is used to obtain the low-accuracy positioning coordinates of the rover station quickly. Then, the ambiguity function method (AFM) is used to search for the coordinates in the corresponding epoch. The real coordinates obtained by AFM can meet the initial accuracy requirement of the LAMBDA method, so that the double difference carrier phase ambiguities can be correctly fixed. Following the above steps, high-precision indoor pseudolite positioning can be realized. Several experiments, including static and dynamic tests, are conducted to verify the feasibility of the new method. According to the results of the experiments, the initial coordinates with the accuracy of decimeter level through the DPL can be obtained. For the AFM part, both a one-meter search scope and two-centimeter or four-centimeter search steps are used to ensure the precision at the centimeter level and high search efficiency. After dealing with the problem of multiple peaks caused by the ambiguity cosine function, the coordinate information of the maximum ambiguity function value (AFV) is taken as the initial value of the LAMBDA, and the ambiguities can be fixed quickly. The new method provides accuracies at the centimeter level for dynamic experiments and at the millimeter level for static ones.

In grassland open-pit mining areas, net primary productivity (NPP) is mainly affected by climate conditions and human activities. The identification and assessment of the influence of human activities on NPP is important for mining production and the implementation of ecological restoration. In this study, we explored the influence of human activities on the NPP in the Shengli mining area in Inner Mongolia, China by using the Carnegie–Ames–Stanford Approach (CASA) model and the Chikugo model, in which a calibration method was applied. An analysis of four representative years showed that the proportion of NPP induced by human activities reached 56.2%, that the percentage of pixels with an inhibitory effect on NPP was 99% in 2011 with the highest intensity of mining activity, and that these two values decreased to 11.9% and 69% in 2020, respectively, with the steady implementation of ecological restoration. Moreover, from the analysis of global and local spatial correlation, mining activities and ecological restoration aggravated and weakened the aggregation of NPP induced by human activities, respectively.