Search

A non-periodic "step-like" variation in displacement is exhibited owing to the repeated instability of expansive soil landslides. The dynamic prediction of deformation for expansive soil landslides has become a challenge in actual engineering for disaster prevention and mitigation. Therefore, a support vector regression prediction (AMPSO-SVR) model based on adaptive mutation particle swarm optimization is proposed, which is suitable for small samples of data. The shallow displacement is decomposed into a trend component and fluctuating component by complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN), and the trend displacement is predicted by cubic polynomial fitting. In this paper, the multiple disaster-inducing factors of expansive landslides and the time hysteresis effect between displacement and its influencing factors are fully considered, and the crucial influencing factors which eliminate the time lag effect and state factors are input into the model to predict the fluctuation displacement. Monitoring data in the Ningming area of China are employed for the model validation. The predicted results are compared with those of the traditional model. The model performance is evaluated through indicators such as the goodness of fit R2 and root mean square error RMSE. The results show that the prediction RMSE of the new model for three monitoring stations can reach 2.6 mm, 6.6 mm, and 2.5 mm, respectively. Compared with the common Grid search support vector regression (GS-SVR), the Particle Swarm Optimization Support Vector Regression (PSO-SVR) and Back Propagation Neural Network (BPNN) models have average improvements of 58.4%, 38.1%, and 25.2% respectively. The goodness of fit R2 is superior to 0.99 in the new method. The proposed model can effectively be deployed for the displacement prediction of non-periodic stepped expansive soil landslides driven by multiple influencing factors, providing a reference idea for the deformation prediction of expansive soil landslides. [ABSTRACT FROM AUTHOR]

For real-time (RT) precise point positioning (PPP), the state space representation (SSR) information is often delayed due to possible communication delays and specific broadcast intervals. In this case, the positioning results will diverge and re-converge due to the increase of SSR products extrapolation errors. In addition, RT orbit and clock offset accuracy, as well as their extrapolation errors, will vary in different systems and satellites. We propose a PPP with ambiguity resolution (PPP-AR) method that combines a time-differenced carrier phase (TDCP) model, in which the characteristics of the orbit and clock are considered. Under normal communication, the PPP-AR solution is obtained by fixing satellites with small SSR product errors. When the communication is abnormal, the TDCP model is utilized to extrapolate user coordinates by considering different extrapolation error characteristics of satellites. The experimental results show that GPS and Galileo SSR products have better accuracy than BDS, with signal-in-space user ranger errors (SISREs) of 2.7, 2.2, and 8.6 cm, respectively. Optimizing the PPP stochastic model based on SISREs can effectively reduce the convergence time. Under 5 min SSR delay, SISREs caused by clock and orbit extrapolation for GPS/Galileo/BDS are 3.5, 1.4, and 2.6 cm, respectively. After optimizing the TDCP stochastic model based on extrapolation errors, the horizontal and vertical positioning accuracies can be maintained at 0.7 cm and 5.0 cm. For multi-GNSS, the combination of the TDCP and PPP-AR can overcome the influence of short delay. After optimizing the stochastic model, the GPS/Galileo/BDS positioning accuracy can be maintained at about 2.4 cm under 3 min delay, showing an accuracy improvement rate of 59.3% compared with the traditional method using only PPP. Additionally, the rapid PPP convergence results can be obtained by inheriting previous filter state information when the communication recovers normally. [ABSTRACT FROM AUTHOR]

Smartphones with the Android operating system can acquire Global Navigation Satellite System (GNSS) raw pseudorange and carrier phase observations, which can provide a new way for the general public to obtain precise position information. However, only postprocessing precise orbit and clock offset products in some older smart devices are applied in current studies. The performances of precise point positioning (PPP) with the smartphone using real-time products and newly smartphones are still unrevealed, which is more valuable for real-time applications. This study investigates the observation data quality and multi-GNSS real-time PPP performance using recent smartphones. Firstly, the observed carrier-to-noise density ratio (C/N0), number of satellites and position dilution of precision (PDOP) of GNSS observations are evaluated. The results demonstrate that the C/N0 received by Huawei Mate40 is better than that of the Huawei P40 for GPS, BDS, QZSS and Galileo systems, while the GLONASS is poorer, and the PDOP of the Huawei P40 is slightly better than that of Mate40. Additionally, a comprehensive analysis of real-time precise orbit and clock offset products performance is conducted. The experiment result expresses that the orbit and clock offset performance of GPS and Galileo is better than that of BDS-3 and GLONASS, and BDS-2 is the worst. Finally, single- and dual-frequency multi-GNSS combined PPP experiments using observations received from smartphones and real-time products are conducted; the results indicate that the real-time static PPP using a smartphone can achieve decimeter-level positioning accuracy, and kinematic PPP can achieve meter-level positioning accuracy after convergence. [ABSTRACT FROM AUTHOR]

Efficient monitoring and early warning are the preconditions of realizing expansive soil landslide hazard prevention and control. Previous early warning of expansive soil landslides was evaluated through soil sampling experiments to analyze the stability coefficient. However, the existing methods lack timeliness and ignore the inconsistent deformation characteristics of different parts of the landslide mass. There are still difficulties in the dynamic numerical early warning of landslides at multiple points. Considering that the degradation of expansive soil landslides' strength is directly reflected by surface displacement, for the Ningming expansive soil demonstration area and based on the GNSS shallow real-time displacement monitoring sequence, a landslide early-warning method based on the GNSS displacement rate combined with the GNSS displacement tangent angle model was proposed, and we thus designed early-warning thresholds for different warning levels. Combined with multi-source data such as soil moisture, soil pressure, and rainfall, the feasibility of accurate early warning of expansive soil landslides based on GNSS real-time surface displacement was verified. The proposed method does not require numerical calculation of internal stress and achieved two successful early warnings of landslides in the test area, which has a certain promotional value. [ABSTRACT FROM AUTHOR]

