Arcabouço geofísico, isostasia e causas do magmatismo cenozóico da província Borborema e de sua margem continental (Nordeste do Brasil)

The Borborema Province (BP) is a geologic domain located in Northeastern Brazil. The BP is limited at the south by the São Francisco craton, at the west by the Parnaíba basin, and both at the north and east by coastal sedimentary basins. Nonetheless the BP surface geology is well known, several key aspects of its evolution are still open, notably: i)its tectonic compartmentalization established after the Brasiliano orogenesis, ii) the architecture of its cretaceous continental margin, iii) the elastic properties of its lithosphere, and iv) the causes of magmatism and uplifting which occurred in the Cenozoic. In this thesis, a regional coverage of geophysical data (elevation, gravity, magnetic, geoid height, and surface wave global tomography) were integrated with surface geologic information aiming to attain a better understanding of the above questions. In the Riacho do Pontal belt and in the western sector of the Sergipano belt, the neoproterozoic suture of the collision of the Sul domain of the BP with the Sanfranciscana plate (SFP) is correlated with an expressive dipolar gravity anomaly. The positive lobule of this anomaly is due to the BP lower continental crust uplifting whilst the negative lobule is due to the supracrustal nappes overthrusting the SFP. In the eastern sector of the Sergipano belt, this dipolar gravity anomaly does not exist. However the suture still can be identified at the southern sector of the Marancó complex arc, alongside of the Porto da Folha shear zone, where the SFP N-S geophysical alignments are truncated. The boundary associated to the collision of the Ceará domain of the BP with the West African craton is also correlated with a dipolar gravity anomaly. The positive lobule of this anomaly coincides with the Sobral-Pedro II shear zone whilst the negative lobule is associated with the Santa Quitéria magmatic arc. Judging by their geophysical signatures, the major BP internal boundaries are: i)the western sector of the Pernambuco shear zone and the eastern continuation of this shear zone as the Congo shear zone, ii) the Patos shear zone, and iii) the Jaguaribe shear zone and its southwestern continuation as the Tatajuba shear zone. These boundaries divide the BP in five tectonic domains in the geophysical criteria: Sul, Transversal, Rio Grande do Norte, Ceará, and Médio Coreaú. The Sul domain is characterized by geophysical signatures associated with the BP and SFP collision. The fact that Congo shear zone is now proposed as part of the Transversal domain boundary implies an important change in the original definition of this domain. The Rio Grande do Norte domain presents a highly magnetized crust resulted from the superposition of precambrian and phanerozoic events. The Ceará domain is divided by the Senador Pompeu shear zone in two subdomains: the eastern one corresponds to the Orós-Jaguaribe belt and the western one to the Ceará-Central subdomain. The latter subdomain exhibits a positive ENE-W SW gravity anomaly which was associated to a crustal discontinuity. This discontinuity would have acted as a rampart against to the N-S Brasiliano orogenic nappes. The Médio Coreaú domain also presents a dipolar gravity anomaly. Its positive lobule is due to granulitic rocks whereas the negative one is caused by supracrustal rocks. The boundary between Médio Coreaú and Ceará domains can be traced below the Parnaíba basin sediments by its geophysical signature. The joint analysis of free air anomalies, free air admittances, and effective elastic thickness estimates (Te) revealed that the Brazilian East and Equatorial continental margins have quite different elastic properties. In the first one 10 km < Te < 20 km whereas in the second one Te ≤ 10 km. The weakness of the Equatorial margin lithosphere was caused by the cenozoic magmatism. The BP continental margin presents segmentations; some of them have inheritance from precambrian structures and domains. The segmentations conform markedly with some sedimentary basin features which are below described from south to north. The limit between Sergipe and Alagoas subbasins coincides with the suture between BP and SFP. Te estimates indicates concordantly that in Sergipe subbasin Te is around 20 km while Alagoas subbasin has Te around 10 km, thus revealing that the lithosphere in the Sergipe subbasin has a greater rigidity than the lithosphere in the Alagoas subbasin. Additionally inside the crust beneath Sergipe subbasin occurs a very dense body (underplating or crustal heritage?) which is not present in the crust beneath Alagoas subbasin. The continental margin of the Pernambuco basin (15 < Te < 25 km) presents a very distinct free air edge effect displaying two anomalies. This fact indicates the existence in the Pernambuco plateau of a relatively thick crust. In the Paraíba basin the free air edge effect is quite uniform, Te ≈ 15 km, and the lower crust is abnormally dense probably due to its alteration by a magmatic underplating in the Cenozoic. The Potiguar basin segmentation in three parts was corroborated by the Te estimates: in the Potiguar rift Te ≅ 5 km, in the Aracati platform Te ≅ 25 km, and in the Touros platform Te ≅ 10 km. The observed weakness of the lithosphere in the Potiguar rift segment is due to the high heat flux while the relatively high strength of the lithosphere in the Touros platform may be due to the existence of an archaean crust. The Ceará basin, in the region of Mundaú and Icaraí subbasins, presents a quite uniform free air edge effect and Te ranges from 10 to 15 km. The analysis of the Bouguer admittance revealed that isostasy in BP can be explained with an isostatic model where combined surface and buried loadings are present. The estimated ratio of the buried loading relative to the surface loading is equal to 15. In addition, the lower crust in BP is abnormally dense. These affirmations are particularly adequate to the northern portion of BP where adherence of the observed data to the isostatic model is quite good. Using the same above described isostatic model to calculate the coherence function, it was obtained that a single Te estimate for the entire BP must be lower than 60 km; in addition, the BP north portion has Te around 20 km. Using the conventional elastic flexural model to isostasy, an inversion of crust thickness was performed. It was identified two regions in BP where the crust is thickened: one below the Borborema plateau (associated to an uplifting in the Cenozoic) and the other one in the Ceará domain beneath the Santa Quitéria magmatic arc (a residue associated to the Brasiliano orogenesis). On the other hand, along the Cariri-Potiguar trend, the crust is thinned due to an aborted rifting in the Cretaceous. Based on the interpretation of free air anomalies, it was inferred the existence of a large magmatism in the oceanic crust surrounding the BP, in contrast with the incipient magmatism in the continent as shown by surface geology. In BP a quite important positive geoid anomaly exists. This anomaly is spatially correlated with the Borborema plateau and the Macaú-Queimadas volcanic lineament. The integrated interpretation of geoid height anomaly data, global shear velocity model, and geologic data allow to propose that and Edge Driven Convection (EDC) may have caused the Cenozoic magmatism. The EDC is an instability that presumably occurs at the boundary between thick stable lithosphere and oceanic thin lithosphere. In the BP lithosphere, the EDC mechanism would have dragged the cold lithospheric mantle into the hot asthenospheric mantle thus causing a positive density contrast that would have generated the main component of the geoid height anomaly. In addition, the compatibility of the gravity data with the isostatic model, where combined surface and buried loadings are present, together with the temporal correlation between the Cenozoic magmatism and the Borborema plateau uplifting allow to propose that this uplifting would have been caused by the buoyancy effect of a crustal root generated by a magmatic underplating in the Cenozoic

