By employing fundamental methods from the theory of curves and surfaces, the concepts of different types of latitude and longitude (reduced latitude and longitude, geometric latitude and longitude, and geodetic latitude and longitude) at an arbitrary point on a triaxial ellipsoidal surface are proposed. Theoretical relationship models are established between Cartesian coordinates (spatial rectangular coordinates) and geodetic coordinates (geodetic latitude, longitude, and height), as well as between Cartesian coordinates and geometric coordinates (geometric latitude, longitude, and height), for any point in space both inside and outside the ellipsoid. Based on the theoretical relationships among the different latitude and longitude types, approximate geodetic coordinates of the study point are derived. Using these approximations as initial values, a novel method for the inverse transformation of coordinates is presented via Newton's iterative approach. Extensive numerical calculations demonstrate that, for any arbitrary point in space, the proposed method accomplishes the transformation from Cartesian to geodetic coordinates nearly instantaneously; even for points near the Earth's surface, only three iterations are required to achieve convergence.

To address the issues of reduced accuracy and misidentification of single water bodies caused by interpolation in large-scale flood inundation mapping using spaceborne GNSS-R technology, this paper proposes a dynamic monitoring method based on the reflectivity change rate using period-matched observation-based interpolation (POBI). Utilizing CYGNSS reflectivity data during the rainy seasons from 2019 to 2023 and China's 1 km resolution monthly precipitation products, we constructed a POBI spatial interpolation model. By combining this with the bistatic radar equation to retrieve surface reflectivity, we systematically analyzed the flood evolution process during the extreme rainfall event in Guilin in June 2024. The results show that, for reflectivity interpolation in unsampled areas, the root mean square error (RMSE) decreased by an average of 20.98% compared to the natural neighbor method. The inundation identification mechanism based on the change rate effectively avoids the subjectivity of traditional threshold methods. The monitoring results indicate that the inundated area in the central urban district expanded by 48.68% from June 19 to June 12, with spatial evolution highly consistent with changes in SMAP soil moisture. This study provides reliable methodological support for refined flood monitoring in highly dynamic, low-coverage areas using spaceborne GNSS-R.

A high-precision seismic dataset was constructed using deep learning methods, and a regional three-dimensional velocity model was obtained through double-difference tomography. The study reveals a correlation between surface velocity structures and geological features: mountainous areas generally exhibit high-velocity P-wave anomalies, whereas depression zones are characterized by low-velocity P-wave anomalies. This pattern may be attributed to the widespread presence of bedrock with relatively high wave velocity in folded regions, in contrast to the Quaternary sedimentary layers with lower velocity in depression areas. Seismic epicenters are densely distributed along the transition zones between high- and low-velocity anomalies. These transitional regions experience significant stress gradients, which facilitate rock fracturing and lead to heightened seismic activity.The Luxi uplift is characterized by several near-vertical high- and low-velocity blocks, with some earthquake epicenters also displaying steeply dipping distributions, suggesting the presence of high-angle faults. This structural configuration may result from multiple phases of compression and extension under the influence of Pacific Plate subduction and the tectonic activity of the Tan-Lu fault zone. These processes have led to cyclic uplift and subsidence in the region, during which progressive compression and extension of rock folds gradually increased their dip angles, forming near-vertical velocity blocks and high-angle faults.Integrated analysis of P- and S-wave velocities suggests the presence of mantle-derived basaltic magmatic intrusions beneath both the Nishan area and the northeastern region of the Lushan. In the Nishan area, the intrusions are predominantly plutonic, whereas in northeastern Lushan, mantle-derived basaltic magma likely reached the surface, forming ancient volcanoes and extrusive rocks.

Based on the methods of low Earth orbit (LEO) satellite orbit determination and inter-satellite single-difference ambiguity resolution, this paper utilizes uncalibrated phase delay (UPD) and observable-specific signal bias (OSB) products to achieve precise orbit determination for eight LEO satellites, namely GRACE-C/D, SWARM-A/B/C, and SENTINEL-3A/3B/6A, and evaluates the impact of these two products on orbit determination accuracy. The results indicate that both products achieve a nearly 100% fixed rate for wide-lane ambiguities; except for the GRACE satellites, which have a narrow-lane fixed rate of approximately 85%, the remaining satellites achieve rates close to 95%. Ambiguity resolution significantly enhances orbit determination accuracy, with the fixed solutions obtained using UPD and OSB products showing comparable three-dimensional accuracy, reaching 1.5 cm, 1.7 cm, and 0.9 cm for the GRACE-FO, SWARM, and SENTINEL series satellites, respectively, representing improvements of up to 50% compared to the float solutions. Additionally, a multi-system analysis of the SENTINEL-6A satellite reveals that the fixed solution accuracy for the combined GPS/Galileo orbit determination is 0.9 cm, representing improvements of 10% and 18% compared to single-GPS and single-Galileo systems, respectively. This study confirms that ambiguity resolution can effectively improve the orbit determination accuracy of LEO satellites and provides a reference for future research on orbit determination for large-scale LEO constellations.

