Magnetotelluric deep soundings, gravity and geoid in the south Sao Francisco craton: Geophysical indicators of cratonic lithosphere rejuvenation and crustal underplating

In the south Sao Francisco craton a circular and 8-m amplitude geoid anomaly coincides with the outcropping terrain of an Archean-Paleoproterozoic basement. Broadband magnetotelluric (MT) data inversions of two radial profiles within the positive geoid and Bouguer gravity anomaly yield geo-electrical crustal sections, whereby the lower crust is locally more conductive (10 to 100 Omega m) in spatial coincidence with a denser lower crust modeled by the gravity data. This anomalous lower crust may have resulted from magmatic underplating, associated with Mesoarchean and Proterozoic episodes of tholeiitic dike intrusion. Long-period MT soundings reveal a low electrical resistivity mantle (20 to 200 Omega m) from depths beyond 120 km. Forward geoid modeling, using the scope of the low electrical resistivity region within the mantle as a constraint, entails a density increase (40 to 50 kg/m(3)) possibly due to Fe enrichment of mantle minerals. However, this factor alone does not explain the observed resistivity. A supplemented presence of small amounts of percolated carbonatite melting (similar to 0.005 vol.%), dissolved water and enhanced oxygen fugacity within the peridotitic mantle are viable agents that could explain the less resistive upper mantle. We propose that metasomatic processes confined in the sub-continental lithospheric mantle foster the conditions for a low degree melting with variable CO(2), H(2)O and Fe content. Even though the precise age of this metasomatism is unknown it might be older than the Early Cretaceous based on the evidence that a high-degree of melting in a lithospheric mantle impregnated with carbonatites originated the tholeiitic dike intrusions dispersed from the southeastern border of the Sao Francisco craton, during the onset of the lithosphere extension and break-up of the western Gondwana. The proxies are the NE Parana and Espinhaco (130 Ma, Ar/Ar ages) tholeiitic dikes, which contain (similar to 3%) carbonatites in their composition. The occurrence of a positive geoid anomaly (+ 10 m) and pre-tholeiites (age > 138 Ma), carbonatites and kimberlites along the west African continental margin (Angola and Namibia) reinforces the presumed age of the Sao Francisco-Congo craton rejuvenation to be prior to its fragmentation in the Lower Cretaceous. (C) 2010 Elsevier B.V. All rights reserved.

Deep crustal structure of the Parana Basin from receiver functions and Rayleigh-wave dispersion: Evidence for a fragmented cratonic root

The deep crustal structure of the Parana Basin of southern Brazil is investigated by analyzing P- and PP-wave receiver functions at 17 Brazilian Lithosphere Seismic Project stations within the basin. The study area can be described as a typical Paleozoic intracratonic basin that hosts one of the largest Large Igneous Province of the world and makes a unique setting for investigating models of basin subsidence and their interaction with mantle plumes. Our study consists of (1) an analysis of the Moho interaction phases in the receiver functions to obtain the thickness and bulk Vp/Vs ratio of the basin`s underlying crust and (2) a joint inversion with Rayleigh-wave dispersion velocities from an independent tomographic study to delineate the detailed S-wave velocity variation with depth. The results of our analysis reveal that Moho depths and bulk Vp/Vs ratios (including sediments) vary between 41 and 48 km and between 1.70 and 1.76, respectively, with the largest values roughly coinciding with the basin`s axis, and that S-wave velocities in the lower crust are generally below 3.8 km/s. Select sites within the basin, however, show lower crustal S-wave velocities slightly above 3.9 km/s suggestive of underplated mafic material. We show that these observations are consistent with a fragmented cratonic root under the Parana basin that defined a zone of weakness for the initial Paleozoic subsidence of the basin and which allowed localized mafic underplating of the crust along the suture zones by Cenozoic magmatism.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

