首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The earthquakes in the seismicity belt extending through Indonesia, New Guinea, Vanuatu and Fiji to the Tonga–Kermadec subduction zone recorded at the 65 portable broad-band stations deployed during the Skippy experiment from 1993–1996 provide good coverage of the lithosphere and mantle under the Australian continent, Coral Sea and Tasman Sea.
The variation in structure in the upper part of the mantle is characterized by deter-mining a suite of 1-D structures from stacked record sections utilizing clear P and S arrivals, prepared for all propagation paths lying within a 10° azimuth band. The azimuth of these bands is rotated by 20° steps with four parallel corridors for each azimuth. This gives 26 separate azimuthal corridors for which 15 independent 1-D seismic velocity structures have been derived, which show significant variation in P and S structure.
The set of 1-D structures is combined to produce a 3-D representation by projecting the velocity values along the ray path using a turning point approximation and stacking into 3-D cells (5° by 50 km in depth). Even though this procedure will tend to underestimate wave-speed perturbations, S -velocity deviations from the ak135 reference model exceed 6 per cent in the lithosphere.
In the uppermost mantle the results display complex features and very high S -wave speeds beneath the Precambrian shields with a significant low-velocity zone beneath. High velocities are also found towards the base of the transition zone, with high S -wave speeds beneath the continent and high P -wave speeds beneath the ocean. The wave-speed patterns agree well with independent surface wave studies and delay time tomography studies in the zones of common coverage.  相似文献   

2.
P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.  相似文献   

3.
Shear-wave splitting from local deep earthquakes is investigated to clarify the volume and the location of two anisotropic bodies in the mantle wedge beneath central Honshu, Japan. We observe a spatial variation in splitting parameters depending on the combination of sources and receivers, nearly N–S fast in the northern region, nearly E–W fast in the southern region and small time delays in the eastern region. Using forward modelling, two models with 30 and 10 per cent anisotropy are tested by means of a global search for the locations of anisotropic bodies with various volumes. The optimum model is obtained for 30 per cent anisotropy, which means a 5 per cent velocity difference between fast and slow polarized waves. The northern anisotropic body has a volume of 1.00° (longitude) × 0.5° (latitude) × 75 km (depth), with the orientation of the symmetry axis being N20°E. The southern anisotropic body has a volume of 1.25° × 1.25° × 100 km with the symmetry axis along N95°E. Our results show that the anisotropic bodies are located in low-velocity and low- Q regions of the mantle. This, together with petrological data and the location of volcanoes in the arc, suggests that the possible cause of the anisotropy is the preferred alignment of cracks filled with melt.  相似文献   

4.
It is well established that the Earth's uppermost mantle is anisotropic, but there are no clear observations of anisotropy in the deeper parts of the mantle. Surface waves are well suited to observe anisotropy since they carry information about both radial and azimuthal anisotropy. Fundamental mode surface waves, for commonly used periods up to 200 s, are sensitive to structure in the first few hundred kilometres, and therefore, do not provide information on anisotropy below. Higher mode surface waves have sensitivities that extend to and beyond the transition zone, and should thus give insight about azimuthal anisotropy at greater depths. We have measured higher mode Love and Rayleigh phase velocities using a model space search approach, which provides us with consistent relative uncertainties from measurement to measurement and from mode to mode. From these phase velocity measurements, we constructed global anisotropic phase velocity maps. Prior to inversion, we determine the optimum relative weighting for anisotropy. We present global azimuthal phase velocity maps for higher mode Rayleigh waves (up to the sixth higher mode) and Love waves (up to the fifth higher mode) with corresponding average model uncertainties. The anisotropy we derive is robust within the uncertainties for all modes. Given the ray theoretical sensitivity kernels of Rayleigh and Love wave modes, the source of anisotropy is complex, but mainly located in the asthenosphere and deeper. Our models show a good correspondence with other studies for the fundamental mode, but we have been able to achieve higher resolution.  相似文献   

5.
Summary. Seismic anisotropy has been previously studied at depths usually not exceeding 100 or 150 km. In this paper we present a method of analysis of seismic records which is very sensitive to azimuthal anisotropy and is applicable at almost any depth range. The idea of the method is to detect and analyse the SH -component of the waves, converted from P to S in the mantle. The procedure of record processing includes frequency filtering, axis rotation, transformation of the record to a standard form, stacking the standardized SH -component records of many seismic events, and the harmonic analysis of amplitude as a function of the direction of wave propagation. When applied to the long-period records of NORSAR the procedure detected a converted wave with the properties implying the possibility of its propagation in a transversely isotropic medium with a horizontal axis of symmetry . Our preferred model postulates anisotropy of ∼ 1 per cent in a layer 50 km thick at the base of the upper mantle.  相似文献   

