The Earth's gravity field

The Earth's gravity field can be determined from gravity measurements made on the surface of the Earth, and through the analysis of the motion of Earth satellites. Gravity data can be used to solve the boundary value problem of gravimetric geodesy in various ways, from the classical formulation using a geoid to the concept of a reference surface interior to the masses of the Earth to a statistical method. We now have gravity information for 1⁰ data blocks over 46% of the Earth's surface and more than several million point measurements available.Satellite observations such as range, range-rate, and optical data have been analyzed to determine potential coefficients used to describe the Earth's gravitational potential field. Coefficients, in a spherical harmonic expansion to degree 12, can be determined from satellite data alone, and to at least degree 20 when the satellite data is combined with surface gravity material. Recent solutions for potential coefficients agree well to degree 4, but with increasing disagreement at higher degrees.     

A comparison of satellite systems for gravity field measurements

A detailed and accurate Earth gravity field model is important both to geophysical progress and to the precise tracking necessary for interpretation of geophysical experiments. Various satellite techniques which may be used to determine the Earth's gravity field are compared and their ability to recover the long wavelength and short wavelength features of the field are described.A high-low configuration satellite-to-satellite tracking mission is recommended for the determination of the long wavelength portion of the gravity field. Satellite altimetry and satellite gradiometry experiments are recommended for determination of the short wavelength portion of the gravity field.     

Interpretation of general geophysical patterns in Iran based on GRACE gradient component analysis

Only with satellites it is possible to cover the entire Earth densely with gravity field related measurements of uniform quality within a short period of time. However, due to the altitude of the satellite orbits, the signals of individual local masses are strongly damped. Based on the approach of Petrovskaya and Vershkov we determine the gravity gradient tensor directly from the spherical harmonic coefficients of the recent EIGEN-GL04C combined model of the GRACE satellite mission. Satellite gradiometry can be used as a complementary tool to gravity and geoid information in interpreting the general geophysical and geodynamical features of the Earth. Due to the high altitude of the satellite, the effects of the topography and the internal masses of the Earth are strongly damped. However, the gradiometer data, which are nothing else than the second order spatial derivatives of the gravity potential, efficiently counteract signal attenuation at the low and medium frequencies. In this article we review the procedure for estimating the gravity gradient components directly from spherical harmonics coefficients. Then we apply this method as a case study for the interpretation of possible geophysical or geodynamical patterns in Iran. We found strong correlations between the cross-components of the gravity gradient tensor and the components of the deflection of vertical, and we show that this result agrees with theory. Also, strong correlations of the gravity anomaly, geoid model and a digital elevation model were found with the diagonal elements of the gradient tensor.     

A new regard on the tectonic map of the Arabian–African region inferred from the satellite gravity analysis

Satellite gravimetry is a powerful and reliable tool for regional tectono-geodynamic zonation. The studied region contains intricate geodynamical features (high seismological indicators, active rift systems and collision processes), richest structural arrangement (existence of mosaic blocks of oceanic and continental Earth’s crust of various age), and a number of high-amplitude gravity anomalies and complex magnetic pattern. The most hydrocarbon reserves of the world and other important economic deposits occur in this region. Comprehensive analysis of satellite gravity data with application of different approaches was used to develop a sequence of maps specifying crucial properties of the region deep structure. Careful examination of numerous geological sources and their combined examination with satellite gravity (main), magnetic, GPS, seismic, seismological and some other geophysical data enabled to develop a new tectonic map of the Arabian–African region. Integrated analysis of series of gravity map transformations and certain geological indicators allowed to reveal significant geodynamic features of the region.     

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.     

The Integrated Global Geodetic Observing System (IGGOS) viewed from the perspective of history

