Journal of Geodesy最新文献

{"title":"An extended w-test for outlier diagnostics in linear models","authors":"Yangkang Yu, Ling Yang, Yunzhong Shen","doi":"10.1007/s00190-024-01855-0","DOIUrl":"https://doi.org/10.1007/s00190-024-01855-0","url":null,"abstract":"<p>The issue of outliers has been a research focus in the field of geodesy. Based on a statistical testing method known as the <i>w</i>-test, data snooping along with its iterative form, iterative data snooping (IDS), is commonly used to diagnose outliers in linear models. However, in the case of multiple outliers, it may suffer from the masking and swamping effects, thereby limiting the detection and identification capabilities. This contribution is to investigate the cause of masking and swamping effects and propose a new method to mitigate these phenomena. First, based on the data division, an extended form of the <i>w</i>-test with its reliability measure is presented, and a theoretical reinterpretation of data snooping and IDS is provided. Then, to alleviate the effects of masking and swamping, a new outlier diagnostic method and its iterative form are proposed, namely data refining and iterative data refining (IDR). In general, if the total observations are initially divided into an inlying set and an outlying set, data snooping can be considered a process of selecting outliers from the inlying set to the outlying set. Conversely, data refining is then a reverse process to transfer inliers from the outlying set to the inlying one. Both theoretical analysis and practical examples show that IDR would keep stronger robustness than IDS due to the alleviation of masking and swamping effect, although it may pose a higher risk of precision loss when dealing with insufficient data.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"13 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141334123","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"IAG newsletter","authors":"Gyula Tóth","doi":"10.1007/s00190-024-01873-y","DOIUrl":"https://doi.org/10.1007/s00190-024-01873-y","url":null,"abstract":"","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"315 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141334418","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Combined algorithms of high-frequency topographical effects for the boundary-value problems based on Helmert's second condensation method","authors":"Jian Ma, Ziqing Wei, Zhenhe Zhai, Xinxing Li","doi":"10.1007/s00190-024-01844-3","DOIUrl":"https://doi.org/10.1007/s00190-024-01844-3","url":null,"abstract":"<p>The Helmert’s second condensation method is usually used to condense the topographical masses outside the boundary surface in the determination of the geoid and quasi-geoid based on the boundary-value theory. The condensation of topographical masses produces direct and indirect topographical effects. Nowadays, the Remove-Compute-Restore (RCR) technique has been widely utilized in the boundary-value problems. In view of spectral consistency, high-frequency direct and indirect topographical effects should be used in the Hotine-Helmert/Stokes–Helmert integral when the Earth gravitational model serves as the reference model in determining the (quasi-) geoid. Thus, the algorithms for high-frequency topographical effects are investigated in this manuscript. First, the prism methods for near-zone direct and indirect topographical effects are derived to improve the accuracies of near-zone effects compared with the traditional surface integral methods. Second, the Molodenskii spectral methods truncated to power <i>H</i>⁴ are put forward for far-zone topographical effects. Next, the \"prism + Molodenskii spectral-spherical harmonic\" combined algorithms for high-frequency topographical effects are further presented. At last, the effectiveness of the combined algorithms for the high-frequency topographical effects are verified in a mountainous test area.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"6 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141315734","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Prospects of GENESIS and Galileo joint orbit and clock determination","authors":"Tomasz Kur, Krzysztof Sośnica, Maciej Kalarus","doi":"10.1007/s00190-024-01869-8","DOIUrl":"https://doi.org/10.1007/s00190-024-01869-8","url":null,"abstract":"<p>The European Space Agency (ESA) is preparing a satellite mission called GENESIS to be launched in 2027 as part of the FutureNAV program. GENESIS co-locates, for the first time, all four space geodetic techniques on one satellite platform. The main objectives of the mission are the realization of the International Terrestrial Reference Frames and the mitigation of biases in geodetic measurements; however, GENESIS will remarkably contribute to the determination of the geodetic parameters. The precise GENESIS orbits will be determined through satellite-to-satellite tracking, employing two GNSS antennas to observe GPS and Galileo satellites in both nadir and zenith directions. In this research, we show results from simulations of GENESIS and Galileo-like constellations with joint orbit and clock determination. We assess the orbit quality of GENESIS based on nadir-only, zenith-only, and combined nadir–zenith GNSS observations. The results prove that GENESIS and Galileo joint orbit and clock determination substantially improves Galileo orbits, satellite clocks, and even ground-based clocks of GNSS receivers tracking Galileo satellites. Although zenith and nadir GNSS antennas favor different orbital planes in terms of the number of collected observations, the mean results for each Galileo orbital plane are improved to a similar extent. The 3D orbit error of Galileo is improved from 27 mm (Galileo-only), 23 mm (Galileo + zenith), 16 mm (Galileo + nadir), to 14 mm (Galileo + zenith + nadir GENESIS observations), i.e., almost by a factor of two in the joint GENESIS + Galileo orbit and clock solutions.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"313 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141264907","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Uncertainties associated with integral-based solutions to geodetic boundary-value problems","authors":"Pavel Novák, Mehdi Eshagh, Martin Pitoňák","doi":"10.1007/s00190-024-01858-x","DOIUrl":"https://doi.org/10.1007/s00190-024-01858-x","url":null,"abstract":"<p>Physical geodesy applies potential theory to study the Earth’s gravitational field in space outside and up to a few km inside the Earth’s mass. Among various tools offered by this theory, boundary-value problems are particularly popular for the transformation or continuation of gravitational field parameters across space. Traditional problems, formulated and solved as early as in the nineteenth century, have been gradually supplemented with new problems, as new observational methods and data are available. In most cases, the emphasis is on formulating a functional relationship involving two functions in 3-D space; the values of one function are searched but unobservable; the values of the other function are observable but with errors. Such mathematical models (observation equations) are referred to as deterministic. Since observed data burdened with observational errors are used for their solutions, the relevant stochastic models must be formulated to provide uncertainties of the estimated parameters against which their quality can be evaluated. This article discusses the boundary-value problems of potential theory formulated for gravitational data currently or in the foreseeable future used by physical geodesy. Their solutions in the form of integral formulas and integral equations are reviewed, practical estimators applicable to numerical solutions of the deterministic models are formulated, and their related stochastic models are introduced. Deterministic and stochastic models represent a complete solution to problems in physical geodesy providing estimates of unknown parameters and their error variances (mean squared errors). On the other hand, analyses of error covariances can reveal problems related to the observed data and/or the design of the mathematical models. Numerical experiments demonstrate the applicability of stochastic models in practice.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"119 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141287231","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Analysis of long-term distributed autonomous orbit determination for BeiDou-3 satellites","authors":"Fengyu Xia, Shanshi Zhou, Ziqiang Li, NaNa Jiang, Xiaogong Hu","doi":"10.1007/s00190-024-01857-y","DOIUrl":"https://doi.org/10.1007/s00190-024-01857-y","url":null,"abstract":"<p>With the support of inter-satellite link technology, GNSS can theoretically achieve the distributed autonomous orbit determination (AOD) function. Traditional AOD operation generally utilizes the forecast ephemeris uploaded by operational control segment (OCS) as the filter reference orbits or to constrain the orbit systematic errors, especially for constellation overall rotation effects in Earth-centered inertial (ECI) coordinate system. To get rid of the dependency on forecast trajectories for saving the OCS workload and also reduce the onboard storage and computation burden, we use a sequential extended Kalman filter to estimate the orbit parameters and consider main perturbation forces acting on satellites in the AOD solution. In particular, for modeling solar radiation pressure (SRP), an empirical prediction function derived by historical SRP estimates is introduced. Using the proposed scheme, the orbit 3D accuracy and user range error (URE) of the first 180-day distributed AOD solution for BeiDou-3 MEOs with precise Earth rotation parameters (ERPs) can reach about 2.10 and 0.43 m, respectively. The constellation rotation errors implied in AOD orbits around the <i>X</i>-, <i>Y</i>- and <i>Z</i>-axis of ECI system are less than 15.0, 11.7 and 15.2 mas, respectively. For real-world AOD scenarios, precise ERP is not available for satellites. With the 180-day prediction ERP, the orbit 3D errors and URE due to the gradually increased UT1-UTC error can be elevated to 14.62 and 2.91 m during our AOD experiments. Result analysis shows if OCS can upload latest prediction ERP at a frequency of once a week, the 180-day distributed AOD is expected to consistently produce real-time orbits preferable to broadcast ephemeris derived by the traditional region L-band tracking network.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"1 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141264856","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"3D large-scale forward modeling of gravitational fields using triangular spherical prisms with polynomial densities in depth","authors":"Fang Ouyang, Long-wei Chen, Leyuan Wu","doi":"10.1007/s00190-024-01863-0","DOIUrl":"https://doi.org/10.1007/s00190-024-01863-0","url":null,"abstract":"<p>To take the sphericity of the Earth into account, tesseroids are often utilized as grid elements in large-scale gravitational forward modeling. However, such elements in a latitude–longitude mesh suffer from degenerating into poorly shaped triangles near poles. Moreover, tesseroids have limited flexibility in describing laterally variable density distributions with irregular boundaries and also face difficulties in achieving completely equivalent division over a spherical surface that may be desired in a gravity inversion. We develop a new method based on triangular spherical prisms (TSPs) for 3D gravitational modeling in spherical coordinates. A TSP is defined by two spherical surfaces of triangular shape, with one of which being the radial projection of the other. Due to the spherical triangular shapes of the upper and lower surfaces, TSPs enjoy more advantages over tesseroids in describing mass density with different lateral resolutions. In addition, such an element also allows subdivisions with nearly equal weights in spherical coordinates. To calculate the gravitational effects of a TSP, we assume the density in each element to be polynomial along radial direction so as to accommodate a complex density environment. Then, we solve the Newton’s volume integral using a mixed Gaussian quadrature method, in which the surface integral over the spherical triangle is calculated using a triangle-based Gaussian quadrature rule via a radial projection that transforms the spherical triangles into linear ones. A 2D adaptive discretization strategy and an extension technique are also combined to improve the accuracy at observation points near the mass sources. The numerical experiments based on spherical shell models show that the proposed method achieves good accuracy from near surface to a satellite height in the case of TSPs with various dimensions and density variations. In comparison with the classical tesseroid-based method, the proposed algorithm enjoys better accuracy and much higher flexibility for density models with laterally irregular shapes. It shows that to achieve the same accuracy, the number of elements required by the proposed method is much less than that of the tesseroid-based method, which substantially speeds up the calculation by more than 2 orders. The application to the tessellated LITHO1.0 model further demonstrates its capability and practicability in realistic situations. The new method offers an attractive tool for gravity forward and inverse problems where the irregular grids are involved.\u0000</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"71 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141287145","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Analysis of the IGS contribution to ITRF2020","authors":"Paul Rebischung, Zuheir Altamimi, Laurent Métivier, Xavier Collilieux, Kevin Gobron, Kristel Chanard","doi":"10.1007/s00190-024-01870-1","DOIUrl":"https://doi.org/10.1007/s00190-024-01870-1","url":null,"abstract":"<p>As its contribution to the latest release of the International Terrestrial Reference Frame, ITRF2020, the International GNSS Service (IGS) provided a 27-year-long series of daily “repro3” terrestrial frame solutions obtained by combining reprocessed solutions from ten Analysis Centers. This contribution represents an improvement over the previous contribution to ITRF2014, not only by the inclusion of more stations with longer and more complete position time series, but also by a general reduction in random and systematic errors. The IGS contribution to ITRF2020 also provided, for the first time, an independent estimate of the terrestrial scale based on the calibration of the Galileo satellite antennas. Despite the various observed improvements, the repro3 station position time series remain affected by a variety of random and systematic errors. This includes spurious periodic variations in several frequency bands, originating mostly from orbit and tide modeling errors, on top of a combination of white and flicker noise, whose origins remain to be precisely understood. These various components should carefully be accounted for when modeling GNSS station position time series and interpreting them in terms of Earth’s surface deformation. The Galileo-based scale of the repro3 solutions is found to be significantly offset (by <span>(+)</span>4.3 mm at epoch 2015.0) and drifting (by <span>(+)</span>0.11 mm/year) from the SLR/VLBI-based scale of ITRF2020. The reasons for this offset and drift remain to be uncovered.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"101 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141251786","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"All-frequency IGS phase clock/bias product combination to improve PPP ambiguity resolution","authors":"Jianghui Geng, Qiang Wen, Guo Chen, Patrick Dumitraschkewitz, Qiyuan Zhang","doi":"10.1007/s00190-024-01865-y","DOIUrl":"https://doi.org/10.1007/s00190-024-01865-y","url":null,"abstract":"<p>Satellite product combination has been a major effort for the International GNSS Service Analysis Center Coordinator to improve the robustness of orbits, clocks and biases over original AC-specific contributions. While the orbit and clock combinations have been well documented, combining phase biases is more of a challenge since they have to be aligned with the clocks precisely to preserve the exactitude of integer ambiguities in precise point positioning (PPP). In the case of dual-frequency signals, frequency-specific phase biases are first translated into an ionosphere-free form to agree with the IGS satellite clocks, and they can then be integrated as integer clocks to facilitate a joint combination. However, regarding multi-frequency phase biases, forming their ionosphere-free counterparts would be cumbersome as they are linearly dependent. We therefore propose a concept of “frequency-specific integer clock” where all third-frequency phase biases are integrated individually with satellite clocks to enable an efficient frequency-wise combination. The resultant combined product will ensure all-frequency PPP ambiguity resolution over any frequency choices and observable combinations. Our combination test based on the GPS/Galileo satellite products from four IGS-ACs in 2020 showed that the mean phase clock/bias consistencies among ACs for all third-frequency signals (i.e., GPS L5, Galileo E6 and E5b) were as high as 10 ps, and the ambiguity fixing rates were all around 95%. Both quantities reached the same levels as those for the baseline frequencies (i.e., GPS L1/L2 and Galileo E1/E5a). The combined products outperformed AC-specific products since outlier contributions were excluded in the combination.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"2013 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141251771","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Determination of the time-variable geopotential by means of orbiting clocks","authors":"Simone Giuliani, Byron D. Tapley, John C. Ries","doi":"10.1007/s00190-024-01868-9","DOIUrl":"https://doi.org/10.1007/s00190-024-01868-9","url":null,"abstract":"<p>Monitoring the time-variable geopotential identifies the mass redistribution across the Earth and reveals, e.g., climate change and availability of water resources. The features of interest are characterized by spatial and temporal scales accessible only through space missions. Among the most important gravity missions are GRACE (2002–2017), its successor GRACE-FO (since 2018), and GOCE (2009–2013), which all sense the Earth’s gravity field via the geopotential derivatives. We investigate the geopotential estimation through frequency comparisons between orbiting clocks by means of the Doppler-canceling technique, describing the clocks’ behavior in the Earth’s gravitational field via Einstein’s general relativity. The novelty of this approach lies in measuring gravity by sensing the geopotential itself. The proof of principle for the measurement is achieved through an innovative mission scenario: for the first time, the observations are collected by a probing clock in LEO. We show gravity solutions obtained by simulating an estimation problem via our proposed architecture. The results suggest that we can conceivably retrieve the geopotential coefficients with accuracy comparable to the GRACE measurement concept by employing clocks with stabilities of order <span>({10}^{-18})</span>. Presently, terrestrial clocks can routinely attain fractional frequency stabilities of <span>({10}^{-18})</span>, whereas spaceborne clocks are still at the <span>({10}^{-15})</span> level. While our findings are promising, further analysis is needed to obtain more realistic indications on the feasibility of an actual mission, whose realization will be possible when clock technology reaches the required performance. The goal is for the technique investigated in this study to become a future staple for gravity field estimation.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"26 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141264820","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}