In response to the imprecise localization of moving radiation sources in complex scenarios and issues related to the accuracy of the observation station's position and the geometric relationship with the target using the two stage weighted least squares(TSWLS) algorithm, we propose a passive radar centroid passive localization algorithm for moving targets based on Taylor expansion and the joint application of TSWLS. This method utilizes a TSWLS algorithm to determine the target's position and velocity. The obtained target parameters serve as initial values for Taylor expansion, constructing the positioning error equation. Through iterative optimization, the algorithm refines and solves for optimal target parameters. The proposed algorithm and TSWLS independently provide estimations, with the final result obtained through centroid localization. Simulation experiments show that, regardless of target speed, in comparison with traditional methods, the presented algorithm significantly improves robustness and positioning accuracy while reducing sensitivity to geometric relationships between the observation station and target.

To objectively evaluate the accuracy of TPXO series ocean tide models for the eastern China sea(ECS), we use the latest FES2014b and EOT20 as reference models. The assessment is performed by computing the root mean square error(RMSE) of the M₂ vertical ocean tide loading(OTL) displacements and the root sum square(RSS) of the 8 principal tidal constituents based on the GPS-observed M₂ vertical OTL displacements and tide gauge observations from the evenly selected 22 GPS sites and 65 tide gauges in the coastal areas, respectively. The results show that: 1) In the whole ECS, EOT20 is the most accurate model with RMSE of 0.41 mm and RSS of 11.1 cm, while the accuracy of FES2014b is slightly lower with RMSE of 0.43 mm and RSS of 11.4 cm. 2) Along the coast of eastern China and the western Korean Peninsula, there are large differences between TPXO models and FES2014b, and the RSS with tide gauge observations are 2-6 times as much as that of FES2014b and EOT20. 3) Both TPXO9_Atlas and TPXO8_Atlas agree better than other models with the M₂ vertical OTL observations in the Ryukyu Islands, and the two models are suggested for the local correction of OTL displacements.

The water level change series of lakes in the Qinghai-Xizang plateau from 2003 to 2022 are extracted using the data of ICESat, Cryosat-2 and ICESat-2. Taking Qinghai lake as an example, the three retracking algorithms of Cryosat-2 LRM(Ocean-CFI, UCL Land-ice and OCOG) are compared and analyzed, and the results show that OCOG algorithm is the best. The accuracy of the time series of lake water level changes extracted from satellite altimetry data is evaluated by using the hydrological station data of Qinghai lake, Nam Co, Mêmar Co and Lumaqangdong Co, and the results show that the RMS of the difference between the water level change time series of Qinghai lake and the measured results of Xiashe station is only 0.092 m, and the correlation coefficient between the water level change time series and the measured results of Nam Co is 0.989, and the accuracy analysis of different measurement techniques in the same time period shows that the accuracy of the Cryosat-2 SARIn measurement mode is slightly higher than that of ICESat-2. The water level changes of 106 lakes in the Qinghai-Xizang plateau from 2003 to 2022 are extracted, and the water level change trend is analyzed, and the results show that the water level of 19 of the 106 lakes show a downward trend, and the water level of 83 lakes show an upward trend.

This paper is based on the AR products from CODE, GFZ, CNES, PRIDE Lab, and in-house AR products for the performance assessment of PPP-AR in terms of convergence time, time to first fix and positioning accuracy. The experiment uses seven-day observation data from 40 IGS stations in 2022 and precise products of each AR product. The results show that under the static solution mode with confidence level of 95%, five products have the most obvious improvement in positioning accuracy at the time of 1 h, which are 46.58%(3.4 cm), 41.10%(3.0 cm), 45.21%(3.3 cm), 34.25%(2.5 cm), and 41.10%(3.0 cm), respectively. In the simulation-kinematic solution mode, the positioning accuracy of the five products in the E and N directions can reach the millimeter level, and the accuracy of the U direction is slightly improved than that of the floating point solution, among which the accuracy of GBM product is the least improved, and 72.73%(2.4 cm), 47.37%(0.9 cm) and 5.41%(0.2 cm) were increased in E, N and U directions. WUM and COM products had better improvement effect. They were 81.82%(2.7 cm), 63.16%(1.2 cm), 24.32%(1.1 cm) and 81.82%(2.7 cm), 63.16%(1.2 cm), 15.58%(0.8 cm).