The orbit maneuver detection is crucial in Global Navigation Satellite System (GNSS) precise orbit determination, which is necessary for adjusting data processing strategies. The frequency of orbit maneuvers for the BeiDou Navigation System is significantly higher than that of other navigation systems, especially for geosynchronous orbit (GEO) and inclined geosynchronous orbit (IGSO) satellites. We propose a novel real-time and postprocessing method for detecting orbit maneuvers for BeiDou satellites based on the orbit differences between the epoch-updated orbit estimated using square root information (SRIF) and the predicted orbit according to the precise orbit estimated during non-maneuver period, as well as the orbital state difference during maneuver and non-maneuver periods. This method has significant advantages over using observation residuals and it is not affected by observation outliers, thus improving the accuracy and timeliness of orbit maneuver detection. We demonstrated that 32 orbit maneuver events of BeiDou satellites were successfully detected in 2022, of which 1 was for medium Earth orbit (MEO), 7 were for IGSO with an average detected maneuvering time of 7–8 min, and 24 were for GEO satellites with an average detected time of 4–5 min. Moreover, our method can be easily integrated into current real-time filter-based precise orbit determination (POD) processing without any extra task line, which simplifies the overall data processing. The data used in this method can be accessed easily, including GNSS observation data, broadcast ephemeris, and other open-source information files. [ABSTRACT FROM AUTHOR]

Real-time satellite clock offset is a crucial element for real-time precise point positioning (RT-PPP). However, the elapsed time for undifferenced (UD) multi-global navigation satellite system (GNSS) real-time satellite clock offset estimation at each epoch is increased with the growth of stations, which may fall short of real-time application requirements. Therefore, a rapid estimation method for UD multi-GNSS real-time satellite clock offset is proposed to improve the computation efficiency, in which both the dimension of the normal equation (NEQ) and the number of redundant observations are calculated before adjustment; if these two values are larger than the predefined thresholds, the elevation mask is gradually increased until they are less than the predefined thresholds. Then, the clock offset estimation is conducted; this method is called clock offset estimation using partial observations. Totals of 50, 60, 70 and 80 stations are applied to perform experiments. Compared to clock offset estimation using all observations, the elapsed times of clock offset estimation using partial observations can be reduced from 6.80 to 3.10 s, 7.93 to 2.97 s, 12.04 to 3.14 s for 60, 70 and 80 stations, respectively. By using the proposed method, the elapsed time of the clock offset estimation at each epoch is less than 5 s. The estimated clock offset accuracy for GPS, BDS-3, Galileo and GLONASS satellites are better than 0.04, 0.05, 0.03 and 0.16 ns when using the partial observations to estimate clock offset with 50, 60, 70 and 80 stations, respectively. For the multi-GNSS kinematic PPP using the estimated clock offset from 50, 60, 70 and 80 stations with partial observations, the positioning accuracy at 95% confidence level in the east, north and up direction are better than 2.70, 2.20 and 5.60 cm, respectively. [ABSTRACT FROM AUTHOR]

The BeiDou Global Satellite Navigation System (BDS-3) broadcast newly developed B1C and B2a signals. To provide a better service for global users, the vertical phase center offset (PCO) and phase center variation (PCV) are estimated for the B1C/B2a ionospheric-free linear combination of the BDS-3 inclined geostationary orbit (IGSO) and medium earth orbit (MEO) satellites in this study. And considering the traditional PCC estimation method needs two Precise orbit determination (POD) processing, based on the correlation between PCO z-offset and PCV, the theoretical analysis and experimental comparison have been made to discuss whether the POD procedure for the PCO estimation can be omitted. The estimated z-offset time series revealed the inadequacy of the solar radiation pressure (SRP) model for the IGSO satellites and the MEO satellites with Pseudo Random Noise code (PRN) C45 and C46. The PCVraws estimated by the traditional method and the PCO estimation omitted method have the same characteristic. The final PCO z-offsets and PCVs calculated by the two schemes agreed very well with differences can be harmlessly ignored, which confirmed that the PCO estimation can be safely omitted to save computation time. The PCC model proposed in this study has been compared with the Test and Assessment Research Center of China Satellite Navigation Office (TARC/CSNO) released model, the qualities of the orbits and BDS-only precise point positioning (PPP) solutions of the new model both show improvements, except for the IGSO orbits. The analysis of the IGSO orbits further verifies the SRP model is not suitable for the IGSO satellites. [ABSTRACT FROM AUTHOR]

Compared with annual crops, perennial crops with longer growing seasons and deeper root systems can fix more sunlight energy, and have advantages in reducing soil erosion and saving water, fertilizer and pesticide inputs. Rice is one of the most important food crops in the world. Perennial rice can be of great significance for protecting the ecological environment and coping with the shortage of young farmers due to urbanization. Oryza longistaminata (OL) is a rhizomatous wild rice with an AA genome and has strong biotic and abiotic resistances. The AA genome makes OL easy to cross with cultivated rice, thus making it an ideal donor material for perennial rice breeding. Sucrose plays an important role in the development and growth of plants. In this study, OL seedlings were cultured in medium with different concentrations of sucrose, and it was found that sucrose of appropriate concentrations can promote the sprout of basal axillary buds and the subsequent development of rhizomes. In order to explore the molecular mechanism, comparative transcriptome analysis was carried out with OL cultured under two concentrations of sucrose, 20 g/L and 100 g/L, respectively. The results showed that the boost of sucrose to rhizome elongation may be due to the glucose and fructose, hydrolyzed from the absorbed sucrose by vacuolar acid invertase. In addition, the consequent increased osmotic pressure of the cells would promote water absorption, which is benefit for the cell elongation, eventually causing the rhizome elongation. These results may provide a reference for elucidating the regulatory mechanism of sucrose on the rhizome development of OL. [ABSTRACT FROM AUTHOR]

Ambiguity resolution (AR) is a core technology that helps to speed up convergence time and increase positioning accuracy for precise point positioning (PPP), and the performance of PPP-AR is based on the quality of ambiguity resolution products. Real-time PPP-AR becomes a reality as users can now obtain publicly accessible real-time observable-specific signal bias (OSB) products from the Centre National d'Etudes Spatiales (CNES). Therefore, an analysis of the quality of OSB products and an evaluation of the performance of PPP-AR are required to promote the application of real-time positioning. For a total of 31 days between day of year (DOY) 121 and 151 in 2021, observation data were collected from 90 stations, and the OSB products were used to assess the experiments. As for the quality of the OSB products, the data availability (DA) of the GPS and Galileo satellites was greater than 97%, whereas that of BDS was less than 60%; the maximum fluctuation value (MAX) and standard deviation (STD) of the GPS, Galileo, and BDS satellites were 0.045 and 0.012; 0.081 and 0.028; and 0.292 and 0.085 cycles, respectively. In terms of ambiguity residuals using the OSB products, the wide-lane (WL) residual percentages within ±0.25 cycles for the GPS, Galileo, BDS-2, and BDS-3 systems were more than 92%, and the narrow-lane (NL) residual percentages within ±0.25 cycles for the four systems were 92%, 89%, 79%, and 60%, respectively. With regard to the performance of PPP-AR, the GPS+Galileo solution showed the best performance in the kinematic positioning mode, in which the mean root mean square (RMS) of positioning accuracy was 1.06, 1.27, and 2.85 cm for the east (E), north (N), and up (U) components, respectively, and the mean convergence time reached 9.6 min. In the static positioning mode, the mean convergence times of the GPS-only and GPS+Galileo solutions were 11.4 min and 8.0 min, respectively, and both of their mean RMS for positioning accuracy reached 0.79, 0.95, and 1.48 cm for the E, N, and U components, respectively. However, the addition of BDS did not further enhance the performance of multi-GNSS PPP-AR in either the kinematic or static positioning mode due to the poor quality of the real-time BDS products. More importantly, a prediction method was proposed to avoid re-convergence and to enhance the reliability of PPP-AR in the event of short-time missing real-time OSB products and to improve the positioning accuracy and the ambiguity fixed rate. [ABSTRACT FROM AUTHOR]

