scholarly journals On inverting gravity changes with the harmonic inversion method: Teide (Tenerife) case study

2015 ◽  
Vol 45 (2) ◽  
pp. 111-134
Author(s):  
Vladimír Pohánka ◽  
Peter Vajda ◽  
Jaroslava Pánisová
Keyword(s):  
Gravity Field ◽  
Integral Transformation ◽  
Gravity Data ◽  
Mass Density ◽  
Surface Gravity ◽  
Inversion Method ◽  
Volcanic Complex ◽  
Size And Shape ◽  
Gravity Changes ◽  
Harmonic Inversion

Abstract Here we investigate the applicability of the harmonic inversion method to time-lapse gravity changes observed in volcanic areas. We carry out our study on gravity changes occured over the period of 2004–2005 during the unrest of the Central Volcanic Complex on Tenerife, Canary Islands. The harmonic inversion method is unique in that it calculates the solution of the form of compact homogeneous source bodies via the mediating 3-harmonic function called quasigravitation. The latter is defined in the whole subsurface domain and it is a linear integral transformation of the surface gravity field. At the beginning the seeds of the future source bodies are introduced: these are quasi-spherical bodies located at the extrema of the quasigravitation (calculated from the input gravity data) and their differential densities are free parameters preselected by the interpreter. In the following automatic iterative process the source bodies change their size and shape according to the local values of quasigravitation (calculated in each iterative step from the residual surface gravity field); the process stops when the residual surface gravity field is sufficiently small. In the case of inverting temporal gravity changes, the source bodies represent the volumetric domains of temporal mass-density changes. The focus of the presented work is to investigate the dependence of the size and shape of the found source bodies on their differential densities. We do not aim here (yet) at interpreting the found solutions in terms of volcanic processes associated with intruding or rejuvenating magma and/or migrating volatiles.

Antarctica crustal model by means of the Bayesian gravity inversion

2020 ◽  
Author(s):  
Martina Capponi ◽  
Daniele Sampietro
Keyword(s):  
Gravity Field ◽  
Gravity Data ◽  
Mass Density ◽  
Geophysical Methods ◽  
Gravity Inversion ◽  
Inversion Algorithm ◽  
Global Gravity ◽  
Uniqueness Of The Solution ◽  
Definition Of ◽  
Regional Gravity

<p>The Antarctica crustal structure is still not completely unveiled due to the presence of thick ice shields all over the continent which prevent direct in situ measurements. In the last decades, various geophysical methods have been used to retrieve information of the upper crust and sediments distribution however there are still regions, especially in central Antarctica, where our knowledge is limited. For these kind of situations, in which the indirect investigation of the subsurface is the only valuable solution, the gravity data are an important source of information. After the recent dedicated satellite missions, like GRACE and GOCE, it has been possible to obtain global gravity field data with spatial resolution and accuracy almost comparable to those of local/regional gravity acquisitions, paving the way to new geophysical applications. The new challenge today is the capability to invert such gravity data on large areas with the aim to obtain a 3D density model of the Earth crust. This is in fact a problem characterized by intrinsic instability and non-uniqueness of the solution that to be solved requires the definition of strong constrains and numerical regularization.</p><p>In this work the authors propose the application of a Bayesian inversion algorithm to the Antarctica continent to infer a model of mass density distribution. The first operation is the definition of an initial geological model in terms of geological horizons and density. These two variables are considered as random variables and, within the iterative procedure based on Markov Chain Monte Carlo methods, they are adjusted in such a way to fit the gravity field on the surface. The test performed show that the method is capable of retrieving an estimated model consistent with the prior information and fitting the gravity observations according to their accuracy.</p>

Highly efficient density inversion of gravity data using nonlinear density polynomial fitting

Geophysics ◽  
2021 ◽  
pp. 1-34
Author(s):  
Guoqing Ma ◽  
Zongrui Li ◽  
Lili Li ◽  
Taihan Wang
Keyword(s):  
Weighting Function ◽  
Gravity Data ◽  
Real Life ◽  
Density Variation ◽  
Density Change ◽  
Inversion Method ◽  
Density Inversion ◽  
Perspective View ◽  
Polynomial Fitting ◽  
Polynomial Coefficients

