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.

GRACE detection of the medium-to far-field coseismic gravity changes caused by the 2004 Mw9.3 Sumatra-Andaman earthquake

Large earthquakes cause observable changes in the Earth’s gravity field, which have been detected by the Gravity Recovery and Climate Experiment (GRACE). Since most previous studies focus on the detection of near-field gravity effects, this study provides the results from the medium- to far-field gravity changes caused by the 2004 Sumatra-Andaman earthquake that are recorded within GRACE monthly solutions. Utilizing a spherical-earth dislocation model we documented that large-scale signals predominate in the global field of the coseismic gravity changes caused by the earthquake. After removing the near-field effects, the coseismic gravity changes show a negative anomaly feature with an average magnitude of -0.18×10-8 m·s-2 in the region ranging ～40° from the epicenter, which is considered as the 'medium ffield' in this study. From the GRACE data released by Center for Space Research from August 2002 to December 2008, we retrieved the large-scale gravity changes smoothed with 3 000 km Gaussian ffilter. The results show that the coseismic gravity changes detected by GRACE in the medium field have an average of (-0.20±0.06)×10-8 m·s-2, which agrees with the model prediction. The detection confirms that GRACE is sensitive to large-scale medium-field coseismic gravitational effects of mega earthquakes, and also validates the spherical-earth dislocation model in the medium field from the perspective of satellite gravimetry.

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.

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.

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.

An Improved Energy Balance Approach and Its Application in CHAMP Gravity Field Recovery

An efficient method for gravity field determination from CHAMP orbits and accelerometer data is referred to as the energy balance approach. A new CHAMP gravity field recovery strategy based on the improved energy balance approach IS developed in this paper. The method simultaneously solves the spherical harmonic coefficients, daily Integration constant, scale and bias parameters. Two 60 degree and order gravitational potential models, XISM-CHAMPO1S from the classical energy balance approach, and XISM-CHAMPO2S from the improved energy balance, are determined using about one year＇s worth of CHAMP kinematic orbits from TUM and accelerometer data from GFZ. Comparisons among XISM-CHAMPO1S, XISM-CHAMPO2S, EIGEN-CGO3C, EIGEN-CHAMPO3S, EIGEN2, ENIGNIS and EGM96 are made. The results show that the XISM-CHAMPO2S model is more accurate than EGM96, EIGENIS, EIGEN2 and XISM-CHAMPO1S at the same degree and order, and has almost the same accuracy as EIGEN-CHAMPO3S.

Effects of rotation and gravity on an electro-magneto-thermoelastic medium with diffusion and voids by using the Lord-Shulman and dual-phase-lag models

The effects of rotation and gravity on an electro-magneto-thermoelastic medium with diffusion and voids in a generalized thermoplastic half-space are studied by using the Lord-Shulman (L-S) model and the dual-phase-lag (DPL) model. The analytical solutions for the displacements, stresses, temperature, diffusion concentration, and volume fraction field with different values of the magnetic field, the rotation, the gravity, and the initial stress are obtained and portrayed graphically. The results indicate that the effects of gravity, rotation, voids, diffusion, initial stress, and electromagnetic field are very pronounced on the physical properties of the material.

High-precision chaotic radial basis function neural network model:Data forecasting for the Earth electromagnetic signal before a strong earthquake

The Earth’s natural pulse electromagnetic field data consists typically of an underlying variation tendency of intensity and irregularities.The change tendency may be related to the occurrence of earthquake disasters.Forecasting of the underlying intensity trend plays an important role in the analysis of data and disaster monitoring.Combining chaos theory and the radial basis function neural network,this paper proposes a forecasting model of the chaotic radial basis function neural network to conduct underlying intensity trend forecasting by the Earth’s natural pulse electromagnetic field signal.The main strategy of this forecasting model is to obtain parameters as the basis for optimizing the radial basis function neural network and to forecast the reconstructed Earth’s natural pulse electromagnetic field data.In verification experiments,we employ the 3 and 6 days’data of two channels as training samples to forecast the 14 and 21-day Earth’s natural pulse electromagnetic field data respectively.According to the forecasting results and absolute error results,the chaotic radial basis function forecasting model can fit the fluctuation trend of the actual signal strength,effectively reduce the forecasting error compared with the traditional radial basis function model.Hence,this network may be useful for studying the characteristics of the Earth’s natural pulse electromagnetic field signal before a strong earthquake and we hope it can contribute to the electromagnetic anomaly monitoring before the earthquake.