Global Navigation Satellite System (GNSS) positioning technology has become the most effective method for real-time three-dimensional landslide monitoring. However, the GNSS observation signal is easily affected by the presence of a complex landslide environment with high occlusion and strong interference, in which case its accuracy and reliability cannot meet the requirements of landslide deformation monitoring. Although the accelerometers have strong autonomous working capacities and can complement the GNSS landslide monitoring technology, regular GNSS/accelerometer coupled deformation monitoring technology relies on high-quality GNSS measurement information in order to obtain high-precision accelerometer-reckoned results, derived by accurately estimating the baseline shift error (BSE). When the GNSS signal suffers severe interference, the GNSS monitoring error will be partially absorbed by the BSE component after Kalman filtering, resulting in the divergence of the deformation solution. In this study, an abnormal observation variance inflation model was used to process the simulated landslide monitoring data (GNSS and accelerometer raw observation) for three typical scenes—GNSS signal normally locked, signal partially lost, and short-term interruption. The results were as follows: (1) When the GNSS signal was normally locked, the accuracy was comparable to that of the coupled solution employing an accelerometer (the Root Mean Square (RMS) values in the East (E), North (N) and Upward (U) directions were 0.11 cm, 0.33 cm, and 0.30 cm, respectively). (2) When the GNSS signal was partially lost, the accelerometer could effectively suppress the low-precision float solution of the GNSS, but here, the accuracy of the coupled solution would also decrease with the duration of the floats (the RMS values were E—1.21 cm, N—0.31 cm, and U: 0.58—cm, respectively, when the floats lasted for 10 s, and increased to E—3.09 cm, N—0.39 cm, and U—1.14 cm when they lasted for 20 s, wherein E was the main simulated sliding direction). (3) When the GNSS signal was interrupted for a short time, the accuracy of the coupled solution gradually decreased during continuous interruption, and decreased more quickly during the sliding period of the landslide (when the interruption persisted for 10 s, the RMS values in the simulated landslide stability period were E—0.61 cm, N—0.24 cm, and U—0.25 cm, respectively, while in the simulated landslide sliding period they reached E—4.10 cm, N—6.84 cm, and U—2.30 cm). However, raw observations of the accelerometer could still effectively be used to assist in identifying the real state of the landslide, thereby providing auxiliary information pertinent to early landslide disaster warning. [ABSTRACT FROM AUTHOR]

Accurate estimation of tropospheric delay is significant for global navigation satellite system's (GNSS) high-precision navigation and positioning. However, due to the random and contingent changes in weather conditions and water vapor factors, the classical tropospheric delay model cannot accurately reflect changes in tropospheric delay. In recent years, with the development of meteorological observation/detection and numerical weather prediction (NWP) technology, the accuracy and resolution of meteorological reanalysis data have been effectively improved, providing a new solution for the inversion and modeling of regional or global tropospheric delays. Here, we evaluate the consistency and accuracy of three different types of reanalysis data (i.e., ERA5, MERRA2, and CRA40) used to invert the zenith tropospheric delay (ZTD) from 436 international GNSS service (IGS) stations in 2020, based on the integral method. The results show that the ZTD inversion of the three types of reanalysis data was consistent with the IGS ZTD, even in heavy rain conditions. Furthermore, the average precision of the ZTD inversion of the ERA5 reanalysis data was higher, where the mean deviation (bias), mean absolute error (MAE), and root mean square (RMS) were –3.39, 9.69, and 12.55 mm, respectively. The ZTD average precisions of the MERRA2 and CRA40 inversions were comparable, showing slightly worse performance than the ERA5. In addition, we further analyzed the global distribution characteristics of the ZTD errors inverted from the reanalysis data. The results show that ZTD errors inverted from the reanalysis data were highly correlated with station latitude and climate type, and they were mainly concentrated in the tropical climate zone at low latitudes. Compared to dividing error areas by latitude, dividing error areas by climatic category could better reflect the global distribution of errors and would also provide a data reference for the establishment of tropospheric delay models considering climate type. [ABSTRACT FROM AUTHOR]

Since BeiDou Navigation Satellite System (BDS) and Japan's Quasi-Zenith Satellite System (QZSS) have more visible satellites in the Asia-Pacific region, and navigation satellites of Global Positioning System (GPS), Galileo satellite navigation system (Galileo), and GLONASS satellite navigation system (GLONASS) are uniformly distributed globally, the service level of multi-mode Global Navigation Satellite System (GNSS) in the Asia-Pacific region should represent the best service capability. Based on the observation data of 10 Multi-GNSS Experiment (MGEX) stations, broadcast ephemeris and precision ephemeris from 13 to 19 October 2021, this paper comprehensively evaluated the service capability of multi-GNSS in the Asia-Pacific region from three aspects of observation data quality, broadcast ephemeris performance, and precision positioning level. The results show that: (1) the carrier-to-noise-density ratio (C/N0) quality of the GPS and Galileo is the best, followed by BDS and GLONASS, and QZSS is the worst. GPS, BDS-2, GLONASS, and QZSS pseudorange multipath values range from 0 to 0.6 m, while Galileo system and BDS-3 pseudorange multipath values range from 0 to 0.8 m. (2) In terms of broadcast ephemeris accuracy, BDS-3 broadcast ephemeris has the best orbit, and the three-dimensional (3D) Root Mean Square (RMS) is 0.21 m; BDS-2 was the worst, with a 3D RMS of 1.99 m. The broadcast ephemeris orbits of GPS, Galileo, QZSS, and GLONASS have 3D RMS of 0.60 m, 0.62 m, 0.83 m, and 1.27 m, respectively. For broadcast ephemeris clock offset: Galileo has the best performance, 0.61 ns, GLONASS is the worst, standard deviation (STD) is 3.10 ns, GPS, QZSS, BDS-3 and BDS-2 are 0.65 ns, 0.75 ns, and 1.72 ns, respectively. For signal-in-space ranging errors (SISRE), the SISRE results of GPS and Galileo systems are the best, fluctuating in the range of 0 m–2 m, followed by QZSS, BDS-3, Galileo, and BDS-2. (3) GPS, BDS, GLONASS, Galileo, GPS/QZSS, and BDS/QZSS were used for positioning experiments. In static PPP, the convergence time and positioning accuracy of GPS show the best performance. The positioning accuracy of GPS/QZSS and BDS/QZSS is improved compared with that of GPS and BDS. In terms of kinematic PPP, the convergence time and positioning accuracy of GPS/QZSS and BDS/QZSS are improved compared with that of GPS and BDS. In addition to GLONASS and Galileo systems, the other combinations outperformed 3 cm, 3 cm, and 5 cm in the east, north, and up directions. [ABSTRACT FROM AUTHOR]

