A High-Resolution Earth’s Gravity Field Model SGG-UGM-2 from GOCE,GRACE,Satellite Altimetry,and EGM2008

This paper focuses on estimating a new high-resolution Earth’s gravity field model named SGG-UGM-2 from satellite gravimetry,satellite altimetry,and Earth Gravitational Model 2008(EGM2008)-derived gravity data based on the theory of the ellipsoidal harmonic analysis and coefficient transformation(EHA-CT).We first derive the related formulas of the EHA-CT method,which is used for computing the spherical harmonic coefficients from grid area-mean and point gravity anomalies on the ellipsoid.The derived formulas are successfully evaluated based on numerical experiments.Then,based on the derived least-squares formulas of the EHA-CT method,we develop the new model SGG-UGM-2 up to degree 2190 and order 2159 by combining the observations of the Gravity Field and Steady-State Ocean Circulation Explorer(GOCE),the normal equation of the Gravity Recovery and Climate Experiment(GRACE),marine gravity data derived from satellite altimetry data,and EGM2008-derived continental gravity data.The coefficients of degrees 251–2159 are estimated by solving the block-diagonal form normal equations of surface gravity anomalies(including the marine gravity data).The coefficients of degrees 2–250 are determined by combining the normal equations of satellite observations and surface gravity anomalies.The variance component estimation technique is used to estimate the relative weights of different observations.Finally,global positioning system(GPS)/leveling data in the mainland of China and the United States are used to validate SGG-UGM-2 together with other models,such as European improved gravity model of the earth by new techniques(EIGEN)-6C4,GECO,EGM2008,and SGG-UGM-1(the predecessor of SGG-UGM-2).Compared to other models,the model SGG-UGM-2 shows a promising performance in the GPS/leveling validation.All GOCE-related models have similar performances both in the mainland of China and the United States,and better performances than that of EGM2008 in the mainland of China.Due to the contribution of GRACE data and the new marine gravity anomalies,SGG-UGM-2 is slightly better than SGG-UGM-1 both in the mainland of China and the United States.

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.

A Solvable Model in Two—Dimensional Gravity Coupled to a Nonlinear Matter Field

The two-dimensional gravity model with a coupling constant and a vanishing cosmological constant coupled to a nonlinear matter field is investigated. We found that the classical equations of motion are exactly solvable and the static solutions of the induced metric and scalar curvature can be obtained analytically. These solutions may be used to describe the naked singularity at the origin.

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.

WHU-Grace01s:A new temporal gravity field model recovered from GRACE KBRR data alone

A new temporal gravity field model called WHU-Grace01s solely recovered from Gravity Recovery and Climate Experiment （GRACE） K-Band Range Rate （KBRR） data based on dynamic integral approach is presented in this paper. After meticulously preprocessing of the GRACE KBRR data, the root mean square of its post residuals is about 0.2 micrometers per second, and seventy-two monthly temporal solutions truncated to degree and order 60 are computed for the period from January 2003 to December 2008. After applying the combi- nation filter in WHU-Grace01s, the global temporal signals show obvious periodical change rules in the large-scale fiver basins. In terms of the degree variance, our solution is smaller at high degrees, and shows a good consistency at the rest of degrees with the Release 05 models from Center for Space Research （CSR）, GeoForschungsZentrum Potsdam （GFZ） and Jet Pro- pulsion Laboratory 0PL）. Compared with other published models in terms of equivalent water height distribution, our solution is consistent with those published by CSR, GFZ, JPL, Delft institute of Earth Observation and Space system （DEOS）, Tongji University （Tongji）, Institute of Theoretical Geodesy （ITG）, Astronomical Institute in University of Bern （AIUB） and Groupe de Recherche de Geodesie Spatiale （GRGS}, which indicates that the accuracy of WHU-Grace01s has a good consistency with the previously published GRACE solutions.

Time-varying gravity field model of Sichuan-Yunnan region based on the equivalent mass source model

High-precision time-varying gravity field is an effective way to study the internal mass movement and understanding the spatio-temporal evolution process of the geodynamic system.Compared to the satellite gravity measurement,the repeated terrestrial gravity observation can provide a more high-order signal related to the shallow crust and subsurface.However,the suitable and unified method for gravity model estimation is a key problem for further applications.In this study,we introduce the spherical hexahedron element to simulate the field source mass and forward model the change of gravity field located at the Sichuan-Yunnan region(99—104°E,23—29°N)in the four epochs from 2015 to 2017.Compared to the experimental results based on Slepian or spherical harmonics frequency domain method,this alternative approach is suitable for constructing the equivalent mass source model of regional-scale gravity data,by introducing the first-order smooth prior condition of gravity time-varying signal to suppress the high-frequency component of the signal.The results can provide a higher spatial resolution reference for regional gravity field modeling in the Sichuan-Yunnan region.

The construction of high precision geostrophic currents based on new gravity models of GOCE and satellite altimetry data