Gravity Field Variations Associated with the Buried Geological Structures: San Marcos Fault (NE Mexico) Case Study

Gravity data are sensitive to local vertical offsets across high-angle faults, where rocks with different densities are juxtaposed. Yet high densities in some Mesozoic sedimentary rocks just above the basement may smear out the subtle gravity signatures of basement faults. At this study the gravity data processing tends to avoid ill-described “black-box” techniques. The study area is situated in the Palomas site, Cuatrociénegas region, Coahuila, NE Mexico. The San Marcos Fault is at least 300 km long and has WNW-ESE trend from the central part of Nuevo León State through Coahuila, and finally to the eastern part of Chihuahua State. Gravimetric data shows that the lowest values of free air and Bouguer anomalies are in the southern part of the area, and the highest values are in the western and central part of the area. Between these parts exists a zone of high horizontal gravity gradient. Configuration of linear elements of gravity field (gradient zones) delimited the San Marcos Fault in the San Marcos valley below thickness of recent sedimentary cover. Two density models were carried out, which showed that the Cretaceous rocks are in discordant contact with the Paleo- zoic rocks that can be related to the San Marcos Fault. The density was determinate using to Nettleton’s method, which results highlight the presence of the San Marcos Fault. Density models showed that the Quaternary sediments are in direct contact with the San Marcos Fault.

Synthetic tidal parameters for gravity over China and its neighbor area

The synthetic tidal parameters with high spatial resolution for gravity over China and its neighbor area are constructed with Earth＇s tidal model and ocean tide loading calculated using TPXO7 global ocean tide model as well as tidal data over China seas. The comparison between synthetic parameters and ones observed by spring gravimeters at some seismic network stations and Hong Kong station and one observed by super-conducting gravimeter at Wuhan station shows that the average differences in amplitude factors and phases are smaller than 0.005 and 0.5° respectively; and that the discrepancies between observational and synthetic parameters are dependent on gravimetric technique in that the synthetic parameters are in well agreement with the superconducting gravimetric observations. This also indicates that the synthetic result is a good estimation for tidal gravity, and the numerical results in the present paper not only can provide ground and space gravimetry such as absolute gravimetry with correction model of tidal gravity, but also provide effective tidal parameters over areas where no observation is carried out.

Magnetic Structure of the Earth’s Crust in the White Sea Region

The geological structure of the White Sea area and the surrounding land areas has been well studied in the framework of individual case studies. There are a number of local models of the deep structure of the Earth’s crust available. We propose a uniform assessment of deep crustal bodies responsible for long-period (regional) magnetic anomalies and consider their correlation with surface structures. The aim of the study is to build a three-dimensional magnetic model of the Earth’s crust in the White Sea region using aeromagnetic data and modeling technologies of the Integro software package. The model is formed on the basis of a digital map of the anomalous magnetic field reduced to the pole. The sources of magnetic anomalies are considered to be located in the Earth’s crust. The 3D distribution of the relative magnetic susceptibility of rocks was obtained by solving the inverse problem of the magnetic survey. To separate the magnetic sources by frequency and depth, it was necessary to continue the magnetic field of the model upward and to calculate the TDR derivatives, which determine the lateral boundaries of the sources of positive magnetic field anomalies. 2D distributions of magnetic sources of the model for vertical and horizontal sections with depths of 10, 15 and 20 km are analyzed. The correlation between the surface and deep structures of magnetic sources of the Earth’s crust in the region is shown.

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.

On the Cosmic Evolution of the Quantum Vacuum Using Two Variable G Models and Winterberg’s Thesis