The density inversion of gravity data is commonly achieved by discretizing the subsurface into prismatic cells and calculating the density of each cell. During this process, a weighting function is introduced to the iterative computation to reduce the skin effect during the inversion. Thus, the computation process requires a significant number of matrix operations, which results in low computational efficiency. We have adopted a density inversion method with nonlinear polynomial fitting (NPF) that uses a polynomial to represent the density variation of prismatic cells in a certain space. The computation of each cell is substituted by the computation of the nonlinear polynomial coefficients. Consequently, the efficiency of the inversion is significantly improved because the number of nonlinear polynomial coefficients is less than the number of cells used. Moreover, because representing the density change of all of the cells poses a significant challenge when the cell number is large, we adopt the use of a polynomial to represent the density change of a subregion with fewer cells and multiple nonlinear polynomials to represent the density changes of all prism cells. Using theoretical model tests, we determine that the NPF method more efficiently recovers the density distribution of gravity data compared with conventional density inversion methods. In addition, the density variation of a subregion with 8 × 8 × 8 prismatic cells can be accurately and efficiently obtained using our cubic NPF method, which can also be used for noisy data. Finally, the NPF method was applied to real gravity data in an iron mining area in Shandong Province, China. Convergent results of a 3D perspective view and the distribution of the iron ore bodies were acquired using this method, demonstrating the real-life applicability of this method.

scholarly journals New insights on the dynamics of the Sumatra and Mariana complexes inferred from the comparative analysis of gravity data and model predictions

2020 ◽  
Author(s):  
Arcangela Bollino ◽  
Anna Maria Marotta ◽  
Federica Restelli ◽  
Alessandro Regorda ◽  
Roberto Sabadini
Keyword(s):  
Comparative Analysis ◽  
Gravity Field ◽  
Wide Spectrum ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Phase Changes ◽  
Wide Range ◽  
Model Predictions ◽  
Dip Angle ◽  
The Masses

<p>Subduction is responsible for surface displacements and deep mass redistribution. This rearrangement generates density anomalies in a wide spectrum of wavelengths which, in turn, causes important anomalies in the Earth's gravity field that are visible as lineaments parallel to the arc-trench systems. In these areas, when the traditional analysis of the deformation and stress fields is combined with the analysis of the perturbation of the gravity field and its slow time variation, new information on the background environment controlling the tectonic loading phase can be disclosed.</p><p>Here we present the results of a comparative analysis between the geodetically retrieved gravitational anomalies, based on the EIGEN-6C4 model, and those predicted by a 2D thermo-chemical mechanical modeling of the Sumatra and Mariana complexes.</p><p>The 2D model accounts for a wide range of parameters, such as the convergence velocity, the shallow dip angle, the different degrees of coupling between the facing plates. The marker in cell technique is used to compositionally differentiate the system. Phase changes in the crust and in the mantle and mantle hydration are also allowed. To be compliant with the geodetic EIGEN-6C4 gravity data, we define a model normal Earth considering the vertical density distribution at the margins of the model domain, where the masses are not perturbed by the subduction process.</p><p>Model predictions are in good agreement with data, both in terms of wavelengths and magnitude of the gravity anomalies measured in the surroundings of the Sumatra and Marina subductions. Furthermore, our modeling supports that the differences in the style of the gravity anomaly observed in the two areas are attributable to the different environments – ocean-ocean or ocean-continental subduction – that drives a significantly different dynamic in the wedge area.</p>

scholarly journals Accuracy Analysis of Gravity Field Changes from Grace RL06 and RL05 Data Compared to in Situ Gravimetric Measurements in the Context of Choosing Optimal Filtering Type

2020 ◽  
Vol 55 (3) ◽  
pp. 100-117
Author(s):  
Viktor Szabó ◽  
Dorota Marjańska
Keyword(s):  
Gravity Field ◽  
Ocean Circulation ◽  
Gravity Data ◽  
The Other ◽  
Subsurface Water ◽  
Satellite Gravity ◽  
Absolute Gravimetry ◽  
The Earth ◽  
Gravity Measurements ◽  
Gravity Recovery