The orbital maneuvers of Global Navigation Satellite System (GNSS) Constellations will decrease the performance and accuracy of positioning, navigation, and timing (PNT). Because satellites in the Chinese BeiDou Navigation Satellite System (BDS) are in Geostationary Orbit (GEO) and Inclined Geosynchronous Orbit (IGSO), maneuvers occur more frequently. Also, the precise start moment of the BDS satellites' orbit maneuvering cannot be obtained by common users. This paper presented an improved real-time detecting method for BDS satellites' orbit maneuvering and anomalies with higher timeliness and higher accuracy. The main contributions to this improvement are as follows: (1) instead of the previous two-steps method, a new one-step method with higher accuracy is proposed to determine the start moment and the pseudo random noise code (PRN) of the satellite orbit maneuvering in that time; (2) BDS Medium Earth Orbit (MEO) orbital maneuvers are firstly detected according to the proposed selection strategy for the stations; and (3) the classified non-maneuvering anomalies are detected by a new median robust method using the weak anomaly detection factor and the strong anomaly detection factor. The data from the Multi-GNSS Experiment (MGEX) in 2017 was used for experimental analysis. The experimental results and analysis showed that the start moment of orbital maneuvers and the period of non-maneuver anomalies can be determined more accurately in real-time. When orbital maneuvers and anomalies occur, the proposed method improved the data utilization for 91 and 95 min in 2017. [ABSTRACT FROM AUTHOR]

The satellite clocks used in the BeiDou-2 satellite navigation System (BDS) are Chinese self-developed Rb atomic clocks, and their performances and stabilities are worse than GPS and Galileo satellite clocks. Due to special periodic noises and nonlinear system errors existing in the BDS clock offset series, the GPS ultra-rapid clock model, which uses a simple quadratic polynomial plus one periodic is not suitable for BDS. Therefore, an improved prediction model for BDS satellite clocks is proposed in order to enhance the precision of ultra-rapid predicted clock offsets. First, a basic quadratic polynomial model which is fit for the rubidium (Rb) clock is constructed for BDS. Second, the main cyclic terms are detected and identified by the Fast Fourier Transform (FFT) method according to every satellite clock offset series. The detected results show that most BDS clocks have special cyclic terms which are different from the orbit periods. Therefore, two main cyclic terms are added to absorb the periodic effects. Third, after the quadratic polynomial plus two periodic fitting, some evident nonlinear system errors also exist in the model residual, and the Back Propagation (BP) neural network model is chosen to compensate for these nonlinear system errors. The simulation results show that the performance and precision using the improved model are better than that of China iGMAS ultra-rapid prediction (ISU-P) products and the Deutsches GeoForschungsZentrum GFZ BDS ultra-rapid prediction (GBU-P) products. Comparing to ISU-P products, the average improvements using the proposed model in 3 h, 6 h, 12 h and 24 h are 23.1%, 21.3%, 20.2%, and 19.8%, respectively. Meanwhile the accuracy improvements of the proposed model are 9.9%, 13.9%, 17.3%, and 21.2% compared to GBU-P products. In addition, the kinematic Precise Point Positioning (PPP) example using 8 Multi-GNSS Experiment MGEX stations shows that the precision based on the proposed clock model has improved about 16%, 14%, and 38% in the North (N), East (E) and Height (H) components. [ABSTRACT FROM AUTHOR]

The weighted mean temperature (Tm) is a crucial parameter for determining the tropospheric delay in transforming precipitable water vapor. We used the reanalysis data provided by European Centre for Medium-Range Weather Forecasts (ECMWF) to analyze the distribution characteristics of Tm in the vertical direction in China. To address the problem that the precision of the traditional linear function model is limited in fitting the Tm profile, a scheme using the linear and Fourier functions to fit the Tm profile was constructed. Based on the least squares principle (LSQ) to fit the change in its coefficients over time, a Tm model for China with nonlinear elevation correction (CTm-h) was constructed. The experimental results show that, using ECMWF and radiosonde data to evaluate the precision of the CTm-h model, the RMS is 3.43 K and 4.64 K, respectively. Compared to GPT2w, the precision of the CTm-h model in China is increased by about 26.8%. The CTm-h model provides a significant improvement in the correction effect of Tm in the vertical direction, and the Tm profile calculated by the model is closer to the reference value. [ABSTRACT FROM AUTHOR]

Variations in weather conditions have a significant impact on thermokarst lakes, such as the sub-lake permafrost thawing caused by global warming. Based on the analysis of Landsat sensor images by ENVI TM 5.3 software, the present study quantitatively determined the area of the thermokarst lakes and the area of the single selected thermokarst lake in the Beilu River Basin from 2000 to 2016. In an effort to explore the reason for changes in the area of thermokarst lakes, this work used Pearson correlation to analyze the relationship between the area of thermokarst lakes and precipitation, wind speed, average temperature, and relative humidity as obtained from the weather station Wudaoliang. Furthermore, this study used multiple linear regression to comprehensively study the correlation between the meteorological factors and changes in the thermokarst lake area. In this case, the total lake-area changes and the single-area changes exhibited unique patterns. The results showed that the total lake area and the single selected lake area increased year by year. Furthermore, the effects of the four meteorological factors defined above on the total area of typical thermokarst lakes are different from the effects of these factors on the single selected thermokarst lake. While the total area of specific thermokarst lakes exhibited a time lag in their response to the four factors, the surface area of the selected thermokarst lake responded to these factors on time. The dominant meteorological factor contributing to total lake area variations of typical thermokarst lakes is the increasing annual average temperature. The Pearson correlation coefficient between the total area and the annual average temperature is 0.717, suggesting a statistically significant correlation between the two factors. For the selected thermokarst lake, the surface area is related to annual average temperature and wind speed. As a result, wind speed and average temperature could infer the variation law on the thermokarst lake due to the linear fitting equation between area and significant meteorological factors. [ABSTRACT FROM AUTHOR]

