Journal of Geodesy最新文献

In this contribution, we introduce the ambiguity-resolved (AR) detector and study its distributional characteristics. The AR-detector is a new detector that lies in between the commonly used ambiguity-float (AF) and ambiguity-known (AK) detectors. As the ambiguity vector can seldomly be known completely, usage of the AK-detector is questionable as reliance on its distributional properties will then generally be incorrect. The AR-detector resolves the shortcomings of the AK-detector by treating the ambiguities as unknown integers. We show how the detector improves upon the AF-detector, and we demonstrate that the, for ambiguity-resolved parameter estimation, commonly required extreme success rates can be relaxed for detection, thus showing that improved model validation is also possible with smaller success rates. As such, the AR-detector is designed to work for mixed-integer GNSS models.

In recent years, significant progress has been in ionospheric modeling research through data ingestion and data assimilation from a variety of sources, including ground-based global navigation satellite systems, space-based radio occultation and satellite altimetry (SA). Given the diverse observing geometries, vertical data coverages and intermission biases among different measurements, it is imperative to evaluate their absolute accuracies and estimate systematic biases to determine reasonable weights and error covariances when constructing ionospheric models. This study specifically investigates the disparities among the vertical total electron content (VTEC) derived from SA data of the Jason and Sentinel missions, the integrated VTEC from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) and global ionospheric maps (GIMs). To mitigate the systematic bias resulting from differences in satellite altitudes, the vertical ranges of various VTECs are mapped to a standardized height. The results indicate that the intermission bias of SA-derived VTEC remains relatively stable, with Jason-1 serving as a benchmark for mapping other datasets. The mean bias between COSMIC and SA-derived VTEC is minimal, suggesting good agreement between these two space-based techniques. However, COSMIC and GIM VTEC exhibit remarkable seasonal discrepancies, influenced by the solar activity variations. Moreover, GIMs demonstrate noticeable hemispheric asymmetry and a degradation in accuracy ranging from 0.7 to 1.7 TECU in the ocean-dominant Southern Hemisphere. While space-based observations effectively illustrate phenomena such as the Weddell Sea anomaly and longitudinal ionospheric characteristics, GIMs tend to exhibit a more pronounced mid-latitude electron density enhancement structure.

The atmospheric phase, which is the sum of vertical stratification and turbulent atmospheric phase, is a major challenge currently faced by small baseline subset interferometric synthetic aperture radar (SBAS-InSAR) measurements. Many previous studies have demonstrated that the former can be separated from the interferogram by establishing a functional model between it and the topography. Due to the high variability of the turbulent atmospheric phase (TAP) in the space and time domains, however, the TAP is difficult to model and remove. Recently, many stochastic models have been developed to reduce the influence of the TAP in SBAS-InSAR. To avoid the rank deficient in stochastic model method, we present a correction method using network-based variance estimation, interferogram stacking and ordinary kriging interpolation (NIO). There are three main steps in the proposed algorithm to ensure the accuracy of the correction result: (1) adaptively identify and select sufficient good-quality interferograms that contain less turbulent atmospheric noise to participate in deformation calculation; (2) further select the short temporal baseline interferogram and mask the corresponding deformation location to avoid the effect of deformation; and 3) take advantage of ordinary kriging interpolation to reduce the effects of TAP from the selected good-quality interferograms. The performance of the proposed method has been validated with a set of simulations and real Sentinel-1A SAR data in Southern California, USA.

A new generation of space-borne LiDAR (Light Detection And Ranging) satellite ICESat-2 (Ice, Cloud, and land Elevation Satellite-2) equipped with ATLAS (Advanced Topographic Laser Altimeter System) can perform earth observation. The main problem is to remove the noise photons from the data. The study proposes a main direction-based noise removal algorithm based on three sets of photon-counting LiDAR data. In order to extract the main direction, features in the spatial neighborhood (k) of photons are calculated, most of the initial noise is removed according to the angle between the main direction of photons and the along-track distance direction. Qualitative and quantitative evaluations are employed to validate the proposed algorithm. The obtained results and the performed analysis reveal that the proposed algorithm can process day and night data with different signal-to-noise ratios, while the accuracy of various surface types exceeds 96%. More specifically, the accuracy of the proposed algorithm for night data can reach 97.43%. Based on quantitative evaluations using SPL (Single photon LiDAR), MATLAS, and airborne LiDAR data, the average R, P, and F values are 0.951, 0.959, and 0.954, respectively. Meanwhile, the result of the proposed algorithm is compatible with the ATL03 photons with low, medium, and high confidence, and its accuracy is superior to ATL08 products. The proposed algorithm had fewer parameters and significantly outperformed the Density-Based Spatial Clustering of Applications with Noise (DBSCAN) and the improved local statistical distance algorithm. This algorithm is expected to provide a reference for subsequent photon-counting LiDAR data processing.

