Journal of Geodesy最新文献

{"title":"IAG Newsletter","authors":"Gyula Tóth","doi":"10.1007/s00190-024-01909-3","DOIUrl":"https://doi.org/10.1007/s00190-024-01909-3","url":null,"abstract":"","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"18 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142597694","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":"Modification methods of the Stokes’ kernel for determining the (quasi-) geoid with the Remove-Compute-Restore technique","authors":"Jian Ma, Ziqing Wei, Zhenhe Zhai, Duan Li, Changqiang Feng, Xiaogang Liu, Bin Guan","doi":"10.1007/s00190-024-01902-w","DOIUrl":"https://doi.org/10.1007/s00190-024-01902-w","url":null,"abstract":"<p>The geoid and quasi-geoid serve as the reference surfaces of the orthometric and normal height systems, respectively. In order to improve the accuracy of the (quasi-) geoid determined by the Stokes integral with use of the Remove-Compute-Restore (RCR) technique, various modification methods for the spherical Stokes’ kernels, including the spheroidal, cosine-, power-, and Molodensky-modified kernels, are studied in this paper. In addition to the traditional Molodensky-modified Stokes’ kernel, a more effective Molodensky-modified Stokes’ kernel is put forward. A general formula for spectral decomposition of the Stokes integral in the RCR mode is derived, followed by the spectral analysis to reveal the transfer principles of gravity data when using different Stokes’ kernels. The spheroidal and modified Stokes integrals can cause spectral leakage phenomenon, and a method to eliminate spectral leakage is presented based on spectral analysis. The research indicates the low truncation degree of the spheroidal Stokes’ kernel and the low modification degrees of the modified Stokes’ kernel affect the accuracy of the (quasi-) geoid significantly. Quantitative methods for estimating the empirical values of the parameters of the low-degree spheroidal and modified Stokes’ kernels are proposed and the effectiveness of the methods is validated through numerical tests.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"77 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142594802","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":"A new method for global ionospheric real-time modeling integrating ionospheric VTEC short-term forecast results","authors":"Peng Chen, Rong Wang, Yibin Yao, Mingzhu Xiong, Yuchen Zhang, Xinyue Yang","doi":"10.1007/s00190-024-01911-9","DOIUrl":"https://doi.org/10.1007/s00190-024-01911-9","url":null,"abstract":"<p>As an important data source for monitoring the behavior and variations of the ionosphere, the accuracy of current real-time global ionospheric maps (RT-GIMs) in low-latitude regions and oceanic regions is usually poor, and the accuracy during geomagnetic storms is not ideal. Therefore, the ionospheric vertical total electron content (VTEC) short-term forecast results were integrated into the global ionospheric real-time modeling process to improve the accuracy of RT-GIMs. Firstly, the preliminary RT-GIMs were established by constructing a virtual grid and determining the number of ionospheric pierce points in the grid. Then, different strategies were used to determine the virtual VTEC observations and filled the preliminary RT-GIMs. Finally, the filled RT-GIMs were modeled using spherical harmonic expansion and generated the final RT-GIMs, XRTG. On this basis, three ways were selected to evaluate the accuracy of XRTG. The GPS dSTEC (differential slant total electron content) assessment results showed that the performance of XRTG was the closest to that of Centre for Orbit Determination in Europe’s final GIMs (CODG), and it outperformed other RT-GIMs during geomagnetic storm periods and low-latitude regions. Compared with Universitat Politècnica de Catalunya’s RT-GIMs (UADG) with better performance in other RT-GIMs, the maximum decrease in root mean square error (RMSE) of XRTG during the geomagnetic storm period exceeds 25%, and the maximum decrease in the overall average RMSE of the 20 stations in low latitudes exceeds 27%. The Jason-3 VTEC assessment results showed that the accuracy of XRTG was closer to that of UADG and CODG, and the performance of XRTG and UADG in the range of 22° N–22° S was significantly better than that of other RT-GIMs. The consistency between XRTG and Universitat Politècnica de Catalunya’s rapid GIMs, Chinese Academy of Sciences’ final GIMs, and CODG was good, and the VTEC deviations from each post-processing GIMs were mainly concentrated in the range of ± 5 TECU.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"127 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142594801","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":"Power function of $${varvec{F}}-$$ distribution: revisiting its computation and solution for geodetic studies","authors":"Cüneyt Aydin, Özge Güneş","doi":"10.1007/s00190-024-01905-7","DOIUrl":"https://doi.org/10.1007/s00190-024-01905-7","url":null,"abstract":"<p>The power function of <span>(F-)</span> distribution is the complementary cumulative distribution function of the non-central <span>(F-)</span> distribution. It is used to evaluate the power of the test based on the <span>(F)</span> or <span>({chi }^{2}-)</span> distributed statistics. This paper revisits its computation and solution for the non-centrality parameter in geodetic studies and shows that the power function related to these studies can be computed efficiently and with minimal effort. To facilitate this, we introduce a novel standalone algorithm that consistently computes the power of the test, even for large non-centrality parameters (e.g., <span>(&gt;{10}^{5})</span>) and for <span>({chi }^{2})</span>-distribution. The solution of the power function for the non-centrality parameter is typically obtained using standard root finding algorithms, such as the bisection or Newton–Raphson methods. However, they may encounter convergence problems, particularly when the non-centrality parameter increases. We demonstrate that a solution can be readily obtained from a logarithmic form of the power function, ensuring convergence and removing the requirement for a precisely defined initial value. Furthermore, we utilize a few geometric relationships during the iteration to expedite the solution process. As a result, we propose a novel solution algorithm that is highly precise, stable, and at least four times faster than standard algorithms, even for the solution interval of <span>(&lt;{0, 10}^{6}&gt;)</span>. This efficient solution is published online as a web-based application for geodetic detectability studies in addition to the given MATLAB and Python codes.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"29 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-11-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142580251","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":"Assessing tropospheric turbulence impact on VGOS telescope placement in the Indian subcontinent for the estimation of earth orientation parameters","authors":"Arnab Laha, Johannes Böhm, Sigrid Böhm, Matthias Schartner, Hana Krásná, Nagarajan Balasubramanian, Onkar Dikshit","doi":"10.1007/s00190-024-01912-8","DOIUrl":"https://doi.org/10.1007/s00190-024-01912-8","url":null,"abstract":"<p>The complete set of five Earth Orientation Parameters (EOP) can only be estimated accurately using geodetic Very Long Baseline Interferometry (VLBI). Their precision and accuracy depends on network geometry and station-dependent properties. Atmospheric turbulence poses one of the largest error sources for geodetic VLBI, impacting the precision of EOP. Thus, it becomes imperative to consider this factor while choosing the optimal locations for geodetic VLBI. The magnitude of tropospheric turbulence is approximated through the refractive index structure constant, <span>(C_textrm{n}^textrm{2})</span>. In this study, we simulate the optimal locations for geodetic VLBI in India, considering individual tropospheric turbulence parameters per telescope location. The study identifies 14 potential VLBI stations, co-located with GPS stations and homogeneously distributed all over India, and computes the <span>(C_textrm{n})</span> values from zenith wet delay variances over 24 h obtained from GPS data. These locations are simulated in addition to three different reference networks, which show the current and future VLBI Global Observing System (VGOS) networks. Multiple schedules have been generated and simulated for each configuration using VieSched++, and the precision of EOP is compared when constant and station-specific tropospheric turbulence parameters are used. The study shows that, for the investigated networks, southern stations are optimal for polar motion and celestial pole offsets estimation, whereas an eastern station is optimal for UT1−UTC estimation. Furthermore, the study highlights that for reference networks with fewer stations, utilizing station-specific <span>(C_textrm{n})</span> values significantly influences the determination of optimal locations. It further demonstrates how station-specific <span>(C_textrm{n})</span> values impact the positioning of VGOS telescopes in each network for each EOP differently. The findings show that higher <span>(C_textrm{n})</span> values generally lead to a degradation in EOP precision. Geometrically, a station might be at a good location, but if the <span>(C_textrm{n})</span> value is too high, that location is not favorable.\u0000</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"242 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-11-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142580250","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":"A data-driven approach for analyzing contributions of individual loading factors to GNSS-measured bridge displacements","authors":"Xuanyu Qu, Xiaoli Ding, Yong Xia, Wenkun Yu","doi":"10.1007/s00190-024-01913-7","DOIUrl":"https://doi.org/10.1007/s00190-024-01913-7","url":null,"abstract":"<p>A bridge may displace due to various loadings (e.g., thermal (Xia et al. in Struct Control Health Monit 28(7):e2738, 2013), winds (Owen et al. in J Wind Eng Ind Aerodyn 206:104389, 2020), and vehicles (Xu et al. in J Struct Eng 133(1):3–11, 2007)) acting upon the bridge. Identifying the contributions of individual loading factors to the measured bridge displacements is important for understanding the structural health conditions of the bridge. There is however no effective method to quantify the contributions when multiple loadings act simultaneously on a bridge. We propose a new data-driven method, termed random forest (RF)-assisted variational mode decomposition (RF-AVMD), for more effective identification of dominant loading factors and for quantifying the contributions of individual loading factors to the measured bridge displacements. The proposed method is applicable to studying the displacements of any bridge structures and allows for the first time to separate the contributions of individual loadings. The effectiveness of the proposed method is validated using data from Tsing Ma Bridge (TMB), a large suspension bridge in Hong Kong recorded during two consecutive strong typhoons. The results reveal that the transverse displacements of TMB mid-span were controlled by the crosswinds, the longitudinal displacements were dominated by the temperature and winds along the bridge centerline, and the vertical displacements were mainly due to the winds along the bridge centerline, temperature, and traffic flows. Displacement time series that responded to each loading factor was derived. The proposed method provides important new insights into the impacts of individual loadings on the displacements of long-span bridges.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"25 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-11-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142574313","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":"Two methods for spherical harmonic analysis of area mean values over equiangular blocks based on exact spherical harmonic analysis of point values","authors":"Rong Sun, Zhicai Luo","doi":"10.1007/s00190-024-01900-y","DOIUrl":"https://doi.org/10.1007/s00190-024-01900-y","url":null,"abstract":"<p>Currently, the least-square estimation method is the mainstream method for recovering spherical harmonic coefficients from area mean values over equiangular blocks. Since the least-square estimation method involves matrix inversion, it requires great computation power when the maximum degree to be solved is large. In comparison, numerical quadrature methods are faster. Recent numerical quadrature methods designed for spherical harmonic analysis of area mean values over blocks delineated by equiangular and Gaussian grids are both fast and exact for band-limited data. However, for band-limited area mean values over an equiangular grid that has <span>(N)</span> blocks along the colatitude direction and <span>(2N)</span> blocks along the longitude direction, the maximum degree that can be recovered by using current exact numerical quadrature methods is no larger than <span>(N/2-1)</span>. In this study, by using Lagrange’s method for polynomial interpolation, recently proposed numerical quadrature methods that employ the recurrence relations for the integrals of the associated Legendre’s functions are modified into two new methods. By using these methods, the maximum degree of recovered spherical harmonic coefficients is <span>(N-1)</span>. The results show that these newly proposed methods are comparable in computation speed with the current numerical quadrature methods and are comparable in accuracy with the least-square estimation method for both band-limited and aliased data. Moreover, solving linear systems is not necessary for these two new methods. The error characteristics of these two new methods are quite different from those of methods that employ least-square methods. The spherical harmonic coefficients recovered using these new methods can effectively supplement those recovered using least-square methods.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"87 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-11-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142566170","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":"On theoretical and practical analyses of BDS-3/Galileo/GPS all-in-view and PPP time transfer with the consideration of satellite clock biases","authors":"Dong Zhang, Fu Zheng, Meng Wang, Yuanhui Lin, Liangcheng Deng, Zhen Wang, Chuang Shi","doi":"10.1007/s00190-024-01904-8","DOIUrl":"https://doi.org/10.1007/s00190-024-01904-8","url":null,"abstract":"<p>Benefited from the advantage of high precision, wide area and low cost, the time transfer method based on precise point positioning (PPP) has become a popular technique for the remote clock comparisons. Although the time reference to which satellite clocks are referred can be eliminated by the difference between stations, the effect of satellite clock biases on the estimation of receiver clock offset is always ignored for PPP time transfer. Considering the PPP technique is extended from the all-in-view (AV) by the full use of precise carrier phase observables, a method to evaluate the effects of satellite clock biases on AV and PPP time transfer is proposed first. Then, the GPS, Galileo and BDS-3 time transfer results with different international GNSS Service (IGS) precise products are compared to verify the negative effects of satellite clock biases on AV and PPP time transfer. In our experiment, precise orbit and clock products provided by GFZ (German Research Centre for Geosciences), ESA (European Space Agency) and COD (Center for Orbit Determination in Europe) are used to obtain the clock comparison results of nine time links, including three short baselines, three medium baselines and three long baselines. The results show that the effects of satellite clock biases on AV and PPP time transfer are related to the magnitude of satellite clock biases and the baseline distance between stations. After removing the satellites with larger satellite clock biases, we assess the negative effects of satellite clock biases on AV and PPP time transfer for GPS, Galileo and BDS-3, respectively. By using GFZ, COD and ESA precise products for AV time transfer, the inconsistency of GPS and Galileo time transfer results caused by satellite clock biases is below 0.2 ns for long baseline. Due to the larger satellite clock biases for BDS-3, the inconsistency of BDS-3 AV time transfer results caused by satellite clock biases could reach 0.3 ns for medium baseline and even reach 0.9 ns for long baseline. For GPS and Galileo PPP time transfer, the inconsistency of time transfer results is below 0.05 ns for long baseline. However, the inconsistency of BDS-3 PPP time transfer results can only achieve 0.11 ns for medium baseline and 0.3 ns for long baseline. Thus, it is concluded that the satellite clock biases of BDS-3 precise satellite clock products need refining to improve the performance of BDS-3 PPP time transfer.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"4 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142519528","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":"A new approach to improve the Earth's polar motion prediction: on the deconvolution and convolution methods","authors":"CanCan Xu, ChengLi Huang, YongHong Zhou, PengShuo Duan, QiQi Shi, XueQing Xu, LiZhen Lian, XinHao Liao","doi":"10.1007/s00190-024-01890-x","DOIUrl":"https://doi.org/10.1007/s00190-024-01890-x","url":null,"abstract":"<p>Combining the Liouville equations for polar motion (PM) with forecasted geophysical effective angular momentum (EAM) functions can significantly improve the accuracy of Earth's PM predictions. These predictions rely on deconvolution and convolution methods. Deconvolution derives the geodetic EAM function from the PM observations, while convolution uses both the geodetic and geophysical EAM functions to reproduce and predict the PM values. However, there are limitations in existing numerical realisations of deconvolution and convolution that must be addressed. These limitations include low-frequency biases, high-frequency errors, and edge errors, which can negatively impact the accuracy of PM prediction. To overcome these concerns, we develop the Convolution Least Squares (Conv-LS) scheme through a multi-perspective analysis in the frequency domain, PM domain, and EAM domain. By comparing representative approaches for reproducing three different PM series (IERS C01, IERS C04, and IGS) with varying sampling intervals (18.25 days, daily, and 6 h), we demonstrate that the Conv-LS scheme can effectively eliminate the usually present spurious signals and also integrate high-accuracy deconvolution algorithms to reduce reproduced errors further. Compared to the traditional approacsh (using a low-accuracy discrete PM equation for deconvolution and numerical integration for convolution), our new approach (utilising a high-accuracy deconvolution algorithm along with the Conv-LS scheme for convolution) reduces the standard deviations of the residuals' x-component by 31.0%, 60.1%, and 13.7% for C01, C04, and IGS PM series, respectively, while also reducing the y-component by 17.3%, 47.0%, and 14.0%, respectively. These results highlight the superiority of the Conv-LS scheme, leading us to recommend it for practical applications.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"34 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-10-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142489994","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":"An accurate and lightweight calculation for the high degree truncation coefficient via asymptotic expansion with applications to spectral gravity forward modeling","authors":"Linshan Zhong, Hongqing Li, Qiong Wu","doi":"10.1007/s00190-024-01895-6","DOIUrl":"https://doi.org/10.1007/s00190-024-01895-6","url":null,"abstract":"<p>The truncation coefficient is widely utilized in non-global coverage computations of geophysics and geodesy and is always altitude dependent. As the two commonly used calculation methods for truncation coefficients, i.e., the spectral form and the recursive formula, both suffer from decreasing precision caused by high-altitude, leading to slow convergence for the former and numerical instability recursion for the latter. The asymptotic expansion mathematically converges with increasing degree and can precisely compensate for the shortcomings of the two methods. This study introduces asymptotic expansion to accurately compute the truncation coefficient for the spectral gravity forward modeling to a high degree. The evaluation at the whole altitudes and whole integral radii indicates that the proposed method has the following advantages: (i) The calculation precision increases with increasing degree and is altitude independent; (ii) the accurate calculation can be supported by a double-precision format; and (iii) the calculation can be conducted nearly without extra time cost with increasing degree. Generally, asymptotic expansion is used to calculate the high degree truncation coefficients, while the truncation coefficients at low degrees can be calculated using spectral form or recursive formulas in multiprecision format as a supplement; and the available range of asymptotic expansion is provided in the appendix.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"64 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142405051","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}