Approximating edges of source bodies from magnetic or gravity anomalies

Geophysics ◽  
1986 ◽  
Vol 51 (7) ◽  
pp. 1494-1498 ◽  
Author(s):  
Richard J. Blakely ◽  
Robert W. Simpson
Keyword(s):  
Gravity Anomaly ◽  
Gravity Data ◽  
Mass Density ◽  
Gravity Anomalies ◽  
Horizontal Gradient ◽  
Plan View ◽  
Crustal Rocks

Cordell and Grauch (1982, 1985) discussed a technique to estimate the location of abrupt lateral changes in magnetization or mass density of upper crustal rocks. The final step of their procedure is to identify maxima on a contoured map of horizontal gradient magnitudes. We attempt to automate their final step. Our method begins with gridded magnetic or gravity anomaly data and produces a plan view of inferred boundaries of magnetic or gravity sources. The method applies to both local surveys and to continent‐wide compilations of magnetic and gravity data (e.g., Zietz, 1982; Simpson et al., 1983a; Kane et al., 1982).

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>

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;

4D gravity changes associated with the 2005 eruption of Sierra Negra volcano, Galápagos

Geophysics ◽  
2008 ◽  
Vol 73 (6) ◽  
pp. WA29-WA35 ◽  
Author(s):  
Nathalie Vigouroux ◽  
Glyn Williams-Jones ◽  
William Chadwick ◽  
Dennis Geist ◽  
Andres Ruiz ◽  
...  
Keyword(s):  
Melt Inclusion ◽  
Gravity Data ◽  
Mass Density ◽  
Remote Sensing Data ◽  
Gravity Anomalies ◽  
Gas Content ◽  
Residual Gravity ◽  
Caldera Floor ◽  
Sierra Negra Volcano ◽  
Sierra Negra

Sierra Negra volcano, the most voluminous shield volcano in the Galápagos archipelago and one of the largest basaltic calderas in the world, erupted on October 22, 2005 after more than [Formula: see text] of quiescence. GPS and satellite radar interferometry (InSAR) monitoring of the deformation of the caldera floor in the months prior to the eruption documented extraordinary inflation rates [Formula: see text]. The total amount of uplift recorded since monitoring began in 1992 approached [Formula: see text] at the center of the caldera over the eight days of the eruption the caldera floor deflated a maximum of 5 m and subsquently renewed its inflation, but at a decelerating rate. To gain insight into the nature of the subsurface mass/density changes associated with the deformation, gravity measurements performed in 2005, 2006, and 2007 are compared to previous measurements from 2001-2002 when the volcano underwent a period of minor deflation and magma withdrawal.The residual gravity decrease between 2001-2002 and 2005 is among the largest ever recorded atan active volcano (−950 μGal) and suggests that inflation was accompanied by a relative density decrease in the magmatic system. Forward modeling of the residual gravity data in 4D (from 2002 to 2005) gives an estimate of the amount of vesiculation in the shallow sill required to explain the observed gravity variations. Geochemical constraints from melt inclusion and satellite remote-sensing data allow us to estimate the pre-eruptive gas content of the magma and place constraints on the thickness of the gas-rich sill necessary to produce the gravity anomalies observed. Results suggest that reasonable sill thicknesses [Formula: see text] and bubble contents (10–50 volume %) can explain the large decrease in residual gravity prior to eruption. Following the eruption (2006 and 2007), the deformation and gravity patterns suggest re-equilibration of the pressure regime in the shallow magma system via a renewed influx of relatively gas-poor magma into the shallow parts of the system.

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 Study of salt structures from gravity and seismic data in Santos Basin, Brazil

2016 ◽  
Vol 55 (3) ◽  
Author(s):  
Renata Regina Constantino Regina Constantino ◽  
Eder Cassola Molina ◽  
Iata Anderson de Souza
Keyword(s):  
Gravity Anomaly ◽  
Crustal Structure ◽  
Gravity Data ◽  
Geophysical Methods ◽  
Gravity Anomalies ◽  
Sedimentary Basins ◽  
Final Conclusion ◽  
Additional Information ◽  
Salt Layer ◽  
Modelling Studies

