Velocity structure and seismic anisotropy in the crust and upper mantle from Receiver Function analysis: three case studies in Italy

The research for this PhD project consisted in the application of the RFs analysis technique to different data-sets of teleseismic events recorded at temporary and permanent stations located in three distinct study regions: Colli Albani area, Northern Apennines and Southern Apennines. We found some velocity models to interpret the structures in these regions, which possess very different geologic and tectonics characteristics and therefore offer interesting case study to face. In the Colli Albani some of the features evidenced in the RFs are shared by all the analyzed stations: the Moho is almost flat and is located at about 23 km depth, and the presence of a relatively shallow limestone layer is a stable feature; contrariwise there are features which vary from station to station, indicating local complexities. Three seismic stations, close to the central part of the former volcanic edifice, display relevant anisotropic signatures­­­ with symmetry axes consistent with the emplacement of the magmatic chamber. Two further anisotropic layers are present at greater depth, in the lower crust and the upper mantle, respectively, with symmetry axes directions related to the evolution of the volcano complex. In Northern Apennines we defined the isotropic structure of the area, finding the depth of the Tyrrhenian (almost 25 km and flat) and Adriatic (40 km and dipping underneath the Apennines crests) Mohos. We determined a zone in which the two Mohos overlap, and identified an anisotropic body in between, involved in the subduction and going down with the Adiratic Moho. We interpreted the downgoing anisotropic layer as generated by post-subduction delamination of the top-slab layer, probably made of metamorphosed crustal rocks caught in the subduction channel and buoyantly rising toward the surface. In the Southern Apennines, we found the Moho depth for 16 seismic stations, and highlighted the presence of an anisotropic layer underneath each station, at about 15-20 km below the whole study area. The moho displays a dome-like geometry, as it is shallow (29 km) in the central part of the study area, whereas it deepens peripherally (down to 45 km); the symmetry axes of anisotropic layer, interpreted as a layer separating the upper and the lower crust, show a moho-related pattern, indicated by the foliation of the layer which is parallel to the Moho trend. Moreover, due to the exceptional seismic event occurred on April 6th next to L’Aquila town, we determined the Vs model for two station located next to the epicenter. An extremely high velocity body is found underneath AQU station at 4-10 km depth, reaching Vs of about 4 km/s, while this body is lacking underneath FAGN station. We compared the presence of this body with other recent works and found an anti-correlation between the high Vs body, the max slip patches and earthquakes distribution. The nature of this body is speculative since such high velocities are consistent with deep crust or upper mantle, but can be interpreted as a as high strength barrier of which the high Vs is a typical connotation.

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.

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.

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.

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

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)

Obtención de la función receptora en las estaciones sísmicas de banda ancha del Instituto Geográfico Nacional.

La estructura de velocidad en los primeros metros debajo de una estación sísmica produce un cierto efecto sobre el sismograma que ésta registra, dicho efecto se denomina función receptora. Si se conoce esta función, es posible, mediante modelización, obtener esa estructura de velocidad. La onda P que genera un telesismo, queda “atrapada” en los últimos kilómetros de su trayectoria en la estructura próxima a la estación que la registra, reverberando entre las discontinuidades de la misma y convirtiéndose, en algunas partes, en onda S. La diferencia de registro existente entre el sismograma del rayo no atrapado (componente vertical del registro) y el producido por este fenómeno (componente radial), sería la función receptora. En este proyecto, se pretende obtener mediante la técnica de Langston (1979),las funciones receptoras de las 55 estaciones que forman la Red Sísmica de banda ancha del Instituto Geográfico Nacional, con el fin de que en trabajos venideros, sea posible conocer de una manera más realista, la estructura de la corteza terrestre bajo el territorio español. Para la obtención de estas funciones, será necesario trabajar con un conjunto de sismogramas registrados por estaciones de banda ancha a partir de telesismos cuyos parámetros focales reúnan una serie de condiciones, definidas y detalladas ampliamente en esta memoria. Tras un proceso de selección, que será constante a lo largo de todo este trabajo,se van a emplear para el cálculo de las funciones receptoras, un total de 5231 sismogramas, generados por 186 telesismos. Para terminar, a modo de ejemplo y con el objetivo de mostrar la finalidad para la que se ha llevado a cabo este proyecto, se procederá a realizar la inversión de la función receptora obtenida para una de las estaciones sísmicas de la red, obteniendo así, un modelo de velocidades aproximado de la estructura de la litosfera existente bajo dicha estación. Abstract: The speed structure in the first layers under a seismic station produces an effect on the seismogram that is registered. This effect is known as receiver function. It is possible to obtain the speed structure if this function is detected through modelling. The P wave, which generates a teleseism, gets trapped in the final kilometres of its trajectory in the nearby station structure, reverberating among its discontinuities and turning into a S wave. The registered difference between the not caught ray seismogram (vertical component of the register) and the one that is produced by this phenomenon (radial component) would be known as receiver function. This project aims to obtain, through the Langston technique (1979), the receiver function of the 55 stations that make up the Broadband Seismic Network of the National Geography Institute to get to know, in a more realistic way, the structure of the Earth’s crust under the Spanish territory. In order to obtain these functions, it will be necessary to work with a number of seismograms registered by the broadband stations whose focal parameters meet the requirements that are well defined and detailed in the development of this memory. After selecting the appropriate data, that would be constant throughout the development of this project, we are going to apply a total of 5231 seismograms obtained by 186 teleseism to calculate the receivers functions. To conclude, and with the aim of demonstrating the purposes and uses that lie behind the development of this project, we will reverse the receiver function for one of the seismic stations, obtaining a approximate speeds’ model of the lithosphere structure that can be found under the previously mentioned station.