Avaliação de dados de altimetria do satélite Geosat a partir de comparação com dados do levantamento gravimétrico marinho Equant

O objetivo principal deste trabalho é avaliar o uso dos dados de altimetria de satélites para mapear a superfície do potencial gravitacional (geóide) no mar. Esta avaliação se faz por comparações da resolução e precisão entre os dados de altimetria processados numa superfície equipotencial (o mar) e dados obtidos a partir de levantamentos convencionais. Uma vez processada a superfície equipotencial, quantidades tais como a anomalia "ar livre" juntamente com o desvio vertical podem ser calculadas. Os dados altimétricos ("altura do mar") utilizados neste trabalho foram coletados pelo satélite GEOSAT. Este satélite rastreou diversas áreas oceânicas do globo processando 44 ciclos em dois anos. Alguns pesquisadores utilizaram os valores médios da "altura do mar" deste satélite para melhoramentos em precisão e resolução dos registros. Estes valores tratados estão disponíveis em NOAA (National Oceanic and Atmospheric Administration) sendo deste modo repassados à UFPa para utilização nesta tese. Os dados de gravimetria marinha utilizados neste trabalho são aqueles obtidos do levantamento "Equatorial Atlantic" (EQUANT I e II) resultantes de uma pesquisa conjunta entre várias instituições com objetivos científicos de conhecer o comportamento da margem equatorial Brasileira. Para comparação e integração entre os dois tipos de dados obtidos através de fontes distintas (medidas de satélite e do navio), poder-se-ia obter a aceleração vertical numa superfície equipotencial partindo-se de um tratamento algébrico dos dados coletados por rastreamento altimétrico do satélite Geosat ou alternativamente poder-se-ia processar transformações dos dados de gravimetria marinha em uma superfície equipotencial equivalente. Em decorrência de diferenças no espaçamento entre as linhas dos levantamentos, ou seja, as linhas das trajetórias do satélite estão largamente espaçadas em comparação com aquelas do levantamento marinho, optou-se por transformar os dados gravimétricos do navio em uma superfície equipotencial. Neste tipo de transformação, consideraram-se vários fatores tais como efeitos "aliasing", nível de ruídos nos levantamentos do navio, redução ao geóide (correção "ar livre"), bem como erros computacionais durante as transformações. Com a supressão parcial desses efeitos (enfatizando o "aliasing") encontrou-se forte correlação entre os dois conjuntos de dados verificando-se um nível de coerência satisfatório a partir do comprimento de onda de 11 km. Comparando-se este resultado com o nível de resolução do satélite Geosat largamente estudado por outros pesquisadores enfatiza-se que de fato a resolução dos valores médios (2 anos) do satélite Geosat aproxima-se da resolução dos levantamentos do Equant I e Equant II.

