Modern continental dynamics and structure geology have determined that shell with control weak interaction between crust and mantle not only restricts the coupling relationships between the spheres of rock, but also affects the interaction plots. The coupling effect can compare the different depth of the deformation and stress characteristics analysis. The widespread low velocity zones at the depth of 20－35 km in the southeastern margin of the Qinghai-Tibetan plateau may have potential decoupling functions in the present seismic and tectonic activities due to melting from the inner part of the crust. This paper, based on the 2015－1999 GPS coordinate time series data to obtain the velocity field, combined with the previous SKS given shear wave data, analyzes the horizontal deformation coupling characteristics; further it collects the Moho and lithosphere bottom deformation data, showing that the deformation characteristics are vertical upward. The use of EGM2008 gravity model 11－36 order spherical harmonic coefficients are calculated in the Sichuan-Yunnan region under the lithosphere mantle convection stress field, collecting and sorting the 1 131 focal mechanism solutions of the data since 2000, the regional stress tensor inversion method of damping the seismogenic layer stress field of the region; and further analysis of the two kinds of stress field correlation. In the ABAQUS software, the 3D lithospheric viscoelastic finite element model of Sichuan-Yunnan region is established, while adding the upper crust boundary displacement constraints and the bottom surface of the mantle convection drag loading conditions. The stress and strain in horizontal and vertical directions are analyzed respectively, and the results of strong earthquake profile are given. The seismogenic mechanism of crust mantle coupling is discussed, and the seismic risk of faults in Sichuan-Yunnan region is analyzed. Through the above research, the dissertation infers that the uneven development of stress is the key of strain energy accumulation, and vertical coupling relationship between the spheres has greatly influenced the block effect. Further, the shell weak zone blocks the extensive development of stress stratification, the edge area of transverse and longitudinal stress changes high volt, crustal seismic hazard and contact.

Taking 11 cities of Hebei province as examples, the correction models of MODIS PWV based GPS are studied. Firstly, we compare GPS PWV and MODIS PWV. There is good correlation and existing deviation between GPS PWV and MODIS PWV; the correlation coefficient is more than 0.879 and the root mean square error amounted to 2.31 mm. Secondly, we construct the cities’ and regional MODIS PWV correction models using linear regression method. Finally, we verify the reliability of city and regional models compared with GPS PWV. It is concluded that city and regional models of MODIS PWV can effectively improve the accuracy of MODIS PWV, while the root mean square error is less than 1 mm. The accuracy of MODIS PWV correction model can meet the application of weather forecasting.

Giving error analysis on three frequently used zenith tropospheric delay correction models, including the Hopfield, Saastamoninen and EGNOS models, combined with the errors of meteorological elements and station position, this study obtains the degree to which the errors affect the assessment value of ZTD. Then, achieving the three models by C++ language, selecting a series of IGS station at different latitudes and elevations, using meteorological data provided by IGS Analysis Center, combined with the time and space information of the station, this paper calculates the real ZTD. Compared with the IGS troposphere products, the estimation of quality of the three models is achieved. The results show, according to the meteorological data, Saastamoninen and Hopfield model could reflect the changes of tropospheric delay in a certain day accurately. The corrective precision of Saastamoninen model is slightly higher than that of Hopfield model. The EGNOS model does not need the measured meteorological data. Meanwhile, its RMS is less than 0.1 m, and can also meet the requirement of GNSS meter level localization accuracy.

As the real-time ionospheric delay grid correct information of the Beidou wide area difference system is calculated based on pseudo-range observation and phase observation with higher observation accuracy, so this paper presents an ionospheric delay corrections algorithm based on carrier phase and pseudo-range combined observations. The experiment is carried out using the GPS data of land network to validate the algorithm. The results show that ionospheric delay corrections based on the new algorithm are smoother, more stable, and better than the pseudo-range observations. The mean inner precision of ionospheric delay grids is improved by 30.2%, and the mean external precision is improved 27.8%.

The vertical total electron content(VTEC) is obtained using uncombined precise point positioning technique, while a series of COSMIC ionospheric occultation products which are near GNSS stations are collected. Both are used as common constraint conditions to retrieve electron density profile above these stations using functional model in near real-time. To verify the feasibility of this new method, GNSS observation data from the last weeks of March, June, September, and December in 2016 at 11 IGS stations in east Asia, as well as corresponding COSMIC ionospheric products(ionprf) occurring in or near these stations are collected. The results are compared with COSMIC occultation electron density profile products, which are seen as true values. The preliminary statistical results show a good consistency between the retrieved electron density profile using the ground-based singular GNSS station and space-based occultation products. The precision of f0F2 and hmF2 parameters are both better than those of IRI2016.

