Repository logo
 

Search Results

Now showing 1 - 10 of 29
  • Sparsity-promoting approach to polarization analysis of seismic signals in the time-frequency domain
    Publication . Mohammadigheymasi, Hamzeh; Crocker, Paul; Fathi, Maryam; Almeida, Eduardo; Silveira, Graça; Gholami, Ali; Schimmel, Martin
    Time - frequency (TF)-domain polarization analysis (PA) methods are widely used as a processing tool to decompose multicomponent seismic signals. However, as a drawback, they are unable to obtain sufficient resolution to discriminate between overlapping seismic phases, as they generally rely on a low-resolution time-frequency representation (TFR) method. In this article, we present a new approach to the TF-domain PA methods. More precisely, we provide an in-detailed discussion on rearranging the eigenvalue decomposition polarization analysis (EDPA) formalism in the frequency domain to obtain the frequency-dependent polarization properties from the Fourier coefficients owing to the Fourier space orthogonality. Then, by extending the formulation to the TF domain and incorporating sparsity promoting TFR (SP-TFR), we improve the resolution when estimating the TF-domain polarization parameters. Finally, an adaptive SP-TFF is applied to extract and filter different phases of the seismic wave. By processing earthquake wave-forms, we show that, by combining amplitude, directivity, and rectilinearity attributes on the sparse TF-domain polarization map of the signal, we are able to extract (or filter) different phases of seismic waves. The SP-TFF method is evaluated on synthetic and real data associated with the source mechanism of the M-w = 8.2 earthquake that occurred in the south-southwest of Tres Picos, Mexico. A discussion on the results is given, verifying the efficiency of the method in separating not only the Rayleigh waves from the Love waves hut also in discriminating them from the body and coda waves. The codes and datasets are available at https://github.com/SigProSeismology/SP-TFF, contributing to the geoscience community.
  • An active seismic zone in intraplate West Iberia inferred from high-resolution geophysical data
    Publication . Custodio, Susana; Batlló, Josep; Zahradnik, Jiri; Arroucau, Pierre; Silveira, Graça; Heimann, Sebastian; Matos, Catarina
    Intraplate Iberia is a region of slow lithopsheric deformation (<1mm/yr) with significant historical earthquake activity. Recent high-quality instrumental data have shown that small-magnitude earthquakes collapse along clusters and lineaments, which however do not bear a clear relationship to geologically mapped active structures. In this article, we investigate the controls of these earthquake clusters. In particular, we study two of the identified clustersthe Arraiolos and the Evora seismic zones (ASZ and ESZ), located in the Western Ossa Morena Zone, southwest Iberia. The ASZ marks a sharp boundary between a seismically active region to its south and a more quiet region to its north. We revise historical earthquakes in order to clarify whether earthquake activity in the region is persistent. We use data from a local network to compute accurate epicenters, focal depth, focal mechanisms, and spatiotemporal clustering, thus characterizing ongoing small-scale fracturing. Finally, we analyze complementary data sets, including tomographic models, Global Navigation Satellite Systems data, magnetic anomalies, and gravity anomalies, in order to discuss the factors that control seismogenesis in the two seismic zones. Consistency between earthquake locations, focal mechanisms and Global Navigation Satellite Systems data suggests that the ASZ is an active right-lateral shear zone, which divides two blocks within the Western Ossa Morena Zone. The ESZ seems to localize microseismicity due to its granitic lithology. These results suggest that high-resolution geophysical data have the potential to reveal blocks with different seismogenic and rheological behaviors, which may be used to improve our understanding of fault systems and the assessment of earthquake hazard in slowly deforming regions. Plain Language Summary Mainland Portugal is a region of slow lithospheric deformation. This means that changes in Earth's outmost layerthe lithosphereoccur at very low rates (<1mm/yr). In such environments, faults producing earthquakes are not easy to identify at the Earth's surface, both because their evidence can be gradually erased by wind and water or simply because they do not reach the surface. Recent studies have shown that small earthquakes in mainland Portugal group together delineating seismically more active regions. In this article we focus in two particular groups of earthquakesthe Arraiolos and the Evora seismic zones (ASZ and ESZ) and we investigate why they occur in these particular locations. We obtain precise maps of earthquake epicenters. When possible, we also analyze the direction of slip during the earthquake and the orientation of the fracture on which it occurred. We compare our results with other data sets, such as images of the Earth's interior, that could give hints about the constitution of crust beneath the ASZ and the ESZ. Earthquakes epicenters show fault sections at depth in the ASZ. These faults separate two crustal blocks with distinct material properties. In the ESZ earthquakes are associated to contrasts in crustal materials.
  • Stratification of the earth beneath the Azores from P and S receiver functions
    Publication . Silveira, Graça; Vinnik, Lev; Stutzmann, E.; Farra, V.; Kiselev, Sergei; Morais, Iolanda
    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
    Publication . Dündar, Süleyman; Dias, Nuno; Silveira, Graça; Kind, Rainer; Vinnik, Lev; Matias, Luís; Bianchi, Marcelo
    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
    Publication . Iolanda Morais; Vinnik, Lev; Silveira, Graça; Kiselev, Sergei; Matias, Luís
    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.
  • Thermal nature of mantle upwellings below the ibero-western Maghreb region inferred from teleseismic tomography
    Publication . Civiero, Chiara; Custodio, Susana; Rawlinson, Nicholas; Strak, Vincent; Silveira, Graça; Arroucau, Pierre; Corela, Carlos
    Independent models of P wave and S wave velocity anomalies in the mantle derived from seismic tomography help to distinguish thermal signatures from those of partial melt, volatiles, and compositional variations. Here we use seismic data from SW Europe and NW Africa, spanning the region between the Pyrenees and the Canaries, in order to obtain a new S-SKS relative arrival-time tomographic model of the upper mantle below Iberia, Western Morocco, and the Canaries. Similar to previous P wave tomographic results, the S wave model provides evidence for (1) subvertical upper-mantle low-velocity structures below the Canaries, Atlas Ranges, and Gibraltar Arc, which are interpreted as mantle upwellings fed by a common lower-mantle source below the Canaries; and (2) two low-velocity anomalies below the eastern Rif and Betics that we interpret as the result of the interaction between quasi-toroidal mantle flow induced by the Gibraltar slab and the mantle upwelling behind it. The analysis of teleseismic P wave and S wave arrival-time residuals and the conversion of the low-velocity anomalies to temperature variations suggest that the upwellings in the upper mantle below the Canaries, Atlas Ranges, and Gibraltar Arc system may be solely thermal in nature, with temperature excesses in the range similar to 100-350 degrees C. Our results also indicate that local partial melting can be present at lithospheric depths, especially below the Atlas Ranges. The locations of thermal mantle upwellings are in good agreement with those of thinned lithosphere, moderate to high heat-flow measurements, and recent magmatic activity at the surface.
  • Response of Fogo volcano (Cape Verde) to lunisolar gravitational forces during the 2014-2015 eruption
    Publication . Dumont, Stéphanie; Silveira, Graça; Custódio, Susana; Lopes, Fernando; Le Mouel, Jean-Louis; Gouhier, Mathieu; GUEHENNEUX, Yannick
    Volcanoes are complex systems that evolve in space and time as a result of their eruptive activity. Volcanic eruptions represent the ultimate expression of a complex interplay between internal and external processes that span across different time scales. Deciphering how internal and external processes interact at the time scale of eruptions may provide key insights on the temporal evolution of eruptions and also help to better evaluate associated volcanic hazards. Studies of the tidal influence on volcanic activity have fallen within this context, although the cause-effect relationship between tides and eruptions is still unclear. In this study, we used Singular Spectrum Analysis to analyze three time-series, namely the seismic tremor, SO2 emission and lava volume flow rate, which cover the first month of effusive activity at Fogo volcano, Cape Verde, in 2014-2015. We detect 9 tidal periodicities and up to 5 in each time-series ranging from semi-diurnal to fortnightly periods. We show that the movement of magma at crustal depths and at surface as well as gas emission during the effusive eruption are all modulated by lunisolar gravitational forces. We highlight the relevance of the volcano location on Earth, which together with the timing of the eruption, associated with a specific astronomical configuration, result in a specific combination of tides that directly influence the volcano eruptive activity. With this data set, we further investigate the response of Fogo volcano to this external forcing. We show that during the 2014-2015 eruption, Fogo volcano acted as a bandpass filter to quasi-permanent tidal oscillations.
  • Characterization of microseismic noise in Cape Verde
    Publication . Carvalho, Joana; Silveira, Graça; Schimmel, Martin; Stutzmann, Eleonore
    The interaction of ocean waves with either the seafloor or other ocean waves generates primary (PM) and secondary microseisms (SM) that propagate through the crust and mantle, predominantly as Rayleigh waves. The horseshoe geometry and surrounding bathymetry of the Cape Verde archipelago play a significant role in the ambient-noise generation in this region. We analyze the microseisms recorded in the region using two different temporary seismic networks, and we determine the number of signals polarized as Rayleigh waves and their back azimuth (BAZ) as a function of time and frequency. The relative number of polarized signals between PM and SM varies between the stations. At most of the stations, the SM can be divided into two frequency bands. At lower frequencies (0.1-0.2 Hz), the number of SM signals is stable throughout the year, whereas at higher frequencies (0.2-0.3 Hz) this number varies with the season, with more polarized signals during the northern hemisphere spring and summer. In both frequency ranges and at most stations, the BAZ does not vary significantly over the year and points toward sources within the archipelago and outside. We compute the source site effect and show that the local bathymetry around the Cape Verde Islands strongly amplifies local SM sources. Finally, we compare the measured BAZ with source areas derived from an ocean-wave model, which confirms that Cape Verde stations mostly record local sources.
  • Educating for earthquake science and risk in a tectonically slowly deforming region
    Publication . Custódio, Susana; Silveira, Graça; Matias, Luís; Mata, Isabel; Matos, Catarina; Palma-Oliveira, José Manuel; Rocha, Francisco; Lopes, Fernando C.
    Over the past decade, scientists have been called to participate more actively in public education and outreach (E&O). This is particularly true in fields of significant societal impact, such as earthquake science. Local earthquake risk culture plays a role in the way that the public engages in educational efforts. In this article, we describe an adapted E&O program for earthquake science and risk. The program is tailored for a region of slow tectonic deformation, where large earthquakes are extreme events that occur with long return periods. The adapted program has two main goals: (1) to increase the awareness and preparedness of the population to earthquake and related risks (tsunami, liquefaction, fires, etc.), and (2) to increase the quality of earthquake science education, so as to attract talented students to geosciences. Our integrated program relies on activities tuned for different population groups who have diferente interests and abilities, namely young children, teenagers, young adults, and professionals.
  • Imaging the crust and uppermost mantle structure of Portugal (West Iberia) with seismic ambient noise
    Publication . Silveira, Graça; Dias, Nuno; Kiselev, Sergey; Stutzmann, Eleonore; Custodio, Susana; Schimmel, Martin
    We present a new high-resolution three-dimensional (3D) shear wave velocity (Vs) model of the crust and uppermost mantle beneath Portugal, inferred from ambient seismic noise tomography. We use broadband seismic data from a dense temporary deployment covering the entire Portuguese mainland between 2010 and 2012 in the scope of the WILAS project. Vertical component data are processed using phase correlation and phase weighted stack to obtain Empirical Green functions (EGF) for 3900 station pairs. Further, we use a random sampling and subset stacking strategy to measure robust Rayleigh wave group velocities in the period range 7-30 s and associated uncertainties. The tomographic inversion is performed in 2 steps: First, we determine group velocity lateral variations for each period. Next, we invert them at each grid point using a new trans-dimensional inversion scheme to obtain the 3D shear wave velocity model. The final 3D model extends from the upper crust (5 km) down to the uppermost mantle (60 km) and has a lateral resolution of similar to 50 km. In the upper and middle crust, the Vs anomaly pattern matches the tectonic units of the variscan massif and alpine basins. The transition between the Lusitanian Basin and the Ossa Morena Zone is marked by a contrast between moderate and high velocity anomalies, in addition to two arched earthquake lineations. Some faults, namely the Manteigas-Vilarica-Braganca fault and the Porto-Tomar-Ferreira do Alentejo fault, have a clear signature from the upper crust down to the uppermost mantle (60 km). Our 3D shear wave velocity model offers new insights into the continuation of the main tectonic units at depth and contributes to better understanding the seismicity of Portugal.