The new gravity field models of gravity field and steady-state ocean circulation explorer(GOCE),TIM_R6 and DIR_R6,were released by the European Space Agency(ESA)in June 2019.The sixth generation of gravity models have the highest possible signal and lowest error levels compared with other GOCE-only gravity models,and the accuracy is significantly improved.This is an opportunity to build high precision geostrophic currents.The mean dynamic topography and geostrophic currents have been calculated by the 5th(TIM_R5 and DIR_R5),6th(TIM_R6 and DIR_R6)release of GOCE gravity field models and ITSG-Grace2018 of GRACE gravity field model in this study.By comparison with the drifter results,the optimal filtering lengths of them have been obtained(for DIR_R5,DIR_R6,TIM_R5 and TIM_R6 models are 1°and for ITSG-Grace2018 model is 1.1°).The filtered results show that the geostrophic currents obtained by the GOCE gravity field models can better reflect detailed characteristics of ocean currents.The total geostrophic speed based on the TIM_R6 model is similar to the result of the DIR_R6 model with standard deviation(STD)of 0.320 m/s and 0.321 m/s,respectively.The STD of the total velocities are 0.333 m/s and 0.325 m/s for DIR_R5 and TIM_R5.When compared with ITSG-Grace2018 results,the STD(0.344 m/s)of total geostrophic speeds is larger than GOCE results,and the accuracy of geostrophic currents obtained by ITSG-Grace2018 is lower.And the absolute errors are mainly distributed in the areas with faster speeds,such as the Antarctic circumpolar circulation,equatorial region,Kuroshio and Gulf Stream areas.After the remove-restore technique was applied to TIM_R6 MDT,the STD of total geostrophic speeds dropped to 0.162 m/s.

Analysis of the Processes of Gravity in the Framework of Curvature of Space and the Substantiation of the New Model

The paper belongs to the sphere of quantum physics, physics of waves and physical fields, in particular—to the gravitation. Their study provides a better understanding of the problems of natural sciences at all levels, from elementary particles, to Universe as a whole. Therefore, the solution of these problems is an urgent and important task, which to the works of many generations of scientists of the world was dedicated. However, they have not been fully resolved. In well-known works, including general relativity, determination of the wave and energy parameters of the gravitational field of the Universe and their numerical values are absent. Solutions found are limited to tensor equations of a general form, which allows their interpretation of over a wide range. Other disadvantages of famous models are: 1) the voluminous world of the Universe reduced to the planes on which space objects and other objects move, sagging planes due to their own mass;2) signs of “top” and “bottom” of the system, which are not in the real Universe, just as they are not on Earth and not in the Solar system;3) the formation of “voids” between the object and the curved space and others. Main goals of the work to identify these contradictions and find ways to resolve them are performed. The main difference and the scientific novelty of the work performed are the justification of the gravity model based on a rigorous determination of the wave and energy parameters of the gravitational field of the Universe and their numerical values. The initial parameters of this worked—is the frequency oscillation νG of the waves of the gravitational field (Nastasenko’s constant) found in 2011. Research Results: Knowing νG can find all wave parameters of the gravitational field and their numerical values. The proposed new spatial-wave model of the action of gravity is based on the wave parameters of the gravitational fields of material objects. In the framework of their unity with electromagnetic fields, it reduces their structures to similar ones and eliminates the drawbacks of the previous model—of replaced gravity on curvature of space.

Gravity field recovery from GOCE orbits using the energy conservation approach

The sum of the dissipative energy and energy constant of the GOCE satellite is found by a priori gravity field model at first, and the GOCE dissipative energy is obtained by computing the adjacent epoch difference via the differential method. Then, a gravity field model GOCE-ECPO1, which up to the degree and order 80, is recovered by the energy conservation approach from the 103-day precise orbital data of the GOCE satellite collected from November 1, 2009 to January 12, 2010. Finally, the model is compared with existing models EGM96, ITG-CHAMP05S, EIGEN-GRACE2010S, EIGEN-6C and GO CONS GCF 2 DIR_R3. The results show that at the same order and degree, the accuracy of model GOCE-EBPO1 is higher than those of models EGM96 and ITG-CHAMP05S, but lower than those of models EIGEN-GRACE2010S, EIGEN-6C and GO_CONS_GCF 2 DIR_R3, which is mainly caused by the pole gap.

Monthly gravity field solution from GRACE range measurements using modified short arc approach

In this paper we present a series of monthly gravity field solutions from Gravity Recovery and Climate Experiment（GRACE） range measurements using modified short arc approach,in which the ambiguity of range measurements is eliminated via differentiating two adjacent range measurements.The data used for developing our monthly gravity field model are same as Tongji-GRACEOl model except that the range measurements are used to replace the range rate measurements,and our model is truncated to degree and order 60,spanning Jan.2004 to Dec.2010 also same as Tongji-GRACE01 model.Based on the comparison results of the C_（2,0）,C_（2,1）,S_（2,1）,and C_（15,15）,S_（15,15）,time series and the global mass change signals as well as the mass change time series in Amazon area of our model with those of Tongji-GRACE01 model,we can conclude that our monthly gravity field model is comparable with Tongji-GRACE01 monthly model.

