Search

Integral transformations of gravitational gradients onto a Gravity Recovery And Climate Experiment (GRACE) type of observable are derived in this article. The gravitational gradients represent components of the gravitational tensor in the local north-oriented frame. The GRACE type of observable corresponds to a difference between two gravitational vectors as projected onto the line of sight between the two GRACE satellites. In total, three integral transformations relating vertical-vertical, vertical-horizontal and horizontal-horizontal gravitational gradients with the GRACE type of observable are provided. Spectral and closed forms of corresponding isotropic kernels are derived for each transformation. Special cases show that the integral transformations are general and relate gravitational gradients to many other quantities of the gravitational field, such as the gravitational vector, and its radial and tangential components. Correctness of the mathematical derivations is validated in a closed-loop simulation using synthetic data. [ABSTRACT FROM AUTHOR]

We derive expressions for computing the gravitational field (potential and its radial derivative) generated by an arbitrary homogeneous or laterally varying density contrast layer with a variable depth and thickness based on methods for a spherical harmonic analysis and synthesis of gravity field. The newly derived expressions are utilised in the gravimetric forward modelling of major known density structures within the Earth's crust (excluding the ocean density contrast) beneath the geoid surface. The gravitational field quantities due to the sediments and crust components density contrasts, shown in numerical examples, are computed using the 2 × 2 arc-deg discrete data from the global crustal model CRUST2.0. These density contrasts are defined relative to the adopted value of the reference crustal density of 2670 kgm. All computations are realised globally on a 1 × 1 arc-deg geographical grid at the Earth's surface. The maxima of the gravitational signal due to the sediments density contrast are mainly along continental shelf regions with the largest sedimentary deposits. The corresponding maxima due to the consolidated crust components density contrast are over areas of the largest continental crustal thickness with variable geological structure. [ABSTRACT FROM AUTHOR]

We developed and applied a novel numerical scheme for a gravimetric forward modelling of the Earth's crustal density structures based entirely on methods for a spherical analysis and synthesis of the gravitational field. This numerical scheme utilises expressions for the gravitational potentials and their radial derivatives generated by the homogeneous or laterally varying mass density layers with a variable height/depth and thickness given in terms of spherical harmonics. We used these expressions to compute globally the complete crust-corrected Earth's gravity field and its contribution generated by the Earth's crust. The gravimetric forward modelling of large known mass density structures within the Earth's crust is realised by using global models of the Earth's gravity field (EGM2008), topography/bathymetry (DTM2006.0), continental ice-thickness (ICE-5G), and crustal density structures (CRUST2.0). The crust-corrected gravity field is obtained after modelling and subtracting the gravitational contribution of the Earth's crust from the EGM2008 gravity data. These refined gravity data mainly comprise information on the Moho interface and mantle lithosphere. Numerical results also reveal that the gravitational contribution of the Earth's crust varies globally from 1,843 to 12,010 mGal. This gravitational signal is strongly correlated with the crustal thickness with its maxima in mountainous regions (Himalayas, Tibetan Plateau and Andes) with the presence of large isostatic compensation. The corresponding minima over the open oceans are due to the thin and heavier oceanic crust. [ABSTRACT FROM AUTHOR]

In geophysical studies investigating the lithosphere structure, the gravitational field generated by the ocean density contrast (i.e., bathymetry-generated gravitational field) represents a significant amount of the signal to be modelled and subsequently removed from the Earth's gravity field. The ocean density contrast is typically calculated as the difference between the mean density values of the Earth's crust and seawater. The approximation of the actual seawater density distribution by its mean value yields relative errors up to about 2% in computed quantities of the gravitational field. To reduce these errors, a more realistic model of the seawater density distribution is utilized based on the analysis of existing oceanographic data of salinity, temperature, and pressure (depth). We study the accuracy of the bathymetry-generated gravitational field quantities formulated for a depth-dependent model of the seawater density distribution. This density distribution approximates the seawater density variations due to an increasing pressure with depth, whereas smaller lateral density variations caused by salinity, temperature, and other oceanographic factors are not taken into consideration. The error analysis reveals that the approximation of the seawater density by the depth-dependent density model reduces the maximum errors to less than 0.6%. The corresponding depth-averaged errors are below 0.1%. The depth-dependent seawater density model is further facilitated in expressions for computing the bathymetry-generated gravitational field quantities by means of the spherical bathymetric (ocean bottom depth) functions. The numerical realization reveals large differences in the results obtained with and without consideration of the depth-dependent seawater density distribution. The maxima of absolute differences reach 201 m/s and 16.5 mGal in computed values of the potential and attraction, respectively. The application of the depth-dependent seawater density model thus significantly improves the accuracy in the forward modelling of the bathymetric gravitational field quantities. [ABSTRACT FROM AUTHOR]

