scholarly journals The Gravity Effect of Topography: A Comparison among Three Different Methods

2021 ◽  
Author(s):  
Carlo Iapige De Gaetani ◽  
Anna Maria Marotta ◽  
Riccardo Barzaghi ◽  
Mirko Reguzzoni ◽  
Lorenzo Rossi
Keyword(s):  
Gravity Data ◽  
Three Dimensional ◽  
Gravity Anomalies ◽  
Gravitational Effect ◽  
Terrain Correction ◽  
Bouguer Anomalies ◽  
Terrain Effect ◽  
Data Points ◽  
The One ◽  
Sharp Corners

In this paper, three different methods for computing the terrain correction have been compared. The terrain effect has been accounted for by using the standard right parallelepiped closed formula, the spherical tesseroid and the flat tesseroid formulas. Particularly, the flat tesseroid approximation is obtained by flattening the top and the bottom sides of the spherical tesseroid. Its gravitational effect can be computed as the gravitational effect of a polyhedron, i.e. a three-dimensional body with flat polygonal faces, straight edges and sharp corners or vertices. These three methods have been applied in the context of a Bouguer reduction scheme. Two tests were devised in the Alpine area in order to quantify possible discrepancies. In the first test, the terrain correction has been evaluated on a grid of points on the DTM. In the second test, Bouguer gravity anomalies were computed on sparse observed gravity data points. The results prove that the three methods are practically equivalent even in an area of rough topography though, in the second test, the Bouguer anomalies obtained by using the tesseroid and the flat tesseroid formulas have slightly smaller RMSs than the one obtained by applying the standard right parallelepiped formula.

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>

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.

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.

On: “Reduction of magnetic and gravity data on an arbitrary surface acquired in a region of high topographic relief” by B. K. Bhattacharyya and K. C. Chan (GEOPHYSICS, December 1977, p. 1411–1430)

Geophysics ◽  
1978 ◽  
Vol 43 (6) ◽  
pp. 1274-1275
Author(s):  
V. Courtillot ◽  
J. Ducruix ◽  
J. L. Le Mouël
Keyword(s):  
Potential Field ◽  
Gravity Data ◽  
Three Dimensional ◽  
Analytic Solutions ◽  
Helmholtz Equations ◽  
Data Points ◽  
The Subject ◽  
Christoffel Transformation ◽  
The Relationship ◽  
Topographic Relief

In their paper, Bhattacharryya and Chan address the problem of reduction of magnetic and gravity data on an arbitrary surface acquired in a region of high topographic relief. In their work, the authors are kind enough to mention our contribution to the solution of the problem of finding the sources responsible for an observed magnetic or gravity anomaly, using the general formalism of inverse problems (Courtillot et al, 1974). Unfortunately, however, the authors seem to be unaware of our other publications which are far more relevant to their subject. Courtillot et al (1973) solved the problem of continuation of a potential field measured on an uneven profile, using the Backus and Gilbert approach. Another reference relevant to this problem (solved by Bhattacharryya and Chan on p. 1424) is Parker and Klitgord (1972), who used the Schwartz‐Christoffel transformation. The work was extended to the case of three‐dimensional potential fields measured on an uneven surface by Ducruix et al (1974). Indeed, the development of our paper is strikingly similar to that of Bhattacharryya and Chan, although the method is quite different. In our paper, we give many illustrations of both theoretical and real cases, in which our method is seen to perform very well. We leave it to the reader to compare the results provided by both methods and to compare the methods themselves. In a third paper (Le Mouël et al, 1975), we generalized the method and showed how one could obtain excellent approximate analytic solutions of the Dirichlet and Neumann problems in the two‐dimensional case for a contour with any arbitrary shape. Finally, let us take the opportunity of this discussion to mention a review of the subject which appears in French in Courtillot (1977) and in English, much expanded, in Courtillot et al (1978). In this last paper, which should be of interest in solving a variety of geophysical problems, we show how our method allows one to continue a potential field measured on an entirely arbitrary set of data points in any number of dimensions for the various coordinate systems in which the Laplace and Helmholtz equations are separable. We also establish the relationship between our method and a generalization of the theory of generalized inverse matrices. One other relevant reference on that subject is Parker (1977). In the case of spherical coordinates, an application can be the continuation of satellite data, a problem studied by Bhattacharryya (1977).