Addressing the under-quantified crustal deformation characteristics and the unclear stress-strain coupling relationship of the Ordos block, this study employs a spherical least squares collocation model to calculate the characteristic strain field of the study area's crust and inverts the principal stress distribution using seismic source mechanism data. The results indicate that crustal deformation is primarily concentrated along plate boundaries and within fault basins, with less deformation within the Ordos block itself. Notably, at the junction of the northeastern Tibetan plateau and the Alxa block, the crustal deformation is particularly intense, characterized by strong compression along the plate boundary, with the main deformation mechanism being a coexistence of NW-SE compression and NE-SW extension, with a maximum principal compressive strain rate exceeding 0.05 μstrain. The crust in the study area is under a relatively uniform tensile stress environment, with the direction of the maximum principal stress being nearly horizontal and pointing towards SEE, aligning with the direction of the principal strain in the corresponding area. The minimum principal stress direction is closer to the direction of the principal compressive strain.

In this paper, we collect 1 Hz high-frequency data from multiple GNSS ground reference stations near the epicenter of the Jishishan earthquake on December 18, 2023, and use PRIDE PPP-AR and TRACK software for dynamic PPP calculation and high-precision static data processing. The coseismic deformation displacement fluctuation and permanent deformation displacement of the GNSS reference stations are obtained. Then, we analyze the deformation impact of the earthquake on the GNSS ground reference station at the epicenter. The results show that PRIDE PPP-AR and TRACK software can both calculate the coseismic deformation displacement impact of the earthquake on the reference station. The static calculation results show that the earthquake caused a permanent deformation displacement of about 18 mm in the horizontal direction at LXJS station (5.3 km away from the epicenter), and the impact range of this earthquake on the stability of the GNSS reference station is about 50 km. The rapid deformation analysis method and results can provide reference and data support for disaster prevention and mitigation, as well as real-time navigation and positioning services.

We comprehensively use InSAR and machine learning technology to evaluate landslide susceptibility in Jinyang county, Sichuan province. The data set is updated by interpreting landslides. Based on 12 evaluation factors, we use three models including random forest(RF), support vector machine(SVM) and extreme gradient boosting(XGBoost) for training in Python environment to complete landslide susceptibility mapping. We use ROC curve to verify prediction performance. We optimize the negative samples, obtain the landslide susceptibility evaluation results after sample optimization using machine learning and update the landslide susceptibility results using surface LOS deformation rate. The results show that the three machine learning models have good zoning effect, and the mapping effect of XGBoost model is best. The accuracy of XGBoost model after sample optimization is highest, and the AUC value reaches 0.95. The surface deformation rate obtained by SBAS-InSAR can reduce zoning errors and give timeliness to landslide susceptibility evaluation.

To study the activity characteristics and patterns of the landslide, this paper uses stacked interferometric synthetic aperture radar (InSAR) based on ascending and descending track C-band Sentinel-1A data to obtain the deformation rates of the Cheyuping landslide from January 2020 to December 2021 along the slope direction and in the vertical direction. We obtain the landslide displacement time series evolution results. The results show that the maximum deformation rate of Cheyuping landslide is more than -90 mm/a along the slope direction and more than -40 mm/a in the vertical direction, and the deformation of the landslide is larger at the leading edge and middle part of the landslide, and smaller at the trailing edge, which is characteristic of the slip-type landslide. According to the analysis of rainfall data and time series results, rainfall has a certain influence on the landslide movement, showing seasonal accelerated deformation characteristics. At the same time, the geological and stratigraphic lithology of the landslide area and the change of the water level of the Lancang river also accelerated the landslide activities to a certain extent. This study can provide technical reference for the interpretation of landslide deformation mechanism and disaster monitoring and early warning in the canyon reservoir area.

Based on the Sentinel-1A image data covering the epicenter, we use D-InSAR to calculate the surface coseismic deformation field of Xinjiang Wushi M_W7.0 earthquake. We use the quadtree method to downsample the coseismic deformation field, and optimize the geometric and physical parameters of the fault by using the downsampling data as constraint conditions. The results show that the maximum uplift values of ascending and descending orbits are 72 cm and 39 cm, respectively. The seismogenic fault is within the Maidan-Shayilamu fault, which is dominated by thrust and has a small amount of left-lateral strike-slip. The fault strike is 230°, the dip is 67°, and the slip angle is about 45°. The maximum slip on the fault plane is 3.53 m, the seismic moment is 3.78×10¹⁹Nm, and the calculated magnitude is M_W7.02. The results of Coulomb stress show that the Tuoshigan fault, Kuokeshale fault, Dashixia fault, Kurukeyuzumu fault, and northern Kuqi fault are all in the stress loading region, on which we need to pay more attention in the future.