High T/P evolution and metamorphic ages of the migmatitic basement of northern Sierras Pampeanas, Argentina: Characterization of a mid-crustal segment of the Famatinian belt

New petrologic, thermobarometric and U-Pb monazite geochronologic information allowed to resolve the metamorphic evolution of a high temperature mid-crustal segment of an ancient subduction-related orogen. The EI Portezuelo Metamorphic-Igneous Complex, in the northern Sierras Pampeanas, is mainly composed of migmatites that evolved from amphibolite to granulite metamorphic facies, reaching thermal peak conditions of 670-820 degrees C and 4.5-5.3 kbar. The petrographic study combined with conventional and pseudosection thermobarometry led to deducing a short prograde metamorphic evolution within migmatite blocks. The garnet-absent migmatites represent amphibolite-facies rocks, whereas the cordierite-garnet-K-feldspar-sillimanite migmatites represent higher metamorphic grade rocks. U-Pb geochronology on monazite grains within leucosome record the time of migmatization between approximate to 477 and 470 Ma. Thus, the El Portezuelo Metamorphic-Igneous Complex is an example of exhumed Early Ordovician anatectic middle crust of the Famatinian mobile belt. Homogeneous exposure of similar paleo-depths throughout the Famatinian back-arc and isobaric cooling paths suggest slow exhumation and consequent longstanding crustal residence at high temperatures. High thermal gradients uniformly distributed in the Famatinian back-arc can be explained by shallow convection of a low-viscosity asthenosphere promoted by subducting-slab dehydration. (C) 2011 Elsevier Ltd. All rights reserved.

Petrogenesis of Early Cretaceous silicic volcanism in SE Uruguay: The role of mantle and crustal sources

Early Cretaceous (similar to 129 Ma) silicic rocks crop out in SE Uruguay between the Laguna Merin and Santa Lucia basins in the Lascano, Sierra Sao Miguel. Salamanca and Minas areas They are mostly rhyolites with minor quartz-trachytes and are nearly contemporaneous with the Parana-Etendeka igneous province and with the first stages of South Atlantic Ocean opening A strong geochemical variability (particularly evident from Rb/Nb, Nb/Y trace element ratios) and a wide range of Sr-Nd isotopic ratios ((143)Nd/(144)Nd((129)) = 0.51178-0.51209, (87)Sr/(86)Sr((129)) = 0.70840-0.72417) characterize these rocks Geochemistry allows to distiniguish two compositional groups, corresponding to the north-eastern (Lascano and Sierra Sao Miguel, emplaced on the Neo-Proterozoic southern sector of the Dom Feliciano mobile belt) and south-eastern localities (Salamanca, Minas, emplace on the much older (Archean) Nico Perez teriane or on the boundary between the Dom Feliciano and Nico Perez termites) These compositional differences between the two groups are explained by variable mantle source and crust contributions. The origin of the silicic magmas is best explained by complex processes involving assimilation and fractional crystallization and mixing of a basaltic magma with upper crustal lithologies, for Lascano and Sierra Sao Miguel rhyolites. In the Salamanea and Minas rocks genesis, a stronger contribution from lower crust is indicated.

Sr-Nd isotopic evidence for crustal contamination in the Niquelandia complex, Goias, Central Brazil

