共查询到20条相似文献，搜索用时 31 毫秒
1.
2.
This paper presents a method to invert underside-reflection ( P d P or S d S arrivals) data for lateral depth variations of upper-mantle discontinuities, combining traveltime and amplitude data. The method greatly improves the resolution of small-scale undulations obtained by existing imaging methods and does not suffer from the long-wavelength biases that are likely to be present in currently available models. Existing inversion methods account for the large size of the Fresnel zone of underside reflections, but not for its complexity, arising from the mini-max traveltime nature of PP- and SS -related waves. This neglect results in long-wavelength artefacts from small-scale undulations of the discontinuities, obscuring true long-wavelength depth variations. The inversion method presented in this paper uses a complex-valued sensitivity kernel, derived from the representation of underside reflections through a Kirchhoff integral formulation. The sensitivity kernel accounts for the varying sensitivity of the waveforms to discontinuity structure over the Fresnel zone. The method is applied to a large, synthetic data set. The data set consists of P d P amplitudes and traveltimes. The results show that the new inversion method resolves depth variations on a lateral scale that is smaller than the size of the Fresnel zone of individual underside reflections (but larger than the dominant wavelength), retaining the resolution of large-scale variations. The results presented here suggest that the discontinuity depth variations induced by slab penetration of the 670 discontinuity could be resolved by current broad-band P 6 7 0 P data sets. 相似文献
3.
Tadashi Kito reas Rietbrock Christine Thomas 《Geophysical Journal International》2007,169(3):1201-1209
We have developed a new array method combining conventional migration with a slowness-backazimuth deviation weighting scheme. All seismic traces are shifted based on the theoretical traveltime of the scattered wave from specific gridpoints in a 3-D volume. Observed slowness and backazimuth are calculated for each raypath and compared with theoretical values in order to estimate slowness and backazimuth deviations. Subsequently, stacked energy calculated by a conventional migration method is weighted by the slowness and backazimuth deviations to suppress any arrival energy whose slowness and backazimuth are inconsistent with the expected theoretical values. This new method was applied to two P- wave data sets which comprise (1) underside reflections at the 410 and 660 km mantle discontinuities and (2) D" reflections as well as their corresponding synthetic data sets. The results show that the weighting scheme dramatically increases the resolution of the migrated images and enables us to obtain well-constrained, focused images, making upper-mantle discontinuities and D" reflections more distinct by reducing their surrounding energy. 相似文献
4.
b
Long-period data of the Global Digital Seismograph Network (GDSN) recorded over the three-year period from 1984 to 1986 were studied for the occurrence of S-P and P-S conversions from the upper mantle transition zone that appear as precursors to teleseismic S arrivals. Conversions of this type were identified on a large number of single-station records. Simple stacking of many records enhanced the appearance of converted phases and demonstrated that no major lateral variations in the nature of the transition zone exist between various tectonic regions. S-P and P-S conversions from the 400 km discontinuity were best observed at distances between 70 and 85 while conversions from the 670 km discontinuity showed up best at distances beyond 87. The analysis of published source mechanisms and comparison with synthetic seismograms suggests that the appearance of converted phases is primarily governed by the earthquake radiation pattern. Phases that have undergone S-P conversions beneath the receiver are best observed from dip-slip events that radiate strong SV - and weak P -waves towards the station. P-S conversions beneath the source area, on the other hand, are frequently observed from events that radiate strong P and little SV energy towards the station, and also from some strike-slip events. Comparison of observed with synthetic seismograms suggests that the PREM model of Dziewonski & Anderson (1981) explains most of the observations. Observed S-P and P-S conversions from the 670 km discontinuity, however, often have larger amplitudes than in the synthetics. Constructive interference of converted waves with the P -wave coda, source radiation effects and a velocity contrast across the 670 km discontinuity which is higher than in PREM may all contribute to the discrepancy. 相似文献
Long-period data of the Global Digital Seismograph Network (GDSN) recorded over the three-year period from 1984 to 1986 were studied for the occurrence of S-P and P-S conversions from the upper mantle transition zone that appear as precursors to teleseismic S arrivals. Conversions of this type were identified on a large number of single-station records. Simple stacking of many records enhanced the appearance of converted phases and demonstrated that no major lateral variations in the nature of the transition zone exist between various tectonic regions. S-P and P-S conversions from the 400 km discontinuity were best observed at distances between 70 and 85 while conversions from the 670 km discontinuity showed up best at distances beyond 87. The analysis of published source mechanisms and comparison with synthetic seismograms suggests that the appearance of converted phases is primarily governed by the earthquake radiation pattern. Phases that have undergone S-P conversions beneath the receiver are best observed from dip-slip events that radiate strong SV - and weak P -waves towards the station. P-S conversions beneath the source area, on the other hand, are frequently observed from events that radiate strong P and little SV energy towards the station, and also from some strike-slip events. Comparison of observed with synthetic seismograms suggests that the PREM model of Dziewonski & Anderson (1981) explains most of the observations. Observed S-P and P-S conversions from the 670 km discontinuity, however, often have larger amplitudes than in the synthetics. Constructive interference of converted waves with the P -wave coda, source radiation effects and a velocity contrast across the 670 km discontinuity which is higher than in PREM may all contribute to the discrepancy. 相似文献
5.
Tadashi Kito Sebastian Rost Christine Thomas Edward J. Garnero 《Geophysical Journal International》2007,169(2):631-645
Broad-band P - and S -waves from earthquakes in South America recorded at Californian network stations are analysed to image lateral variations of the D"-discontinuity beneath the Cocos plate. We apply two array processing methods to the data set: a simplified migration method to the P -wave data set and a double-array method to both the P - and S -wave data sets, allowing us to compare results from the two methods. The double-array method images a dipping reflector at a depth range from 2650 to 2700 km in the southern part of the study area. We observe a step-like topography of 100 km to a shallower reflector at about 2600 km depth to the north, as well as evidence for a second (deeper) reflector at a depth range from 2700 to 2750 km in the north. Results from the simplified migration agree well with those from the double-array method, similarly locating a large step in reflector depth in a similar location (about 2650 km depth in the south and about 2550 km in the north) as well as the additional deeper reflector at the depth of about 2750 km in the north. Waveform modelling of the reflected waves from both methods suggests a positive velocity contrast for S waves, but a negative velocity contrast for P waves for the upper reflector in agreement with predictions from mineral physical calculations for a post-perovskite phase transition. The data also show some evidence for the existence of another deeper reflector that could indicate a double intersection of the geotherm with the post-perovskite stability field, that is, the back-transformation of post-perovskite to perovskite close to the core–mantle boundary. 相似文献
6.
Jordi Julià 《Geophysical Journal International》2007,171(1):286-301
This paper investigates the ability of P -wave receiver functions to constrain both the velocity and density contrasts across the Moho. Expressions for receiver function amplitudes corresponding to a layer over a half-space are extended to media with depth-dependent properties by explicitly accounting for free-surface reflection coefficients. Forward numerical computations show that receiver function amplitudes become frequency-dependent for depth-dependent structures and that, for a given frequency, wavelengths associated to Ps converted phases are larger than wavelengths associated to multiply reverberated phases. When consistent measurements of the amplitudes of Ps phase and multiples can be obtained, the extended expressions suffice to account for the observed amplitudes. Receiver function amplitudes are sensitive to both velocity and density contrasts across the Moho, and a two-step grid-search procedure is proposed to recover the contrasts from their measurement. The near-surface velocity is recovered in the first step from the amplitude of the direct P wave and then used in the second step to work out the relative density and S -wave velocity contrast from the amplitude of the secondary phases. Examples from central Spain (PAB) and the Indian Shield (HYB) are discussed and demonstrate that receiver function amplitudes can constrain contrasts across the Moho accurately enough to be utilized in geological interpretation. 相似文献
7.
8.
9.
Analysis of the crustal velocity structure of the British Isles using teleseismic receiver functions
J. P. Tomlinson P. Denton P. K. H. Maguire D. C. Booth 《Geophysical Journal International》2006,167(1):223-237
The onshore crustal and upper mantle velocity structure of the British Isles has been investigated by teleseismic receiver function analysis. The results of the study augment the dense offshore and sparse onshore models of the velocity structure beneath the area. In total almost 1500 receiver functions have been analysed, which have been calculated using teleseismic data from 34 broadband and short-period, three-component seismic recording instruments. The crustal structure has primarily been investigated using 1-D grid search and forward modelling techniques, returning crustal thicknesses, bulk crustal Vp / Vs ratio and velocity-depth models. H −κ stacking reveals crustal thicknesses between 25 and 36 km and Vp / Vs ratios between 1.6 and 1.9. The crustal thicknesses correlate with the results of previous seismic reflection and refraction profiles to within ±2 km. The significant exceptions are the stations close to the Iapetus Suture where the receiver function crustal thicknesses are up to 5 km less than the seismic refraction Moho. This mismatch could be linked to the presence of underplated magmatic material at the base of the crust. 1-D forward modelling has revealed subcrustal structures in northern Scotland. These correlate with results from other UK receiver function studies, and correspond with the Flannan and W-reflectors. The structures are truncated or pinch out before they reach the Midland Valley of Scotland. The isolated subcrustal structure at station GIM on the Isle of Man may be related to the closure of the Iapetus Ocean. 相似文献
10.
The causes of the major seismic discontinuities in the transition zone of the Earth's mantle 总被引：1，自引：0，他引：1
Lin-gun Liu 《Geophysical Journal International》1995,121(2):404-408
11.
12.
Convective coupling between ductile lower crust and upper mantle, and its tectonic implications 总被引：1，自引：0，他引：1
Masao Nakada 《Geophysical Journal International》1994,118(3):579-603
13.
We present a neural network approach to invert surface wave data for a global model of crustal thickness with corresponding uncertainties. We model the a posteriori probability distribution of Moho depth as a mixture of Gaussians and let the various parameters of the mixture model be given by the outputs of a conventional neural network. We show how such a network can be trained on a set of random samples to give a continuous approximation to the inverse relation in a compact and computationally efficient form. The trained networks are applied to real data consisting of fundamental mode Love and Rayleigh phase and group velocity maps. For each inversion, performed on a 2°× 2° grid globally, we obtain the a posteriori probability distribution of Moho depth. From this distribution any desired statistic such as mean and variance can be computed. The obtained results are compared with current knowledge of crustal structure. Generally our results are in good agreement with other crustal models. However in certain regions such as central Africa and the backarc of the Rocky Mountains we observe a thinner crust than the other models propose. We also see evidence for thickening of oceanic crust with increasing age. In applications, characterized by repeated inversion of similar data, the neural network approach proves to be very efficient. In particular, the speed of the individual inversions and the possibility of modelling the whole a posteriori probability distribution of the model parameters make neural networks a promising tool in seismic tomography. 相似文献
14.
15.
We present a new 2-D analytical solution of the fourth-order differential equation, which describes the flexure of a thin elastic plate.
The new analytical solution allows the differential equation for an elastic plate to be solved for any irregular shaped topography with a high spatial resolution. We apply the new method to the Central Andes. The flexural rigidity distribution calculated by this technique correlates well with tectonic units and the location of fault zones, for example, the Central Andean Gravity High correlates with the presence of a rigid, high-density body. 相似文献
The new analytical solution allows the differential equation for an elastic plate to be solved for any irregular shaped topography with a high spatial resolution. We apply the new method to the Central Andes. The flexural rigidity distribution calculated by this technique correlates well with tectonic units and the location of fault zones, for example, the Central Andean Gravity High correlates with the presence of a rigid, high-density body. 相似文献
16.
17.
18.
19.