The Mu Us dune field is one of China's four major dune fields, which are ecologically vulnerable areas of northwest semiarid land across Shaanxi, Ningxia, and Inner Mongolia, also very sensitive to the global temperature rise and environmental changes. This paper uses data on the temperature, precipitation, and precipitable water vapor (PWV) in the Mu Us dune field and its surrounding areas to analyze and discuss the time series and spatial distribution characteristics of these three factors in this area. The results of the study show that, in recent years, the trend of temperature increase in the Mu Us dune field has been higher than the average level in China, but this trend has gradually subsided since 2000. The spatial distribution of temperature presents an obvious characteristic of gradual increase from north to south and is affected by latitude, altitude, and topography. The annual cumulative precipitation of the Mu Us dune field is lower than the average level in China. However, in recent years, the rate of the increase in precipitation in this area has been significantly higher than that of the average rate of increase in China. The eastern part of the dune field has the most precipitation, which gradually decreases to the west. The spatial distribution of precipitation is greatly affected by monsoon factors in the region and the distribution of rivers. In the research field, PWV has been rising in recent years, which is greatly related to the increase of vegetation coverage in this region. This demonstrates that the Mu Us dune field has experienced a "warmer and wetter" trend in recent years. [ABSTRACT FROM AUTHOR]

The positioning, navigation, and timing (PNT) service of the Global Navigation Satellite System (GNSS) is developing in the direction of real time and high precision. However, there are some problems that restrict the development of real-time and high-precision PNT technology. Satellite orbit maneuvering is one of the factors that reduce the reliability of real-time navigation products, especially the high-frequency orbit maneuvering of geostationary earth orbit (GEO) and inclined geosynchronous orbit (IGSO) satellites. The BeiDou Navigation Satellite System (BDS) constellation is designed to contain GEO, IGSO, and medium earth orbit (MEO). These orbit maneuvers bring certain difficulties for data processing, especially for BeiDou satellites, such as decreased real-time service performance, which results in real-time navigation products including unusable maneuvered satellites. Additionally, the performance of real-time navigation products will decrease because the orbit maneuvers could not be known in advance, which diminishes the real-time PNT service performance of BDS for users. Common users cannot obtain maneuvering times and strategies owing to confidentiality, which can lead to a decline in the BDS real-time service performance. Thus, we propose a method to predict orbit maneuvers. BDS data from the broadcast ephemeris were analyzed to verify the availability of the proposed method. In addition, the results of real-time positioning were analyzed by using ultra-rapid orbit products, demonstrating that their reliability is improved by removing maneuvered satellites in advance. This is vital to improve the reliability of real-time navigation products and BDS service performance. [ABSTRACT FROM AUTHOR]

In this study, an uncombined precise point positioning (PPP) model was established and was used for estimating fractional cycle bias (FCB) products and for achieving ambiguity resolution (AR), using GPS, BDS-2, and Galileo raw observations. The uncombined PPP model is flexible and efficient for positioning services and generating FCB. The FCBs for GPS, BDS-2, and Galileo were estimated using the uncombined PPP model with observations from the Multi-GNSS Experiment (MGEX) stations. The root mean squares (RMSs) of the float ambiguity a posteriori residuals associated with all of the three GNSS constellations, i.e., GPS, BDS-2, and Galileo, are less than 0.1 cycles for both narrow-lane (NL) and wide-lane (WL) combinations. The standard deviation (STD) of the WL combination FCB series is 0.015, 0.013, and 0.006 cycles for GPS, BDS-2, and Galileo, respectively, and the counterpart for the NL combination FCB series is 0.030 and 0.0184 cycles for GPS and Galileo, respectively. For the BDS-2 NL combination FCB series, the STD of the inclined geosynchronous orbit (IGSO) satellites is 0.0156 cycles, while the value for the medium Earth orbit (MEO) satellites is 0.073 cycles. The AR solutions produced by the uncombined multi-GNSS PPP model were evaluated from the positioning biases and the success fixing rate of ambiguity. The experimental results demonstrate that the growth of the amount of available satellites significantly improves the PPP performance. The three-dimensional (3D) positioning accuracies associated with the PPP ambiguity-fixed solutions for the respective only-GPS, GPS/BDS-2, GPS/Galileo, and GPS/BDS-2/Galileo models are 1.34, 1.19, 1.21, and 1.14 cm, respectively, and more than a 30% improvement is achieved when compared to the results related to the ambiguity-float solutions. Additionally, the convergence time based on the GPS/BDS-2/Galileo observations is only 7.5 min for the ambiguity-fixed solutions, and the results exhibit a 53% improvement in comparison to the ambiguity-float solutions. The values of convergence time based on the only-GPS observations are estimated as 22 and 10.5 min for the ambiguity-float and ambiguity-fixed solutions, respectively. Lastly, the success fixing rate of ambiguity is also dramatically raised for the multi-GNSS PPP AR. For example, the percentage is approximately 99% for the GPS/BDS-2/Galileo solution over a 10 min processing period. In addition, the inter-system bias (ISB) between GPS, BDS-2, and Galileo, which is carefully considered in the uncombined multi-GNSS PPP method, is modeled as a white noise process. The differences of the ISB series between BDS-2 and Galileo indicate that the clock datum bias of the satellite clock offset estimation accounts for the variation of the ISB series. [ABSTRACT FROM AUTHOR]

