Searching for Regional Crustal Velocity-Density Relations with the Use of 2-D Gravity Modelling – Central Europe Case

In this paper we search for a reference relation between seismic P-wave velocity V and density ρ _ref for the continental crust. Based on the results of modern seismic experiments, we compiled 2-D seismic models into a network of four, each about 1100–1400 km long, continental-scale seismic transects cutting all main tectonic units in Central Europe. The Moho depth (about 52 km beneath the TESZ in SE Poland, to about 25 km beneath the Pannonian Basin) and the crustal structure of this area are characterised by a large variation. This structural variation provides an interesting basis for gravity studies and especially for analysing the difference of the density structure between two major tectonic provinces of distinctive age difference: Precambrian and Phanerozoic. The 2-D gravity modelling applied for crustal cross-sections representing the regional structure, based on a unified gravity anomaly map of the area, allows for a stable determination of some general features of the regional reference velocity-density relation for the continental crust. In general three major seismo-petrological types of rocks can be distinguished: sediments, crystalline crust and mantle. In compacted sediments the reference velocity-density relation is well described by the Gardner or Nafe-Drake model. Calculated gravity anomalies, using unified velocity-density relation for the whole crystalline crust, well describe observed anomalies, with an average difference of 14 mGal. However, calculated gravity anomalies, using separated velocity-density relations for the crystalline crust of Precambrian and Phanerozoic Europe, describe observed anomalies better than for the entire crust, with an average difference 12 mGal. The most important feature of these relations is the large differentiation of the derivative dρ _ref /dV in the crystalline crust, being about 0.3 g s/m⁴ for Precambrian, and about 0.1 g s/m⁴ for the Phanerozoic crystalline crust. The modelling suggests a very small density value in the uppermost mantle ρ = 3.11 g/cm³ below the younger area, while for the older area it is ρ = 3.3 g/cm³.     

On the accuracy of the bathymetry-generated gravitational field quantities for a depth-dependent seawater density distribution

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.     

Using ground gravity to improve ice mass change estimation from GOCE gravity gradients in mid-west Greenland

Vertical gravity gradient anomalies from the Gravity and steady-state Ocean Circulation Explorer (GOCE) DIR-3 model have been used to determine gravity anomalies in mid-west Greenland by using Least-Squares Collocation (LSC) and the Reduced Point Mass (RPM) method. The two methods give nearly identical results. However, compared to LSC, the RPM method needs less computational time as the number of equations to be solved in LSC equals the number of observations. The advantage of the LSC, however, is the acquired error estimates. The observation periods are winter 2009 and summer 2012. In order to enhance the accuracy of the calculated gravity anomalies, ground gravity data from West Greenland is used over locations where the gravity change resulting from ice mass changes is negligible, i.e. over solid rock. In the period considered, the gravity anomaly change due to changes in ice mass varies from ?5 mGal to 4 mGal. It is negative over the outlet glacier Jacobshavn Isbræ, where the mass loss corresponds to a gravity change of approximately ?4 mGal. When using only GOCE vertical gravity gradients, the error estimates range from 5 mGal at the coast to 17 mGal over the ice sheet. Introducing the ground gravity data from West Greenland in the prediction reduces the errors to range from 2 to 10 mGal.     

Error sources and data limitations for the prediction of surface gravity: a case study using benchmarks

Gravity-based heights require gravity values at levelled benchmarks (BMs), which sometimes have to be predicted from surrounding observations. We use the Earth Gravitational Model 2008 (EGM2008) and the Australian National Gravity Database (ANGD) as examples of model and terrestrial observed data respectively to predict gravity at Australian National Levelling Network (ANLN) BMs. The aim is to quantify errors that may propagate into the predicted BM gravity values and then into gravimetric height corrections (HCs). Our results indicate that an approximate ±1 arc-min horizontal position error of the BMs causes maximum errors in EGM2008 BM gravity of ~22 mGal (~55 mm in the HC at ~2200 m elevation) and ~18 mGal for ANGD BM gravity because the values are not computed at the true location of the BM. We use RTM (residual terrain modelling) techniques to show that ~50% of EGM2008 BM gravity error in a moderately mountainous region can be accounted for by signal omission. Non-representative sampling of ANGD gravity in this region may cause errors of up to 50 mGals (~120 mm for the Helmert orthometric correction at ~2200 m elevation). For modelled gravity at BMs to be viable, levelling networks need horizontal BM positions accurate to a few metres, while RTM techniques can be used to reduce signal omission error. Unrepresentative gravity sampling in mountains can be remedied by denser and more representative re-surveys, and/or gravity can be forward modelled into regions of sparser gravity.     

