Description: Multi-constellation GNSS (multi-GNSS) and multi-frequency signals open new prospects for fast ambiguity resolution (AR) of precise point positioning (PPP). Currently, all the BDS and Galileo satellites are capable of transmitting signals on three or more frequencies. In this contribution, we investigate the triple-frequency PPP ambiguity resolution with B1, B2 and B3 observations from BDS satellites and E1, E5a and E5b observations from Galileo satellites and evaluate the contribution of BDS + Galileo combination to triple-frequency PPP AR. The uncalibrated phase delay (UPD) products are estimated based on triple-frequency observations, and the temporal characteristic as well as the residual distributions are analyzed. Our results show that the extra-wide-lane (EWL) and wide-lane (WL) UPDs for BDS and Galileo satellites are both stable during the 30 days and the daily narrow-lane (NL) UPD series are also steady with no obvious fluctuation. The Galileo UPDs exhibit better performance than BDS UPDs due to the high-quality observations. It is also interesting to find that the EWL UPD corrections for all Galileo satellites are very close to the zero. With the precise UPD products, the triple-frequency PPP AR with BDS and Galileo observations was implemented in both static and kinematic modes. Compared to the ambiguity-float solution, the performance can be significantly improved by triple-frequency PPP AR with the positioning accuracy improved by 30–70% in both static and kinematic modes. Moreover, the triple-frequency PPP fixed solutions also present better performance than the dual-frequency PPP fixed solutions in terms of time to the first fix and positioning accuracy, especially for the Galileo-only and BDS + Galileo solutions. And the fusion of multi-GNSS (BDS and Galileo) can further improve the position estimations compared to the single system with more satellites and better spatial geometry.

Description: The latest generation of GPS satellites, termed Block IIF, provides a new L5 signal. Multi-frequency signals open new prospects for precise positioning and fast ambiguity resolution and have become the trend in Global Navigation Satellite System (GNSS) development. However, a new type of inter-frequency clock bias (IFCB), i.e., the difference between the current clock products computed with L1/L2 and the satellite clocks computed with L1/L5, was noticed. Consequently, the L1/L2 clock products cannot be used for L1/L5 precise point positioning (PPP). In order to solve this issue, the IFCB should be estimated with a high accuracy. Datasets collected at 129 globally distributed Multi-GNSS Experiment (MGEX) stations from 2015 are employed to investigate the IFCB. The results indicate that the IFCB is satellite dependent and varies with the relative sun–spacecraft–earth geometry. Other factors, however, may also contribute to the IFCB variations according to the harmonic analysis of the single-day IFCB time series. In addition, the results show that the IFCB exhibits periodic signal with a notable period of 43,080 s and the peak-to-peak amplitude is 0.023–0.269 m. After considering a time lag of 240 s, the average cross-correlation coefficient between the IFCB series of two consecutive days is 0.943, and the prediction accuracy of IFCB is 0.006 m. A triple-frequency PPP model that takes the IFCB into account is proposed. When using 3-h datasets, the positioning accuracy of triple-frequency PPP can be improved by 19, 13 and 21 % compared with the L1/L2-based PPP in the east, north and up directions, respectively.