spectral representation of the topographic corrections to gravity field quantities is formulated in terms of spherical height functions. When computing the far-zone contributions to the topographic corrections, various types of the truncation coefficients are applied to a spectral representation of Newton's integral. In this study we utilise Molodensky's truncation coefficients in deriving the expressions for the far-zone contributions to topographic corrections. The expressions for computing the far-zone gravity field contributions corrected for the effect of topography are then obtained by combining the expressions for the far-zone contributions to the gravity field quantities and to the respective topographic corrections, both expressed in terms of Molodensky's truncation coefficients. The numerical examples of the far-zone contributions to the topographic corrections and to the topography-corrected gravity field quantities are given over the study area situated in the Canadian Rocky Mountains with adjacent planes. Coefficients of the global elevation and geopotential models are used as the input data. [ABSTRACT FROM AUTHOR]

The geoid is an equipotential surface of the static Earth's gravity field which plays a fundamental role in definition of physical heights related to the mean sea level (orthometric heights) in geodesy and which represents a reference surface in many geoscientific studies. Its determination with the cm-level accuracy or better, in particular over dry land, belongs to major tasks of modern geodesy. Traditional data and underlined theory have significantly been affected in recent years by rapid advances in observation techniques. This study reviews gradients of the disturbing gravity potential, both currently available and foreseen, and systematically discusses mathematical models for geoid determination based on gradient data. Fundamentals required for geoid definition and its estimation from measured potential gradients are shortly reviewed at the beginning of the text. Then particular mathematical models based on solutions to boundary-value problems of the potential theory, which include both integral transforms and integral equations, are formulated. Properties of respective integral kernel functions are demonstrated and discussed. With the new mathematical models introduced, new research topics are opened which must be resolved in order to allow for their full-fledged applicability in geoid modelling. Stochastic modelling is also discussed which estimates gradient spatial resolution and accuracy required for geoid modelling with the cm-level accuracy. Results of stochastic modelling suggest that the cm-geoid can be estimated using available gradient data if related problems, namely reduction of gradient data for gravitational effects of all masses outside the geoid and their downward continuation, are solved at the same level of accuracy. [ABSTRACT FROM AUTHOR]

The uncertainties of the geoidal heights estimated from ground gravity data caused by their spatial distribution and noise are investigated in this study. To test these effects, the geoidal heights are estimated from synthetic ground gravity data using the Stokes- Helmert approach. Five different magnitudes of the random noise in ground gravity data and three types of their spatial distribution are considered in the study, namely grid, semigrid and random. The noise propagation is estimated for the two major computational steps of the Stokes-Helmert approach, i.e., the downward continuation of ground gravity and Stokes's integration. Numerical results show that in order to achieve the 1-cm geoid, the ground gravity data should be distributed on the grid or semi-grid with the average angular distance less than 2′. If they are randomly distributed (scattered gravity points), the 1-cm geoid cannot be estimated if the average angular distance between scattered gravity points is larger than 1′. Besides, the noise of the gravity data for the tree types of their spatial distribution should be below 1 mGal to estimate the 1-cm geoid. The advantage of interpolating scattered gravity points onto the regular grid, rather than using them directly, is also investigated in this study. Numerical test shows that it is always worth interpolating the scattered points to the regular grid except if the scattered gravity points are sparser than 5′. [ABSTRACT FROM AUTHOR]