Data requirements for a 5-mm quasi-geoid in the Netherlands

We assess the surface gravity data requirements for a 5-mm quasi-geoid model for the Netherlands mainland and continental shelf in terms of omission and commission errors. The omission error critically depends on the roughness of the topography and bathymetry. For the Netherlands continental shelf, Central and Northern Netherlands, the omission error is well described by the model 0.32d mm, where d is the data spacing in km. For the more hilly Southern Netherlands, the omission error model is 0.92d mm. The commission error depends on the kernel modification, the data spacing, and the data accuracy. When using the spheroidal Stokes kernel, it is well described by 0.277 d σ_Δg mm, where σ_Δg is the noise standard deviation of surface gravity data in mGal. An upper bound of the commission error of the state-of-the-art satellite-only gravity model GOCO05S over the Netherlands is e^{0.03676L–11.419} mm, where L is the maximum degree up to which this model is used. Only if this model is truncated at a sufficiently low degree, e.g., at degree 100, its contribution to the total commission error can be neglected. We determine the total error as the sum of commission and omission errors. Hence, to realize a 5-mm quasi-geoid model for the Netherlands mainland and continental shelf, a data spacing of 3.5 km is needed when assuming a noise standard deviation of 1.5 mGal for surface gravity data. The currently available land-based gravity data fulfill this requirement. This does not apply to the situation at sea, where the density of the shipboard gravity data and the accuracy of the radar altimeter-derived data do not allow the realization of a 5-mm quasi-geoid model.     

Upward continuation of Dome-C airborne gravity and comparison with GOCE gradients at orbit altitude in east Antarctica

An airborne gravity campaign was carried out at the Dome-C survey area in East Antarctica between the 17th and 22nd of January 2013, in order to provide data for an experiment to validate GOCE satellite gravity gradients. After typical filtering for airborne gravity data, the cross-over error statistics for the few crossing points are 11.3 mGal root mean square (rms) error, corresponding to an rms line error of 8.0 mGal. This number is relatively large due to the rough flight conditions, short lines and field handling procedures used. Comparison of the airborne gravity data with GOCE RL4 spherical harmonic models confirmed the quality of the airborne data and that they contain more high-frequency signal than the global models. First, the airborne gravity data were upward continued to GOCE altitude to predict gravity gradients in the local North-East-Up reference frame. In this step, the least squares collocation using the ITGGRACE2010S field to degree and order 90 as reference field, which is subtracted from both the airborne gravity and GOCE gravity gradients, was applied. Then, the predicted gradients were rotated to the gradiometer reference frame using level 1 attitude quaternion data. The validation with the airborne gravity data was limited to the accurate gradient anomalies (T_XX, T_YY, T_ZZ and T_XZ) where the long-wavelength information of the GOCE gradients has been replaced with GOCO03s signal to avoid contamination with GOCE gradient errors at these wavelengths. The comparison shows standard deviations between the predicted and GOCE gradient anomalies T_XX, T_YY, T_ZZ and T_XZ of 9.9, 11.5, 11.6 and 10.4 mE, respectively. A more precise airborne gravity survey of the southern polar gap which is not observed by GOCE would thus provide gradient predictions at a better accuracy, complementing the GOCE coverage in this region.     

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.     

Ability of the EGM2008 high degree geopotential model to calculate a local geoid model in Valencia,Eastern Spain

A new generation of global geopotential models (GGM) is being developed. These solutions offer a file with fully-normalized spherical harmonic coefficients of the Earth’s gravitational potential up to a degree greater than 2000 with very low commission errors. This paper analyses the recent Earth Gravitational Model EGM2008, developed up to degree and order 2159 with additional coefficients to degree 2190 and order 2159, which means recovering the gravitational field up to approximately 20 km wavelengths. 223 GPS/levelling points of the new Spanish High Precision Levelling Network in the Valencia region (Eastern Spain) are used as external tool for evaluation in that particular region. The same evaluation has been performed to other different global (EGM96 and EIGENCG03C), continental (EGG97), regional (IGG2005 and IBERGEO2006) and local (GCV07) geoid models for comparison purposes only. These comparisons show that EGM2008 is the geoid model that best fits to the GPS/levelling data in that region.     

捷联式航空重力测量算法比较

