A contrastive study on the influences of radial and three-dimensional satellite gravity gradiometry on the accuracy of the Earth's gravitational field recovery

The accuracy of the Earth＇s gravitational field measured from the gravity field and steady-state ocean circulation explorer（GOCE）,up to 250 degrees,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij from the satellite gravity gradiometry（SGG） are contrastively demonstrated based on the analytical error model and numerical simulation,respectively.Firstly,the new analytical error model of the cumulative geoid height,influenced by the radial gravity gradient V zz and three-dimensional gravity gradient V ij are established,respectively.In 250 degrees,the GOCE cumulative geoid height error measured by the radial gravity gradient V zz is about 2 1/2 times higher than that measured by the three-dimensional gravity gradient V ij.Secondly,the Earth＇s gravitational field from GOCE completely up to 250 degrees is recovered using the radial gravity gradient V zz and three-dimensional gravity gradient V ij by numerical simulation,respectively.The study results show that when the measurement error of the gravity gradient is 3×10 12 /s 2,the cumulative geoid height errors using the radial gravity gradient V zz and three-dimensional gravity gradient V ij are 12.319 cm and 9.295 cm at 250 degrees,respectively.The accuracy of the cumulative geoid height using the three-dimensional gravity gradient V ij is improved by 30%-40% on average compared with that using the radial gravity gradient V zz in 250 degrees.Finally,by mutual verification of the analytical error model and numerical simulation,the orders of magnitude from the accuracies of the Earth＇s gravitational field recovery make no substantial differences based on the radial and three-dimensional gravity gradients,respectively.Therefore,it is feasible to develop in advance a radial cold-atom interferometric gradiometer with a measurement accuracy of 10 13 /s 2-10 15 /s 2 for precisely producing the next-generation GOCE Follow-On Earth gravity field model with a high spatial resolution.

Dependence of Gravity Induced Absorption Changes on the Earth’s Magnetic Field as Measured during Parabolic Flight Campaigns

Various spectroscopic experiments performed on the AIRBUS ZERO G—located in Bordeaux, France—in the years 2002 to 2012 exhibit minute optical reflection/absorption changes (GIACs) as a result of gravitational changes between 0 and 1.8 g in various biological species such as maize, oats, Arabidopsis and particularly Phycomyces sporangiophores. During a flight day, the AIRBUS ZERO G conducts 31 parabolas, each of which lasts about three minutes including a period of 22 s of weightlessness. So far, we participated in 11 parabolic flight campaigns including more than 1000 parabolas performing various kinds of experiments. During our campaigns, we observed an unexplainable variability of the measuring signals (GIACs). Using GPS-positioning systems and three dimensional magnetic field sensors, these finally were traced back to the changing earth’s magnetic field associated with the various flight directions. This is the first time that the interaction of gravity and the Earth’ magnetic field in the primary induction process in living system has been observed.

Discovering Crustal Deformation Bands by Processing Regional Gravity Field

Objectives： This article presents a new computational procedure to discover scratches buried in the earth＇s crust. We also validate this new interdisciplinary analysis method with regional gravity data located in a well-known Dabie orogenic zone for test. Methods： Based on the scratch analysis method evolved with mathematical morphology of surfaces, we present a procedure that extracts information of the crustal scratches from regional gravity data. Because the crustal scratches are positively and highly correlated to crustal deformation bands, it can be used for delineation of the crustal deformation belts. The scratches can be quantitatively characterized by calculation of the ridge coefficient function, whose high value traces delineate the deformation bands hidden in the regional gravity field. In addition, because the degree of crustal deformation is an important indicator of tectonic unit divisions, so the crust can be further divided according to the degree of crustal deformation into some tectonic units by using the ridge coefficient data, providing an objective base map for earth scientists to build tectonic models with quantitative evidence. Results： After the ridge coefficients are calculated, we can further enhance the boundary of high ridge-coefficient blocks, resulting in the so-called ridge-edge coefficient function. The high-value ridge-edge coefficients are well correlated with the edge faults of tectonic units underlay, providing accurate positioning of the base map for compilation of regional tectonic maps. In order to validate this new interdisciplinary analysis method, we select the Dabie orogenic zone as a pilot area for test, where rock outcrops are well exposed on the surface and detailed geological and geophysical surveys have been carried out. Tests show that the deformation bands and the tectonic units, which are conformed by tectonic scientists based on surface observations, are clearly displayed on the ridge and ridge-edge coefficient images obtained in this article. Moreover, these computer-generated images provide more accurate locations and geometric details. Conclusions： This work demonstrates that application of modern mathematical tools can promote the quantitative degree in research of modern geosciences, helping to open a door to develop a new branch of mathematical tectonics.

