Journal of Geodesy最新文献

In this paper, we analyze the general errors-in-variables (EIV) model, allowing both the uncertain coefficient matrix and the dispersion matrix to be rank-deficient. We derive the weighted total least-squares (WTLS) solution in the general case and find that with the model consistency condition: (1) If the coefficient matrix is of full column rank, the parameter vector and the residual vector can be uniquely determined independently of the singularity of the dispersion matrix, which naturally extends the Neitzel/Schaffrin rank condition (NSC) in previous work. (2) In the rank-deficient case, the estimable functions and the residual vector can be uniquely determined. As a result, a unified approach for WTLS is provided by using generalized inverse matrices (g-inverses) as a principal tool. This method is unified because it fully considers the generality of the model setup, such as singularity of the dispersion matrix and multicollinearity of the coefficient matrix. It is flexible because it does not require to distinguish different cases before the adjustment. We analyze two examples, including the adjustment of the translation elimination model, where the centralized coordinates for the symmetric transformation are applied, and the unified adjustment, where the higher-dimensional transformation model is explicitly compatible with the lower-dimensional transformation problem.

The resolution and above all the stability of the geodetic reference frames is crucially important when global change, such as the sea level rise is observed. In this context systematic errors are still presenting a significant challenge to the measurement techniques of space geodesy. In order to overcome this unfortunate situation for the satellite laser ranging technique, we have utilized the injection of a mode-locked laser to provide a stable low-noise link between the optical domain, where the measurements are carried out, and the microwave regime in which the station clock is defined. We obtained a considerably enhanced measurement delay stability by 10–20 ps over several days, albeit with some experimental challenges. The implementation of waveform scans required us to revisit the issue of target structure and intensity variation in satellite laser ranging.

To support real-time global navigation satellite systems (GNSS) precise applications, satellite clock corrections need to be precisely estimated at a high-rate update interval, which remains a challenge due to the rapid development of multi-GNSS constellations. In this study, we developed an undifferenced (UD) ambiguity resolution (AR) procedure to improve both the accuracy and computational efficiency for real-time multi-GNSS clock estimation realized by a square root information filter. In the proposed method, UD ambiguities are resolved after correcting the simultaneously estimated uncalibrated phase delays (UPD) and the fixed UD ambiguity parameters are eliminated immediately from the filter, so that the computational burden is significantly reduced. Moreover, based on the linear relationship between double-differenced (DD) and UD ambiguities, we investigated the difference between DD and UD AR in clock estimation. We found that the major reason why DD AR contributes little to the clock estimation while UD AR can speed up the convergence remarkably is that UD AR additionally provides a stable clock datum compared with DD AR. GNSS observations from about 100 globally distributed stations were processed with the proposed method to generate simulated real-time clocks and UPDs for GPS, Galileo, and BDS satellites over a one-month period. The results show that the percentage of wide-lane (WL) UPD residuals within ± 0.25 cycles and narrow-lane (NL) UPD residuals within ± 0.15 cycles are over 97.0% and 90.0%, respectively, which contributes to an ambiguity fixing rate of more than 90% for three systems. The mean daily standard deviation (STD) of the clocks of the UD-fixed solution with respect to Center for Orbit Determination in Europe 30 s final products is 0.021, 0.020, and 0.035 ns for GPS, Galileo, and BDS satellite, respectively, which is improved by 78.1%, 58.3%, and 79.8% compared to the float solution. Benefiting from the removal of fixed ambiguities, the average computation time per epoch was reduced from 3.88 to 1.05 s with a remarkable improvement of 72.9%. The quality of the satellite clock and UPD products was also evaluated by the performance of kinematic precise point positioning (PPP). The results show that fast and reliable multi-GNSS PPP-AR can be achieved with the derived UD-fixed clocks and UPDs, which outperforms that using DD-fixed clock and off-line UPD products with an average improvement of 7.9% and 19.9% in terms of convergence time and positioning accuracy, respectively. Furthermore, we demonstrated the effectiveness of the proposed UD AR method through a 7-day real-time clock estimation experiment.