Description: During 2016–2018, satellite metadata/information including antenna parameters, attitude laws and physical characteristics such as mass, dimensions and optical properties were released for Galileo and QZSS (except for the QZS-1 optical coefficients). These metadata are critical for improving the accuracy of precise orbit and clock determination. In this contribution, we evaluate the benefits of these new metadata to orbit and clock in three aspects: the phase center offsets and variations (PCO and PCV), the yaw-attitude model and solar radiation pressure (SRP) model. The updating of Galileo PCO and PCV corrections, from the values estimated by Deutsches Zentrum für Luft- und Raumfahrt and Deutsches GeoForschungsZentrum to the chamber calibrations disclosed by new metadata, has only a slight influence on Galileo orbits, with overlap differences within only 1 mm. By modeling the yaw attitude of Galileo satellites and QZS-2 spacecraft (SVN J002) according to new published attitude laws, the residuals of ionosphere-free carrier-phase combinations can be obviously decreased in yaw maneuver seasons. With the new attitude models, the 3D overlap RMS in eclipse seasons can be decreased from 12.3 cm, 14.7 cm, 16.8 cm and 34.7 cm to 11.7 cm, 13.4 cm, 15.8 cm and 32.9 cm for Galileo In-Orbit Validation (IOV), Full Operational Capability (FOC), FOC in elliptical orbits (FOCe) and QZS-2 satellites, respectively. By applying the a priori box-wing SRP model with new satellite dimensions and optical coefficients, the 3D overlap RMS are 5.3 cm, 6.2 cm, 5.3 cm and 16.6 cm for Galileo IOV, FOCe, FOC and QZS-2 satellites, with improvements of 11.0%, 14.7%, 14.0% and 13.8% when compared with the updated Extended CODE Orbit Model (ECOM2). The satellite laser ranging (SLR) validation reveals that the a priori box-wing model has smaller mean biases of − 0.4 cm, − 0.4 cm and 0.6 cm for Galileo FOCe, FOC and QZS-2 satellites, while a slightly larger mean bias of − 1.0 cm is observed for Galileo IOV satellites. Moreover, the SLR residual dependencies of Galileo IOV and FOC satellites on the elongation angle almost vanish when the a priori box-wing SRP model is applied. As for satellite clocks, a visible bump appears in the Modified Allan deviation at integration time of 20,000 s for Galileo Passive Hydrogen Maser with ECOM2, while it almost vanishes when the a priori box-wing SRP model and new metadata are applied. The standard deviations of clock overlap can also be significantly reduced by using new metadata.

Description: Azimuthal anisotropy derived from multimode Rayleigh wave tomography in China exhibits depth-dependent variations in Tibet, which can be explained as induced by the Cenozoic India-Eurasian collision. In west Tibet, the E-W fast polarization direction at depths 〈100 km is consistent with the accumulated shear strain in the Tibetan lithosphere, whereas the N-S fast direction at greater depths is aligned with Indian Plate motion. In northeast Tibet, depth-consistent NW-SE directions imply coupled deformation throughout the whole lithosphere, possibly also involving the underlying asthenosphere. Significant anisotropy at depths of 225 km in southeast Tibet reflects sublithospheric deformation induced by northward and eastward lithospheric subduction beneath the Himalaya and Burma, respectively. The multilayer anisotropic surface wave model can explain some features of SKS splitting measurements in Tibet, with differences probably attributable to the limited back azimuthal coverage of most SKS studies in Tibet and the limited horizontal resolution of the surface wave results.

Description: Power-to-Gas (PtG) is a chemical energy storage technology that converts electrical energy into a high-energy density combustible gas. This technology alleviates the contradiction between power supply and demand due to the intermittent electricity production from environment-friendly renewable energy. Hydrogen (H2), produced by electrolysis, can be used to produce synthetic methane (CH4) by reacting with carbon dioxide (CO2) that is captured from carbon emission sources. Subsurface gas storage is one of the most important processes in a PtG system. However, nearly a half of the total stored gas is used as cushion gas to maintain a suitable reservoir pressure, indicating large amounts of CH4 might be wasted and trapped in geological formations. Based on a PtG system, CO2 can be a promising choice as a cushion gas due to its high effective compressibility near its critical conditions. In this paper, a numerical model is established based on the theory of the fluid flow and molecular diffusion to study the role of CO2 as a cushion gas in increasing the gas storage capacity. The accuracy of this model is verified by comparing with that of Curtis M. Oldenburg. However, because of the declining purity of recovered gas induced by the convection and diffusion of two kinds of gases in the same reservoir. The influences of reservoir thickness on the distribution of mixed region are discussed. The results show that thicker reservoir would be a better choice for the geological storage of CH4 with CO2 as a cushion gas.

