Exploring viable geologic interpretations of gravity models using distance-based global sensitivity analysis and kernel methods

Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. G79-G92 ◽  
Author(s):  
Geoff Phelps ◽  
Celine Scheidt ◽  
Jef Caers
Keyword(s):  
Sensitivity Analysis ◽  
Gravity Anomaly ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Gravity Models ◽  
Normal Faults ◽  
Reverse Fault ◽  
Density Values ◽  
Generalized Sensitivity

ABSTRACT We have explored ways to integrate alternative geologic interpretations into the modeling of gravity data. These methods are applied to the Vaca Fault east of Fairfield, California, USA, where the structure across the fault is in question, and the Vaca Fault is used as a case study to demonstrate the method. The Vaca Fault is modeled using gravity data collected along a 10 km line perpendicular to the strike of the fault. Of particular interest is how the gravity data might inform on the dip of the Vaca Fault and thickness of the nonmarine section and whether spatial autocorrelation of density internal to the geologic units significantly influences the resulting gravity anomaly. We approach these questions by creating a suite of structural geologic models, which we then populate with geostatistically generated densities and from which the respective synthetic gravity anomalies are calculated. We perform distance-based generalized sensitivity analysis to identify which model inputs most leverage the calculated gravity anomaly. We then use multidimensional scaling to transform the gravity anomalies into a metric space and estimate the posterior probabilities of each structural geologic model using a Bayesian approach. We find that the gravity anomalies are particularly sensitive to zones of autocorrelated density values generated from geostatistical modeling. The structural geologic models most likely to produce gravity anomalies that match the observed data are the moderately dipping normal faults, 45° and 60°, although the probability that the fault dips more steeply, including in a strike slip or reverse fault orientation, is approximately 30%. The probability of a thicker nonmarine unit is 67%, more probable than a thinner nonmarine unit. This suggests that the Vaca Fault dips moderately to the east and truncates a thicker nonmarine unit, but that any further process modeling should include alternatives of the geologic structures.

Geophysical constraints on the nature of the Highland Boundary Fault Zone in western Scotland

1992 ◽  
Vol 129 (4) ◽  
pp. 411-419 ◽  
Author(s):  
M. C. Dentith ◽  
A. Trench ◽  
B. J. Bluck
Keyword(s):  
Fault Zone ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Thrust Fault ◽  
Gravity Models ◽  
Reverse Fault ◽  
Red Sandstone ◽  
Boundary Fault ◽  
Models Of Gravity ◽  
Surface Geology

AbstractPreviously published models of gravity anomalies across the Highland Boundary Fault in western Scotland interpret this structure as a high-angle reverse fault. These gravity anomalies have been re-interpreted in the light of more extensive gravity data now available, and new density data from the Highland Border Complex. The new data suggest that earlier interpretations have overestimated the fault anomaly and used over-simplified density models. New gravity models of the Highland Boundary Fault Zone are presented which show that the interface between the Dalradian and Highland Border Complex dips to the northwest at an angle of about 20°. We interpret the contact between these two formations as a thrust fault. The interface between the Highland Border Complex and the Lower Old Red Sandstone is shown to be vertical as suggested by surface geology, with the latter rocks a few hundred metres thick.

scholarly journals Marquardt inverse modeling of the residual gravity anomalies due to simple geometric structures: A case study of chromite deposit

2019 ◽  
Vol 49 (2) ◽  
pp. 153-180 ◽  
Author(s):  
Ata Eshaghzadeh ◽  
Alireza Dehghanpour ◽  
Sanaz Seyedi Sahebari
Keyword(s):  
Gravity Anomaly ◽  
Gravity Data ◽  
Synthetic Data ◽  
Vertical Cylinder ◽  
Gravity Anomalies ◽  
Horizontal Cylinder ◽  
Inversion Method ◽  
Chromite Deposit ◽  
Different Shapes

Abstract In this paper, an inversion method based on the Marquardt’s algorithm is presented to invert the gravity anomaly of the simple geometric shapes. The inversion outputs are the depth and radius parameters. We investigate three different shapes, i.e. the sphere, infinite horizontal cylinder and semi-infinite vertical cylinder for modeling. The proposed method is used for analyzing the gravity anomalies from assumed models with different initial parameters in all cases as the synthetic data are without noise and also corrupted with noise to evaluate the ability of the procedure. We also employ this approach for modeling the gravity anomaly due to a chromite deposit mass, situated east of Sabzevar, Iran. The lowest error between the theoretical anomaly and computed anomaly from inverted parameters, determine the shape of the causative mass. The inversion using different initial models for the theoretical gravity and also for real gravity data yields approximately consistent solutions. According to the interpreted parameters, the best shape that can imagine for the gravity anomaly source is the vertical cylinder with a depth to top of 7.4 m and a radius of 11.7 m.