Analyses on gravity variation before and after the Lijiang earthquake based on a finite rectangular dislocation model

The methods were discussed to calculate the gravity variation due to crustal deformation based on a model of dis-location on a finite rectangular plane. Taking the Lijiang MS=7.0 earthquake as an example the calculating princi-ple of fault parameters were determined, and the results were given. Of particular interests were the characteristics of the gravity variations in different dislocation types. With comparison between the calculated results and the practical measurements, it was found that the model could to some extent account for the observations. But it failed to give explanations to the more far spatial gravity variation.

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.

Monthly gravity field recovery from GRACE orbits and K-band measurements using variational equations approach

The Gravity Recovery and Climate Experiment（GRACE） mission can significantly improve our knowledge of the temporal variability of the Earth＇s gravity field.We obtained monthly gravity field solutions based on variational equations approach from GPS-derived positions of GRACE satellites and K-band range-rate measurements.The impact of different fixed data weighting ratios in temporal gravity field recovery while combining the two types of data was investigated for the purpose of deriving the best combined solution.The monthly gravity field solution obtained through above procedures was named as the Institute of Geodesy and Geophysics（IGG） temporal gravity field models.IGG temporal gravity field models were compared with GRACE Release05（RL05） products in following aspects：（i） the trend of the mass anomaly in China and its nearby regions within 2005-2010; （ii） the root mean squares of the global mass anomaly during 2005-2010; （iii） time-series changes in the mean water storage in the region of the Amazon Basin and the Sahara Desert between 2005 and 2010.The results showed that IGG solutions were almost consistent with GRACE RL05 products in above aspects（i）-（iii）.Changes in the annual amplitude of mean water storage in the Amazon Basin were 14.7 ± 1.2 cm for IGG,17.1 ± 1.3 cm for the Centre for Space Research（CSR）,16.4 ± 0.9 cm for the GeoForschungsZentrum（GFZ） and 16.9 ± 1.2 cm for the Jet Propulsion Laboratory（JPL） in terms of equivalent water height（EWH）,respectively.The root mean squares of the mean mass anomaly in Sahara were 1.2 cm,0.9 cm,0.9 cm and 1.2 cm for temporal gravity field models of IGG,CSR,GFZ and JPL,respectively.Comparison suggested that IGG temporal gravity field solutions were at the same accuracy level with the latest temporal gravity field solutions published by CSR,GFZ and JPL.

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.

Determination of Antarctic geoid by using global gravity field

With Chinese latest global gravity field model WDM94, the authors provide the geoid height and mean free air gravity anomaly of Antarctica (The range of latitude is from -60° to -90°). In order to conclude and analyze the characters of Antarctic geoid roundly, the authors collect the latest oversea global gravity field model OSU91 (to degree and order 360) and JGMOSU (to degree and order 360), get the corresponding geoid height and mean free air gravity anomaly. The results are compared with the results got from WDM94, thus we get the difference. The standard deviation of geoid height between WDM94 and OSU91 is ±1.90 m;the deviation of geoid between WDM94 and JGMOSU is ±2.09 m. The standard deviation of mean gravity anomaly are ±8.97 mGal and ±9.32 mGal respectively. [WT9.HZ]

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.

Forward Modeling Research of Gravity Data Offset in North China Region

Due to the complexity and overlap effects of gravity anomalies,the gravity high point always deviates from the real position of the target.Researching the offset characteristics is helpful to determine the right location of gravity anomalies.In this paper,a series of forward models were designed to calculate the offset features under different assumptions.Different density models were established based on the measured gravity data in north China Sunhu area and a variety of methods were tried to eliminate offset effects under different conditions.The results indicate that the gravity anomalies of different density bodies can be separated effectively and the real position of target can be correctly located by using the“layer stripping method”.

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.

A data-adaptive network design for the regional gravity field modelling using spherical radial basis functions

A high-precision regional gravity field model is significant in various geodesy applications.In the field of modelling regional gravity fields,the spherical radial basis functions(SRBFs)approach has recently gained widespread attention,while the modelling precision is primarily influenced by the base function network.In this study,we propose a method for constructing a data-adaptive network of SRBFs using a modified Hierarchical Density-Based Spatial Clustering of Applications with Noise(HDBSCAN)algorithm,and the performance of the algorithm is verified by the observed gravity data in the Auvergne area.Furthermore,the turning point method is used to optimize the bandwidth of the basis function spectrum,which satisfies the demand for both high-precision gravity field and quasi-geoid modelling simultaneously.Numerical experimental results indicate that our algorithm has an accuracy of about 1.58 mGal in constructing the gravity field model and about 0.03 m in the regional quasi-geoid model.Compared to the existing methods,the number of SRBFs used for modelling has been reduced by 15.8%,and the time cost to determine the centre positions of SRBFs has been saved by 12.5%.Hence,the modified HDBSCAN algorithm presented here is a suitable design method for constructing the SRBF data adaptive network.