Accurate orientation of geodetic instruments is fundamental for understanding deformation processes within the Earth's interior. Misalignment can lead to significant errors in data interpretation, affecting various geophysical applications. However, accurate alignment of standalone instruments like seismometers, strainmeters and tiltmeters remains a challenge in field geodesy. While numerous seismic-wave-based orientation methods have been successfully applied to seismometers, they are often inapplicable to tiltmeters due to their high-frequency filtering behavior and the requirement for a neighboring, pre-oriented instrument. In response to these challenges, we propose a novel orientation calibration method for borehole tiltmeters based on maximizing the correlation between recorded tilt data and theoretical tides by adjusting azimuthal angles. Our study encompasses two kinds of borehole tiltmeters and four datasets from three different field sites. Using solid and ocean tides modeling together with local topography and cavity disturbances, we obtain coefficient correlations ranging between 0.831 and 0.963, and 95% confidence intervals of azimuthal angles below 3.3°. The correlation-based method demonstrates robustness across various tidal-signal extraction techniques, including different averaging window sizes and band-pass filters. Moreover, it yields azimuthal results in agreement with direct compass measurements for known orientations, while exhibiting a moderate sensitivity to factors such as ocean tides and site-specific topography for the studied cases. This method appears to be advantageous when direct measurements are either unavailable or challenging, and emerges as an accurate tool for determining borehole tiltmeter orientation. Its potential applicability may extend beyond tiltmeters to other instruments that can also record tidal phenomena, such as strainmeters and broadband seismometers. Additionally, its utility could be extended to environments like the seafloor, in order to refine the precision of azimuthal angle estimation and simplify methods for azimuthal angle determination.

The combination of satellite laser ranging (SLR) observations to various low earth orbit (LEO) satellites can enhance the accuracy and robustness of SLR-derived geodetic parameters, benefiting the realization of the International terrestrial reference frames. Observation stochastic models play a critical role in the integrated processing of SLR observations to multiple LEO satellites. The consideration of precision in heterogeneous SLR observations from various satellites is essential. In this study, we aim to improve the combination of multi-LEO SLR observations for geodetic parameters determination by optimizing the stochastic model using variance component estimation (VCE). We perform weekly estimates of the geodetic parameters, including station coordinates, Earth rotation parameters, and geocenter coordinates (GCC), using three years of SLR observations to seven LEO satellites at different orbits. The satellite-dependent, station-dependent, and satellite–station-dependent variance components are separately estimated through VCE processing to refine the stochastic models. Given the fact that the precision of SLR observations significantly differs in satellites and stations, the multiple LEO combination can be significantly improved with the implementation of VCE. Satellite–station-pair-dependent variance components are more suitable to the SLR VCE and the accuracy of station coordinates, pole coordinates, and length of day can be averagely improved by 8.4, 22.6, and 21.9%, respectively, compared to the equal-weight solution. Our result also indicates that the observation insufficiency for some stations may result in an unreliable VCE estimation, and eventually leads to an accuracy degradation for station coordinates. To overcome this deficiency, we adopt the variance components derived from the monthly solutions to build the stochastic model in the weekly solutions. The application of monthly weights can effectively mitigate the accuracy deterioration of station coordinates, improving the repeatability of the station coordinates by 15.9, 14.6, and 9.2% with respect to the equal-weight solution in E, N, and U components. The global geodetic parameters also benefit from this processing. The import of monthly weight decreases the outliers in the GCC series, especially in the X and Y components.

In the geocenter motion determination using the Global Navigation Satellite Systems (GNSS), satellite clock offsets are usually estimated as white noise process. The correlation between geocenter coordinates (GCC) and the epoch-wise satellite clocks brings inferior GCC estimates, especially for the Z component. In this contribution, satellite clock offsets are described by the polynomial model, and the deviation of the model from the truth is estimated as a random parameter whose process noise is described by the variogram. Based on 3.7 years of BDS, Galileo and GPS observations from 98 global stations, we investigate the impact of the atomic clock model on GCC estimates. After employing the proposed model, the formal errors of GCC-Z component are reduced by 23–46%, 15–31% and 3–9% for BDS, Galileo and GPS, respectively. When the 7-parameter extended empirical CODE orbit model with the a priori box-wing model (BE7) is used, the atomic clock model reduces the correlation of the B_1C parameter and GCC-Z component by 0.28, 0.23 and 0.07 for BDS, Galileo and GPS, respectively. Besides, a mitigation of about 60% is obtained at the 3rd and 5th BDS draconitic harmonics and a mitigation of 55% at the 3rd Galileo draconitic harmonic for the GCC-Z component. The proposed model also contributes to reduce the annual amplitudes of single BDS, Galileo and GPS solutions, improving the agreement with the Satellite Laser Ranging solutions. As an additional verification, the resulting satellite orbits are also improved by satellite clock modeling. When the BE7 model is applied, the day boundary discontinuities of daily orbits are reduced by 3.4–3.6%, and the RMS of orbit differences relative to the ESA precise orbits is reduced by 8.2–8.5% for BDS and Galileo.

