Helmert Variance Component Estimation for Multi-GNSS Relative Positioning

The Multi-constellation Global Navigation Satellite System (Multi-GNSS) has become the standard implementation of high accuracy positioning and navigation applications. It is well known that the noise of code and phase measurements depend on GNSS constellation. Then, Helmert variance component estimation (HVCE) is usually used to adjust the contributions of di¿erent GNSS constellations by determining their individual variances of unit weight. However, HVCE requires a heavy computation load. In this study, the HVCE posterior weighting was employed to carry out a kinematic relative Multi-GNSS positioning experiment with six short-baselines from day of year (DoY) 171 to 200 in 2019. As a result, the HVCE posterior weighting strategy improved Multi-GNSS positioning accuracy by 20.5%, 15.7% and 13.2% ineast-north-up(ENU) components, compared to an elevation-dependent (ED) priori weighting strategy. We observed that the weight proportion of both code and phase observations for each GNSS constellation were consistent during the entire 30 days, which indicates that the weight proportions of both code and phase observations are stable over a long period of time. It was also found that the quality of a phase observation is almost equivalent in each baseline and GNSS constellation, whereas that of a code observation is different. In order to reduce the time consumption off the HVCE method without sacrificing positioning accuracy, the stable variances of unit weights of both phase and code observations obtained over 30 days were averaged and then frozen as a priori information in the positioning experiment. The result demonstrated similar ENU improvements of 20.0%, 14.1% and 11.1% with respect to the ED method but saving 88% of the computation time of the HCVE strategy. Our study concludes with the observations that the frozen variances of unit weight (FVUW) could be applied to the positioning experiment for the next 30 days, that is, from DoY 201 to 230 in 2019, improving the positioning ENU accuracy of the ED method by 18.1%, 13.2% and 10.6%, indicating the effectiveness of the FVUW.Peer ReviewedPostprint (published version

The Mechanism for GNSS-Based Kinematic Positioning Degradation at High-Latitudes Under the March 2015 Great Storm

In this study, we focus on the kinematic precise point positioning (PPP) solutions at high-latitudes during the March 2015 great geomagnetic storm. We aim to discover the mechanism behind the positioning degradation from the perspective of the impacts of the storm-induced ionospheric disturbance on the global navigation satellite system (GNSS) data processing. We observed that the phase scintillation dominated the amplitude scintillation at high-latitudes and the variation pattern of the rate of total electron content index (ROTI) was consistent with that of the phase scintillation during the storm. The kinematic PPP errors at high-latitudes were almost three times larger than those at the middle- and low-latitude, which were accompanied by large ROTI variations. From the perspective of GNSS data processing, the large positioning errors were also found to be related to the large number of satellites experiencing cycle slips (CSs). Based on the lock time from the ionospheric scintillation monitoring receiver, we found that a large amount of the CSs was falsely detected under the conventional threshold of the CS detector. By increasing such threshold, the kinematic positioning accuracy at high-latitudes can be improved to obtain similar magnitude as at middle- and low-latitude. The improved positioning accuracy may suggest that the ionospheric disturbance induced by the geomagnetic storm at high-latitudes has minor effects on triggering the CSs. Therefore, precise positioning can be achieved at high-latitudes under geomagnetic storms, given that the CS problem is well addressed.The study is funded by the National Natural Science Foundation of China (No.42004012, 42004025), the Natural Science Foundation of Shandong Province, China (No.ZR2020QD048), the State Key Laboratory of GeoInformation Engineering (No.SKLGIE2019-Z-2-2), the State Key Laboratory of Geodesy and Earth's Dynamics (No. SKLGED-2021-3-4) and by the project RTI2018-094295-B-I00 funded by the MCIN/AEI 10.13039/501100011033, which is co-funded by the FEDER programme.Peer ReviewedPostprint (published version

Evaluation of precipitable water vapor from five reanalysis products with ground-based GNSS observations

At present, the global reliability and accuracy of Precipitable Water Vapor (PWV) from different reanalysis products have not been comprehensively evaluated. In this study, PWV values derived by 268 Global Navigation Satellite Systems (GNSS) stations around the world covering the period from 2016 to 2018 are used to evaluate the accuracies of PWV values from five reanalysis products. The temporal and spatial evolution is not taken into account in this analysis, although the temporal and spatial evolution of atmospheric flows is one of the most important information elements available in numerical weather prediction products. The evaluation results present that five reanalysis products with PWV accuracy from high to low are in the order of the fifth generation of European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA5), ERA-Interim, Japanese 55-year Reanalysis (JRA-55), National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR), and NCEP/DOE (Department of Energy) according to root mean square error (RMSE), bias and correlation coefficient. The ERA5 has the smallest RMSE value of 1.84 mm, while NCEP/NCAR and NCEP/DOE have bigger RMSE values of 3.34 mm and 3.51 mm, respectively. The findings demonstrate that ERA5 and two NCEP reanalysis products have the best and worst performance, respectively, among five reanalysis products. The differences in the accuracy of the five reanalysis products are mainly attributed to the differences in the spatial resolution of reanalysis products. There are some large absolute biases greater than 4 mm between GNSS PWV values and the PWV values of five reanalysis products in the southwest of South America and western China due to the limit of terrains and fewer observations. The accuracies of five reanalysis products are compared in different climatic zones. The results indicate that the absolute accuracies of five reanalysis products are highest in the polar regions and lowest in the tropics. Furthermore, the effects of different seasons on the accuracies of five reanalysis products are also analyzed, which indicates that RMSE values of five reanalysis products in summer and in winter are the largest and the smallest in the temperate regions. Evaluation results from five reanalysis products can help us to learn more about the advantages and disadvantages of the five released water vapor products and promote their applications.Peer ReviewedPostprint (published version

Comprehensive analysis reveals CCDC60 as a potential biomarker correlated with prognosis and immune infiltration of head and neck squamous cell carcinoma

BackgroundCoiled-coil domain containing 60 (CCDC60) is a member of the CCDC family, which participates in the progression of many types of cancer. However, the prognostic value of CCDC60 in head and neck squamous cell carcinoma (HNSC) and its function in tumor immunity remain unclear.MethodsCCDC60 expression and its prognostic potential in HNSC were evaluated by bioinformatics approaches, which was validated in human HNSC samples. Genetic alteration analysis of CCDC60 and the underlying biological function of CCDC60 related co-expressed genes in HNSC were analyzed. The impact of CCDC60 on the regulation of immune infiltration in HNSC was comprehensively investigated. In vitro, a series of functional assays on CCDC60 were performed in HNSC cells.ResultsOur study has indicated that compared with the adjacent normal tissues, CCDC60 expression was considerably downregulated in HNSC tissues. High CCDC60 expression was connected with favorable outcome of HNSC patients, and its prognostic significance was examined by distinct clinical characteristics. We identified the CCDC60-related co-expression genes, which were mainly enriched in the NOD-like receptor signaling pathway associated with the inhibition of tumor growth, leading to a better prognosis of HNSC patients. In vitro, CCDC60 overexpression significantly inhibited the growth, migration and invasiveness but regulated cell cycle progression, and promoted cell adhesion of Fadu and Cal27 cells. Additionally, high CCDC60 expression had strong connections with the infiltrating levels of immune cells, immune marker sets, immunomodulators and chemokines in HNSC, suggesting that targeting CCDC60 could be a promising strategy to enhance the efficacy of immunotherapy for HNSC patients.ConclusionTumor suppressor CCDC60 may be identified as a prognostic and immune-related indicator in HNSC, which had the potential functions in regulating the immune infiltration of HNSC and improving the response to immunotherapy for HNSC patients

Middle- and Long-Term UT1-UTC Prediction Based on Constrained Polynomial Curve Fitting, Weighted Least Squares and Autoregressive Combination Model

Universal time (UT1-UTC) is a key component of Earth orientation parameters (EOP), which is important for the study of monitoring the changes in the Earth’s rotation rate, climatic variation, and the characteristics of the Earth. Many existing UT1-UTC prediction models are based on the combination of least squares (LS) and stochastic models such as the Autoregressive (AR) model. However, due to the complex periodic characteristics in the UT1-UTC series, LS fitting produces large residuals and edge distortion, affecting extrapolation accuracy and thus prediction accuracy. In this study, we propose a combined prediction model based on polynomial curve fitting (PCF), weighted least squares (WLS), and AR, namely, the PCF+WLS+AR model. The PCF algorithm is used to obtain accurate extrapolation values, and then the residuals of PCF are predicted by the WLS+AR model. To obtain more accurate extrapolation results, annual and interval constraints are introduced in this work to determine the optimal degree of PCF. Finally, the multiple sets prediction experiments based on the International Earth Rotation and Reference Systems Service (IERS) EOP 14C04 series are carried out. The comparison results indicate that the constrained PCF+WLS+AR model can efficiently and precisely predict the UT1-UTC in the mid and long term. Compared to Bulletin A, the proposed model can improve accuracy by up to 33.2% in mid- and long-term UT1-UTC prediction

Middle- and Long-Term UT1-UTC Prediction Based on Constrained Polynomial Curve Fitting, Weighted Least Squares and Autoregressive Combination Model

Universal time (UT1-UTC) is a key component of Earth orientation parameters (EOP), which is important for the study of monitoring the changes in the Earth’s rotation rate, climatic variation, and the characteristics of the Earth. Many existing UT1-UTC prediction models are based on the combination of least squares (LS) and stochastic models such as the Autoregressive (AR) model. However, due to the complex periodic characteristics in the UT1-UTC series, LS fitting produces large residuals and edge distortion, affecting extrapolation accuracy and thus prediction accuracy. In this study, we propose a combined prediction model based on polynomial curve fitting (PCF), weighted least squares (WLS), and AR, namely, the PCF+WLS+AR model. The PCF algorithm is used to obtain accurate extrapolation values, and then the residuals of PCF are predicted by the WLS+AR model. To obtain more accurate extrapolation results, annual and interval constraints are introduced in this work to determine the optimal degree of PCF. Finally, the multiple sets prediction experiments based on the International Earth Rotation and Reference Systems Service (IERS) EOP 14C04 series are carried out. The comparison results indicate that the constrained PCF+WLS+AR model can efficiently and precisely predict the UT1-UTC in the mid and long term. Compared to Bulletin A, the proposed model can improve accuracy by up to 33.2% in mid- and long-term UT1-UTC prediction

Earth Rotation Parameters Derived from BDS-3 New Signals B1C/B2a Dual-Frequency Combination Observations

The Earth rotation parameters (ERP) play a crucial role in defining the global reference frame and the Global Navigation Satellite System (GNSS) is one of the important tools used to obtain ERP, including polar motion (PM), its rates, and length of day (LOD). The latest IGS Repro3 ERP products, which provided the IGS contribution to the latest ITRF2020, were generated without consideration of the Beidou Navigation Satellite System (BDS) observations. The global BDS, namely the BDS-3 constellation, has been completely constructed from July 2020 and the observing stations are evenly distributed globally now. Two couple dual-frequency combinations, including the B1I/B3I and B1C/B2a combinations, are commonly used for BDS-3 ionosphere-free combination usage. With the goal of identifying the optimal dual-frequency combination for BDS-3 ERP estimates for the future ITRF definition with a consideration of BDS-3, this research evaluated the performance of ERP estimation using B1I/B3I and B1C/B2a combinations. Firstly, we conducted a comparison of the ambiguity resolutions. The mean percentage of successfully resolved ambiguities for the BDS-3 B1C/B2a combination is 86.5%, surpassing that of B1I/B3I. The GNSS satellite orbits and ERP were estimated simultaneously, thus the accuracy of orbits could also reflect the performance of the ERP estimates. Subsequently, we validated the orbits of 22 BDS-3 Medium Earth Orbit (MEO) satellites using Satellite Laser Ranging (SLR), and the root mean square error (RMS) of the SLR residuals for the 3-day arc orbit with B1C/B2a signals was 5.72 cm, indicating superior accuracy compared with the B1I/B3I combination. Finally, we compared the performance of ERP estimation, considering both internal and external accuracy. For the internal accuracy, B1C/B2a-based solutions demonstrated a reduction in mean formal errors of approximately 17% for PM, 22% for LOD, and 21% for PM rates compared with B1I/B3I-based solutions. In terms of external accuracy, we compared BDS-3-derived ERP estimates with the IERS 20C04 products. The B1C/B2a combination exhibited a slightly better standard deviation performance and a significant reduction in mean bias by 56%, 54%, 39%, 64%, and 23% for X, Y polar motion, X, Y polar motion rates, and LOD, respectively, compared with B1I/B3I solutions. In conclusion, the results highlight the excellent signal quality for BDS-3 B1C/B2a and its superiority in ERP estimation when compared with the B1I/B3I combination

Evaluation of Precipitable Water Vapor from Five Reanalysis Products with Ground-Based GNSS Observations

At present, the global reliability and accuracy of Precipitable Water Vapor (PWV) from different reanalysis products have not been comprehensively evaluated. In this study, PWV values derived by 268 Global Navigation Satellite Systems (GNSS) stations around the world covering the period from 2016 to 2018 are used to evaluate the accuracies of PWV values from five reanalysis products. The temporal and spatial evolution is not taken into account in this analysis, although the temporal and spatial evolution of atmospheric flows is one of the most important information elements available in numerical weather prediction products. The evaluation results present that five reanalysis products with PWV accuracy from high to low are in the order of the fifth generation of European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA5), ERA-Interim, Japanese 55-year Reanalysis (JRA-55), National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR), and NCEP/DOE (Department of Energy) according to root mean square error (RMSE), bias and correlation coefficient. The ERA5 has the smallest RMSE value of 1.84 mm, while NCEP/NCAR and NCEP/DOE have bigger RMSE values of 3.34 mm and 3.51 mm, respectively. The findings demonstrate that ERA5 and two NCEP reanalysis products have the best and worst performance, respectively, among five reanalysis products. The differences in the accuracy of the five reanalysis products are mainly attributed to the differences in the spatial resolution of reanalysis products. There are some large absolute biases greater than 4 mm between GNSS PWV values and the PWV values of five reanalysis products in the southwest of South America and western China due to the limit of terrains and fewer observations. The accuracies of five reanalysis products are compared in different climatic zones. The results indicate that the absolute accuracies of five reanalysis products are highest in the polar regions and lowest in the tropics. Furthermore, the effects of different seasons on the accuracies of five reanalysis products are also analyzed, which indicates that RMSE values of five reanalysis products in summer and in winter are the largest and the smallest in the temperate regions. Evaluation results from five reanalysis products can help us to learn more about the advantages and disadvantages of the five released water vapor products and promote their applications