Modelos geoidais híbridos para o Estado de São Paulo: desenvolvimento e avaliação

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Validation of the MAPGEO2010 and comparison with recent geopotential models

This paper presents the offorts to calculate the geoid model for Brazil. It is limited by 6 degrees N and 35 degrees S in latitude and 30 degrees W and 75 degrees W in longitude. The terrestrial gravity data for the continent have been updated by means of the most recent surveys in Brazil and in the neighbour countries. The complete Bouguer and Helmert gravity anomalies have been derived through the Canadian package SHGEO. The short wavelength component was estimated via FFT. The geopotential model EGM2008 was used as a reference field restricted to degree and order 150. The model was validated over 844 GPS observations on Bench Marks of the spirit leveling network. The height anomalies plus a topography dependent correction term derived from EGM2008 (degree 2190 and order 2159), GO_CONS_GCF_2_DIR_R2 (degree and order 240), GOCO02S (degree and order 250), EIGEN 51C (degree and order 359) and EIGEN 6C (degree and order 1420), geoidal height derived from MAPGEO2004 (old official geoid model in Brazil) have also been compared to the GPS points on Bench Marks.

Sea Level Variations from Decades to Millennia: A few Case Studies

The main goals of this Ph.D. study are to investigate the regional and global geophysical components related to present polar ice melting and to provide independent cross validation checks of GIA models using both geophysical data detected by satellite mission, and geological observations from far field sites, in order to determine a lower and upper bound of uncertainty of GIA effect. The subject of this Thesis is the sea level change from decades to millennia scale. Within ice2sea collaboration, we developed a Fortran numerical code to analyze the local short-term sea level change and vertical deformation resulting from the loss of ice mass. This method is used to investigate polar regions: Greenland and Antarctica. We have used mass balance based on ICESat data for Greenland ice sheet and a plausible mass balance for Antarctic ice sheet. We have determined the regional and global fingerprint of sea level variations, vertical deformations of the solid surface of the Earth and variations of shape of the geoid for each ice source mentioned above. The coastal areas are affected by the long wavelength component of GIA process. Hence understanding the response of the Earth to loading is crucial in various contexts. Based on the hypothesis that Earth mantle materials obey to a linear rheology, and that the physical parameters of this rheology can be only characterized by their depth dependence, we investigate the Glacial Isostatic Effect upon the far field sites of Mediterranean area using an improved SELEN program. We presented new and revised observations for archaeological fish tanks located along the Tyrrhenian and Adriatic coast of Italy and new RSL for the SE Tunisia. Spatial and temporal variations of the Holocene sea levels studied in central Italy and Tunisia, provided important constraints on the melting history of the major ice sheets.