Efficient and rapid accuracy estimation of the Earth's gravitational field from next-generation GOCE Follow-On by the analytical method

Firstly, a new analytical error model of the cumulative geoid height using the three-dimensional diagonal tensors of satellite gravity gradiometry （SGG） is introduced based on the variance-covariance matrix principle. Secondly, a study for the requirements demonstration on the next-generation GOCE Follow-On satellite gravity gradiometry system is developed using different satellite orbital altitudes and measurement accuracies of satellite gravity gradiometer by the new analytical error model of SGG. The research results show that it is preferable to design satellite orbital altitudes of 300 km–400km and choose the measurement accuracies of 10-13/s2 –10-15/s2 from satellite gravity gradiometer. Finally, the complementarity of the four-stage satellite gravity missions, including past CHAMP, current GRACE, and GOCE, and next-generation GOCE Follow-On, is contrastively demonstrated for precisely recovering the Earth’s full-frequency gravitational field with high spatial resolution.

Construction of Earth's gravitational field model from CHAMP,GRACE and GOCE data

The basic principle of spectral combination method is discussed,and the general expressions of the spectral weight and spectral combination of the united-processing of various types of gravimetric data are shown.What＇s more,based on degree error RMS of potential coefficients,the detailed expressions of spectral combination formulae and the corresponding spectral weights in the Earth＇s gravitational field model（EGM） determination using GOCE ＋ GRACE and CHAMP ＋ GRACE ＋ GOCE are derived.The fundamental situation that ulux-champ2013 s,tongji-GRACE01,go-cons-gcf-2-tim-r5 constructed respectively by CHAMP,GRACE,GOCE data and go-cons-gcf-2-dir-r5 constructed by syncretic processing of GRACE,GOCE and LAGEOS data are explained briefly,the degree error RMS,cumulative geoid height error and cumulative gravity anomaly error of these models are calculated.A syncretic model constructed from CHAMP,GRACE and GOCE data,which is expressed by champ ＋ grace ＋ goce,is obtained based on spectral combination method.Experimentation results show that the precision of CHAMP data model is the lowest in satellite-only models,so it is not needed in the determination of syncretic models.The GRACE data model can improve the GOCE data model in medium-long wavelength,so the overall precision of syncretic model can be improved.Consequently,as many types of gravimetric data as possible should be combined together in the data processing in order to strengthen the quality and reliability with widening scope and improve the precision and spatial resolution of the computational results.

A GOCE only gravity model GOSG01S and the validation of GOCE related satellite gravity models

