首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 734 毫秒
1.
The velocity spectrum stacking method is applied to receiver functions from stations ATD and AAE to image P -to- S converted phases originating at the 410 and 660 km discontinuities beneath Afar. A transition zone thickness of 244 ± 19 km is obtained, similar to the global average transition zone thickness. This result suggests that any broad thermal anomaly beneath Afar probably does not extend as far down as the transition zone. However, because of the 19 km uncertainty in the thickness estimate, a small thermal anomaly of ~100–150 K at mantle transition zone depths cannot be ruled out.  相似文献   

2.
Global mapping of 410 and 660 km discontinuity topography and transition zone thickness has proven to be a powerful tool for constraining mantle chemistry, dynamics and mineralogy. Numerous seismic and mineral physics studies suggest that the 410 km discontinuity results from the phase change of olivine to wadsleyite and the 660 km discontinuity results from the phase change of ringwoodite to perovskite and magnesiowustite. Underside reflections of the 410 and 660 km discontinuities arrive as precursors to SS . With the recent development of a semi-automated method of determining SS arrivals, we have more than tripled the Flanagan and Shearer (1998a) data set of handpicked SS waveforms. We are able to increase resolution by stacking waveforms in 5° rather than 10° radius bins as well as increasing data coverage significantly in the southern hemisphere. The resulting SS-S410S and SS-S660S times are heavily influenced by upper-mantle velocity structure. We perform a joint inversion for discontinuity topography and velocity heterogeneity as well as performing a simple velocity correction to the precursor differential times and find little difference between the two methods. The 660 km discontinuity topography and transition zone thickness are correlated with velocities in the transition zone whereas the 410 km discontinuity topography is not. In addition, the 410 km discontinuity topography is not correlated with the 660 km discontinuity topography, rather anticorrelated, as expected due to the opposite signs of the Clapeyron slopes of their respective phase changes. These results suggest that, whereas the topography of 660 km discontinuity could be dominated by thermal effects, the topography of the 410 km discontinuity is likely dominated by compositional effects. In addition, unlike previous studies which find less topography on the 410 km discontinuity than on the 660 km discontinuity, our 410 and 660 km topography have similar amplitudes.  相似文献   

3.
We study the crustal structure of eastern Marmara region by applying the receiver function method to the data obtained from the 11 broad-band stations that have been in operation since the 1999 İzmit earthquake. The stacked single-event receiver functions were modelled by an inversion algorithm based on a five-layered crustal velocity model to reveal the first-order shear-velocity discontinuities with a minimum degree of trade-off. We observe crustal thickening from west (29–32 km) to east (34–35 km) along the North Anatolian Fault Zone (NAFZ), but we observe no obvious crustal thickness variation from north to south while crossing the NAFZ. The crust is thinnest beneath station TER (29 km), located near the Black Sea coast in the west and thickest beneath station TAR (35 km), located inland in the southeast. The average crustal thickness and S -wave velocity for the whole regions are  31 ± 2  km and  3.64 ± 0.15 km s−1  , respectively. The eastern Marmara region with its average crustal thickness, high heat flow value (101 ± 11 mW m−2) and with its remarkable extensional features seems to have a Basin and Range type characteristics, but the higher average shear velocities (∼3.64 km s−1) and crustal thickening from 29 to 35 km towards the easternmost stations indicate that the crustal structure shows a transitional tectonic regime. Therefore, we conclude that the eastern Marmara region seems to be a transition zone between the Marmara Sea extensional domain and the continental Anatolian inland region.  相似文献   