6.
The deployment of temporary arrays of broadband seismological stations over dedicated targets is common practice. Measurement of surface wave phase velocity across a small array and its depth-inversion gives us information about the structure below the array which is complementary to the information obtained from body-wave analysis. The question is however: what do we actually measure when the array is much smaller than the wave length, and how does the measured phase velocity relates to the real structure below the array? We quantify this relationship by performing a series of numerical simulations of surface wave propagation in 3-D structures and by measuring the apparent phase velocity across the array on the synthetics. A principal conclusion is that heterogeneities located outside the array can map in a complex way onto the phase velocities measured by the array. In order to minimize this effect, it is necessary to have a large number of events and to average measurements from events well-distributed in backazimuth. A second observation is that the period of the wave has a remarkably small influence on the lateral resolution of the measurement, which is dominantly controlled by the size of the array. We analyse if the artefacts created by heterogeneities can be mistaken for azimuthal variations caused by anisotropy. We also show that if the amplitude of the surface waves can be measured precisely enough, phase velocities can be corrected and the artefacts which occur due to reflections and diffractions in 3-D structures greatly reduced.  相似文献   

7.
A 3-D P -velocity map of the crust and upper mantle beneath the southeastern part of India has been reconstructed through the inversion of teleseismic traveltimes. Salient geological features in the study region include the Archean Dharwar Craton and Eastern Ghat metamorphic belt (EGMB), and the Proterozoic Cuddapah and Godavari basins. The Krishna–Godavari basin, on the eastern coastal margin, evolved in response to the Indo–Antarctica breakup. A 24-station temporary network provided 1161 traveltimes, which were used to model 3-D P -velocity variation. The velocity model accounts of 80 per cent of the observed data variance. The velocity picture to a depth of 120 km shows two patterns: a high velocity beneath the interior domain (Dharwar craton and Cuddapah basin), and a lower velocity beneath the eastern margin region (EGMB and coastal basin). Across the array velocity variations of 7–10 per cent in the crust (0–40 km) and 3–5 per cent in the uppermost mantle (40–120 km) are observed. At deeper levels (120–210 km) the upper-mantle velocity differences are insignificant among different geological units. The presence of such a low velocity along the eastern margin suggests significantly thin lithosphere (<100 km) beneath it compared to a thick lithosphere (>200 km) beneath the eastern Dharwar craton. Such lithospheric thinning could be a consequence of Indo–Antarctica break-up.  相似文献   

8.
Borehole logs often show considerable variation over small lengths, consistent with fine layering of the formation. Elastic upscaling allows layering on a fine length scale to be replaced by thicker, homogeneous equivalent transversely isotropic layers, thus creating a model with far fewer layers. In this paper, the effect of inclusions of finite lateral extent on elastic upscaling is examined. For inclusions with a small thickness to width ratio, Backus averaging gives a good approximation to the long-wavelength elastic properties of the medium. However, for larger thickness to width ratios, the anisotropy of the medium can be significantly reduced. For inclusions that are long in comparison to their width, azimuthal anisotropy may result if the inclusions show a preferential orientation. This could result from the presence of flow during deposition. The azimuthal anisotropy increases with increasing thickness to width ratio of the inclusions, and vanishes when this ratio is zero. For larger thickness to width ratios, the magnitude of the azimuthal anisotropy may be of similar magnitude to that commonly seen in sedimentary basins. Azimuthal anisotropy is usually attributed to the presence of aligned fractures within the formation. The present results suggest a further source of azimuthal anisotropy that may be significant.  相似文献   