Seismic is one of the main methods used for the identification of structures and stratigraphic studies in sedimentary basins. In the Santos Basin, numerous 2D and 3D seismic surveys are being conducted in order to get a better ima-ge of the geological section to depths beyond the base of salt layer. Crustal modelling studies that make joint interpretation of seismic and gravity data are found in the literature, however there are few studies that relate gravity anomalies directly to salt structures. This work aims to associate gravity anomalies with salt structures from seismic and gravimetric interpretation. For studies aimed to model the crustal structure from gravity field data, the knowledge of two major discontinuities is required, the basement and the Moho. Such interfaces are often not easily seen by seismic and so, during this study, they were found by different methods involving analysis of gravity anomalies. The other interfaces involving density contrasts were analyzed based on seismic interpretation. The results showed that the obtained seismic geological interpretations may provide additional information when compared to gravity anomaly data. In all the modelled profiles, some geological information of the Santos Basin that are not visible in the seismic, could be interpreted according to the geological model and the adjustment of gravity anomaly curves. As a final conclusion of this work, it is suggested that the combined analysis of the two cited geophysical methods, can provide important information about the crustal structure and to assist in modelling the salt layer.

Stable inversion of gravity anomalies of sedimentary basins with nonsmooth basement reliefs and arbitrary density contrast variations

Geophysics ◽  
1999 ◽  
Vol 64 (3) ◽  
pp. 754-764 ◽  
Author(s):  
Valéria C. F. Barbosa ◽  
João B. C. Silva ◽  
Walter E. Medeiros
Keyword(s):  
Gravity Anomaly ◽  
Gravity Data ◽  
A Priori ◽  
Bouguer Anomaly ◽  
Gravity Anomalies ◽  
Sedimentary Basins ◽  
Upper And Lower Bounds ◽  
Density Contrast ◽  
Good Precision ◽  
Rift Basins

We present a new, stable method for interpreting the basement relief of a sedimentary basin which delineates sharp discontinuities in the basement relief and incorporates any law known a priori for the spatial variation of the density contrast. The subsurface region containing the basin is discretized into a grid of juxtaposed elementary prisms whose density contrasts are the parameters to be estimated. Any vertical line must intersect the basement relief only once, and the mass deficiency must be concentrated near the earth’s surface, subject to the observed gravity anomaly being fitted within the experimental errors. In addition, upper and lower bounds on the density contrast of each prism are introduced a priori (one of the bounds being zero), and the method assigns to each elementary prism a density contrast which is close to either bound. The basement relief is therefore delineated by the contact between the prisms with null and nonnull estimated density contrasts, the latter occupying the upper part of the discretized region. The method is stabilized by introducing constraints favoring solutions having the attributes (shared by most sedimentary basins) of being an isolated compact source with lateral borders dipping either vertically or toward the basin center and having horizontal dimensions much greater than its largest vertical dimension. Arbitrary laws of spatial variations of the density contrast, if known a priori, may be incorporated into the problem by assigning suitable values to the nonnull bound of each prism. The proposed method differs from previous stable methods by using no smoothness constraint on the interface to be estimated. As a result, it may be applied not only to intracratonic sag basins where the basement relief is essentially smooth but also to rift basins whose basements present discontinuities caused by faults. The method’s utility in mapping such basements was demonstrated in tests using synthetic data produced by simulated rift basins. The method mapped with good precision a sequence of step faults which are close to each other and present small vertical slips, a feature particularly difficult to detect from gravity data only. The method was also able to map isolated discontinuities with large vertical throw. The method was applied to the gravity data from Reco⁁ncavo basin, Brazil. The results showed close agreement with known geological structures of the basin. It also demonstrated the method’s ability to map a sequence of alternating terraces and structural lows that could not be detected just by inspecting the gravity anomaly. To demostrate the method’s flexibility in incorporating any a priori knowledge about the density contrast variation, it was applied to the Bouguer anomaly over the San Jacinto Graben, California. Two different exponential laws for the decrease of density contrast with depth were used, leading to estimated maximum depths between 2.2 and 2.4 km.

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 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 A least-squares minimization approach to interpret gravity anomalies caused by a 2D thick, vertically faulted slab

2019 ◽  
Vol 49 (3) ◽  
pp. 229-247
Author(s):  
El-Sayed Abdelrahman ◽  
Mohamed Gobashy
Keyword(s):  
Least Squares ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Central Valley ◽  
Real Data ◽  
Horizontal Gradient ◽  
Vertical Fault ◽  
Independent Observations ◽  
Minimization Approach ◽  
Least Squares Minimization

Abstract We present a least-squares minimization approach to estimate simultaneously the depth to and thickness of a buried 2D thick, vertically faulted slab from gravity data using the sample spacing – curves method or simply s-curves method. The method also provides an estimate for the horizontal location of the fault and a least-squares estimate for the density contrast of the slab relative to the host. The method involves using a 2D thick vertical fault model convolved with the same finite difference second horizontal gradient filter as applied to the gravity data. The synthetic examples (noise-free and noise affected) are presented to illustrate our method. The test on the real data (Central Valley of Chile) and the obtained results were consistent with the available independent observations and the broader geological aspects of this region.

Sign in / Sign up

Export Citation Format

Share Document