The Southeast Indian Ridge between 88°E and 118°E: Gravity anomalies and crustal accretion at intermediate spreading rates

Although slow spreading ridges characterized by a deep axial valley and fast spreading ridges characterized by an axial bathymetric high have been extensively studied, the transition between these two modes of axial morphology is not well understood. We conducted a geophysical-survey of the intermediate spreading rate Southeast Indian Ridge between 88 degrees E and 118 degrees E, a 2300-km-long section of the ridge located between the Amsterdam hot spot and the Australian-Antarctic Discordance where satellite gravity data suggest that the Southeast Indian Ridge (SEIR) undergoes a change from an axial high in the west to an axial valley in the east. A basic change in axial morphology is found near 103 degrees 30'E in the shipboard data; the axis to the west is marked by an axial high, while a valley is found to the east. Although a well-developed axial high, characteristic of the East Pacific Rise (EPR), is occasionally present, the more common observation is a rifted high that is lower and pervasively faulted, sometimes with significant (> 50 m throw) faults within a kilometer of the axis. A shallow axial valley (< 700 m deep) is observed from 104 degrees E to 114 degrees E with a sudden change to a deep (>1200 m deep) valley across a transform at 114 degrees E. The changes in axial morphology along the SEIR are accompanied by a 500 m increase in near-axis ridge flank depth from 2800 m near 88 degrees E to 3300 m near 114 degrees E and by a 50 mGal increase in the regional level of mantle Bouguer gravity anomalies over the same distance, The regional changes in depth and mantle Bouguer anomaly (MBA) gravity can be both explained by a 1.7-2.4 km change in crustal thickness or by a mantle temperature change of 50 degrees C-90 degrees C. In reality, melt supply (crustal thickness) and mantle temperature are linked, so that changes in both may occur simultaneously and these estimates serve as upper bounds. The along-axis MBA gradient is not uniform. Pronounced steps in the regional level of the MBA gravity occur at 103 degrees 30'E-104 degrees E and at 114 degrees E-116 degrees E and correspond to the changes in the nature of the axial morphology and in the amplitude of abyssal hill morphology suggesting that the different forms of morphology do not grade into each other but rather represent distinctly different forms of axial (s)tructure and tectonics with a sharp transition between them. The change from an axial high to an axial valley requires a threshold effect in which the strength of the lithosphere changes quickly. The presence or absence of a quasi-steady state magma chamber may provide such a mechanism. The different forms of axial morphology are also associated with different intrasegment MBA gravity patterns. Segments with an axial high have an MBA low located at a depth minimum near the center of the segment, At EPR-like segments, the MBA low is about 10 mGal with along-axis gradients of 0.15-0.25 mGal/km, similar to those observed at the EPR, Rifted highs have a shallower low and lower gradients suggesting an attenuated composite magma chamber and a reduced and perhaps episodic melt supply. Segments with a shallow axial valley have very flat along-axis MBA profiles with little correspondence between axial depth and axial MBA gravity.

Gravity Modelling of the Crustal Structure Below Bavali Shear Zone And Adjacent Major Plutons Of Northern Kerala

The primary aim of the present study is to acquire a large amount of gravity data, to prepare gravity maps and interpret the data in terms of crustal structure below the Bavali shear zone and adjacent regions of northern Kerala. The gravity modeling is basically a tool to obtain knowledge of the subsurface extension of the exposed geological units and their structural relationship with the surroundings. The study is expected to throw light on the nature of the shear zone, crustal configuration below the high-grade granulite terrain and the tectonics operating during geological times in the region. The Bavali shear is manifested in the gravity profiles by a steep gravity gradient. The gravity models indicate that the Bavali shear coincides with steep plane that separates two contrasting crustal densities extending beyond a depth of 30 km possibly down to Moho, justifying it to be a Mantle fault. It is difficult to construct a generalized model of crustal evolution in terms of its varied manifestations using only the gravity data. However, the data constrains several aspects of crustal evolution and provides insights into some of the major events.

