{"title":"The statistical testing of regularized mathematical models in geodetic data processing","authors":"Artur Fischer, Krzysztof Nowel, Slawomir Cellmer","doi":"10.1007/s00190-025-01934-w","DOIUrl":"https://doi.org/10.1007/s00190-025-01934-w","url":null,"abstract":"<p>The geodetic community commonly challenges the composite hypotheses in the statistical testing of mathematical models. Since the composite hypotheses are not specified as opposed to their simple counterparts, they require <i>a prior</i> estimation of the model parameters. However, if the mathematical models are ill-conditioned, the regularized estimation is often applied for the parameters of interest. Due to the biased property, the regularized estimation does not rigorously originate in the principle of maximum likelihood (ML) estimation, which was the base for developing the theory of the generalized likelihood ratio (GLR) test. Since the regularized estimator of the parameters of interest is consequently inconsistent with the ML one, one cannot construct the GLR test, which is the uniformly most powerful invariant (UMPI) test. So far, only the bias correction approach has been suggested to solve this problem. In this contribution, an implicit representation of the regularized mathematical model is proposed. It eliminates the complete impact of regularized estimation on a mathematical model and delivers the misclosures analytically free from the influence of regularization. Thus, one can construct the GLR test, which belongs to the UMPI family, and then formulate the test statistic in terms of misclosures.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"15 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143393048","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":"Spatially enhanced interpolating vertical adjustment model for precipitable water vapor","authors":"Hao Yang, Vagner Ferreira, Xiufeng He, Wei Zhan, Xiaolei Wang, Shengyue Ji","doi":"10.1007/s00190-025-01936-8","DOIUrl":"https://doi.org/10.1007/s00190-025-01936-8","url":null,"abstract":"<p>As a critical parameter in meteorological monitoring, precipitable water vapor (PWV) is widely used in short-term extreme weather forecasting and long-term climate change research. However, as PWV exhibits significant vertical attenuation, especially within 2 km, achieving accurate vertical interpolation is essential for comparisons and fusion across different measurement techniques, such as sampling water vapor at different heights. PWV vertical adjustment relies only on the empirical or time-varying lapse rate models (e.g., GPWV-H). The non-uniform vertical distribution of PWV and the uncertain variation trend in the low-latitude region still limit the accuracy. To address these issues, we propose the Spatially enhanced Vertical Adjustment Model for PWV (SPWV-H), taking into account the non-uniform distribution in the vertical direction based on the fifth-generation European Centre for Medium-Range Weather Forecasts Atmospheric Reanalysis (ERA5) products. The assessment, validated against the ERA5 benchmark, highlights the SPWV-H model’s superior performance with an RMSE of 1 mm and a bias of 0.03 mm, especially pronounced in the low-latitude region. Compared to global radiosonde datasets, the SPWV-H model achieves notable reductions in RMSE of 12%, 11%, and 18% when evaluated against the EPWV-H, GPWV-H, and GPT3-1 models, respectively. In spatial interpolation, the SPWV-H model achieves an RMSE of 1.22 mm, indicating an improvement of 10%, 9%, and 14% compared to the EPWV-H, GPWV-H, and GPT3-1 models, respectively. Therefore, the SPWV-H model can provide a reliable service for multi-source PWV fusion and real-time PWV monitoring by GNSS.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"142 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2025-02-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143367391","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}
Yongqiang Yuan, Xingxing Li, Hongjie Zheng, Chutian Gao, Xia Yao
{"title":"Incorporating Satellite Laser Ranging observations into BDS analysis: from the perspectives of orbit validation, precise orbit determination, and geodetic parameters estimation","authors":"Yongqiang Yuan, Xingxing Li, Hongjie Zheng, Chutian Gao, Xia Yao","doi":"10.1007/s00190-025-01939-5","DOIUrl":"https://doi.org/10.1007/s00190-025-01939-5","url":null,"abstract":"<p>In February 2023, the International Laser Ranging Service started the tracking of additional medium Earth orbit satellites from the global BeiDou navigation satellite system (BDS) constellation, increasing the total number of tracked BDS satellites to 27. As an optical space geodesy technique, the Satellite Laser Ranging (SLR) provides another important measurement for BDS other than the microwave (L-band) one. Based on three years of data from June 2021 to May 2024, the potential benefits of introducing SLR data into BDS processing and analysis are investigated from three key aspects: orbit validation, precise orbit determination, and geodetic parameters estimation. The independent SLR validations of BDS precise orbit products from four analysis centers show that using the a priori box-wing model for solar radiation pressure (SRP) modeling can achieve superior performance than purely empirical models. The results also indicate the existence of SRP modeling deficiencies for some satellites such as C45 and C46 with Search and Rescue payloads. Given a sparse ground network with 5 stations, the introduction of SLR significantly stabilizes the SRP parameter estimates and improves the orbit accuracy by 44.4%. In terms of geodetic parameter estimation, the scatter of the Z-component geocenter motion can be effectively reduced with the inclusion of SLR data, presenting 10.9–15.3% smaller root mean square (RMS) values during February 2023 and May 2024, depending on the SRP models. In addition, the annual amplitudes of the Z-component geocenter motion are reduced by 7.2–48.2%. The improvement is more pronounced with a limited number of microwave stations, due to the greater strength of SLR observations in geocenter motion estimation. On the other hand, since the SLR observations are unhomogeneously distributed in both space and time, the incorporation of SLR does not evidently enhance the accuracy of Earth rotation parameters, and may even to some extent contaminate the results when the number of microwave stations is limited.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"62 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2025-02-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143367392","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":"Trends in the M $$_2$$ ocean tide observed by satellite altimetry in the presence of systematic errors","authors":"Richard D. Ray, Michael Schindelegger","doi":"10.1007/s00190-025-01935-9","DOIUrl":"https://doi.org/10.1007/s00190-025-01935-9","url":null,"abstract":"<p>Trends in the deep-ocean M<span>(_2)</span> barotropic tide, deduced from nearly three decades of satellite altimetry and recently presented by Opel et al. (Commun Earth Environ 5:261, https://doi.org/10.1038/s43247-024-01432-5, 2024), are here updated with a slightly longer time series and with a focus on potential systematic errors. Tidal changes are very small, of order 0.2 mm/year or less, with a tendency for decreasing amplitudes, which is evidently a response to the ocean’s increasing stratification and an increasing energy loss to baroclinic motion. A variety of systematic errors in the satellite altimeter system potentially corrupt these small trend estimates. The Dynamic Atmosphere Correction (DAC), derived from an ocean model and used for de-aliasing, introduces a spurious trend (exceeding 0.1 mm/year in places) caused by changes in ECMWF atmospheric tides. Both operational and reanalysis atmospheric tides have spurious trends over the altimeter era. Tidally coherent errors in satellite orbits, including from use of inconsistent tidal geocenter models, are more difficult to bound, although differences between two sets of satellite ephemerides are found to reach 0.1 mm/year for M<span>(_2)</span>. Orbit errors are more deleterious for some other constituents, including the annual cycle. Tidal leakage in the “mesoscale correction,” needed here to suppress non-tidal ocean variability, is a known potential problem, and if the leakage changes over time, it impacts ocean-tide trend estimation. Tests show the error is likely small in the open ocean (<span>(<0.04)</span> mm/year) but large in some marginal seas (<span>(>0.2)</span> mm/year). Potential contamination from other altimeter corrections (e.g., ionospheric path delay) is likely negligible for M<span>(_2)</span> but can be difficult to bound.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"1209 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2025-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143125096","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":"Effect of the Earth’s triaxiality on the tide-generating potential","authors":"Sergey Kudryavtsev","doi":"10.1007/s00190-024-01926-2","DOIUrl":"https://doi.org/10.1007/s00190-024-01926-2","url":null,"abstract":"<p>Latest harmonic developments of the Earth tide-generating potential (TGP), e.g., HW95 (Hartmann and Wenzel in Geoph Res Lett 22:3553, 1995), RATGP95 (Roosbeek in Geophys J Int 126:197, 1996), KSM03 (Kudryavtsev in J Geodesy 77:829, 2004), include a number of terms caused by the joint effect of the Earth’s <i>polar</i> flattening (that can be numerically described by the <span>({J}_{2})</span> geopotential coefficient) and the Moon/the Sun gravitational attraction. In the present study, we additionally consider the effect of the Earth’s <i>equatorial</i> flattening due to the Earth’s triaxiality. Explicit analytical expressions for the relevant part of the TGP are derived. New terms of the TGP development due to the Earth’s triaxial figure are found. Amplitudes of nineteen of them exceed the threshold level of 10<sup>–8</sup> m<sup>2</sup>s<sup>−2</sup> used by the modern tidal potential catalogs. Three of the new terms have the frequency sign opposite to that of the Earth rotation. It is not the case for any previously known term of the Earth TGP development. Every term has a new feature that an integer multiplier of the mean local lunar time used in its argument is not equal to the order of the spherical harmonic associated with the term. It necessitates a relevant modification of the standard HW95 format for representing the Earth TGP. The new terms are suggested for including in the current and future tidal potential catalogs.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"56 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2025-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142990032","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 machine learning-based partial ambiguity resolution method for precise positioning in challenging environments","authors":"Zhitao Lyu, Yang Gao","doi":"10.1007/s00190-024-01932-4","DOIUrl":"https://doi.org/10.1007/s00190-024-01932-4","url":null,"abstract":"<p>Partial ambiguity resolution (PAR) has been widely adopted in real-time kinematic (RTK) and precise point positioning with augmentation from continuously operating reference station (PPP-RTK). However, current PAR methods, either in the position domain or the ambiguity domain, suffer from high false alarm and miss detection, particularly in challenging environments with poor satellite geometry and observations contaminated by non-line-of-sight (NLOS) effects, gross errors, biases, and high observation noise. To address these issues, a PAR method based on machine learning is proposed to significantly improve the correct fix rate and positioning accuracy of PAR in challenging environments. This method combines two support vector machine (SVM) classifiers to effectively identify and exclude ambiguities those are contaminated by bias sources from PAR without relying on satellite geometry. The proposed method is validated with three vehicle-based field tests covering open sky, suburban, and dense urban environments, and the results show significant improvements in terms of correct fix rate and positioning accuracy over the traditional PAR method that only utilizes ambiguity covariances. The fix rates achieved with the proposed method are 93.9%, 81.9%, and 93.1% with the three respective field tests, with no wrong fixes, compared to 72.8%, 20.9%, and 16.0% correct fix rates using the traditional method. The positioning error root mean square (RMS) is 0.020 m, 0.035 m, and 0.056 m in the east, north, and up directions for the first field test, 0.027 m, 0.080 m, and 0.126 m for the second field test, and 0.035 m, 0.042 m, and 0.071 m for the third field test. In contrast, only decimeter- to meter-level accuracy was obtainable with these datasets using the traditional method due to the high wrong fix rate. The proposed method provides a correct and fast time-to-first-fix (TTFF) of 3–5 s, even in challenging environments. Overall, the proposed method offers significant improvements in positioning accuracy and ambiguity fix rate with high reliability, making it a promising solution for PAR in challenging environments.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"83 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2025-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142987635","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}
Shuyun Zheng, Hao Zhou, Zhiyu Ma, Xiang Guo, Zhicai Luo
{"title":"On the feasibility of retrieving the temporal gravity field via improved optical clocks","authors":"Shuyun Zheng, Hao Zhou, Zhiyu Ma, Xiang Guo, Zhicai Luo","doi":"10.1007/s00190-024-01930-6","DOIUrl":"https://doi.org/10.1007/s00190-024-01930-6","url":null,"abstract":"<p>The development of optical clocks has experienced significant acceleration in recent years, positioning them as one of the most promising quantum optical sensors for next-generation gravimetric missions (NGGMs). This study investigates the feasibility of retrieving the temporal gravity field via improved optical clocks through a closed-loop simulation. It evaluates optical clock capabilities in temporal gravity field inversion by considering the clock noise characteristics, designing satellite formations, and simulating the performance of optical clocks. The results indicate that optical clocks exhibit higher sensitivity to low-degree gravity field signals. However, when the optical clock noise level drops below 1 × 10<sup>−19</sup><span>(/sqrt{uptau })</span> (τ being the averaging time in seconds) in the satellite-to-ground (SG) mode or below 1 × 10<sup>−20</sup><span>(/sqrt{uptau })</span> in the satellite-to-satellite (SS) mode, atmospheric and oceanic (AO) errors become the dominant source of error. At this noise level, optical clocks can detect time-variable gravity signals up to approximately degree 30. Compared to existing gravity measurement missions such as GRACE-FO, optical clocks exhibit consistent results in detecting signals below degree 20. If the orbital altitude is reduced to 250 km, the performance of optical clocks across all degrees aligns with the results of GRACE-FO. Furthermore, the study reveals that lowering the orbital altitude of satellite-based optical clocks from 485 to 250 km improves results by an average of 33%. Switching from the SS mode to the SG mode results in an average improvement of 51%, while each order-of-magnitude improvement in clock precision enhances results by an average of 60%. In summary, these findings highlight the tremendous potential of optical clock technology in determining Earth’s temporal gravity field and provide crucial technological support for NGGMs.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"4 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2025-01-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142917199","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":"From one-dimensional to three-dimensional: effect of lateral inhomogeneity on tidal gravity and its implications for lithospheric strength","authors":"Zhenyu Wang, Qian Zhao, Zhigang Shao, Wuxing Wang","doi":"10.1007/s00190-024-01907-5","DOIUrl":"https://doi.org/10.1007/s00190-024-01907-5","url":null,"abstract":"<p>Lateral inhomogeneity in the Earth’s mantle affects the tidal response. The current study reformulates the expressions for estimating the lateral inhomogeneity effects on tidal gravity with respect to the unperturbed Earth and supplements some critical derivation process to enhance the methodology. The effects of lateral inhomogeneity are calculated using several real Earth models. By considering the collective contributions of seismic wave velocity disturbances and density disturbance, the global theoretical changes of semidiurnal gravimetric factor are obtained, which vary from − 0.22 to 0.22% compared to those in a layered Earth model, about 1/2 of the ellipticity’s effect. The gravity changes caused by lateral-inhomogeneous disturbances are also computed and turn out to be up to 0.16% compared to the changes caused by tide-generating potential. The current study compares the influences of lateral inhomogeneity with rotation and ocean tide loading. The results indicate that the rotation and ellipticity on tidal gravity are the most dominant factors, the ocean tide loading is the moderate one, and the lateral inhomogeneity in the mantle has the least significant influence. Moreover, an anti-correlation between the effective elastic thickness and gravimetric factor change caused by lateral inhomogeneity is found, implying that it is difficult to generate tidal response at regions with high rigidity. We argue that the gravimetric factor change can be used as an effective indicator for lithospheric strength.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"31 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-12-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142886765","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":"Regional sea level budget around Taiwan and Philippines over 2002‒2021 inferred from GRACE, altimetry, and in-situ hydrographic data","authors":"Wen-Hau Lan, Chi-Ming Lee, Chung-Yen Kuo, Li-Ching Lin, Eko Yuli Handoko","doi":"10.1007/s00190-024-01928-0","DOIUrl":"https://doi.org/10.1007/s00190-024-01928-0","url":null,"abstract":"<p>The regional sea level budget and interannual sea level changes around Taiwan and Philippines are studied using altimetry, GRACE, and <i>in-situ</i> hydrographic data during 1993‒2021. Results show that the average sea level trend around Taiwan and Philippines during 1993–2021 derived from the altimetric data is 3.6 ± 0.2 mm/yr. Over 2002–2021, the study shows closure of sea level budget in the eastern ocean of Taiwan and Philippines within the observed data uncertainties, and the ocean mass accounts for 88%–100% of the observed sea level rise. In contrast, the sea level budget is not closed in the western ocean of Taiwan and Philippines, probably due to the lack of complete coverage by <i>in-situ</i> ocean observing systems. In addition, both regional sea level anomalies and their steric component around Taiwan and Philippines exhibit pronounced interannual and decadal variabilities. The trade wind stress associated with El Niño–Southern Oscillation and Pacific Decadal Oscillation offers a compelling explanation for the interannual and decadal signals of sea level anomalies in the southern ocean of Taiwan, with negative correlations of − 0.78 to − 0.64, indicating that trade wind stress makes a negative contribution to interannual-to-decadal sea level variability. In the northwestern ocean of Taiwan, the sea level variation is strongly influenced by the local monsoon system and shallow bathymetry with an annual amplitude of 90.3 ± 2.9 mm, larger than those in other regions around Taiwan and Philippines, where ocean mass is dominant with a high correlation with the sea level (+ 0.75 to + 0.78).</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"31 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-12-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142858557","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}