We compile the GOCE-only satellite model GOSG01S complete to spherical harmonic degree of 220 using Satellite Gravity Gradiometry （SGG） data and the Satellite-to-Satellite Tracking （SST） observations along the GOCE orbit based on applying a least-squares analysis. The diagonal components （Vxx, Vyy, Vzz） of the gravitational gradient tensor are used to form the system of observation equations with the band-pass ARMA filter. The point-wise acceleration observations （ax, ay, az） along the orbit are used to form the system of observation equations up to the maximum spherical harmonic degree/order 130. The analysis of spectral accuracy characteristics of the newly derived gravitational model GOSG01S and the existing models GOTIM04S, GODIR04S, GOSPW04S and JYY_GOCE02S based on their comparison with the ultrahigh degree model EIGEN-6C2 reveals a significant consistency at the spectral window approximately between 80 and 190 due to the same period SGG data used to compile these models. The GOCE related satellite gravity models GOSG01S, GOTIM05S, GODIR05S, GOTIM04S, GODIR04S, GOSPW04S, JYY_- GOCE02S, EIGEN-6C2 and EGM2008 are also validated by using GPS-leveling data in China and USA. According to the truncation at degree 200, the statistic results show that all GGMs have very similar differences at GPS-leveling points in USA, and all GOCE related gravity models have better performance than EGM2008 in China. This suggests that all these models provide much more information on the gravity field than EGM2008 in areas with low terrestrial gravity coverage. And STDs of height anomaly differences in China for the selected truncation degrees show that GOCE has improved the accuracy of the global models beyond degree 90 and the accuracies of the models improve from 24 cm to 16 cm. STDs of geoid height differences in USA show that GOSG01S model has best consistency comparing with GPSleveling data for the frequency band of the degree between 20 and 160.

Gravity gradient distribution in China's Mainland from GOCE satellite gravity gradiometry data

At present, gravity field and steady-state ocean circulation explorer（GOCE） gravity data are always used to compute regional gravity anomaly and geoid height. In this study, the latest GOCE gravity field model data（from Oct. 2009 to Jul. 2010） are used to compute the gravity gradient of China's Mainland according to a rigorous recursion formula（in all the six directions）. The results show that the numerical values of the gravity gradients are larger in the T rr direction than those in the other directions. They reflect the terrain characteristics in detail and correlate with the regional tectonics; however, in the T ql and T r l directions,the numerical values are relatively smaller and the gravity gradients in the T r l direction do not reflect the terrain characteristics in detail.

GPS水准方法比较与分析

把GPS水准方法分为3类:数值逼近法、重力法、综合法。在此基础上分别总结了3类方法各自的方法原理,运用特点。最后通过实例对各类方法用于GPS高程转换时的精度进行了比较分析。

Chinese Gravimetry Augment and Mass Change Exploring Mission Status and Future

The satellite gravimetry technology effectively recovers the global Earth’s gravity field.Since 2000s,HL-SST satellite CHAMP,LL-SST satellite GRACE,Gravity Gradient Measurement(GGM)satellite GOCE have been launched successfully,producing some Earth’s gravity models solely from satellites data.However,the space and time resolution of the Earth’s gravity fields do not adequately satisfy scientific objectives.The main reason is that the gravimetry satellites are not enough and observation data insufficient.The paper outlines the current and future status of Chinese gravity satellite missions.The Chinese gravimetry satellite system,named Chinese Gravimetry augment and Mass change exploring mission(ChiGaM),successfully launched in Dec.2021 after four years of production and over a year of calibration and valiation.The accelerometer,K-band ranging system and the three stellar sensors,among others,were precisely calibrated and trimmed.The satellite mass center was determined and coordinated with the proof center of accelerometer with an accuracy 100μm.The inter-satellite ranging system and BDS/GPS receiver operate together seamlessly.The range and range rate noise is less than 3μm/Hz^(1/2) and 1μm/s/Hz^(1/2),respectively,in band of 0.025~0.1 Hz.The electrostatic suspension accelerometer is working well.Its high-sensitive axis noise level is 3×10^(-10) m/s^(2)/Hz^(1/2)near the frequency 0.1 Hz,and 1×10^(-9) m/s^(2)/Hz^(1/2) for the less-sensitive axis.Meanwhile the BDS/GPS receiver is used to achieve data for precise orbit determination,yielding an orbit result with accuracy better than 2 cm.When compared with KBR observations,the RMS of the bias is less than 1 mm.Besides above mission,next gravimetric satellite is being developed now.TQ-2 mission is designed as a totally experimental satellite for gravitational wave detection at low Earth orbit,which can detect the Earth’s gravity field simultaneously.The Bender-type mission is considered the most promising configuration for TQ-2 and consists of two polar satellites and two inclined satellites.Due to the extra observations at the east-west direction derived from the inclined satellite pair,significant improvements has been made in detecting temporal signals with higher accuracy and spatial-temporal resolution.To achieve the scientific goal,the ACC MBW can shift from 0.001~0.1 Hz to 0.004~0.1 Hz for ACC,and the LRI MBW can shift from 0.01~1 Hz to 0.1~1 Hz.For future research,a gravimetric potential survey using cold-atomic-clock based on the general relativity theory,cold atom gradiometer should be pursued.Gravimetric technologies should be mined and researched,and the gravimetry satellite constellation should be developed,so as to improve the time resolution and space resolution for meeting the requirements of geophysics,geodesy,earthquake,water resources environment,oceanography,etc.