Integrated marine gravity field in the Brazilian coast from altimeter-derived sea surface gradient and shipborne gravity

Sea surface gradients derived from the Geosat and ERS-1 satellite altimetry geodetic missions were integrated with marine gravity data from the National Geophysical Data Center and Brazilian national surveys. Using the least squares collocation method, models of free-air gravity anomaly and geoid height were calculated for the coast of Brazil with a resolution of 2` x 2`. The integration of satellite and shipborne data showed better statistical results in regions near the coast than using satellite data only, suggesting an improvement when compared to the state-of-the-art global gravity models. Furthermore, these results were obtained with considerably less input information than was used by those reference models. The least squares collocation presented a very low content of high-frequency noise in the predicted gravity anomalies. This may be considered essential to improve the high resolution representation of the gravity field in regions of ocean-continent transition. (C) 2010 Elsevier Ltd. All rights reserved.

Investigation of the Goias Alkaline Province, Central Brazil: Application of gravity and magnetic methods

We investigate the strong magnetic and gravity anomalies of the Goias Alkaline Province (GAP), a region of Late Cretaceous alkaline magmatism along the northern border of the Parana Basin, Brazil. The alkaline complexes (eight of which are present in outcrops, two others inferred from magnetic signals) are characterized by a series of small intrusions forming almost circular magnetic and gravimetric anomalies varying from -4000 to +6000 nT and from -10 to +40 mGal, respectively. We used the Aneuler method and Analytical Signal Amplitude to obtain depth and geometry for mapped sources from the magnetic anomaly data. These results were used as the reference models in the 3D gravity inversion. The 3D inversion results show that the alkaline intrusions have depths of 10-12 km. The intrusions in the northern GAP follow two alignments and have different sizes. In the anomaly magnetic map, dominant guidelines correlate strongly with the extensional regimes that correlate with the rise of alkaline magmatism. The emplacement of these intrusions marks mechanical discontinuities and zones of weakness in the upper crust. According to the 3D inversion results, those intrusions are located within the upper crust (from the surface to 18 km depth) and have spheres as the preferable geometry. Such spherical shapes are more consistent with magmatic chambers instead of plug intrusions. The Registro do Araguaia anomaly (similar to 15 by 25 km) has a particular magnetic signature that indicates that the top is deeper than 1500 m. North of this circular anomaly are lineaments with structural indices indicating contacts on their edges and dikes/sills in the interiors. Results of 3D inversion of magnetic and gravity data suggest that the Registro do Araguaia is the largest body in the area, reaching 18 km depth and indicating a circular layered structure. (C) 2011 Elsevier Ltd. All rights reserved.

Gravity mapping in the southern Weddell Sea region, Antarctica (Scale 1:2 000 000)

New maps of free-air and the Bouguer gravity anomalies on the Weddell Sea sector (70-81° S, 6-75° W) of Antarctica are presented. These maps are based on the first computer compilation of available gravity data collected by ''Sevmorgeologia'' in 1976-89 in the southern Weddell Sea and adjacent coasts of western Dronning Maud Land (WDML) and Coats Land. The accomplished gravity studies comprise airborne observations with a line spacing of about 20 km and conventional measurements at over-the-ice points, which were spaced at 10-30 km and supplemented by seismic soundings. Hence, anomalies on the maps represent mainly large-scale and deep crustal features. The dominant feature in free-air gravity map is a large dipolar gravity anomaly stretching along the continental margin. Following the major grain of seabed morphology this shelf-edge/slope anomaly (SESA) is clearly divided into three segments characterized by diverse anomaly amplitudes, wavelengths and trends. They are associated with continental margins of different geotectonic provinces of Antarctica surrounding the Weddell Sea. Apparent distinctions in the SESA signatures are interpreted as the gravity expression of tectonic, deep crustal structure segmentation of the continental margin. The prominent gravity highs (100-140 mGal) of the shelf edge anomaly mapped along WDML are assumed to represent high-density mantle injections intruded into the middle/lower crust during initial rifting of continental breakup. Enlarged wavelengths and diminished amplitudes of the gravity anomaly westwards, along the Weddell Sea embayment (WSE) margin, reflect a widening of the continental slope and a significant increase in thickness of underlying sediment strata. Low amplitude, negative free-air anomalies in the Filchner-Ronne Ice Shelves (FRIS) contrast sharply with the dominating positive anomalies offshore. This indicates a greater sedimentary thickness of the basin in this area. Crustal response to the enlarged sediment load is impressed in mostly positive features of the Bouguer gravity field observed here. Two pronounced positive Bouguer anomalies of 50-70 mGal and an average widths of 200 km dominate the Weddell Sea embayment margins towards the Antarctic Peninsula and the East Antarctic craton. They correlate well with very deep seabed troughs (> 1000 m below sea level). The gravity highs are most likely caused by a shallow upper mantle underneath graben-rift structures evolved at the margins of the WSE basin. A regional zone (> 100 km in width) of the prominent Bouguer and free-air negative anomalies (-40 to -60 mGal) adjacent Coats Land to the north of the ice shelf edge may indicate the presence of the thick old cratonic crust far offshore beneath the Weddell Sea Embayment.