9.
Finite-frequency sensitivity kernels for head waves   总被引：2，自引：0，他引：2  
Head waves are extremely important in determining the structure of the predominantly layered Earth. While several recent studies have shown the diffractive nature and the 3-D Fréchet kernels of finite-frequency turning waves, analogues of head waves in a continuous velocity structure, the finite-frequency effects and sensitivity kernels of head waves are yet to be carefully examined. We present the results of a numerical study focusing on the finite-frequency effects of head waves. Our model has a low-velocity layer over a high-velocity half-space and a cylindrical-shaped velocity perturbation placed beneath the interface at different locations. A 3-D finite-difference method is used to calculate synthetic waveforms. Traveltime and amplitude anomalies are measured by the cross-correlation of synthetic seismograms from models with and without the velocity perturbation and are compared to the 3-D sensitivity kernels constructed from full waveform simulations. The results show that the head wave arrival-time and amplitude are influenced by the velocity structure surrounding the ray path in a pattern that is consistent with the Fresnel zones. Unlike the 'banana–doughnut' traveltime sensitivity kernels of turning waves, the traveltime sensitivity of the head wave along the ray path below the interface is weak, but non-zero. Below the ray path, the traveltime sensitivity reaches the maximum (absolute value) at a depth that depends on the wavelength and propagation distance. The sensitivity kernels vary with the vertical velocity gradient in the lower layer, but the variation is relatively small at short propagation distances when the vertical velocity gradient is within the range of the commonly accepted values. Finally, the depression or shoaling of the interface results in increased or decreased sensitivities, respectively, beneath the interface topography.  相似文献   

10.
A new formulation for the propagation of surface waves in three-dimensionally varying media is developed in terms of modal interactions. A variety of assumptions can be made about the nature of the modal field: a single set of reference modes, a set of local modes for the structure beneath a point, or a set of local modes for a laterally varying reference structure. Each modal contribution is represented locally as a spectrum of plane waves propagating in different directions in the horizontal plane. The influence of 3-D structure is included by allowing coupling between different modal branches and propagation directions. For anisotropic models, with allowance for attenuation, the treatment leads to a set of coupled 2-D partial differential equations for the weight functions for different modal orders.
The representation of the guided wavefield requires the inclusion of a full set of modes, so that, even for isotropic models, both Love and Rayleigh modes appear as different polarization states of the modal spectrum. The coupling equations describe the interaction between the different polarizations induced by the presence of the 3-D structure.
The level of lateral variation within the 3-D model is not required to be small. Horizontal refraction or reflection of the surface wavefield can be included by allowing for transfer between modes travelling in different directions. Approximate forms of the coupled equation system can be employed when the level of heterogeneity is small, for example the coupling between the fundamental mode and higher modes can often be neglected, or forward propagation can be emphasized by restricting the interaction to a limited band of plane waves covering the expected direction of propagation.  相似文献   

11.
In contrast to previous work, which treats the Earth's lateral heterogeneity as an infinitesimal perturbation to a spherically symmetrical starting model, we conduct iterative linearized waveform inversion for the Earth's laterally heterogeneous structure. We use the Direct Solution Method (DSM) (Geller et al. 1990a) to calculate synthetic seismograms and their partial derivatives for a laterally heterogeneous earth model. We invert surface-wave data from the IDA and GEOSCOPE networks. We expand the lateral heterogeneity of rigidity in spherical harmonics up to angular order number 8 and use three parameters to specify the depth dependence of each harmonic, giving us a total of 240 unknowns. The short-wavelength lateral heterogeneity (s = 4, 6 and 8) in the deeper part of the upper mantle obtained by this study differs significantly from M84A. The relative improvement in the variance reduction as compared with model M84A is about 20 per cent for the IDA data and more than 100 per cent for the GEOSCOPE data.  相似文献   

12.
P-wave velocity anisotropy in crystalline rocks   总被引：1，自引：0，他引：1  
Summary. Compressional wave velocities and anisotropy coefficients determined at high hydrostatic pressures are compiled from the data published for the main types of crystalline rocks. The crack-free elastic anisotropy of igneous crustal rocks is generally very low, between 1 and 3 per cent on average. The anisotropy of metamorphic rocks is higher (up to 22 per cent), but very variable. The average anisotropy coefficients in schists and amphibolites are about 10 per cent, in gneisses between 3 and 7 per cent, and in granulites less than 3 per cent. The average anisotropy of olivine ultramafites is between 7 and 12 per cent, whereas in pyroxenites and eclogites it is usually less than 4 per cent. A comparison of ranges of average velocities and average anisotropies for the individual rock groups suggests that, whereas in the crust the lateral velocity variations are mainly due to compositional changes, in the olivine of the uppermost mantle the velocity variations due to anisotropic structures could be of the same magnitude as the variations due to inhomogeneities.  相似文献   