Comparison of GOCE-GPS gravity fields derived by different approaches

Several techniques have been proposed to exploit GNSS-derived kinematic orbit information for the determination of long-wavelength gravity field features. These methods include the (i) celestial mechanics approach, (ii) short-arc approach, (iii) point-wise acceleration approach, (iv) averaged acceleration approach, and (v) energy balance approach. Although there is a general consensus that—except for energy balance—these methods theoretically provide equivalent results, real data gravity field solutions from kinematic orbit analysis have never been evaluated against each other within a consistent data processing environment. This contribution strives to close this gap. Target consistency criteria for our study are the input data sets, period of investigation, spherical harmonic resolution, a priori gravity field information, etc. We compare GOCE gravity field estimates based on the aforementioned approaches as computed at the Graz University of Technology, the University of Bern, the University of Stuttgart/Austrian Academy of Sciences, and by RHEA Systems for the European Space Agency. The involved research groups complied with most of the consistency criterions. Deviations only occur where technical unfeasibility exists. Performance measures include formal errors, differences with respect to a state-of-the-art GRACE gravity field, (cumulative) geoid height differences, and SLR residuals from precise orbit determination of geodetic satellites. We found that for the approaches (i) to (iv), the cumulative geoid height differences at spherical harmonic degree 100 differ by only ≈10 % ; in the absence of the polar data gap, SLR residuals agree by ≈96 % . From our investigations, we conclude that real data analysis results are in agreement with the theoretical considerations concerning the (relative) performance of the different approaches.

Time Variable Gravity: Contributions of GOCE Satellite Data to Monthly and Bi-monthly GRACE Gravity Estimates

A feasibility study by Pail et al. (Can GOCE help to improve temporal gravity field estimates? In: Ouwehand L (ed) Proceedings of the 4th International GOCE User Workshop, ESA Publication SP-696, 2011b) shows that GOCE (‘Gravity field and steady-state Ocean Circulation Explorer’) satellite gravity gradiometer (SGG) data in combination with GPS derived orbit data (satellite-to-satellite tracking: SST-hl) can be used to stabilize and reduce the striping pattern of a bi-monthly GRACE (‘Gravity Recovery and Climate Experiment’) gravity field estimate. In this study several monthly (and bi-monthly) combinations of GRACE with GOCE SGG and GOCE SST-hl data on the basis of normal equations are investigated. Our aim is to assess the role of the gradients (solely) in the combination and whether already one month of GOCE observations provides sufficient data for having an impact in the combination. The estimation of clean and stable monthly GOCE SGG normal equations at high resolution ( >  d/o 150) is found to be difficult, and the SGG component, solely, does not show significant added value to monthly and bi-monthly GRACE gravity fields. Comparisons of GRACE-only and combined monthly and bi-monthly solutions show that the striping pattern can only be reduced when using both GOCE observation types (SGG, SST-hl), and mainly between d/o 45 and 60.

Improving Surveying Accuracy and Efficiency in Connecticut: An Accuracy Assessment of GEOID03 and GEOID09