scholarly journals Basin Depth Mapping Beneath Wellington City Based on Residual Gravity Anomalies

2021 ◽  
Author(s):  
◽  
Alistair Stronach
Keyword(s):  
Gravity Anomaly ◽  
Sedimentary Basin ◽  
Gravity Data ◽  
Three Dimensional ◽  
Maximum Amplitude ◽  
Central Business District ◽  
Gravity Anomalies ◽  
Depth Map ◽  
Capital City ◽  
Detailed Knowledge

<p><b>New Zealand’s capital city of Wellington lies in an area of high seismic risk, which is further increased by the sedimentary basin beneath the Central Business District (CBD). Ground motion data and damage patterns from the 2013 Cook Strait and 2016 Kaikōura earthquakes indicate that two- and three-dimensional amplification effects due to the Wellington sedimentary basin may be significant. These effects are not currently accounted for in the New Zealand Building Code. In order for this to be done, three-dimensional simulations of earthquake shaking need to be undertaken, which requires detailed knowledge of basin geometry. This is currently lacking, primarily because of a dearth of deep boreholes in the CBD area, particularly in Thorndon and Pipitea where sediment depths are estimated to be greatest.</b></p> <p>A new basin depth map for the Wellington CBD has been created by conducting a gravity survey using a modern Scintrex CG-6 gravity meter. Across the study area, 519 new high precision gravity measurements were made and a residual anomaly map created, showing a maximum amplitude anomaly of -6.2 mGal with uncertainties better than ±0.1 mGal. Thirteen two-dimensional geological profiles were modelled to fit the anomalies, then combined with existing borehole constraints to construct the basin depth map. </p> <p>Results indicate on average greater depths than in existing models, particularly in Pipitea where depths are interpreted to be as great as 450 m, a difference of 250 m. Within 1 km of shore depths are interpreted to increase further, to 600 m. The recently discovered basin bounding Aotea Fault is resolved in the gravity data, where the basement is offset by up to 13 m, gravity anomaly gradients up to 8 mGal/km are observed, and possible multiple fault strands identified. A secondary strand of the Wellington Fault is also identified in the north of Pipitea, where gravity anomaly gradients up to 18 mGal/km are observed.</p>

The gravity field over the Ungava Bay region from satellite altimetry and new land-based data: implications for the geology of the area

1999 ◽  
Vol 36 (1) ◽  
pp. 75-89 ◽  
Author(s):  
Hamid Telmat ◽  
Jean-Claude Mareschal ◽  
Clément Gariépy
Keyword(s):  
Satellite Altimetry ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Gravity Models ◽  
Crustal Thickening ◽  
Seismic Line ◽  
Crustal Thinning ◽  
Gravity Measurements ◽  
George River ◽  
Superior Craton

Gravity data were obtained along two transects on the southern coast of Ungava Bay, which provide continuous gravity coverage between Leaf Bay and George River. The transects and the derived gravity profiles extend from the Superior craton to the Rae Province across the New Quebec Orogen (NQO). Interpretation of the transect along the southwestern coast of Ungava Bay suggests crustal thickening beneath the NQO and crustal thinning beneath the Kuujjuaq Terrane, east of the NQO. Two alternative interpretations are proposed for the transect along the southeastern coast of the bay. The first model shows crustal thickening beneath the George River Shear Zone (GRSZ) and two shallow bodies correlated with the northern extensions of the GRSZ and the De Pas batholith. The second model shows constant crustal thickness and bodies more deeply rooted than in the first model. The gravity models are consistent with the easterly dipping reflections imaged along a Lithoprobe seismic line crossing Ungava Bay and suggest westward thrusting of the Rae Province over the NQO. Because no gravity data have been collected in Ungava Bay, satellite altimetry data have been used as a means to fill the gap in data collected at sea. The satellite-derived gravity data and standard Bouguer gravity data were combined in a composite map for the Ungava Bay region. The new land-based gravity measurements were used to verify and calibrate the satellite data and to ensure that offshore gravity anomalies merge with those determined by the land surveys in a reasonable fashion. Three parallel east-west gravity profiles were extracted: across Ungava Bay (59.9°N), on the southern shore of the bay (58.5°N), and onshore ~200 km south of Ungava Bay (57.1°N). The gravity signature of some major structures, such as the GRSZ, can be identified on each profile.

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

