Stratification of the earth beneath the Azores from P and S receiver functions

Seismic recordings of IRIS/IDA/GSN station CMLA and of several temporary stations in the Azores archipelago are processed with P and S receiver function (PRF and SRF) techniques. Contrary to regional seismic tomography these methods provide estimates of the absolute velocities and of the Vp/Vs ratio up to a depth of similar to 300 km. Joint inversion of PRFs and SRFs for a few data sets consistently reveals a division of the subsurface medium into four zones with a distinctly different Vp/Vs ratio: the crust similar to 20 km thick with a ratio of similar to 1.9 in the lower crust, the high-Vs mantle lid with a strongly reduced VpNs velocity ratio relative to the standard 1.8, the low-velocity zone (LVZ) with a velocity ratio of similar to 2.0, and the underlying upper-mantle layer with a standard velocity ratio. Our estimates of crustal thickness greatly exceed previous estimates (similar to 10 km). The base of the high-Vs lid (the Gutenberg discontinuity) is at a depth of-SO km. The LVZ with a reduction of S velocity of similar to 15% relative to the standard (IASP91) model is terminated at a depth of similar to 200 km. The average thickness of the mantle transition zone (TZ) is evaluated from the time difference between the S410p and SKS660p, seismic phases that are robustly detected in the S and SKS receiver functions. This thickness is practically similar to the standard IASP91 value of 250 km. and is characteristic of a large region of the North Atlantic outside the Azores plateau. Our data are indicative of a reduction of the S-wave velocity of several percent relative to the standard velocity in a depth interval from 460 to 500 km. This reduction is found in the nearest vicinities of the Azores, in the region sampled by the PRFs, but, as evidenced by SRFs, it is missing at a distance of a few hundred kilometers from the islands. We speculate that this anomaly may correspond to the source of a plume which generated the Azores hotspot. Previously, a low S velocity in this depth range was found with SRF techniques beneath a few other hotspots.

Deep structure of crust and mantle beneath Iberia and Western Mediterranean from P and S Receiver functions and SKS Waveforms

We have conducted a P and S receiver functions [PRFs and SRFs] analysis for 19 seismic stations on the Iberia and western Mediterranean. In the transition zone [TZ] the PRFs analysis reveals a band [from Gibraltar to Balearic] increased by 10-20 km relative to the standard 250 km. The TZ thickness variations are strongly correlated with the P660s times in PRFs. We interpret the variable depth of the 660-km discontinuity as an effect of subduction. Over the anomalous TZ we found a reduced velocity zone in the upper mantle. Joint inversion of PRFs and SRFs reveals a subcrustal high S velocity lid and an underlying LVZ. A reduction of the S velocity in the LVZ is less than 10%. The Gutenberg discontinuity is located at 65±5 km, but in several models sampling the Mediterranean, the lid is missing or its thickness is reduced to ~30 km. In the Gibraltar and North Africa this boundary is located at ~100 km. The lid Vp/Vs beneath Betics is reduced relative to the standard 1.8. Another evidence of the Vp/Vs anomaly is provided by S410p phase late arrivals in the SRFs. The azimuthal anisotropy analysis with a new technology was conducted at 5 stations and at 2 groups of stations. The fast direction in the uppermost mantle layer is ~90º in Iberian Massif. In Balearic is in the azimuth of ~120º. At a depth of ~60 km the direction becomes 90º. Anisotropy in the upper layer can be interpreted as frozen, whereas anisotropy in the lower layer is active, corresponding to the present-day or recent flow. The effect of the asthenosphere in the SKS splitting is much larger than the effect of the subcrustal lithosphere.

Estimation of the Crustal Bulk Properties Beneath Mainland Portugal from P-Wave Teleseismic Receiver Functions