4.
Data from 90 permanent broad-band stations spread over central and eastern Europe were analysed using Ps receiver functions to study the crustal and upper-mantle structure down to the mantle transition zone. Receiver functions provide valuable information on structural features, which are important for the resolution of European lithospheric dynamics. Moho depths vary from less than 25 km in extensional areas in central Europe to more than 50 km at stations in eastern Europe (Craton) and beneath the Alpine–Carpathian belt. A very shallow Moho depth can be observed at stations in the Upper Rhine Graben area ( ca. 25 km), whereas, for example, stations in the SW Bohemian Massif show a significantly deeper Moho interface at a depth of 38 km. Vp / Vs ratios vary between 1.60 and 1.96, and show no clear correlation to the major tectonic units, thus probably representing local variations in crustal composition. Delayed arrivals of converted phases from the mantle transition zone are observed at many stations in central Europe, whereas stations in the cratonic area show earlier arrivals compared with those calculated from the IASP91 Earth reference model. Differential delay times between the P410s and P660s phases indicate a thickened mantle transition zone beneath the eastern Alps, the Carpathians and the northern Balkan peninsula, whereas the transition zone thickness in eastern and central Europe agrees with the IASP91 value. The thickening of the mantle transition zone beneath the eastern Alps and the Carpathians could be caused by cold, deeply subducted oceanic slabs.  相似文献   

5.
We present a 3-D radially anisotropic S velocity model of the whole mantle (SAW642AN), obtained using a large three component surface and body waveform data set and an iterative inversion for structure and source parameters based on Non-linear Asymptotic Coupling Theory (NACT). The model is parametrized in level 4 spherical splines, which have a spacing of ∼ 8°. The model shows a link between mantle flow and anisotropy in a variety of depth ranges. In the uppermost mantle, we confirm observations of regions with   VSH > VSV   starting at ∼80 km under oceanic regions and ∼200 km under stable continental lithosphere, suggesting horizontal flow beneath the lithosphere. We also observe a   VSV > VSH   signature at ∼150–300 km depth beneath major ridge systems with amplitude correlated with spreading rate for fast-spreading segments. In the transition zone (400–700 km depth), regions of subducted slab material are associated with   VSV > VSH   , while the ridge signal decreases. While the mid-mantle has lower amplitude anisotropy (<1 per cent), we also confirm the observation of radially symmetric   VSH > VSV   in the lowermost 300 km, which appears to be a robust conclusion, despite an error in our previous paper which has been corrected here. The 3-D deviations from this signature are associated with the large-scale low-velocity superplumes under the central Pacific and Africa, suggesting that   VSH > VSV   is generated in the predominant horizontal flow of a mechanical boundary layer, with a change in signature related to transition to upwelling at the superplumes.  相似文献   

6.
Global mapping of upper mantle reflectors from long-period SS precursors   总被引:1,自引:0,他引:1  
Long-period precursors to SS resulting from underside reflections off upper mantle discontinuities ( SdS where d is the discontinuity depth) can be used to map the global distribution and depth of these reflectors. We analyse 5,884 long-period seismograms from the Global Digital Seismograph Network (1976-1987, shallow sources, transverse component) in order to identify SdS arrivals. Corrections for velocity dispersion, topography and crustal thickness at the SS bounce point, and lateral variation in mantle velocity are critical for obtaining accurate estimates of discontinuity depths. The 410 and 660 km discontinuities are observed at average depths of 413 and 653 km, and exhibit large-scale coherent patterns of topography with depth variations up to 40 km. These patterns are roughly correlated with recent tomographic models, with fast anomalies in the transition zone associated with highs in the 410 km discontinuity and lows in the 660 km discontinuity, a result consistent with laboratory measurements of Clapeyron slopes for the appropriate phase changes. The best resolved feature in these maps is a trough in the 660 km discontinuity in the northwest Pacific, which appears to be associated with the subduction zones in this region. Amplitude variations in SdS arrivals are not correlated with discontinuity depths and probably result from focusing and defocusing effects along the ray paths. The SdS arrivals suggest the presence of regional reflectors in the upper mantle above 400 km. However, only the strongest of these features are above probable noise levels due to sampling inadequacies.  相似文献   