捷联式航空重力测量系统与平台式系统相比具有体积小、重量轻、功耗低等许多优点,近些年来取得了显著的研究进展.本文给出了捷联式航空重力测量的两种算法模型:捷联式惯性标量重力测量(SISG)和旋转不变式标量重力测量(RISG)模型,并对其误差模型作了初步讨论.利用我国首套捷联式航空重力仪SGA-WZ01在某海域的部分试验数据,对两种算法模型进行了比较分析,表明其差值之标准差对于200s的滤波长度小于0.5mGal.同时,利用两组重复测线数据估算了不同滤波尺度下的两种算法的内符合精度,表明SISG算法略优于RISG算法.对于200s和300s的滤波长度,SISG的内符合精度分别为1.06mGal和0.80mGal.     

Gravity control network at airports in nigeria

A network of gravity of control bases located at most of the airports and airfields in Nigeria has been established as a means of controlling future gravimetric studies in Nigeriaand of integrating existing gravity surveys into a unified whole for the purpose of producing gravity anomaly maps for the country. The network will also provide a means of checking and adjusting gravimeter calibration differences to a common standard.The LaCoste and Romberg gravimeters, models G446, G464 and G468, used for the observations were transported between the base stations by means of a six-seater fixed-wing aircraft. The scale adopted of the network was the same as that of gravimeter G464 to which the other two meters were referenced. The absolute gravity values observed ranged from 978184.882 ± 0.014 mGal at Nguru to 977844.379 ± 0.020 mGal at Jos, yielding a gravity range of 340.503 mGal. The accuracy of the individual absolute gravity value is better than 0.030 mGal using a datum value of 978120.933 mGal for site 103401 estlished at the Kano International Airport. The 34 observation sites are illustrated with sufficient detail to permit reoccupation within 0.50 m to 1.0 m.     

The Rock–Water–Ice Topographic Gravity Field Model RWI_TOPO_2015 and Its Comparison to a Conventional Rock-Equivalent Version

RWI_TOPO_2015 is a new high-resolution spherical harmonic representation of the Earth’s topographic gravitational potential that is based on a refined Rock–Water–Ice (RWI) approach. This method is characterized by a three-layer decomposition of the Earth’s topography with respect to its rock, water, and ice masses. To allow a rigorous separate modeling of these masses with variable density values, gravity forward modeling is performed in the space domain using tesseroid mass bodies arranged on an ellipsoidal reference surface. While the predecessor model RWI_TOPO_2012 was based on the \(5'\times 5'\) global topographic database DTM2006.0 (Digital Topographic Model 2006.0), the new RWI model uses updated height information of the \(1'\times 1'\) Earth2014 topography suite. Moreover, in the case of RWI_TOPO_2015, the representation in spherical harmonics is extended to degree and order 2190 (formerly 1800). Beside a presentation of the used formalism, the processing for RWI_TOPO_2015 is described in detail, and the characteristics of the resulting spherical harmonic coefficients are analyzed in the space and frequency domain. Furthermore, this paper focuses on a comparison of the RWI approach to the conventionally used rock-equivalent method. For this purpose, a consistent rock-equivalent version REQ_TOPO_2015 is generated, in which the heights of water and ice masses are condensed to the constant rock density. When evaluated on the surface of the GRS80 ellipsoid (Geodetic Reference System 1980), the differences of RWI_TOPO_2015 and REQ_TOPO_2015 reach maximum amplitudes of about 1 m, 50 mGal, and 20 mE in terms of height anomaly, gravity disturbance, and the radial–radial gravity gradient, respectively. Although these differences are attenuated with increasing height above the ellipsoid, significant magnitudes can even be detected in the case of the satellite altitudes of current gravity field missions. In order to assess their performance, RWI_TOPO_2015, REQ_TOPO_2015, and RWI_TOPO_2012 are validated against independent gravity information of current global geopotential models, clearly demonstrating the attained improvements in the case of the new RWI model.     

Atmospheric effects in the derivation of geoid-generated gravity anomalies

Parameters of the gravity field harmonics outside the geoid are sought in solving the Stokes boundary-value problem while harmonics outside the Earth in solving the Molodensky boundary-value problem. The gravitational field generated by the atmosphere is subtracted from the Earth’s gravity field in solving either the Stokes or Molodensky problem. The computation of the atmospheric effect on the ground gravity anomaly is of a particular interest in this study. In this paper in particular the effect of atmospheric masses is discussed for the Stokes problem. In this case the effect comprises two components, specifically the direct and secondary indirect atmospheric effects. The numerical investigation is conducted at the territory of Canada. Numerical results reveal that the complete effect of atmosphere on the ground gravity anomaly varies between 1.75 and 1.81 mGal. The error propagation indicates that precise determination of the atmospheric effect on the gravity anomaly depends mainly on the accuracy of the atmospheric mass density distribution model used for the computation.     

