Search

In terrestrial studies, Bouguer gravity data is routinely computed by adopting various numerical schemes, starting from the most basic concept of approximating the actual topography by an infinite Bouguer plate, through adding a planar terrain correction to account for a local/regional terrain geometry, to more advanced schemes that involve the computation of the topographic gravity correction by taking into consideration a gravitational contribution of the whole topography while adopting a spherical (or ellipsoidal) approximation. Moreover, the topographic density information has significantly improved the gravity forward modeling and interpretations, especially in polar regions (by accounting for a density contrast of polar glaciers) and in regions characterized by a complex geological structure. Whereas in geodetic studies (such as a gravimetric geoid modeling) only the gravitational contribution of topographic masses above the geoid is computed and subsequently removed from observed (free-air) gravity data, geophysical studies focusing on interpreting the Earth's inner structure usually require the application of additional stripping gravity corrections that account for known anomalous density structures in order to reveal an unknown (and sought) density structure or density interface. In planetary studies, numerical schemes applied to compile Bouguer gravity maps might differ from terrestrial studies due to two reasons. While in terrestrial studies the topography is defined by physical heights above the geoid surface (i.e., the geoid-referenced topography), in planetary studies the topography is commonly described by geometric heights above the geometric reference surface (i.e., the geometric-referenced topography). Moreover, large parts of a planetary surface have negative heights. This obviously has implications on the computation of the topographic gravity correction and consequently Bouguer gravity data because in this case the application of this correction not only removes the gravitational contribution of a topographic mass surplus, but also compensates for a topographic mass deficit. In this study, we examine numerically possible options of computing the topographic gravity correction and consequently the Bouguer gravity data in planetary applications. In agreement with a theoretical definition of the Bouguer gravity correction, the Bouguer gravity maps compiled based on adopting the geoid-referenced topography are the most relevant. In this case, the application of the topographic gravity correction removes only the gravitational contribution of the topography. Alternative options based on using geometric heights, on the other hand, subtract an additional gravitational signal, spatially closely correlated with the geoidal undulations, that is often attributed to deep mantle density heterogeneities, mantle plumes or other phenomena that are not directly related to a topographic density distribution. [ABSTRACT FROM AUTHOR]

The static Earth’s gravitational field has traditionally been described in geodesy and geophysics by the gravitational potential (geopotential for short), a scalar function of 3-D position. Although not directly observable, geopotential functionals such as its first- and second-order gradients are routinely measured by ground, airborne and/or satellite sensors. In geodesy, these observables are often used for recovery of the static geopotential at some simple reference surface approximating the actual Earth’s surface. A generalized mathematical model is represented by a surface integral equation which originates in solving Dirichlet’s boundary-value problem of the potential theory defined for the harmonic geopotential, spheroidal boundary and globally distributed gradient data. The mathematical model can be used for combining various geopotential gradients without necessity of their re-sampling or prior continuation in space. The model extends the apparatus of integral equations which results from solving boundary-value problems of the potential theory to all geopotential gradients observed by current ground, airborne and satellite sensors. Differences between spherical and spheroidal formulations of integral kernel functions of Green’s kind are investigated. Estimated differences reach relative values at the level of 3% which demonstrates the significance of spheroidal approximation for flattened bodies such as the Earth. The observation model can be used for combined inversion of currently available geopotential gradients while exploring their spectral and stochastic characteristics. The model would be even more relevant to gravitational field modelling of other bodies in space with more pronounced spheroidal geometry than that of the Earth. [ABSTRACT FROM AUTHOR]

Downward continuation of the potential data (gravity or magnetic) helps to interpret contributing sources by transforming the physical signal to their neighborhood. This paper applies three continuation approaches (Landweber's iteration, the direct method and the inverse-matrix approach with Tikhonov's regularization), based on the Poisson integral equation, and four integration schemes to the second vertical derivative of the anomalous potential $$T_{\rm zz}$$ obtained from GOCE data. In the experiments, $$T_{\rm zz}$$ was downward continued for 250 km in the area of Central Europe with special attention to edge effects. For the integration schemes, which use prior information outside the region of interest to reduce edge effects, the best agreement with TIM-r4 was achieved by the inverse-matrix approach with Tikhonov's regularization (RMS = 1.15 eotvos). On the contrary, without prior information the iterative approach performed with RMS = 1.17 eotvos. [ABSTRACT FROM AUTHOR]

The global crustal model CRUST1.0 (refined using additional global datasets of the solid topography, polar ice sheets and geoid) is used in this study to estimate the average densities of major crustal structures. We further use this refined model to compile the gravity field quantities generated by the Earth's crustal structures and to investigate their spatial and spectral characteristics and their correlation with the crustal geometry in context of the gravimetric Moho determination. The analysis shows that the average crustal density is 2,830 kg/m, while it decreases to 2,490 kg/m when including the seawater. The average density of the oceanic crust (without the seawater) is 2,860 kg/m, and the average continental crustal density (including the continental shelves) is 2,790 kg/m. The correlation analysis reveals that the gravity field corrected for major known anomalous crustal density structures has a maximum (absolute) correlation with the Moho geometry. The Moho signature in these gravity data is seen mainly at the long-to-medium wavelengths. At higher frequencies, the Moho signature is weakening due to a noise in gravity data, which is mainly attributed to crustal model uncertainties. The Moho determination thus requires a combination of gravity and seismic data. In global studies, gravimetric methods can help improving seismic results, because (1) large parts of the world are not yet sufficiently covered by seismic surveys and (2) global gravity models have a relatively high accuracy and resolution. In regional and local studies, the gravimetric Moho determination requires either a detailed crustal density model or seismic data (for a combined gravity and seismic data inversion). We also demonstrate that the Earth's long-wavelength gravity spectrum comprises not only the gravitational signal of deep mantle heterogeneities (including the core-mantle boundary zone), but also shallow crustal structures. Consequently, the application of spectral filtering in the gravimetric Moho determination will remove not only the gravitational signal of (unknown) mantle heterogeneities, but also the Moho signature at the long-wavelength gravity spectrum. [ABSTRACT FROM AUTHOR]