While the theory of random isotropic scalar fields on the sphere is well established, it has not been fully extended to the case of vector fields yet. In this contribution, several theoretical results are thus generalized to random isotropic vector fields on the sphere, including an equivalent of the Wiener–Khinchin theorem, which relates the distance-dependent covariance of the field’s components in a particular rotationally invariant basis to the covariance of the vector spherical harmonic coefficients of the field, i.e., its angular power spectrum. A parametric model, based on a stochastic partial differential equation, is proposed to represent the spatial covariance and angular power spectrum of such fields. Such a model is adjusted, with minor modifications, to empirical spatial correlations of the white noise and flicker noise components of 3D displacement time series of ground global navigation satellite system (GNSS) tracking stations. The obtained spatial correlation model may find several applications such as enhanced detection of offsets in GNSS station position time series, enhanced estimation of long-term ground deformation (i.e., station velocities), enhanced isolation of station-specific displacements (i.e., spatial filtering) and more realistic assessment of uncertainties in all GNSS network-based applications (e.g., estimation of crustal strain rates, of glacial isostatic adjustment models or of tectonic plate motion models).

Precise knowledge of geocenter motion, i.e., the relative motion between the Earth’s center of mass (CM) and the center of figure of the Earth’s surface (CF), is crucial to high-stakes geodetic applications such as sea-level rise monitoring with satellite altimetry or the establishment of regional and global mass budgets with satellite gravimetry. The computation of the latest release of the International Terrestrial Reference Frame, ITRF2020, involved the estimation of a field of seasonal motions for a global network of geodetic stations, expressed with respect to CM, as sensed by satellite laser ranging, from which the translational part represents seasonal geocenter motion. This paper presents two different methods to isolate seasonal geocenter motion from the field of ITRF2020 seasonal station motions, among which a new method based on a direct weighted average of seasonal station motions, with station-specific weights chosen so as to provide a better approximation of CF than the standard network shift approach. The ITRF2020 annual geocenter motion model thus obtained is then compared with other recent geodetic and geophysical estimates. Although different sub-groups of estimates with relatively good internal consistency may be identified, the overall scatter of recent geodetic estimates remains at the level of several mm, i.e., close to the amplitude of annual geocenter motion itself. Efforts toward reconciling seasonal geocenter motion estimates therefore still appear necessary. Meanwhile, it would seem safe to assume that seasonal geocenter motion models, in particular those currently used in satellite altimetry and satellite gravimetry, are still uncertain.

Spherical geodetic satellites tracked by satellite laser ranging (SLR) stations provide indispensable scientific products that cannot be replaced by other sources. For studying the time-variable gravity field, two low-degree coefficients C₂₀ and C₃₀ derived from GRACE and GRACE Follow-On missions are replaced by the values derived from SLR tracking of geodetic satellites, such as LAGEOS-1/2, LARES-1/2, Starlette, Stella, and Ajisai. The subset of these satellites is used to derive the geocenter motion which is fundamental in the realization of the origin of the terrestrial reference frames. LAGEOS satellites provide the most accurate standard gravitational product GM of the Earth. In this study, we use the Kaula theorem of gravitational perturbations to find the best possible satellite height, inclination, and eccentricity for a future geodetic satellite to maximize orbit sensitivity in terms of the recovery of low-degree gravity field coefficients, geocenter, and GM. We also derive the common station-satellite visibility-coverability coefficient as a function of the inclination angle and satellite height. We found that the best inclination for a future geodetic satellite is 35°–45° or 135°–145° with a height of about 1500–1700 km to support future GRACE/MAGIC missions with C₂₀ and C₃₀. For a better geocenter recovery and derivation of the standard gravitational product, the preferable height is 2300–3500 km. Unfortunately, none of the existing geodetic satellites has the optimum height and inclination angle for deriving GM, geocenter, and C₂₀ because there are no spherical geodetic satellites at the heights between 1500 (Ajisai and LARES-1) and 5800 km (LAGEOS-1/2, LARES-2).

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.