The Niquelandia complex is a Neoproterozoic mafic-ultramafic intrusion resulting from fractional crystallization of primary picritic basalt intrusions. It consists of two layered sequences: a lower and larger one (LS), where four stratigraphic units exhibit an upward decrease of ultramafic layers and increase of gabbroic layers; an upper, smaller sequence (US), separated from LS by a high-temperature shear zone and consisting of two stratigraphic units (gabbros + anorthosites and amphibolites). Nd and Sr isotopic analyses and rare earth element (REE) profiles provide evidence that the complex suffered important crustal contamination. The LS isotopic array trends from a DM region with positive epsilon Nd and moderately positive epsilon Sr towards a field occupied by crustal xenoliths, especially abundant in the upper LS (negative epsilon Nd and large, positive E:Sr). Each LS stratigraphic unit is distinct from the next underlying unit, showing lower epsilon Nd and higher epsilon Sr, suggesting inputs of fresh magma and mixing with the contaminated, residual magma. The US is characterised by a relatively high variation of epsilon Nd and constant epsilon Sr. REE patterns vary within each unit from LREE depleted to LREE enriched in the samples having lower epsilon Nd and higher epsilon Sr. The contamination process has been modelled by using the EC-AFC algorithms from [Spera, F.J., Bohrson, W.A., 2001. Energy-constrained open-system magmatic processes 1: general model and energy-constrained assimilation and fractional crystallization (EC-AFC) formulation. J. Petrology 42, 999-1018]. The differences between the LS and US isotopic arrays are consistent with contamination by the same crustal component, provided that its melting degree was higher in LS than in US. The different degrees of anatexis are explained by the heat budget released from the magma, higher in LS (because of its larger mass) than in US. Comparison of the correlations between isotopes and incompatible trace element ratios of the models and of the gabbros shows some differences, which are demonstrably related with the variable amount of cumulus phases and trapped melt in the gabbros. (c) 2007 Elsevier Ltd. All rights reserved.

Cordierite-bearing lavas from Jaguarao, southern Brazil: Petrological evidence for crustal melts during early rifting of Gondwana

The Jaguarao stratoid dacites (Rio Grande do Sul, Brazil) are limited in areal extent, are comprised of about 3.2 km(3) of preserved erupted material, and outcrop only in areas of the region underlain by mylonitic and ultramylonitic rocks. They are S-type volcanic rocks containing cordierite, orthopyroxene, plagioclase, and ilmenite as liquidus phases, and partially melted granite, gneiss, and migmatite enclaves that are very similar to the Precambrian basement rocks. The Jaguarao lavas have distinct geochemical signatures and Sr-Nd isotopes with respect to other volcanic rocks of the region. Available geochronological data for Jaguarao dacites range between 157 +/- 5 Ma and 139.6 +/- 7.4 Ma. Considering the errors, the younger ages obtained for Jaguarao lavas overlap the 138-128 Ma age of rocks of the Serra Geral Group, and thus indicate that the dacites were erupted prior to the break-up of Gondwana in this region. Petrographic, mineralogical, and petrochemical data, as well as the tectonic context of the Jaguarao lavas, suggest that magma genesis was linked, at least in part, to friction melts. The dacitic magma was generated by partial melting reactions involving biotite breakdown in a dominantly quartz-feldspathic source terrane, leaving a granulite facies residue in subsurface. These melts were probably generated as a consequence of crustal thinning linked to simple shear extension just prior to Gondwana break-up and rifting of the southern Atlantic Ocean. (C) 2009 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

(40)Ar/(39)Ar ages and Sr-Nd-Pb-Os geochemistry of CAMP tholeiites from Western Maranhao basin (NE Brazil)