This article discusses methods for selecting criteria used for stopping iterative regularization methods applied to solving continuation of gravity observed at or outside the Earth's surface down to the geoid. This computational step is required for determination of the gravimetric geoid by the Stokes approach. For surface gravity data measured with spatial resolutions at the kilometer level and higher, their downward continuation becomes ill posed and its solution numerically unstable. Two iterative methods often used for solving this problem, namely the Landweber and conjugate gradient least-squares (CGLS) methods, are investigated using a sample of surface gravity data synthesized from a global gravitational model over a mountainous test area in Colorado, USA. Five different commonly used stopping criteria are applied to both iterative methods and neither of them, except L-curve, succeed to stop the iterative process at the reasonable solution. A simple approach is proposed to stop the iterative methods exploiting stochastic properties of surface gravity data. Results of numerical tests suggest that a rather simple stopping criterion based on stochastics parameters of surface gravity data estimated using a leave-one-out cross-validation procedure yields results superior to those based on the L-curve approach which is best performing among the five standard methods for selecting the stopping criterion. Results based on the Landweber method outperform those based on the CGLS method. Iterative methods based on optimum numbers of iterations represent an alternative to commonly used regularization techniques for solving ill-posed inverse problems, such as Tikhonov, that require optimal estimates of regularization parameters. [ABSTRACT FROM AUTHOR]

The geoid-to-quasigeoid separation is often computed only approximately as a function of the simple planar Bouguer gravity anomaly and the height of the computation point while disregarding the contributions of terrain geometry and anomalous topographic density as well as the sub-geoid masses. In this study we demonstrate that these contributions are significant and, therefore, should be taken into consideration when investigating the relation between the normal and orthometric heights particularly in the mountainous, polar and geologically complex regions. These contributions are evaluated by applying the spectral expressions for gravimetric forward modelling and using the EIGEN-6C4 gravity model, the Earth2014 datasets of terrain, ice thickness and inland bathymetry and the CRUST1.0 sediment and (consolidated) crustal density data. Since the global crustal density models currently available (e.g. CRUST1.0) have a limited accuracy and resolution, the comparison of individual density contributions is-for consistency-realized with a limited spectral resolution up to a spherical harmonic degree 360 (or 180). The results reveal that the topographic contribution globally varies between $$-$$ 0.33 and 0.57 m, with maxima in Himalaya and Tibet. The contribution of ice considerably modifies the geoid-to-quasigeoid separation over large parts of Antarctica and Greenland, where it reaches $$\sim $$ 0.2 m. The contributions of sediments and bedrock are less pronounced, with the values typically varying only within a few centimetres. These results, however, have still possibly large uncertainties due to the lack of information on the actual sediment and bedrock density. The contribution of lakes is mostly negligible; its maxima over the Laurentian Great Lakes and the Baikal Lake reach only several millimetres. The contribution of the sub-geoid masses is significant. It is everywhere negative and reaches extreme values of $$-$$ 4.43 m. According to our estimates, the geoid-to-quasigeoid separation globally varies within $$-$$ 4.19 and 0.26 m while the corresponding values computed according to a classical definition are only negative and reach extreme values of $$-$$ 3.5 m. A comparison of these results reveals that inaccuracies caused by disregarding the terrain geometry and mass density heterogeneities distributed within the topography and below the geoid surface can reach $$\pm $$ 2 m or more in the mountainous regions. [ABSTRACT FROM AUTHOR]

We compute globally the consolidated crust-stripped gravity disturbances/anomalies. These refined gravity field quantities are obtained from the EGM2008 gravity data after applying the topographic and crust density contrasts stripping corrections computed using the global topography/bathymetry model DTM2006.0, the global continental ice-thickness data ICE-5G, and the global crustal model CRUST2.0. All crust components density contrasts are defined relative to the reference crustal density of 2,670 kg/m. We demonstrate that the consolidated crust-stripped gravity data have the strongest correlation with the crustal thickness. Therefore, they are the most suitable gravity data type for the recovery of the Moho density interface by means of the gravimetric modelling or inversion. The consolidated crust-stripped gravity data and the CRUST2.0 crust-thickness data are used to estimate the global average value of the crust-mantle density contrast. This is done by minimising the correlation between these refined gravity and crust-thickness data by adding the crust-mantle density contrast to the original reference crustal density of 2,670 kg/m. The estimated values of 485 kg/m (for the refined gravity disturbances) and 481 kg/m (for the refined gravity anomalies) very closely agree with the value of the crust-mantle density contrast of 480 kg/m, which is adopted in the definition of the Preliminary Reference Earth Model (PREM). This agreement is more likely due to the fact that our results of the gravimetric forward modelling are significantly constrained by the CRUST2.0 model density structure and crust-thickness data derived purely based on methods of seismic refraction. [ABSTRACT FROM AUTHOR]

