NON-STANDARD GRAVITY ANOMALIES OVER SEDIMENTARY STRUCTURES

1976 ◽  
Vol 16 (1) ◽  
pp. 117
Author(s):  
A. J. Flavelle ◽  
W. Anfiloff

Calculations of the gravity effect of sedimentary structures can now be easily and rapidly made using digital computers. When density variations are normal, they generate gravity profiles which are easily predictable. However, if the general rule regarding increase of density with depth is not fulfilled then anomaly variations can be unexpected. The standard profile for a fault is a change in Bouguer anomaly level, but in some circumstances it can be marked by a positive (or negative) anomaly and no overall change in level. The rather unexpected type of gravity anomaly generated from fairly standard geological sections emphasises the need to carry out modelling calculations based on the known or suspected geology.A gravity anomaly across the boundary between the Galilee and Drummond Basins in central Queensland is analysed. Part of the structure has been determined seismically, and a thickening sedimentary sequence appears to have caused a reverse density contrast with the basement resulting in a positive gravity anomaly. However, the flank of this anomaly has been affected by the negative contribution of a granite body, resulting in a displaced gradient, a minor low which is also displaced, and another minor high and low which could be mistakenly attributed to near-surface structure. A further complication is that a topographic feature occurs at this point of interaction, strongly affecting the multiple density Bouguer profiles, and thus the interpretation. The successful interpretation of this type of situation requires accurate gravity data, its accurate positioning relative to the seismic section, and gravity modelling facilities which allow a rapid trial-and-error approach to modelling.

Geophysics ◽  
1999 ◽  
Vol 64 (3) ◽  
pp. 754-764 ◽  
Author(s):  
Valéria C. F. Barbosa ◽  
João B. C. Silva ◽  
Walter E. Medeiros

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.


2021 ◽  
Author(s):  
◽  
Alistair Stronach

<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>


2021 ◽  
Author(s):  
Mirko Scheinert ◽  
Philipp Zingerle ◽  
Theresa Schaller ◽  
Roland Pail ◽  
Martin Willberg

&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;


2014 ◽  
Vol 4 (1) ◽  
Author(s):  
L.E. Sjöberg

AbstractTraditional gravimetric geoid determination relies on Stokes’ formula with removal and restoration of the topographic effects. It is shown that this solution is in error of the order of the quasigeoid-to-geoid difference, which is mainly due to incomplete downward continuation (dwc) of gravity from the Earth’s surface to the geoid. A slightly improved estimator, based on the surface Bouguer gravity anomaly, is also biased due to the imperfect harmonic dwc the Bouguer anomaly. Only the third estimator,which uses the (harmonic) surface no-topography gravity anomaly, is consistent with the boundary condition and Stokes’ formula, providing a theoretically correct geoid height. The difference between the Bouguer and no-topography gravity anomalies (on the geoid or in space) is the “secondary indirect topographic effect”, which is a necessary correction in removing all topographic signals.


Geophysics ◽  
1993 ◽  
Vol 58 (12) ◽  
pp. 1785-1791 ◽  
Author(s):  
El‐Sayed M. Abdelrahman ◽  
Hesham M. El‐Araby

The gravity anomaly expression produced by most geologic structures can be represented by a continuous function in both shape (shape factor) and depth variables with an amplitude coefficient related to the mass. Correlation factors between successive least‐squares residual gravity anomalies from a buried vertical cylinder, horizontal cylinder, and sphere are used to determine the shape and depth of the buried geologic structure. For each shape factor value, the depth is determined automatically from the correlation value. The computed depths are plotted against the shape factor representing a continuous correlation curve. The solution for the shape and depth of the buried structure is read at the common intersection of correlation curves. This method can be applied to a Bouguer anomaly profile consisting of a residual component caused by local structure and a regional component. This is a powerful technique for automatically separating the Bouguer data into residual and regional polynomial components. This method is tested on theoretical examples and a field example. In both cases, the results obtained are in good agreement with drilling results.


Author(s):  
A. Tugi ◽  
A. H. M. Din ◽  
K. M. Omar ◽  
A. S. Mardi ◽  
Z. A. M. Som ◽  
...  

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.


2016 ◽  
Vol 55 (3) ◽  
Author(s):  
Renata Regina Constantino Regina Constantino ◽  
Eder Cassola Molina ◽  
Iata Anderson de Souza

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.


1972 ◽  
Vol 9 (8) ◽  
pp. 942-959 ◽  
Author(s):  
J. M. Woodside

Detailed maps of free-air, Bouguer, and residual gravity anomalies for a survey area 250 km wide across the Mid-Atlantic Ridge between 45° and 46 °N have been compiled. The Bouguer anomaly was terrain-corrected to a radius of 40 km. The residual anomaly was computed from the terrain-corrected Bouguer anomaly using an empirical linear relationship between the Bouguer anomaly and the bathymetry to predict a 'regional' Bouguer anomaly from the depth data. North–south and east–west trends in the gravity data are enhanced in the residual anomaly; and it is suggested that at least one short east–west transform fault may offset the ridge in a right-lateral sense. The offset is presumably a response to a change in sea-floor spreading direction from west–northwest/east–southeast to west/east about 10 m.y. ago. A change in spreading rate may have occurred at the same time. A difference in accretion rate on either side of the ridge axis is indicated by asymmetry in the gravity data and by differences in the topographic compensation across the axis. Variations in the relationship of terrain-corrected Bouguer anomaly to bathymetry within the survey area suggest that a density deficiency or buoyant forces in the upper mantle are responsible for the overall elevation of the crestal mountain region but that the topography of the high-fractured plateau may be partially compensated by undulations of the crust–mantle interface.


2021 ◽  
Author(s):  
Alistair Stronach

<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>


Sign in / Sign up

Export Citation Format

Share Document