By analyzing the problems existing in the cycle slip detection method of wide-lane phase minus narrow-lane pseudorange and geometry-free combination, we propose a cycle slip detection method that considers the GNSS satellite elevation factor and is suitable for universal GNSS observation data. After detecting the epoch position where the cycle slip occurs, we use the method based on spatial search and minimum criterion of objective function to repair cycle slip. We then validate the proposed cycle slip processing method by selecting satellite data of different elevations observed from the universal GNSS receiver located at a landslide hazard near the Linxia Highway in Gansu Province. The results show that the proposed method is effective. the cycle slip detection and repair method can realize the accurate detection and repair of cycle slips in carrier phase data with different elevations and can reduce the influence of GNSS data observation error on the cycle slip test. It is feasible for the cycle slip detection and repair in universal GNSS observation data.

Abnormal measurement noise leads to lower filtering accuracy and even divergence of GNSS/SINS integrated navigation system. Thus, firstly we propose the Sage-Husa adaptive Kalman filtering method(SHAKF) suitable for integrated navigation systems. Then, we construct a control factor according to the theoretical and actual estimates of filtered innovation covariance, propose a detection method to measure the start and end time of measurement noise RMS mutation, construct a forgetting factor based on the exponential function variation rule, and propose a Sage-Husa adaptive Kalman filtering method based on exponential weighted average(EWASHAKF). Finally, we apply EWASHAKF to GNSS/SINS integrated navigation system and carry out the simulation experiment. The results show that compared with SHAKF, EWASHAKF can accurately track the variation of measurement noise RMS, and improve the filtering accuracy of integrated navigation system.

Based on the seismic ambient noise cross-correlation method, we calculate the relative velocity variation within the underground medium. We utilize continuous waveform data from the CCH and BDS stations at Dongtan mine, spanning from 2 a.m. on February 8, 2022, to 7 a.m. on March 8, 2023. The analysis focuses on the correlation between these velocity changes and both the mining operation and seismic activity at the 63^up06 working face. The findings reveal a notable impact of the mining operation at the 63^up06 working face on the wave velocity within the underground medium. Additionally, a clear correlation exists between the changes in wave velocity and the incidence of mine earthquakes, highlighting the interplay between variations in underground stress levels and seismic events. Furthermore, the release of underground stress during mine earthquakes can lead to slight fluctuations in wave velocity.

Based on the digital seismic observation report recorded by the Inner Mongolia seismological network, using Pg and Sg wave data of 3 398 earthquakes with M≥1.0, we deploy the double-difference tomography method to study the underground structure in the northwest margin of the Ordos block(37°~43°N, 104°~110°E) and conduct the earthquake relocation. At the same time, we analyze in detail the velocity structure of the upper crust in the area. The results show, from the velocity structure of Pg wave, there is a continuous irregular high-speed body from northeast to southwest at 5 to 10 km, and the high-speed body at 15 to 20 km is no longer continuous but scattered. From the velocity structure of the Sg wave, the velocity structure at 5 km and 15 km is basically consistent with that of the Pg wave. There are more low-speed areas at 10 km and 20 km, and earthquakes mostly occured at the junction of high and low velocity bodies and around faults, which proves the fragility of the crust in this area.

We propose a method for ground deformation monitoring based on distributed scatterer InSAR(DS-InSAR). The method employs the HTCI algorithm for homogenous pixel identification, utilizes eigenvalue decomposition for phase optimization, determines distributed scatterer(DS) based on spatiotemporal coherence, and establishes a relationship model between DS phase and deformation for time series surface deformation inversion. We use 50 scenes of Sentinel-1A images to obtain spatiotemporal distribution information of ground deformation in Jiangbei region of Nanjing from October 2021 to October 2023. The results show that: 1) The correlation coefficient between DS-InSAR results and PS-InSAR results is 0.731, and the point density increases by 6.22 times. DS-InSAR can more completely reflect the spatiotemporal surface deformation of the study area. 2) DS-InSAR is in good agreement with the measured leveling data, with a maximum error of 5.6 mm, an average error of 2.6 mm, and a root mean square error(RMSE) of 3.0 mm. 3) The geological structures, urban development, and rainfall are the main influencing factors causing ground subsidence.

We lay a shallow seismic exploration line L1 with a length of 10.77 km extended southward across Serteng mountain piedmont fault. The results show that the Sertengshan piedmont fault is clearly visible on the profile, and the main fault has bifurcations in the upper part, consisting of faults F_P1.1 and F_P1, both of which are normal faults with an apparent inclination towards the south. Combined with geological and geomorphological surveys, the latest active age of the F_P1.1 fault is Holocene. We carry out a composite drilling geological section for the fault F_P1 of the Serteng mountain piedmont fault and distinguish two faults, Fa and Fb. Based on the results of optical stimulated luminescence dating, the faults Fa and Fb are Holocene faults and the latest active age is later than 2.5±0.6 ka BP. The Serteng mountain piedmont fault has a forked structure near the surface, which reflects the multiple activities of the fault in the late Quaternary.