Description: In recent times increasing numbers of high-rate GNSS stations have been installed around the world and set-up to provide data in real-time. These networks provide a great opportunity to quickly capture surface displacements, which makes them important as potential constituents of earthquake/tsunami monitoring and warning systems. The appropriate GPS real-time data analysis with sufficient accuracy for this purpose is a main focus of the current GNSS research. The objective of this thesis is to develop high-precision GNSS algorithms for better seismological applications. The core research and the contributions of this thesis are summarized as following: With the availability of real-time high-rate GNSS observations and precise satellite orbit and clock products, the interest in the real-time Precise Point Positioning (PPP) technique has greatly increased to construct displacement waveforms and to invert for source parameters of earthquakes in real time. Furthermore, PPP ambiguity resolution approaches, developed in the recent years, overcome the accuracy limitation of the standard PPP float solution and achieve comparable accuracy with relative positioning. In this thesis, we introduce the real-time PPP service system and the key techniques for real-time PPP ambiguity resolution. We assess the performance of the ambiguity-fixed PPP in real-time scenarios and confirm that positioning accuracy in terms of root mean square (RMS) of 1.0~1.5 cm can be achieved in horizontal components. For the 2011 Tohoku-Oki (Japan) and the 2010 El Mayor-Cucapah (Mexico) earthquakes, the displacement waveforms, estimated from ambiguity-fixed PPP and those provided by the accelerometer instrumentation are consistent in the dynamic component within few centimeters. The PPP fixed solution not only can improve the accuracy of coseismic displacements, but also provides a reliable recovery of earthquake magnitude and of the fault slip distribution in real time. We propose an augmented point positioning method for GPS based hazard monitoring, which can achieve fast or even instantaneous precise positioning without relying on data of a specific reference station. The proposed method overcomes the limitations of the currently mostly used GPS processing approaches of relative positioning and global precise point positioning. The advantages of the proposed approach are demonstrated by using GPS data, which was recorded during the 2011 Tohoku-Oki earthquake in Japan. We propose a new approach to quickly capture coseismic displacements with a single GNSS receiver in real-time. The new approach can overcome the convergence problem of precise point positioning (PPP), and also avoids the integration process of the variometric approach. Using the results of the 2011 Tohoku-Oki earthquake, it is demonstrated that the proposed method can provide accurate displacement waveforms and permanent coseismic offsets at an accuracy of few centimeters, and can also reliably recover the moment magnitude and fault slip distribution. We investigate three current existing single-receiver approaches for real-time GNSS seismology, comparing their observation models for equivalence and assessing the impact of main error components. We propose some refinements to the variometric approach and especially consider compensating the geometry error component by using the accurate initial coordinates before the earthquake to eliminate the drift trend in the integrated coseismic displacements. We propose an approach for tightly integrating GPS and strong motion data on raw observation level to increase the quality of the derived displacements. The performance of the proposed approach is demonstrated using 5 Hz high-rate GPS and 200 Hz strong motion data collected during the El Mayor-Cucapah earthquake (Mw 7.2, 4 April, 2010) in Baja California, Mexico. The new approach not only takes advantages of both GPS and strong motion sensors, but also improves the reliability of the displacement by enhancing GPS integer-cycle phase ambiguity resolution, which is very critical for deriving displacements with highest quality. We also explore the use of collocated GPS and seismic sensors for earthquake monitoring and early warning. The GPS and seismic data collected during the 2011 Tohoku-Oki (Japan) and the 2010 El Mayor-Cucapah (Mexico) earthquakes are analyzed by using a tightly-coupled integration. The performance of the integrated results are validated by both time and frequency domain analysis. We detect the P-wave arrival and observe small-scale features of the movement from the integrated results and locate the epicenter. Meanwhile, permanent offsets are extracted from the integrated displacements highly accurately and used for reliable fault slip inversion and magnitude estimation.

Description: The tropospheric horizontal gradients with high spatiotemporal resolutions provide important information to describe the azimuthally asymmetric delays and significantly increase the ability of ground-based GNSS (Global Navigation Satellite Systems) within the field of meteorological studies, like the nowcasting of severe rainfall events. The recent rapid development of multi-GNSS constellations has potential to provide such high-resolution gradients with a significant degree of accuracy. In this study, we develop a multi-GNSS process for the precise retrieval of high-resolution tropospheric gradients. The tropospheric gradients with different temporal resolutions, retrieved from both single-system and multi-GNSS solutions, are validated using independent numerical weather models (NWM) data and water vapor radiometer (WVR) observations. The benefits of multi-GNSS processing for the retrieval of tropospheric gradients, as well as for the improvement of precise positioning, are demonstrated. The multi-GNSS high-resolution gradients agree well with those derived from the NWM and WVR, especially for the fast-changing peaks, which are mostly associated with synoptic fronts. The multi-GNSS gradients behave in a much more stable manner than the single-system estimates, especially in cases of high temporal resolution, benefiting from the increased number of observed satellites and improved observation geometry. The high-resolution multi-GNSS gradients show higher correlation with the NWM and WVR gradients than the low-resolution gradients. Furthermore, the precision of station positions can also be noticeably improved by multi-GNSS fusion, and enhanced results can be achieved if the high-resolution gradient estimation is performed, instead of the commonly used daily gradient estimation in the multi-GNSS data processing.