Towards the end of the 19th century, geodetic observation techniques allowed it to create geodetic networks of continental size. The insight that big networks can only be set up through international collaboration led to the establishment of an international collaboration called “Central European Arc Measurement”, the predecessor of the International Association of Geodesy (IAG), in 1864. The scope of IAG activities was extended already in the 19th century to include gravity.At the same time, astrometric observations could be made with an accuracy of a few tenths of an arcsecond. The accuracy stayed roughly on this level, till the space age opened the door for milliarcsecond (mas) astrometry. Astrometric observations allowed it at the end of the 19th century to prove the existence of polar motion. The insight that polar motion is almost unpredictable led to the establishment of the International Latitude Service (ILS) in 1899.The IAG and the ILS were the tools (a) to establish and maintain the terrestrial and the celestial reference systems, including the transformation parameters between the two systems, and (b) to determine the Earth's gravity field.Satellite-geodetic techniques and astrometric radio-interferometric techniques revolutionized geodesy in the second half of the 20th century. Satellite Laser Ranging (SLR) and methods based on the interferometric exploitation of microwave signals (stemming from Quasars and/or from satellites) allow it to realize the celestial reference frame with (sub-)mas accuracy, the global terrestrial reference frame with (sub-)cm accuracy, and to monitor the transformation between the systems with a high time resolution and (sub-)mas accuracy. This development led to the replacement of the ILS through the IERS, the International Earth Rotation Service in 1989.In the pre-space era, the Earth's gravity field could “only” be established by terrestrial methods. The determination of the Earth's gravitational field was revolutionized twice in the space era, first by observing geodetic satellites with optical, Laser, and Doppler techniques, secondly by implementing a continuous tracking with spaceborne GPS receivers in connection with satellite gradiometry. The sequence of the satellite gravity missions CHAMP, GRACE, and GOCE allow it to name the first decade of the 21st century the “decade of gravity field determination”.The techniques to establish and monitor the geometric and gravimetric reference frames are about to reach a mature state and will be the prevailing geodetic tools of the following decades. It is our duty to work in the spirit of our forefathers by creating similarly stable organizations within IAG with the declared goal to produce the geometric and gravimetric reference frames (including their time evolution) with the best available techniques and to make accurate and consistent products available to wider Earth sciences community as a basis for meaningful research in global change. IGGOS, the Integrated Global Geodetic Observing System, is IAG's attempt to achieve these goals. It is based on the well-functioning and well-established network of IAG services.     

Contribution of satellite altimetry data in the geophysical investigation of the Red Sea Region, Egypt

Successful development of geodetic satellite missions has aroused new interest in determining global and regional gravity field based on satellite data. Satellite altimetry data enable direct determination of the geoid over sea regions. In Egypt, where land and marine geophysical data are inadequate because of rough topography and economic reasons, the use of satellite altimetry data is of special importance. The northern Red Sea region has been selected as a site for case study of the current research, after applying spectral analysis to reveal near-surface structure, the residual geoid of the studied region shows a good correlation with the known geologic features. Moreover, satellite-based gravity data enhance small-scale features and agrees well with land and marine gravity data. Thus, geoid undulation and satellite gravity data can be a complementary source of data to determine near-surface and deep structures.     

On The Determination Of The Earth's Model – The Mean Equipotential Surface

The sea surface cannot be used as reference for Major Vertical Datum definition because its deviations from the ideal equipotential surface are very large compared to rms in the observed quantities. The quasigeoid is not quite suitable as the surface representing the most accurate Earth's model without some additional conditions, because it depends on the reference field. The normal Earth's model represented by the rotational level ellipsoid can be defined by the geocentric gravitational constant, the difference in the principal Earth's inertia moments, by the angular velocity of the Earth's rotation and by the semimajor axis or by the potential (U ₀ ) on the surface of the level ellipsoid. After determining the geopotential at the gauge stations defining Vertical Datums, gravity anomalies and heights should be transformed into the unique vertical system (Major Vertical Datum). This makes it possible to apply Brovar's (1995) idea of determining the reference ellipsoid by minimizing the integral, introduced by Riemann as the Dirichlet principle, to reach a minimum rms anomalous gravity field. Since the semimajor axis depends on tidal effects, potential U ₀ should be adopted as the fourth primary fundamental geodetic constant. The equipotential surface, the actual geopotential of which is equal to U ₀ , can be adopted as reference for realizing the Major Vertical Datum.     

Gravitational Gradients at Satellite Altitudes in Global Geophysical Studies