Comparing published NAVD 88 Helmert orthometric heights of First-Order bench marks against GPS-determined orthometric heights showed that GEOID03 and GEOID09 perform at their reported accuracy in Connecticut. GPS-determined orthometric heights were determined by subtracting geoid undulations from ellipsoid heights obtained from a network least-squares adjustment of GPS occupations in 2007 and 2008. A total of 73 markers were occupied in these stability classes: 25 class A, 11 class B, 12 class C, 2 class D bench marks, and 23 temporary marks with transferred elevations. Adjusted ellipsoid heights were compared against OPUS as a check. We found that: the GPS-determined orthometric heights of stability class A markers and the transfers are statistically lower than their published values but just barely; stability class B, C and D markers are also statistically lower in a manner consistent with subsidence or settling; GEOID09 does not exhibit a statistically significant residual trend across Connecticut; and GEOID09 out-performed GEOID03. A "correction surface" is not recommended in spite of the geoid models being statistically different than the NAVD 88 heights because the uncertainties involved dominate the discrepancies. Instead, it is recommended that the vertical control network be re-observed.

What Does Height Really Mean? Part III: Height Systems

This is the third paper in a four-part series considering the fundamental question, “what does the word “height” really mean?” The first paper reviewed reference ellipsoids and mean sea level datums. The second paper reviewed the physics of heights culminating in a simple development of the geoid and explained why mean sea level stations are not all at the same orthometric height. This third paper develops the principle notions of height, namely measured, differentially deduced changes in elevation, orthometric heights, Helmert orthometric heights, normal orthometric heights, dynamic heights, and geopotential numbers. We conclude with a more in-depth discussion of current thoughts regarding the geoid.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Sediment volume variations inferred from DGPS measurements within Bossons glacier proglacial area (Mont-Blanc massif, France)

This dataset presents Differential Global Positioning System data (DGPS) acquired within the Bossons glacier proglacial area. Bossons glacier is a rapidly retreating glacier and its proglacial area is deglaciated for ~30 years. Bossons stream is one of the outlets of the subglacial drainage system. It starts as a 800 m steep cascade reach, then flows through an area with gentler slope : the Plan des Eaux (PdE). PdE is a 300 m long, 50 m wide proglacial alluvial plain with an increasing channel mobility in the downstream direction but decreasing slope gradient and incision. As it may act a sediment trap, studying periglacial and proglacial erosion processes in the Bossons catchment requires to quantify PdE sediment volume evolution. A several meter-sized block located within Bossons proglacial area was set up as GPS base : its location was measured by one antenna (Topcon Hyper Pro) by performing 600 consecutive measurements throughout one day. A second antenna (Topcon Hyper Pro) was then used to measure XYZ location of points in the proglacial area with a ~2 m grid. Radio communication between the two antennas allowed differential calculations to be automatically carried out on field using the Topcon FC-250 hand controller. This methodology yields 3 cm XY and 1.5 cm Z uncertainties. DGPS data have been acquired through 10 campaigns from 2004 to 2014; campaigns from 2004 to 2008 cover a smaller area than those from 2010 to 2014. Digital Elevation Model (DEM) have been interpolated from DGPS data and difference between two DEMs yields deposited and eroded volume within PdE. Maps of PdE volume variation between two campaigns show that incision mainly occurs in the upper and lower sections where as deposition dominates in the middle section. Deposition, denudation and net rate (deposition rate - denudation rate) are calculated by normalizing volumes by DEM areas. Deposition dominates results with a mean net rate of 29 mm/yr. However, strong inter-annual variability exists and some years are dominated by denudation : -36 mm/yr and -100 mm/yr for 2006 and 2011, respectively. Nonetheless, oldest campaigns (2004 to 2008) were carried out on the lower part part of the alluvial plain and ruling them out to keep only complete DEM (2010 to 2014) yields a mean net rate of ~15 mm/yr. This results is coherent with field observations of both strong deposition (e.g. flood deposits) and strong erosion (e.g. 30 cm incision) evidences. Bossons glacier proglacial area is thus dynamic with year-to-year geormorphological changes but may leans toward increasing its mean elevation through a deposition dominated system.