Gravity Field Models Research Articles

Impacts of Digital Elevation Model Elevation Error on Terrain Gravity Field Calculations: A Case Study in the Wudalianchi Airborne Gravity Gradiometer Test Site, China

Accurate Digital Elevation Models (DEMs) are essential for precise terrain gravity field calculations, which are critical in gravity field modeling, airborne gravimeter and gradiometer calibration, and geophysical inversion. This study evaluates the accuracy of various satellite DEMs by comparing them with a LiDAR DEM at the Wudalianchi test site, a location requiring ultra-accurate terrain gravity fields. Major DEM error sources, particularly those related to vegetation, were identified and corrected using a least squares method that integrates canopy height, vegetation cover, NDVI, and airborne LiDAR DEM data. The impact of DEM vegetation errors on terrain gravity anomalies and gravity gradients was quantified using a partitioned adaptive gravity forward-modeling method at different measurement heights. The results indicate that the TanDEM-X DEM and AW3D30 DEM exhibit the highest vertical accuracy among the satellite DEMs evaluated in the Wudalianchi area. Vegetation significantly affects DEM accuracy, with vegetation-related errors causing an impact of approximately 0.17 mGal (RMS) on surface gravity anomalies. This effect is more pronounced in densely vegetated and volcanic regions. At 100 m above the surface and at an altitude of 1 km, vegetation height affects gravity anomalies by approximately 0.12 mGal and 0.07 mGal, respectively. Additionally, vegetation height impacts the vertical gravity gradient at 100 m above the surface by approximately 4.20 E (RMS), with errors up to 48.84 E over vegetation covered areas. The findings underscore the critical importance of using DEMs with vegetation errors removed for high-precision terrain gravity and gravity gradient modeling, particularly in applications such as airborne gravimeter and gradiometer calibration.

Evaluation of GOCE/GRACE and combined global geopotential models using GNSS/levelling data over Nigeria

AbstractGlobal Geopotential Models (GGMs) provide valuable information about Earth’s gravity field functionals, such as geoid heights and gravity anomalies. However, ground-based datasets are required to validate these GGMs at the regional and local scales. In this study, we validated the accuracy of GGMs by comparing them with ground-based Global Navigation Satellite System (GNSS)/levelling data for the first time in Nigeria. We employed two validation scenarios: with and without considering spectral consistency using the spectral enhancement method (SEM) to incorporate high and very high frequencies of the gravity field spectrum from the combined global gravity field model (XGM2019e_2159) and the residual terrain model (RTM) derived from the Shuttle Radar Topography Mission (SRTM) data, respectively. The results of this evaluation confirmed that the application of SEM improved the assessment of the GGM solutions in an unbiased manner. Integrating XGM2019e_2159 and SRTM data to constrain the high-frequency component of geoid heights in Gravity Field and Steady-State Ocean Circulation Explorer (GOCE)-based GGMs led to an improvement of approximately 10% in reducing the standard deviation (STD) relative to when SEM was not applied. GO_CONS_GCF_2_TIM_R6 at spherical harmonics (SH) of up to degree and order 260 demonstrated the lowest STD when compared to GO_CONS_GCF_2_DIR_R6 and GO_CONS_GCF_2_SPW_R5, with a reduction from 0.380 m without SEM application to 0.342 m with SEM implementation. In addition, four transformation models, namely, linear, four-parameter, five-parameter, and seven-parameter models, were evaluated. The objective is to mitigate the reference system offsets between the GNSS/levelling data and the GGMs and to identify the particular parametric model with the smallest STD across all GGMs. This effort reduced the GGMs misfits to GNSS/levelling to 0.30 m, representing a 15.3% decrease in STD. Notably, the XGM2019e_2159 model provides this improvement.

Planning, emphasis and proposal for implementation of gravimetric control network in Israel