Modeling of Earth’s Gravity Fields Visualization Based on Quad Tree

The problems of the earth's gravity fields' visualization are both focus and puzzle currently. Aiming at multiresolution rendering, modeling of the Earth's gravity fields' data is discussed in the paper by using LOD algorithm based on Quad Tree. First, this paper employed the method of LOD based on Quad Tree to divide up the regional gravity anomaly data, introduced the combined node evaluation system that was composed of viewpoint related and roughness related systems, and then eliminated the T-cracks that appeared among the gravity anomaly data grids with different resolutions. The test results demonstrated that the gravity anomaly data grids' rendering effects were living, and the computational power was low. Therefore, the proposed algorithm was a suitable method for modeling the gravity anomaly data and has potential applications in visualization of the earth's gravity fields.

A rapid compensation method for launch data of long-range rockets under influence of the Earth＇s disturbing gravity field

Regarding the rapid compensation of the influence of the Earth＇ s disturbing gravity field upon trajectory calculation,the key point lies in how to derive the analytical solutions to the partial derivatives of the state of burnout point with respect to the launch data.In view of this,this paper mainly expounds on two issues：one is based on the approximate analytical solution to the motion equation for the vacuum flight section of a long-range rocket,deriving the analytical solutions to the partial derivatives of the state of burnout point with respect to the changing rate of the finalstage pitch program;the other is based on the initial positioning and orientation error propagation mechanism,proposing the analytical calculation formula for the partial derivatives of the state of burnout point with respect to the launch azimuth.The calculation results of correction data are simulated and verified under different circumstances.The simulation results are as follows：（1） the accuracy of approximation between the analytical solutions and the results attained via the difference method is higher than 90%,and the ratio of calculation time between them is lower than 0.2%,thus demonstrating the accuracy of calculation of data corrections and advantages in calculation speed;（2） after the analytical solutions are compensated,the longitudinal landing deviation of the rocket is less than 20 m and the lateral landing deviation of the rocket is less than 10 m,demonstrating that the corrected data can meet the requirements for the hit accuracy of a long-range rocket.

GPS and the Search for Axions

GPS, an excellent tool for geodesy, may serve also particle physics. In the presence of Earth’s magnetic field, a GPS photon may be transformed into an axion. The proposed experimental setup involves the transmission of a GPS signal from a satellite to another satellite, both in low orbit around the Earth. To increase the accuracy of the experiment, we evaluate the influence of Earth’s gravitational field on the whole quantum phenomenon. There is a significant advantage in our proposal. While the geomagnetic field B is low, the magnetized length L is very large, resulting into a scale (BL)2 orders of magnitude higher than existing or proposed reaches. The transformation of the GPS photons into axion particles will result in a dimming of the photons and even to a “light shining through the Earth” phenomenon.

高分辨率地球重力场模型DQM99

