scholarly journals Long-period modelling of MEDNET waveforms for the December 13, 1990 Eastern Sicily earthquake

1995 ◽  
Vol 38 (2) ◽  
Author(s):  
D. Giardini ◽  
B. Palombo ◽  
N. A. Pino

The availability of broad-band digital data allow "the analysis" of the seismic signal in the low-noise frequency band for structural and seismic source studies, We model complete seismograms -surface and body waves - for the December 13, 1 990, Eastern Sicily earthquake, recorded at regional distances on the MEDNET stations AQU, BNI and KEG, The inversion for the moment tensor is carried out following two approaches: a) a model-independent stralegy to fit complete seismograms in the 100-130 s period range; b) the calibration of phase velocity curves along each path to fit surface waves in the 40-100 s range, Both methodologies yield stable and consistent results: the 1990 Eastern Sicily earthquake had a seismic moment of M/o = 37 x 1024 dyne.cm, corresponding lo magnitude values of mb = 5,5, Ms = 5,7, Mw = ML = 5.8. In a second stage, synthetic seismograms arc generated by full reflectivity in the 5-50 s period range, with a new code based on the "minors" integration. By modelling surface waves in the 30-50 s range and body waves to 5 s periods, we derive average velocity models for the Central and Eastern Mediterranean, and constrain the hypocenter of the 1990 Eastern Sicily earthquake to be in the 13-17 km depth range.

1968 ◽  
Vol 58 (6) ◽  
pp. 1899-1932
Author(s):  
A. Ben-Menahem ◽  
H. Jarosch ◽  
M. Rosenman

ABSTRACT A composite compute program has been devised for a fast reduction of multistation seismic data in the period range 50–500 sec for mantle surface waves and 20–100 sec for body waves. The analysis aims at the reconstruction of the seismic source from the spectrum of its far radiation field and the correlation of its parameters with its depth, size and regional environment. The capability of the computational procedure has been demonstated in two studies of WWNSS records: one includes a spectral analysis of surface waves from a shallow shock in the Kurile Islands; the other includes a spectral analysis of P waves from 9 shocks in the depth range 550–700 km at Fiji, Mariana, Java, Japan, Peru and Brazil. Other applications of the proposed data processing routine are foreseen; a tsunami warning system and focal depth determination from spectral modal ratios. It is believed that a persistent search for stress patterns, based on the processing of a sufficiently large sample of seismic events, is essential to any future program of earthquake prediction.


1967 ◽  
Vol 57 (1) ◽  
pp. 55-81
Author(s):  
E. J. Douze

abstract This report consists of a summary of the studies conducted on the subject of short-period (6.0-0.3 sec period) noise over a period of approximately three years. Information from deep-hole and surface arrays was used in an attempt to determine the types of waves of which the noise is composed. The theoretical behavior of higher-mode Rayleigh waves and of body waves as measured by surface and deep-hole arrays is described. Both surface and body waves are shown to exist in the noise. Surface waves generally predominate at the longer periods (of the period range discussed) while body waves appear at the shorter periods at quiet sites. Not all the data could be interpreted to define the wave types present.


2021 ◽  
Author(s):  
Christina Tsarsitalidou ◽  
Pierre Boué ◽  
Gregor Hillers ◽  
Bruno Giammarinaro ◽  
Michel Campillo ◽  
...  