AbstractThis study outlines the planning, emphases, and proposed implementation of a gravimetric control network for Israel. The establishment of a comprehensive and accurate gravimetric control network is crucial for various geodetic applications, including geoid determination, gravity field modelling, vertical datum establishment, and geophysical studies. This study aims to provide a framework for the development and implementation of such a network in Israel. The proposed network is divided into two orders. The first order, also known as the absolute order or zero order, comprises four stations situated in geologically stable and tranquil locations. The second order, referred to as the first order, is established through relative gravimetric measurements. This order will encompass 29 stations positioned in areas that are both stable and secure. The zero and first order stations will be evenly distributed throughout the country to ensure homogeneous coverage. The proposed gravimetric network is in conjunction with the requirements to improve the geoid undulation model of Israel. Through the establishment of a small gravimetric network and the implementation of relative gravimetric measurements, we were able to determine the optimal control point density and establish a well-structured methodology for measurement, analysis, and the adjustment process.

Performance assessment of sentinel-3/6 altimeter data for marine gravity recovery

High-precision sea surface height is crucial for determining the marine gravity field. The Sentinel-3/6 altimetry missions, equipped with SRAL and Poseidon-4 altimeters, provide this essential data. However, there is a lack of comprehensive assessment of the Sentinel-3/6 altimeters for inverting marine gravity anomalies (MGA). In this study, we employ the inverse Venning-Meinsz method to derive nine sets of 1’×1’ MGAs in the South China Sea (SCS) and the Ross Sea (RS). Specifically, MGAs from Sentinel-3A, Sentinel-3B, Sentinel-6 SARM, Sentinel-6 LRM, HY-2A, ICESat-2, and CryoSat-2 are denoted as S3A, S3B, S6S, S6L, H2A, IS2, and CS2, respectively. MGA from the combined HY-2A, ICESat-2, and CryoSat-2 is referred to as HIC, while 3SHIC denotes the MGA from the combination of Sentinel-3/6 SARM, HY2A, ICESat-2, and CryoSat-2. We assess the performance of these MGAs using the EGM08, DTU17, SIO V32.1, and SDUST2021 gravity field models, as well as shipboard gravity across different ocean regions. Among the Sentinel-3/6 MGAs, S3B exhibits the highest accuracy in the SCS, with a root mean square error (RMSE) of 5.277 mGal, followed closely by S3A. Conversely, S3A demonstrates the highest accuracy with an RMSE of 4.635 mGal, followed by S3B in the RS. The inversion accuracy of MGAs from S6S and S6L are comparable, though S6S outperforms S6L in the open sea. The performance of MGAs from Sentinel-3/6 matches or surpasses that of other altimetry missions during the same period. In the SCS, the best-performing MGA is 3SHIC, with an RMSE of 4.585 mGal, closely matching DTU17. However, 3SHIC exhibits superior performance in the RS with an RMSE of 4.263 mGal compared to DTU17 and SDUST 2021. Furthermore, the performance of 3SHIC, which integrates Sentinel-3/6 data, improves that of HIC by 0.74% and 3.37% in the SCS and RS, respectively. These results underscore the contribution of Sentinel-3/6 altimeters to the MGA, particularly in coastal and high-latitude regions. Integration of Sentinel-3/6 data with other altimetry satellites is expected to enhance the spatial resolution and accuracy of the global marine gravity field, especially with the successful establishment of the network of Sentinel-6 in the future.

Residual and Unmodeled Ocean Tide Signal From 20+ Years of GRACE and GRACE‐FO Global Gravity Field Models

AbstractWe analyze remaining ocean tide signal in K/Ka‐band range‐rate (RR) postfit residuals, obtained after estimation of monthly gravity field solutions from 21.5 years of Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow‐On sensor data. Low‐pass filtered and numerically differentiated residuals are assigned to grids and a spectral analysis is performed using Lomb‐Scargle periodograms. We identified enhanced amplitudes at over 30 ocean tide periods. Spectral replicas revealed several tides from sub‐semidiurnal bands. Increased ocean tide amplitudes are located in expected regions, that is, in high‐latitude, coastal and shallow water regions, although some tides also show distinct patterns over the open ocean. While most identified tides are considered during processing, and therefore the amplitudes represent residual signal w.r.t. the ocean tide model, several unmodeled tides were found, including astronomical degree‐3 tides , , , , and radiational and/or compound tides , , , and . The astronomical degree‐3 tides were observed on a global level for the first time a few years ago in altimeter data. We are unaware of any global data‐constrained solutions for the other tides. The amplitude patterns of these tides exhibit similarities to purely hydrodynamic solutions, and altimeter observations (astronomical degree‐3 only). The sensitivity of the satellites to these rather small tidal effects demands their inclusion into the gravity field recovery processing to reduce orbit modeling errors and a possible aliasing. The conducted study shows enormous potential of RR postfit residuals analysis for validating ocean tide models and improving gravity field recovery processing strategies.