Topography of the crust–mantle interface under the Western Superior craton from gravity data

2003 ◽  
Vol 40 (10) ◽  
pp. 1307-1320 ◽  
Author(s):  
B Nitescu ◽  
A R Cruden ◽  
R C Bailey
Keyword(s):  
Gravity Data ◽  
Three Dimensional ◽  
Gravity Anomalies ◽  
Gravity Inversion ◽  
Constant Density ◽  
Inversion Algorithm ◽  
Data Set ◽  
Mantle Boundary ◽  
Refraction Seismic ◽  
Superior Craton

The Moho undulations beneath the western part of the Archean Superior Province have been investigated with a three-dimensional gravity inversion algorithm for a single interface of constant density contrast. Inversion of the complete gravity data set produces unreal effects in the solution due to the ambiguity in the possible sources of some crustal gravity anomalies. To avoid these effects a censored gravity data set was used instead. The inversion results are consistent with reflection and refraction seismic data from the region and, therefore, provide a basis for the lateral correlation of the Moho topography between parallel seismic lines. The results indicate the existence of a major linear east–west-trending rise of the Moho below the metasedimentary English River subprovince, which is paralleled by crustal roots below the granite–greenstone Uchi and Wabigoon subprovinces. This correlation between the subprovincial structure at the surface and deep Moho undulations suggests that the topography of the crust–mantle boundary is related to the tectonic evolution of the Western Superior belts. Although certain features of the crust–mantle boundary are likely inherited from the accretionary and collisional stages of the Western Superior craton, gravity-driven processes triggered by subsequent magmatism and crustal softening may have played a role in both the preservation of those features, as well as in the development of new ones.

Correlation of gravity anomalies with Yellowknife Supergroup rocks, North Arm, Great Slave Lake

1980 ◽  
Vol 17 (11) ◽  
pp. 1506-1516 ◽  
Author(s):  
R. A. Gibb ◽  
M. D. Thomas
Keyword(s):  
Greenstone Belt ◽  
Gravity Data ◽  
Volcanic Rocks ◽  
Three Dimensional ◽  
Gravity Anomalies ◽  
Rock Density ◽  
The North ◽  
Great Slave Lake ◽  
Western Shore ◽  
Gravity Map

A gravity map compiled from observations made on the frozen surface of Great Slave Lake shows that the positive gravity anomaly associated with the Yellowknife greenstone belt extends offshore into the North Arm of the lake. On the western shore of Yellowknife Bay the axis of the anomaly coincides with mafic volcanic rocks of the Kam Formation. Offshore the axis continues southwards for about 10 km to the West Mirage Islands where it takes a dramatic turn to the southeast and continues for a further 60 km to the Outer Whaleback Rocks. Using the geology and rock density determinations on land for control, a three-dimensional geological model comprising a large number of prismatic blocks was derived from the gravity anomalies. In the model the simplifying assumption has been made that the greenstone belt is everywhere floored by granodiorite similar to the adjacent Western and South-east granodiorites. According to the model, mafic volcanic rocks of the Kam Formation are generally 1–3 km thick with a maximum thickness of 7 km at the mouth of Yellowknife Bay. Greywacke and mudstone of the Burwash Formation vary in thickness from 1 to 3 km. Locally these sedimentary rocks attain a thickness of 8 km but this is probably an overestimated value as they may very well be underlain by volcanic rocks of the Kam Formation. The presence of a third pluton of granodiorite flanking the belt to the southwest is also inferred from the gravity data. Previous seismic work indicated a greenstone basin with an average thickness of about 10 km. However, reexamination of the seismic records suggests that weak arrivals interpreted as originating from the base of the greenstone belt are more likely to be pulses associated with earlier arrivals.

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>

scholarly journals Usporedba različitih pristupa Bouguer-ove redukcije na temelju satelitskih gravimetrijskih podataka

Geofizika ◽  
2020 ◽  
Vol 37 (2) ◽  
pp. 237-261
Author(s):  
Fan Luo ◽  
Xin Tao ◽  
Guangming Fu ◽  
Chong Zhang ◽  
Kun Zhang ◽  
...  
Keyword(s):  
Gravity Anomaly ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Gravitational Effect ◽  
Seismic Profile ◽  
Bouguer Gravity ◽  
Seismic Profiles ◽  
Satellite Gravity ◽  
Free Air ◽  
Free Air Gravity