SUBSURFACE GRAVITY MEASUREMENTS AND DENSITY MAPPING IN THE MAILARAM COPPER MINES,ANDHRA PRADESH,INDIA*

Gravity measurements were made in the Mailaram copper mines, Andhra Pradesh. The observations were distributed between the two shafts situated about 220 m apart and in the three levels up to a maximum depth of 100 m. Assuming a normal free-air gradient, average densities for the three layers were determined as 2.631, 2.604, and 2.823 g cm^-3, respectively. Upon incorporating the weighted mean density values from measurements on samples, the free-air gradients were found to be 0.315 mGal m^-1 for the second layer (i.e. between the first and second levels) and 0.2978 mGal m^-1 for the third layer (i.e. between the second and third levels). The density variation map obtained from the gravity data, the deduced anomalies, the weighted mean density values from measurements on rock samples, and the varying free-air coefficients all suggest correspondence with the concentration of ore lodes.     

The first airborne scalar gravimetry system based on SINS/DGPS in China

China has developed an airborne gravimetry system based on SINS/DGPS named SGA-WZ,the first system in which a strapdown inertial navigation system(SINS)has been used for airborne gravimetry in China.This gravity measurement system consists of a strap-down inertial navigation system and a differential global positioning system(DGPS).In April 2010,a flight test was carried out in Shandong Province of China to test the accuracy of this system.The test was designed to assess the repeatability and accuracy of the system.Two repeated flights and six grid flights were made.The flying altitude was about 400 m.The average flying speed was about 60 m/s,which corresponds to a spatial resolution of 4.8 km when using 160-s cutoff low-pass filter.This paper describes the data processing of the system.The evaluation of the internal precision is based on repeated flights and differences in crossover points.Gravity results in this test from the repeated flight lines show that the repeatability of the repeat lines is 1.6 mGal with a spatial resolution of 4.8 km,and the internal precision of grid flight data is3.2 mGal with a spatial resolution of 4.8 km.There are some systematic errors in the gravity results,which can be modeled using trigonometric function.After the systematic errors are compensated,the precision of grid flight data can be better than1 mGal.     

Assessment of some global solar radiation parameterizations

In spite of their practicability, most classical models are not versatile but rather restrictive in their application. Consequently, their applicability for a particular location depends largely on validation against actual measurements. Global solar radiation parameterizations have been evaluated in this study for a lowland and a mountain site. Tested models were broadly categorised as cloud-based (Kasten) and sunshine-based (Ångström–Prescott, Garg and Garg, Sivkov). Data sets utilised for the evaluation extended from 1991 to 1994. Adjustable parameters in the models were determined. Observed monthly mean values of solar radiation G and those estimated using Kasten model agreed within 2.5% for the lowland site and 13% for the mountain site. Root mean square errors of estimated hourly values of G using Kasten model appreciated significantly with fractional cloud cover N (particularly for N>4 octals). For the study sites as well as other locations examined here, Ångström–Prescott coefficients did not show a distinctive trend with respect to season, geographical co-ordinate or altitude. Monthly mean values of G estimated using Ångström–Prescott model agreed with observation within 2.5% for the lowland site and 3.4% for the mountain site. The effect of air mass, latitude and water vapour terms on the Ångström–Prescott relation has also been investigated. In general, Ångström–Prescott as well as Garg and Garg models yielded the least RMSE (<0.047) for the study sites and are thus recommended for estimating G for an arbitrary location.     

Studying the Representation Accuracy of the Earth’s Gravity Field in the Polar Regions Based on the Global Geopotential Models

The development of studies on estimating the accuracy of the Earth’s modern global gravity models in terms of the spherical harmonics of the geopotential in the problematic regions of the world is discussed. The comparative analysis of the results of reconstructing quasi-geoid heights and gravity anomalies from the different models is carried out for two polar regions selected within a radius of 1000 km from the North and South poles. The analysis covers nine recently developed models, including six high-resolution models and three lower order models, including the Russian GAOP2012 model. It is shown that the modern models determine the quasi-geoid heights and gravity anomalies in the polar regions with errors of 5 to 10 to a few dozen cm and from 3 to 5 to a few dozen mGal, respectively, depending on the resolution. The accuracy of the models in the Arctic is several times higher than in the Antarctic. This is associated with the peculiarities of gravity anomalies in every particular region and with the fact that the polar part of the Antarctic has been comparatively less explored by the gravity methods than the polar Arctic.     

EGM2008地球重力模型数据在中国大陆地区的精度分析