In this article, a new method called IEEMD (improved ensemble empirical mode decomposition) is introduced to build the correction of multipath error. Compared to EEMD (ensemble empirical mode decomposition), the new method can effectively solve the problem of model mixing and incomplete decomposition. It can also decompose the original coordinate sequence into different scales of modal components. The product of Gaussian white noise energy density and the average period is constant. Taking this condition into account, a new method is designed for automatic selection of scale and reconstruction to build the multipath error model. We can use it to correct the later coordinate series by the first day’s multipath model by the reason of the strong correlation between two successive days. The experiment results show that IEEMD can get more accurate coordinate sequence than EEMD, based on the same way of scale selection and stellar day filtering. To a certain extent, this study solidifies the theoretical basis of GPS high-precision real-time dynamic deformation monitoring.

Beidou lacks effective methods of navigation performance assessment in its representative applying scenes. This paper puts forward a way to evaluate the performance of Beidou user terminals based on the modeling of the scenes. This paper analyzes the possible challenge environment that Beidou might face and lists the elements to construct it. Then, the hilly scene is built following the steps. On this condition, the visibility of satellites and the precision of navigation are listed as major parameters to evaluate the navigation performance of user terminals. In the end, we construct a specific scene based on DEM and analyze the parameters mentioned above. The research demonstrates that the scene we have constructed can show the navigation performance in real life. The navigation performance assessment is practical in these circumstances.

GPS/BDS dual-frequency RTK positioning is fast and reliable, and is most widely used in high-precision dynamic positioning. However, its receiver is expensive and difficult to be widely used in urban environments. It is easy to promote GNSS single-frequency RTK because of low-price GNSS receiver. In this paper, we introduce a single-frequency BDS+GPS RTK algorithm, taking account of velocity information to assist in ambiguity fix, the inheritance of ambiguity, and the fusion of Doppler observation. Through an actual on-board test in urban environment, the single-frequency RTK positioning performance of BDS, GPS and BDS+GPS are compared and analyzed.

Galileo satellite navigation system has begun to provide initial services. GAL-201 and GAL-202, which did not enter the scheduled orbit, have broadcast navigation messages as well. Based on the measured data of MGEX, the positioning performance of the current Galileo system is evaluated and analyzed. The results show that: 1) Galileo standard point positioning (SPP) can reach positioning accuracy 2 m in horizontal, 4 m in vertical; 2) GAL-201 and GAL-202 can now be used as SPP solvers. After their introduction, the average positioning availability increases by 10% and reaches more than 80%; 3) when the observation time is 3 hours, Galileo precise point positioning (PPP) can only reach an accuracy of 5 cm in horizontal,7 cm in vertical. However, with the increase in observation time, accuracy gradually improves. When the observation time is 24 hours, the horizontal accuracy is better than 1 cm, and the vertical accuracy is better than 3 cm.

Receiver antenna phase center variation is one of important error sources influencing the precision orbit determination of low earth orbit satellites. In this paper, based on the GPS data of ZY3-01 star, the PCV models are estimated using the reduced-dynamic method and the residual method, and the influence of PCV on the precision orbit determination is analyzed. After using the PCV model, 3D accuracy in overlapping arcs comparison can increase 4.5mm, and SLR residual can increase 1.2 mm. All SLR station checking results also increase.

This study is concerned with the difficulty of parameter selection in the LSSVM model and the selection of parameters for high dimensional reconstruction of single variable sequences. To this end, we introduce the phase space reconstruction theory, and fruit fly algorithm into the LSSVM model. The deformation prediction model of phase space reconstruction based on FOA-GLSSVM is established. In order to validate the effectiveness and reliability of the proposed model, we compare practical examples and the GLSSVM model, least square support vector machine model and support vector machine model. The experimental results show that the proposed model is more accurate and stable.

Three InSAR products and DEM data are involved in detecting potential landslides over Wudongde Hydropower station section, Jinshajiang. Several known and unknown landslides are successfully detected, and their spatial patterns and stabilities are detected, showing that this procedure is an efficient means to detect active landslides in a large area. Meanwhile, small baseline subsets InSAR technique is applied to monitor the Jinpingzi landslide. Both spatial zonation features and time-series deformation results in key landslide zone are achieved. The InSAR precision in slide direction is 1.8 cm by comparing with in-situ total station measurements. Different InSAR techniques are illustrated in this paper for different scale landslide applications.

In spatial linear fitting, it is not possible to directly adopt the method of total least squares algorithm and mixed total least squares. This paper proposes a space linear fitting algorithm based on the PEIV model of the total least squares method. The standard equation of space straight line is transformed and the equation is rewritten as the overall least squares EIV model. Then, according to the characteristics of the coefficient matrix, the model is transformed into a more reasonable PEIV model, linearly similar to the least squares indirect adjustment form, and the iterative method is used to solve the fitting parameters. The adjustment process ensures that the correction number of the elements of the coefficient matrix is the same, and the correction number of the constant element is zero, which is in accordance with the actual theory. Finally, the feasibility and superiority of the method are verified by an example.

