Linear inversion of gravity data using the spectral expansion method

Geophysics ◽  
1985 ◽  
Vol 50 (5) ◽  
pp. 820-824 ◽  
Author(s):  
Rene E. Chavez ◽  
George D. Garland
Keyword(s):  
Inverse Problem ◽  
Gravity Data ◽  
Linear Equations ◽  
Thin Sheet ◽  
Gravity Anomalies ◽  
Expansion Method ◽  
Spectral Expansion ◽  
Uniform Density ◽  
Linear Inversion ◽  
Prismatic Structure

Inversion of gravity anomalies in terms of an anomalous mass distribution with irregular outline but uniform density leads to a nonlinear inverse problem. An alternative approach based on the thin‐sheet approximation can, however, be formulated as a linear inverse problem provided the structures are two‐dimensional. The anomalous mass is represented by a thin sheet, which is located at a depth [Formula: see text] and is divided into M strips with N < M data points. Thus, an underdetermined system of linear equations is obtained which is solved by the spectral expansion method for the surface density distribution of each segment. This set of parameters is then transformed into a prismatic structure with variable depth but uniform density. The modeling procedure involves a noniterative method. A gravity problem is investigated, and the solution obtained compares well with previous interpretations.

Linear inversion of gravity data for 3-D density distributions

Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1354-1364 ◽  
Author(s):  
Glenn W. Bear ◽  
Haydar J. Al‐Shukri ◽  
Albert J. Rudman
Keyword(s):  
Density Distribution ◽  
Gravity Data ◽  
Linear Equations ◽  
Large Data ◽  
Linear Inversion ◽  
Solution Quality ◽  
Data Set ◽  
User Input ◽  
Statistical Measures ◽  
Unknown Density

We have developed an improved Levenburg‐Marquart technique to rapidly invert Bouguer gravity data for a 3-D density distribution as a source of the observed field. This technique is designed to replace tedious forward modeling with an automatic solver that determines density models constrained by geologic information supplied by the user. Where such information is not available, objective models are generated. The technique estimates the density distribution within the source volume using a least‐squares inverse solution that is obtained iteratively by singular value decomposition using orthogonal decomposition of matrices with sequential Householder transformations. The source volume is subdivided into a series of right rectangular prisms of specified size but of unknown density. This discretization allows the construction of a system of linear equations relating the observed gravity field to the unknown density distribution. Convergence of the solution to the system is tightly controlled by a damping parameter which may be varied at each iteration. The associated algorithm generates statistical measures of solution quality not available with most forward methods. Along with the ability to handle large data sets within reasonable time constraints, the advantages of this approach are: (1) the ease with which pre‐existing geological information can be included to constrain the solution, (2) its minimization of subjective user input, (3) the avoidance of difficulties encountered during wavenumber domain transformations, and (4) the objective nature of the solution. Application to a gravity data set from Hamilton County, Indiana, has yielded a geologically reasonable result that agrees with published models derived from interpretation of gravity, magnetic, seismic, and drilling data.

Inversion of the power spectrum from gravity anomalies of prismatic bodies

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1698-1703 ◽  
Author(s):  
Juan García‐Abdeslem
Keyword(s):  
Inverse Problem ◽  
Power Spectrum ◽  
Ridge Regression ◽  
Field Data ◽  
Gravity Data ◽  
Source Parameters ◽  
Gravity Anomalies ◽  
Independent Functions ◽  
Prismatic Bodies ◽  
And Inversion

I develop methods for modeling and inversion of gravity anomalies produced by prismatic bodies. The forward problem is solved in the wavenumber domain, where the power spectrum of the gravity anomaly is given by the product of independent functions that describe depth, thickness, horizontal dimensions, and density of the source body. The solution of the inverse problem is based on an iterative ridge‐regression algorithm, starting from an initial trial of the geometry and density of the source body. This procedure is assessed in a number of examples with both numeric and field data spectra. The method is first tested on the spectrum obtained from simple gravity anomalies and is found to be capable of recovering the source parameters. Next the method is applied to gravity data from a survey made near Noranda, Quebec, Canada. This interpretation compares favorably with drillhole data and provides an estimation of the mass of the source body which is similar to previous estimates.

Nonlinear inversion of isostatic residual gravity data from Montage Basin, northern Gulf of California

Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. G45-G55 ◽  
Author(s):  
Juan García-Abdeslem
Keyword(s):  
Inverse Problem ◽  
Gravity Anomaly ◽  
Gulf Of California ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Nonlinear Inversion ◽  
Residual Gravity ◽  
Sedimentary Sequences ◽  
Residual Gravity Anomaly ◽  
Northwestern Mexico

