Although satellite gravity gradient data plays a great role in determining short-wavelength part of static gravity field model,accuracy of the long-wavelength part of gravity field model recovered by them are poor,whi...Although satellite gravity gradient data plays a great role in determining short-wavelength part of static gravity field model,accuracy of the long-wavelength part of gravity field model recovered by them are poor,which leads to only a few applications in time-variable gravity investigation.The reason is that some factors limit the accuracy of the gravity field recovered using gradient data,including accuracy of the gravity gradient observations,measurement bandwidth(MBW)of gradiometer,satellite inclination,etc.This paper aims at analyzing the influence of these limitations on gravity field recovery and discusses the possibility of time-variable gravity field detection by using gravity gradient observations.Firstly,for arbitrary satellite orbit inclination,we give the frequency distributions of all the components of gravity gradients(i.e.Txx;Tyy;Tzz;Txy;Txz and Tyz,).The results show that the maximum frequency of each component of the gravity gradients is the same,i.e.l=Ts(l is degree of the gravity field model,Ts is the orbital periods),and it is not influenced by the inclination of the satellite orbits.Secondly,the paper gives a theory proof to explain why only the low orders of the coefficients are influenced by polar gaps.Big polar gaps are experimented by a numerical test with inclination of 45°.Finally,considering that the measurement bandwidth can be expanded and accuracy of gradient observations can be improved by superconducting gravity gradiometer(SGG)compared to gradiometer used in Gravity field and steadystate Ocean Circulation Explorer(GOCE),the possibility of detecting time-variable gravity using gravity gradient observations is discussed.The results show that the SGG creates errors in MBW with magnitude of 0.014 m E,which is smaller than the magnitude of the time-variable gravity gradient signals(i.e.,0.02 m E)derived from Gravity Recovery and Climate Experiment(GRACE)gravity field models.This indicates the potential of SGG in time-variable gravity detection.展开更多
In this paper we studied gravitational potential, gravity and its gradient changes, which are caused by a pointdislocation, and gave the concise mathematical deduction with definite physical implication in dealing wit...In this paper we studied gravitational potential, gravity and its gradient changes, which are caused by a pointdislocation, and gave the concise mathematical deduction with definite physical implication in dealing with thesingular integral at a seismic source. We also analysed the features of the fields of gravity and gravity gradient,gravity-vertical-displacement gradient. The conclusions a re: (1) Gravity and gravit y gradient changes are verysmall with the change of vertical positionl (2) Gravity change is much greater than the gravity gradient changewhich is not so distinct; (3) The gravity change due to redistribution of mass accounts for 10-50 percent of thetotal gravity change caused by dislocation. The signs (positive or negative) of total gravity change and verticaldisplacement are opposite each other at the same point for strike slip and dip slip 1 (4) Gravity-vertical displacement-gradient is not constantl it manifests a variety of patterns for different dislocation modelsl (5) Gravityvertical-displacement-gradient is approximately equal to a pparent gra vit y-- vert ical- displacement -gradient.展开更多
When the computational point is approaching the poles, the variance and covariance formulae of the disturbing gravity gradient tensors tend to be infinite, and this is a singular problem. In order to solve the problem...When the computational point is approaching the poles, the variance and covariance formulae of the disturbing gravity gradient tensors tend to be infinite, and this is a singular problem. In order to solve the problem, the authors deduced the practical non-singular computational formulae of the first- and second-order derivatives of the Legendre functions and two kinds of spherical harmonic functions, and then constructed the nonsingular formulae of variance and eovarianee function of disturbing gravity gradient tensors.展开更多
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 calculati...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.展开更多
How to deal with colored noises of GOCE （Gravity field and steady - state Ocean Circulation Explorer） satellite has been the key to data processing. This paper focused on colored noises of GOCE gradient data and the...How to deal with colored noises of GOCE （Gravity field and steady - state Ocean Circulation Explorer） satellite has been the key to data processing. This paper focused on colored noises of GOCE gradient data and the frequency spectrum analysis. According to the analysis results, gravity field model of the optima] degrees 90-240 is given, which is recovered by COCE gradient data. This paper presents an iterative Wiener filtering method based on the gravity gradient invariants. By this method a degree-220 model was calculated from GOCE SGG （Satellite Gravity Gradient） data. The degrees above 90 of ITG2010 were taken as the prior gravity field model, replacing the low degree gravity field model calculated by GOCE orbit data. GOCE gradient colored noises was processed by Wiener filtering. Finally by Wiener filtering iterative calculation, the gravity field model was restored by space-wise harmonic analysis method. The results show that the model＇s accuracy matched well with the ESA＇s （European Space Agency） results by using the same data,展开更多
Gravity Recovery and Climate Experiment（GRACE） observations have been used to de-tect the co-seismic and post-seismic gravity field variations due to the Mw=9.3 Sumatra-Andaman earthquake that occurred on December 2...Gravity Recovery and Climate Experiment（GRACE） observations have been used to de-tect the co-seismic and post-seismic gravity field variations due to the Mw=9.3 Sumatra-Andaman earthquake that occurred on December 26,2004.This article focuses on investigating some gravita-tional effects caused by this huge earthquake.We computed the geoid height changes,the equivalent water height（EWH） changes,and the gravity changes using the GRACE Level-2 monthly spherical harmonic（SH） solutions released by University of Texas Center for Space Research（UTCSR）.The GRACE results agree well with the prediction by a dislocation model and are consistent with the results obtained by some previous scholars.In particular,we calculated the three components of the gravity gradient variations and found that they can recover the seismic-related signature more sensitively due to a certain degree of amplification of the signals.A positive-negative-positive mode predominates in the spatial distribution of the horizontal components of the gravity gradient variations,which is possibly attributed to the anomalies in the crustal density distribution caused by the uplift-subduction effect of the dip-slip earthquake.Moreover,the latitude components of the gravity gradient changes show strong suppression of the north-south stripes,which is due to the along-orbit measurements of the two GRACE satellites.We conclude that the posi-tive-negative-positive mode in latitude gravity gradient changes would be a more sensitive fea-ture to detect the deformations of some major dip-slip earthquakes by GRACE data.展开更多
On the basis of the results of improved analytical expression of computation of gravity anomalies due to a homogeneous polyhedral body composed of polygonal facets, and applying the forward theory with the coordinate ...On the basis of the results of improved analytical expression of computation of gravity anomalies due to a homogeneous polyhedral body composed of polygonal facets, and applying the forward theory with the coordinate transformation of vectors and tensors, we deduced both the analytical expressions for gravity gradient tensors and for magnetic anomalies of a polygon, and obtained new analytical expressions for computing vertical gradients of gravity anomalies and vertical component of magnetic anomalies caused by a polyhedral body. And also we developed explicitly the complete unified expressions for the calculation of gravity anomalies, gravity gradient, and magnetic anomalies due to the homogeneous polyhedron. Furthermore, we deduced new analytical expressions for computing vertical gradients of gravity anomalies due to a finite rectangular prism by applying the newly obtained expressions for gravity gradient tensors due to a polyhedral target body. Comparison with forward calculation of models shows the correctness of these new expressions. It will reduce forward calculation time of gravity-magnetic anomalies and improve computational efficiency by applying our unified expressions for joint forward modeling of gravity-magnetic anomalies due to homogeneous polyhedral bodies.展开更多
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 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.展开更多
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 sme...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.展开更多
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-dimension...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.展开更多
This study adopted the Euler deconvolution method to conduct an inversion and interpretation of the depth and spatial distribution pattern of field source that lead to gravity variation. For this purpose, mobile gravi...This study adopted the Euler deconvolution method to conduct an inversion and interpretation of the depth and spatial distribution pattern of field source that lead to gravity variation. For this purpose, mobile gravity data from four periods in the Hexi region between 2011 and 2015 were obtained from an observation network. With a newly established theoretical model, we acquired the optimum inversion parameters and conducted calculation and analysis with the actual data. The results indicate that one is the appropriate value of the structure index for the inversion of the mobile gravity data. The inversion results of the actual data showed a comparable spatial distribution of the field source and a consistent structural trend with observations from the Qilian-Haiyuan Fault zone between 2011 and 2015. The distribution was in a blocking state at the epicenter of the Menyuan earthquake in 2016. Our quantitative study of the field source provides new insights into the inversion and interpretation of signals of mobile gravity variation.展开更多
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 GOC...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.展开更多
Based on the absolute and relative gravity observations in North China from 2009 to 2014,spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theore...Based on the absolute and relative gravity observations in North China from 2009 to 2014,spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theoretical model to get the best estimates of parameters. Gravity field change caused by the depth and distribution in North China is calculated by back analysis. The results show the structural index that equals 1 is suitable for inversion of the gravity variation data. The inversion results indicate that the depths of anomaly field sources are spread over the Hetao fault. The research method of this paper can be used in the quantitative study on the field source and may shed new light on the interpretations of gravity change, and also provide quantitative basis for earthquake prediction index criterions based on the gravity change.展开更多
The average acceleration approach was applied to recover a gravity field model Model;CA from GOCE precise science orbits from September 2 to November 2, 2010, and furthermore a so called sequential least square adjust...The average acceleration approach was applied to recover a gravity field model Model;CA from GOCE precise science orbits from September 2 to November 2, 2010, and furthermore a so called sequential least square adjustment was used. The model was compared with other gravity field models based on CHAMP, GRACE and GOCE. The result shows that the model is superior to gravity field based on CHAMP, and with higher accuracy than other international gravity field models based on only GOCE data before 80 degree. The degree geoid height of Model;CA reaches 3cm up to 90 degree and order.展开更多
Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to...Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to search over the entire solution space for a more refined result. However, the inversion will be difficult with the increased parameters in the large search space and the number of computations increases exponentially. ｜n this paper, we propose a novel approach based on the frequency characteristics of the density distribution over the mesh. The purposes of our study are to reduce the parameters of three- dimensional gravity inversion and to lighten the image quality of the inversion result. The results show that the new method can expedite the inversion processing and get a better geological interpretation than tradition methods.展开更多
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 locat...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”.展开更多
Asteroid research is of global security interest, we are more afraid about any impact of these asteroids on Earth. We seek through this work to provide a study on the effect of the asteroid’s shape on the gravity asp...Asteroid research is of global security interest, we are more afraid about any impact of these asteroids on Earth. We seek through this work to provide a study on the effect of the asteroid’s shape on the gravity aspect. Knowledge of the potential is crucial to bodies approaching the asteroid. There is a range of asteroid shapes to consider. Some well-known asteroids such Ceres, Vesta, Iris and Oumuamua are considered in this study. After determining the moment of inertia of the asteroids depending on their materials, the gravity fields and the potential gravity of the asteroids are established when varying their shapes. A representation of the gravity field is given in three dimensional coordinate systems. Also, the behavior of the potential gravity is drawing in the function of the object’s location. The second part is dealing with the interpretation of all the obtained results in order to deduce some rules and features which would be useful for the identification of the asteroids. Thus, by the best knowledge of the effect of the asteroid’s shape, we would be more informed in the survey of the hazardous near earth objects.展开更多
Theoretical horizontal displacements caused by the 2004 Sumatra earthquake in the Sichuan-Yunnan area have been calculated according to a spherical dislocation theory and an earthquake-fault model. The results show th...Theoretical horizontal displacements caused by the 2004 Sumatra earthquake in the Sichuan-Yunnan area have been calculated according to a spherical dislocation theory and an earthquake-fault model. The results show that the theoretical displacements are basically consistent with the GPS observations in situ. On this basis,we have calculated the co-seismic displacements, strains, changes of gravity and geoid of the whole Earth, including China mainland and vicinity, caused by this earthquake. Key wards：展开更多
基金funded by National Natural Science Foundation of China(No.41674026,41404019,41774089)Fundamental Research Funds for the Central University(No.2652018027)+2 种基金China Geological Survey(DD20191006)Open Research Fund of Qian Xuesen Laboratory of Space Technology,CAST(No.GZZKFJJ2020006)Open Research Fund of Key Laboratory of Space Utilization,Chinese Academy of Sciences(LSU-KFJJ201902)
文摘Although satellite gravity gradient data plays a great role in determining short-wavelength part of static gravity field model,accuracy of the long-wavelength part of gravity field model recovered by them are poor,which leads to only a few applications in time-variable gravity investigation.The reason is that some factors limit the accuracy of the gravity field recovered using gradient data,including accuracy of the gravity gradient observations,measurement bandwidth(MBW)of gradiometer,satellite inclination,etc.This paper aims at analyzing the influence of these limitations on gravity field recovery and discusses the possibility of time-variable gravity field detection by using gravity gradient observations.Firstly,for arbitrary satellite orbit inclination,we give the frequency distributions of all the components of gravity gradients(i.e.Txx;Tyy;Tzz;Txy;Txz and Tyz,).The results show that the maximum frequency of each component of the gravity gradients is the same,i.e.l=Ts(l is degree of the gravity field model,Ts is the orbital periods),and it is not influenced by the inclination of the satellite orbits.Secondly,the paper gives a theory proof to explain why only the low orders of the coefficients are influenced by polar gaps.Big polar gaps are experimented by a numerical test with inclination of 45°.Finally,considering that the measurement bandwidth can be expanded and accuracy of gradient observations can be improved by superconducting gravity gradiometer(SGG)compared to gradiometer used in Gravity field and steadystate Ocean Circulation Explorer(GOCE),the possibility of detecting time-variable gravity using gravity gradient observations is discussed.The results show that the SGG creates errors in MBW with magnitude of 0.014 m E,which is smaller than the magnitude of the time-variable gravity gradient signals(i.e.,0.02 m E)derived from Gravity Recovery and Climate Experiment(GRACE)gravity field models.This indicates the potential of SGG in time-variable gravity detection.
文摘In this paper we studied gravitational potential, gravity and its gradient changes, which are caused by a pointdislocation, and gave the concise mathematical deduction with definite physical implication in dealing with thesingular integral at a seismic source. We also analysed the features of the fields of gravity and gravity gradient,gravity-vertical-displacement gradient. The conclusions a re: (1) Gravity and gravit y gradient changes are verysmall with the change of vertical positionl (2) Gravity change is much greater than the gravity gradient changewhich is not so distinct; (3) The gravity change due to redistribution of mass accounts for 10-50 percent of thetotal gravity change caused by dislocation. The signs (positive or negative) of total gravity change and verticaldisplacement are opposite each other at the same point for strike slip and dip slip 1 (4) Gravity-vertical displacement-gradient is not constantl it manifests a variety of patterns for different dislocation modelsl (5) Gravityvertical-displacement-gradient is approximately equal to a pparent gra vit y-- vert ical- displacement -gradient.
基金supported by the National 973 Foundation of China(61322201)the National Natural Science Foundation of China(41304022,41174026,41104047)Key Laboratory Foundation of Geo-space Environment and Geodesy,Ministry of Education(11-01-03)
文摘When the computational point is approaching the poles, the variance and covariance formulae of the disturbing gravity gradient tensors tend to be infinite, and this is a singular problem. In order to solve the problem, the authors deduced the practical non-singular computational formulae of the first- and second-order derivatives of the Legendre functions and two kinds of spherical harmonic functions, and then constructed the nonsingular formulae of variance and eovarianee function of disturbing gravity gradient tensors.
基金Joint Seismological Science Foundation of China (No.101005).
文摘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.
基金supported by the National Natural Science Foundation of China(41404020)
文摘How to deal with colored noises of GOCE （Gravity field and steady - state Ocean Circulation Explorer） satellite has been the key to data processing. This paper focused on colored noises of GOCE gradient data and the frequency spectrum analysis. According to the analysis results, gravity field model of the optima] degrees 90-240 is given, which is recovered by COCE gradient data. This paper presents an iterative Wiener filtering method based on the gravity gradient invariants. By this method a degree-220 model was calculated from GOCE SGG （Satellite Gravity Gradient） data. The degrees above 90 of ITG2010 were taken as the prior gravity field model, replacing the low degree gravity field model calculated by GOCE orbit data. GOCE gradient colored noises was processed by Wiener filtering. Finally by Wiener filtering iterative calculation, the gravity field model was restored by space-wise harmonic analysis method. The results show that the model＇s accuracy matched well with the ESA＇s （European Space Agency） results by using the same data,
基金supported by the National Natural Science Foundation of China (Nos. 40974015,40637034)the Fund of Key Laboratory of Geodynamic Geodesy, Chinese Academy of Sciences (No. 09-18)the Fund of Key Laboratory of Geospace Environment and Geodesy,Ministry of Education,China (No. 07-12)
文摘Gravity Recovery and Climate Experiment（GRACE） observations have been used to de-tect the co-seismic and post-seismic gravity field variations due to the Mw=9.3 Sumatra-Andaman earthquake that occurred on December 26,2004.This article focuses on investigating some gravita-tional effects caused by this huge earthquake.We computed the geoid height changes,the equivalent water height（EWH） changes,and the gravity changes using the GRACE Level-2 monthly spherical harmonic（SH） solutions released by University of Texas Center for Space Research（UTCSR）.The GRACE results agree well with the prediction by a dislocation model and are consistent with the results obtained by some previous scholars.In particular,we calculated the three components of the gravity gradient variations and found that they can recover the seismic-related signature more sensitively due to a certain degree of amplification of the signals.A positive-negative-positive mode predominates in the spatial distribution of the horizontal components of the gravity gradient variations,which is possibly attributed to the anomalies in the crustal density distribution caused by the uplift-subduction effect of the dip-slip earthquake.Moreover,the latitude components of the gravity gradient changes show strong suppression of the north-south stripes,which is due to the along-orbit measurements of the two GRACE satellites.We conclude that the posi-tive-negative-positive mode in latitude gravity gradient changes would be a more sensitive fea-ture to detect the deformations of some major dip-slip earthquakes by GRACE data.
基金This paper is supported by the National Natural Science Foundation of China (No.40374039)Program for New Century Excellent Talents in University (No. NCET-04-0726)the Focused Subject Program of Beijing (No. XK104910598).
文摘On the basis of the results of improved analytical expression of computation of gravity anomalies due to a homogeneous polyhedral body composed of polygonal facets, and applying the forward theory with the coordinate transformation of vectors and tensors, we deduced both the analytical expressions for gravity gradient tensors and for magnetic anomalies of a polygon, and obtained new analytical expressions for computing vertical gradients of gravity anomalies and vertical component of magnetic anomalies caused by a polyhedral body. And also we developed explicitly the complete unified expressions for the calculation of gravity anomalies, gravity gradient, and magnetic anomalies due to the homogeneous polyhedron. Furthermore, we deduced new analytical expressions for computing vertical gradients of gravity anomalies due to a finite rectangular prism by applying the newly obtained expressions for gravity gradient tensors due to a polyhedral target body. Comparison with forward calculation of models shows the correctness of these new expressions. It will reduce forward calculation time of gravity-magnetic anomalies and improve computational efficiency by applying our unified expressions for joint forward modeling of gravity-magnetic anomalies due to homogeneous polyhedral bodies.
基金supported by the Fundamental Research Funds for the Central Universities(SWJTU12BR012)
文摘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.
文摘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.
基金Project supported by the Main Direction Program of Knowledge Innovation of the Chinese Academy of Sciences for Distinguished Young Scholars (Grant No. KZCX2-EW-QN114)the National Natural Science Foundation of China for Young Scholars (GrantNos. 41004006,41131067,and 11173049)+3 种基金the Merit-Based Scientific Research Foundation of the State Ministry of Human Resources and Social Security of China for Returned Overseas Chinese Scholars (Grant No. 2011)the Open Research Fund Programof the Key Laboratory of Computational Geodynamics of the Chinese Academy of Sciences (Grant No. 2011-04)the Frontier Field Program of Knowledge Innovation of Institute of Geodesy and Geophysics of the Chinese Academy of Sciencesthe Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Grant No. PLN1113)
文摘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.
基金funded by Science and Technology Project of Shanxi Province (2014K13-04)the National Science Foundation of China (41274083)+1 种基金the Special Fund for Earthquake Scientific Research of China (201308009)the Youth Seismic Regime Tracking Project in the Year of 2016, China Earthquake Administration (2016010222)
文摘This study adopted the Euler deconvolution method to conduct an inversion and interpretation of the depth and spatial distribution pattern of field source that lead to gravity variation. For this purpose, mobile gravity data from four periods in the Hexi region between 2011 and 2015 were obtained from an observation network. With a newly established theoretical model, we acquired the optimum inversion parameters and conducted calculation and analysis with the actual data. The results indicate that one is the appropriate value of the structure index for the inversion of the mobile gravity data. The inversion results of the actual data showed a comparable spatial distribution of the field source and a consistent structural trend with observations from the Qilian-Haiyuan Fault zone between 2011 and 2015. The distribution was in a blocking state at the epicenter of the Menyuan earthquake in 2016. Our quantitative study of the field source provides new insights into the inversion and interpretation of signals of mobile gravity variation.
基金financially supported by the National Key Basic Research Program of China(973 program,grant no.:2013CB733302,2013CB733301)the Major International(Regional) Joint Research Project(grant no.:41210006)+1 种基金DAAD Thematic Network Project(grant no.:57173947)the National Natural Science Foundation of China(grant No.41374022)
文摘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.
基金funded by the Natural Science Foundation of China(61627824,41274083)the Youth Foundation of Earthquake Prediction(2017010227)
文摘Based on the absolute and relative gravity observations in North China from 2009 to 2014,spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theoretical model to get the best estimates of parameters. Gravity field change caused by the depth and distribution in North China is calculated by back analysis. The results show the structural index that equals 1 is suitable for inversion of the gravity variation data. The inversion results indicate that the depths of anomaly field sources are spread over the Hetao fault. The research method of this paper can be used in the quantitative study on the field source and may shed new light on the interpretations of gravity change, and also provide quantitative basis for earthquake prediction index criterions based on the gravity change.
基金supported by the Fundamental Research Funds for the Central Universities(SWJTU12BR012)
文摘The average acceleration approach was applied to recover a gravity field model Model;CA from GOCE precise science orbits from September 2 to November 2, 2010, and furthermore a so called sequential least square adjustment was used. The model was compared with other gravity field models based on CHAMP, GRACE and GOCE. The result shows that the model is superior to gravity field based on CHAMP, and with higher accuracy than other international gravity field models based on only GOCE data before 80 degree. The degree geoid height of Model;CA reaches 3cm up to 90 degree and order.
基金supported by the Key Project Fund of the Chinese Academy of Sciences under grant number (kzcx2-yw-203-01)the Major State Basic Research Development Program of China(973 Program,Grant No.2007CB41170404)
文摘Three-dimensional gravity inversion based on the mass property model is very popular in recent years. The time and efficiency of inversion algorithms is relative to the magnitude of the target mesh. One approach is to search over the entire solution space for a more refined result. However, the inversion will be difficult with the increased parameters in the large search space and the number of computations increases exponentially. ｜n this paper, we propose a novel approach based on the frequency characteristics of the density distribution over the mesh. The purposes of our study are to reduce the parameters of three- dimensional gravity inversion and to lighten the image quality of the inversion result. The results show that the new method can expedite the inversion processing and get a better geological interpretation than tradition methods.
基金National Natural Science Foundation of China(41674107,41874119,41574064)National Key Research and Development Program of China(2017YFB0202904).
文摘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”.
文摘Asteroid research is of global security interest, we are more afraid about any impact of these asteroids on Earth. We seek through this work to provide a study on the effect of the asteroid’s shape on the gravity aspect. Knowledge of the potential is crucial to bodies approaching the asteroid. There is a range of asteroid shapes to consider. Some well-known asteroids such Ceres, Vesta, Iris and Oumuamua are considered in this study. After determining the moment of inertia of the asteroids depending on their materials, the gravity fields and the potential gravity of the asteroids are established when varying their shapes. A representation of the gravity field is given in three dimensional coordinate systems. Also, the behavior of the potential gravity is drawing in the function of the object’s location. The second part is dealing with the interpretation of all the obtained results in order to deduce some rules and features which would be useful for the identification of the asteroids. Thus, by the best knowledge of the effect of the asteroid’s shape, we would be more informed in the survey of the hazardous near earth objects.
基金supported by Basic Research Foundation from Institute of Earthquake Science,CEA(0210240101)
文摘Theoretical horizontal displacements caused by the 2004 Sumatra earthquake in the Sichuan-Yunnan area have been calculated according to a spherical dislocation theory and an earthquake-fault model. The results show that the theoretical displacements are basically consistent with the GPS observations in situ. On this basis,we have calculated the co-seismic displacements, strains, changes of gravity and geoid of the whole Earth, including China mainland and vicinity, caused by this earthquake. Key wards：