The distribution and ages of regional lithologies in the lunar maria

A research program was designed (1) to map regional lithological units of the lunar surface based on measurements of spatial variations in spectral reflectance, and, (2) to establish the sequence of the formation of such lithological units from measurements of the accumulated affects of impacting bodies.

Spectral reflectance data were obtained by scanning luminance variations over the lunar surface at three wavelengths (0.4µ, 0.52µ, and 0.7µ). These luminance measurements were reduced to normalized spectral reflectance values relative to a standard area in More Serenitotis. The spectral type of each lunar area was identified from the shape of its reflectance spectrum. From these data lithological units or regions of constant color were identified. The maria fall into two major spectral classes: circular moria like More Serenitotis contain S-type or red material and thin, irregular, expansive maria like Mare Tranquillitatis contain T-type or blue material. Four distinct subtypes of S-type reflectances and two of T-type reflectances exist. As these six subtypes occur in a number of lunar regions, it is concluded that they represent specific types of material rather than some homologous set of a few end members.

The relative ages or sequence of formation of these more units were established from measurements of the accumulated impacts which have occurred since more formation. A model was developed which relates the integrated flux of particles which hove impacted a surface to the distribution of craters as functions of size and shape. Erosion of craters is caused chiefly by small bodies which produce negligible individual changes in crater shape. Hence the shape of a crater can be used to estimate the total number of small impacts that have occurred since the crater was formed. Relative ages of a surface can then be obtained from measurements of the slopes of the walls of the oldest craters formed on the surface. The results show that different maria and regions within them were emplaced at different times. An approximate absolute time scale was derived from Apollo 11 crystallization ages under an assumption of a constant rote of impacting for the last 4 x 10^9 yrs. Assuming, constant flux, the period of mare formation lasted from over 4 x 10^9 yrs to about 1.5 x 10^9 yrs ago.

A synthesis of the results of relative age measurements and of spectral reflectance mapping shows that (1) the formation of the lunar maria occurred in three stages; material of only one spectral type was deposited in each stage, (2) two distinct kinds of maria exist, each type distinguished by morphology, structure, gravity anomalies, time of formation, and spectral reflectance type, and (3) individual maria have complicated histories; they contain a variety of lithic units emplaced at different times.

Evolution of the tectogene concept, 1930-1965