13.
b
Rayleigh-wave phase velocities are investigated in the period range 17–100 s by the two-station method over several paths covering most of French Polynesia. Our results confirm the validity of theoretical models obtained through regionalization of data pertaining to longer paths. They also exhibit a 2–3.5 per cent anisotropy, with the axis of maximum velocity oriented in the direction of spreading of the plate. Part of this anisotropy is, however, due to the presence of the Tuamotu archipelago; when this is removed, the remaining anisotropy (about 1.5 per cent) correlates with the present direction of spreading, indicating that a relaxation of the anisotropy has taken place since the East Pacific ridge jump. Finally, the presence of the Tuamotu Islands explains anomalous waveshapes for surface waves travelling in their vicinity, due to multipathing through their faster structure.  相似文献   

14.
We develop an approach that allows us to invert for the mantle velocity structure within a finely parametrized region as a perturbation with respect to a low-resolution, global tomographic model. We implement this technique to investigate the upper-mantle structure beneath Eurasia and present a new model of shear wave velocity, parametrized laterally using spherical splines with ∼2.9° spacing in Eurasia and ∼11.5° spacing elsewhere. The model is obtained from a combined data set of surface wave phase velocities, long-period waveforms and body-wave traveltimes. We identify many features as narrow as few hundred kilometres in diameter, such as subducting slabs in eastern Eurasia and slow-velocity anomalies beneath tectonically active regions. In contrast to regional studies in which these features have been identified, our model encompasses the structure of the entire Eurasian continent. Furthermore, including mantle- and body-wave waveforms helped us constrain structures at depths larger than 250 km, which are poorly resolved in earlier models. We find that up to +9 per cent faster-than-average anomalies within the uppermost ∼200 km of the mantle beneath cratons and some orogenic regions are separated by a sharp gradient zone from deeper, +1 to +2 per cent anomalies. We speculate that this gradient zone may represent a boundary separating the lithosphere from the continental root, which might be compositionally distinct from the overlying lithosphere and remain stable either due to its compositional buoyancy or due to higher viscosity compared with the suboceanic mantle. Our regional model of anisotropy is not significantly different from the global one.  相似文献   

15.
Summary. Available seismic refraction data from three different continental areas, northern Britain and the eastern and western United States, has been studied for possible Pn , velocity anisotropy using the methods described by Bamford. There are various deficiencies in the time—distance data used in each case but, while the uppermost mantle beneath northern Britain and the eastern United States seems to be isotropic within the limits of measurement error, there is a small but significant anisotropy beneath the western United States.
Both the amount (up to 3 per cent) and the direction (70–80° east of north) of this anisotropy are very similar to the results obtained in the Pacific Ocean off California. We tentatively conclude that this anisotropy is present as a consequence of the subduction of oceanic lithosphere beneath the western United States.  相似文献   

16.
Summary. Asymptotic ray theory is applied to surface waves in a medium where the lateral variations of structure are very smooth. Using ray-centred coordinates, parabolic equations are obtained for lateral variations while vertical structural variations at a given point are specified by eigenfunctions of normal mode theory as for the laterally homogeneous case. Final results on wavefields close to a ray can be expressed by formulations similar to those for elastic body waves in 2-D laterally heterogeneous media, except that the vertical dependence is described by eigenfunctions of 'local' Love or Rayleigh waves. The transport equation is written in terms of geometrical-ray spreading, group velocity and an energy integral. For the horizontal components there are both principal and additional components to describe the curvature of rays along the surface, as in the case of elastic body waves. The vertical component is decoupled from the horizontal components. With complex parameters the solutions for the dynamic ray tracing system correspond to Gaussian beams: the amplitude distribution is bell-shaped along the direction perpendicular to the ray and the solution is regular everywhere, even at caustics. Most of the characteristics of Gaussian beams for 2-D elastic body waves are also applicable to the surface wave case. At each frequency the solution may be regarded as a set of eigenfunctions propagating over a 2-D surface according to the phase velocity mapping.  相似文献   