Satellite gravity data are widely used in the field of geophysics to study deep structures at the regional and global scales. These data comprise free-air gravity anomaly data, which usually need to be corrected to a Bouguer gravity anomaly for practical application. Bouguer reduction approaches can be divided into two methods based on the coordinate system: the spherical coordinates method (SBG) and the Cartesian coordinates method; the latter is further divided into the CEBG and CBG methods, which do and do not include the Earth’s curvature correction. In this paper, free-air gravity anomaly data from the eastern Tibetan Plateau and its adjacent areas were used as the basic data to compare the CBG, CEBG, and SBG Bouguer gravity correction methods. The comparison of these three Bouguer gravity correction methods shows that the effect of the Earth’s curvature on the gravitational effect increases with increasing elevation in the study area. We want to understand the inversion accuracy for the data obtained by different Bouguer gravity reduction approaches. The depth distributions of the Moho were obtained by the interface inversion of the Bouguer gravity anomalies obtained by the CBG, CEBG, and SBG, and active seismic profiles were used as references for comparison and evaluation. The results show that the depths of the Moho obtained by the SBG inversion are more consistent with the measured seismic profile depths. Therefore, the SBG method is recommended as the most realistic approach in the process of global or regional research employing gravity data.

scholarly journals Subsurface structure investigation of the United Arab Emirates using gravity data

2021 ◽  
Vol 13 (1) ◽  
pp. 262-271
Author(s):  
Hakim Saibi ◽  
Diab Bakri Hag ◽  
Mohammed Saeed Mohammed Alamri ◽  
Hamdan Abdo Ali
Keyword(s):  
United Arab Emirates ◽  
Crustal Structure ◽  
Gravity Data ◽  
Three Dimensional ◽  
Petroleum Exploration ◽  
Gravity Inversion ◽  
Subsurface Structure ◽  
Bouguer Anomalies ◽  
Groundwater Aquifers ◽  
3D Gravity

Abstract The crustal structure beneath the United Arab Emirates (UAE) is still relatively unknown. Here, we use regional gravity data to constrain the subsurface density distribution and structure of the crust of the UAE by applying diverse gravity derivatives methods such as horizontal derivative (HDR), analytic signal (AS), and tilt angle (TA) to analyze the subsurface structure and perform three-dimensional (3D) gravity inversion for imaging crustal structure from the surface down to 35 km depth. The results are compared with known geological regional structures and the location of the petroleum fields. The Bouguer anomalies range from −100.8 to 113.5 mGal. The 3D gravity inversion results and the maximum Bouguer values coincide with the ophiolitic Hajar mountains in the east and the successive anticlines (uplifted basement rocks) and synclines in different parts of UAE, which could be promising sites for future mining and petroleum exploration. Also, the 3D density model results and the minimum Bouguer anomalies are located over the Aruma Basin, eastern UAE Platform, and Low Central UAE Platform, which can be the places for deep groundwater aquifers. These new results from HDR, AS, and TA successfully identify known geological structures, especially in the eastern part of UAE.

scholarly journals Simple Formulae for Estimating the Depth and Dimension of Some Geological Bodies From Gravity Anomalies

2001 ◽  
Vol 6 (1) ◽  
pp. 71
Author(s):  
Sayyadul Arafin
Keyword(s):  
Gravity Anomalies ◽  
Horizontal Cylinder ◽  
Radial Symmetry ◽  
The Body ◽  
Simple Method ◽  
Depth Estimate ◽  
Data Points ◽  
The Mean ◽  
The One ◽  
Body Parameters

A simple method  for determining the depth and radius of some simple geological bodies such as a sphere and a horizontal cylinder  is presented. The method utilizes all of the data points to calculate the  body parameters.  The depth is estimated  from the mean of the depths determined from various fractions of the maximum anomaly and their corresponding widths. The radii of the two bodies are  calculated from the areas under the respective gravity profiles. Additionally , the radial symmetry of the gravity anomaly due to a sphere has been utilized to calculate the excess mass from a single profile. The radius determined from the excess mass is independent of depth estimate. This method of determining the radius is more robust than the one based on calculating the area.  An error analysis showing the effects of truncation and zero-level errors on the estimated radius  is  also discussed.  

Sign in / Sign up

Export Citation Format

Share Document