The Central Atlantic Magmatic Province (CAMP), emplaced at the Triassic-Jurassic (T-J) boundary (-200 Ma), is among the largest igneous provinces on Earth. The Maranhao basin in NE Brazil is located around 700 km inland and 2000 km from the site of the earliest Pangea disruption. The CAMP tholeiites occur only in the western part of the basin and have been described as low and high-Ti. Here we document the occurrence of two sub-groups among the high-Ti tholeiites in the Western Maranhao basin. The major and trace elements and the Sr-Nd-Pb isotopic ratios define three chemical groups corresponding to the low-Ti (TiO(2)< 1.3 wt.%), high-Ti (TiO(2)-2.0 wt.%) and evolved high-Ti (TiO(2 >)3 wt.%) western Maranhao basin tholeiites (WMBT). The new (40)Ar/(39)Ar plateau ages obtained on plagioclase separates for high-Ti (199.7 +/- 2.4 Ma) and evolved high-Ti WMBT (197.2 +/- 0.5 Ma and 198.2 +/- 0.6 Ma) are indistinguishable and identical to those of previously analyzed low-Ti WMBT (198.5 +/- 0.8 Ma) and to the mean (40)Ar/(39)Ar age of the CAMP (199 +/- 2.4 Ma). We also present the first Re-Os isotopic data for CAMP basalts. The low and high-Ti samples display mantle-like initial ((187)Os/(188)Os)(i) ranging from 0.1267 to 0.1299, while the evolved high-Ti samples are more radiogenic (((187)Os/ (188)Os)(i) up to 0.184) We propose that the high-Ti WMBT were derived from the sub-lithospheric asthenosphere, and contaminated during ascent by interaction with the subcontinental lithospheric mantle (SCLM). The evolved high-Ti WMBT were derived from the same asthenospheric source but experienced crustal contamination. The chemical characteristics of the low-Ti group can be explained by partial melting of the most fertile portions of the SCLM metasomatized during paleo-subduction. Alternatively, the low-Ti WMBT could be derived from the sub-lithospheric asthenosphere but the resulting melts may have undergone contamination by the SCLM. The occurrences of high-Ti basalts are apparently not restricted to the area of initial continental disruption which may bring into question previous interpretations such as those relating high-Ti CAMP magmatism to the initiation of Atlantic ridge spreading or as the expression of a deep mantle plume. We propose that the CAMP magmatism in the Maranhao basin may be attributed to local hotter mantle conditions due to the combined effects of edge-driven convection and large-scale mantle warming under the Pangea supercontinent. The involvement of a mantle-plume with asthenosphere-like isotopic characteristics cannot be ruled out either as one of the main source components of the WMBT or as a heat supplier. (C) 2010 Elsevier BM. All rights reserved.

Moho map of South America from receiver functions and surface waves