An Improved Average Acceleration Approach of Modelling Earth Gravity Field Based on K-Band Range-Rate Observations

The conventional average acceleration approach relies on K-band range observation, containing an unknown bias, which leads to possible degradation of the precision of Earth’s gravity field modelling. It also suffers from correlated errors caused by three-point numerical differentiation. In this study, an improved approach is proposed that makes use of K-band range-rate observations instead and overcoming the influence of correlated errors by introducing a whitening filter. GRACE-Follow On data spanning the period from January 2019 to December 2022 were processed by the proposed approach and a series of time-varying gravity field models was derived, referred to as SSM-AAA-GFO in this paper. This model series is compared comprehensively with three official model series. Results demonstrate that all model series are highly coincident below degree 30 and reflect similar time-varying gravity field signals in both large and small basins. After filtering, SSM-AAA-GFO shows uncertainty, in the form of equivalent water height below 2.5 cm, which is comparable with three official model series. The comparison results confirm the effectiveness of the proposed approach for precisely modelling a time-varying gravity field based on K-band range-rate observations.

Multidimensional Evaluation of Altimetry Marine Gravity Models with Shipborne Gravity Data from a New Platform Marine Gravimeter

With the development of satellite altimetry technology and the application of new altimetry satellites, the accuracy and resolution of altimeter-derived gravity field models have improved over the last decades. Nowadays, they are close enough to shipborne gravimetry. In this paper, multi-source shipborne gravity data in the South China Sea were taken to evaluate the accuracies of two high-precision altimeter-derived marine gravity field models (SS V30.1, DTU17). In these shipborne gravity data, there are dozens of routes’ ship gravimetry data, obtained from the National Geophysical Data Center (NGDC); data were tracked from a marine survey with a commercial marine gravimeter (type KSS31M), and data were tracked from a marine gravimetry campaign that was conducted with a newly developed platform gravimeter (type JMG) in the South China Sea in September 2020. After various data filtering, processing, and calibrations, the shipborne gravity data were validated with crossover points analysis. Then, the processed shipborne data were employed to evaluate the accuracy of the altimeter-derived marine gravity field models. During this procedure, the quality of JMG shipborne gravity data was compared with the results of KSS31M and NGDC data. Analysis and evaluation results show that the crossover points verification accuracies of KSS31M and JMG are 0.70 mGal and 1.61 mGal, which are much better than the accuracy of NGDC, which is larger than 8.0 mGal. In the area where the bathymetry changes slowly, the root mean square error values between altimetry gravity models and KSS31M data are respectively 3.28 mGal and 4.54 mGal, and those of the JMG data are respectively 2.94 mGal and 2.60 mGal. According to the above results, we can conclude that the JMG has the same 1–2 mGal accuracy level as KSS31M and can meet the measurement requirements of marine gravity.

Constellation design and performance of future quantum satellite gravity missions

Temporal aliasing is currently the largest error contributor to time-variable satellite gravity field models. Therefore, the evolution of sensor technologies has to be complemented by strategies to reduce temporal aliasing errors. The most straightforward way to improve temporal aliasing is through extended satellite constellations because they improve the observation geometry and increase the achievable temporal resolution. Therefore, strategies to optimize the design of larger satellite constellations are investigated in this contribution. A complete constellation modeling procedure is presented, starting from primary design variables (such as the required targeted resolutions) and concluding with concrete orbital elements for the individual satellites. In parallel, it is evaluated if improved instrument sensitivities based on quantum technologies (cold atom interferometry) can be fully exploited in the case of larger constellations. For this, future quantum satellite gravity missions adopting the gradiometry concept (similar to the GOCE mission) and the low-low satellite-to-satellite tracking concept (similar to GRACE/-FO) are simulated on optimized constellations with up to 6 satellites/pairs. The retrieval performance of a 6-pair mission in terms of the global equivalent water height RMS can be improved by a factor of roughly 3 compared to an inclined double-pair mission. 3D-gradiometry intrinsically has a better de-aliasing behavior but has extremely high accuracy requirements for the gradiometer (about 10 µEotvos) and the attitude reconstruction to be of any benefit. All simulations show that when incorporating improved sensor technologies, such as future quantum sensing instruments in extended constellations, temporal aliasing will remain the dominant error source by far, up to five orders of magnitude larger than the instrument errors. Therefore, improving sensor technologies has to go hand in hand with larger satellite constellations and improved space–time parameterization strategies to further reduce temporal aliasing effects.Graphical