We analyse spatial and spectral characteristics of various refined gravity data used for modelling and gravimetric interpretation of the crust-mantle interface and the mantle-lithosphere structure. Depending on the purpose of the study, refined gravity data have either a strong or weak correlation with the Moho depths (Moho geometry). The compilation of the refined gravity data is purely based on available information on the crustal density structure obtained from seismic surveys without adopting any isostatic hypothesis. We demonstrate that the crust-stripped relative-to-mantle gravity data have a weak correlation with the CRUST2.0 Moho depths of about 0.02. Since gravitational signals due to the crustal density structure and the Moho geometry are subtracted from gravity field, these refined gravity data comprise mainly the information on the mantle lithosphere and sub-lithospheric mantle. On the other hand, the consolidated crust-stripped gravity data, obtained from the gravity field after applying the crust density contrast stripping corrections, comprise mainly the gravitational signal of the Moho geometry, although they also contain the gravitational signal due to anomalous mass density structures within the mantle. In the absence of global models of the mantle structure, the best possible option of computing refined gravity data, suitable for the recovery/refinement of the Moho interface, is to subtract the complete crust-corrected gravity data from the consolidated crust-stripped gravity data. These refined gravity data, that is, the homogenous crust gravity data, have a strong absolute correlation of about 0.99 with the CRUST2.0 Moho depths due to removing a gravitational signal of inhomogeneous density structures within the crust and mantle. Results of the spectral signal decomposition and the subsequent correlation analysis reveal that the correlation of the homogenous crust gravity data with the Moho depths is larger than 0.9 over the investigated harmonic spectrum up to harmonic degree 90. The crust-stripped relative-to-mantle gravity data correlate substantially with the Moho depths above harmonic degree 50 where the correlation exceeds 0.5. [ABSTRACT FROM AUTHOR]

Until recently, the information about the topographic density distribution has been limited to only certain regions and some countries, while missing in the global context. The UNB_TopoDens is the first model that provides the information about a lateral topographic density globally. The analysis of this model also reveals that the average topographic density for the entire continental landmass (excluding polar glaciers) is 2247 kg m−3. This density differs significantly from the value of 2670 kg m−3 that is typically adopted to represent the continental upper crustal density. In this study, we use the UNB_TopoDens density model to inspect how the topographic density variations affect interpretational properties of Bouguer gravity maps. Since this model provides also the information about density uncertainties of individual lithologies (main rock types), we estimate the corresponding errors in the Bouguer gravity data. Despite a new estimate of the average topographic density corresponds to relative changes of ∼16 per cent in values of the topographic gravity correction, these changes do not affect interpretational properties of Bouguer gravity maps. The anomalous topographic density distribution (taken with respect to the average density of 2247 kg m−3), however, modifies the Bouguer gravity pattern. We demonstrate that the gravitational contribution of anomalous topographic density is globally mostly within ±25 mGal, but much large values are detected in Himalaya, Tibet, central Andes and along the East African Rift System. Our estimates also indicate that errors in the Bouguer gravity data attributed to topographic density uncertainties are mostly less than ±15 mGal, but in mountainous regions could reach large values exceeding even ±50 mGal. Unarguably, the UNB_TopoDens model provides an improved information about the global topographic density variations and their uncertainties. Nevertheless, much more in situ measurements of rock density samples together with detailed 3-D geological models are still necessary to understand better the actual density distribution within the whole topography, particularly to mention a density change with depth. [ABSTRACT FROM AUTHOR]

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]

