The Bouguer-Rudzki-Hotine scheme for geoid computations

<p>The Rudzki inversion gravimetric reduction maps the Earth’s topographic masses inside the geoid in such a way that the inverted masses produce exactly the same potential as the topographic masses on the geoid. In other words, the indirect effect to the geoid is zero so that its computation is not needed. This paper proposes a geoid computation scheme that combines the Bouguer reduction and Rudzki inversion reduction under the spherical approximation and constant density assumption. The proposed computation scheme works with the Bouguer gravity field that is smooth and theoretically legitimate for the harmonic downward continuation. Then the Bouguer potential is compensated by the potential of the inverted masses, ensuring zero indirect effect to the geoid. The direct effect of the Rudzki inversion gravimetric reduction is added to the Bouguer gravity disturbance, resulting in the reduced gravity disturbance for geoid computation. A spherical harmonic reference gravity model is also developed so that the kernel modification/truncation can be applied to the Hotine integral. If the density of the topographic masses becomes available, the effect of density anomalies can be computed separately and added to the geoid computed under the constant density assumption. The combined ellipsoidal effect of the Bouguer and Rudzki inversion reduction should be insignificant because of the canceling effect between them.</p>

A New Moho Depth Model for Fennoscandia with Special Correction for the Glacial Isostatic Effect

In this study, we present a new Moho depth model in Fennoscandia and its surroundings. The model is tailored from data sets of XGM2019e gravitationl field, Earth2014 topography and seismic crustal model CRUST1.0 using the Vening Meinesz-Moritz model based on isostatic theory to a resolution of 1° × 1°. To that end, the refined Bouguer gravity disturbance is determined by reducing the observed field for gravity effect of topography, density heterogeneities related to bathymetry, ice, sediments, and other crustal components. Moreover, stripping of non-isostatic effects of gravity signals from mass anomalies below the crust due to crustal thickening/thinning, thermal expansion of the mantle, Delayed Glacial Isostatic Adjustment (DGIA), i.e., the effect of future GIA, and plate flexure has also been performed. As Fennoscandia is a key area for GIA research, we particularly investigate the DGIA effect on the gravity disturbance and gravimetric Moho depth determination in this area. One may ask whether the DGIA effect is sufficiently well removed in the application of the general non-isostatic effects in such an area, and to answer this question, the Moho depth is determined both with and without specific removal of the DGIA effect prior to non-isostatic effect and Moho depth determinations. The numerical results yield that the RMS difference of the Moho depth from our model HVMD19 vs. the seismic CRUST19 and GRAD09 models are 3.8/4.2 km and 3.7/4.0 km when the above strategy for removing the DGIA effect is/is not applied, respectively, and the mean value differences are 1.2/1.4 km and 0.98/1.4 km, respectively. Hence, our study shows that the specific correction for the DGIA effect on gravity disturbance is slightly significant, resulting in individual changes in the gravimetric Moho depth up to − 1.3 km towards the seismic results. On the other hand, our study shows large discrepancies between gravimetric and seismic Moho models along the Norwegian coastline, which might be due to uncompensated non-isostatic effects caused by tectonic motions.

Sensitivity of Algorithms for Estimating the Gravity Disturbance Vector to Its Model Uncertainty

The work considers the results of filtering and smoothing of the gravity disturbance vector horizontal components and focuses on the sensitivity of these results to the model parameters in the case when the inertial-geodesic method is applied in the framework of a marine survey on a sea vessel.

Research on the Gravity Disturbance Compensation Terminal for High-Precision Position and Orientation System

The Position and Orientation System (POS) is the core device of high-resolution aerial remote sensing systems, which can obtain the real-time object position and collect target attitude information. The goal of exceeding 0.015°/0.003° of its real-time heading/attitude measurement accuracy is unlikely to be achieved without gravity disturbance compensation. In this paper, a high-precision gravity data architecture for gravity disturbance compensation technology is proposed, and a gravity database with accuracy better than 1 mGal is constructed in the test area. Based on the “Block-Time Variation” Markov Model (B-TV-MM), a gravity disturbance compensation device is developed. The gravity disturbance compensation technology is applied to POS products for the first time, and is applied in the field of aerial remote sensing. Flight test results show that the heading accuracy and attitude accuracy of POS products are improved by at least 6% and 16%, respectively. The device can be used for the gravity disturbance compensation of various inertial technology products.

Plate and faults boundary detection using gravity disturbance and Bouguer gravity anomaly from space geodesy

Nowadays, satellite technology has developed significantly. Geodesy satellites such as Grace and Grace-FO can be used for subsurface mapping. The mapping is in the form of detection of the plate details, faults, and regional geodynamic conditions. This study aims to detect plate and faults from space geodesy using the gravity disturbance and Bouguer gravity anomaly parameter. The study area is in the Sunda Strait. Gravity disturbance is one of the gravity model parameters. Gravity disturbance is the gravitational potential of the topography expressed by the spherical harmonic model and the topographic effect by Barthelmes's calculations. Gravity disturbance can visualize subsurface conditions. Bouguer gravity anomaly is needed to get the condition on subsurface objects. This parameter visualizes subsurface conditions in the form of rocks and non-rocks. These conditions can distinguish oceanic crust and continental crust. Gravity contours are needed to obtain plate and faults predictions. The results obtained are validated patterns and shapes with plate and faults secondary data. The tolerance used in this validation is 80%. The gravity disturbance parameter obtained a value of 83% in verifying the accuracy of assessment in plate and faults detection. The Bouguer gravity disturbance parameter obtained a verification value of accuracy assessment in plate detection but 65% in faults detection. This accuracy assessment uses pattern and texture parameters in detecting the similarity of two or more images. This plate and faults detection results are more detailed and can be used for geophysical, geological, earthquake, and earth dynamics applications.

Vertical Deflections and Gravity Disturbances Derived from HY-2A Data

The first Chinese altimetry satellite, Haiyang-2A (HY-2A), which was launched in 2011, has provided a large amount of sea surface heights which can be used to derive marine gravity field. This paper derived the vertical deflections and gravity disturbances using HY-2A observations for the major area of the whole Earth’s ocean from 60°S and 60°N. The results showed that the standard deviations (STD) of vertical deflections differences were 1.1 s and 3.5 s for the north component and the east component between HY-2A’s observations and those from EGM2008 and EIGEN-6C4, respectively. This indicates the accuracy of the east component was poorer than that of the north component. In order to clearly demonstrate contribution of HY-2A’s observations to gravity disturbances, reference models and the commonly used remove-restore method were not adopted in this study. Therefore, the results can be seen as ‘pure’ signals from HY-2A. Assuming the values from EGM2008 were the true values, the accuracy of the gravity disturbances was about −1.1 mGal in terms of mean value of the errors and 8.0 mGal in terms of the STD. This shows systematic errors if only HY-2A observations were used. An index of STD showed that the accuracy of HY-2A was close to the theoretical accuracy according to the vertical deflection products. To verify whether the systematic errors of gravity field were from the long wavelengths, the long-wavelength parts of HY-2A’s gravity disturbance with wavelengths larger than 500 km were replaced by those from EGM2008. By comparing with ‘pure’ HY-2A version of gravity disturbance, the accuracy of the new version products was improved largely. The systematic errors no longer existed and the error STD was reduced to 6.1 mGal.