Repository logo
 

Search Results

Now showing 1 - 4 of 4
  • 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.
  • 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.
  • 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.
  • Ambient seismic noise tomography of SW Iberia integrating seafloor- and land-based data
    Publication . Corela, Carlos; Silveira, Graça; Matias, Luis; Schimmel, Martin; Geissler, Wolfram
    We used ambient seismic noise recorded by 24 Broadband Ocean Bottom Seismometers (OBS) deployed in the Gulf of Cadiz during the EC funded NEAREST project and seven broadband land stations located in the South of Portugal to image the sedimentary and crustal structure beneath the Eastern Atlantic and SW Iberia. We computed ambient noise cross-correlations to obtain empirical Green's functions (EGFs) between all station pairs using land seismometers and both OBS sensors, seismometers and hydrophones. Despite the great difference in the recording conditions and local crustal structure between the OBSs and land stations, we could compute EGFs, by applying a linear cross-correlation with running absolute mean average time normalization, followed by a time-frequency phase weighted stack. Dispersion analysis was then applied to the EGFs, between 4 and 20s period. The obtained dispersion curves allowed mapping the lateral variation of Rayleigh-wave group velocities, as a function of period. Finally, dispersion curves extracted from each cell of the 2D group velocity maps were inverted, as a function of depth, to obtain the 3D distribution of the shear-wave velocities. The 3-D shear wave velocity model, computed from joint inversion of OBSs and land stations data allowed to estimate the thickness of sediments and crust and the Moho depth. Despite the gap that exists between the OBSs and land station locations, our model displays a good correlation with the known geological structure. The derived sedimentary layer and crustal thicknesses and the obtained Moho depth are locally in agreement with the models proposed by other studies using near vertical, refraction and wide-angle seismic profiling. We conclude that ambient noise tomography could be a valuable tool to image oceanic domains, and also that it is possible to integrate seafloor- and land-based stations to derive a structure model in the transition domain between continent and ocean.