17.
Summary. Numerical modelling is one of the most efficient methods for an investigation of the relationship between structural features and peculiarities of observed wavefields. It is practically the only method for 2-D and 3-D inhomogeneous media.
An algorithm based on ray theory has been developed for calculations of travel times and amplitudes of seismic waves in 3-D inhomogeneous media with curved interfaces. It was applied for numerical modelling of kinematic and dynamic characteristics of seismic waves propagating in laterally inhomogeneous media.
Travel-time and amplitude patterns were studied in the 2-D and 3-D models of a geosyncline, in which velocity distribution was given by an analytical function of the coordinates. For a more complicated model representing a subducting high-velocity lithospheric plate in a transition zone between oceanic and continental upper mantle, the velocity distribution was given by discrete values on a 2-D non-rectangular grid. It was shown that when a source was placed above the lithospheric plate, a shadow zone appeared along a strike of the structure, i.e. in the direction which is perpendicular to a strong lateral velocity gradient. Travel-time residuals were calculated along the seismological profile for a 3-D velocity distribution in the upper mantle beneath Central Asia, obtained as a result of inversion of travel times by the Backus-Gilbert method. They were found to be in a good agreement with the observed data.  相似文献   

18.
Amplitude measurements of the transverse component of SKS waves, the so-called splitting intensity, can be used to formulate a non-linear inverse problem to image the 3-D variations of upper mantle anisotropy. Assuming transverse isotropy (or hexagonal symmetry), one can parametrize anisotropy by two anisotropic parameters and two angles describing the orientation of the symmetry axis. These can also be written as two collinear pseudo-vectors. The tomographic process consists of retrieving the spatial distribution of these pseudo-vectors, and thus resembles surface wave vectorial tomography. Spatial resolution results from the sensitivity of low-frequency SKS waves to seismic anisotropy off the ray path. The expressions for the 3-D sensitivity kernels for splitting intensity are derived, including the near-field contributions, and validated by comparison with a full wave equation solution based upon the finite element method. These sensitivity kernels are valid for any orientation of the symmetry axis, and thus generalize previous results that were only valid for a horizontal symmetry axis. It is shown that both lateral and vertical subwavelength variations of anisotropy can be retrieved with a dense array of broad-band stations, even in the case of vertically propagating SKS waves.  相似文献   

19.
Seismic imaging of the laterally varying D" region beneath the Cocos Plate   总被引：1，自引：0，他引：1  
We use an axisymmetric, spherical Earth finite difference algorithm to model SH -wave propagation through cross-sections of laterally varying lower mantle models beneath the Cocos Plate derived from recent data analyses. Synthetic seismograms with dominant periods as short as 4 s are computed for several models: (1) a D" reflector 264 km above the core–mantle boundary with laterally varying S -wave velocity increases of 0.9–2.6 per cent, based on localized structures from a 1-D double-array stacking method; (2) an undulating D" reflector with large topography and uniform velocity increase obtained using a 3-D migration method and (3) cross-sections through the 3-D mantle S -wave velocity tomography model TXBW. We apply double-array stacking to assess model predictions of data. Of the models explored, the S -wave tomography model TXBW displays the best overall agreement with data. The undulating reflector produces a double Scd arrival that may be useful in future studies for distinguishing between D" volumetric heterogeneity and D" discontinuity topography. Synthetics for the laterally varying models show waveform variability not observed in 1-D model predictions. It is challenging to predict 3-D structure based on localized 1-D models when lateral structural variations are on the order of a few wavelengths of the energy used, particularly for the grazing geometry of our data. Iterative approaches of computing synthetic seismograms and adjusting model characteristics by considering path integral effects are necessary to accurately model fine-scale D" structure.  相似文献   

20.
Teleseismic P waves passing through low-wave-speed bodies in the mantle are refracted, causing anomalies in their propagation directions that can be measured by seismometer arrays. Waves from earthquakes in the eastern Pacific and western North America arriving at the NORSAR array in Norway and at seismic stations in Scotland pass beneath the Iceland region at depths of ∼ 1000–2000 km. Waves arriving at NORSAR have anomalous arrival azimuths consistent with a low-wave-speed body at a depth of ∼ 1500 km beneath the Iceland–Faeroe ridge with a maximum diameter of ∼250 km and a maximum wave-speed contrast of ∼ 1.5 per cent. This agrees well with whole-mantle tomography results, which image a low-wave-speed body at this location with a diameter of ∼ 500 km and a wave-speed anomaly of ∼ 0.5 per cent, bearing in mind that whole-mantle tomography, because of its limited resolution, broadens and weakens small anomalies. The observations cannot resolve the location of the body, and the anomaly could be caused in whole or in part by larger bodies farther away, for example by a body imaged beneath Greenland by whole-mantle tomography.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号