<div> <div> <div> <p>Dense seismic arrays are equivalent to medical ultrasound transducers in the sense that both “devices” allow the reconstruction of refocusing wave fields at near-field distances. In this work we explore the imaging potential of refocusing surface waves constructed from USArray noise correlation functions using sensors located between the US west coast and -90 degrees West. So-called focal spots---a term adopted from elastography---are constructed from the noise correlation amplitude field at zero lag time around the origin, i.e., each sensor in the array. Similar to the related SPAC method, properties of the Bessel-function-shaped focal spot are controlled by the local medium properties, which underpins the local imaging approach. Unlike USArray SPAC applications in the 5 – 40 s period range, however, we proceed in the spirit of elastographic local measurements and demonstrate the possibility to estimate properties of Rayleigh wave propagation between 80 – 300 s period using the vertical-vertical and vertical-radial focal spot components of the Green’s tensor. Clearly, the up to five-fold extension of the period range compared to noise- based USArray surface wave tomography studies are an intriguing asset of the approach that suggests a significantly increased depth resolution. In addition to demonstrating the general applicability of the focal spot method using dense array data, we address the biasing effects of less-than-ideal ambient wave field properties on our measurements. Impinging body wave energy and non- isotropic surface wave energy flux contribute to focal spot shapes and properties that are not compatible with the theoretical assumptions and used model functions and parametrizations. We show the space and period dependent distributions of these biasing components based on the focal spot representation in the wave number domain. Numerical and theoretical work discussed in an accompanying abstract is used to assess the impact on the dispersion measurements, and to test the effectiveness of filtering strategies for making improved estimates.</p> </div> </div> </div>


2020 ◽  
Vol 222 (1) ◽  
pp. 69-85 ◽  
Author(s):  
Y Lu ◽  
L Stehly ◽  
R Brossier ◽  
A Paul ◽  

SUMMARY We present an improved crustal Vs model and Moho depth map using ambient noise wave-equation tomography. The so-called ‘ambient noise wave-equation tomography’ is a method to invert seismic ambient noise phase dispersion data based on elastic waveform simulation, which accounts for 3-D and finite-frequency effects. We use cross-correlations of up to 4 yr of continuous vertical-component ambient seismic noise recordings from 304 high-quality broad-band stations in the Alpine region. We use model LSP_Eucrust1.0 obtained from traditional ambient noise tomography as initial model, and we iteratively improve the initial model by minimizing frequency-dependent phase traveltime differences between the observed and synthetic waveforms of Rayleigh waves in the period range 10–50 s. We obtain the final model after 15 iterations with ∼65 per cent total misfit reduction compared to the initial model. At crustal depth, the final model significantly enhances the amplitudes and adjusts the shapes of velocity anomalies. At Moho and upper-mantle depth, the final model corrects an obvious systematic velocity shift of the initial model. The resulting isovelocity Moho map confirms a Moho step along the external side of the external crystalline massifs of the northwestern Alps and reveals underplated gabbroic plutons in the lower most crust of the central and eastern Alps. Ambient noise wave-equation tomography turns out to be a useful tool to refine shear wave velocity models obtained by traditional ambient noise tomography based on ray theory.


Author(s):  
Fumiaki Nagashima ◽  
Hiroshi Kawase

Summary P-wave velocity (Vp) is an important parameter for constructing seismic velocity models of the subsurface structures by using microtremors and earthquake ground motions or any other geophysical exploration data. In order to reflect the ground survey information in Japan to the Vp structure, we investigated the relationships among Vs, Vp, and depth by using PS-logging data at all K-NET and KiK-net sites. Vp values are concentrated at around 500 m/s and 1,500 m/s when Vs is lower than 1,000 m/s, where these concentrated areas show two distinctive characteristics of unsaturated and saturated soil, respectively. Many Vp values in the layer shallower than 4 m are around 500 m/s, which suggests the dominance of unsaturated soil, while many Vp values in the layer deeper than 4 m are larger than 1,500 m/s, which suggests the dominance of saturated soil there. We also investigated those relationships for different soil types at K-NET sites. Although each soil type has its own depth range, all soil types show similar relationships among Vs, Vp, and depth. Then, considering the depth profile of Vp, we divided the dataset into two by the depth, which is shallower or deeper than 4 m, and calculated the geometrical mean of Vp and the geometrical standard deviation in every Vs bins of 200 m/s. Finally, we obtained the regression curves for the average and standard deviation of Vp estimated from Vs to get the Vp conversion functions from Vs, which can be applied to a wide Vs range. We also obtained the regression curves for two datasets with Vp lower and higher than 1,200 m/s. These regression curves can be applied when the groundwater level is known. In addition, we obtained the regression curves for density from Vs or Vp. An example of the application for those relationships in the velocity inversion is shown.