A New Combination Approach for Gibbs Phenomenon Suppression in Regional Validation of Global Gravity Field Model: A Case Study in North China

A global gravity field model (GGM) is essential to be validated with ground-based or airborne observational data for the accurate application of the GGM at a regional scale. Furthermore, accurately understanding the commission errors between the GGM and observational data are crucial for improving regional gravity fields. Taking the North China region as an example, to circumvent the omission errors, it is necessary to unify the spatial resolutions of the EIGEN-6C4 model and terrestrial gravity observational data to 110 km (determined by the distribution of gravity stations) by employing the spherical harmonic function for the EIGEN-6C4 model and the Slepian basis function for the gravity data, respectively. However, the application of spherical harmonic function expansions in the gravity model results in the Gibbs phenomenon, which may be a primary factor contributing to commission errors and impedes the accurate validation of the EIGEN-6C4 model with terrestrial gravity data. To effectively mitigate this issue, this study proposes a combination approach of window function filtering and regional eigenvalue constraint (based on the Slepian basis). Utilizing the EIGEN-6C4 gravity model to derive the gravity disturbance field at a resolution of 110 km (with spherical harmonic expansion up to the 180th degree and order), the combination approach effectively suppresses over 90% of high-degree (above the 120th degree) Gibbs phenomena. This approach also reduces signal leakage outside the region, thus enhancing the spatial accuracy of the regional gravity disturbance field. A subsequent comparison of the regional gravity disturbance field derived from the true model and terrestrial gravity data in North China indicates excellent consistency, with a root mean squared error (RMSE) of 0.80 mGal. This validation confirms that the combined approach of window function filtering and regional eigenvalue constraints effectively mitigates the Gibbs phenomenon and yields precise regional gravity fields. This approach is anticipated to significantly benefit scientific applications such as improving the accuracy of regional elevation benchmarks and accurately inverting the Earth’s internal structure.

HUST-Grace2024: a new GRACE-only gravity field time series based on more than 20 years of satellite geodesy data and a hybrid processing chain

Abstract. To improve the accuracy of monthly temporal gravity field models for the Gravity Recovery and Climate Experiment (GRACE) and the GRACE Follow-On (GRACE-FO) missions, a new series named HUST-Grace2024 is determined based on the updated L1B datasets (GRACE L1B RL03 and GRACE-FO L1B RL04) and the newest atmosphere and ocean de-aliasing product (AOD1B RL07). Compared to the previous HUST temporal gravity field model releases, we have made the following improvements related to updating the background models and the processing chain: (1) during the satellite onboard events, the inter-satellite pointing angles are calculated to pinpoint outliers in the K-band ranging (KBR) range-rate and accelerometer observations. To exclude outliers, the advisable threshold is 50 mrad for KBR range rates and 20 mrad for accelerations. (2) To relieve the impacts of KBR range-rate noise at different frequencies, a hybrid data-weighting method is proposed. Kinematic empirical parameters are used to reduce the low-frequency noise, while a stochastic model is designed to relieve the impacts of random noise above 10 mHz. (3) A fully populated scale factor matrix is used to improve the quality of accelerometer calibration. Analyses in the spectral and spatial domains are then implemented, which demonstrate that HUST-Grace2024 yields a noticeable reduction of 10 % to 30 % in noise level and retains consistent amplitudes of signal content over 48 river basins compared with the official GRACE and GRACE-FO solutions. These evaluations confirm that our aforementioned efforts lead to a better temporal gravity field series. This data set is identified with the following DOI: https://doi.org/10.5880/ICGEM.2024.001 (Zhou et al., 2024).

Contribution of landmark-tracking data to the estimation of the gravity field and rotation parameters of 311P/PANSTARRS for the incoming Tianwen-2 space mission