Ambiguity resolution (AR) is critical for achieving a fast, high-precision solution in precise point positioning (PPP). In the standard uncombined PPP (S-UPPP) method, ionosphere-free code biases are superimposed by ambiguity and receiver clock offsets to be estimated. However, besides the time-constant part of the receiver code bias, the complex and time-varying term in receivers destroy the stability of ambiguities and degrade the performance of the UPPP AR. The variation of receiver code bias can be confirmed by the analysis in terms of ionospheric observables, code multipath (MP) of the Melbourne–Wübbena (MW) combination and the ionosphere-free combination. Therefore, the effect of receiver code biases should be rigorously mitigated. We introduce a modified UPPP (M-UPPP) method to reduce the effects of receiver code biases in ambiguities and to decouple the correlation between receiver clock parameters, code biases, and ambiguities parameters. An extra receiver code bias is set to isolate the code biases from ambiguities. The more stable ambiguities without code biases are expected to achieve a higher success rate of ambiguity resolution and a shortened convergence time. The variations of the receiver code biases, which are the unmodeled errors in measurement residuals of the S-UPPP method, can be estimated in the M-UPPP method. The maximum variation of the code biases is up to 16 ns within two-hour data. In the M-UPPP method, the averaged epoch residuals for code and phase measurements recover their zero-mean features. For the ambiguity-fixed solutions in the M-UPPP method, the convergence times are 14 and 43 min with 17.7% and 69.2% improvements compared to that in the S-UPPP method which are 17 and 90 min under the 68% and 95% confidence levels. [ABSTRACT FROM AUTHOR]

The real-time kinematic (RTK) positioning technique of global navigation satellite systems (GNSS) has been widely used for deformation monitoring in the past several decades. The RTK technique can provide relative displacements in a local reference frame defined by a highly stable reference station. However, the traditional RTK solution does not account for reference stations that experience displacement. This presents a challenge for establishing a near real-time GNSS monitoring system, as since the displacement of a reference station can be easily misinterpreted as a sign of rapid movement at the monitoring station. In this study, based on the reference observations in different time domains, asynchronous and synchronous RTK are proposed and applied together to address this issue, providing more reliable displacement information. Using the asynchronously generated time difference of a reference frame, the proposed approach can detect whether a measured displacement has occurred in the reference or the monitoring station in the current epoch. This allows for the separation of reference station movements from monitoring station movements. The results based on both simulated and landslide monitoring data demonstrate that the proposed method can provide reliable displacement determinations, which are critical in deformation monitoring applications, such as the early warning of landslides. [ABSTRACT FROM AUTHOR]

The determination of slant total electron content (STEC) between satellites and receivers is the first step for establishing an ionospheric model. However, the leveling errors, caused by the smoothed ambiguity solutions in the carrier-to-code leveling (CCL) method, degrade the performance of ionosphere modeling and differential code bias (DCB) estimation. To reduce the leveling errors, an uncombined and undifferenced precise point positioning (PPP) method with ambiguity resolution (AR) was used to directly extract the STEC. Firstly, the ionospheric observables were estimated with CCL, PPP float-ambiguity solutions, and PPP fixed-ambiguity solutions, respectively, to analyze the short-term temporal variation of receiver DCB in zero or short baselines. Then, the global ionospheric map (GIM) was modeled using three types of ionospheric observables based on the single-layer model (SLM) assumption. Compared with the CCL method, the slight variations of receiver DCBs can be obviously distinguished using high precise ionospheric observables, with a 58.4% and 71.2% improvement of the standard deviation (STD) for PPP float-ambiguity and fixed-ambiguity solutions, respectively. For ionosphere modeling, the 24.7% and 27.9% improvements for posteriori residuals were achieved for PPP float-ambiguity and fixed-ambiguity solutions, compared to the CCL method. The corresponding improvement for residuals of the vertical total electron contents (VTECs) compared with the Center for Orbit Determination in Europe (CODE) final GIM products in global accuracy was 9.2% and 13.7% for PPP float-ambiguity and fixed-ambiguity solutions, respectively. The results show that the PPP fixed-ambiguity solution is the best one for the GIM product modeling and satellite DCBs estimation. [ABSTRACT FROM AUTHOR]

Satellite antenna phase center offsets (PCOs) and phase variations (PVs) for BeiDou-3 satellites are estimated based on the tracking data of the Multi-GNSS Experiment (MGEX) and the international GNSS Monitoring and Assessment System (iGMAS) network. However, when estimating the (PCOs) of BeiDou-3 medium Earth orbit (MEO) satellites by pure Extending the CODE Orbit Model (ECOM1), the x-offset estimations of the PCOs have a systematic variation of about 0.4 m with the elevation of the Sun above the orbital plane (β-angle). Thus, a priori box-wing solar radiation pressure (SRP) model of BeiDou-3 MEO was assisted with ECOM1. Then, the satellite type-specific PCOs and common PVs were obtained. The estimations of PCOs and PVs were compared with the MGEX PCOs from the precise orbit and clock offset. When the MGEX PCOs were used, the root mean square (RMS) of 24 h overlap was 6.76, 4.36, 1.46 cm, in along-track, cross-track, and radial directions, respectively; the RMS and standard deviations (STD) of the 24 h clock offset overlap were 0.28 and 0.15 ns; the fitting RMS of the 72 h clock offset of the quadratic polynomial was 0.243 ns. After comparing this with the estimated PCOs and PVs, the RMS of the 24 h orbit overlap was decreased by 6.5 mm (10.54%), 1.8 mm (4.4%), and 1.1 mm (8.03%) in the along-track, cross-track, and radial directions, respectively; the RMS and STD of the 24 h clock offset overlap were decreased by 0.024 ns (8.6%) and 0.020 ns (13.1%), respectively; the fitting RMS of the 72 h clock offset of the quadratic polynomial was reduced by about 0.016 ns (6.5%). [ABSTRACT FROM AUTHOR]

In the Global Navigation Satellite System (GNSS) community, the Quasi-Zenith Satellite System (QZSS) is an augmentation system for users in the Asia-Pacific region. However, the characteristics and performance of four QZSS satellite clocks in a long-term scale are unknown at present. However, it is crucial to the positioning, navigation and timing (PNT) services of users, especially in Asia-Pacific region. In this study, the characteristics and performance variation of four QZSS satellite clocks, which including the phase, frequency, frequency drift, fitting residuals, frequency accuracy, periodic terms, frequency stability and short-term clock prediction, are revealed in detail for the first time based on the precise satellite clock offset products of nearly 1000 days. The important contributions are as follows: (1) It is detected that the times of phase and frequency jump are 2.25 and 1.5 for every QZSS satellite clock in one year. The magnitude of the frequency drift is about 10−18. The periodic oscillation of frequency drift of J01 and J02 satellite clocks is found. The clock offset model precision of QZSS is 0.33 ns. (2) The two main periods of QZSS satellite clock are 24 and 12 hours, which is the influence of the satellite orbit; (3) The frequency stability of 100, 1000 and 10,000 s are 1.98 × 10−13, 6.59 × 10−14 and 5.39 × 10−14 for QZSS satellite clock, respectively. The visible "bump" is found at about 400 s for J02 and J03 satellite clocks. The short-term clock prediction accuracy of is 0.12 ns. This study provides a reference for the state monitoring and performance variation of the QZSS satellite clock. [ABSTRACT FROM AUTHOR]

