15 May 2026, Volume 46 Issue 5
-
Select all|
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 522-529. https://doi.org/10.14075/j.jgg.2025.06.210
Based on the MSATSI method with 1°×1° grid resolution, we inverted the present-day tectonic stress field of the Qinghai-Xizang plateau. The results show that: 1) The focal mechanisms are predominantly strike-slip and thrust types, with 72.70% of P-axes and 58.97% of T-axes plunging less than 30°, indicating that the Qinghai-Xizang plateau is dominated by horizontal compression and shear stress. 2) The stress shape ratio R is significantly higher at the plateau margins than in the interior, with 76.21% of maximum principal compressive stress axes (σ1) plunging less than 15° and trending overall NNE-SSW. 3) The maximum horizontal stress (SHmax) of Qinghai-Xizang plateau rotates clockwise around the eastern syntaxis. In the shallow crust (0 to 10 km), SHmax direction aligns with GPS velocity fields and is controlled by active tectonics. SHmax of the northeast margin of Qinghai-Xizang plateau, northeast of the Longmenshan, rotates counterclockwise by 20° to 30° at 10 to 20 km depth, indicating decoupling between shallow and deep stress fields, likely associated with differential deformation mechanisms of upper strike-slip and deep thrust thickening.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 530-540. https://doi.org/10.14075/j.jgg.2025.07.255
We analyze the geochemical characteristics of five hot spring fluids in Eryuan and its surrounding areas, and explore the relationship among geothermal distribution, mantle-derived fluid degassing, hydrochemical anomalies, and seismic activity. The results show that hydrochemical type of hot springs in the study area include HCO3·SO4-Na·Ca, HCO3·SO4-Na, HCO3·Cl-Na, SO4-Ca·Na, and SO4·Cl-Na·Ca. The main sources of hot spring gases are atmospheric saturated water and air, with contributions from volatile components in the crust and mantle. The heat storage temperatures of hot springs in the study area ranges from 164 to 275 ℃, and circulation depths is from 6 369.78 to 11 143.98 meters. The content of He is from 12.90 to 821.62 ppm. Large-scale deep faults in Eryuan area provide favorable channels for the upward migration of deep-sourced fluids, and the deep-sourced fluids are the direct sources of geothermal energy and He. There is a correlation between the geochemical characteristics of hot spring fluids and seismic activities, which is characterized by the deeper the circulation depth, the higher the geothermal reservoir temperature, and the stronger the seismic activity. Moreover, most of the M≥4 earthquakes are located in the transition zone between high and low geothermal reservoir temperatures. The Cl- and SO42- concentrations of Wenmiao hot spring show abnormal increases before and after earthquakes, mainly due to the contribution of the mixing of deep-sourced fluids and different aquifers. The temperature at spring outlet shows an abnormal decrease, which may be related to the fact that multiple fractures provide heat dissipation channels after the earthquake. The hydrochemical values of hot springs were restored after the earthquake, indicating that the deep-sourced fluid activity in Eryuan area weakened.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 541-550. https://doi.org/10.14075/j.jgg.2025.06.212
Based on multi-period ground gravity and geomagnetic observation data from 2016 to 2021 in western Yunnan, this study investigates the 2021 Yangbi M6.4 earthquake using a joint analysis method of gravity and magnetic fields. It analyzes the coupling characteristics of the gravity field and geomagnetic field during the seismogenic process and proposes a new short-term and imminent earthquake prediction method. The synchronously observed gravity and total magnetic intensity variation data were normalized, and the kernel density distribution and correlation coefficient changes of both within the same region were analyzed to study the epicenter location and the imminence of the earthquake occurrence. The results show that: 1) The gravity field changes and total magnetic intensity changes in the epicentral area are coupled, both exhibiting minor variations at the epicenter but significant changes in the surrounding areas. As the time of the earthquake approaches, the trends of both changes become consistent, while they rapidly diverge after the earthquake. 2) For the Yangbi earthquake, within a 1°×1° grid, the region with the largest overlap area of the annual variation kernel density curves of the gravity and magnetic fields is located between 25°-26°N and 99°-100°E, which is highly consistent with the actual epicenter of the Yangbi earthquake. 3) The joint analysis method also demonstrates practical applicability for earthquakes below magnitude 6 in the study area.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 551-560. https://doi.org/10.14075/j.jgg.2025.07.225
Based on the gravity observation data from the Hubei gravity network in 2019, we study the impact of datum points selection on the adjustment accuracy of gravity network from three aspects: the spatial distribution of datum points, the number of datum points, and the maximum gravity segment difference of selected datum points using the classical adjustment method, and obtain the optimal scheme for datum points selection in Hubei gravity network. Additionally, we conduct a forward simulation using priori mean error to evaluate the adjustment effects controlled by different datum points, and explore the optimal selection schemes for single-datum, dual-datum, triple-datum, and quadruple-datum configurations. The results show that: 1) The selection of datum points should prioritize their spatial distribution to maximize the control effect of each datum point. The number of datum points should be considered next, after a certain number is reached, the adjustment accuracy of gravity network will tend to saturate. Finally, the maximum gravity segment difference of selected datum points should be considered, with preference given to datum points with larger segment difference. 2) For the Hubei gravity network with fixed gravity datums, the average adjustment accuracy improvement rates from optimal single-datum point to optimal quintuple-datum points are 16.81%, 8.31%, 4.84%, and 2.98%, respectively. After statistical and simulation analysis, it is suggested that Hubei gravity network use four datum points of Wanzhou, Zigui, Xiangyang, and Wuhan. 3) Using priori mean error to simulate the selection of optimal datum points, the optimal single-datum point is Jingmen, the optimal dual-datum points are Fengjie and Wuhan, the optimal triple-datum points are Wanzhou, Buzhuanghe, and Wuhan, and the optimal quadruple-datum points are Wanzhou, Shiyan, Buzhuanghe, and Wuhan. Compared with the optimal datum points obtained from the actual data analysis, the simulated optimal datum points have a more rational spatial distribution and a slight improvement in adjustment accuracy. The priori mean error method used in this article to simulate the optimal datum points of gravity network has certain reference value for datum points optimization of the existing national network, regional network, or the planning gravity measurements.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 561-568. https://doi.org/10.14075/j.jgg.2025.06.220
To investigate the correlation between precipitable water vapor (PWV) and aerosol optical depth (AOD), a novel AOD estimation model was developed based on PWV. The model incorporates not only a linear relationship between PWV and AOD but also accounts for the periodic characteristics observed in the residuals of the linear fit. Utilizing high-precision PWV and AOD data from 20 stations in Yunnan from 2016 to 2019, the model was constructed and validated. The results reveal a strong lagged correlation between PWV and AOD in the region, with an average correlation coefficient of 0.82. The proposed model achieved excellent performance in AOD estimation, both in internal and external validations, yielding a mean standard deviation (STD) of 0.11 mm and a root mean square error (RMSE) of 0.12 mm. This model effectively paves the way for estimating high-precision AOD from PWV, thereby opening up new avenues for atmospheric environmental monitoring and research using GNSS remote sensing.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 569-578. https://doi.org/10.14075/j.jgg.2025.06.215
An empirical atmospheric weighted average temperature (Tm) model for Qinghai was established based on data from 7 radiosonde stations in Qinghai from 2019 to 2022. This model was then used to retrieve precipitable water vapor (PWV) from 42 GNSS stations in the province for analyzing the spatiotemporal characteristics of water vapor during multiple heavy rainfall processes in August 2022. The results show that: 1) The newly established Tm model has a mean bias of 0.4 K and a root mean square error (RMSE) of 3.41 K, indicating its accuracy is significantly superior to the Bevis model and China-regional models, while being comparable to that of the single-station models. 2) The GNSS PWV retrieved based on the localized Tm model highly agrees with radiosonde results (correlation coefficient of 0.98). 3) Before heavy rainfall occurs, PWV exhibits a rapid increase sustained for several hours, and its 3-hour variation (PWV*) remains positive, whereas before weak rainfall, PWV shows fluctuating changes. 4) Short-term heavy rainfall is likely to occur when the regional PWV value exceeds 34 mm and PWV* is positive; conversely, rainfall intensity is lower.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 579-586. https://doi.org/10.14075/j.jgg.2025.06.223
To investigate the impact of receiver IFCB on GLONASS undifferenced precise time transfer, different IFCB modeling schemes were proposed and their PPP time transfer performance was evaluated, namely the IFCB0 scheme ignoring IFCB, the IFCB1 scheme modeling IFCB as a linear function of frequency numbers, and the IFCB2 scheme modeling IFCB as a quadratic polynomial. The results indicate that IFCB modeling is of great significance for improving the accuracy of PPP time transfer. Compared to the IFCB0 scheme, the IFCB modeling approach demonstrates the most notable improvement in long-term frequency stability for time transfer. The IFCB modeling significantly enhances time transfer stability for hybrid and heterogeneous receiver/antenna types, while its improvement for homogeneous equipment is limited. Additionally, compatibility tests with different satellite precision products show that the IFCB2 scheme aligns best with the precision products from the European Orbit Determination Centre, whereas the IFCB1 scheme is more compatible with the precision products provided by Wuhan University.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 587-598. https://doi.org/10.14075/j.jgg.2025.07.249
Given the issue that multi-system satellite clock bias estimation is susceptible to interference from various noise types, an estimation model that fully considers the characteristics of satellite clock noise components is proposed. This model employs two strategies, namely sliding window and overall modeling, and utilizes the Hadamard variance to solve for noise parameters. Based on this, clock bias estimation and accuracy assessment are conducted. In the experimental section, multiple representative satellites from the four major satellite navigation systems—GPS, BDS-3, Galileo, and GLONASS are selected as research subjects. A systematic comparative analysis is performed on the performance differences of the two modeling strategies across different systems and satellites. The experimental results reveal that the BDS-3 exhibits significant inter-satellite variability. Specifically, most satellites in MEO demonstrate more pronounced time-varying noise characteristics, making them more suitable for the sliding modeling strategy. Conversely, for some satellites equipped with MEO rubidium clocks and stable IGSO satellites, the overall modeling strategy offers greater advantages. In the case of the Galileo, the global modeling strategy performs better. For the GPS and GLONASS, after adopting the sliding window strategy, the positioning accuracy in E, N, U directions improves by 10.34%, 17.07%, and 4.30%, respectively, and by 6.27%, 10.34%, and 3.76%, respectively, compared to the overall modeling strategy.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 599-609. https://doi.org/10.14075/j.jgg.2025.06.219
Accurate estimation of phase biases is an essential prerequisite for achieving PPP ambiguity resolution (AR). We focus on the estimation methods of observable-specific signal bias (OSB) and multi-frequency multi-GNSS PPP-AR technique, and conduct experimental analysis on the combined observations from GPS, Galileo and BDS at 30 MGEX stations during March 2024. The results show that the phase and pseudorange OSB have good stability, the STD of phase OSB in different systems is less than 0.05 ns, and the consistency of different types of inter-frequency phase and pseudorange OSB of each system is good. The wide-lane/narrow-lane ambiguity resolution rates reach 98.5%/96.7% for GPS, 99.3%/98.8% for Galileo, and 95.4%/94.9% for BDS. Compared with float solutions, the convergence times of triple-frequency PPP fixed solutions are reduced by 28.37% in kinematic scenario, shortened by 31.61% in static scenario. Under static condition, average single-day solution positioning accuracy in E, N and U directions improved by 24.91%(0.71 cm), 31.22%(0.50 cm), and 18.77%(1.10 cm), respectively, demonstrating significant enhancement in positioning performance.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 610-617. https://doi.org/10.14075/j.jgg.2025.07.248
Currently, the accuracy of medium- and long-term predictions of Earth rotation parameters (ERP) (prediction duration exceeding 30 days) can only reach the milliarcsecond level. This level of accuracy is significantly lower than the microarcsecond level required by astrometry and the millimeter level needed by space geodesy. This paper systematically evaluated the performance of four deep learning models (BP, RNN, LSTM, and Transformer) in short-, medium-, and long-term ERP forecasts, with a focus on long-term predictions. The experimental results show that in 180-day predictions, the precision of polar motion forecasts by the Transformer model is 43.8% higher than that of traditional methods, and the prediction error of length of day (LOD) is reduced to 0.27 ms. In UT1-TAI predictions, the BP model achieved the best accuracy with an RMSE of 7.83 ms. This study can provide reference for solving the medium- and long-term ERP forecasting challenges in the autonomous operation of the Beidou system.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 618-627. https://doi.org/10.14075/j.jgg.2025.12.462
Conventional SBAS-InSAR methods are limited by their inability to capture the full dimensionality of deformation, making it difficult to effectively represent three-dimensional (3D) surface displacements in mining areas. To address this issue, this study proposes a 3D deformation inversion method that integrates ascending and descending SAR data by combining SBAS-InSAR with Offset-Tracking techniques. Using the Liuhuanggou coal mine as a case study, experiments were conducted to derive the 3D surface deformation field. The inversion results were systematically analyzed in terms of subsidence magnitude, affected area, and spatial distribution. Furthermore, leveling data were employed to compare the accuracy of the proposed 3D deformation inversion method with conventional ascending and descending SBAS-InSAR 2D deformation results. This comparison aimed to evaluate the advantages of the 3D decomposition approach over existing 2D methods. The findings indicate that the 3D inversion method aligns more closely with actual leveling data and achieves higher inversion accuracy than the 2D approach. The proposed method enhances the monitoring precision of SBAS-InSAR to a certain extent and expands its measurement dimensionality, thereby providing a novel and effective technical reference for surface deformation monitoring in mining areas.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 628-636. https://doi.org/10.14075/j.jgg.2025.06.209
Existing landslide research is mostly limited to single deformation data, ignoring the influence of multidimensional environmental characteristics factors on landslide deformation, such as temperature, dew point temperature, precipitation, and spatiotemporal correlation. To address this issue, we propose a multivariate Transformer prediction model based on time series InSAR technology. The model integrates multidimensional environmental features and spatiotemporal correlation data, and uses wavelet transform to decompose deformation sequences into periodic and trend components. Different Transformer structures are constructed to predict each component. The results show that the prediction accuracy of this model is significantly better than that of traditional methods, with RMSE of 0.72 mm and R2 of 0.96. This research provides a new approach for landslide deformation prediction and is of great significance for improving early warning capabilities.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 637-644. https://doi.org/10.14075/j.jgg.2025.07.244
Taking the Yingpan landslide in the Lancang river basin as the study object, this study employed the SBAS-InSAR technique to process 82 scenes of Sentinel-1 descending orbit images acquired from January 2021 to December 2023 for deformation calculation. The spatiotemporal evolution characteristics of its overall deformation were analyzed. Subsequently, the local deformation features along the road and elevation profiles were thoroughly examined, and the response relationship between time-series deformation and rainfall was investigated. Finally, the fast independent component analysis (FastICA) method was applied to decompose the deformation displacement, revealing the spatial distribution differences and temporal variation characteristics of different deformation components of the Yingpan landslide. The results indicate that during the monitoring period, the maximum deformation rate along the line of sight (LOS) direction reached -65 mm/a, with a maximum cumulative deformation of -175 mm. The deformation area extended over 1 000 m in both east-west and north-south directions, and the deformation characteristics varied across different land cover types. A significant correlation was observed between the magnitude of deformation and rainfall, with deformation being markedly greater during the rainy season than the dry season. FastICA analysis revealed that the overall deformation of the landslide exhibits periodicity.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 645-654. https://doi.org/10.14075/j.jgg.2025.07.233
The reliability of SWOT satellite altimetry-based water level data was validated through extraction and analysis, with the July 2024 flood event in the Dongting lake basin serving as a case study to explore its application potential in dynamic flood monitoring. The results demonstrate that the water level data provided by the SWOT satellite exhibits high accuracy, with an optimal correlation coefficient of 0.99 and a RMSE of 0.55 m when compared to in-situ measurements. The derived water extents also show high spatial consistency with optical imagery.
-
Journal of Geodesy and Geodynamics. 2026, 46(5): 655-661. https://doi.org/10.14075/j.jgg.2025.07.229
To improve the interpretation accuracy of aeromagnetic data at different altitudes, this study addresses issues in upward continuation such as the loss of weak signals, excessive suppression of shallow anomalies, and error accumulation from multi-source superposition effects. An integrated compensation method combining height attenuation and local anomaly factors is proposed. Results from theoretical simulations indicate that the error between compensated data and forward-model data is significantly reduced, with the improvement ratio reaching 217.02% at an upward continuation of 5 m (equal to one model height) and still achieving 54.36% at 25 m (five times the model height). For field data, when upward continuing by 10 m from 130 m altitude, the RMSE decreases from 2.955 2 nT to 0.993 2 nT, an improvement of 66.39%. Even at 120 m continuation, the error is further reduced by 2.64%, demonstrating clear enhancement in accuracy. While preserving the structural characteristics of the original data, the compensation method significantly enhances the clarity and spatial resolution of anomaly responses. This provides theoretical support and technical means for refined processing and quantitative interpretation of aeromagnetic data, with promising potential for practical applications.
