scholarly journals A NOVEL REGIONAL-RESIDUAL SEPARATION APPROACH FOR GRAVITY DATA THROUGH CRUSTAL MODELING

2018 ◽  
Vol 36 (4) ◽  
pp. 1
Author(s):  
Nelson Ribeiro Filho ◽  
Cristiano Mendel Martins ◽  
Renata de Sena Santos

ABSTRACT. Gravity anomalies normally contain information of all sources beneath Earth’s surface. Once residual anomalies exhibit information about the main target, the knowledge of this specific residual signal is extremely important to interpretation. To find this signal, it’s necessary to perform regional-residual separation. We present here a new approach of separation by using gravity crustal modeling. We divide the surface in prisms, with density given by GEMMA. We calculate the regional signal, assuming Earth’s crust can be the source of observed anomaly. This methodology was applied on Barreirinhas basin-Brazil. Its formation is related to geologic events in South America-Africa break. Besides, the complex geology is the main obstacle on finding the residual anomaly. We compare our methodology with robust-polynomial fitting and spectral-analysis. They were not able to identify the residual anomaly. Main trouble relies on absence of crust information. Those kind of environment usually requires forward modeling and/or gravity inversion. On the other hand, our approach considers all crust’s parameters. Then the difficulty on choosing the residual no longer exists. The residual anomaly follows a geologic pattern. The crustal depocenter was mapped between structural faults. Therefore, our results satisfies the main expectation and are extremely linked to Barreirinhas basin’s geological background. We recommend this separation procedure, once Earth’s crustal model and gravity data are available for all planet.Keywords: Gravity modeling; GEMMA model; Barreirinhas basin; residual anomaly. RESUMO. Anomalias gravimétricas contêm informações de todas as fontes na superfície terrestre. Uma vez que anomalias residuais exibem informações sobre alvos principais, o conhecimento desse específico sinal residual é extremamente importante para interpretação. Para encontrá-lo, é necessário realizar separação regional-residual. Apresentamos aqui uma nova abordagem de separação utilizando a modelagem gravimétrica crustal. Discretizamos a superfície em prismas, com densidade fornecida pelo modelo GEMMA. Calculamos o sinal regional, assumindo que a crosta terrestre é a fonte da anomalia observada. Aplicamos esta metodologia na bacia de Barreirinhas - Brasil, que tem sua formação relacionada aos eventos geológicos de separação da América do Sul e África. Além disso, a complexidade geológica é considerada o principal obstáculo para encontrar esta anomalia residual. Comparamos nossa metodologia com Ajuste Polinomial Robusto e Análise Espectral. Essas técnicas não foram capazes de identificar a anomalia residual. O principal problema se dá pela ausência de informações acerca da crosta. Para esse ambiente, geralmente requer modelagem direta e/ou inversão geofísica. Por outro lado, nossa abordagem considera todos os parâmetros crustais e a dificuldade em escolher o residual deixa de existir. A anomalia residual apresenta um padrão geológico. O depocentro crustal foi mapeado entre falhas estruturais. Nossos resultados satisfazem a expectativa principal e estão extremamente ligados ao cenário geológico da bacia. Recomendamos este procedimento de separação, uma vez que os modelos crustais e dados gravimétricos estão disponíveis para todo o planeta.Palavras-chave: Modelagem gravimétrica; modelo GEMMA; bacia de Barreirinhas; anomalia residual

Geosciences ◽  
2021 ◽  
Vol 11 (10) ◽  
pp. 398
Author(s):  
Federico Cella ◽  
Rosa Nappi ◽  
Valeria Paoletti ◽  
Giovanni Florio

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.


Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. G15-G23
Author(s):  
Andrea Vitale ◽  
Domenico Di Massa ◽  
Maurizio Fedi ◽  
Giovanni Florio