本文介绍了5′×5 ′的EGM2008地球重力模型及其在全球的精度评价.按照地形变化规律,将中国大陆大致分为7个区域,在10 km网度上,将EGM2008地球重力模型数据与中国地面实测空间重力网格数据进行了对比.由于数据源的问题,中国大陆的模型数据精度普遍低于北美和欧洲.二种数据在地形平坦的东部地区差别较小,向西随着地形复杂程度的增加,二种数据之间的标准差从小于10 mGal增大到50多mGal.畸变点分析表明精度极低的网格点均分布在地形起伏大的地区.总体而言,5′×5′的EGM2008地球重力模型数据在中国大陆将近80％的面积上的精度可达10 mGal之内,可用于小比例尺重力编图和构造研究.在地形起伏较大的中国西部以青藏高原为例进一步比较了EGM2008重力模型和重力测点数据,结果表明在重力点分布稀疏不均匀的地区,平面网格数据难以准确表达重力场信息.由于缺少地面重力数据控制,EGM2008重力模型数据在中国西部精度较低,但模型数据依然在很大程度上提高了空间重力异常信息的丰富程度.将中国区域重力调查成果数据应用于地球模型的构建是一项有意义的工作.     

Improving airborne strapdown vector gravimetry using stabilized horizontal components

Integrating the deflections of the vertical along the flight line can yield geoid profiles which are valuable in the study of geodesy and geophysics, fortunately, the deflections can be measured directly by vector gravimetry. Airborne vector gravimetry using a Strapdown Inertial Navigation System and the Global Navigation Satellite System (SINS/GNSS) has shown promising results in previous studies. However, the quality of the SINS and GNSS is a major limitation; in particular, the attitude errors induced by the gyros will result in large measurement errors to the horizontal components of the gravity disturbance, and these measurement errors represent the behavior of low-frequency trend. An airborne vector gravimetry method used to remove the bias and low-frequency trends in the gravity disturbance estimated for each survey line has been developed. This method uses the horizontal components of the gravity disturbance computed from EGM2008 (Earth Gravitational Model 2008) as a reference. Firstly, the horizontal measurement results obtained from the gravimeter are divided into high- and low-frequency components according to the resolution of the EGM2008, and then, the bias and low-frequency trends of the low-frequency components are corrected using a linear fit to the EGM2008 reference data. Finally, the ultimate results can be acquired after combining the high-frequency components and the corrected low-frequency components. The data used was obtained from the SGA-WZ, which is the first strapdown airborne gravimeter developed in China. The results of this method are promising. The internal accuracy of the gravity disturbance's horizontal components for repeated survey lines exceeds 3.5 mGal, and the corresponding resolution is approximately 4.8 km based on 160-s data smoothing and an airplane averaging speed of approximately 216 km/h. After applying the WCF (Wavenumber Correlation Filter), the internal accuracy of the horizontal components exceeds 2 mGal. This can satisfy the requirement of the application in geodesy and solid earth geophysics.     

利用解析法有效快速估计将来GRACE Follow-On地球重力场的精度

本文首次利用解析法有效快速估计了将来GRACE（Gravity Recovery and Climate Experiment） Follow-On地球重力场的精度. 第一,基于功率谱原理分别建立了新的GRACE Follow-On卫星激光干涉星间测量系统星间速度、GPS接收机轨道位置和轨道速度以及加速度计非保守力误差影响累计大地水准面的单独和联合解析误差模型. 第二,利用提出的GRACE卫星关键载荷匹配精度指标和美国喷气推进实验室（JPL）公布的GRACE Level 1B实测精度指标的一致性,以及估计的GRACE累计大地水准面精度和德国波兹坦地学研究中心（GFZ）公布的EIGEN-GRACE02S地球重力场模型实测精度的符合性,验证了本文建立的解析误差模型是可靠的. 第三,论证了GRACE Follow-On卫星不同关键载荷匹配精度指标和轨道高度对地球重力场精度的影响. 在360阶处,利用轨道高度250 km、星间距离50 km、星间速度误差1×10^-9m/s、轨道位置误差3×10^-5m、轨道速度误差3×10^-8m/s和非保守力误差3×10^-13m/s²,基于联合解析误差模型估计累计大地水准面的精度为1.231×10^-1 m. 本文的研究不仅为当前GRACE和将来GRACE Follow-On地球重力场精度的有效快速确定提供了理论基础和计算保证,同时对国际将来GRAIL（Gravity Recovery and Interior Laboratory）月球卫星重力测量计划的成功实施具有重要的参考意义.