&lt;p&gt;In the frame of the IAG Subcommission 2.4f &amp;#8220;Gravity and Geoid in Antarctica&amp;#8221; (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).&lt;/p&gt;&lt;p&gt;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).&lt;/p&gt;&lt;p&gt;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&amp;#8217;s interior (cf. the contribution of Schaller et al. in session G4.3).&lt;/p&gt;

Three-dimensional gravity modelling applied to the exploration of uranium unconformity-related basement-hosted deposits: the Contact prospect case study, Kiggavik, northeast Thelon region (Nunavut, Canada)

2017 ◽  
Vol 54 (8) ◽  
pp. 869-882 ◽  
Author(s):  
Régis Roy ◽  
Antonio Benedicto ◽  
Alexis Grare ◽  
Mickaël Béhaegel ◽  
Yoann Richard ◽  
...  
Keyword(s):  
Gravity Data ◽  
Three Dimensional ◽  
Gravity Anomalies ◽  
3D Models ◽  
Low Density ◽  
Uranium Deposits ◽  
Geological Interpretation ◽  
Thelon Basin ◽  
Dimensional Gravity ◽  
Density Values

In unconformity-related uranium deposits, mineralization is associated with hydrothermal clay-rich alteration haloes that decrease the density of the host rock. In the Kiggavik uranium project, located in the eastern Thelon Basin, Nunavut (Canada), basement-hosted shallow deposits were discovered by drilling geophysical anomalies in the 1970s. In 2014, gravity data were inverted for the first time using the Geosoft VOXI Earth ModellingTM system to generate three-dimensional (3D) models to assist exploration in the Contact prospect, the most recent discovery at Kiggavik. A 3D unconstrained inversion model was calculated before drilling, and a model constrained by petrophysical data was computed after drilling. The unconstrained inversion provided a first approximation of the geometry and depth of a low-density body and helped to collar the discovery holes of the Contact mineralization. The constrained inversion was computed using density values measured on 315 core samples collected from 21 drill holes completed between 2014 and 2015. The constrained modelling highlights three shallower and smaller low-density bodies that match the geological interpretation and refines the footprint of the gravity anomalies in relation to the current understanding of the deposit. The 3D inversion of gravity data is a valuable tool to guide geologists in exploration of shallow basement-hosted uranium deposits associated with alteration haloes and to assess the deposit gravity geometry.

scholarly journals GRAVITY ANOMALY ASSESSMENT USING GGMS AND AIRBORNE GRAVITY DATA TOWARDS BATHYMETRY ESTIMATION

2016 ◽  
Vol XLII-4/W1 ◽  
pp. 287-297
Author(s):  
A. Tugi ◽  
A. H. M. Din ◽  
K. M. Omar ◽  
A. S. Mardi ◽  
Z. A. M. Som ◽  
...  
Keyword(s):  
Gravity Anomaly ◽  
Gravity Field ◽  
Ocean Circulation ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Airborne Gravity ◽  
Satellite Gravity ◽  
Earth’S Gravity Field ◽  
Explorer Mission ◽  
Satellite Gravity Missions

The Earth’s potential information is important for exploration of the Earth’s gravity field. The techniques of measuring the Earth’s gravity using the terrestrial and ship borne technique are time consuming and have limitation on the vast area. With the space-based measuring technique, these limitations can be overcome. The satellite gravity missions such as Challenging Mini-satellite Payload (CHAMP), Gravity Recovery and Climate Experiment (GRACE), and Gravity-Field and Steady-State Ocean Circulation Explorer Mission (GOCE) has introduced a better way in providing the information on the Earth’s gravity field. From these satellite gravity missions, the Global Geopotential Models (GGMs) has been produced from the spherical harmonics coefficient data type. The information of the gravity anomaly can be used to predict the bathymetry because the gravity anomaly and bathymetry have relationships between each other. There are many GGMs that have been published and each of the models gives a different value of the Earth’s gravity field information. Therefore, this study is conducted to assess the most reliable GGM for the Malaysian Seas. This study covered the area of the marine area on the South China Sea at Sabah extent. Seven GGMs have been selected from the three satellite gravity missions. The gravity anomalies derived from the GGMs are compared with the airborne gravity anomaly, in order to figure out the correlation (R<sup>2</sup>) and the root mean square error (RMSE) of the data. From these assessments, the most suitable GGMs for the study area is GOCE model, GO_CONS_GCF_2_TIMR4 with the R<sup>2</sup> and RMSE value of 0.7899 and 9.886 mGal, respectively. This selected model will be used in the estimating the bathymetry for Malaysian Seas in future.

scholarly journals Regional Recovery of Gravity Anomaly from the Inversion of Diagonal Components of GOCE Gravitational Tensor: A Case Study in Ethiopia