7.
We present velocity constraints for the upper-mantle transition zones beneath Central Siberia based on observations of the 1982 RIFT Deep Seismic Sounding (DSS) profile. The data consist of seismic recordings of a nuclear explosion in north-western Siberia along a 2600 km long seismic profile extending from the Yamal Peninsula to Lake Baikal. We invert seismic data from the mantle transition zones using a non-linear inversion scheme using a genetic algorithm for optimization and the WKBJ method to compute the synthetic seismograms. A statistical error analysis using a graph-binning technique was performed to provide uncertainty values in the velocity models.
Our best model for the upper-mantle velocity discontinuity near 410 km depth has a two-stage velocity-gradient structure, with velocities increasing from 8.70–9.25 km s−1 over a depth range of 400–415 km, a gradient of 0.0433 s−1, and from 9.25–9.60 km s−1 over a depth range of 415–435 km, a gradient of 0.0175 s−1. This derived model is consistent with other seismological observations and mineral-physics models. The model for the velocity discontinuity near 660 km depth is simple, sharp and includes velocities increasing from 10.15 km s−1 at 655 km depth to 10.70 km s−1 at 660 km depth, a gradient of 0.055 s−1.  相似文献   

8.
Seismic phase conversions provide important constraints on the layered nature of subduction zone structures. Recordings from digital stations in North Island, New Zealand, have been examined for converted ScS ‐to‐ p ( ScSp ) arrivals from deep (>150 km) Tonga–Kermadec earthquakes to image layering in the underlying Hikurangi subduction zone. Consistent P ‐wave energy prior to ScS has been identified from stations in eastern and southern North Island, where the subducted plate interface is at a depth of between 15 and 30 km. Two ScS precursors are observed. Ray tracing indicates that the initial precursor ( ScSp 1) corresponds to conversion from the base of an 11–14 km thick subducting Pacific crust. The second precursor is interpreted as a conversion from the top of the subducting plate. The amplitude ratio, ScSp 1: ScS , increases from 0.10 to 0.19 from northern to southern North Island. This is within the range expected from a simple first‐order velocity discontinuity at an oceanic Moho. A 1–2 km thick layer of low‐velocity sediment at the top of the subducting plate is required to explain the remaining ScSp waveform. Our results imply that the abnormally thick Hikurangi–Chatham Plateau has been subducting beneath New Zealand for at least 2.9 Myr, thus explaining the high uplift rates observed across eastern North Island.  相似文献   

9.
We present a stepwise inversion procedure to assess the focal depth and model earthquake source complexity of seven moderate-sized earthquakes  (6.2 > M w > 5.1)  that occurred in the Afar depression and the surrounding region. The Afar depression is a region of highly extended and intruded lithosphere, and zones of incipient seafloor spreading. A time-domain inversion of full moment tensor was performed to model direct P and SH waves of teleseismic data. Waveform inversion of the selected events estimated focal depths in the range of 17–22 km, deeper than previously published results. This suggests that the brittle–ductile transition zone beneath parts of the Afar depression extends more than 22 km. The effect of near-source velocity structure on the moment tensor elements was also investigated and was found to respond little to the models considered. Synthetic tests indicate that the size of the estimated, non-physical, non-isotropic source component is rather sensitive to incorrect depth estimation. The dominant double couple part of the moment tensor solutions for most of the events indicates that their occurrence is mainly due to shearing. Parameters associated with source directivity (rupture velocity and azimuth) were also investigated. Re-evaluation of the analysed events shows predominantly normal faulting consistent with the relative plate motions in the region.  相似文献   