The tectogene, or crustal downbuckle, was proposed in the early 1930s by F.A. Vening Meinesz to explain the unexpected belts of negative gravity anomalies in island arcs. He attributed the isostatic imbalance to a deep sialic root resulting from the action of subcrustal convection currents. Vening Meinesz's model was initially corroborated experimentally by P.H. Kuenen, but additional experiments by D.T. Griggs and geological analysis by H.H. Hess in the late 1930s led to substantial revision in detail. As modified, the tectogene provided a plausible model for the evolution of island arcs into alpine mountain belts for another two decades. Additional revisions became necessary in the early 1950s to accommodate the unexpected absence of sialic crust in the Caribbean and the marginal seas of the western Pacific. By 1960 the cherished analogy between island arcs and alpine mountain belts had collapsed under the weight of the detailed field investigations by Hess and his students in the Caribbean region. Hess then incorporated a highly modified form of the tectogene into his sea-floor spreading hypothesis. Ironically, this final incarnation of the concept preserved some of the weaker aspects of the 1930s original, such as the ad hoc explanation for the regular geometry of island arcs.

均衡重力异常计算与综合解释——以黄海地区为例

Petroleum and Natural Gas is an important strategic resources. The reserves of Petroleum and Natural Gas can’t meet the need of our country, which also blocks the development of economy and threatens the safety of national. Therefore, it makes a great sense to bring “the second round of oil & gas exploration” into effect and study the exploration of oil and gas of Pre-Cenozoic residual basins in China. The integrated geophysical exploration is the main way to research the Pre-Cenozoic residual basins. Gravity exploration is one of the most important exploration methods, which has played an important role in oil and gas prospecting, such as compartmentalizing geotectonic elements, delineating the distribution range of sedimentary basins, searching oil and gas structure, abstracting oil and gas information, and so on, from its naissance. The isostatic gravity anomalies is significant for exploration, which can help us research deep crustal structure, the equilibrium state of earth, the geologic structure of shallow crust, the basement shape of sedimentary basins and the genetic evolution of sedimentary basins. In the paper, we stress the implication and physical meanings systemically, and discuss the calculation theory. On the basis of previous work, we test different isostatic compensation models and parameters to find out their influences to the result of isostatic gravity anomalies. In addition, we improve the method of isostatic gravity anomalies calculation and give a system of isostatic gravity anomalies calculation which is proved has satisfying effect. From the research above, we find that the results of Platt model and Airy model are consistent, which have similar form and almost the same value. However, by contrast, the Airy model is proved has better adaptability than Platt model. The two main parameters——crust thickness and density difference of crust and mantle, both have influence to the isostatic gravity anomalies, but the latter have more. Finally, we adopt the regional field extending edge method to make the result more of actual geologic condition. On the methods above, we calculate the isostatic gravity anomalies field in Yellow Sea area from the Bouguer gravity anomalies and the water depth and altitude data. And then the isostatic gravity anomalies character is analyzed and the integrated geological-geophysical interpretation is made on the basis of summarizing the previous research result systemically and analyzing other geophysical data and geological information. From the research, we find that the Yellow Sea area belongs to continental type crust equilibrium regions, where the isostatic gravity anomalies field is placid and has less fluctuation values, which implies that the area is in equilibrium state to different extends.

Geophysical Constraints on Structure and Tectonics of the Eastern Arabian Sea and the Adjoining West Coast of India with Special Reference to the Kerala Basin

The main objective of the present study is to model the gravity fields in terms of lithospheric structure below the western continental margin of India (WCMI) identify zones of crustal mass anomalies and attempt to infer the location of Ocean Continent transition in the Arabian Sea. In this study, the area starting from the western shield margin to the region covering the deep oceanic parts of the Arabian Sea which is bounded by Carlsberg and Cerg and Central Indian ridges in the south, eastern part of the Indus Cone in the west and falling between 630E and 800E longitudes, and 50N - 240N latitudes has been considered. The vast amount of seismic reflection and refraction data in the form of crustal velocities, basement configuration and crustal thicknesses available for the west coast as well as the eastern Arabian Sea has been utilized for this purpose

Relics of eclogite facies assemblages in the Ceara Central Domain, NW Borborema Province, NE Brazil: Implications for the assembly of West Gondwana