Recent applications in Earth sciences require geoid models to be determined with a sub-centimetre internal error. Regional models of the geoid are usually determined using discrete gravity values measured at and/or outside the Earth, and global models of the Earth gravity field and topographic surface. In this article, we review previous studies that (to some extent) discuss the estimation of the geoid internal error, and provide formulations and methodologies required for a comprehensive formal propagation of errors of gravity data and global models through a mathematical model used for regional geoid determination. The mathematical model is based on combining the inverse Poisson integral equation and the Hotine integral transform in the Helmert harmonic space; also called the one-step integration method. Calculations and tests are performed in one of the most challenging test areas ("the Colorado test area") using ground and airborne gravity observations, a global digital terrain model (DTM) for topographic effects on gravity and the geoid, and a global Earth gravitational model (EGM) for the long-wavelength components of gravity and the geoid. There are three main contributors to the total internal error of the geoid height, namely those associated with the EGM (for estimating the long-wavelength geoid height), DTM heights (for evaluation of the topographic effects on observed gravity and the geoid height), and gravity observations (for determining the short-wavelength components of the geoid height). The geoid errors stemming from the EGM formal variances of its spherical harmonic coefficients amount to 0.3 cm of the total internal error budget of the geoid. The mean value of the standard deviation of the geoid height stemming from the topographic effects (due to uncertainties of DTM heights) is 1.0 cm. The observation errors and spatial distribution (resolution) of regional ground and airborne gravity observations, i.e., the design of the project, are the dominant contributors to the internal error of the geoid height. The internal error estimate of the geoid model in the Colorado test area computed on a 1′ × 1′ grid is 2.7 cm which agrees with the differences between various geoid models that contributed to the Colorado 1-cm geoid experiment. Determining a regional geoid model with a sub-centimetre internal error in areas of rough topography, such as that of Colorado, requires high-accuracy and high-resolution gravity observations. To assess the accuracy of regional gravity measurements required for sub-centimetre geoid models, we simulate airborne gravity measurements in the Colorado test area at a practical lower flight altitude, and improve the spatial distribution of existing ground gravity observations by filling in gaps using synthetic (EGM-based) gravity disturbances. Results show that carrying out airborne gravity surveys with a gentle drape approach at an average flight altitude between 300 and 500 m above the Earth's surface provides airborne gravity measurements with a mean standard deviation of 0.75 mGal at 2.2 km spatial resolution. Thus, a sub-centimetre regional gravimetric geoid model in the Colorado test area would be achievable using the proposed configuration of the ground and airborne gravity observations and more accurate DTMs. [ABSTRACT FROM AUTHOR]

In geodesy, the geoid and the quasigeoid are used as a reference surface for heights. Despite some similarities between these two concepts, the differences between the geoid and the quasigeoid (i.e. the geoid-to-quasigeoid correction) have to be taken into consideration in some specific applications which require a high accuracy. Over the world's oceans and marginal seas, the quasigeoid and the geoid are identical. Over the continents, however, the geoid-to-quasigeoid correction could reach up to several metres especially in the mountainous, polar and geologically complex regions. Various methods have been developed and applied to compute this correction regionally in the spatial domain using detailed gravity, terrain and crustal density data. These methods utilize the gravimetric forward modelling of the topographic density structure and the direct/inverse solutions to the boundary-value problems in physical geodesy. In this article, we provide a brief summary of existing theoretical and numerical studies on the geoid-to-quasigeoid correction. We then compare these methods with the newly developed procedure and discuss some numerical and practical aspects of computing this correction. In global applications, the geoid-to-quasigeoid correction can conveniently be computed in the spectral domain. For this purpose, we derive and present also the spectral expressions for computing this correction based on applying methods for a spherical harmonic analysis and synthesis of global gravity, terrain and crustal structure models. We argue that the newly developed procedure for the regional gravity-to-potential conversion, applied for computing the geoid-to-quasigeoid correction in the spatial domain, is numerically more stable than the existing inverse models which utilize the gravity downward continuation. Moreover, compared to existing spectral expressions, our definition in the spectral domain takes not only the terrain geometry but also the mass density heterogeneities within the whole Earth into consideration. In this way, the geoid-to-quasigeoid correction and the respective geoid model could be determined more accurately. [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]

New integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients are derived in this article. They provide more options for continuation of gravitational gradient combinations and extend available mathematical apparatus formulated for this purpose up to now. The starting point represents the analytical solution of the spherical gradiometric boundary value problem in the spatial domain. Applying corresponding differential operators on the analytical solution of the spherical gradiometric boundary value problem, a total of 18 integral formulas are provided. Spatial and spectral forms of isotropic kernels are given and their behaviour for parameters of a GOCE-like satellite is investigated. Correctness of the new integral formulas and the isotropic kernels is tested in a closed-loop simulation. The derived integral formulas and the isotropic kernels form a theoretical basis for validation purposes and geophysical applications of satellite gradiometric data as provided currently by the GOCE mission. They also extend the well-known Meissl scheme. [ABSTRACT FROM AUTHOR]

