Ellipsoidal vertical deflections and ellipsoidal gravity disturbance: Case studies

First, we present three different definitions of the vertical which relate to (i) astronomical longitude and astronomical latitude as spherical coordinates in gravity space, (ii) Gauss surface normal coordinates (also called geodetic coordinates) of type ellipsoidal longitude and ellipsoidal latitude and (iii) Jacobi ellipsoidal coordinates of type spheroidal longitude and spheroidal latitude in geometry space. Up to terms of second order those vertical deflections agree to each other. Vertical deflections and gravity disturbances relate to a reference gravity potential. In order to refer the horizontal and vertical components of the disturbing gravity field to a reference gravity field, which is physically meaningful, we have chosen the Somigliana-Pizzetti gravity potential as well as its gradient. Second, we give a new closed-form representation of Somigliana-Pizzetti gravity, accurate to the sub Nano Gal level. Third, we represent the gravitational disturbing potential in terms of Jacobi ellipsoidal harmonics. As soon as we take reference to a normal potential of Somigliana-Pizzetti type, the ellipsoidal harmonics of degree/order (0,0), (1,0), (1, − 1), (1,1) and (2,0) are eliminated from the gravitational disturbing potential. Fourth, we compute in all detail the gradient of the gravitational disturbing potential, in particular in orthonormal ellipsoidal vector harmonics. Proper weighting functions for orthonormality on the International Reference Ellipsoid are constructed and tabulated. In this way, we finally arrive at an ellipsoidal harmonic representation of vertical deflections and gravity disturbances. Fifth, for an ellipsoidal harmonic Gravity Earth Model (SEGEN: http://www.uni-stuttgart.de/gi/research/paper/coefficients/coefficients.zip) up to degree/order 360/360 we compute the global maps of ellipsoidal vertical deflections and ellipsoidal gravity disturbances which transfer a great amount of geophysical information in a properly chosen equiareal ellipsoidal map projection.     

用GRACE卫星跟踪数据反演地球重力场

利用141天GRACE卫星观测资料,包括K波段、星载加速度和卫星轨道数据,反演了80阶地球重力场模型IGGGRACE01S,该模型在半波长为500km的空间分辨率上,确定大地水准面的精度约为0012m,中长波(<80阶)精度优于重力卫星发射以前研制的重力场模型. 与EIGEN_GRACE02S、EIGEN_CHAMP03S和EGM96模型的位系数相比,该模型系数最接近于EIGEN_GRACE02S,与另两个模型差异较大. 比较几种模型确定的全球重力异常和大地水准面起伏,结果发现IGGGRACE01S与EIGEN_GRACE02S模型的计算结果比较接近,与EGM96模型结果差异较大,差别较大地区主要在南极等地区. 对于中国大陆,比较IGGGRACE01S模型（前72阶）计算的重力异常和NIMA重力异常数据（25°×25°网格）,两者之间的标准偏差为48mGal.     

Structure and State of Stress of the Chilean Subduction Zone from Terrestrial and Satellite-Derived Gravity and Gravity Gradient Data

It is well known that the quality of gravity modelling of the Earth’s lithosphere is heavily dependent on the limited number of available terrestrial gravity data. More recently, however, interest has grown within the geoscientific community to utilise the homogeneously measured satellite gravity and gravity gradient data for lithospheric scale modelling. Here, we present an interdisciplinary approach to determine the state of stress and rate of deformation in the Central Andean subduction system. We employed gravity data from terrestrial, satellite-based and combined sources using multiple methods to constrain stress, strain and gravitational potential energy (GPE). Well-constrained 3D density models, which were partly optimised using the combined regional gravity model IMOSAGA01C (Hosse et al. in Surv Geophys, 2014, this issue), were used as bases for the computation of stress anomalies on the top of the subducting oceanic Nazca plate and GPE relative to the base of the lithosphere. The geometries and physical parameters of the 3D density models were used for the computation of stresses and uplift rates in the dynamic modelling. The stress distributions, as derived from the static and dynamic modelling, reveal distinct positive anomalies of up to 80 MPa along the coastal Jurassic batholith belt. The anomalies correlate well with major seismicity in the shallow parts of the subduction system. Moreover, the pattern of stress distributions in the Andean convergent zone varies both along the north–south and west–east directions, suggesting that the continental fore-arc is highly segmented. Estimates of GPE show that the high Central Andes might be in a state of horizontal deviatoric tension. Models of gravity gradients from the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission were used to compute Bouguer-like gradient anomalies at 8 km above sea level. The analysis suggests that data from GOCE add significant value to the interpretation of lithospheric structures, given that the appropriate topographic correction is applied.     

Combined Regional Gravity Model of the Andean Convergent Subduction Zone and Its Application to Crustal Density Modelling in Active Plate Margins

The Central Andean subduction system is one of the most active geological structures on Earth. Although there have been a few previous studies, the structure and dynamics of the system are still not well understood. In the present study, we determine a combined regional gravity model of the Andean convergent subduction region for constraining lithospheric models. After a thorough validation and cleaning of the terrestrial gravity and height databases, the method of Least Squares Collocation was applied to consistently combine terrestrial and satellite gravity data, putting much emphasis on the stochastic modelling of the individual data components. As a result, we computed the first high-resolution regional gravity model of the study region that includes GOCE satellite gravity information. The inclusion of GOCE is an essential distinction from the independent global gravity model EGM2008. Validation against EGM2008 reveals that our regional solution is very consistent in regions where terrestrial gravity data are available, but shows systematic differences in areas with terrestrial data gaps. Artefacts in the EGM2008 of up to 150 mGal could be identified. The new combined regional model benefits from the very homogeneous error characteristics and accuracy of GOCE gravity data in the long-to-medium wavelengths down to 80–100 km. Reliable density modelling became possible also in the region of Central Andes, which lacks terrestrial gravity data. Finally, density models were adapted to fit the new regional gravity field solution. The results clearly demonstrate the capabilities of GOCE to better constrain lithospheric models.     

Layer-Based Modelling of the Earth’s Gravitational Potential up to 10-km Scale in Spherical Harmonics in Spherical and Ellipsoidal Approximation

Global forward modelling of the Earth’s gravitational potential, a classical problem in geophysics and geodesy, is relevant for a range of applications such as gravity interpretation, isostatic hypothesis testing or combined gravity field modelling with high and ultra-high resolution. This study presents spectral forward modelling with volumetric mass layers to degree 2190 for the first time based on two different levels of approximation. In spherical approximation, the mass layers are referred to a sphere, yielding the spherical topographic potential. In ellipsoidal approximation where an ellipsoid of revolution provides the reference, the ellipsoidal topographic potential (ETP) is obtained. For both types of approximation, we derive a mass layer concept and study it with layered data from the Earth2014 topography model at 5-arc-min resolution. We show that the layer concept can be applied with either actual layer density or density contrasts w.r.t. a reference density, without discernible differences in the computed gravity functionals. To avoid aliasing and truncation errors, we carefully account for increased sampling requirements due to the exponentiation of the boundary functions and consider all numerically relevant terms of the involved binominal series expansions. The main outcome of our work is a set of new spectral models of the Earth’s topographic potential relying on mass layer modelling in spherical and in ellipsoidal approximation. We compare both levels of approximations geometrically, spectrally and numerically and quantify the benefits over the frequently used rock-equivalent topography (RET) method. We show that by using the ETP it is possible to avoid any displacement of masses and quantify also the benefit of mapping-free modelling. The layer-based forward modelling is corroborated by GOCE satellite gradiometry, by in-situ gravity observations from recently released Antarctic gravity anomaly grids and degree correlations with spectral models of the Earth’s observed geopotential. As the main conclusion of this work, the mass layer approach allows more accurate modelling of the topographic potential because it avoids 10–20-mGal approximation errors associated with RET techniques. The spherical approximation is suited for a range of geophysical applications, while the ellipsoidal approximation is preferable for applications requiring high accuracy or high resolution.     

On High-Frequency Topography-Implied Gravity Signals for a Height System Unification Using GOCE-Based Global Geopotential Models

National height reference systems have conventionally been linked to the local mean sea level, observed at individual tide gauges. Due to variations in the sea surface topography, the reference levels of these systems are inconsistent, causing height datum offsets of up to ±1–2 m. For the unification of height systems, a satellite-based method is presented that utilizes global geopotential models (GGMs) derived from ESA’s satellite mission Gravity field and steady-state Ocean Circulation Explorer (GOCE). In this context, height datum offsets are estimated within a least squares adjustment by comparing the GGM information with measured GNSS/leveling data. While the GNSS/leveling data comprises the full spectral information, GOCE GGMs are restricted to long wavelengths according to the maximum degree of their spherical harmonic representation. To provide accurate height datum offsets, it is indispensable to account for the remaining signal above this maximum degree, known as the omission error of the GGM. Therefore, a combination of the GOCE information with the high-resolution Earth Gravitational Model 2008 (EGM2008) is performed. The main contribution of this paper is to analyze the benefit, when high-frequency topography-implied gravity signals are additionally used to reduce the remaining omission error of EGM2008. In terms of a spectral extension, a new method is proposed that does not rely on an assumed spectral consistency of topographic heights and implied gravity as is the case for the residual terrain modeling (RTM) technique. In the first step of this new approach, gravity forward modeling based on tesseroid mass bodies is performed according to the Rock–Water–Ice (RWI) approach. In a second step, the resulting full spectral RWI-based topographic potential values are reduced by the effect of the topographic gravity field model RWI_TOPO_2015, thus, removing the long to medium wavelengths. By using the latest GOCE GGMs, the impact of topography-implied gravity signals on the estimation of height datum offsets is analyzed in detail for representative GNSS/leveling data sets in Germany, Austria, and Brazil. Besides considerable changes in the estimated offset of up to 3 cm, the conducted analyses show that significant improvements of 30–40% can be achieved in terms of a reduced standard deviation and range of the least squares adjusted residuals.     

Optimal Model for Geoid Determination from Airborne Gravity

Two different approaches for transformation of airborne gravity disturbances, derived from gravity observations at low-elevation flying platforms, into geoidal undulations are formulated, tested and discussed in this contribution. Their mathematical models are based on Green's integral equations. They are in these two approaches defined at two different levels and also applied in a mutually reversed order. While one of these approaches corresponds to the classical method commonly applied in processing of ground gravity data, the other approach represents a new method for processing of gravity data in geoid determination that is unique to airborne gravimetry. Although theoretically equivalent in the continuous sense, both approaches are tested numerically for possible numerical advantages, especially due to the inverse of discretized Fredholm's integral equation of the first kind applied on different data. High-frequency synthetic gravity data burdened by the 2-mGal random noise, that are expected from current airborne gravity systems, are used for numerical testing. The results show that both approaches can deliver for the given data a comparable cm-level accuracy of the geoidal undulations. The new approach has, however, significantly higher computational efficiency. It would be thus recommended for real life geoid computations. Additional errors related to regularization of gravity data and the geoid, and to accuracy of the reference field, that would further deteriorate the quality of estimated geoidal undulations, are not considered in this study.     

Regional gravity field modeling based on rectangular harmonic analysis

Regional gravity field modeling with high-precision and high-resolution is one of the most important scientific objectives in geodesy,and can provide fundamental information for geophysics,geodynamics,seismology,and mineral exploration.Rectangular harmonic analysis(RHA)is proposed for regional gravity field modeling in this paper.By solving the Laplace’s equation of gravitational potential in local Cartesian coordinate system,the rectangular harmonic expansions of disturbing potential,gravity anomaly,gravity disturbance,geoid undulation and deflection of the vertical are derived,and so are the formula for signal degree variance and error degree variance of the rectangular harmonic coefficients(RHC).We also present the mathematical model and detailed algorithm for the solution of RHC using RHA from gravity observations.In order to reduce the edge effects caused by periodic continuation in RHA,we propose the strategy of extending the size of computation domain.The RHA-based modeling method is validated by conducting numerical experiments based on simulated ground and airborne gravity data that are generated from geopotential model EGM2008 and contaminated by Gauss white noise with standard deviation of 2 mGal.The accuracy of the 2.5′×2.5′geoid undulations computed from ground and airborne gravity data is 1 and 1.4cm,respectively.The standard error of the gravity disturbances that downward continued from the flight height of 4 km to the geoid is only 3.1 mGal.Numerical results confirm that RHA is able to provide a reliable and accurate regional gravity field model,which may be a new option for the representation of the fine structure of regional gravity field.     

On Solving the Forward Problem of Gravimetry in Curvilinear and Cartesian Coordinates: Krasovskii’s Ellipsoid and Plane Modeling

Correcting the effects of the sphericity of the Earth in the results of the interpretation of gravimetric data is a topical issue in modern gravimetry. Estimating the error of the gravity field calculations due to the replacement of the spherical Earth model by the plane model is an important part of this problem. In this paper, a method is proposed for transforming the plane density models into spherical ones and vice versa. Algorithms for calculating the vertical component of gravity field for both model types are presented. For two extensive plane models of the Earth’s density, their transformation into spherical models is carried out and the resulting gravity fields are compared. The relative root mean square residuals between the fields calculated with this replacement are at most 5%.     

新疆库尔勒地区秋季下平流层重力波特性

利用2011年秋季无线电探空数据,采用矢端曲线法首次分析了新疆库尔勒地区下平流层重力波特征参量,得到36组准单色重力波的结果.结果统计显示：库尔勒秋季下平流层重力波垂直波长、水平波长平均值分别为2.8 km和580 km,固有频率平均值为1.74f.垂直传播方向以上传为主,约占78%,其中下传重力波水平波长较短,固有频率较高.水平传播方向以西北和东南为主,各占1/3,其中上传（下传）重力波水平传播方向以西北（东南）居多,这与热带低纬站点和其他中纬站点观测结果不同.与其他站点比较,库尔勒地区ŵ/f最小,中高纬地区水平波长、垂直波长随纬度增加大致有减小的趋势,库尔勒地区偏离这一趋势,波长偏大.     

Pacific Plate motion and undulations in geoid and bathymetry

Previous studies have shown that the Pacific geoid and gravity fields exhibit lineated anomalies, trending approximately in the direction of absolute plate motion over the underlying mantle. Because the undulations obliquely cross fracture zones they have often been attributed a convective origin. Recently, lithospheric boudinage caused by diffuse extension has been proposed as a possible mechanism. We have examined the undulations in the free-air anomalies, geoid and bathymetry over a portion of the Pacific Plate to determine quantitatively how the undulations are related to plate motion. We compare the observed data to an axisymmetric, sinusoidal undulation defined in an arbitrary frame of reference; in particular, we seek the north pole of this reference frame that maximizes the correlation between data and model. Poles that are close to the Pacific hotspot pole represent copolar undulations possibly related to plate motion. The distance between the best-fitting poles and the hotspot pole is determined as a function of undulation wavelength and reveals several minima (with distance < 10°) for discrete geoid wavebands centered on wavelengths of 160 km, 225 km, 287 km, 400 km, 660 km, 850 km, 1000 km and 1400 km. Bathymetry data have copolar bathymetric expressions as well, giving an implied admittance of 2–3 m/km. The most co-polar geoid/bathymetry undulations (with poles within 2–3° of the average Pacific Euler pole) have wavelengths of 280 km and 1050 km, respectively. The latter could have a convective origin or be related to the spacing of hotspot swells. The former may reflect lithospheric boudinage formed in response to diffuse extension, but could also have a dynamic origin since flexural dampening may only have attenuated the bathymetric amplitude by 50% or less. Radiometric dating of volcanic ridges found in the troughs of prominent gravity lineations gives ages that correlate well with documented changes in Pacific and Indo/Australian Plate motion, suggesting the ridges formed in response to intermittent plate boundary stresses and not as a direct consequence of small-scale convection or diffuse extension.     

Calculation of atmospheric electric fields penetrating from the ionosphere

The spatial distributions of electric fields and currents in the Earth’s atmosphere are calculated. Electric potential distributions typical of substorms and quiet geomagnetic conditions are specified in the ionosphere. The Earth is treated as a perfect conductor. The atmosphere is considered as a spherical layer with a given height dependence of electrical conductivity. With the chosen conductivity model and an ionospheric potential of 300 kV with respect to the Earth, the electric field near the ground is vertical and reaches 110 Vm⁻¹. With the 60-kV potential difference in the polar cap of the ionosphere, the electric field disturbances with a vertical component of up to 13 V m⁻¹ can occur in the atmosphere. These disturbances are maximal near the ground. If the horizontal scales of field nonuniformity are over 100 km, the vertical component of the electric field near the ground can be calculated with the one-dimensional model. The field and current distributions in the upper atmosphere can be obtained only from the three-dimensional model. The numerical method for solving electrical conductivity problems makes it possible to take into account conductivity inhomogeneities and the ground relief.     

What Can be Expected from the GRACE-FO Laser Ranging Interferometer for Earth Science Applications?

The primary objective of the gravity recovery and climate experiment follow-on (GRACE-FO) satellite mission, due for launch in August 2017, is to continue the GRACE time series of global monthly gravity field models. For this, evolved versions of the GRACE microwave instrument, GPS receiver, and accelerometer will be used. A secondary objective is to demonstrate the effectiveness of a laser ranging interferometer (LRI) in improving the satellite-to-satellite tracking measurement performance. In order to investigate the expected enhancement for Earth science applications, we have performed a full-scale simulation over the nominal mission lifetime of 5 years using a realistic orbit scenario and error assumptions both for instrument and background model errors. Unfiltered differences between the synthetic input and the finally recovered time-variable monthly gravity models show notable improvements with the LRI, on a global scale, of the order of 23 %. The gain is realized for wavelengths smaller than 240 km in case of Gaussian filtering but decreases to just a few percent when anisotropic filtering is applied. This is also confirmed for some typical regional Earth science applications which show randomly distributed patterns of small improvements but also degradations when using DDK4-filtered LRI-based models. Analysis of applied error models indicates that accelerometer noise followed by ocean tide and non-tidal mass variation errors are the main contributors to the overall GRACE-FO gravity model error. Improvements in these fields are therefore necessary, besides optimized constellations, to make use of the increased LRI accuracy and to significantly improve gravity field models from next-generation gravity missions.     

A precise gravimetric geoid model in a mountainous area with scarce gravity data: a case study in central Turkey

In mountainous regions with scarce gravity data, gravimetric geoid determination is a difficult task that needs special attention to obtain reliable results satisfying the demands, e.g., of engineering applications. The present study investigates a procedure for combining a suitable global geopotential model and available terrestrial data in order to obtain a precise regional geoid model for Konya Closed Basin (KCB). The KCB is located in the central part of Turkey, where a very limited amount of terrestrial gravity data is available. Various data sources, such as the Turkish digital elevation model with 3 ?? × 3?? resolution, a recently published satellite-only global geopotential model from the Gravity Recovery and Climate Experiment satellite (GRACE) and the ground gravity observations, are combined in the least-squares sense by the modified Stokes?? formula. The new gravimetric geoid model is compared with Global Positioning System (GPS)/levelling at the control points, resulting in the Root Mean Square Error (RMS) differences of ±6.4 cm and 1.7 ppm in the absolute and relative senses, respectively. This regional geoid model appears to be more accurate than the Earth Gravitational Model 2008, which is the best global model over the target area, with the RMS differences of ±8.6 cm and 1.8 ppm in the absolute and relative senses, respectively. These results show that the accuracy of a regional gravimetric model can be augmented by the combination of a global geopotential model and local terrestrial data in mountainous areas even though the quality and resolution of the primary terrestrial data are not satisfactory to the geoid modelling procedure.     

Experiences with the use of mass-density maps in residual gravity forward modelling

In many modern local and regional gravity field modelling concepts, the short-wavelength gravitational signal modeled by the residual terrain modelling (RTM) technique is used to augment global geopotential models, or to smooth observed gravity prior to data gridding. In practice, the evaluation of RTM effects mostly relies on a constant density assumption, because of the difficulty and complexity of obtaining information on the actual distribution of density of topographic masses. Where the actual density of topographic masses deviates from the adopted value, errors are present in the RTM mass-model, and hence, in the forward-modelled residual gravity field. In this paper we attempt to overcome this problem by combining the RTM technique with a high-resolution mass-density model. We compute RTM gravity quantities over New Zealand, with different combinations of elevation models and mass-density assumptions using gravity and GPS/levelling measurements, precise terrain and bathymetry models, a high-resolution mass-density model and constant density assumptions as main input databases. Based on gravity observations and the RTM technique, optimum densities are detected for North Island of ~2500 kg m^?3, South Island of ~2600 kg m^?3, and the whole New Zealand of ~2590 kg m^?3. Comparison among the three sets of residual gravity disturbances computed from different mass-density assumptions show that, together with a global potential model, the high-resolution New Zealand density model explains ~89.5% of gravitational signals, a constant density assumption of 2670 kg m^?3 explains ~90.2%, while a regionally optimum mass-density explains ~90.3%. Detailed comparison shows that the New Zealand density model works best over areas with small residual heights. Over areas with larger residual heights, subsurface density variations appear to affect the residual gravity disturbance. This effect is found to reach about 30 mGal over Southern Alpine Fault. In order to improve the RTM modelling with mass-density maps, a higher-quality mass-density model that provides radially varying mass-density data would be desirable.     

A new ellipsoidal gravimetric, satellite altimetry and astronomic boundary value problem, a case study: The geoid of Iran

A new gravimetric, satellite altimetry, astronomical ellipsoidal boundary value problem for geoid computations has been developed and successfully tested. This boundary value problem has been constructed for gravity observables of the type (i) gravity potential, (ii) gravity intensity (i.e. modulus of gravity acceleration), (iii) astronomical longitude, (iv) astronomical latitude and (v) satellite altimetry observations. The ellipsoidal coordinates of the observation points have been considered as known quantities in the set-up of the problem in the light of availability of GPS coordinates. The developed boundary value problem is ellipsoidal by nature and as such takes advantage of high precision GPS observations in the set-up. The algorithmic steps of the solution of the boundary value problem are as follows:

- Application of the ellipsoidal harmonic expansion complete up to degree and order 360 and of the ellipsoidal centrifugal field for the removal of the effect of global gravity and the isostasy field from the gravity intensity and the astronomical observations at the surface of the Earth.

- Application of the ellipsoidal Newton integral on the multi-cylindrical equal-area map projection surface for the removal from the gravity intensity and the astronomical observations at the surface of the Earth the effect of the residual masses at the radius of up to 55 km from the computational point.

- Application of the ellipsoidal harmonic expansion complete up to degree and order 360 and ellipsoidal centrifugal field for the removal from the geoidal undulations derived from satellite altimetry the effect of the global gravity and isostasy on the geoidal undulations.

- Application of the ellipsoidal Newton integral on the multi-cylindrical equal-area map projection surface for the removal from the geoidal undulations derived from satellite altimetry the effect of the water masses outside the reference ellipsoid within a radius of 55 km around the computational point.

- Least squares solution of the observation equations of the incremental quantities derived from aforementioned steps in order to obtain the incremental gravity potential at the surface of the reference ellipsoid.

- The removed effects at the application points are restored on the surface of reference ellipsoid.

- Application of the ellipsoidal Bruns’ formula for converting the potential values on the surface of the reference ellipsoid into the geoidal heights with respect to the reference ellipsoid.

- Computation of the geoid of Iran has successfully tested this new methodology.

A study of regional gravity field recovery from GOCE vertical gravity gradient data in the Auvergne test area using collocation

Gravity field and steady-state Ocean Circulation Explorer (GOCE) is the first satellite mission that observes gravity gradients from the space, to be primarily used for the determination of high precision global gravity field models. However, the GOCE gradients, having a dense data distribution, may potentially provide better predictions of the regional gravity field than those obtained using a spherical harmonic Earth Geopotential Model (EGM). This is investigated in Auvergne test area using Least Squares Collocation (LSC) with GOCE vertical gravity gradient anomalies (T_zz), removing the long wavelength part from EGM2008 and the short wavelength part by residual terrain modelling (RTM). The results show that terrain effects on the vertical gravity gradient are significant at satellite altitude, reaching a level of 0.11 E?tv?s unit (E.U.) in the mountainous areas. Removing the RTM effects from GOCE T_zz leads to significant improvements on the LSC predictions of surface gravity anomalies and quasigeoid heights. Comparison with ground truth data shows that using LSC surface free air gravity anomalies and quasi-geoid heights are recovered from GOCE T_zz with standard deviations of 11 mGal and 18 cm, which is better than those obtained by using GOCE EGMs, demonstrating that information beyond the maximal degree of the GOCE EGMs is present. Investigation of using covariance functions created separately from GOCE T_zz and terrestrial free air gravity anomalies, suggests that both covariance functions give almost identical predictions. However, using covariance function obtained from GOCE T_zz has the effect that the predicted formal average error estimates are considerably larger than the standard deviations of predicted minus observed gravity anomalies. Therefore, GOCE T_zz should be used with caution to determine the covariance functions in areas where surface gravity anomalies are not available, if error estimates are needed.     

Gravity Spectra from the Density Distribution of Earth’s Uppermost 435 km

The Earth masses reside in a near-hydrostatic equilibrium, while the deviations are, for example, manifested in the geoid, which is nowadays well determined by satellite gravimetry. Recent progress in estimating the density distribution of the Earth allows us to examine individual Earth layers and to directly see how the sum approaches the observed anomalous gravitational field. This study evaluates contributions from the crust and the upper mantle taken from the LITHO1.0 model and quantifies the gravitational spectra of the density structure to the depth of 435 km. This is done without isostatic adjustments to see what can be revealed with models like LITHO1.0 alone. At the resolution of 290 km (spherical harmonic degree 70), the crustal contribution starts to dominate over the upper mantle and at about 150 km (degree 130) the upper mantle contribution is nearly negligible. At the spatial resolution \(<150\,\hbox {km},\) the spectra behavior is driven by the crust, the mantle lid and the asthenosphere. The LITHO1.0 model was furthermore referenced by adding deeper Earth layers from ak135, and the gravity signal of the merged model was then compared with the observed satellite-only model GOCO05s. The largest differences are found over the tectonothermal cold and old (such as cratonic), and over warm and young areas (such as oceanic ridges). The misfit encountered comes from the mantle lid where a velocity–density relation helped to reduce the RMS error by 40%. Global residuals are also provided in terms of the gravitational gradients as they provide better spatial localization than gravity, and there is strong observational support from ESA’s satellite gradiometry mission GOCE down to the spatial resolution of 80–90 km.     

Grim4-C1, -C2p: combination solutions of the global earth gravity field

On the basis of the GRIM4-S1 satellite-only Earth gravity model, being accomplished in a common effort by DGFI and GRGS, a combination solution, called GRIM4-C1, has been derivcd using 1° × 1° mean gravity anomalies and 1° × 1° Seasat altimeter derived mean geoid undulations. In the meantime improvements could be achieved by incorporating more tracking data (GEOSAT, SPOT2-DORIS) into the solution, resulting in the two new parallel versions, the satellite-only gravity model GRIM4-S2 and the combined solution GRIM4-C2p (preliminary). All GRIM4 Earth gravity models cover the spectral gravitational constituents complete up to degree and order 50.In this report the emphasis is on the discussion of the combined gravity models: combination and estimation techniques, capabilities for application in precise satellite orbit computation and accuracies in long wavelength geoid representation. It is shown that with the new generation of global gravity models general purpose satellite-only models are no longer inferior to combination solutions if applied to satellite orbit restitution.     

Ellipsoidal Neumann geodetic boundary-value problem based on surface gravity disturbances: case study of Iran

An ellipsoidal Neumann type geodetic boundary-value problem (GBVP) for the computation of disturbing potential on the surface of the Earth based on the surface gravity disturbance as the boundary data is formulated. The solution methodology of the GBVP can be algorithmically summarized as follows: (i) using global navigation satellite systems (GNSS) coordinates of the gravity stations, the surface gravity disturbances are generated as the boundary data. (ii) Applying the deflection correction to the gravity disturbances to arrive at the derivative of the surface disturbing potential along the ellipsoidal normal. (iii) Removing the low frequencies part of the gravity field using harmonic expansion to degree and order 110. (iv) Using the short wavelength part of the corrected gravity disturbances derived in the previous section as the boundary data within the constructed GBVP to derive the short wavelength disturbing potential over the Earth surface. (v) The computed shortwave length signals of disturbing potentials are converted to disturbing potential values by restoring the removed effects.