The Borborema Province, in the NE of Brazil, is a rather complex piece in the Brazil-Africa puzzle as it represents the junction of the Dahomeyide/Pharusian, Central African, Aracuai and Brasilia fold belts located between the West-African/Sao Luis, Congo/Sao Francisco and Amazonas craton. The correlation between the Dahomeyides from W-Africa (Ghana, Benin, Togo, and Mali) and the Borborema Province involves the Medio Coreau and Central Ceara domains. The inferred continuation of the main oceanic suture zone exposed in the Dahomeyides of W Africa is buried beneath the Phanerozoic Parnaiba Basin in Brazil (northwest of the Medio Coreau domain) where some high density gravity anomalies may represent hidden remnants of an oceanic suture. In addition to this major suture a narrow, nearly continuous strip composed of mainly mafic pods containing relics of eclogite-facies assemblages associated with partially migmatized granulite-facies metapelitic gneisses has been found further east in the NW Borborema Province. These high pressure mafic rocks, interpreted as retrograded eclogites, are located between the Transbrasiliano Lineament and the Santa Quiteria continental arc and comprise primitive to evolved arc-related rocks with either arc- or MORB-type imprints that can indicate either deep subduction of oceanic lithosphere or roots of continental and oceanic magmatic arcs. Average peak P-T conditions under eclogite-facies metamorphism (T=770 degrees C and P = 17.3 kbar) were estimated using garnet-clinopyroxene thermometry and Jd content in clinopyroxene. Transition to granulite-facies conditions, as well as later widespread re-equilibration under amphibolite facies, were registered both in the basic and the metapelitic rocks and suggest a clockwise P-T path characterized by an increase in temperature followed by strong decompression. A phenomenon possibly related to the exhumation of a highly thickened crust associated with the suturing of the Medio Coreau and Central Ceara domains, two distinct crustal blocks separated by the Transbrasiliano Lineament. (C) 2009 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Relics of eclogite facies assemblages in the Ceara Central Domain, NW Borborema Province, NE Brazil: Implications for the assembly of West Gondwana

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Uma nova abordagem para interpretação de anomalias gravimétricas regionais e residuais aplicada ao estudo da organização crustal: exemplo da Região Norte do Piauí e Noroeste do Ceará