In order to provide better service for the Asia-Pacific region, the BeiDou navigation satellite system (BDS) is designed as a constellation containing medium earth orbit (MEO), geostationary earth orbit (GEO), and inclined geosynchronous orbit (IGSO). However, the multi-orbit configuration brings great challenges for orbit determination. When orbit maneuvering, the orbital elements of the maneuvered satellites from broadcast ephemeris are unusable for several hours, which makes it difficult to estimate the initial orbit in the process of precise orbit determination. In addition, the maneuvered force information is unknown, which brings systematic orbit integral errors. In order to avoid these errors, observation data are removed from the iterative adjustment. For the above reasons, the precise orbit products of maneuvered satellites are missing from IGS (international GNSS (Global Navigation Satellite System) service) and iGMAS (international GNSS monitoring and assessment system). This study proposes a method to determine the precise orbits of maneuvered satellites for BeiDou GEO and IGSO. The initial orbits of maneuvered satellites could be backward forecasted according to the precise orbit products. The systematic errors caused by unmodeled maneuvered force are absorbed by estimated pseudo-stochastic pulses. The proposed method for determining the precise orbits of maneuvered satellites is validated by analyzing data of stations from the Multi-GNSS Experiment (MGEX). The results show that the precise orbits of maneuvered satellites can be estimated correctly when orbit maneuvering, which could supplement the precise products from the analysis centers of IGS and iGMAS. It can significantly improve the integrality and continuity of the precise products and subsequently provide better precise products for users. [ABSTRACT FROM AUTHOR]

The unscented Kalman filter (UKF) can effectively reduce the linearized model error and the dependence on initial coordinate values for indoor pseudolite (PL) positioning unlike the extended Kalman filter (EKF). However, PL observations are prone to various abnormalities because the indoor environment is usually complex. Standard UKF (SUKF) lacks resistance to frequent abnormal observations. This inadequacy brings difficulty in guaranteeing the accuracy and reliability of indoor PL positioning, especially for phase-based high-precision positioning. In this type of positioning, the ambiguity resolution (AR) will be difficult to achieve in the presence of abnormal observations. In this study, a robust UKF (RUKF) and partial AR (PAR) algorithm are introduced and applied in indoor PL positioning. First, the UKF is used for parameter estimation. Then, the anomaly recognition statistics and optimal ambiguity subset of PAR are constructed on the basis of the posterior residuals. The IGGIII scheme is adopted to weaken the influence of abnormal observation, and the PAR strategy is conducted in case of failure of the conventional PL-AR. The superiority of our proposed algorithm is validated using the measured indoor PL data for code-based differential PL (DPL) and phase-based real-time kinematic (RTK) positioning modes. Numerical results indicate that the positioning accuracy of RUKF-based indoor DPL is higher with a decimeter-level improvement compared that of the SUKF, especially in the presence of large gross errors. In terms of high-precision RTK positioning, RUKF can correctly identify centimeter-level anomalous observations and obtain a corresponding positioning accuracy improvement compared with the SUKF. When relatively large gross errors exist, the conventional method cannot easily realize PL-AR. By contrast, the combination of RUKF and the PAR algorithm can achieve PL-AR for the selected ambiguity subset successfully and can improve the positioning accuracy and reliability significantly. In summary, our proposed algorithm has certain resistance ability for abnormal observations. The indoor PL positioning of this algorithm outperforms that of the conventional method. Thus, the algorithm has some practical application value, especially for kinematic positioning. [ABSTRACT FROM AUTHOR]

Real-time kinematic (RTK) positioning is a satellite navigation technique that is widely used to enhance the precision of position data obtained from global navigation satellite systems (GNSS). This technique can reduce or eliminate significant correlation errors via the enhancement of the base station observation data. However, observations received by the base station are often interrupted, delayed, and/or discontinuous, and in the absence of base station observation data the corresponding positioning accuracy of a rover declines rapidly. With the strategies proposed till date, the positioning accuracy can only be maintained at the centimeter-level for a short span of time, no more than three min. To address this, a novel asynchronous RTK method (that addresses asynchronous errors) that can bridge significant gaps in the observations at the base station is proposed. First, satellite clock and orbital errors are eliminated using the products of the final precise ephemeris during post-processing or the ultra-rapid precise ephemeris during real-time processing. Then the tropospheric error is corrected using the Saastamoinen model and the asynchronous ionospheric delay is corrected using the carrier phase measurements from the rover receiver. Finally, a straightforward first-degree polynomial function is used to predict the residual asynchronous error. Experimental results demonstrate that the proposed approach can achieve centimeter-level accuracy for as long as 15 min during interruptions in both real-time and post-processing scenarios, and that the accuracy of the real-time scheme can be maintained for 15 min even when a large systematic error is projected in the U direction. [ABSTRACT FROM AUTHOR]

Due to the cuboid satellite body of BeiDou-3 satellites, the accuracy of their orbit showed a trend of systematic variation with the sun-satellite-earth angle (ε) using the Extend CODE Orbit Model (ECOM1). Therefore, an a priori cuboid box-wing model (named the cuboid model) is necessary to compensate ECOM1. Considering that the body-dimensions and optical properties of the BeiDou-3 satellites used to construct the box-wing model have not yet been fully released, the adjustable box-wing model (ABW) was used for precise orbit determination (POD). The a priori cuboid box-wing model was directly estimated by the precision radiation accelerations, obtained from ABW POD. When using ECOM1 model, for 14 < β < 40°, a linear systematic variation of D0 related to the elevation of the sun above the orbital plane (β-angle) with a slope of 0.048 nm/s2/°, was found for C30. After adding the cuboid model to assist ECOM1 (named Cuboid + ECOM1), the slope was reduced to 0.005 nm/s2/°, and for C20 satellite, the standard deviation (STD) of D0 was improved, from 1.28 to 0.85 nm/s2 (34%). For satellite laser ranging (SLR) validation, when using the ECOM1 model, the systematic variation with the ε angle was about 14 cm for C20 and C30. After using the Cuboid + ECOM1 model, the variation was significantly reduced to about 5 cm. For C20 and C21, compared with the ECOM1 model, the root mean square (RMS) of the ECOM2 and Cuboid + ECOM1 model was improved by about 0.54 (10.3%) and 0.43 cm (8.7%). For C29 and C30, the RMS of ECOM2 and Cuboid + ECOM1 model was improved for about 0.7 (10.9%) and 1.6 cm (25.6%). Finally, the RMS of the SLR residuals of 4.37 to 4.88 cm was achieved for BeiDou-3 POD. [ABSTRACT FROM AUTHOR]