10.
We invert differential SdS-SS traveltime residuals measured from stacked waveforms and finite-frequency sensitivity kernels for topography on the 410- and 660-km discontinuities. This approach yields higher resolution images of transition zone thickness than previous stacking methods, which simply average/smooth over topographic features. Apparent structure measured using simple stacking is highly dependent upon the bin size of each stack. By inverting for discontinuity topography with a variety of bin sizes, we can more accurately calculate the true structure. The inverted transition zone model is similar to simple stack models with an average thickness of 242 km, but the lateral variations in thickness are larger in amplitude and smaller in scale. Fast seismic velocities in 3-D mantle models such as SB4L18 correlate with areas of thicker transition zone. The elongated curvilinear regions of thickened transition zone that occur near subduction zones are narrow and high amplitude, which suggests relatively little lateral spreading and warming of subducted lithosphere within the transition zone. The anomalously thin transition zone regions are laterally narrow, and not broadly continuous. If these variations in transition zone thickness are interpreted as thermal in nature, then this model suggests significant temperature variations on small lateral scales.  相似文献   

11.
The Ethiopian side of central Afar was struck in August 1989 by the largest seismic sequence (three 6.1 ≤ M s ≤ 6.3 events, 15 with M s or m b ≥ 5.0) since that of Serdo in 1969. Using the Djibouti seismological network, we relocated 297 of the events of that sequence. As most of the large events took place outside the network, we assessed the accuracy and stability of earthquake relocations by using three different velocity models and two relocation codes to try to relate individual shocks to distinct faults and surface breaks. A majority of the events apparently occurred underneath the floor of the Dôbi graben, an area about 45  km long and 15  km wide, rupturing boundary and inner floor faults, in agreement with the surface cracks and scarps that we mapped in the area. The relocation shows that the principal events propagated about 50  km northwestwards along the graben in the first 40  hr. A day and a half after the beginning of the sequence, smaller events ( M ≤ 4) started to propagate more than 55  km eastwards, towards Asal Lake. Using two three-component stations installed near the Ethiopian border, we could determine reliable depths for 21 events. The depths are compatible with a seismogenic crust about 14  km thick in the Dôbi and Hanle graben area. Although the Dôbi sequence ruptured about 50  km of the fault array extending from Serdo to Asal, where the regional stress was released by earthquakes in 1969 and 1978, respectively, a seismic gap about 50  km long still subsists along the northern part of the Gaggade region (Der'êla half-graben).  相似文献   

12.
We have analysed the fundamental mode of Love and Rayleigh waves generated by 12 earthquakes located in the mid-Atlantic ridge and Jan Mayen fracture zone. Using the multiple filter analysis technique, we isolated the Rayleigh and Love wave group velocities for periods between 10 and 50  s. The surface wave propagation paths were divided into five groups, and average group velocities calculated for each group. The average group velocities were inverted and produced shear wave velocity models that correspond to a quasi-continental oceanic structure in the Greenland–Norwegian Sea region. Although resolution is poor at shallow depth, we obtained crustal thickness values of about 18  km in the Norwegian Sea area and 9  km in the region between Svalbard and Iceland. The abnormally thick crust in the Norwegian Sea area is ascribed to magmatic underplating and the thermal blanketing effect of sedimentary layers. Maximum crustal shear velocities vary between 3.5 and 3.9  km  s−1 for most paths. An average lithospheric thickness of 60  km was observed, which is lower than expected for oceanic-type structure of similar age. We also observed low shear wave velocities in the lower crust and upper mantle. We suggest that high heat flow extending to depths of about 30  km beneath the surface can account for the thin lithosphere and observed low velocities. Anisotropy coefficients of 1–5 per cent in the shallow layers and >7 per cent in the upper mantle point to the existence of polarization anisotropy in the region.  相似文献   

13.
We image the Hikurangi subduction zone using receiver functions derived from teleseismic earthquakes. Migrated receiver functions show a northwest dipping low shear wave feature down to 60 km depth, which we associate with the crust of the subducted Pacific Plate. Receiver functions (RF) at several stations also show a pair of negative and positive polarity phases with associated conversion depths of ∼20–26 km, where the subducted Pacific Plate is at a depth of ∼40–50 km beneath the overlying Australian Plate. RF inversion solutions model these phases with a thin low S -wave velocity zone less than 4 km thick, and an S -wave velocity contrast of more than ∼0.5 km s−1 with the overlying crust. We interpret this phase pair as representing fluids near the base of the lower crust of the Australian Plate, directly overlying the forearc mantle wedge.  相似文献   

14.
Rayleigh wave phase velocity maps in southern Africa are obtained at periods from 6 to 40 s using seismic ambient noise tomography applied to data from the Southern Africa Seismic Experiment (SASE) deployed between 1997 and 1999. These phase velocity maps are combined with those from 45 to 143 s period which were determined previously using a two-plane-wave method by Li & Burke. In the period range of overlap (25–40 s), the ambient noise and two-plane-wave methods yield similar phase velocity maps. Dispersion curves from 6 to 143 s period were used to estimate the 3-D shear wave structure of the crust and uppermost mantle on an 1°× 1° grid beneath southern Africa to a depth of about 100 km. Average shear wave velocity in the crust is found to vary from 3.6 km s–1 at 0–10 km depths to 3.86 km s–1 from 20 to 40 km, and velocity anomalies in these layers correlate with known tectonic features. Shear wave velocity in the lower crust is on average low in the Kaapvaal and Zimbabwe cratons and higher in the surrounding Proterozoic terranes, such as the Limpopo and the Namaqua-Natal belts, which suggests that the lower crust underlying the Archean cratons is probably less mafic than beneath the Proterozoic terranes. Crustal thickness estimates agree well with a previous receiver function study of Nair et al. . Archean crust is relatively thin and light and underlain by a fast uppermost mantle, whereas the Proterozoic crust is thick and dense with a slower underlying mantle. These observations are consistent with the southern African Archean cratons having been formed by the accretion of island arcs with the convective removal of the dense lower crust, if the foundering process became less vigorous in arc environments during the Proterozoic.  相似文献   

15.
Absolute S-velocity estimation from receiver functions   总被引:2,自引:0,他引:2  
We present a novel method to recover absolute S velocities from receiver functions.
For a homogeneous half-space the S velocity can be calculated from the horizontal slowness and the angle of surface particle motion for an incident P wave. Generally, the calculated S velocity is an apparent half-space value which depends on model inhomogeneity and P -waveform. We therefore, suggest to calculate such apparent half-space S velocities from low-pass filtered (smoothed) receiver functions using a suite of filter-parameters, T . The use of receiver functions neutralize the influence of the P -waveform, and the successive low-pass filterings emphasize the variation of S velocity with depth.
We apply this   V S ,app.( T )  technique to teleseismic data from three stations: FUR, BFO and SUM, situated on thick sediments, bedrock and the Greenland ice cap, respectively. The observed   V S ,app.( T )  curves indicate the absolute S velocities from the near surface to the uppermost mantle beneath each station, clearly revealing the different geological environments. Application of linearized, iterative inversion quantify these observations into   V S ( z )  models, practically independent of the S -velocity starting model. The obtained models show high consistency with independent geoscientific results. These cases provide also a general validation of the   V S ,app.( T )  method.
We propose the computation of   V S ,app.( T )  curves for individual three-component broad-band stations, both for direct indication of the S velocities and for inverse modelling.  相似文献   

16.
Crust and upper mantle structure of the central Iberian Meseta (Spain)   总被引:2,自引:0,他引:2  
Summary. Quarry blasts recorded along three lines on the central Iberian Meseta are used in an attempt to interpret the crustal structure. The results of the interpretation of the data, together with published surface wave and earthquake data, suggest a layered structure of the crust having the following features: the basement, in some areas covered by up to 4 km of sediments, has a P -velocity of 6.1 km s−1; a low-velocity layer, between 7 and 11 km depth, seems to exist on the basis of both P and S interpretation of seismic data; a thick middle crust of 12 km has a P -velocity of 6.4 km s−1 and overlies a lower crust with a mean P -velocity of 6.9 km s−1 and a possible slight negative gradient; the mean v p/ v s ratio for the crust is about 1.75; the Moho is reached at about 31 km depth and consists of a transition zone at least 1.5 km thick. The P -velocity of the upper mantle is close to 8.1 km s−1 and the S -velocity about 4.5 km s−1, which gives a v p /v s ratio of 1.8 for the uppermost mantle. A tentative petrological interpretation of the velocities and composition of the layers is given.  相似文献   

17.
Summary. Four seismic refraction lines, three of which had shots every 250 m, were shot across, along and parallel to the median valley of the Mid-Atlantic Ridge at 37° N. A method has been developed for calculating the effect on the travel times of the rough sea-floor relief beneath the profiles and has been used to correct all the travel times for this effect. Most arrivals were from a main refractor of apparent velocity 5·4 to 6·3 km s−1; only beyond 35 km were faster arrivals observed from an 8·09 ± 36 km s−1 refractor. The main refractor corresponds in depth, at least approximately, to the top of Layer 3 of the ocean basins but its velocity is significantly less than normal for Layer 3, perhaps due to dip. A study of time residuals along two profiles across the median valley indicates the presence of a 2 to 3 km wide low velocity zone (about 3·2 km s−1) beneath the median valley floor. This zone extends over the upper 2·5 km of the crust and is believed to represent a zone of intrusion through which magma passes on its way to the sea floor.  相似文献   

18.
Upper mantle shear structure of North America   总被引:5,自引:0,他引:5  
Summary. The waveforms and travel times of S and SS phases in the range 10°–60° have been used to derive upper mantle shear velocity structures for two distinct tectonic provinces in North America. Data from earthquakes on the East Pacific Rise recorded at stations in western North America were used to derive a tectonic upper mantle model. Events on the north-west coast of North America and earthquakes off the coast of Greenland provided the data to investigate the upper mantle under the Canadian shield. All branches from the triplications due to velocity jumps near 400 and 660 km were observed in both areas. Using synthetic seismograms to model these observations placed tight constraints on heterogeneity in the upper mantle and on the details of its structure. SS–S travel-time differences of 30 s along with consistent differences in waveforms between the two data sets require substantial heterogeneity to at least 350 km depth. Velocities in the upper 170 km of the shield are about 10 per cent higher than in the tectonic area. At 250 km depth the shield velocities are still greater by about 4.5 per cent and they gradually merge near 400 km. Below 400 km no evidence for heterogeneity was found. The two models both have first-order discontinuities of 4.5 per cent at 405 km and 7.5 per cent at 695 km. Both models also have lids with lower velocities beneath. In the western model the lid is very thin and of relatively low velocity. In the shield the lid is 170 km thick with very high elocity (4.78 km s-1); below it the velocity decreases to about 4.65 km s-1. Aside from these features the models are relatively smooth, the major difference between them being a larger gradient in the tectonic region from 200 to 400 km.  相似文献   

19.
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.  相似文献   

20.
The Narmada zone in central India is a zone of weakness that separates the region of Vindhyan (Meso-Neoproterozoic) deposition to the north from Gondwana (Permo-Carboniferous–lower Cretaceous) deposits to the south. The reinterpretation of analogue seismic refraction data, acquired during the early 1980s, using 2-D ray-tracing techniques reveals a basement (velocity 5.8–6.0 km s−1 ) topography suggesting that the Narmada zone, bounded by the Narmada North and Narmada South faults is a region of basement uplift. A layer of anomalously high velocity (6.5–6.7 km s−1 ) at depths between 1.5 and 9.0 km appears to be present in the entire region. Within the Narmada zone this layer occurs at shallower depths than outside the Narmada zone. At two places within the Narmada zone this layer is at a depth of about 1.5 km. This layer cannot be considered as the top of the lower crust because in this case it should have produced large positive gravity anomalies at the shallowest parts. Instead, these parts correspond to Bouguer gravity lows. Furthermore, lower crust at such shallow depths has not been reported from any other part of the Indian shield. Therefore, this layer is likely to represent the top of a high-velocity mafic body that has different thicknesses in different places.  相似文献   

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

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