A interpretação de anomalias gravimétricas é de grande importância no estudo de feições geológicas que ocorrem na crosta terrestre. Esta interpretação é, no entanto, dificultada pelo fato das anomalias gravimétricas serem resultantes da soma total dos efeitos produzidos por todos os contrastes de densidades de subsuperfície. Desse modo, com o intuito de separar efeitos de feições mais profundas de efeitos de feições mais rasas, bem como a caracterização da geometria desses dois conjuntos de feições, apresentamos um método de separação das componentes regional e residual do campo e a subsequente interpretação de cada componente. A separação regional-residual de dados gravimétricos é efetuada através da aproximação do campo regional por um polinômio ajustado ao campo observado por um método robusto. Este método é iterativo e usa como aproximação inicial a solução obtida através do ajuste polinomial pelo método dos mínimos quadrados. O método empregado minimiza a influência de observações contendo forte contribuição do campo residual no ajuste do campo regional. A componente regional obtida a partir da separação regional-residual é transformada em um mapa de distâncias verticais em relação a um nível de referência. Esta transformação compreende duas etapas. A primeira consiste na obtenção da continuação para baixo da componente regional, que é pressuposta ser causada por uma interface suave separando dois meios homogêneos, representando a interface crosta-manto, cujo contraste de densidade é supostamente conhecido. A segunda consiste na transformação do mapa de continuação para baixo em um mapa de distâncias verticais entre o nível de continuação (tomado como nível de referência) e a interface. Este método apresenta duas dificuldades. A primeira está ligada à instabilidade, havendo portanto a necessidade do emprego de um estabilizador o que acarreta a perda de resolução das feições que se desejam mapear. A segunda, inerente ao método gravimétrico, consiste na impossibilidade da determinação das profundidades absolutas da interface em cada ponto, bastando entretanto o conhecimento da profundidade absoluta em um ponto, através de informação independente, para que todas as outras profundidades absolutas sejam conhecidas. A componente residual obtida a partir da separação regional-residual é transformada em um mapa de contrastes de densidade aparente. Esta transformação consiste no cálculo do contraste de densidade de várias fontes prismáticas através de uma inversão linear pressupondo que as fontes reais estejam das a uma placa horizontal, com contrastes de densidade variando apenas nas direções horizontais. O desempenho do método de separação regional-residual apresentado foi avaliado, através de testes empregando dados sintéticos, fornecendo resultados superiores em relação aos métodos dos mínimos quadrados e da análise espectral. O método de interpretação da componente regional teve seu desempenho avaliado em testes com dados sintéticos onde foram produzidos mapeamentos de interfaces bem próximas das estruturas reais. O limite de resolução das feições que se desejam mapear depende não só do grau do polinômio ajustante, como também da própria limitação inerente ao método gravimétrico. Na interpretação da componente residual é necessário que se postule ou tenha informação a priori sobre a profundidade do topo e espessura da placa onde as fontes estão supostamente confinadas. No entanto, a aplicação do método em dados sintéticos, produziu estimativas razoáveis para os limites laterais das fontes, mesmo na presença de fontes interferentes, e pressupondo-se valores para profundidade do topo e espessura da placa, diferentes dos valores verdadeiros. A ambiguidade envolvendo profundidade do topo, espessura e densidade pode ser visualizada através de gráficos de valores de densidade aparente contra profundidade do topo presumida para a placa para vários valores postulados para a espessura da placa. Estes mesmos gráficos permitem, pelo aspecto das curvas, a elaboração de uma interpretação semi-quantitativa das profundidades das fontes reais. A seqüência dos três métodos desenvolvidos neste trabalho foi aplicada a dados gravimétricos da região norte do Piauí e noroeste do Ceará levando a um modelo de organização crustal que compreende espessamentos e adelgaçamentos crustais associados a um evento compressivo que possibilitou a colocação de rochas densas da base da crosta a profundidades rasas. Este modelo ê compatível com os dados geológicos de superfície. É ainda sugerida a continuidade, por mais 200 km em direção a sudoeste, do Cinturão de Cisalhamento Noroeste do Ceará por sob os sedimentos da Bacia do Parnaíba, com base nas evidências fornecidas pela interpretação da anomalia residual. Embora esta seqüência de métodos tenha sido desenvolvida com vistas ao estudo de feições crustais de porte continental, ela também pode ser aplicada ao estudo de feições mais localizadas como por exemplo no mapeamento do relevo do embasamento de/bacias sedimentares onde os sedimentos são cortados por rochas intrusivas mais densas.

Aplicação do gravímetro de poço em problemas de avaliação de formação no Brasil

No período de outubro a novembro de 1990, foi realizado na Bacia Sergipe-Alagoas o primeiro levantamento gravimétrico de poço no Brasil, através de um programa com a participação da Universidade Federal do Pará, a PETROBRÁS e o U.S.G.S. (United States Geological Survey). Este levantamento teve como objetivos o teste do equipamento, a resposta dos perfis gravimétricos em relação aos problemas encontrados nas bacias sedimentares brasileiras e a comparação com os dados de densidade obtidos com o perfil de densidade compensada (CDL). Os levantamentos foram realizados em três poços e os dados obtidos passaram por um processo de redução, onde os valores foram transformados para miligals e corrigidos dos efeitos de maré, deriva e terreno, para obter as densidades preliminares. Porém, a região onde foram realizados os levantamentos apresenta além do embasamento muito raso, uma seqüência evaporítica com contraste de densidade grande em relação ao resto do pacote sedimentar, gerando gradientes verticais anômalos. Estes efeitos foram então corrigidos e obtidos os valores finais de densidade. Com base em uma caracterização litológica prévia, os valores de densidade foram então comparados com o perfil CDL. Verificou-se a boa qualidade dos dados, o que permitiu o cálculo de valores de porosidade, mostrando assim possíveis intervalos de interesse à exploração de hidrocarbonetos, tanto no embasamento, reservatório mais importante, quanto no intervalo sedimentar.