In order to discuss the impact on the adjustment of gravity network structure, the latest four groups gravity data of the upgraded northern Tianshan gravity network are used to analyze the impact of the long open line in gravity network. The results show that the existence of long open line in the gravity network brings the maximum gravity change is up to 260 μGal. Meanwhile, when the open lines are converted into a closed loop, the maximum is 60 μGal. Therefore, accumulation of error on the open line causes unreliable gravity value on the station at the end of the open line. Consequently, the identification of anomalous gravity change will be seriously affected. In order to improve the reliability of adjustment result effectively, proper data processing method is adopted to convert the open lines into a closed loop, when some long lines exist in gravity network owing to the observation work, helping to identify anomaly information.

The estimation method and accuracy of gravity disturbance in local horizontal coordinate system and inertial coordinate system are compared for the strapdown GNSS/INS combined airborne gravity measurement system. The results show that the accuracy of the Kalman filter method based on the acceleration update in the inertial coordinate system is equivalent to the Kalman filter method based on the position and velocity update in the local horizontal coordinate system. However, the former algorithm is simpler. By using the two methods, the mean standard deviation of the gravity disturbance difference of the six repeatlines are 1.21 mGal and 1.27 mGal respectively, with half-wavelength resolution of 7.5 km. For the same lines, the standard deviation of the gravity disturbance difference calculated by the two methods is 0.65 mGal. The gravity disturbance of the repeatlines obtained by the two methods have evident systematic error, which can be obviously reduced by linear compensation according to time. The maximum difference between the lines is reduced from 8.20 mGal, 8.19 mGal to 2.70 mGal, 2.54 mGal respectively, and the mean internal accuracy is increased from 1.21 mGal, 1.27 mGal to 1.06 mGal, 1.10 mGal respectively.

This paper introduces a method using recorded timing pulses to estimate the time tag accuracy and system characteristic of a seismic digitizer. The estimated time tag accuracy of a digitizer with internal timing error log is consistent with its error log record. The estimated empirical transfer functions in two other experiments are consistent with their frequency responses calculated from FIR filter coefficients provided by the manufacturer,and their filter delay times are compensated. Our experiment’s results show that the non-linearity of time tag accuracy drift without timing synchronization could differ very much between different digitizers, even using the same model; additionally, that there could be about 1 sample point timing difference among different digitizers and different decimation filters due to the inconsistency in their filter delay time compensation. Such timing inconsistency could be unacceptable in work that requires accurate timing service.

Power spectral density analysis of the VP vertical pendulum tiltmeter second sampling data on normal time and interference periods helps to further understand the normal observation background and interference characteristics of observation data, and it provides a more objective reference for identifying the seismic anomalies information and revealing the process of crustal change.

This paper introduces the test and calibration of the LCR-G gravimeter after having been returned to the factory for repair. Comparison of each capability with the observational data of the gravimeter before and after its repair shows that the fine and targeted inspection and calibration technique is very important for the objective evaluation of the instrument's maintenance and performance.

Temperature gradients of several wells, including Baoshan well in Yunnan， Chuan No.03 well in Sichuan, and Jingle well in Shanxi, are carefully measured. The reason for the co-seismic response of the location of main water in the borehole, the aquifer distribution, and water temperature and water level are inverted. The method can even analyze the effect of different lithologies and wellhead casing temperatures on the wellbore.

The main objective of this research is to investigate the response characteristics of water level to atmospheric loading and solid Earth tide in different frequency bands in the Shiquan well, Shaanxi. Using the recorded water level, barometric pressure and the corresponding theoretical gravity tide from 1 November 2015 to 29 February 2016, we obtain the barometric and tidal responses of the water level by fast Fourier transform (FFT), coherence function, as well as cross-correlation. The results show that barometric response is relatively weak and fluctuates greatly, and the tidal response is extremely weak in low frequency band (f＜0.5 cpd); in the intermediate frequency band (0.5－8 cpd), tidal response is excellent; furthermore, the barometric responses are relatively strong at the frequency of K1, S2, and S3 waves, while the barometric and tidal responses are both relatively weak. The coherence functions between the water level and the barometric pressure are relatively low in high frequency (f＞8 cpd), reflecting a bad correlation between the two signals, which may be caused by the low signal-to-noise ratio in the water level and the relatively weak barometric signals. Moreover, the time lag of barometric response increases from 1 minute to 720 minutes with increasing periods of barometric pressure over all frequency bands. Lead effects and anomalous fluctuations appear in barometric response at some frequency points or ranges over medium-low frequency bands.