We work within a Winterberg framework where space, i.e., the vacuum, consists of a two component superfluid/super-solid made up of a vast assembly (sea) of positive and negative mass Planck particles, called planckions. These material particles interact indirectly, and have very strong restoring forces keeping them a finite distance apart from each other within their respective species. Because of their mass compensating effect, the vacuum appears massless, charge-less, without pressure, net energy density or entropy. In addition, we consider two varying G models, where G, is Newton’s constant, and G-1, increases with an increase in cosmological time. We argue that there are at least two competing models for the quantum vacuum within such a framework. The first follows a strict extension of Winterberg’s model. This leads to nonsensible results, if G increases, going back in cosmological time, as the length scale inherent in such a model will not scale properly. The second model introduces a different length scale, which does scale properly, but keeps the mass of the Planck particle as, ± the Planck mass. Moreover we establish a connection between ordinary matter, dark matter, and dark energy, where all three mass densities within the Friedman equation must be interpreted as residual vacuum energies, which only surface, once aggregate matter has formed, at relatively low CMB temperatures. The symmetry of the vacuum will be shown to be broken, because of the different scaling laws, beginning with the formation of elementary particles. Much like waves on an ocean where positive and negative planckion mass densities effectively cancel each other out and form a zero vacuum energy density/zero vacuum pressure surface, these positive mass densities are very small perturbations (anomalies) about the mean. This greatly alleviates, i.e., minimizes the cosmological constant problem, a long standing problem associated with the vacuum.

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.

多类重力场模型的精度分析及联合确定GPS点正常高的方法

针对模型尺度参数ae对高程异常计算的影响进行数值分析;利用各模型位系数方差信息及区域高精度GPS/水准数据分析国际上新近公布的ITG-GRACE2010S、AIUB-GRACE03S、EIGEN-6C、GOCO02S、DIR_R3、TIM_R3、SPW_R2、gif48以及EGM2008多类重力场模型对应不同谱域的位系数精度;提出采用简单谱组合方法和加权谱组合方法来提高高程异常的计算精度。在试验区域,计算结果表明,各模型的中低阶(约135阶)位系数整体精度要优于EGM2008模型,gif48和EIGEN-6C模型的整体精度均优于同阶次的EGM2008模型,简单谱组合方法和加权谱组合方法均能提高区域模型高程异常精度。

地球重力场模型在线路GPS高程转换中的应用

在GPS高程转换中 ,为了减少对已知高程点的需求量 ,本文引用了地球重力场模型。采取的方法是 :选择一种拟合函数 ,将每一个GPS向量用地球重力场模型计算的高程异常差值作为观测值 ,已知高程点上用GPS求得的高程异常作为固定值 ,进行拟合处理。实例表明 ,当测段长度小于5 0km ,仅用两个已知高程点 ,高程转换精度可以满足

利用GOCE模拟观测反演重力场的Torus法

在介绍Torus方法反演地球重力场模型的基本原理和方法的基础上,基于圆环面上均匀分布的卫星引力梯度模拟观测值解算了200阶次的地球重力场模型,在无误差情况下,Torus方法解算模型的阶误差RMS小于10-16,验证了该方法的严密性。利用61dGOCE卫星轨道上无误差的模拟引力梯度观测值解算了200阶次的地球重力场模型,分析了格网化误差、极空白对解算精度的影响,迭代3次后,在不考虑低次系数情况下,模型的大地水准面阶误差和累积误差均较小,最大值仅为0.022mm和0.099mm。在沿轨卫星引力梯度模拟数据中加入5mE/Hz1/2的白噪声,基于Torus方法和空域最小二乘法解算了200阶次的地球重力场模型,Torus方法的精度略低于空域最小二乘法的精度,在不考虑低次项的情况下,两种方法解算模型的大地水准面阶误差最大值分别为1.58cm和1.45cm,累积误差最大值分别为6.37cm和5.55cm。但由于采用了二维快速傅里叶技术和块对角最小二乘法,极大地提高了计算效率。本文数值结果说明Torus方法是一种独立有效的方法,可用于GOCE任务海量卫星引力梯度观测值反演重力场的快速解算。

复杂艰险山区铁路GNSS高程转换技术研究

复杂艰险山区国家水准点非常稀少,如何通过联测尽量少的水准点进行GNSS高程转换并满足铁路勘测要求是亟需解决的问题。本文开展了观测高差法、地球重力场模型法和似大地水准面模型精化法进行GNSS高程转换技术研究,试验结果表明:(1)采用观测高差法和地球重力场模型法可大大减少联测已知水准点的数量,并满足复杂艰险山区铁路勘测要求;(2)采用研制的川西似大地水准面精化模型可将GNSS测量得到的国家2000椭球大地高直接转换为1985国家高程基准高程,成果精度优于10 cm。