Geophysics ◽  
2021 ◽  
pp. 1-84
Author(s):  
Chunying Yang ◽  
Wenchuang Wang

Irregular acquisition geometry causes discontinuities in the appearance of surface wave events, and a large offset causes seismic records to appear as aliased surface waves. The conventional method of sampling data affects the accuracy of the dispersion spectrum and reduces the resolution of surface waves. At the same time, ”mode kissing” of the low-velocity layer and inhomogeneous scatterers requires a high-resolution method for calculating surface wave dispersion. This study tested the use of the multiple signal classification (MUSIC) algorithm in 3D multichannel and aliased wavefield separation. Azimuthal MUSIC is a useful method to estimate the phase velocity spectrum of aliased surface wave data, and it represent the dispersion spectra of low-velocity and inhomogeneous models. The results of this study demonstrate that mode-kissing affects dispersion imaging, and inhomogeneous scatterers change the direction of surface-wave propagation. Surface waves generated from the new propagation directions are also dispersive. The scattered surface wave has a new dispersion pattern different to that of the entire record. Diagonal loading was introduced to improve the robustness of azimuthal MUSIC, and numerical experiments demonstrate the resultant effectiveness of imaging aliasing surface waves. A phase-matched filter was applied to the results of azimuthal MUSIC, and phase iterations were unwrapped in a fast and stable manner. Aliased surface waves and body waves were separated during this process. Overall, field data demonstrate that azimuthal MUSIC and phase-matched filters can successfully separate aliased surface waves.


1995 ◽  
Vol 85 (3) ◽  
pp. 716-735 ◽  
Author(s):  
John F. Cassidy ◽  
Garry C. Rogers

Abstract On 6 April 1992, a magnitude 6.8 (MS) earthquake occurred in the triple-junction region at the northern end of the Cascadia subduction zone. This was the largest earthquake in at least 75 yr to occur along the 110-km-long Revere-Dellwood-Wilson (RDW) transform fault and the first large earthquake in this region recorded by modern broadband digital seismic networks. It thus provides an opportunity to examine the rupture process along a young (<2 Ma) oceanic transform fault and to gain better insight into the tectonics of this triple-junction region. We have investigated the source parameters and the rupture process of this earthquake by modeling broadband body waves and long-period surface waves and by accurately locating the mainshock and the first 10 days of aftershocks using a well-located “calibration” event recorded during an ocean-bottom seismometer survey. Analysis of P and SH waveforms reveals that this was a complex rupture sequence consisting of three strike-slip subevents in 12 sec. The initial rupture occurred 5 to 6 km to the SW of the seafloor trace of the RDW fault at 50.55° N, 130.46° W. The dominant subevent occurred 2 to 3 sec later and 4.3 km beneath the seafloor trace of the RDW fault, and a third subevent occurred 5 sec later, 18 km to the NNW, suggesting a northwestward propagating rupture. The aftershock sequence extended along a 60- to 70-km-long segment of the RDW fault, with the bulk of the activity concentrated ∼30 to 40 km to the NNW of the epicenter, consistent with this interpretation. The well-constrained mechanism of the initial rupture (strike/dip/slip 339°/90°/−168°) and of the largest aftershock (165°/80°/170°) are rotated 15° to 20° clockwise relative to the seafloor trace of the RDW fault but are parallel to the Pacific/North America relative plate motion vector. In contrast, the mechanisms of the dominant subevent (326°/87°/−172°), and the long-period solution derived from surface waves aligns with the RDW fault. This suggests that small earthquakes (M < 6) in this area occur along faults that are optimally aligned with respect to the regional stress field, whereas large earthquakes, involving tens of kilometers of rupture, activate the RDW fault. For the mainshock, we estimate a seismic moment (from surface waves) of 1.0 × 1026 dyne-cm, a stress drop of 60 bars, and an average slip of 1.2 m. This represents only 21 yr of strain accumulation, implying that there is either a significant amount of aseismic slip along the RDW fault or that much of the strain accumulation manifests itself as deformation within the Dellwood and Winona blocks or along the continental margin.