Early research on broadband seismometer azimuth mainly based on the P-wave particle polarization principle and calculates the angle corresponding to the minimum value of P-wave energy on the lateral component to obtain the angle between the north-south direction of seismometer and the geographic north pole. In addition, the elliptical motion characteristics of Rayleigh waves can be used to estimate the seismometer azimuth of stations with azimuth deviation. We estimate the seismometer azimuth deviation of 20 stations of CENC based on the two different methods. The results show that the azimuth deviation calculated by the two methods is in good agreement, indicating that Rayleigh wave analysis can be appropriately added in future seismometer azimuth detection.

In view of the problem that GNSS observation data in complex environments such as landslide is affected by factors such as high occlusion and multipath, it is easy for the success rate of ambiguity resolution to decrease and the positioning reliability to deteriorate. We propose an improved algorithm of GNSS partial ambiguity resolution based on IGGⅢ robust estimation. The robust estimation algorithm is introduced to weaken the influence of gross error on parameter estimation and improve the accuracy of floating-point ambiguity. In addition, to avoid the influence of residual error on the subsequent integer ambiguity resolution, and to further improve the success rate of ambiguity resolution, we adopt a partial ambiguity resolution improvement strategy based on conditional variance matrix method. The experimental results show that compared with the full ambiguity resolution method, the algorithm can significantly improve ambiguity resolution rate, and the results of three stations are increased by 5.9%, 52.7% and 48.0%, respectively. At the same time, compared with the traditional Kalman filter algorithm, the improved robust estimation algorithm significantly improves the situation of incorrect ambiguity resolution and greatly improves the GNSS positioning accuracy. The horizontal RMS is better than 1 cm and the vertical RMS is better than 2 cm.

Based on the observation data of the reservoir network in the lower reaches of the Jinsha River, we utilize the PhaseNet model, a deep neural network, to pick up seismic phase arrival times. By combining the methods of group triggering and isochronous octree searching, we associate seismic phases and use the NLLoc non-linear location method to automatically locate the microseismic sequence in Yuanbaoshan, thereby quickly constructing a high-precision machine learning catalog. The seismic numbers provided by the machine learning catalog are 3.89 times that of the manual catalog, with a matching rate of 94.76%. 89.8% of the events have a time deviation of less than 0.5 seconds, 98.7% of the events have an epicenter location deviation of less than 3km, 90.9% of the events have a source depth deviation of less than 5km, and 91.5% of the events have a magnitude deviation of less than 0.2. Based on the focal mechanism solution and the seismic distribution characteristics provided by the machine learning catalog, we infer that the seismic events of this sequence are caused by the sub-branch fault of the Ebian-Jinyang fault. The earthquake catalog construction method of this study can provide rapid data support for determining earthquake trends in the reservoir area downstream of the Jinsha River, precise aftershock locating, determination of focal mechanism nodal planes, detection of hidden faults, seismic event structures, and body wave travel time tomography imaging.

We focus on the complex characteristics of groundwater level data, including nonlinear trends, seasonal fluctuations, and random disturbances, and introduce the Prophet time series prediction model developed by Facebook. The aim is to use its nonlinear trend capture, seasonal fluctuation analysis, and flexible response ability to outliers and data missing to significantly improve the accuracy of groundwater level anomaly identification. Through observation data from Beilin district seismic station in Suihua city, Heilongjiang province, it is shown that the Prophet model performs well in capturing dynamic characteristics of time series data and can effectively identify anomalies. The high fitting accuracy and predictive ability of adjusted model have been confirmed, with low prediction error and high determination coefficient. In addition, the model identifies water level anomalies related to earthquakes in earthquake prediction, providing a new perspective for earthquake precursor research. This study demonstrates the effectiveness of Prophet model in processing complex time series data, providing a new tool for earthquake prediction.

Large earthquakes can simultaneously excite numerous normal modes of Earth's free oscillations. Due to the Earth's ellipticity, rotation, and internal anisotropy, these normal modes undergo splitting, with frequencies of splitting modes becoming very close(only a few μHz apart). This imposes higher demands on the detection of Earth's free oscillation modes. This article derives and verifies a new method for detecting free oscillation modes based on normal time-frequency transform(NTFT). Taking the detection of ₃S₁ splitting modes as an example, compared with the classical FT spectrum method and the latest OSE method, our method has higher frequency resolution for detecting free oscillation modes.