We estimate crustal structure and thickness of South America north of roughly 40 degrees S. To this end, we analyzed receiver functions from 20 relatively new temporary broadband seismic stations deployed across eastern Brazil. In the analysis we include teleseismic and some regional events, particularly for stations that recorded few suitable earthquakes. We first estimate crustal thickness and average Poisson`s ratio using two different stacking methods. We then combine the new crustal constraints with results from previous receiver function studies. To interpolate the crustal thickness between the station locations, we jointly invert these Moho point constraints, Rayleigh wave group velocities, and regional S and Rayleigh waveforms for a continuous map of Moho depth. The new tomographic Moho map suggests that Moho depth and Moho relief vary slightly with age within the Precambrian crust. Whether or not a positive correlation between crustal thickness and geologic age is derived from the pre-interpolation point constraints depends strongly on the selected subset of receiver functions. This implies that using only pre-interpolation point constraints (receiver functions) inadequately samples the spatial variation in geologic age. The new Moho map also reveals an anomalously deep Moho beneath the oldest core of the Amazonian Craton.

Geoelectric evidence for centripetal resurge of impact melt and breccias over central uplift of Araguainha impact structure

We present five profiles from electrical resistivity tomography (ERT), with surface constraints and gravity data, in the central uplift of the Araguainha impact structure in central Brazil. The central uplift, the overlying polymict breccias, and decameter-scale impact melt rocks are characterized by contrasting ranges of electrical resistivity. Our resistivity model provides empirical evidence that supports the existing model in which impact melt and breccias resurged toward the crater center in the final stages of the cratering process. On the basis of our results from the first use of ERT in impact cratering studies, we conclude that the deposition and flow of impact melt and breccias over the central uplift were influenced by the geometry of the lithologic boundaries in the central uplift.

Induced seismicity in the Castarthao reservoir, NE Brazil - Preliminary results

The Castanhao reservoir was built in the state of Ceara, a dry region in Northeastern Brazil, to regulate the flow of the Jaguaribe River, for irrigation, and for power generation. It is an earth-filled dam, 60 m high, with a water capacity of 4.5 x 10(9) m(3). The seismicity in the area has been monitored since 1998, with a few interruptions, using one analog or one digital station and, during a few periods, a three-station network. The first earthquakes likely to be induced events were detected in 2003, when the water level was about 20 in high. In early 2004 a very heavy rainfall season quickly filled the reservoir. Shortly after, an increase in the seismic activity occurred and many micro-earthquakes were recorded. We suggest that this activity resulted from an increase in pore pressure due to undrained response. Therefore, we may classify this cluster of microearthquakes as ""initial seismicity."" We deployed a network with four analog stations in the area, following this activity, to determine the epicentral zone. At least three epicentral areas under the reservoir were detected. The spatio-temporal analysis of the available data revealed that the seismicity occurs in clusters and that these were activated at different periods. We identified four sets of faults (N-S-, E-W-, NW-SE-, and NE-SW-oriented), some of which moved in shallow crustal levels and as recently as the Quaternary (1.8 Ma). Under the present-day stress regime, the last two sets moved as strike-slip structures. We suggest a possible correlation between dormant faults and the observed induced seismicity. (c) 2008 Elsevier B.V. All rights reserved.

Basement depths in the Parecis Basin (Amazon) with receiver functions from small local earthquakes in the Porto dos Gauchos seismic zone

Receiver functions from small local earthquakes were used to determine sediment thicknesses in Porto dos Gauchos seismic zone (PGSZ), Parecis basin, Amazonian craton, Brazil. The high velocity contrast between basement and sediments (P-wave velocities of 6.1 and 3.0 km/s, respectively) favors the generation of clear P-to-S converted phases (Ps) seen in the radial component, and also S-to-P conversions (Sp) seen in the vertical component. A reference 10 velocity model determined with shallow refraction experiment in PGSZ helped to convert Ps P time differences to basement depths at 15 stations deployed for aftershocks studies. The results of receiver function integrated with the shallow refraction reveal that the basement depths in the PGSZ increases from the basin border in the north up to about 600 m depth in the south. The basement topography, however, does not vary smoothly and a basement high with a steep topography was detected near the epicentral area. A 400 m elevation difference within 1.7 km distance suggests a possible border fault of a buried graben. This feature seems to be oriented roughly WSW-ENE and could indicate basement structures related to the seismicity of the Porto dos Gauchos Seismic Zone. (C) 2011 Elsevier Ltd. All rights reserved.

Programs to compute magnetization to density ratio and the magnetization inclination from 3-D gravity and magnetic anomalies

Vector field formulation based on the Poisson theorem allows an automatic determination of rock physical properties (magnetization to density ratio-MDR-and the magnetization inclination-MI) from combined processing of gravity and magnetic geophysical data. The basic assumptions (i.e., Poisson conditions) are: that gravity and magnetic fields share common sources, and that these sources have a uniform magnetization direction and MDR. In addition, the previously existing formulation was restricted to profile data, and assumed sufficiently elongated (2-D) sources. For sources that violate Poisson conditions or have a 3-D geometry, the apparent values of MDR and MI that are generated in this way have an unclear relationship to the actual properties in the subsurface. We present Fortran programs that estimate MDR and MI values for 3-D sources through processing of gridded gravity and magnetic data. Tests with simple geophysical models indicate that magnetization polarity can be successfully recovered by MDR-MI processing, even in cases where juxtaposed bodies cannot be clearly distinguished on the basis of anomaly data. These results may be useful in crustal studies, especially in mapping magnetization polarity from marine-based gravity and magnetic data. (c) 2007 Elsevier Ltd. All rights reserved.

Neoproterozoic geodynamic evolution of SW-Gondwana: a southern African perspective

Our current understanding of the tectonic history of the principal Pan-African orogenic belts in southwestern Africa, reaching from the West Congo Belt in the north to the Lufilian/Zambezi, Kaoko, Damara, Gariep and finally the Saldania Belt in the south, is briefly summarized. On that basis, possible links with tectono-stratigraphic units and major structures on the eastern side of the Rio de la Plata Craton are suggested, and a revised geodynamic model for the amalgamation of SW-Gondwana is proposed. The Rio de la Plata and Kalahari Cratons are considered to have become juxtaposed already by the end of the Mesoproterozoic. Early Neoproterozoic rifting led to the fragmentation of the northwestern (in today`s coordinates) Kalahari Craton and the splitting off of several small cratonic blocks. The largest of these ex-Kalahari cratonic fragments is probably the Angola Block. Smaller fragments include the Luis Alves and Curitiba microplates in eastern Brazil, several basement inliers within the Damara Belt, and an elongate fragment off the western margin, named Arachania. The main suture between the Kalahari and the Congo-So Francisco Cratons is suspected to be hidden beneath younger cover between the West Congo Belt and the Lufilian/Zambezi Belts and probably continues westwards via the Cabo Frio Terrane into the Goias magmatic arc along the Brasilia Belt. Many of the rift grabens that separated the various former Kalahari cratonic fragments did not evolve into oceanic basins, such as the Northern Nosib Rift in the Damara Belt and the Gariep rift basin. Following latest Cryogenian/early Ediacaran closure of the Brazilides Ocean between the Rio de la Plata Craton and the westernmost fragment of the Kalahari Craton, the latter, Arachania, became the locus of a more than 1,000-km-long continental magmatic arc, the Cuchilla Dionisio-Pelotas Arc. A correspondingly long back-arc basin (Marmora Basin) on the eastern flank of that arc is recognized, remnants of which are found in the Marmora Terrane-the largest accumulation of oceanic crustal material known from any of the Pan-African orogenic belts in the region. Corresponding foredeep deposits that emerged from the late Ediacaran closure of this back-arc basin are well preserved in the southern areas, i.e. the Punta del Este Terrane, the Marmora Terrane and the Tygerberg Terrane. Further to the north, present erosion levels correspond with much deeper crustal sections and comparable deposits are not preserved anymore. Closure of the Brazilides Ocean, and in consequence of the Marmora back-arc basin, resulted from a change in the Rio de la Plata plate motion when the Iapetus Ocean opened between the latter and Laurentia towards the end of the Ediacaran. Later break-up of Gondwana and opening of the modern South Atlantic would have followed largely along the axis of the Marmora back-arc basin and not along major continental sutures.