In this work, we present results from teleseismic P-wave receiver functions (PRFs) obtained in Portugal, Western Iberia. A dense seismic station deployment conducted between 2010 and 2012, in the scope of the WILAS project and covering the entire country, allowed the most spatially extensive probing on the bulk crustal seismic properties of Portugal up to date. The application of the H-κ stacking algorithm to the PRFs enabled us to estimate the crustal thickness (H) and the average crustal ratio of the P- and S-waves velocities V p/V s (κ) for the region. Observations of Moho conversions indicate that this interface is relatively smooth with the crustal thickness ranging between 24 and 34 km, with an average of 30 km. The highest V p/V s values are found on the Mesozoic-Cenozoic crust beneath the western and southern coastal domain of Portugal, whereas the lowest values correspond to Palaeozoic crust underlying the remaining part of the subject area. An average V p/V s is found to be 1.72, ranging 1.63-1.86 across the study area, indicating a predominantly felsic composition. Overall, we systematically observe a decrease of V p/V s with increasing crustal thickness. Taken as a whole, our results indicate a clear distinction between the geological zones of the Variscan Iberian Massif in Portugal, the overall shape of the anomalies conditioned by the shape of the Ibero-Armorican Arc, and associated Late Paleozoic suture zones, and the Meso-Cenozoic basin associated with Atlantic rifting stages. Thickened crust (30-34 km) across the studied region may be inherited from continental collision during the Paleozoic Variscan orogeny. An anomalous crustal thinning to around 28 km is observed beneath the central part of the Central Iberian Zone and the eastern part of South Portuguese Zone.

Mantle beneath the Gibraltar Arc from receiver functions

P and S receiver functions (PRF and SRF) from 19 seismograph stations in the Gibraltar Arc and the Iberian Massif reveal new details of the regional deep structure. Within the high-velocity mantle body below southern Spain the 660-km discontinuity is depressed by at least 20 km. The Ps phase from the 410-km discontinuity is missing at most stations in the Gibraltar Arc. A thin (similar to 50 km) low-S-velocity layer atop the 410-km discontinuity is found under the Atlantic margin. At most stations the S410p phase in the SRFs arrives 1.0-2.5 s earlier than predicted by IASP91 model, but, for the propagation paths through the upper mantle below southern Spain, the arrivals of S410p are delayed by up to +1.5 s. The early arrivals can be explained by elevated Vp/Vs ratio in the upper mantle or by a depressed 410-km discontinuity. The positive residuals are indicative of a low (similar to 1.7 versus similar to 1.8 in IASP91) Vp/Vs ratio. Previously, the low ratio was found in depleted lithosphere of Precambrian cratons. From simultaneous inversion of the PRFs and SRFs we recognize two types of the mantle: 'continental' and 'oceanic'. In the 'continental' upper mantle the S-wave velocity in the high-velocity lid is 4.4-4.5 km s(-1), the S-velocity contrast between the lid and the underlying mantle is often near the limit of resolution (0.1 km s(-1)), and the bottom of the lid is at a depth reaching 90 100 km. In the 'oceanic' domain, the S-wave velocities in the lid and the underlying mantle are typically 4.2-4.3 and similar to 4.0 km s(-1), respectively. The bottom of the lid is at a shallow depth (around 50 km), and at some locations the lid is replaced by a low S-wave velocity layer. The narrow S-N-oriented band of earthquakes at depths from 70 to 120 km in the Alboran Sea is in the 'continental' domain, near the boundary between the 'continental' and 'oceanic' domains, and the intermediate seismicity may be an effect of ongoing destruction of the continental lithosphere.

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.

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.

Estimation of the Crustal Bulk Properties Beneath Mainland Portugal from P-Wave Teleseismic Receiver Functions