GNSS technology is widely used in landslide monitoring, but it is difficult to install monitoring equipment in complex high-risk landslide environment. UAV-dropped is expected to realize unmanned deployment of monitoring equipment. In view of the problems of UAV-dropped endurance and terrain threat in complex environment, the path planning of UAV-dropped is particularly critical as the basic mission. In this paper, the path planning problem of UAV-dropped in complex mountainous area is studied in three aspects: 3D real map construction of geological disaster area, cost function design and flight path planning algorithm. In addition, whale optimization algorithm based on adaptive weight and Levy flight strategy is applied to the flight path planning of UAV-dropped. Taking Hongyanzi landslide in Hanyuan county, Sichuan province as the study area, the path planning of UAV-dropped GNSS monitoring equipment was realized.

On the basis of effectively handling the gross errors of landslide monitoring data and fully considering the characteristics of landslide monitoring data, we develop a deep learning landslide displacement prediction model combining time series decomposition and similar component reorganization. First, we deal with the significant gross errors of landslide time series monitoring data using the isolation forest algorithm, and then comprehensively analyze its smoothness, autocorrelation, and normality to determine the optimal length of input feature sequence. Second, the non-stationary landslide monitoring data are decomposed into multiple smooth time series using the ensemble empirical mode decomposition(EEMD) method, which is then classified into three types combining the sample entropy and K-means algorithm, namely high, medium, and low frequency. Finally, comparing the prediction accuracy of different neural networks, the prediction models suitable for three types of time components are constructed respectively, and then the prediction results are superimposed to realize the high-precision prediction of landslide displacement. The testing results of Beidou/GNSS monitoring data of typical landslide body in experimental area show that the combination prediction model proposed in this paper has a better applicability to the landslide monitoring data containing significant gross errors, and can significantly improve the prediction accuracy of landslide displacement compared with single and existing combination models.

By utilizing satellite time-monitoring data to mitigate the impact of satellite-side and transmission path errors in standard single-point positioning, the timing stability and positioning accuracy at the user end are improved. At the same time, the timing deviation of the user relative to the satellite navigation system is corrected to be relative to the national standard time. Taking GPS as an example for a short baseline experiment, the results show that in terms of timing, the mean timing deviation of the L1 frequency band was improved from -42.8 ns to 5.4 ns, and from -53.6 ns to 2.8 ns for the L2 frequency band; in terms of positioning, the mean positioning errors in the X, Y, and Z directions for the L1 frequency band were improved by more than 55%, and by more than 77% for the L2 frequency band. This proves that by using satellite time-monitoring data, users of the satellite navigation system can not only enhance the pseudo-code positioning accuracy to the sub-meter level but also obtain the deviation relative to the national standard time, achieving ns-level timekeeping of the national standard time.

Tonga volcano violently erupted at 12:26:30 on January 15, 2022. The atmospheric pressure observation station of Yunnan Geophysical Network clearly recorded the intense pressure change process caused by the eruption. The atmospheric pressure change forms and characteristics of all observation stations were highly consistent, with the average peak amplitude of 2.16 hPa and the average propagation velocity of 314.30 m/s. The wave train has the properties of Lamb wave propagation, and the atmospheric pressure observation data are basically reliable. The water level observation in Yunnan deep well clearly recorded the water level change caused by the short-time pressure change. The static pressure efficiency of the observation well was estimated by the ratio of peak water level change to the amplitude of pressure change. We analyze by regression the short-time pressure change and well water level data, and take the determination coefficient, regression line slope and static pressure efficiency as indicators. Then we analyze the bearing property and stress-strain sensitivity of the borehole water-bearing system. The results show that the water level and pressure observation data of 25 wells, such as Kunming station, have high goodness of fit, small lag time, and good permeability of aquifer. The atmospheric pressure efficiency of 7 wells in Kunming station, Daguanxijiagou, Tengchong station, Ludianciyuan, Nanhua, Shiping and Zhaotong SK3 is slightly higher, and the sensitivity of well water level to stress and strain may be low. The larger the coefficient of determination for an observation well, the higher the probability that its monitoring and forecasting performance will be assessed as qualified.

Aiming the issue of reduced accuracy in sea surface height inversion using global navigation satellite system interferometric reflectometry (GNSS-IR) due to elevation angle deviation, a correction method for elevation angle deviation based on the ionospheric region is proposed. To verify the reliability of the correction method for elevation angle deviation, a comparison is made between the inversion results before and after applying the correction to the elevation angle deviation. The results indicate that the inversion results, after correction for elevation angle deviation, exhibit higher accuracy. The method for correcting elevation angle deviation can effectively rectify the deviation caused by the ionosphere.

