15 September 2025, Volume 45 Issue 9
-
Select all|
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 881-887. https://doi.org/10.14075/j.jgg.2025.02.040
Based on the PWV time series data inverted from GPS during 2010-2018 and the PWV data calculated from the ERA5 reanalysis data for the same period, a combination of the generative adversarial networks(GAN) model and the Transformer neural network model in deep learning was employed to achieve short-term forecasting of PWV in Greenland using GPS-PWV data. The prediction results were evaluated using the ERA5 data from 2019. The results show that the model performed well in most areas, with a root mean square error(RMSE) better than 4.5 mm and a correlation coefficient greater than 0.7. The correlation coefficients in spring, autumn, and winter were all above 0.5, while in summer, the correlation coefficients on some dates were slightly lower due to severe weather changes. This method is capable of predicting the spatial distribution and temporal variations of PWV in Greenland.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 888-893. https://doi.org/10.14075/j.jgg.2024.10.465
With the widespread use of GNSS on civil aviation and UAV platforms, scenarios with large elevation differences are increasingly emerging, which puts higher demands on the accuracy of vertical corrections for tropospheric delays. As the most advanced ZTD model currently available, GPT3 has limitations in its vertical corrections for temperature and pressure, making it only suitable for near-surface applications. To address this issue, this paper proposes an improved method of GPT3 by integrating a global temperature lapse rate grid model and a multi-variable atmospheric pressure-altitude formula, constructing the IGPT3 model, and evaluating its adaptability. The results show that compared with UNB3 and GPT3, IGPT3 has improved accuracy by 25.0%, 18.2%, 16.1%, and 72.2%, 50.0%, 69.9% in predicting temperature lapse rate, temperature profile, and pressure profile, respectively. Moreover, the precision improvements in predicting ZHD and ZTD profiles using IGPT3 reached 29.3 mm and 29.0 mm, respectively. Furthermore, the precision of GNSS water vapor inversion at most IGS stations was improved with the use of IGPT3, with the maximum improvement reaching 3.07 mm.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 894-898, 904. https://doi.org/10.14075/j.jgg.2024.09.436
To assess the accuracy of Earth rotation parameters(ERP) derived from the current VLBI global observing system(VGOS) observational data, this paper processes data from VGOS observations made between 2019 and 2023. Simultaneously, traditional S-band and X-band VLBI observational data, which were obtained during the same periods as the VGOS data, were also included in the processing. The EOP 14 C04 series provided by the IERS were used as reference values for accuracy assessment. The results show that, during the selected experimental data periods, the average weighted root mean square(WRMS) of the post-fit residuals of time delay from VGOS observational data is 16.77 ps, which is significantly better than the 30.97 ps from the traditional S/X VLBI measurements during the same periods. In terms of ERP determination, there are still differences between the results from global network station observational data of VGOS and traditional S/X VLBI. Specifically, the ERP polar motion component accuracy derived from VGOS is inferior to that of traditional S/X VLBI, while the UT1-UTC accuracy is comparable to that of traditional S/X VLBI measurements. Further analysis using experimental stations centered on co-located stations indicates that VGOS measurements can provide ERP determination results with higher accuracy.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 899-904. https://doi.org/10.14075/j.jgg.2024.10.487
In order to improve the medium and long term prediction accuracy of polar motion, the LS+Informer(least-squares extrapolation and Informer, LS+Informer) prediction model is proposed based on the traditional LS+AR prediction model combined with deep learning method. The mean absolute error(MAE) was used as the evaluation index based on LS+AR forecast model and Bulletin A forecast results. After LS+Informer method is adopted, the medium and long term prediction accuracy of PMX and PMY components of polar motion can be significantly improved, and the amplitude can reach 61.94% and 64.86% respectively. At the same time, the MAE value can be controlled within 10 mas in the medium and long term 365-day phase. The results show that the prediction accuracy of LS+Informer model is better than that of LS+AR model, and better than that of Bulletin A. The results demonstrate that the LS+Informer model can be effectively used in medium and long term polar motion prediction.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 905-909. https://doi.org/10.14075/j.jgg.2024.09.437
This paper introduces the data preprocessing and interpolation methods for universal time(UT1) and analyzes the factors affecting the interpolation accuracy of UT1, including different interpolation methods and differences in data preprocessing. The study combines international products with domestic co-observation data to analyze the interpolation accuracy of UT1 under different observation frequencies of very long baseline interferometry(VLBI). The results show that among different interpolation methods, cubic spline interpolation is superior. The interpolation process of UT1 must consider the influence of Earth's solid tides, while high-frequency corrections can be neglected. Compared with daily observations, when the observation interval is less than 2 days, the deviation of UT1 can be stabilized at around 0.05 ms. Increasing the observation interval further will lead to a more pronounced increase of UT1 deviation. Verification using domestic observational data shows that with an observation interval of 2 days, the accuracy of UT1 can reach around 0.05 ms.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 910-914. https://doi.org/10.14075/j.jgg.2024.08.409
To address the issues of large pseudorange measurement noise, discontinuous observation data, weak geometric strength in smartphone GNSS observations, and high noise levels in the built-in micro inertial measurement unit(MIMU) of smartphones, a tightly coupled GNSS/MIMU navigation method enhanced by low-Earth orbit(LEO) satellites is proposed. Based on the characteristics of the Walker constellation and the Keplerian orbital parameters of LEO satellites, LEO pseudorange observation data are simulated. The positioning performance of smartphone GNSS single-point positioning, GNSS/MIMU tight integration, and LEO-enhanced GNSS/MIMU tight integration are compared and analyzed. Real smartphone data are used for validation. The experimental results show that after incorporating LEO satellites, the average number of visible satellites increases by about 7. The faster movement speed of LEO satellites further optimizes the GNSS geometric strength, reducing the average PDOP value by 50%. Compared with GNSS single-point positioning, the GNSS/MIMU tight integration improves positioning accuracy in the east, north, and up directions by 18%, 23%, and 10%, respectively. The LEO-enhanced GNSS/MIMU tight integration improves positioning accuracy in the east, north, and up directions by 73%, 81%, and 60% compared with GNSS single-point positioning, and by 67%, 75%, and 56% compared with GNSS/MIMU tight integration.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 915-921. https://doi.org/10.14075/j.jgg.2024.10.476
The presence of differential code bias(DCB) can affect the accuracy of total electron content(TEC) calculations and pseudorange observations, thereby impacting navigation, positioning, timing, and meteorological research results.In order to accurately estimate DCB, this paper analyzes the temporal variations of GPS inter-frequency DCB data provided by the Chinese Academy of Sciences(CAS) from 2021 to 2022 and proposes an LR-KF-LSTM combined model for precise prediction and analysis of DCB. Experimental results indicate that the average absolute percentage error of this method is less than 1.9%, the average absolute error is less than 0.03 ns, and the root mean square error is less than 0.04 ns.Compared with the LSTM model, BP neural network model, and the DCB values from the CAS product, the combined model shows better prediction performance under different solar activity and geomagnetic conditions. This combined network model can effectively predict satellite DCB and also provides a reference for addressing the issue of single-day or multi-day missing data in the DCB data CAS products.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 922-928, 936. https://doi.org/10.14075/j.jgg.2024.10.484
The accuracy of gravity matching navigation is influenced by the adaptability of the sailing sea area. To ensure the effectiveness of matching, it is crucial to conduct path planning for the submersible. In response to the issues of low efficiency, susceptibility to local optimal solutions, and unstable performance of existing algorithms, this paper proposes a hybrid algorithm-particle swarm optimization and grey wolf optimization(PSO-GWO) to enhance the path planning efficiency of submersibles. Building on the traditional grey wolf optimizer(GWO), the proposed algorithm incorporates the memory mechanism of the best position information of particle movement from particle swarm optimization(PSO), reconstructing the update functions for velocity and position, thereby improving the optimization capability and convergence speed of the GWO algorithm. Experiments were conducted using PSO-GWO, GWO, and PSO algorithms on a two-dimensional grid map discretized based on adaptability. The results show that compared with the GWO and PSO algorithms, the PSO-GWO algorithm reduces the planned path length by 10.99% and the number of iterations by 82.43%. This algorithm reduces redundant travel and accelerates convergence, significantly improving the underwater navigation efficiency of the submersible.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 929-936. https://doi.org/10.14075/j.jgg.2024.08.414
We select the hourly data of two components of extensometer at Yichang station from 2013 to 2023, calculate the M2 wave tidal factors of the two components, perform variational mode decomposition(VMD) on them, and compute the power spectral density(PSD) for each decomposition results. Based on this, a Hilbert time-frequency analysis is conducted. The results show that NS component contains information with a frequency of 0.002 7/d, which may be related to the annual variation of water level in the Three Gorges reservoir(with a frequency of 0.002 7/d) and the current geological stress state. The EW component contains information with a frequency of 0.004-0.012/d, which may be related to the rapid changes of water level during the flood season in Three Gorges reservoir(with frequencies of 0.004-0.007/d and 0.010-0.013/d). Combining the VMD decomposition results of tidal factors and earthquake events with magnitude above 4.0 within 200 km of Yichang station, abnormal high value is observed in the detailed information of EW component tidal factor before the Badong M5.1 earthquake in 2013. Additionally, in terms of meteorological factors, the tidal factors do not exhibit a regular annual variation similar to temperature. However, trend changes in air temperature lead to cave temperature changes, which has a certain impact on the variation trend of tidal factors.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 937-943, 949. https://doi.org/10.14075/j.jgg.2024.09.419
The linear or nonlin drift of the spring gravimeter must be addressed for the existing observation system based on spring gravimeters to monitor the long-term change signal of gravity. Simple linear fitting is usually difficult to simulate the real drift of the instrument over a long period of time, and too high-order fitting will also remove the trend change signal of the signal. This study introduces linear and nonlinear drift correction models to more reasonably remove the drift. This study uses the data of the gPhone gravimeter at Zigui station in 2017 and the gPhone gravimeter at Wushi station in 2023 for analysis and verification. First, the correlation between each channel of Zigui station and the original gravity channel is analyzed, and the theoretical gravity effect of the Three Gorges water storage on the Zigui station is calculated with the results that reach more than 100 μGal. Then, the model coefficients are obtained using a variety of regression analysis methods, and the residual results are studied after deducting various correction values and drift terms. The research results show that the changes in the internal pressure and temperature of the instrument are the main causes of gravity data drift. In the regression model used in this paper, the nonlinear regression model(pure quadratic regression) drift removal results can better reflect the gravity changes caused by water storage, and the residual effect obtained after smoothing is closer to theory than the conventional polynomial fitting method for drift removal. Similar conclusions were obtained by verifying the data from Wushi station.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 944-949. https://doi.org/10.14075/j.jgg.2024.09.439
We conduct experiments on the parallel operation of DC and AC geoelectrical resistivity observation systems at Chengdu earthquake monitoring center station by setting different power supply parameters, and summarize the interference characteristics of different geoelectrical resistivity observation systems on geomagnetic records at the same station. The results show that the interference caused by DC geoelectrical resistivity observation on geomagnetic records is in the form of positive and negative large amplitude square waves, and the output current amplitude is positively correlated with the interference amplitude. The interference caused by AC geoelectrical resistivity observation on geomagnetic records is in the form of large amplitude sine waves, and will not interfere geomagnetic records when the output current frequency is 1-10 Hz and the output current amplitude is 0.1-0.5 A. In addition, based on the experimental results of AC geoelectrical resistivity observation, we simulate the amplitude frequency characteristic curve of the FGM-01 fluxgate magnetometer, which conforms to the frequency response characteristics of a low-pass filter.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 950-953, 963. https://doi.org/10.14075/j.jgg.2024.09.453
Based on the earthquakes recorded by Gansu early warning system since its operation and the output results of early warning system, a comparative study is conducted on the output magnitude of earthquake early warning system and the officially cataloged magnitude, and the distribution of magnitude errors is analyzed. We use the back-propagation neural network to learn and predict the warning magnitude parameters, and obtain the warning magnitude. The results indicate that for earthquakes with smaller actual magnitudes, the estimated magnitude is larger, while for earthquakes with larger magnitudes, the estimated magnitude is smaller. The neural network propagation algorithms can simulate and predict warning magnitude, and reduce the deviation between warning magnitude and cataloged magnitude.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 954-963. https://doi.org/10.14075/j.jgg.2024.09.425
An adaptive method for removing power frequency interference is proposed by combining seasonal trend dispersion(STD) decomposition with recursive least squares(RLS). This method utilizes STD decomposition to extract the seasonal component containing power frequency interference from microseismic waves, which serves as the reference signal for RLS, and dynamically updates the algorithm coefficients to make the calculated signal approach the power frequency interference signal. Simulation experiments were designed by superimposing 1 set of microseismic waves without power frequency interference with 3 different types of power frequency interference. The finite impulse response(FIR) filtering, wavelet threshold(WT) filtering, and the method presented in this paper were used for processing. The results indicate that the method in this paper can effectively remove power frequency interference while fully preserving the key time-frequency characteristics of microseismic waves. Furthermore, the application of this method to the microseismic monitoring data processing of Taoyuan coal mine has verified its feasibility and reliability in engineering applications.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 964-971, 982. https://doi.org/10.14075/j.jgg.2025.03.062
Based on the P-wave travel-time data of 2 538 earthquake events(ML≥1.0) recorded by 198 fixed stations of CENC from 2019 to 2023, the three-dimensional velocity structure of the crust in this area is inversed by double-difference tomography, earthquake relocation via double-difference tomography and waveform cross-correlation for enhanced accuracy and reliability, it reveals the coupling relationship between crustal structural heterogeneity and seismic activity under the background of craton destruction. The results show that after relocation, the root-mean-square(RMS) of seismic travel time residuals is reduced from 1.09 s to 0.52 s, and the residual average is optimized from －0.33 s to 0.02 s, the accuracy of the velocity model has been improved. Set four profile lines distributed in the Shanxi Jinzhong fault zone, Yanshan structural zone, Tangshan fault zone, Xingtai-Handan fault zone and other fault zones where earthquake events are more concentrated. Velocity structure analysis reveals that most seismic events are clustered within transitional zones between high- and low-velocity anomalies, the crustal thickness exhibits significant lateral heterogeneity. In the volcanic rift basin of Datong, Shanxi province, enhanced thermal-chemical upwelling corresponds to prominent low-velocity anomalies. The Yanshan tectonic zone displays alternating low-velocity and high-velocity anomalies, reflecting unstable stress conditions and structural anisotropy. Multiple strong earthquakes along the Tangshan fault zone predominantly occur within transitional zones between high-velocity and low-velocity anomalies. The Xingtai-Handan fault zone demonstrates concentrated seismic events likely triggered by lithospheric delamination processes associated with high-velocity anomalies.The dynamic mechanism of the destruction of the North China Craton is jointly dominated by lithospheric delamination, thermochemical erosion, and mantle convective erosion induced by the subduction of the Pacific plate, providing new geophysical constraints for understanding the stability evolution of the craton and the mechanism of strong earthquake propagation.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 972-982. https://doi.org/10.14075/j.jgg.2024.10.472
This study employs time-series InSAR technology to extract the interseismic surface deformation field along the Anninghe-Zemuhe fault zone on the eastern margin of the Sichuan-Yunnan rhombic block. We further invert the interseismic slip rates and locking depths of the fault system. Key findings include: the cross-fault InSAR deformation velocity gradient along the Anninghe fault increases progressively from ~2 mm/a in the northern segment to ~2.5 mm/a in the southern segment, whereas the Zemuhe fault exhibits a relatively smaller gradient of ~1.8 mm/a.Slip rates along the Anninghe fault vary spatially: ~2.8 mm/a at the northern Shimian segment, ~4.0 mm/a at the central Mianning segment, and ~3.2 mm/a at the southern Xichang segment. The Zemuhe fault shows a decreasing slip rate from ~3.1 mm/a in the northern segment to ~2.6 mm/a in the southern Ningnan segment.Locking depth analysis reveals maximum values (~37.8 km) beneath the Xichang segment, gradually shallowing northward to ~5.0 km at Shimian and southward to ~13.2 km at Ningnan.These results demonstrate that the Xichang segment, characterized by higher slip rates and deeper locking depths, exhibits enhanced seismic potential, consistent with its historically dense seismicity. This spatial correlation between geodetic deformation patterns and seismic activity provides critical insights into strain accumulation processes within this active tectonic boundary.
-
Journal of Geodesy and Geodynamics. 2025, 45(9): 983-990. https://doi.org/10.14075/j.jgg.2024.10.471
To investigate the response of L-band ionospheric scintillation to equatorial anomaly characteristics in low-latitude regions, this study utilizes L-band scintillation observation data from three Australian stations(Darwin, Weipa, and Willis Island) combined with global ionospheric VTEC data provided by the Chinese Academy of Sciences(CAS).The study examines the relationship between ionospheric scintillation and VTEC equatorial anomalies under geomagnetically quiet conditions. Results show a significant positive correlation between ionospheric scintillation and equatorial anomaly characteristics at the three stations. On scintillation days, the ionospheric VTEC peak-to-valley difference, peak-to-valley ratio increase, and gradient changes are more pronounced. This paper reveals the response mechanism and investigates the correlation between VTEC and its gradients with scintillation, which is significant for the rapid identification and judgment of scintillation phenomena and for GNSS engineering applications.