In this work, we present results from teleseismic P-wave receiver functions (PRFs) obtained in Portugal, Western Iberia. A dense seismic station deployment conducted between 2010 and 2012, in the scope of the WILAS project and covering the entire country, allowed the most spatially extensive probing on the bulk crustal seismic properties of Portugal up to date. The application of the H-kappa stacking algorithm to the PRFs enabled us to estimate the crustal thickness (H) and the average crustal ratio of the P- and S-waves velocities V (p)/V (s) (kappa) for the region. Observations of Moho conversions indicate that this interface is relatively smooth with the crustal thickness ranging between 24 and 34 km, with an average of 30 km. The highest V (p)/V (s) values are found on the Mesozoic-Cenozoic crust beneath the western and southern coastal domain of Portugal, whereas the lowest values correspond to Palaeozoic crust underlying the remaining part of the subject area. An average V (p)/V (s) is found to be 1.72, ranging 1.63-1.86 across the study area, indicating a predominantly felsic composition. Overall, we systematically observe a decrease of V (p)/V (s) with increasing crustal thickness. Taken as a whole, our results indicate a clear distinction between the geological zones of the Variscan Iberian Massif in Portugal, the overall shape of the anomalies conditioned by the shape of the Ibero-Armorican Arc, and associated Late Paleozoic suture zones, and the Meso-Cenozoic basin associated with Atlantic rifting stages. Thickened crust (30-34 km) across the studied region may be inherited from continental collision during the Paleozoic Variscan orogeny. An anomalous crustal thinning to around 28 km is observed beneath the central part of the Central Iberian Zone and the eastern part of South Portuguese Zone.

Cape Verde hotspot from the upper crust to the top of the lower mantle

We investigate the crust, upper mantle and mantle transition zone of the Cape Verde hotspot by using seismic P and S receiver functions from several tens of local seismograph stations. We find a strong discontinuity at a depth of similar to 10 km underlain by a similar to 15-km thick layer with a high (similar to 1.9) Vp/Vs velocity ratio. We interpret this discontinuity and the underlying layer as the fossil Moho, inherited from the pre-hotspot era, and the plume-related magmatic underplate. Our uppermost-mantle models are very different from those previously obtained for this region: our S velocity is much lower and there are no indications of low densities. Contrary to previously published arguments for the standard transition zone thickness our data indicate that this thickness under the Cape Verde islands is up to similar to 30 km less than in the ambient mantle. This reduction is a combined effect of a depression of the 410-km discontinuity and an uplift of the 660-km discontinuity. The uplift is in contrast to laboratory data and some seismic data on a negligible dependence of depth of the 660-km discontinuity on temperature in hotspots. A large negative pressure-temperature slope which is suggested by our data implies that the 660-km discontinuity may resist passage of the plume. Our data reveal beneath the islands a reduction of S velocity of a few percent between 470-km and 510-km depths. The low velocity layer in the upper transition zone under the Cape Verde archipelago is very similar to that previously found under the Azores and a few other hotspots. In the literature there are reports on a regional 520-km discontinuity, the impedance of which is too large to be explained by the known phase transitions. Our observations suggest that the 520-km discontinuity may present the base of the low-velocity layer in the transition zone. (C) 2011 Elsevier B.V. All rights reserved.

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.

Estimativa de espessura crustal na Província Borborema (NE/Brasil) através de função do recpetor

Crustal thickness and VP/VS estimates are essential to the studies of subsurface geological structures and also to the understanding of the regional tectonic evolution of a given area. In this dissertation, we use the Langston´s (1979) Receiver Function Method using teleseismic events reaching the seismographic station with angles close to the vertical. In this method, the information of the geologic structures close to the station is isolated so that effects related to the instrument response and source mechanics are not present. The resulting time series obtained after the deconvolution between horizontal components contains the larger amplitude referring to the P arrival, followed by smaller arrival caused by the reverberation and conversion of the P-wave at the base of the crust. We also used the HK-Stacking after Zhu & Kanamori (2000) to obtain crustal thickness and Vp/VS estimates. This method works stacking receiver functions so that the best estimates of crustal thickness and Vp/VS are found when the direct P, the Ps wave and the first multiple are coherently stacked. We used five broadband seismographic stations distributed over the Borborema Province, NE Brazil. Crustal thickness and Vp/VS estimates are consistent with the crust-mantle interface obtained using gravity data. We also identified crutal thickening in the NW portion of the province, close to Sobral/CE. Towards the center-north portion of the province, there is an evident crustal thinning which coincides with a geological feature consisting of an alignment of sedimentary basins known as the Cariris-Potiguar trend. Towards the NE portion of the province, in Solânea/PB and Agrestina/PE regions, occurs a crustal thickening and a systematic increase in the VP/VS values which suggest the presence of mafic rocks in the lower crust also consistent with the hypothesis of underplating in the region