We perform a statistical sensitivity analysis on a parametric fit to vertical daily displacement time series of 244 European Permanent GNSS stations, with a focus on linear vertical land motion (VLM), i.e., station velocity. We compare two independent corrections to the raw (uncorrected) observed displacements. The first correction is physical and accounts for non-tidal atmospheric, non-tidal oceanic and hydrological loading displacements, while the second approach is an empirical correction for the common-mode errors. For the uncorrected case, we show that combining power-law and white noise stochastic models with autoregressive models yields adequate noise approximations. With this as a realistic baseline, we report improvement rates of about 14% to 24% in station velocity sensitivity, after corrections are applied. We analyze the choice of the stochastic models in detail and outline potential discrepancies between the GNSS-observed displacements and those predicted by the loading models. Furthermore, we apply restricted maximum likelihood estimation (RMLE), to remove low-frequency noise biases, which yields more reliable velocity uncertainty estimates. RMLE reveals that for a number of stations noise is best modeled by a combination of random walk, flicker noise, and white noise. The sensitivity analysis yields minimum detectable VLM parameters (linear velocities, seasonal periodic motions, and offsets), which are of interest for geophysical applications of GNSS, such as tectonic or hydrological studies.

In this paper, we propose a solution of designing a topside broadcast ionospheric model to enable the future low earth orbit (LEO) navigation augmentation (LEO-NA) services. Considering the lack of global station observations to develop the LEO-NA ionosphere model, we utilize abundant global navigation satellite system (GNSS) data from LEO satellites to determine the topside global broadcast ionospheric delay. This delay can be combined with existing GNSS broadcast ionospheric delay correction models to determine LEO-NA ionospheric delay. First, the performance of the different-order spherical harmonic (SH) model is evaluated in generating a global topside ionospheric map. The results indicate that by increasing the order from 1 to 2, the internal and external accuracy of the model improves significantly. However, increasing the order from 2 to 8 leads to a decrease in accuracy of 0.10 and 0.11 TECU (total electron content unit) for the internal and external root mean square error. Taking into account compatibility with the Beidou global ionospheric delay correction model, limited data capacity in the navigation message, ionospheric model accuracy, and computational efficiency, we select the second-order SH model as the topside ionosphere broadcast model and outline the strategy for calculating broadcast coefficients. Finally, the accuracy of the topside global broadcast ionospheric delay correction model is evaluated during periods of high and low solar activity. The mean values of root mean square in 2009 and 2014 are 1.49 and 1.88 TECU, respectively. The model in 2009 and 2014 can correct for 67.30% and 72.49% of the ionospheric delay, respectively.

Abstract

The strong noise of satellite-based Time-Variable Gravity (TVG) field is often suppressed by applying the averaging filters. However, how to appropriately compromise the data blurring and de-noising remains as a challenge. In our hypothesis, the optimum spatial averaging filter expects to contain averaging kernels that capture the same amount of orbital samples everywhere, to avoid introducing excessive data blurring. To achieve the goal, we take advantages of the spherical convolution and introduce extra spatial constraints into a Gaussian kernel: (1) its half-width radius adapts to the global inhomogeneity of satellite orbit, and (2) the kernel is reshaped as an ellipsoid to adapt to the regional anisotropy. In this way, we designed optimal filters that contain a spatially-Varying non-isotropic Gaussian-based Convolution (VGC) kernel. The VGC-based filter is compared against three most popular filters through real TVG fields and another closed-loop simulation. In both scenarios, VGC-based filters retain more realistic secular trend and seasonal characteristics, in particular at high latitudes. The spatial correlation between the VGC estimates and the simulated ground truth is found to be 0.95 and 0.86 over Greenland and Antarctica, which is found to be 10% better than other tested filters. Temporal correlations with the ground truth are also found to be considerably better than the other filters over 90% of the globally distributed river basin. Besides, the VGC-based filters provide tolerable efficiency (3.5 s per month) and sufficient accuracy (integral error less than 3%). The method can be extended to the next generation gravity mission as well.

Plain Language Summary

Time-Variable Gravity (TVG) fields of the Gravity Recovery and Climate Experiment (GRACE) and its Follow-On mission (GRACE-FO) need proper filtering to suppress the noise before being applied for intended geophysical studies. Existing filters are generally designed in the spectral domain. Though they are numerically efficient, they can hardly treat the noise in fairness, globally. As a result, the TVG fields may get over-smoothed after applying those filters, particularly in regions with high-latitudes. However, it would be mathematically simple to design a filter by applying a spherical convolution, whose kernels can be easily constrained and tuned in the spatial domain. This study introduces filters with spatially-Varying non-isotropic Gaussian-based Convolution kernel (VGC) that is enforced to comply with the spatial distribution of the TVG noise. The proposed filter is found to preserve a finer spatial resolution of TVG fields, and at the same time, to be able to de-noise them at a comparable level as the existing techniques. Geophysical applications that use GRACE-like TVG fields might have benefits from this practical filtering technique.