Author(s):  
L. Borcea ◽  
J. Garnier ◽  
K. Sølna

We derive a radiative transfer equation that accounts for coupling from surface waves to body waves and the other way around. The model is the acoustic wave equation in a two-dimensional waveguide with reflecting boundary. The waveguide has a thin, weakly randomly heterogeneous layer near the top surface, and a thick homogeneous layer beneath it. There are two types of modes that propagate along the axis of the waveguide: those that are almost trapped in the thin layer, and thus model surface waves, and those that penetrate deep in the waveguide, and thus model body waves. The remaining modes are evanescent waves. We introduce a mathematical theory of mode coupling induced by scattering in the thin layer, and derive a radiative transfer equation which quantifies the mean mode power exchange. We study the solution of this equation in the asymptotic limit of infinite width of the waveguide. The main result is a quantification of the rate of convergence of the mean mode powers toward equipartition.


Author(s):  
Tianshi Liu ◽  
Haiming Zhang

The cross-correlations of ambient noise or earthquake codas are massively used in seismic tomography to measure the dispersion curves of surface waves and the travel times of body waves. Such measurements are based on the assumption that these kinematic parameters in the cross-correlations of noise coincide with those in Green's functions. However, the relation between the cross-correlations of noise and Green's functions deserves to be studied more precisely. In this paper, we use the asymptotic analysis to study the dispersion relations of surface waves and the travel times of body waves, and come to the conclusion that for the spherically symmetric Earth model, when the distribution of noise sources is laterally uniform, the dispersion relations of surface waves and the travel times of SH body-wave phases in noise correlations should be exactly the same as those in Green's functions.


1994 ◽  
Vol 37 (3) ◽  
Author(s):  
R. G. North ◽  
C. R. D. Woodgold

An algorithm for the automatic detection and association of surface waves has been developed and tested over an 18 month interval on broad band data from the Yellowknife array (YKA). The detection algorithm uses a conventional STA/LTA scheme on data that have been narrow band filtered at 20 s periods and a test is then applied to identify dispersion. An average of 9 surface waves are detected daily using this technique. Beamforming is applied to determine the arrival azimuth; at a nonarray station this could be provided by poIarization analysis. The detected surface waves are associated daily with the events located by the short period array at Yellowknife, and later with the events listed in the USGS NEIC Monthly Summaries. Association requires matching both arrival time and azimuth of the Rayleigh waves. Regional calibration of group velocity and azimuth is required. . Large variations in both group velocity and azimuth corrections were found, as an example, signals from events in Fiji Tonga arrive with apparent group velocities of 2.9 3.5 krn/s and azimuths from 5 to + 40 degrees clockwise from true (great circle) azimuth, whereas signals from Kuriles Kamchatka have velocities of 2.4 2.9 km/s and azimuths off by 35 to 0 degrees. After applying the regional corrections, surface waves are considered associated if the arrival time matches to within 0.25 km/s in apparent group velocity and the azimuth is within 30 degrees of the median expected. Over the 18 month period studied, 32% of the automatically detected surface waves were associated with events located by the Yellowknife short period array, and 34% (1591) with NEIC events; there is about 70% overlap between the two sets of events. Had the automatic detections been reported to the USGS, YKA would have ranked second (after LZH) in terms of numbers of associated surface waves for the study period of April 1991 to September 1992.


Sign in / Sign up

Export Citation Format

Share Document