analise da dispersao de ondas de superficie na provincia borborema, Nordeste do Brasil

The Borborema Province, Northeastern Brazil, had its internal structure investigated by different geophysical methods like gravity, magnetics and seismics. Additionally, many geological studies were also carried out to define the structural domains of this province. Despite the plethora of studies, there are still many important open aspects about its evolution. Here, we study the velocity structure of S-wave in the crust using dispersion of surface waves. The dispersion of surface waves allows an estimate of the average thickness of the crust across the region between the stations. The inversion of the velocity structure was carried out using the inter-station dispersion of surface waves of Rayleigh and Love types. The teleseismic events are mainly from the edges of the South and North American plates. The period of data collection occurred between 2007 and 2010 and we selected 7 events with magnitude above 5.0 MW and up to 40 km depth. The difference between the events back-azimuths and the interstation path was not greater than 10. We also know the depth of the Moho, results from Receiver Functions (Novo Barbosa, 2008), and use those as constrains in inversion. Even using different parameterizations of models for the inversion, our results were very similar the mean profiles velocity structure of S-wave. In pairs of stations located in the Cear´a Central Domain Borborema the province, there are ranges of depths for which the velocities of S are very close. Most of the results in the profile near the Moho complicate their interpretation at that depth, coinciding with the geology of the region, where there are many shear zones. In particular, the profile that have the route Potiguar Bacia in inter-station, had low velocities in the crust. We combine these results to the results of gravimetry and magnetometry (Oliveira, 2008) and receptor function (Novo Barbosa, 2008). We finally, the first results on the behavior of the velocity structure of S-wave with depth in the Province Borborema

Feições crustais determinadas pela análize azimutal da função do receptor, na região da estação sismológica de Rio Claro (RCLB)

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

Espessura da zona de transição do manto sob a Província Borborema, NE Brasil, através de conversões Ps telessísmicas

The mantle transition zone is defined by two seismic discontinuities, nominally at 410 and 660 km depth, which result from transformations in the mineral olivine. The topography of these discontinuities provides information about lateral temperature changes in the transition zone. In this work, P-to-S conversions from teleseismic events recorded at 32 broadband stations in the Borborema Province were used to determine the transition zone thickness beneath this region and to investigate whether there are lateral temperature changes within this depth range. For this analysis, stacking and migration of receiver functions was performed. In the Borborema Province, geophysical studies have revealed a geoid anomaly which could reflect the presence of a thermal anomaly related to the origin of intraplate volcanism and uplift that marked the evolution of the Province in the Cenozoic. Several models have been proposed to explain these phenomena, which include those invoking the presence of a deep-seated mantle plume and those invoking shallower sources, such as small-scale convection cells. The results of this work show that no thermal anomalies are present at transition zone depths, as significant variations in the transition zone thickness were not observed. However, regions of depressed topography for both discontinuities (410 and 660 km) that approximately overlap in space were identified, suggesting that lower-thanaverage, lateral variations in seismic velocity above 410 km depth may exist below the the Borborema Province. This is consistent with the presence of a thermally-induced, low-density body independently inferred from analysis of geoid anomalies. Therefore, the magma source responsible for the Cenozoic intraplate volcanism and related uplift in the Province, is likely to be confined above the upper mantle transition zone.

Mapeamento 3D de descontinuidades sísmicas na província Borborema com funções de receptor

The crustal architecture of the Borborema Province was investigated through migration and stacking of receiver functions (phase-weighted-stack). The stacks were developed from teleseismic