We have developed a method to interpret potential fields, which obtains 1D models by inverting vertical soundings of potential field data. The vertical soundings are built through upward continuation of potential field data, measured on either a profile or a surface. The method assumes a forward problem consisting of a volume partitioned in layers, each of them homogeneous and horizontally finite, but with the density changing versus depth. The continuation errors, increasing with the altitude, are automatically handled by determining the coefficients of a third-order polynomial function of the altitude. Due to the finite size of the source volume, we need a priori information about the total horizontal extent of the volume, which is estimated by boundary analysis and optimized by a Markov chain process. For each sounding, a 1D inverse problem is independently solved by a nonnegative least-squares algorithm. Merging of the several inverted models finally yields approximate 2D or 3D models that are, however, shown to generate a good fit to the measured data. The method is applied to synthetic models, producing good results for either perfect or continued data. Even for real data, i.e., the gravity data of a sedimentary basin in Nevada, the results are interesting, and they are consistent with previous interpretation, based on 3D gravity inversion constrained by two gamma-gamma density logs.


2003 ◽  
Vol 40 (10) ◽  
pp. 1307-1320 ◽  
Author(s):  
B Nitescu ◽  
A R Cruden ◽  
R C Bailey

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.


Author(s):  
S. H. Anikeyev ◽  
S. M. Bahriy ◽  
B. B. Hablovskiy

In accordance with the purpose of geophysical exploration, the gravity data interpretation is aimed at prospecting mineral resources which is based on the study of the geological cross-section structure. The task of quantitative interpretation, which uses methods of gravity modeling and gravity inversion, is the modelling of a gravity field (gravity modeling) and of a density structure of geological environments (gravity inversion). The article presents the definition and steps of the gravity data modelling technique. This technique is based on the construction of an informal sequence of equivalent solutions. The technological and geological features of methods for modelling the density structure of complex geological environments are given; among them geological content, consistency with a priori data and the subordination of modelling to geological hypotheses are important. The topicality and methods of simulation modelling are outlined. The purpose of simulation modelling is to study the properties of gravity inversion in the general formulation, as well as to evaluate the degree of detail and reliability of the methods and technologies of gravity modelling, which claim to be an effective solution to geological problems. The example of structural simulation testing of the methods of informal sequence of equivalent solutions and its computer technologies shows that a complex interpretation of seismic and gravity measurements data enables the creation of detailed density models of structural cross-sections. The ways of increasing the veracity of gravity data modelling of structural cross-sections have been studied. It is revealed that the best approximation of the regional background is an inclined plane which approximates the observed field of gravity according to characteristic pickets over the research areas that are better studied. The increase in the veracity of modeling can also be achieved by rebuilding the near side zones in the structural type models in an interactive process of solving structural gravity inversion problems. Substantive modeling depends primarily on the experience of the interpreter since computer technologies for gravity modeling and gravity inversion are merely an interpretation tool.


2021 ◽  
Author(s):  
Igor Ognev ◽  
Jörg Ebbing ◽  
Peter Haas

<p>A new crustal model of the Volga-Uralian subcraton was built. The compilation of the model was subdivided in two steps: (1) inverse gravity modeling followed by (2) thorough forward gravity modeling.</p><p>For inverse gravity modeling GOCE gravity gradients were used. The effect of the Earth sphericity was taken into account by using tesseroids. Density contrasts between crust and mantle were varied laterally according to the tectonic units present in the region.  The model is constrained by the available seismic data including receiver function studies, and deep reflection and refraction profiles.</p><p>The Moho discontinuity obtained during the gravity inversion was consequently modified, and complemented by the sedimentary cover, upper crust, lower crust, and lithospheric mantle layers in the process of forward gravity modeling. Obtained model showed crustal thickness variation from 34 to more than 55 km in some areas. The thinnest crust with the thickness below 40 km appeared on the Pericaspian basin with the thickest sedimentary column. A relatively thin crust was found along the central Russia rift system, while the thickest crust is located underneath Ural Mountains as well as in the center of the Volga-Uralian subcraton. In both areas the crustal thickness exceeds 50 km. At the same time, the gravity misfit of ca. 95 mGal between the measured Bouguer gravity anomaly and forward calculated gravity field was revealed in the central area of the Volga-Uralian subcraton. This misfit was interpreted and modeled as high-density lower crust which can possibly represent an underplated material.</p><p>In the end, the new crustal model of Volga-Uralian subcraton respects the gravity and seismic constraints, and reflects the main geological features of the region. This model will be used for further geothermal analysis of the area.</p>