Estrutura crustal e mantélica da província Borborema através de funções do receptor e dispersão de ondas superficiais

The Borborema Province, located in northeastern Brazil, has a basement of Precambrian age and a tectonic framework structured at the Neoproterozoic (740-560 Ma). After separation between South America and Africa during the Mesozoic, a rift system was formed, giving rise to a number of marginal and inland basins in the Province. After continental breakup, episodes of volcanism and uplift characterized the evolution of the Province. Plateau uplift was initially related to magmatic underplating of mafic material at the base of the crust, perhaps related to the generation of young continental plugs (45-7 Ma) along the Macau-Queimadas Alignment (MQA), due to a small-scale convection at the continental edge. The goal of this study is to investigate the causes of intra-plate uplift and its relationship to MQA volcanism, by using broadband seismology and integrating our results with independent geophysical and geological studies in the Borborema Province. The investigation of the deep structure of the Province with broadband seismic data includes receiver functions and surface-wave dispersion tomography. Both the receiver functions and surface-wave dispersion tomography are methods that use teleseismic events and allow to develop estimates of crustal parameters such as crustal thickness, Vp/Vs ratio, and S-velocity structure. The seismograms used for the receiver function work were obtained from 52 stations in Northeast Brazil: 16 broadband stations from the RSISNE network (Rede Sismográfica do Nordeste do Brasil), and 21 short-period and 6 broadband stations from the INCT-ET network (Instituto Nacional de Ciência e Tecnologia – Estudos Tectônicos). These results add signifi- cantly to previous datasets collected at individual stations in the Province, which include station RCBR (GSN - Global Seismic Network), stations CAUB and AGBL (Brazilian Lithosphere Seismic Project IAG/USP), and 6 other broadband stations that were part of the Projeto Milênio - Estudos geofísicos e tectônicos na Província Borborema/CNPq. For the surface-wave vii tomography, seismograms recorde at 22 broadband stations were utilized: 16 broadband stations from the RSISNE network and 6 broadband stations from the Milênio project. The new constraints developed in this work include: (i) estimates of crustal thickness and bulk Vp/Vs ratio for each station using receiver functions; (ii) new measurements of surfassewave group velocity, which were integrated to existing measurementes from a continental-scale tomography for South America, and (iii) S-wave velocity models (1D) at various locations in the Borborema Province, developed through the simultaneous inversion of receiver functions and surface-wave dispersion velocities. The results display S-wave velocity structure down to the base of the crust that are consistent with the presence of a 5-7.5 km thick mafic layer. The mafic layer was observed only in the southern portion of the Plateau and absent in its northern portion. Another important observation is that our models divide the plateau into a region of thin crust (northern Plateau) and a region of thick crust (southern Plateau), confirming results from independent refraction surveys and receiver function analyses. Existing models of plateau uplift, nonetheless, cannot explain all the new observations. It is proposed that during the Brazilian orogeny a layer of preexisting mafic material was delaminated, as a whole or in part, from the original Brasiliano crust. Partial delamination would have happened in the southern portion of the plateau, where independent studies found evidence of a more resistant rheology. During Mesozoic rifting, thinning of the crust around the southern Plateau would have formed the marginal basins and the Sertaneja depression, which would have included the northern part of the Plateau. In the Cenozoic, uplift of the northern Plateau would have occurred, resulting in a northern Plateau without mafic material at the base of the crust and a southern Plateau with partially delaminated mafic layer.

Crustal-scale shear zones and heterogeneous structure beneath the North Anatolian Fault Zone, Turkey, revealed by a high-density seismometer array

Peer reviewed

Crustal-scale shear zones and heterogeneous structure beneath the North Anatolian Fault Zone, Turkey, revealed by a high-density seismometer array

Peer reviewed

EMA analysis of tongue function in children with dysarthria following traumatic brain injury

Primary objective : To investigate the speed and accuracy of tongue movements exhibited by a sample of children with dysarthria following severe traumatic brain injury (TBI) during speech using electromagnetic articulography (EMA). Methods and procedures : Four children, aged between 12.75-17.17 years with dysarthria following TBI, were assessed using the AG-100 electromagnetic articulography system (Carstens Medizinelektronik). The movement trajectories of receiver coils affixed to each child's tongue were examined during consonant productions, together with a range of quantitative kinematic parameters. The children's results were individually compared against the mean values obtained by a group of eight control children (mean age of 14.67 years, SD 1.60). Main outcomes and results : All four TBI children were perceived to exhibit reduced rates of speech and increased word durations. Objective EMA analysis revealed that two of the TBI children exhibited significantly longer consonant durations compared to the control group, resulting from different underlying mechanisms relating to speed generation capabilities and distances travelled. The other two TBI children did not exhibit increased initial consonant movement durations, suggesting that the vowels and/or final consonants may have been contributing to the increased word durations. Conclusions and clinical implications : The finding of different underlying articulatory kinematic profiles has important implications for the treatment of speech rate disturbances in children with dysarthria following TBI.