The flexural isostatic response to surface loads is used to estimate the crustal thickness in northwestern Mexico and Southwestern USA. This estimate is used to compute an isostatic regional gravity, which, subtracted from Bouguer gravity anomalies, led to the isostatic residual gravity anomaly at Montage Basin. This basin is located between the southern portion of the Mexicali Valley and the northern Gulf of California, it roughly has an extension of [Formula: see text] wide, and it shows a gravity minimum reaching approximately [Formula: see text]. Montage Basin is within the extensional province of the Gulf of California, where rifting is currently an ongoing geologic process, and deep exploratory wells drilled by Petróleos Mexicanos have shown that the basin accommodates thick sedimentary sequences greater than 5 km. The interpretation of the isostatic residual gravity anomaly is considered as a nonlinear inverse problem, constrained using density as a function of depth derived from Gardner’s equation applied to dual time [Formula: see text]-logs, assuming isostatic equilibrium and considering the basin as a subsurface load that is compensated at depth by a mass of unknown shape and density. The outcome of the inverse problem suggests that Montage Basin accommodates as much as 7.5 km thick sedimentary sequences and a compensating mass at a minimum depth of 13 km.

scholarly journals Joint Gravity and Seismic Reflection Methods to Characterize the Deep Aquifers in Arid Ain El Beidha Plain (Central Tunisia, North Africa)

Water ◽  
2021 ◽  
Vol 13 (9) ◽  
pp. 1310
Author(s):  
Hajer Azaiez ◽  
Hakim Gabtni ◽  
Mourad Bédir
Keyword(s):  
Seismic Reflection ◽  
Large Scale ◽  
Seismic Analysis ◽  
Gravity Data ◽  
Geophysical Methods ◽  
Gravity Anomalies ◽  
Integrated Approach ◽  
Central Tunisia ◽  
Deep Aquifers ◽  
Advanced Analysis

Electric resistivity sounding and tomography, as well as electromagnetic sounding, are the classical methods frequently used for hydrogeological studies. In this work, we propose the development and implementation of an original integrated approach using the unconventional hydro–geophysical methods of gravity and seismic reflection for the fast, large–scale characterization of hydrogeological potential using the Ain El Beidha plain (central Tunisia) as an analogue. Extending the values of vintage petroleum seismic reflection profiles and gravity data, in conjunction with available geological and hydrogeological information, we performed an advanced analysis to characterize the geometry of deep tertiary (Oligocene and Eocene) aquifers in this arid area. Residual and tilt angle gravity maps revealed that most gravity anomalies have a short wavelength. The study area was mainly composed of three major areas: the Oued Ben Zitoun and Ain El Beidha basins, which are both related to negative gravity trends corresponding to low–density subsiding depocenters. These basins are separated by an important NE–SW trend called “El Gonna–J. El Mguataa–Kroumet Zemla” gravity high. Evaluation of the superposition of detected lineaments and Euler deconvolution solutions’ maps showed several NE–SW and N–S relay system faults. The 3D density inversion model using a lateral and vertical cutting plane suggested the presence of two different tectonic styles (thin VS thick). Results from the gravity analysis were in concordance with the seismic analysis. The deep Oligocene and Eocene seismic horizons were calibrated to the hydraulic wells and surrounding outcrops. Oligocene and Eocene geological reservoirs appear very fractured and compartmented. The faulting network also plays an important role in enhancing groundwater recharge process of the Oligocene and Eocene aquifers. Finally, generated isochron maps provided an excellent opportunity to develop future comprehensive exploration surveys over smaller and more favorable areas’ sub–basins.

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 Basement Mapping of the Fucino Basin in Central Italy by ITRESC Modeling of Gravity Data

Geosciences ◽  
2021 ◽  
Vol 11 (10) ◽  
pp. 398
Author(s):  
Federico Cella ◽  
Rosa Nappi ◽  
Valeria Paoletti ◽  
Giovanni Florio
Keyword(s):  
Seismic Activity ◽  
Site Response ◽  
Gravity Data ◽  
A Priori ◽  
Central Italy ◽  
Gravity Anomalies ◽  
Gravity Modeling ◽  
Ground Shaking ◽  
Sedimentary Evolution ◽  
Fucino Basin

