Upper-mantle seismic structure beneath SE and Central Brazil from P- and S-wave regional traveltime tomography

We present models for the upper-mantle velocity structure beneath SE and Central Brazil using independent tomographic inversions of P- and S-wave relative arrival-time residuals (including core phases) from teleseismic earthquakes. The events were recorded by a total of 92 stations deployed through different projects, institutions and time periods during the years 1992-2004. Our results show correlations with the main tectonic structures and reveal new anomalies not yet observed in previous works. All interpretations are based on robust anomalies, which appear in the different inversions for P-and S-waves. The resolution is variable through our study volume and has been analyzed through different theoretical test inversions. High-velocity anomalies are observed in the western portion of the Sao Francisco Craton, supporting the hypothesis that this Craton was part of a major Neoproterozoic plate (San Franciscan Plate). Low-velocity anomalies beneath the Tocantins Province (mainly fold belts between the Amazon and Sao Francisco Cratons) are interpreted as due to lithospheric thinning, which is consistent with the good correlation between intraplate seismicity and low-velocity anomalies in this region. Our results show that the basement of the Parana Basin is formed by several blocks, separated by suture zones, according to model of Milani & Ramos. The slab of the Nazca Plate can be observed as a high-velocity anomaly beneath the Parana Basin, between the depths of 700 and 1200 km. Further, we confirm the low-velocity anomaly in the NE area of the Parana Basin which has been interpreted by VanDecar et al. as a fossil conduct of the Tristan da Cunha Plume related to the Parana flood basalt eruptions during the opening of the South Atlantic.

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.

Upper-mantle seismic anisotropy from SKS splitting in the South American stable platform: A test of asthenospheric flow models beneath the lithosphere

Upper-mantle seismic anisotropy has been extensively used to infer both present and past deformation processes at lithospheric and asthenospheric depths. Analysis of shear-wave splitting (mainly from core-refracted SKS phases) provides information regarding upper-mantle anisotropy. We present average measurements of fast-polarization directions at 21 new sites in poorly sampled regions of intra-plate South America, such as northern and northeastern Brazil. Despite sparse data coverage for the South American stable platform, consistent orientations are observed over hundreds of kilometers. Over most of the continent, the fast-polarization direction tends to be close to the absolute plate motion direction given by the hotspot reference model HS3-NUVEL-1A. A previous global comparison of the SKS fast-polarization directions with flow models of the upper mantle showed relatively poor correlation on the continents, which was interpreted as evidence for a large contribution of ""frozen"" anisotropy in the lithosphere. For the South American plate, our data indicate that one of the reasons for the poor correlation may have been the relatively coarse model of lithospheric thicknesses. We suggest that improved models of upper-mantle flow that are based on more detailed lithospheric thicknesses in South America may help to explain most of the observed anisotropy patterns.

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.

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.

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.

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.

Age and geochemistry of mantle peridotites and diorite dykes from the Baldissero body: Insights into the Paleozoic-Mesozoic evolution of the Southern Alps

Trace element and isotopic data obtained for mantle spinel Iherzolites and diorite dykes from the Baldissero massif (Ivrea-Verbano Zone, Western Italy) provide new, valuable constraints on the petrologic and geodynamic evolution of the Southern Alps in Paleozoic to Mesozoic times. Whole rock and mineral chemistry indicates that Baldissero Iherzolites can be regarded as refractory mantle residues following limited melt extraction. In particular, the Light Rare Earth Elements (LREE)-depleted and fractionated compositions of whole rock and clinopyroxene closely match modelling results for refractory residues after low degrees (similar to 4-5%) of near-fractional melting of depleted mantle, possibly under garnet-facies conditions. Following this, the peridotite sequence experienced subsolidus re-equilibration at lithospheric spinel-facies conditions and intrusion of several generations of dykes. However, Iherzolites far from dykes show very modest metasomatic changes, as evidenced by the crystallisation of accessory titanian pargasite and the occurrence of very slight enrichments in highly incompatible trace elements (e.g. Nb). The Re-Os data for Iherzolites far from the dykes yield a 376 Ma (Upper Devonian) model age that is considered to record a partial melting event related to the Variscan orogenic cycle s.l. Dioritic dykes cutting the mantle sequence have whole rock, clinopyroxene and plagioclase characterised by high radiogenic Nd and low radiogenic Sr, which point to a depleted to slightly enriched mantle source. Whole rock and mafic phases of diorites have high Mg# values that positively correlate with the incompatible trace element concentrations. The peridotite at the dyke contact is enriched in orthopyroxene, iron and incompatible trace elements with respect to the Iherzolites far from dykes. Numerical simulations indicate that the geochemical characteristics of the diorites can be explained by flow of a hydrous, silica-saturated melt accompanied by reaction with the ambient peridotite and fractional crystallisation. The composition of the more primitive melts calculated in equilibrium with the diorite minerals show tholeiitic to transitional affinity. Internal Sm-Nd, three-point isochrons obtained for two dykes suggest an Upper Triassic-Lower Jurassic emplacement age (from 204 31 to 198 29 Ma). Mesozoic igneous events are unknown in the southern Ivrea-Verbano Zone (IVZ), but the intrusion of hydrous melts, mostly silica-saturated, have been well documented in the Finero region, i.e. the northernmost part of IVZ and Triassic magmatism with calc-alkaline to shoshonitic affinity is abundant throughout the Central-Eastern Alps. The geochemical and chronological features of the Baldissero diorites shed new light on the geodynamic evolution of the Southern Alps before the opening of the Jurassic Tethys. (C) 2010 Elsevier B.V. All rights reserved.