Due to the ESA’s satellite mission GOCE launched in March 2009, gravitational gradients sampled along the orbital trajectory approximately 250 km above the Earth’s surface have become available. Since 2010, gravitational gradients have routinely been applied in geodesy for the derivation of global Earth’s gravitational models provided in terms of fully normalized coefficients in a spherical harmonic series representation of the Earth’s gravitational potential. However, in geophysics, gravitational gradients observed by spaceborne instruments have still been applied relatively seldom. This contribution describes their possible geophysical applications in structural studies where gravitational gradients observed at satellite altitudes are compared with those derived by a spectral forward modeling technique using available global models of selected Earth’s mass components as input data. In particular, GOCE gravitational gradients are interpreted in terms of a superposition principle of gravitation as combined gravitational effects generated by a homogeneous reference ellipsoid of revolution, mean topographic and ice mass density distributions, depth-dependent mass density contrasts within bathymetry and lateral mass density anomalies with sediments and crustal layers. Respective gravitational effects are one by one removed from gravitational gradients observed at approximately 250 km elevation above ground. Removing respective gravitational gradients from observed gravitational gradients gradually reveals problematic geographic areas with model deficiencies. For the full interpretation of observed gravitational gradients, deficiencies of CRUST2.0 must be corrected and effects of deeper laying mass anomalies not included in the study considered. These findings are confirmed by parameters describing spectral properties of the gravitational gradients. The methodology can be applied for validating Earth’s gravitational models and for constraining crustal models in the development phase.     

中国岩石层底面地幔对流状态的探讨

本文应用卫星测量地球重力场的球谐函数系数和地幔流体力学方程,计算了中国地区岩石层下的地幔对流状态。计算中选用了三种不同的模式:模式Ⅰ,完全场（2-30阶）;模式Ⅱ,高阶场（13-25阶）;模式ⅢI,低阶场（4-18阶）。其目的是探求一种与我国地质和地球物理实测资料相符合的最佳模式。考虑到岩石层厚度与所假定的地幔对流区厚度相比是很薄的,因而可以把地幔流对岩石层底面产生的拖曳力近似地看作岩石层内构造力。将上述三个模式分别与中国及邻区的131个地震（1961-1977,M≥5）震源机制解和近代地质构造形变等资料进行比较。结果表明,模式Ⅰ与我国大多数地区的实测资料符合得最好,该模式所显示的地幔对流特征与相应各地区近代构造运动特点相符。模式Ⅰ为探讨近代构造运动的力源问题,提供了新的线索。     

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.     

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.     

Global,Regional and National Geodetic Reference Frames for Geodesy and Geodynamics

In July 2003 the International Association of Geodesy (IAG) established the Global Geodetic Observing System (GGOS). The GGOS is integrating the three basic components: geometry, the earth rotation and gravity. The backbone of this integration is the existing global ground network, based on the geodetic space techniques: very long baseline interferometry, satellite laser ranging, global navigation satellite systems and Doppler orbitography and radiopositioning integrated by satellite. These techniques have to operate as one global entity and in one global reference frame. The global reference frame in the GGOS is a realization of the International Terrestrial Reference System (ITRS). The ITRS is a world spatial reference system co-rotating with the Earth in its diurnal motion in the space. The IAG Subcommision for the European Reference Frame (EUREF) in 1991 recommended that the terrestrial reference system for Europe should be coincident with ITRS at the epoch t ₀ = 1989.0 and fixed to the stable part of the Eurasian Plate. It was named the European Terrestrial Reference System 89 (ETRS89). On the 2nd of June 2008, the Head Office of Geodesy and Cartography in Poland commenced operating the ASG-EUPOS multifunctional precise satellite positioning system. The ASG-EUPOS network defines the European Terrestrial Reference System ETRS89 in Poland. A close connection between the ASG-EUPOS stations and 15 out of 18 Polish EUREF permanent network stations controls the realization of the ETRS89 on Polish territory. This paper is a review of the global ITRS, as well as a regional and a national geodetic reference systems ETRS89.     

Use of data of satellite system for gravity recovery and climate experiment (GRACE) for studying and assessment of hydrological-geohydrological characteristics of large river basins