AbstractGlobal satellite gravity measurements provide unique information regarding gravity field distribution and its variability on the Earth. The main cause of gravity changes is the mass transportation within the Earth, appearing as, e.g. dynamic fluctuations in hydrology, glaciology, oceanology, meteorology and the lithosphere. This phenomenon has become more comprehensible thanks to the dedicated gravimetric missions such as Gravity Recovery and Climate Experiment (GRACE), Challenging Minisatellite Payload (CHAMP) and Gravity Field and Steady-State Ocean Circulation Explorer (GOCE). From among these missions, GRACE seems to be the most dominating source of gravity data, sharing a unique set of observations from over 15 years. The results of this experiment are often of interest to geodesists and geophysicists due to its high compatibility with the other methods of gravity measurements, especially absolute gravimetry. Direct validation of gravity field solutions is crucial as it can provide conclusions concerning forecasts of subsurface water changes. The aim of this work is to present the issue of selection of filtration parameters for monthly gravity field solutions in RL06 and RL05 releases and then to compare them to a time series of absolute gravimetric data conducted in quasi-monthly measurements in Astro-Geodetic Observatory in Józefosław (Poland). The other purpose of this study is to estimate the accuracy of GRACE temporal solutions in comparison with absolute terrestrial gravimetry data and making an attempt to indicate the significance of differences between solutions using various types of filtration (DDK, Gaussian) from selected research centres.

Towards an updated, enhanced regional gravity field solution for Antarctica

2021 ◽  
Author(s):  
Mirko Scheinert ◽  
Philipp Zingerle ◽  
Theresa Schaller ◽  
Roland Pail ◽  
Martin Willberg
Keyword(s):  
Gravity Anomaly ◽  
Gravity Field ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Data Set ◽  
Gravity Field Model ◽  
Terrestrial Gravity ◽  
Field Solution ◽  
Gravity Disturbances ◽  
Regional Gravity

<p>In the frame of the IAG Subcommission 2.4f “Gravity and Geoid in Antarctica” (AntGG) a first Antarctic-wide grid of ground-based gravity anomalies was released in 2016 (Scheinert et al. 2016). That data set was provided with a grid space of 10 km and covered about 73% of the Antarctic continent. Since then a considerably amount of new data has been made available, mainly collected by means of airborne gravimetry. Regions which were formerly void of any terrestrial gravity observations and have now been surveyed include especially the polar data gap originating from GOCE satellite gravimetry. Thus, it is timely to come up with an updated and enhanced regional gravity field solution for Antarctica. For this, we aim to improve further aspects in comparison to the AntGG 2016 solution: The grid spacing will be enhanced to 5 km. Instead of providing gravity anomalies only for parts of Antarctica, now the entire continent should be covered. In addition to the gravity anomaly also a regional geoid solution should be provided along with further desirable functionals (e.g. gravity anomaly vs. disturbance, different height levels).</p><p>We will discuss the expanded AntGG data base which now includes terrestrial gravity data from Antarctic surveys conducted over the past 40 years. The methodology applied in the analysis is based on the remove-compute-restore technique. Here we utilize the newly developed combined spherical-harmonic gravity field model SATOP1 (Zingerle et al. 2019) which is based on the global satellite-only model GOCO05s and the high-resolution topographic model EARTH2014. We will demonstrate the feasibility to adequately reduce the original gravity data and, thus, to also cross-validate and evaluate the accuracy of the data especially where different data set overlap. For the compute step the recently developed partition-enhanced least-squares collocation (PE-LSC) has been used (Zingerle et al. 2021, in review; cf. the contribution of Zingerle et al. in the same session). This method allows to treat all data available in Antarctica in one single computation step in an efficient and fast way. Thus, it becomes feasible to iterate the computations within short time once any input data or parameters are changed, and to easily predict the desirable functionals also in regions void of terrestrial measurements as well as at any height level (e.g. gravity anomalies at the surface or gravity disturbances at constant height).</p><p>We will discuss the results and give an outlook on the data products which shall be finally provided to present the new regional gravity field solution for Antarctica. Furthermore, implications for further applications will be discussed e.g. with respect to geophysical modelling of the Earth’s interior (cf. the contribution of Schaller et al. in session G4.3).</p>

scholarly journals Performance Comparison of Geoid Refinement between XGM2016 and EGM2008 Based on the KTH and RCR Methods: Jilin Province, China