Description: The multiple global navigation satellite systems (multi-GNSS) bring great opportunity for the real-time retrieval of high-quality zenith tropospheric delay (ZTD), which is a critical quality for atmospheric science and geodetic applications. In this contribution, a multi-GNSS precise point positioning (PPP) ambiguity resolution (AR) analysis approach is developed for real-time tropospheric delay retrieval. To validate the proposed multi-GNSS ZTD estimates, we collected and processed data from 30 Multi-GNSS Experiment (MGEX) stations; the resulting real-time tropospheric products are evaluated by using standard post-processed troposphere products and European Centre for Medium-Range Weather Forecasts analysis (ECMWF) data. An accuracy of 4.5 mm and 7.1 mm relative to the Center for Orbit Determination in Europe (CODE) and U.S. Naval Observatory (USNO) products is achievable for real-time tropospheric delays from multi-GNSS PPP ambiguity resolution after an initialization process of approximately 5 min. Compared to Global Positioning System (GPS) results, the accuracy of retrieved zenith tropospheric delay from multi-GNSS PPP-AR is improved by 16.7% and 31.7% with respect to USNO and CODE final products. The GNSS-derived ZTD time-series exhibits a great agreement with the ECMWF data for a long period of 30 days. The average root mean square (RMS) of the real-time zenith tropospheric delay retrieved from multi-GNSS PPP-AR is 12.5 mm with respect to ECMWF data while the accuracy of GPS-only results is 13.3 mm. Significant improvement is also achieved in terms of the initialization time of the multi-GNSS tropospheric delays, with an improvement of 50.7% compared to GPS-only fixed solutions. All these improvements demonstrate the promising prospects of the multi-GNSS PPP-AR method for time-critical meteorological applications.

Description: Permeability is known as a key factor affecting the effectiveness of gas production from a natural gas hydrate reser-voir. In addition to the lithological effects of the hydrate-bearing sediment itself, the formation and dissociation processes of gas hydrate in the sediments may also have an influence on the permeability. In this study, samples from a gas hydrate reservoir in the Qilian Mountain permafrost (borehole DK-8 and SK-2 in northern Mulicoalfield) were taken for the permeability experiments. Permeabilities were measured before hydrate formation,with hydrate and also after hydrate dissociation. The presence of solid methane hydrate in the pores lowers the permeability depending on hydrate saturation. However an unexpected high permeability decrease was observed after the dissociation of methane hydrate. Six kinds of permeability tests were carried out to detect the reason for formation damage after hydrate dissociation. The results indicate that the fresh water released from the hydratedissociation may cause an activation and following the migration of fine particles which block the pore throats and finally result in a decrease of permeability. Scanning Electron Microscopy (SEM) analysis on the filter papers which attached on the inlet and outlet of the core sample provides visible evidence on fine migration. In our study we present the experimental results of the permeability tests under different conditions and discuss the potential reasons for our observations. One possible explanation for this phenomenon may be that the release of fresh water causes an increase in the double layer thickness at the water mineral surfaces and therefore increase the repulsion forces between rock particles. This process can release small particles which were attached to the surface of bigger sediment grains.

Description: In this study, we evaluated the SAR information extraction of oil spilled at sea and the effectiveness of combining the fuzzy C-means (FCM) and distance regularized level set evolution (DRLSE) models to extract SAR oil-spill information. In light of the inability of this approach to prevent small-noise effects and its poor thin-oil- film extraction performance, we propose a method for extracting oil-spill information that combines threshold data and the DRLSE model. With this method, the initial contour of the oil-spill region is constructed based on the threshold, which overcomes the influence of small noises on the oil extraction, and the extraction of thin-oil-film information is facilitated. Our method demonstrates better oil-extraction precision than the H/A/alpha-Wishart unsupervised classification method and the combined FCM and DRLSE models.