For the GPS modernization by adding L₅ frequency observable values, we propose a new algorithm of GPS single epoch double-difference integer ambiguity resolution called GPS Main and Auxiliary Frequency Correlation Algorithm (MAXCOM). We use only GPS single epoch data, including pseudo-range observation values and phase observation values of main-frequency and auxiliary-frequency for the algorithm. We set L₁ is set as the main-frequency and L₂ or L₅ as the auxiliary-frequency by expanding the idea of dual frequency correlation method to establish an internal-correlation model between double-difference phase observation values of GPS main-frequency and auxiliary-frequency. We conduct a cross-search between the adjustment value-domain of observation values and the value-domain of ambiguity, and then use the ratio significance test method to determine the double-difference integer ambiguity of GPS main-frequency. The example shows that for the three original phase observations of L₁, L₂ and L₅ frequency after GPS modernization, selecting L₅ as auxiliary-frequency of the algorithm is slightly better than L₂ in terms of efficiency and success-rate of ambiguity resolution. The proposed algorithm is effective.

On August 23, 2023, a M4.6 earthquake hit the Pulandian, Liaoning province. To describe the source characteristics of this earthquake, we explore its incubation and seismogenic mechanisms. In this paper we determined the focal depth of the Pulandian M4.6 earthquake, and calculate the focal mechanism solution, moment tensor solution and centroid depth of the earthquake, obtain the center of focal mechanisms of the earthquake. In the same, we analyze the relationship between focal mechanism and tectonic stress field, and fit the fault plane according to the results of small earthquake relocation. We preliminarily determine the seismogenic fault plane of this earthquake. The results show that the initial rupture depth of the Pulandian M4.6 earthquake is 12.0 km. The focal mechanism solution is strike: 50°, dip: 75° and rake: －169° for nodal planeⅠ, strike: 317°, dip: 80°and rake: －15° for nodal plane Ⅱ, the moment magnitude is M_W4.8, while the centroid is 12 km. The seismic moment M₀ is 1.796×10¹⁶ Nm, the moment tensor solution is M_rr: －0.004, M_tt: 0.946, M_pp: －0.942, M_rt: 0.017, M_rp: －0.305, M_tp: －0.125. The center of focal mechanisms solution is strike: 47.03°, dip: 79.04° and rake: －168.15° for nodal planeⅠ, strike: 314.75°, dip: 78.37°and rake: －11.19°for nodal plane Ⅱ. The relative shear stress and relative normal stress generated by the tectonic stress field acting on the nodal plane Ⅰ are 0.877 and －0.544; on nodal plane Ⅱ they are 0.911 and 0.161. The fitted seismogenic fault plane is strike: 148.91°, dip: 89.85° and rake: 26.47°. The analysis shows that the Pulandian M4.6 earthquake occurred on the Pulandian-Changhai NW tectonic belt. It is a natural earthquake that after stress accumulation occurs along the optimal nodal plane of the stress field and takes left-lateral strike-slip as the dislocation mode.

Using the satellite fire monitoring data accumulated by an independent development geostationary meteorological satellite system in the western Pacific from 2020 to 2023 and the VOTW Holocene volcanoes data released by the global volcanism program(GVP), we propose a volcanic activity monitoring algorithm based on geostationary meteorological satellite fire observation data. We use cluster recognition technology to extract volcanic fire from satellite hotspot data and calculate the volcanic normalized thermal index. The comparison of volcanic monitoring results with the GVP volcanic activity report shows that the identification rate of volcanic activity is about 90%. Due to the type of volcanic activity and the impact of tropical convective cloud systems, the missed reported volcanoes are mainly located in Indonesia. The spreading direction, distance and speed of lava flow can be estimated according to the geographical distribution and time variation of the volcanic fire. The satellite volcano monitoring data can be used as input variables in real-time mode to drive the volcanic activity forecast model to issue volcano warnings in real-time. The accumulated data can also be used as high time resolution data sources to support the study of volcanic activity.

In order to investigate the applicability of InSAR technology in CO₂ geological storage safety monitoring, we use 43 Envisat images from October 2003 to September 2010, and obtain the deformation rate field around CO₂ injection well during the implementation of In Salah CCS project based on StaMPS-InSAR and SBAS-InSAR technologies. Furthermore, the underground CO₂ reservoir parameters of KB-501 injection well are inverted using the geodetic Bayesian method. The results show that: 1) With the injection of CO₂, deformation of 2-4 mm/a occurred around three injection wells, and the maximum cumulative deformation during monitoring period reached 32 mm. 2) Through the joint analysis of surface displacement time series and CO₂ injection data, we find that CO₂ injection rate significantly affects surface uplift rate. 3) Based on the horizontal Sill model, we use the Bayesian method to invert the depth of underground CO₂ reservoir in KB-501 injection well, and the best fitting value is 1 897 meters. The strike angle of horizontal edge is 37°, and the tensor of dislocation plane is 0.04 meters, which proves the applicability of horizontal Sill model and Bayesian method in geophysical processes analysis related to CO₂ geological storage.