311P/PANSTARRS is one of the two targets of the Chinese asteroid exploration mission Tianwen-II. The estimation of the gravity field and rotation parameters is an important part of the mission. Optical landmark observations directly record the relative relationship between the asteroid and probe. For this purpose, we investigate the role of the addition of landmark observations to radio-tracking data in achieving precise orbit determination and deriving the physical parameters of 311P. More specifically, we use a linear re-weighting fusion scheme to aggregate these two heterogeneous data to determine 311P's gravitational field and rotational parameters. Our results show that landmark observations enhance parameter estimation drastically. We conduct extensive experiments to understand how to optimize the use of landmark observations. Besides, we find that the upper bound error on 311P ephemeris for estimating a 2nd-degree gravity field model is 5 km.

A new spherical harmonic approach to residual terrain modeling: a case study in the central European Alps

For decades, the residual terrain model (RTM) concept (Forsberg and Tscherning in J Geophys Res Solid Earth 86(B9):7843–7854, https://doi.org/10.1029/JB086iB09p07843, 1981) has been widely used in regional quasigeoid modeling. In the commonly used remove-compute-restore (RCR) framework, RTM provides a topographic reduction commensurate with the spectral resolution of global geopotential models. This is usually achieved by utilizing a long-wavelength (smooth) topography model known as reference topography. For computation points in valleys this neccessitates a harmonic correction (HC) which has been treated in several publications, but mainly with focus on gravity. The HC for the height anomaly only recently attracted more attention, and so far its relevance has yet to be shown also empirically in a regional case study. In this paper, the residual spherical-harmonic topographic potential (RSHTP) approach is introduced as a new technique and compared with the classic RTM. Both techniques are applied to a test region in the central European Alps including validation of the quasigeoid solutions against ground-truthing data. Hence, the practical feasibility and benefits for quasigeoid computations with the RCR technique are demonstrated. Most notably, the RSHTP avoids explicit HC in the first place, and spectral consistency of the residual topographic potential with global geopotential models is inherently achieved. Although one could conclude that thereby the problem of the HC is finally solved, there remain practical reasons for the classic RTM reduction with HC. In this regard, both intra-method comparison and ground-truthing with GNSS/leveling data confirms that the classic RTM (Forsberg and Tscherning 1981; Forsberg in A study of terrain reductions, density anomalies and geophysical inversion methods in gravity field modeling. Report 355, Department of Geodetic Sciences and Surveying, Ohio State University, Columbus, Ohio, USA, https://earthsciences.osu.edu/sites/earthsciences.osu.edu/files/report-355.pdf, 1984) provides reasonable results also for a high-resolution (degree 2160) RTM, yet neglecting the HC for the height anomaly leads to a systematic bias in deep valleys of up to 10–20 cm.

Improved gravity field estimation by incorporating the Tiangong Space Station and next-generation CAI satellite gravity mission

The next-generation gravity satellite mission equipped with the Cold Atom Interferometry (CAI) gradiometer has great potential for the Earth's gravity field estimation. Deploying a CAI gradiometer on the Chinese Tiangong Space Station launched for long-term Earth science research not only reduces the cost compared to a dual-satellite constellation but also enhances interdisciplinary collaboration in the Earth's gravity field detection. In this study, we conducted gravity gradient-based simulations to assess the contribution of deploying a CAI gradiometer on the Tiangong Space Station to collaboratively observe the Earth's gravity field with a polar-orbit gravity satellite. The simulation results demonstrate that whether utilizing Vyy component, three diagonal components or full components, the derived gravity field models show significant improvements within 100 degree and above 200 degree after incorporating Tiangong Space Station. In particular, the gravity field solution recovered from three diagonal components achieves the best accuracy. In the case of using diagonal components, the collaboration observation scheme effectively reduced the cumulative geoid height error by approximately 5.3cm (300 d/o). In the spatial domain, the incorporation of the Tiangong Space Station primarily impacts the estimated gravity field within the orbital coverage area of the space station, and this effect is particularly pronounced when just employing Vyy component. However, due to the limitation of angular velocity observation inaccuracy associated with the CAI gradiometer in nadir mode, there is no substantial accuracy improvement observed above 200 degree when adding gradient components.

Benefit of classical leveling for geoid-based vertical reference frames

Classically, vertical reference frames were realized as national or continent-wide networks of geopotential differences derived from geodetic leveling, i.e., from the combination of spirit leveling and gravimetry. Those networks are affected by systematic errors in leveling, leading to tilts in the order of decimeter to meter in larger networks. Today, there opens the possibility to establish a worldwide unified vertical reference frame based on a conventional (quasi)geoid model. Such a frame would be accessible through GNSS measurements, i.e., physical heights would be derived by the method of GNSS-leveling. The question arises, whether existing geodetic leveling data are abolished completely for the realization of vertical reference frames, are used for validation purposes only, or whether existing or future geodetic leveling data can still be of use for the realization of vertical reference frames. The question is mainly driven by the high quality of leveled potential differences over short distances. In the following we investigate two approaches for the combination of geopotential numbers from GNSS-leveling and potential differences from geodetic leveling. In the first approach, both data sets are combined in a common network adjustment leading to potential values at the benchmarks of the leveling network. In the second approach, potential differences from geodetic leveling are used as observable for regional gravity field modeling. This leads to a grid of geoid heights based on classical observables like gravity anomalies and now also on leveled potential differences. Based on synthetic data and a realistic stochastic model, we show that incorporating leveled potential differences improves the quality of a continent-wide network of GNSS-heights (approach 1) by about 40% and that formal and empirical errors of a regional geoid model (approach 2) are reduced by about 20% at leveling benchmarks. While these numbers strongly depend on the chosen stochastic model, the results show the benefit of using leveled potential differences for the realization of a modern geoid-based reference frame. Independent of the specific numbers of the improvement, an additional benefit is the consistency (within the error bounds of each observation type) of leveling data with vertical coordinates from GNSS and a conventional geoid model. Even though we focus on geodetic leveling, the methods proposed are independent of the specific technique used to observe potential (or equivalently height) differences and can thus be applied also to other techniques like chronometric or hydrodynamic leveling.

Connecting the Brazilian Vertical System to the International Height Reference Frame by estimating the vertical datum parameters

One of the challenges in the unification of the global vertical datum is to determine the vertical offset between the local vertical datum and the global vertical datum. Since 2015, the scientific community has been working on the unification of the global vertical datum, following the publication of the resolution for the definition and realization of an International Height Reference System. This paper aims to estimate the vertical offset in Brazil (Imbituba and Santana Data) by combining a regional gravity field model, using the Geodetic Boundary Value Problem approach applying the Remove-Compute-Restore procedure, and GNSS/leveling data. In the Imbituba vertical datum 1271 stations were used and 66 stations were connected to the Santana tide gauge. The results showed that the geopotential value of Imbituba was 62,636,849.87 ± 0.224 m2s-2, and the vertical offset concerning the reference surface (W0) was 0.358 ± 0.023 m. In Santana, the geopotential value was 62,636,836.88 ± 4.108 m2s-2, and the vertical offset was 1.689 ± 0.420 m. The offset between Imbituba and Santana was computed once the physical connection between Brazil's two local vertical data is a challenge. The estimated value was 1.331 ± 0.421 m. This means the Santana vertical datum is 1.331 m higher than the Imbituba vertical datum.

Assessment of Callisto Gravity-field Determination Using the Inter-satellite Range-rate Link

China will launch the “Tianwen-IV” mission around 2030, focusing on the orbiting exploration of Jupiter and Callisto, a moon of Jupiter. As part of this ambitious mission, a main satellite will carry another satellite that will be released in the Jupiter system to continue its journey toward Uranus. Considering the current mission planning, we propose an inter-satellite radio-observation mode that differs from the conventional observation mode of tracking from Earth to precisely determine the orbit of the satellites. Given the significance of the Callisto gravity field model in both science objectives and satellite navigation, we have conducted a series of simulation experiments to evaluate the potential of this inter-satellite range-rate data for accurately estimating the Callisto gravity field. The results obtained from the analysis demonstrate that by utilizing 40 days of ground station observations, it is possible to estimate the gravity field model of Callisto up to a degree of 70. Remarkably, when combining these ground station observations with inter-satellite observations, a comparable level of accuracy can be achieved with just 10 days of observations. Furthermore, with reduced inter-satellite observation noise, accuracy improves, enabling estimation up to 80 degrees or higher. Initial inter-satellite distance selection impacts estimation accuracy. These findings serve as a valuable test bed for the future “Tianwen-IV” mission to perform precise orbit determination and gravity field model estimation to reduce reliance on deep space stations.

A Direct Approach for Local Quasi-Geoid Modeling Based on Spherical Radial Basis Functions Using a Noisy Satellite-Only Global Gravity Field Model

A Direct Approach for Local Quasi-Geoid Modeling Based on Spherical Radial Basis Functions Using a Noisy Satellite-Only Global Gravity Field Model

Gravity field modeling in mountainous areas based on band-limited SRBFs

Gravity field modeling in mountainous areas based on band-limited SRBFs

Orbit Determination of High Orbit Spacecraft based on BDS-3 using Extended Kalman Filter

The positioning of highly elliptical orbit (HEO) spacecraft is very important in deep space exploration and other missions. The Beidou Navigation Satellite System (BDS-3), which was completed in 2020, has put forward new ideas for using GNSS technology to make up for the lack of earth-based TT&C systems. At the same time, the selection of the gravity field order and the perturbation forces must be considered when using dynamical methods for HEO spacecraft positioning. In this paper, simulation experiments are carried out based on the Two-Line Orbital Element (TLE) data of BDS. Two HEO satellites with different semimajor-axis (67509.6 km, 212857.3 km) are selected as the target spacecraft orbits. The single point positioning (SPP) method and Extended Kalman Filter (EKF) method are used for HEO spacecrafts positioning respectively. The experiments show that a 20×20 order gravity field model can be chosen for the orbital integral of HEO spacecraft orbits, and the integration error caused by solar gravity and lunar gravity should be taken into account with the change of orbital altitude. Adding Geostationary Orbit (GEO) and Inclined GeoSynchronous Orbit (IGSO) navigation satellites improved the visibility of these two orbits by 19.45% and 16.64%, respectively, and improved the GDOP value by 14.1% and 3.8%, respectively, compared to observing MEO satellites only. The use of the EKF method can effectively improve the positioning accuracy of the HEO spacecraft compared to the SPP method. For the spacecraft in HEO-1, the RMS value of the positioning errors of the EKF method was 67.82 m, and the accuracy of each direction was improved by 44.28%, 38.85%, and 38.62%, respectively. For the spacecraft in HEO-2, the positioning errors in the x direction of EKF method was within 12 km and within 2 km in the y and z directions. The accuracy of these three directions was improved by 84.91%, 79.24%, and 58.64%, respectively.

Improved GRACE‐FO Gravity Field Solution by Combining Different Accelerometer Transplant Products

AbstractGravity field solutions determined from the Gravity Recovery and Climate Experiment Follow‐on (GRACE‐FO) are affected by the unexpected performance degradation of one accelerometer onboard, which is circumvented by the accelerometer data transplant (ACT) technique. Three operational ACT products are presently available: JPL‐ACT and JPL‐ACH from NASA Jet Propulsion Laboratory (JPL) and TUG‐ACT from Graz University of Technology (TUG). Considering the heterogeneous qualities of individual ACT products, we propose to combine them to improve the gravity field models after a comprehensive assessment of individual products using a consistent Level‐1B data processing. The combination is carried out on the normal equation of gravity field solutions. Spectral, temporal, and spatial evaluations over the period October 2018 to December 2022 demonstrate that the combined solution significantly reduces the noise relative to solutions using individual ACT products, and the average noise reduction rates over JPL‐ACT range from 5% to 12% with varying post‐processing filters and C20 and C30 coefficient treatments, which are double those of the second‐best JPL‐ACH. Moreover, the C20 estimate exhibits significant improvement in the combined solution with its deviation from the satellite laser ranging (SLR) solution showing a root mean square (RMS) value of 7.6 × 10−11, compared to similar RMS values of 15.9 × 10−11, 12.1 × 10−11, and 14.9 × 10−11 for JPL‐ACT, JPL‐ACH, and TUG‐ACT, respectively. This improvement helps greatly reduce the amplitude of the 161‐day spurious signal over polar regions. For C30, the JPL‐ACT gives degraded estimates with an RMS of differences to SLR of 6.2 × 10−11, while the differences are greatly reduced to 2.6 × 10−11, 2.7 × 10−11 and 3.6 × 10−11 for the combined solution, JPL‐ACH and TUG‐ACT, respectively, approaching the formal error of SLR estimates.