Author(s):  
P. Gryshchuk

The article is devoted to the description of the new algorithm for inversion of gravitational anomalies.The principle of the geological objects mass displacement is used for definition of its geometry. In the grid model, the movement of blocks is performed in four directions for a two-dimensional model. The motion and changing of density is determined in a random style. This algorithm defined the geometry of three rectangular bodies of model. The conditions for determining the exact form of gravity sources by the proposed approach were clarified. An analysis of the inversion of gravity acceleration anomalies from three bodies for sections, which consisted of three, four and five layers, was performed. The bodies of identical density were determined exactly for three and four layers, and with errors for the five layers. The model with two densities had some errors in determining the geometry of the bodies. Scientific novelty is related to the development of new algorithm, which is based on the principle of movingfor density between blocks. The new approach fulfils gravity inversion much faster than the genetic algorithm. The practical significance is determination of the number of layers needed for a proper inversion of gravity anomalies. The main factors influencing the accuracy of geometry are the data of body density and the number of layers. The method was used for the interpretation of gravity data over the coal deposits of the Donbas. The model consisted of five layers and a range of minimum to their maximum densities for the geological section. Data of the inversion was determined by the depth to the surface of the limestone, which is confirmed by the data of the drilling. The developed approach determines the geometry of section from density, which is important for geological, search and environmental tasks.


Geophysics ◽  
1983 ◽  
Vol 48 (6) ◽  
pp. 713-721 ◽  
Author(s):  
B. J. Last ◽  
K. Kubik

We present a new criterion for the inversion of gravity data. The principle employed is to minimize the volume of the causative body, which is equivalent to maximizing its compactness. The anomalous density distribution is obtained using an iterative technique which is numerically stable and rapidly convergent. The principle can also be adapted to include modeling of gravity anomalies by single‐density sources. The advantage of this approach is that desirable geologic characteristics are automatically incorporated into the model with a minimum of subjective judgments on the part of the interpreter. The treatment of noise in the data fits naturally into the formulation of the inversion procedure. The method is illustrated by the inversion of noise‐free and noisy data generated from a two‐dimensional model consisting of a regular array of identical rectangular blocks whose densities can be individually specified. In every case the algorithm successfully recovers the correct density distribution from the data. In the case of noise‐contaminated data, a complete separation of the noise from the signal is achieved. The practical effectiveness of the method is demonstrated by the inversion of published gravity data. The results obtained are compared with existing models and with available drilling information.


2019 ◽  
Vol 2019 ◽  
pp. 1-15 ◽  
Author(s):  
Sévérin Nguiya ◽  
Willy Lemotio ◽  
Philippe Njandjock Nouck ◽  
Marcelin M. Pemi ◽  
Alain-Pierre K. Tokam ◽  
...  