2018 ◽  
Vol 53 (2) ◽  
pp. 55-74 ◽  
Author(s):  
Mehdi Eshagh ◽  
Andenet A. Gedamu ◽  
Tulu B. Bedada
Keyword(s):  
Gravity Anomaly ◽  
Gravity Field ◽  
Ocean Circulation ◽  
Opposite Sign ◽  
Joint Inversion ◽  
Gravity Anomalies ◽  
The Earth ◽  
Tikhonov Regularisation ◽  
The Difference

Abstract The tensor of gravitation is traceless as the gravitational field of the Earth is harmonic outside the Earth’s surface. Therefore, summation of the 2nd-order horizontal derivatives on its diagonal components should be equal to the radial one but with the opposite sign. The gravity field can be recovered locally from either of them, or even their combination. Here, we use the in-orbit diagonal components of the gravitational tensor measured by the gravity field and steady state ocean circulation explorer (GOCE) mission for recovering gravity anomaly with a resolution of 1°×1° at sea level in Ethiopia. In order to solve the system of equations, derived after discretisation of integral equations, the Tikhonov regularisation is applied and the bias of this regularisation is estimated and removed from the estimated gravity anomalies. The errors of the anomalies are estimated and their significance of recovery from these diagonal components is investigated. Statistically, the difference between the recovered anomalies from each scenario is not significant comparing to their errors. However, their joint inversion of the diagonal components improved the solution by about 1 mGal. Furthermore, the inversion processes are better stabilised when using errors of the input data compared with its exclusion, but at the penalty of degradation in accuracy of the estimates.

Bathymetry Inversion Using Vertical Deflections: a Case Study in South China Sea

2021 ◽  
Author(s):  
Xiaoyun Wan ◽  
Bo Liu ◽  
Xiaohong Sui ◽  
Richar Fiifi Annan ◽  
Yijun Min
Keyword(s):  
South China Sea ◽  
Gravity Anomaly ◽  
South China ◽  
Additional Data ◽  
Gravity Anomalies ◽  
Local Area ◽  
China Sea ◽  
The North ◽  
East Component

Abstract As an alternative method, an algorithm for bathymetry inversion using vertical deflections is proposed. Firstly, the formulas for the bathymetry inversion from north and east components of vertical deflections are derived and the data processing is introduced. Then a local area in the South China Sea is selected as an example to experiment the method. The bathymetry inversion based on gravity anomaly was also conducted for a comparison. The results show that the bathymetry derived from the north component of the vertical deflections have almost the same accuracy as that derived from gravity anomalies and the results derived from the east component have the poorest accuracy. The experiment’s results also show that accuracy of the derived bathymetry can be improved if the fitting parameters are adjusted according to the water depths. In summary, among the gravity field products used in this study, although the gravity anomaly yielded the best performance in the bathymetry inversion, the vertical defections can still be used as supplements, especially in areas where accurate vertical deflections exist. This is because deriving gravity anomaly from altimetry observations needs additional data and calculation efforts.

Forward modeling of gravity data using geostatistically generated subsurface density variations

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. G81-G94 ◽  
Author(s):  
Geoff Phelps
Keyword(s):  
Gravity Data ◽  
Solution Space ◽  
Cumulative Effect ◽  
Normal Fault ◽  
Gravity Anomalies ◽  
Reverse Fault ◽  
Heterogeneous Distribution ◽  
Geostatistical Models ◽  
Air Force Base ◽  
Models Of Gravity

Using geostatistical models of density variations in the subsurface, constrained by geologic data, forward models of gravity anomalies can be generated by discretizing the subsurface and calculating the cumulative effect of each cell (pixel). The results of such stochastically generated forward gravity anomalies can be compared with the observed gravity anomalies to find density models that match the observed data. These models have an advantage over forward gravity anomalies generated using polygonal bodies of homogeneous density because generating numerous realizations explores a larger region of the solution space. The stochastic modeling can be thought of as dividing the forward model into two components: that due to the shape of each geologic unit and that due to the heterogeneous distribution of density within each geologic unit. The modeling demonstrates that the internally heterogeneous distribution of density within each geologic unit can contribute significantly to the resulting calculated forward gravity anomaly. Furthermore, the stochastic models match observed statistical properties of geologic units, the solution space is more broadly explored by producing a suite of successful models, and the likelihood of a particular conceptual geologic model can be compared. The Vaca Fault near Travis Air Force Base, California, can be successfully modeled as a normal or strike-slip fault, with the normal fault model being slightly more probable. It can also be modeled as a reverse fault, although this structural geologic configuration is highly unlikely given the realizations we explored.

Sign in / Sign up

Export Citation Format

Share Document