2020 ◽  
Vol 12 (2) ◽  
pp. 324
Author(s):  
Qiong Wu ◽  
Hongyao Wang ◽  
Bin Wang ◽  
Shengbo Chen ◽  
Hongqing Li
Keyword(s):  
Gravity Field ◽  
Mean Square Error ◽  
Gravity Data ◽  
Jilin Province ◽  
Mountainous Area ◽  
Mean Square ◽  
Geoid Model ◽  
Plain Area ◽  
Gravity Field Models ◽  
Global Mean

The selection of an appropriate global gravity field model and refinement method can effectively improve the accuracy of the refined regional geoid model in a certain research area. We analyzed the accuracy of Experimental Geopotential Model (XGM2016) based on the GPS-leveling data and the modification parameters of the global mean square errors in the KTH geoid refinement in Jilin Province, China. The regional geoid was refined based on Earth Gravitational Model (EGM2008) and XGM2016 using both the Helmert condensation method with an RCR procedure and the KTH method. A comparison of the original two gravity field models to the GPS-leveling benchmarks showed that the accuracies of XGM2016 and EGM2008 in the plain area of Jilin Province are similar with a standard deviation (STD) of 5.8 cm, whereas the accuracy of EGM2008 in the high mountainous area is 1.4 cm better than that of XGM2016, which may be attributed to its low resolution. The modification parameters between the two gravity field models showed that the coefficient error of XGM2016 model was lower than that of EGM2008 for the same degree of expansion. The different modification limits and integral radii may produce a centimeter level difference in global mean square error, while the influence of the truncation error caused by the integral was at the millimeter level. The terrestrial gravity data error accounted for the majority of the global mean square error. The optimal least squares modification obtained the minimum global mean square error, and the global mean square error calculated based on XGM2016 model was reduced by about 1~3 cm compared with EGM2008. The refined geoid based on the two gravity field models indicated that both KTH and RCR method can effectively improve the STD of the geoid model from about six to three centimeters. The refined accuracy of EGM2008 model using RCR and KTH methods is slightly better than that of XGM2016 model in the plain and high mountain areas after seven-parameter fitting. EGM2008 based on the KTH method was the most precise at ± 2.0 cm in the plain area and ± 2.4 cm in the mountainous area. Generally, for the refined geoid based on the two Earth gravity models, KTH produced results similar to RCR in the plain area, and had relatively better performance for the mountainous area where terrestrial gravity data is sparse and unevenly distributed.

scholarly journals Evaluation of Gravity Data Derived from Global Gravity Field Models Using Terrestrial Gravity Data in Enugu State, Nigeria

2018 ◽  
Vol 8 (1) ◽  
pp. 145-153 ◽  
Author(s):  
O.I. Apeh ◽  
E.C. Moka ◽  
V.N. Uzodinma
Keyword(s):  
Gravity Field ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Mean Square ◽  
Bouguer Gravity ◽  
Bouguer Anomalies ◽  
The Earth ◽  
Global Gravity ◽  
Enugu State ◽  
Gravity Field Models

Abstract Spherical harmonic expansion is a commonly applied mathematical representation of the earth’s gravity field. This representation is implied by the potential coeffcients determined by using elements/parameters of the field observed on the surface of the earth and/or in space outside the earth in the spherical harmonic expansion of the field. International Centre for Gravity Earth Models (ICGEM) publishes, from time to time, Global Gravity Field Models (GGMs) that have been developed. These GGMs need evaluation with terrestrial data of different locations to ascertain their accuracy for application in those locations. In this study, Bouguer gravity anomalies derived from a total of eleven (11) recent GGMs, using sixty sample points, were evaluated by means of Root-Mean-Square difference and correlation coeficient. The Root-Mean-Square differences of the computed Bouguer anomalies from ICGEMwebsite compared to their positionally corresponding terrestrial Bouguer anomalies range from 9.530mgal to 37.113mgal. Additionally, the correlation coe_cients of the structure of the signal of the terrestrial and GGM-derived Bouguer anomalies range from 0.480 to 0.879. It was observed that GECO derived Bouguer gravity anomalies have the best signal structure relationship with the terrestrial data than the other ten GGMs. We also discovered that EIGEN-6C4 and GECO derived Bouguer anomalies have enormous potential to be used as supplements to the terrestrial Bouguer anomalies for Enugu State, Nigeria.

Sign in / Sign up

Export Citation Format

Share Document