Interpretação quantitativa de anomalias gravimétricas ar-livre da margem continental norte brasileira, setor nordeste

Anomalias gravimétricas ar-livre de perfis perpendiculares a margem continental do tipo passiva apresentam uma configuração padrão. Esta configuração é, satisfatoriamente, explicada por um modelo geofísico formado por uma distribuição de descontinuidades horizontais bidimensionais. Um processo automático de busca aleatória é proposto para a interpretação quantitativa dos dados. Através do método de poliedros flexíves (Simplex), os parâmetros principais do modelo - o contraste de densidade, a profundidade, o rejeito e a localização de cada descontinuidade, puderam ser encontrados, admitindo uma relação número de pontos/número de parâmetros, a determinar, conveniente. Sobre a região do talude, as anomalias ar-livre da margem continental podem ser explicadas por uma única descontinuidade horizontal (degrau simples); e tendo que a resposta dos dados gravimétricos no domínio do número de onda contém informações sobre esta anomalia, foi proposto um procedimento gráfico iterativo para a análise espectral deste sinal. Aplicando a transformada de Fourier é possível determinar a profundidade e o rejeito da descontinuidade, e conhecendo estes parâmetros a densidade é calculada unicamente. O objetivo básico do uso destes procedimentos seria combinar os dois métodos de interpretação nos domínios do espaço e do número de onda, com a finalidade de obter soluções vinculadas mais plausíveis quanto ao contexto geológico esperado para a área estudada. Os dois procedimentos de interpretação foram aplicados nas anomalias gravimétricas ar-livre da margem continental norte brasileira, setor nordeste, abrangendo os estados do Maranhão ao Rio Grande do Norte. As respectivas capacidade de resolução de cada procedimento foram então analisadas. Demonstrou-se que a inversão realizada diretamente no domínio do espaço é mais favorável na interpretação das anomalias ar-livre, embora o tratamento espectral seja relativamente mais simples.

Inversão gravimétrica estável do relevo do embasamento e da variação da densidade com a profundidade em bacias sedimentares

O presente método postula uma variação hiperbólica para o contraste de densidade de uma bacia sedimentar em função da profundidade, e tem dois objetivos: (1) delinear o relevo do embasamento de uma bacia, conhecendo-se a anomalia gravimétrica, o contraste de densidade na superfície da bacia e o fator de decaimento do contraste de densidade com a profundidade; (2) estimar, além do relevo, o valor do contraste de densidade na superfície de uma bacia sedimentar e o fator de decaimento do contraste de densidade com a profundidade, sendo fornecida a anomalia gravimétrica e a profundidade do embasamento em alguns pontos da bacia. Nos dois casos o modelo interpretativo é um conjunto de prismas retangulares verticais justapostos, cujas espessuras, que são parâmetros a serem estimados, representam a profundidade da interface de separação entre os sedimentos e o embasamento. As soluções obtidas nos dois problemas acima formulados são estáveis devido à incorporação de informações adicionais sobre a suavidade do relevo estimado, e o conhecimento da profundidade do relevo do embasamento em alguns pontos, fornecido por furos de sondagem. O método foi testado em anomalias gravimétricas sintéticas produzidas pela simulação de bacias sedimentares com relevos suaves. Os resultados mostraram relevos com boa resolução e valores estimados do contraste de densidade na superfície da bacia e do fator de decaimento do contraste de densidade com a profundidade, próximos aos verdadeiros, indicando dessa maneira o potencial do método em interpretações gravimétricas de bacias sedimentares. O método foi aplicado à anomalia Bouguer da Bacia do Recôncavo, Brasil, delineando um relevo com um valor para a profundidade máxima de cerca de 6 km, semelhante ao estimado em interpretações sísmicas. As estimativas para o contraste de densidade na superfície da Bacia e o fator de decaimento com a profundidade foram -0,30 g/cm³ e 30 km, respectivamente, produzindo uma estimativa para a compactação máxima dos sedimentos da ordem de 4%.