Because there are different types of BeiDou constellations with participating geostationary orbit (GEO) and inclined geosynchronous orbit (IGSO) satellites, the maneuvering frequency of BeiDou satellites is higher than that of other navigation systems. The satellite orbital maneuvers lead to orbital parameter failure for several hours from broadcast ephemeris. Due to the missing initial orbit, the maneuvering thrust, and the period of orbital maneuvering, the orbit products of maneuvering satellites cannot be provided by the International Global Navigation Satellite System (GNSS) Service (IGS) and International GNSS Monitoring and Assessment System (iGMAS). In addition, the period of unhealthy status and the orbital parameters of maneuvering satellites in broadcast ephemeris are unreliable, making the detection of orbital maneuver periods more difficult. Here, we develop a method to detect orbital maneuver periods involving two key steps. The first step is orbit prediction of maneuvering satellites based on precise orbit products. The second step is time period detection of orbit maneuvering. The start time detection factor is calculated by backward prediction orbit and pseudo-range observations, and the end time detection factor is calculated by forward prediction orbit and pseudo-range observations. Data of stations from the Multi-GNSS Experiment (MGEX) and iGMAS were analyzed. The results show that the period of orbit maneuvering could be detected accurately for BeiDou GEO and IGSO satellites. In addition, the orbital maneuver period of other GNSS medium Earth orbit (MEO) satellites could also be determined by this method. The results of period detection for orbit maneuvering provide important reference information for precision orbit and clock offset determination during satellite maneuvers. [ABSTRACT FROM AUTHOR]

The emergence of multiple global navigation satellite systems (multi-GNSS), including global positioning system (GPS), global navigation satellite system (GLONASS), Beidou navigation satellite system (BDS), and Galileo, brings not only great opportunities for real-time precise point positioning (PPP), but also challenges in quality control because of inevitable data anomalies. This research aims at achieving the real-time quality control of the multi-GNSS combined PPP using additional observations with opposite weight. A robust multiple-system combined PPP estimation is developed to simultaneously process observations from all the four GNSS systems as well as single, dual, or triple systems. The experiment indicates that the proposed quality control can effectively eliminate the influence of outliers on the single GPS and the multiple-system combined PPP. The analysis on the positioning accuracy and the convergence time of the proposed robust PPP is conducted based on one week's data from 32 globally distributed stations. The positioning root mean square (RMS) error of the quad-system combined PPP is 1.2 cm, 1.0 cm, and 3.0 cm in the east, north, and upward components, respectively, with the improvements of 62.5%, 63.0%, and 55.2% compared to those of single GPS. The average convergence time of the quad-system combined PPP in the horizontal and vertical components is 12.8 min and 12.2 min, respectively, while it is 26.5 min and 23.7 min when only using single-GPS PPP. The positioning performance of the GPS, GLONASS, and BDS (GRC) combination and the GPS, GLONASS, and Galileo (GRE) combination is comparable to the GPS, GLONASS, BDS and Galileo (GRCE) combination and it is better than that of the GPS, BDS, and Galileo (GCE) combination. Compared to GPS, the improvements of the positioning accuracy of the GPS and GLONASS (GR) combination, the GPS and Galileo (GE) combination, the GPS and BDS (GC) combination in the east component are 53.1%, 43.8%, and 40.6%, respectively, while they are 55.6%, 48.1%, and 40.7% in the north component, and 47.8%, 40.3%, and 34.3% in the upward component. [ABSTRACT FROM AUTHOR]

Global navigation satellite systems (GNSS) have provided an excellent way to monitor micro-deformation in real-time. However, at local sites where landslides frequently occur, the environment can include complex surroundings with mountains, dense vegetation, and human settlements, which can severely degrade the accuracy of positioning with the GNSS technique. In this study, we propose an azimuth-dependent elevation weight (ADEW) model using an azimuth-dependent elevation mask (ADEM) to reduce the effects of multipath errors and improve the accuracy of real-time deformation monitoring in such environments. We developed an adaptive fixed-elevation mask to serve as the outlier of low precision observations at lower elevations for the ADEM, and then, we applied the weighted phase observations into the mitigation process for the effects of multipath errors. The real numerical results indicate that the ADEM model performs better than the conventional model, and the average improvements were 18.91% and 34.93% in the horizontal and vertical direction, respectively. The ADEW model further improved upon the ADEM model results by an additional 21.9% and 29.8% in the horizontal and vertical direction, respectively. Therefore, we propose that the ADEW model can significantly mitigate the effects of multipath errors and improve the accuracy of micro-deformation monitoring via GNSS receivers. [ABSTRACT FROM AUTHOR]

The frequent maneuvering of BeiDou Geostationary Orbit (GEO) and Inclined Geosynchronous Orbit (IGSO) satellites affects the availability of real-time orbit, and decreases the accuracy and performance of positioning, navigation and time (PNT) services. BeiDou satellite maneuver information cannot be obtained by common users. BeiDou broadcast ephemeris is the only indicator of the health status of satellites, which are broadcast on an hourly basis, easily leading to ineffective observations. Sometimes, identification errors of satellite abnormity also appear in the broadcast ephemeris. This study presents a real-time robust detection method for a satellite orbital maneuver with high frequency and high reliability. By using the broadcast ephemeris and pseudo-range observations, the time discrimination factor and the satellite identification factor were defined and used for the real-time detection of start time and the pseudo-random noise code (PRN) of satellites was used for orbital maneuvers. Data from a Multi-GNSS Experiment (MGEX) was collected and analyzed. The results show that the start time and the PRN of the satellite orbital maneuver could be detected accurately in real time. In addition, abnormal start times and satellite abnormities caused by non-maneuver factors also could be detected using the proposed method. The new method not only improves the utilization of observations for users with the data effective for about 92 min, but also promotes the reliability of real-time PNT services., Competing Interests: The authors declare no conflicts of interest.