Analysis is given of a relatively new method for studying the distribution and dynamics of land water resources, which is based on measuring the anomalies of Earth gravity field with the use of GRACE satellite system. The international satellite experiment on determining the gravity and climate changes has been carried out since 2002 with the aim of making high-frequency (within the frequency range of 10900?C36000 GHz) measurements of time variations in Earth??s gravity field. The measurement method and the procedure of data processing and evaluating the hydrological-geohydrological characteristics of large river basins and regions based on GRACE data.     

Glacio-geophysical implication of an NNSS positioning in and around Syowa Station, East Antarctica

Utilization of the two-wave NNSS receiver drastically improved the positioning accuracy on the ice sheet of Antarctica. The NNSS positioning gives us 3 m three-dimensional convergence with 25 accepted satellite passes, and is most useful for the measurement of ice flow velocity. The flow velocity vectors along Route S-H-Z on Mizuho plateau, East Antarctica were obtained by estimating positional change of glaciological traverse stations after 7 years' interval. The transformation of the coordinate system was necessary before the comparison of the positioning in 1973 by the geodetic traverse method on the Bessel Reference Ellipsoid with that in 1980 by the satellite Doppler method on the NWL-8E Ellipsoid. NNSS receiving experiments at astronomical datum point in Syowa Station enabled us to estimate the translation of the coordinate origin of the Bessel Reference Ellipsoid against the geocentric NWL-8E Ellipsoid and correction terms for coordinate transformation. The obtained flow velocity is 15 m/a at H17 of around 1000 m a.s.l. and 70 m/a at Z2 of around 2000 m a.s.l. with ±6% uncertainty, and is too large to be explained only by the laminar flow of the ice sheet. The obtained velocity vectors are found to be mostly parallel to the maximum slope of the free-air gravity anomaly contours and can be interpreted as the ice sheet sliding down the slope of the subglacial mound of 2400 m relative height from the average subglacial bedrock topography. If such bedrock sliding occurrs over the whole region of Mizuho Plateau, the related thinning of the ice sheet may be detected by the precise measurement of the height change of the same marker station. By dynamically modelling the ice sheet and substituting the observed parameters such as precipitation, principal strain rate, etc., into the equation of ice thickness change the submergence velocity of around −1 m/a is expected and will be detected by carefully designed repetitive NNSS receiving experiments after several year's interval.     

武汉台重力潮汐长期观测结果

采用武汉台超导重力仪(SG C032)14年多的长期连续观测资料,研究了固体地球对二阶和三阶引潮力的响应特征,精密测定了重力潮汐参数,系统研究了最新的固体潮模型和海潮模型在中国大陆的有效性.采用最新的8个全球海潮模型计算了海潮负荷效应,从武汉台SG C032的观测中成功分离出63个2阶潮汐波群和15个3阶潮汐波群信号,3阶潮波涵盖了周日、半日和1/3日三个频段.重力潮汐观测的精度非常高,标准偏差达到1.116 nm·s^-2,系统反映了非流体静力平衡、非弹性地球对2阶和3阶引潮力的响应特征.结果表明,现有的武汉国际重力潮汐基准在半日频段非常精确,但在周日频段存在比较明显的偏差,需要进一步精化.对于中国大陆的大地测量观测,固体潮可以采用Dehant等考虑地球内部介质非弹性和非流体静力平衡建立的固体潮理论模型或Xu 等基于全球SG观测建立的重力潮汐全球实验模型作为参考和改正模型,海潮负荷效应应该采用Nao99作为改正模型.     

卫星重力梯度测量与地球引力场的精度研究

本文根据地球引力位的球谐函数展开式,利用重力梯度张量各分量导出了位系数模型的精度估计公式.从三方面进行了研究:假定卫星重力梯度仪测量精度,探讨用重力梯度数据确定地球重力场模型的精度;求出位系数模型和大气阻力引起的重力梯度卫星的轨道误差;最后,反求轨道误差和位系数误差对重力梯度测量值的影响.数值计算表明,与地面技术和常规卫星方法相比,卫星梯度测量可使重力场模型的精度至少提高3-5倍;利用重力梯度张量全分量求得的重力值精度比单用径向分量V_rr的结果提高40％以上;若仅顾及位系数模型和大气阻力误差,则轨道误差对梯度测量值的影响△V_i3（i=3,2,1）至少可分别在1/4和1/3弧圈内达到△V_i3≤σ（仪器精度）.     

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.