Major shear zones of southern Brazil and Uruguay: escape tectonics in the eastern border of Rio de La plata and Paranapanema cratons during the Western Gondwana amalgamation

The Mantiqueira Province represents a series of supracrustal segments of the South-American counterpart formed during the Gondwana Supercontinent agglutination. In this crustal domain, the process of escape tectonics played a conspicuous role, generating important NE-N-S-trending lineaments. The oblique component of the motions of the colliding tectonic blocks defined the transpressional character of the main suture zones: Lancinha-Itariri, Cubato-Arcadia-Areal, Serrinha-Rio Palmital in the Ribeira Belt and Sierra Ballena-Major Gercino in the Dom Feliciano Belt. The process as a whole lasted for ca. 60 Ma, since the initial collision phase until the lateral escape phase predominantly marked by dextral and subordinate sinistral transpressional shear zones. In the Dom Feliciano Belt, southern Brazil and Uruguay, transpressional event at 630-600 Ma is recognized and in the Ribeira Belt, despite less coevally, the transpressional event occurred between 590 and 560 Ma in its northern-central portion and between ca. 625 and 595 Ma in its central-southern portion. The kinematics of several shear zones with simultaneous movement in opposite directions at their terminations is explained by the sinuosity of these lineaments in relation to a predominantly continuous westward compression.