Abstract: In classical isostatic models for a gravimetric recovery of the Moho parameters (i.e., Moho depths and density contrast) the isostatic gravity anomalies are usually defined based on the assumption that the topographic mass surplus and the ocean mass deficiency are compensated within the Earth''s crust. As acquired in this study, this assumption yields large disagreements between isostatic and seismic Moho models. To assess the effects not accounted for in classical isostatic models, we conduct a number of numerical experiments using available global gravity and crustal structure models. First, we compute the gravitational contributions of mass density contrasts due to ice and sediments, and subsequently evaluate respective changes in the Moho geometry. Residual differences between the gravimetric and seismic Moho models are then used to predict a remaining non-isostatic gravity signal, which is mainly attributed to unmodeled density structures and other geophysical phenomena. We utilize three recently developed computational schemes in our numerical studies. The apparatus of spherical harmonic analysis and synthesis is applied in forward modeling of the isostatic gravity disturbances. The Moho depths are estimated globally on a 1 arc-deg equiangular grid by solving the Vening-Meinesz Moritz inverse problem of isostasy. The same estimation model is applied to evaluate the differences between the isostatic and seismic models. We demonstrate that the application of the ice and sediment density contrasts stripping gravity corrections is essential for a more accurate determination of the Moho geometry. We also show that the application of the additional non-isostatic correction further improves the agreement between the Moho models derived based on gravity and seismic data. Our conclusions are based on comparing the gravimetric results with the CRUST2.0 global crustal model compiled using results of seismic surveys. [Copyright &y& Elsevier]

Shipborne gravity can be used to refine altimeter-derived gravity whose accuracy is low in shallow waters and areas with complex submarine topography. As altimeter-derived gravity only within a small radius around the shipborne data can be corrected by traditional methods, a new method based on multi-layer perceptron (MLP) neural network is proposed to refine the altimeter-derived gravity. Input variables of MLP include the positional information at observation points and geophysical information (from our own South China Sea gravity anomaly model (SCSGA) V1.0 and bathymetry model ETOPO1) at grid points around observation points. Output variables of MLP are the refined residual gravity anomalies at observation points. Training shipborne data are classified into four cases to train four MLP models, which are used to predict the refined gravity anomaly model SCSGA V1.1. Then all of the training shipborne data are used for training an MLP model to predict the refined gravity anomaly model SCSGA V1.2. Assessed by testing shipborne data, the accuracy of SCSGA V1.2 is 0.14 mGal higher than that of SCSGA V1.0, and similar to that of SCSGA V1.1. Compared with the original gravity anomaly model (SCSGA V1.0), the accuracy of the refined gravity anomaly model (SCSGA V1.2) by MLP is improved by 4.4% in areas where the training data are concentrated, and also improved by 2.2% in other areas. Therefore, the method of MLP can be used to refine the altimeter-derived gravity model by shipborne gravity, overcoming the problem of limited correction radius for traditional methods. [ABSTRACT FROM AUTHOR]

Local gravity field modelling demands high-quality gravity data as well as an appropriate mathematical model. Particularly in coastal areas, there may be different types of gravity observations available, for instance, terrestrial, aerial, marine gravity, and satellite altimetry data. Thus, it is important to develop a proper tool to merge the different data types for local gravity field modelling and determination of the geoid. In this study, radial basis functions, as a commonly useful tool for gravity data integration, are employed to model the gravity potential field of the southern part of Iran using terrestrial gravity anomalies, gravity anomalies derived from re-tracked satellite altimetry, marine gravity anomalies, and gravity anomalies synthesized from an Earth gravity model. Reference GNSS/levelling (geometric) geoidal heights are used to evaluate the accuracy of the estimated local gravity field model. The gravimetric geoidal heights are in acceptable agreement with the geometric ones in terms of the standard deviation and the mean value which are 4.1 and 12 cm, respectively. Besides, the reference benchmark of the national first-order levelling network of Iran is located in the study area. The derived gravity model was used to compute the gravity potential difference at this point and then transformed into a height difference which results in the value of the shift of this benchmark with respect to the geoid. The estimated shift shows a good agreement with previously published studies. [ABSTRACT FROM AUTHOR]