Sediments infilling in intermontane basins in areas with high seismic activity can strongly affect ground-shaking phenomena at the surface. Estimates of thickness and density distribution within these basin infills are crucial for ground motion amplification analysis, especially where demographic growth in human settlements has implied increasing seismic risk. We employed a 3D gravity modeling technique (ITerative RESCaling—ITRESC) to investigate the Fucino Basin (Apennines, central Italy), a half-graben basin in which intense seismic activity has recently occurred. For the first time in this region, a 3D model of the Meso-Cenozoic carbonate basement morphology was retrieved through the inversion of gravity data. Taking advantage of the ITRESC technique, (1) we were able to (1) perform an integration of geophysical and geological data constraints and (2) determine a density contrast function through a data-driven process. Thus, we avoided assuming a priori information. Finally, we provided a model that honored the gravity anomalies field by integrating many different kinds of depth constraints. Our results confirmed evidence from previous studies concerning the overall shape of the basin; however, we also highlighted several local discrepancies, such as: (a) the position of several fault lines, (b) the position of the main depocenter, and (c) the isopach map. We also pointed out the existence of a new, unknown fault, and of new features concerning known faults. All of these elements provided useful contributions to the study of the tectono-sedimentary evolution of the basin, as well as key information for assessing the local site-response effects, in terms of seismic hazards.

scholarly journals DETERMINATION AND MAPPING OF THE TERRESTRIAL GRAVITY ANOMALIES IN THE MOUNTAINOUS AREAS OF IRAQ

2021 ◽  
Vol 6 (24) ◽  
pp. 213-225
Author(s):  
Shazad Jamal Jalal ◽  
Tajul Ariffin Musa ◽  
Ami Hassan Md Din ◽  
Wan Anom Wan Aris
Keyword(s):  
Gravity Data ◽  
Physical Geodesy ◽  
Gravity Anomalies ◽  
Gravity Survey ◽  
Mountainous Areas ◽  
Entire Study ◽  
Terrestrial Gravity ◽  
The Past ◽  
Regional Gravity ◽  
Entire Study Area

Gravity data and computing gravity anomalies are regarded as vital for both geophysics and physical geodesy fields. The mountainous areas of Iraq are characterized by the lack of regional gravity data because gravity surveys are rarely performed in the past four decades due to the Iraq-Iran war and the internal unstable political situation of this particular region. In addition, the formal map of the available terrestrial gravity which was published by the French Database of Bureau Gravimetrique International (International Gravimetric Bureau-in English) (BGI), introduces Iraq and the study area as a remote area and in white color because of the unavailability of gravity data. However, a dense and local (not regional) gravity data is available which was conducted by geophysics researchers 13 years ago. Therefore, the regional gravity survey of 160 gravity points was performed by the authors at an average 11 km apart, which was covers the whole area of Sulaymaniyah Governorate (part of the mountainous areas of Iraq). In spite of Although the risk of mine fields within the study area, suitable safe routes as well as a helicopter was used for the gravity survey of several points on the top of mountains. The survey was conducted via Lacoste and Romberg geodetic gravimeter and GPS handheld. The objective of the study is to determine and map the gravity anomalies for the entire study area, the data of which would assist different geosciences applications.

scholarly journals Structural model of the Northern Latium volcanic area constrained by MT, gravity and aeromagnetic data

1997 ◽  
Vol 40 (5) ◽  
Author(s):  
P. Capuano ◽  
G. Florio ◽  
P. Gasparini
Keyword(s):  
Structural Model ◽  
Gravity Data ◽  
Gravity Anomalies ◽  
Volcanic Area ◽  
Aeromagnetic Data ◽  
Geothermal Exploration ◽  
Magnetotelluric Soundings ◽  
Aeromagnetic Anomalies ◽  
Deep Boreholes

The results of about 120 magnetotelluric soundings carried out in the Vulsini, Vico and Sabatini volcanic areas were modeled along with Bouguer and aeromagnetic anomalies to reconstruct a model of the structure of the shallow (less than 5 km of depth) crust. The interpretations were constrained by the information gathered from the deep boreholes drilled for geothermal exploration. MT and aeromagnetic anomalies allow the depth to the top of the sedimentary basement and the thickness of the volcanic layer to be inferred. Gravity anomalies are strongly affected by the variations of morphology of the top of the sedimentary basement, consisting of a Tertiary flysch, and of the interface with the underlying Mesozoic carbonates. Gravity data have also been used to extrapolate the thickness of the neogenic unit indicated by some boreholes. There is no evidence for other important density and susceptibility heterogeneities and deeper sources of magnetic and/or gravity anomalies in all the surveyed area.

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

&lt;p&gt;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.&lt;/p&gt;&lt;p&gt;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.&lt;/p&gt;&lt;p&gt;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.&lt;/p&gt;&lt;p&gt;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 &amp;#8211; ocean-ocean or ocean-continental subduction &amp;#8211; that drives a significantly different dynamic in the wedge area.&lt;/p&gt;

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.

Sign in / Sign up

Export Citation Format

Share Document