We conduct analysis on Sentinel-1A/B datasets(2015 to 2022) across the Xianshuihe fault zone(XSHFZ) and obtain the interseismic line-of-sight(LOS) deformation fields. Under the constraint of the fault parallel horizontal velocities derived from the decomposition of the ascending and descending InSAR-derived LOS deformations, we invert the interseismic kinematic parameters of the XSHFZ. Inversion results show a spatial variation of the tectonic loading rate(3.5-13.0 mm/a), locking depth(14.7-26.2 km), shallow creep rate(3.2-12.5 mm/a), and the corresponding top depth(0.3-2 km) and bottom depth(3.4-13.7 km). Combining the historic seismicity and the inferred interseismic kinematic parameters of the XSHFZ, we calculate the seismic moment deficit rate, accumulated seismic moment and moment magnitude of the Luhuo, Daofu, Qianning, Kangding, and Moxi segments during the interseismic stage, and then quantitatively evaluate the seismic risk of different segments. Results show that all segments have the potential to produce an earthquake with magnitude larger than M_W6.5, and there is a risk of earthquakes with magnitude of ~M_W6.9 on the Luhuo and Qianning segments. We should pay more attention to the risk on these faults.

Using the BIC_tp noise model identification criterion, we analyzed the temporal noise characteristics of sea surface height (SSH) and their influencing factors at 350 virtual coastal satellite altimetry stations worldwide. Furthermore, we explored the impact of environmental loads on SSH time series and the characteristics of sea level changes across different time spans. The results indicate that the assumptions of FN+WN, FN+RW+WN, and PL+WN noise models can lead to errors in estimating the velocity parameters of SSH time series. Furthermore, there is a diversity of SSH time series noise models, and the velocity uncertainties under different noise models vary significantly. Identifying their noise characteristics is crucial for accurately obtaining sea level change parameters. The average global sea level rise rates for the different time spans of 1993 to 2006, 2007 to 2020, and 1993 to 2020 are 2.59 mm/a, 4.41 mm/a, and 3.08 mm/a, respectively, indicating that sea level is accelerating its rise. The optimal noise models and velocity changes for the SSH time series before and after environmental loading correction are not significant, indicating that the load effect has a relatively minor impact on the characteristics of SSH time series.

We calculate local terrain and isostatic correction using spectral combination and nested different topographic resolutions and determine the free air gravity anomaly grid by the gravity reduction method. The results show that the accuracy of terrain correction calculated by the spectral combination method is the best in the test area, and the accuracy of isostatic correction calculated by nested different resolutions in the central and edge regions is best. When the integral radius is 1° and 2°, there is little difference between the local terrain correction, and the isostatic correction results differ significantly. The accuracy of free air gravity anomaly grid results obtained from the complete Bouguer and isostatic gravity anomalies by gravity reduction is the same, and the accuracy of the results is much higher than that obtained from the direct grid of free air gravity anomalies at ground points.

Based on the observation data from 163 GPS continuous stations in Qinghai-Xizang plateau and adjacent areas from 2010 to 2020, we solve the GPS vertical coordinate time series. Using the global loading grid data provided by GFZ, the environmental loading displacements are calculated by double cubic linear interpolation method, and the GPS vertical time series are corrected. The results show that the hydrological loading has the greatest influence on GPS station displacement in the study area, which is larger in the southern stations, with a RMS of 9.45 mm(MUET site). By comparing the changes of weighted root mean square(WRMS), annual term amplitude and velocity uncertainty before and after environmental loading correction, it is found that the environmental loading correction is able to weaken the nonlinear variations of vertical time series of most GPS stations in Qinghai-Xizang plateau and adjacent areas and it can accurately improve the velocity accuracy of most GPS stations.

Based on the images from Gaofen-1 satellite, taking the first section of the Three Gorges reservoir as an example, 160 landslide samples are visually interpreted and divided into training and validation samples according to the ratio of 9:1. Then, landslide automatic identification models are constructed using Mask R-CNN and Faster R-CNN algorithms. To further compare and analyze the performance of the two models under different sample ratios, calculations are performed using sample ratios of 8:2, 7:3, and 6:4. The results show that the recognition results based on Mask R-CNN model are better than Faster R-CNN in three indicators, including precision, recall and F1 score. Moreover, cross-validation proves that the performance of Mask R-CNN model is more stable.