Far-field earthquake triggering mechanisms represent a critical frontier in earthquake prediction research. To explore the spatiotemporal relationship between strong seismic activity in Japan and moderate-to-strong earthquakes in northeast China, this study develops a data-driven prediction model. Based on the USGS earthquake catalog from 1980 to 2024, and using M≥6.0 earthquakes in Japan as potential triggers, the model predicts the probability of earthquakes with M≥4 occurring in northeast China within the next 60 days. A daily-resolution time series dataset is constructed, incorporating 32-dimensional features covering earthquake statistics, spatial distribution, energy release, and aftershock sequences. Deep learning models including long short-term memory (LSTM) neural networks, attention mechanism-enhanced LSTM (Attention-LSTM), and Transformer are systematically compared against traditional machine learning methods such as logistic regression and random forest. Results indicate that the Attention-LSTM model performs optimally, achieving an F1 score of 0.915 and an AUC value of 0.661, and shows significant advantages in Molchan error diagram analysis, a method specific to earthquake prediction evaluation. The model enables both binary classification prediction and generation of spatial probability distribution maps at 1°×1° grid resolution. This study demonstrates the potential of deep learning for revealing cross-regional far-field earthquake triggering mechanisms and provides new insights for short-to-medium-term regional seismic hazard assessment.

To systematically evaluate the accuracy performance of the next-generation VLBI global observing system (VGOS) in determining UT1-UTC, we process synchronous observation data from the 2024 VGOS-INT-A and traditional S/X-band IVS-INT-1 sessions using identical strategies, with the international IERS C04 series serving as the reference for accuracy assessment. The results indicate that, under identical station counts and geometric configurations, the number of observations per session for VGOS-INT-A (approximately 65) is significantly greater than that for IVS-INT-1 (approximately 30). The weighted root mean square (WRMS) of the post-fit delay residuals for both systems is comparable. Regarding UT1-UTC solution accuracy, VGOS-INT-A demonstrates a clear advantage over IVS-INT-1, the RMS of its deviation with the C04 series is ±0.046 ms, superior to the ±0.060 ms of the latter. Its mean formal error is 14.26 μs, while that of IVS-INT-1 is 30.38 μs. Furthermore, preliminary analysis of data from China's VGOS stations shows good formal error (22.11 μs) and high stability in the solutions, although external agreement accuracy still lags behind international levels. This research provides quantitative evidence for the precision efficacy of VGOS technology in practical applications and highlights its potential for high-frequency Earth orientation parameter monitoring.

Under an elastic Earth model, clear analytical solutions exist for surface gravity changes induced by a single strike-slip fault. However, the quantitative relationship between multiple fault zones and their geometric parameters with the resulting surface gravity change patterns remains unclear. Based on dislocation theory in an elastic half-space, this paper systematically investigates the quantitative connection between the geometric parameters of parallel strike-slip dislocation sources and the characteristics of the surface gravity change patterns through numerical simulations. Using the control variable method, we separately analyze the influence of fault length (L), width (W), slip amount (K), number of faults (N), dip angle (δ), and top depth (dep) on the four-quadrant distribution pattern and amplitude of the gravity anomaly. The results show that: 1) Fault length L is linearly positively correlated with the distance between gravity extremal points along the fault strike (X-direction). 2) Increasing fault width W elevates the gravity anomaly amplitude, but the growth trend diminishes, without altering the four-quadrant structure. 3) The strike-slip width K is linearly positively correlated with the distance between gravity extremal points along the Y-direction, but negatively correlated with the gravity anomaly amplitude. 4) An increase in the number of strike-slip faults N leads to a linear increase in the gravity change amplitude. 5) Variations in the dip angle δ affect the symmetry of the gravity change pattern. 6) An increase in the top depth dep causes a significant increase in the distance between extremal points in the Y-direction. 7) The fault center position (s_x, s_y) coincides with the four-quadrant symmetric center of the gravity change pattern, and the line connecting a pair of extremal points in the X-direction aligns with the fault strike angle α. In summary, the main geometric parameters of multiple dislocation sources can be determined based on observed gravity change patterns.