地球扰动位的球谐展开式表示是地球重力场模型应用最广泛的一种表示方法。目前通用的位系数模型的分辨率是 0 5° ,最高完全阶次为 360。由于全球重力数据覆盖的密度和数据精度差别较大 ,所以已知的重力场模型的实际分辨率及其精度也因地区而异 ,故重力场模型的精化是一个漫长的过程 ,局部积分谱权综合法是改善已知重力场模型的有效而简便的方法。基于重力场模型OSU91A及EGM 96,利用我国 1 5′× 1 5′及 5′× 5′重力异常 ,我们解算了更高阶次的地球重力场模型DQM99A、B、C、D四个模型 。

高分辨率区域重力场模型DQM2000

基于局部积分谱权综合法 ,选择国外重力场模型OSU91A和EGM96作为基础模型 ,利用中国地域2 0′× 2 0′、1 5′× 1 5′、1 0′× 1 0′和 5′× 5′平均重力异常 ,解算了一组高分辨率地球重力场系列模型DQM2 0 0 0 ,即DQM2 0 0 0A、DQM2 0 0 0B、DQM2 0 0 0C和DQM2 0 0 0D ,其完整阶次分别为 5 40、72 0、1 0 80和 2 1 60 ,采用内部符合和外部检核的方法以评价新模型的精度。

卫星重力梯度数据的模拟研究

推导了运用地球重力场模型计算单点、格网点以及格网平均的扰动重力梯度复组合分量的公式;提出了广义球谐函数及其定积分的新算法,并利用EGM96地球重力场模型试算了全球地区卫星轨道面上的重力梯度分量的格网平均观测值;通过对角线分量满足Laplace方程的精度,验证了该算法的有效性和实用性。

卫星重力梯度恢复地球重力场的时域法研究

由于卫星重力梯度技术能够有效补偿重力场信号随卫星高度增加快速衰减的影响,因此GOCE卫星的发射使得利用重力梯度测量数据恢复高精度重力场模型成为可能。以轨道摄动理论为基础,利用时域最小二乘法,分析重力梯度卫星的轨道高度、轨道倾角、卫星寿命(数据长度)、重力梯度测量精度以及数据采样频率等对其恢复地球重力场误差的影响。结果表明:时域最小二乘法适合于卫星重力梯度测量的误差分析和方案设计。

2009年7月22日武汉日全食重力效应研究

为研究2009年7月22日发生在长江流域的日全食对地球重力场的影响,处理和分析了武汉九峰国家重力基准点的GWR超导重力仪和gPhone弹簧重力仪的观测数据,结果表明:超导重力仪的观测数据与气压存在显著的负相关性;3台弹簧重力仪残差时间序列之间有较好的相关性,噪声水平在其正常精度范围之内;超导重力仪和弹簧重力仪在日全食发生期间均未观测到重力场的异常效应。

弹道主动段全射向扰动引力快速逼近方法

目前快速逼近方法都必须在已知主动段弹道的条件下才能进行计算,这对于弹道导弹的机动发射(指随时的打击目标的改变)重力场保障带来困难。这里提出了一种能够实现对某种标准弹道的弹道导弹主动段全射向扰动引力快速逼近方法。试验表明:在3×10-5ms-2的精度要求下,在该实验区域只需要建立6组分频多项式模型就可以实现弹道导弹主动段全射向扰动引力的快速逼近。

高低卫-卫跟踪模式主要技术指标的模拟分析与验证

基于重力场信号的频谱特性建立了SST-hl加速度计和星载GPS等有效载荷的解析误差分析模型。以CHAMP卫星的相关技术指标为例,对卫星高度、加速度计精度和GPS定轨精度等技术指标与恢复重力场的性能作了模拟分析。将模拟分析结果与EGM96和EIGEN-CHAMP03S模型进行了比较,验证了解析模拟分析模型的有效性。

清远市清城区似大地水准面的精密确定

利用GPS技术可获得精度优于1 cm的大地高,若能够获得相应点的高精度似大地水准面(高程异常),则可将大地高转换为相应精度的正常高,从而代替相应等级的水准测量。文中介绍了清远市清城区似大地水准面的计算方法,并对计算结果进行了分析和评价。结果表明,清远市清城区1 km格网似大地水准面的精度为1.3 cm。