Aiming at the problem of missing data and gross errors in GRACE Level1B observations, we propose an optimization method to complete the missing data of SCA1B. At the same time, we use an optimization strategy to eliminate gross errors for KBR1B and kinematic orbit data, then use these data sets to invert the time-varying gravity field model. In addition, based on synthesized data and real data respectively, we analyze the impact of the error from ACC1B in the Y direction on the inversion results and propose an optimization correction strategy. The optimized completion method can recover the missing data of SCA1B effectively, and fully considers the variational characteristics of the observations throughout the arc. The accuracy of the result derived from ACC1B data after correction using the optimization strategy is 3.7 mm, which is higher than the result derived from the uncorrected data. It can be seen from the results that the overall accuracy of the gravity field model solved by the data processed using the optimization method is like that of the three official centers. However, the detailed signals performance of the results from different institutions in regionals are different. It shows that the optimized data preprocessing strategy is effective and feasible.

This paper uses small baseline subset InSAR(SBAS-InSAR) technology to carry out long-term monitoring of landslide body deformation, and studies the optimization of gated recurrent unit(GRU), long short-term memory network(LSTM), support vector machine(SVM), differential autoregressive moving average model(ARIMA) prediction periodic and trend deformation through long-term monitoring data. Based on the results of the optimization analysis, this paper designs a landslide body displacement prediction model based on the time series InSAR inversion results and the GRU-ARIMA combination. Taking a landslide body in Zhaotong, Yunnan province as an example, the paper verifies and analyzes the prediction model. The results show that the combined displacement prediction model can take into account the accuracy of the results of the periodic term and the trend term. The error between the predicted value of landslide displacement and the root mean square of InSAR inversion data is 1.15 mm and R2 is 0.996 mm. The experimental results showed that the combined prediction method had high consistency with the actual observation.

Based on the analysis of the daily occultation data provided by FY-3E and its global distribution characteristics, we evaluate the accuracy of foF2, hmF2, and TEC derived from FY-3E using the critical frequency(foF2) and peak height(hmF2) obtained by the ionosonde and the total electron content(TEC) provided by the global ionospheric map(GIM) and the international reference ionosphere model(IRI2016). The results show that the correlation coefficients of foF2 and hmF2 between ionosonde and FY-3E are 0.86 and 0.75, and the root mean square errors are 0.99 MHz and 19.26 km, respectively. Taking the GIM TEC data as a reference, the accuracy of TEC retrieved from FY-3E GPS and BDS data is approximately equivalent. Moreover, the accuracy of FY-3E TEC is negatively correlated with solar activity. The correlation coefficient between FY-3E TEC and GIM TEC is 0.90, and the root mean square error is 5.23 TECu. The correlation coefficient between FY-3E TEC and IRI2016 TEC in the same height range is 0.88 and the root mean square error is 3.70 TECu.

We use the ESPAC method to detect the geological structure below the Mogao Grottoes of Dunhuang using observed data obtained by linear array. The detection results show that the layers determined by microtremor survey method are clearly visible, with significant fluctuations of each layer on the two-dimensional microtremor apparent S-wave velocity profile. The covering layer is thick at both ends and thin in the middle. There is a bedrock intrusion in the middle of profile, and there is a solitary stone with a width of 30 m×40 m. Using geological survey data calibration of the relevant region, we show that the thickness of Quaternary cover layer represented by mixed backfill, fine sand, and gravel sand is about 8 m, and the burial depth of complete bedrock surface is about 30 m. The research results can provide geophysical basis for the protection of cultural relics in the Mogao Grottoes.

Current robust fitting algorithms are prone to failure in point cloud data with a high contamination rate. To solve this problem, we introduce a weight matrix in the plane fitting method based on eigenvalue decomposition, where the initial weight is determined based on the normal vector information of point cloud, thus reducing the influence of outliers. We further use the scale-adaptive method to enhance the robustness. The scale factor of this method gradually decreases as the number of iterations increases, effectively solving the problem of low breakdown ration of traditional M estimation. The experimental results show that the breakdown ratio of the proposed algorithm is higher than that of existing robust point cloud fitting methods, and it can accurately estimate plane parameters in the point cloud data with a high outlier ratio.

We use the precise clock offset data released by IGS to analyze the key performance indicators of six GPS-Ⅲ satellite clocks that have been launched and operated in orbit, including frequency accuracy, frequency drift rate, and frequency stability, and compare them with the previous GPS clocks. The results show that rubidium clocks carried by GPS-Ⅲ are generally stable in operation and have good data quality, which is related to the improvement of anti-interference ability. There is no significant change or difference in time frequency performance, which is basically consistent with the previous generation Block ⅡF rubidium clocks. The frequency stability of GPS-Ⅲ rubidium clocks are better, which helps to improve the accuracy of PNT services. The research can provide reference and guidance for the upgrading and transformation of other GNSS system satellites.