The structure of the transition zone between the north-western boundary of the Congo Craton and the Kribi-Campo sedimentary basin is still a matter of scientific debate. In this study, the existing gravity data are interpreted in order to better understand the geodynamics of the area. Qualitatively, results show that the major gravity highs are associated with long-wavelength shallow sources of the coastal sedimentary basin, while large negative anomalies trending E-W correlate to low dense intrusive bodies found along the northern limit of the Congo Craton. For the delineation of the causative sources, the gravity anomalies have been inverted based on the Parker-Oldenburg iterative process. As inputs, we used a reference depth of 20 km obtained by spectral analysis and successively, the density contrasts 0.19 g/cm3 and 0.24 g/cm3, deduced from available 1D shear wave velocity models. The results reveal an irregular topography of the mafic interface characterized by a sequence of horst and graben structures with mafic depths varying between 15.6 km and 23.4 km. The shallower depths (15.6-17 km) are associated with the uprising of the mafic interface towards the upper crust. This intrusion may have been initiated during the extension of the Archean Ntem crust resulting in a thinning of the continental crust beneath the coastal sedimentary basin. The subsidence of the mafic interface beneath the craton is materialized by 2 similar graben structures located beneath both Matomb and Ebolowa at a maximum depth of 23.4 km. The intermediate depths (18-22 km) are correlated to the suture zone along the Pouma-Bipindi area. The location of some landslides across the area matches within the northern margin of the Congo Craton and suggests that this margin may also impact on their occurrence. This work provides new insights into the geodynamics, regional tectonics, and basin geometry.


2018 ◽  
Vol 80 (2) ◽  
Author(s):  
Nurul Fairuz Diyana Bahrudin ◽  
Umar Hamzah

Major fault structures of Peninsular Malaysia were interpreted by satellite gravity data obtained from EGM2008. Filtering including THD, TVD, TDR, TDX and Euler Deconvolution inversion techniuques were applied to the data and successfully delineated the major faults especially located in the area separating the granite and sedimentary rock such as Bok Bak, Kuala Lumpur, Bukit Tinggi and Lebir faults. The main finding of this research is the boundary separating the western and eastern belt of Peninsular Malaysia namely the Bentong-Raub Suture by the abrupt changes of gravity anomalies between the two belts. The average depths and dips of Kuala Lumpur, Bukit Tinggi and Seremban faults were estimated by gravity inversion 3D Euler Deconvolution. 


2021 ◽  
Author(s):  
Igor Ognev ◽  
Jörg Ebbing ◽  
Peter Haas

Abstract. Volgo-Uralia is a Neoarchean easternmost part of the East European craton. Recent seismic studies of the Volgo-Uralian region provided new insights into the crustal structure of this area. In this study, we combine satellite gravity and seismic data in a common workflow to perform a complex study of Volgo-Uralian crustal structure which is useful for further basin analysis of the area. In this light, a new crustal model of the Volgo-Uralian subcraton is presented from a step-wise approach: (1) inverse gravity modeling followed by (2) 3D forward gravity modeling. First, inversion of satellite gravity gradient data was applied to determine the Moho depth for the area. Density contrasts between crust and mantle were varied laterally according to the tectonic units present in the region, and the model is constrained by the available active seismic data. The Moho discontinuity obtained from the gravity inversion was consequently modified and complemented in order to define a complete 3D crustal model by adding information on the sedimentary cover, upper crust, lower crust, and lithospheric mantle layers in the process of forward gravity modeling where both seismic and gravity constraints were respected. The obtained model shows crustal thickness variations from 32 to more than 55 km in certain areas. The thinnest crust with a thickness below 40 km is found beneath the Pericaspian basin, which is covered by a thick sedimentary layer. The thickest crust is located underneath the Ural Mountains as well as in the center of the Volga-Uralian subcraton. In both areas the crustal thickness exceeds 50 km. At the same time, initial forward gravity modeling has shown a gravity misfit of ca. 95 mGal between the measured Bouguer gravity anomaly and the forward calculated gravity field in the central area of the Volga-Uralian subcraton. This misfit was interpreted and modeled as a high-density lower crust which possibly represents underplated material. Our preferred crustal model of the Volga-Uralian subcraton respects the gravity and seismic constraints and reflects the main geological features of the region with Moho thickening in the cratons and under the Ural Mountains and thinning along the Paleoproterozoic rifts, Pericaspian sedimentary basin, and Pre-Urals foredeep.


Sign in / Sign up

Export Citation Format

Share Document