Based on the observation data of VP-type broadband vertical pendulum inclinometers from 5 stations in Inner Mongolia, this study adopted a method combining the wave propagation theoretical model and the deep learning prediction framework to preliminarily reveal the multi-frequency quantitative influence mechanism of the overburden thickness of observation caves on the observation signals. The results show that the overburden thickness is significantly negatively correlated with background noise (especially in the north-south direction), but has a weak impact on low-frequency signals. It exerts a directional regulatory effect on the amplitude of M₂ tidal waves: the amplitude of the north-south component increases linearly with the increase of thickness, while the east-west component shows a unique resonance amplification effect at the thickness of 20 m (the peak value is about 3.6 times higher than that at 10 m). When the thickness is greater than or equal to 20 m, it can effectively suppress the seasonal fluctuation of data, significantly reduce the variation amplitude, delay the peak value by 1-2 seasons, and simultaneously decrease the barometric admittance coefficient remarkably. Verified by the random forest model, it is found that the optimal thickness range of the cave overburden is 18-22 m (with 20 m as the optimal value). Lithology is a secondary influencing factor, while the overburden thickness and its interaction with the density of mountain lithologic media constitute the dominant control parameters. In this study, 20 m of overburden thickness was selected as the engineering threshold, which provides an important design basis for effectively suppressing background noise, ensuring the fidelity of tidal signals and resisting seasonal interference, and gradually promotes the construction of cave-based observation at stations from an empirical mode to a quantitative and engineering-oriented direction.

In this study, 1 216 consecutive days of ESA precise clock products (January 2022 to April 2025) were analyzed by combining clock-fit parameters, frequency-stability metrics, and spectral analysis to assess the long-term behavior of the onboard clocks. Results show that E14 experienced three phase jumps and two frequency jumps during this period, more frequently than nominal circular-orbit satellites. Following the frequency jump in April 2024, its frequency drift and fitting accuracy deteriorated markedly: residual peaks exceeded 10 ns, and an anomalous periodicity of about 0.63 h emerged with increasing amplitude and no evident secondary component. The frequency stability simultaneously worsened, showing oscillatory behavior near 100 s. After another frequency jump in March 2025, both residuals and stability improved, with stability reaching 2.50×10^-13, 6.52×10^-14, and 9.70×10^-14 at 100 s, 1 000 s, and 10 000 s, respectively. No dominant periodic term was observed afterward, but the drift exhibited a systematic bias, and a "dip" feature appeared near 4 000 s. In contrast, the co-orbital E18 satellite maintained stable performance throughout the same period, with fitting residuals confined within ±0.2 ns, no significant high-frequency periodic components, and frequency stability comparable to that of nominal MEO satellites, showing no anomalous behavior.

Landslide disasters pose a continuous threat to the safety of life and property of mountain residents. Especially under heavy rainfall, landslide areas are often accompanied by severe surface deformation, resulting in serious incoherence of synthetic aperture radar interference data, which restricts the deformation monitoring capability of traditional timing InSAR technology in this type of area. To address this problem, this paper takes the large landslide that occurred in Nanyu township, Zhouqu county, Gansu province in July 2018 as the research object, and prensents the distributed scatterer InSAR (DS-InSAR) deformation monitoring method based on statistical homogeneity detection and time phase link optimization. This method introduces a statistical homogeneity determination model under the traditional DS-InSAR framework, and identifies clusters of pixels with stable scattering characteristics through non-parametric testing, thereby expanding the monitoring coverage of low-coherence areas. At the same time, the shortest path constraint strategy is used to optimize the interferogram network structure to improve connection density and temporal coherence. It is also combined with a time-Phase Linking joint estimation algorithm to optimize the interference phase as a whole, enhancing the temporal continuity and stability of deformation information. The research results show that this method significantly improves the distribution density and temporal coherence of monitoring points in the collapse area of Nanyu landslide. The number of effective coherent points in creased from 23 232 to 43 463, an increase of about 87%, and the average posterior coherence coefficient was improved by 0.3, which effectively revealing the dynamic evolution characteristics of the landslide main body. The optimized annual average deformation rate chart shows that the maximum sliding rate of the main deformation area of the landslide exceeds 70 mm/a, and the timing curve of the characteristic point is highly consistent with the rainfall process. This study verifies the applicability and effectiveness of the timing phase optimization method in high-voltage landslide areas, and provides strong technical support for InSAR monitoring of landslides in complex mountainous areas.