In order to improve the accuracy of ZTD model, a fine-scale stratification ZTD (ZTDFS) model in China is proposed. The new model is finely stratified below 10 km and takes into account the seasonal and daily variation characteristics of each layer. At the same time, the model also fully considers the spatiotemporal variation of the vertical lapse rate of third-order exponential function. The radiosonde data shows that the mean RMSE of ZTDFS model is 31.43 mm, which is 36.6% smaller than that of GPT3 model. The GNSS ZTD data shows that above 500 m, the mean RMSE of ZTDFS model is 28.64 mm, which is improved by 9% compared to GPT3. With the increase in altitude, the performance advantage of ZTDFS model gradually emerges.
To address the issue of deviations in amplitude extraction caused by unmodeled errors in global navigation satellite system interferometric reflectometry(GNSS-IR) soil moisture inversion, this study proposes a method that utilizes a semi-parametric model to introduce nonparametric vectors compensating for unmodeled errors, thereby optimizing amplitude extraction. The soil moisture inversion results based on parametric models are compared with those derived from semi-parametric models. The findings demonstrate that the amplitude extracted by the semi-parametric model achieves higher accuracy in both single satellite and multi-satellite fusion scenarios. The proposed approach effectively improves inversion precision by mitigating unmodeled errors through semi-parametric compensation.
This study investigates the issue of using different gravity field model orders for the varying altitudes of highly elliptical orbits(HEOs). Using the gravity field optimization recursive algorithm, 54 circular orbits with different orbital altitudes were integrated with gravity field models of different orders for 6 hours and their errors were statistically analyzed. Four gravity field order selection models, including polynomial, power function, rational number approximation and piecewise selection, were constructed.The four HEOs in the lunar exploration project were simulated for experimental verification. The integration time for the four models decreased from 13-15 seconds to 2-3 seconds, significantly enhancing computational efficiency, with the integration times of the four selection models being close to each other.The polynomial fitting model achieved the highest integration accuracy among the four selection models, with the accuracy of the four HEOs using the polynomial fitting model being 55.25%, 95.88%, 77.33%, and 76.23% higher than that of the piecewise selection model, respectively.
Using the real-time orbit and clock products from CNES over a one-year period, with URE(user range error) as the evaluation metric, this study analyzes the fault rates, fault types, fault sources, as well as the post-fault accuracy and precise point positioning(PPP) performance of GPS, Galileo, BDS-2, and BDS-3 by comparing them with post-processed products. The results indicate that the fault rates of GPS and Galileo are lower than those of BDS, with significant differences observed between the fault rates of BDS-3 and BDS-2. Regarding fault types, abrupt faults are the least common, while clock offsets are the primary source of faults(accounting for over 65%). After fault elimination, the mean URA(user range accuracy) values for GPS, Galileo, and BDS-3 are 0.02 m, 0.03 m, and 0.05 m, respectively, whereas there is considerable variation in URA among different satellites in BDS-2. Global station simulated kinematic PPP tests demonstrate that excluding satellite fault periods can effectively ensure positioning accuracy.
In inertial pedestrian navigation, traditional fixed-threshold detection methods often significantly increase the positioning error of pedestrians due to false positives and false negatives. To address this issue, this paper proposes an improved zero-velocity interval detection algorithm. The algorithm first constructs an adaptive threshold detection model based on the generalized likelihood ratio test(GLRT). It then uses the threshold values obtained from multi-constraint detection to eliminate false positives in the adaptive detection process, thereby improving the detection accuracy of zero-velocity intervals. To verify the performance of the algorithm, two experimental scenarios were designed: constant-speed and fast walking along L-shaped and rectangular routes. The experimental results show that compared with acceleration magnitude variance(AMV) detection and GLRT detection, the proposed algorithm is not only simple to implement but also reduces the root mean square error(RMSE) of reference points by more than 40%, significantly enhancing the accuracy of inertial pedestrian navigation.
A combined positioning method that integrates ultra wide band (UWB) with low-cost micro-electro-mechanical systems (MEMS) inertial measurement unit (IMU) is proposed. Initially, UWB ranging data is optimized using ranging bias correction and Kalman filtering. Subsequently, the complementary strengths of multiple sensors are leveraged through a UWB/IMU fusion model. Experimental results demonstrate that the optimized UWB achieves an improvement in average positioning accuracy from 15.70 cm to 13.75 cm, with the RMSE decreasing from 22.45 cm to 19.70 cm. In non-line-of-sight (NLOS) scenarios, even when UWB signals experience abrupt changes or are lost, the fused system can still provide continuous and reliable positional information, with an average accuracy of 13.23 cm and an RMSE of 17.90 cm.
Existing research on GNSS signal classification primarily targets professional-grade receivers, distinguishing between line-of-sight(LOS) and non-line-of-sight(NLOS) signals through the extraction of signal features to achieve high-accuracy positioning. However, the signal features and classification methods used in these studies may not be applicable to the signal characteristics of low-cost GNSS chips in smartphones. This paper proposes a machine learning-based method for classifying smartphone GNSS reception signals into LOS and NLOS. The method combines five key indicators: carrier-to-noise ratio, elevation angle, pseudorange residuals, normalized pseudorange residuals, and pseudorange residual percentage. It also integrates deep learning-based zenith view semantic segmentation and satellite calibration techniques to construct an accurate training dataset. Evaluations of various machine learning models in typical urban scenarios show that the proposed method achieves an average classification accuracy of 85.8% across different smartphone models.
To simulate the impact of terrain on co-seismic deformation, this paper established three different terrain models. The results show that compared with flat terrain, convex terrain reduces the amount of horizontal displacement, while concave terrain significantly increases the amount of horizontal displacement. The terrain difference in the epicentral area of the 2024 Wushi, Xinjiang MW7.0 earthquake is significant, and the impact of the real terrain on the co-seismic deformation of this earthquake is not yet clear. This paper uses the spectral element method to simulate the co-seismic deformation of the Wushi earthquake and compares it with the co-seismic displacement field of D-InSAR. The results indicate that the impact of terrain factors on co-seismic deformation can reach 8.5%, 26.9%, and 6.8% in the E, N, and U directions, respectively.
We analyzed the Dingri MS6.8 earthquake as a case study. A multivariate dataset comprising land surface temperature, net surface solar radiation, and other variables in the study area was constructed. After decomposing the time series of each variable, the generalized extreme studentized deviate(GESD) algorithm was applied to extract anomaly signals. The weighted averaging method was then employed to integrate multi-variable anomaly extraction results. The results indicate that the thermal infrared anomalies in the half-year period before the Dingri MS6.8 earthquake exhibited a cyclic change pattern of "emergence-enhancement-weakening". The Dingri MS6.8 earthquake occurred 100 days after the maximum value of the anomalies appeared, and the main direction of the anomalies was NW-SE, forming an approximately 45° angle with the Dengmo Co fault, which is the seismogenic structure of this earthquake.
Using Sentinel-1A data, this study employed D-InSAR technology to obtain the line-of-sight (LOS) co-seismic deformation field of the MS6.8 earthquake in Dingri, Xizang, on January 7, 2025. The geometric parameters of the seismogenic fault were derived through a Bayesian approach. Fault slip distribution was inverted using the SDM, and co-seismic Coulomb stress changes were calculated. The inversion results indicate that the seismogenic fault has a strike of 187.34°, a dip angle of 56.40°, with an epicenter located at 87.50°E and 28.64°N, at a depth of 6.81 km. The maximum slip in the co-seismic slip distribution is 4.24 m, with an average rake angle of-58.82°, suggesting that the earthquake was predominantly a normal faulting event with a minor component of left-lateral strike-slip. The inverted moment magnitude is MW7.0, and the seismogenic fault is preliminarily identified as the Dengmo Co fault. Additionally, regions along the Zanda-Lhazê-Qiongduojiang fault, Darjeeling-Ngamring-Rinbung fault, and parts of the south xizang detachment system fault, as well as the southern sections of the Xainza-Dingjie fault and Dengmo Co fault, exhibit Coulomb stress increases exceeding 0.1 bar, indicating that these areas warrant greater attention for future seismic hazards.
In this paper, the geometry, kinematics and chronology of Shenjia-Guanzhuang ductile shear zone on the eastern margin of Dabie orogenic belt are studied in detail by means of field investigation, laboratory study and zircon U-Pb dating, and the geological significance of shear activity duration is discussed. The results showed that: The Shenjia-Guanzhuang ductile shear zone changes from NNW-SSE to SWW-NEE from west to east, and its tendency rotates counterclockwise from NEE to NNW, and its inclination gradually increases from about 40° to about 55°. Its kinematic characteristics change from forward sliding overlying shear to right lateral translacing shear. The overall performance is a "bow" shape with SW protrating, and the upper plate moves towards NEE. The start-up time is not earlier than 129.3±2.5 Ma, and the end is before 124.7±1.7 Ma. The end of ductile shear activity marks the end of intracontinental orogenic evolution of Dabie orogenic belt, and the beginning of post-orogenic uneven block uplift stage. The tectonic regime transformation of Dabie mountains and even eastern China has been completed, and the Dabie mountains have entered the coastal Pacific tectonic domain, which may be closely related to the major adjustment of the movement mode of the Izanagi plate in the western Pacific.
By comparing the differences between various drift models and adjustment methods through simulated gravity data adjustment, certain optimizations were made to the gravity adjustment algorithm. Methods for multiple segment difference adjustment and constrained linear programming adjustment were proposed. The results indicate that the multiple segment adjustment method outperforms other adjustment methods in some aspects, while the constrained linear programming adjustment method demonstrates overall superiority over other adjustment methods.
A multipath error mitigation method suitable for dynamic landslide deformation monitoring points is proposed. This method replaces the fixed position parameters used in the residual extraction process with dynamic displacement parameters that reflect actual deformation to eliminate systematic errors caused by deformation. The method employs sidereal filtering(SF) and multipath hemispherical mapping(MHM) to mitigate multipath effects in dynamic landslide scenarios. Experimental results show that the proposed multipath mitigation method based on dynamic displacement parameters can be effectively applied to dynamic deformation monitoring scenarios. Compared with the original unfiltered results, during the constant velocity deformation stage of the landslide, the positioning accuracy in the east direction is improved by 13% and 27%, in the north direction by 15% and 30%, and in the height direction by 11% and 22% for SF and MHM based on dynamic displacement parameters. During the accelerated deformation stage of the landslide, the positioning accuracy in the east direction is improved by 12% and 24%, in the north direction by 23% and 32%, and in the height direction by 5% and 19% for SF and MHM based on dynamic displacement parameters.
To address the issues of abnormal and non-smooth fitting results in sparse monitoring station areas and block margins when using the traditional least squares collocation (TLSC) algorithm for high-precision fitting of crustal motion velocity fields in large-scale regions, a K-means clustering-based least squares collocation method (KLSC) is developed by integrating the K-means clustering algorithm with the TLSC algorithm. The effectiveness of this method is validated using the GNSS-measured crustal motion velocity field on the Qinghai-Xizang plateau. The results show that: 1) Compared to the TLSC algorithm, the KLSC algorithm leverages the advantages of the K-means algorithm in unsupervised classification to first divide the study area into multiple sub-regions with similar velocities based on the inherent characteristics of the GNSS velocity field, and then applies TLSC for velocity field fitting within each sub-region, thereby avoiding the impact of local complex geological environments on the accuracy of regional velocity field fitting; 2) The KLSC algorithm selects fitting parameters based on the proximity of each grid point to the cluster centers, resolving the issue of poor fitting results in data-sparse areas; 3) By utilizing secondary nearest-neighbor fitting combined with convolution filtering, the KLSC algorithm effectively addresses the non-smooth fitting results at block margins; 4) The RMSE accuracy and correlation of the velocity fields fitted by the KLSC algorithm are superior to those of the TLSC algorithm, with improvements in RMSE accuracy of 37% to 48.2% and 52.1% to 67.2%, and increases in correlation of 24.1% to 24.7% and 4.7% to 5.2% for the eastward and northward fitted velocity fields, respectively.
We utilize long-term observation data from five superconductor gravimeters in China and Europe within the IGETS network, and apply various atmospheric pressure correction strategies to analyze the impact of time-varying admittance on the determination of Earth tide parameters. The results show that the temporal variability of admittance significantly influences Earth tide parameters, with effects on time-varying characteristics of the Earth's free core nutation period reaching up to 10 sidereal days. Taking into account the time-varying characteristics of admittance and global atmospheric pressure effects correction, the tidal wave amplitude factor is found to be closer to theoretical values. This has important implications for the construction of regional tidal models and subtle dynamic signals analysis.
We quantitatively analyzed the crustal strain changes in Xuzhou and Gaochun areas of Jiangsu province from 2019 to 2023 using four-component borehole strainmeters, and explored the characteristics of regional strain field and its relationship with seismic activity. According to the strain parameter results, both Xuzhou and Gaochun areas were in a stable state during the study period. The response values of tidal factors at Gaochun seismic station exhibited an abnormal change of first dispersion, then stabilization before and after the M5.0 earthquake in Dafeng sea area of Yancheng, Jiangsu province. Through the analysis of its variation characteristics, the relationship with the strike of faults and the focal mechanism solution of earthquakes, it can be inferred that in the process of earthquake preparation, the stress and strain in this area experienced a process from accumulation to release, and the tidal factors can reflect this change to a certain extent.