首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 187 毫秒
1.
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.  相似文献   

2.
Magnetotelluric data from the backarc of the Central Andes in NW Argentinawere re-examined by employing impedance tensor decomposition and 2-D inversion and modelling techniques. The data in the period range of 50–15 000 s were collected on a profile of 220 km length reaching from the Eastern Cordillera across the Santa Barbara System to the Andean foreland of the Argentinean Chaco.
After a dimensionality analysis, data from most sites were treated as regional 2-D. The exception was the eastern section of the profile, where the magnetotelluric transfer functions for periods ≤ 1000 s reflect a 3-D earth. Application of two tensor decomposition schemes yielded a regional strike direction of N–S, which is the azimuth of the Central Andean mountain chains. Several 2-D models were obtained by pseudo- and full 2-D Occam inversion schemes. Special emphasis was placed on the inversion of phase data to reduce the influence of static shifts in the apparent resistivity data. The smooth inversion models all show a good conductor at depth. A final model was then calculated using a finite element forward algorithm.
The most prominent feature of the resulting model is a conductor which rises from depths of 180 km below the Chaco region to 80 km beneath the Santa Barbara System and the Eastern Cordillera. Its interpretation as a rise of the electrical asthenosphere is supported by seismic attenuation studies. Magnetotelluric results, surface heat-flow distribution in the area, and the electrical properties of crustal and mantle rocks suggest that the upper mantle is predominantly ductile beneath the Eastern Cordillera and the western Santa Barbara System. This generally agrees with anelastic seismic attenuation models of the area and is useful in discriminating between models of Q quality factor distribution.  相似文献   

3.
The conductivity structure of the Earth's mantle was estimated using the induction method down to 2100  km depth for the Europe–Asia region. For this purpose, the responses obtained at seven geomagnetic observatories (IRT, KIV, MOS, NVS, HLP, WIT and NGK) were analysed, together with reliable published results for 11  yr variations. 1-D spherical modelling has shown that, beneath the mid-mantle conductive layer (600–800  km), the conductivity increases slowly from about 1  S  m−1 at 1000  km depth to 10  S  m−1 at 1900  km, while further down (1900–2100  km) this increase is faster. Published models of the lower mantle conductivity obtained using the secular, 30–60  yr variations were also considered, in order to estimate the conductivity at depths down to the core. The new regional model of the lower mantle conductivity does not contradict most modern geoelectrical sounding results. This model supports the idea that the mantle base, situated below 2100  km depth, has a very high conductivity.  相似文献   

4.
Earthquake arrival time data from a 36-station deployment of portable seismographs on the Raukumara Peninsula have been used to determine the 3-D Vp and Vp/Vs structure of this region of shallow subduction. A series of inversions have been performed, starting with an inversion for 1-D structure, then 2-D, and finally 3-D. This procedure ensures a smooth regional model in places of low resolution. The subducted plate is imaged as a northwest-dipping feature, with Vp consistently greater than 8.5  km  s−1 in the uppermost mantle of the plate. Structure in the overlying plate changes significantly along strike. In the northeast, there is an extensive low-velocity zone in the lower crust underlying the most rapidly rising part of the Raukumara Range. It is bounded on its arcward side by an upwarp of high velocity. A viable explanation for the low-velocity zone is that it represents an accumulation of underplated subducted sediment, while serpentinization of the uppermost mantle may be responsible for the adjacent high-velocity region. The low-velocity zone decreases and the adjacent high-velocity region is less extensive in the southwest. This change is interpreted to be related to a change in the thickness of the crust of the overlying plate. In the northeast the crust is thinner, and subducted sediment ponds against relatively strong uppermost mantle, while in the southwest the crust is thicker, and the relatively weak lower crust allows sediment subduction to greater depths. A narrow zone of high Vp/Vs parallels the shallow part of the plate interface. This suggests elevated fluid pressures, with the distribution of earthquakes about this zone further suggesting that these pressures may be close to lithostatic. The plate interface at 20  km depth beneath the Raukumara Peninsula may thus be a closed system for fluid flow, similar to that seen at much shallower depths in other subduction décollements.  相似文献   

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

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

7.
The Southern Andes differ significantly from the Central Andes with respect to topography and crustal structures and are, from a geophysical point of view, less well known. In order to provide insight into the along-strike segmentation of the Andean mountain belt, an integrated 3-D density model was developed for the area between latitudes 36°S and 42°S. The model is based on geophysical and geological data acquired in the region over the past years and was constructed using forward density modelling. In general, the gravity field of the South American margin is characterized by a relatively continuous positive anomaly along the coastline and the forearc region, and by negative anomalies along the trench and the volcanic arc. However, in the forearc region of the central part of the study area, located just to the south of the epicentre of the largest ever recorded earthquake (Valdivia, 1960), the trench-parallel positive anomaly is disrupted. The forearc gravity anomaly differences thus allow the study area to be divided into three segments, the northern Arauco-Lonquimay, the middle Valdivia-Liquiñe, and the southern Bahía-Mansa-Osorno segment, which are also evident in geology. In the proposed model, the observed negative gravity anomaly in the middle segment is reproduced by an approximately 5 km greater depth to the top of the slab beneath the forearc region. The depth to the slab is, however, dependent upon the density of the upper plate structures. Therefore, both the upper and lower plates and their interaction have a significant impact on the subduction-zone gravity field.  相似文献   

8.
A whole mantle SH velocity model is obtained by using a unique data set and techniques. Body and surface waveforms including major and multi-orbit phases are used as a data set and are inverted by using 3-D Born kernels. The resultant model, SH18CE, reveals the different natures of the two major upwelling systems: the strong low velocity anomalies beneath Africa extend for more than 1000 km from the core–mantle boundary (CMB), whereas those beneath the Pacific are restricted to 300–400 km from the CMB. The results also show the variable natures of stagnant slabs on the 670 discontinuity around Japan: the depths of the strongest high velocity anomalies within the stagnant slabs are different region by region, which is consistent with the detailed delay time tomography model in this area.  相似文献   

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

10.
According to the theory of isostasy, the Earth has a tendency to deform its surface in order to reach an equilibrium state. The land-uplift phenomenon in the area of the Fennoscandian Shield is thought to be a process of this kind. The geoid, as an equipotential surface of the Earth's gravity field, contains information on how much the Earth's surface departs from the equilibrium state. In order to study the isostatic process through geoidal undulations, the structural effects of the crust on the geoid have to be investigated.
  The structure of the crust of the Fennoscandian Shield has been extensively explored by means of deep seismic sounding (DSS). The data obtained from DSS are used to construct a 3-D seismic-velocity structure model of the area's crust. The velocity model is converted to a 3-D density model using the empirical relationship that holds between seismic velocities and crustal mass densities. Structural effects are then estimated from the 3-D density model.
  The structural effects computed from the crustal model show that the mass deficiency of the crust in Fennoscandia has caused a geoidal depression twice as deep as that observed from the gravimetric geoid. It proves again that the crust has been isostatically compensated by the upper mantle. In other words, an anomalously high-density upper mantle must exist beneath Fennoscandia.  相似文献   

11.
The response of a viscoelastic Earth to the melting of the Late Pleistocene ice sheets has been the subject of a number of investigations employing PREM. In PREM, a non-adiabatic density gradient (NADG) exists in the upper mantle, and to understand the implications of this model it is thus important to examine the effects of this NADG on the Earth's response to surface loads. This paper is based on the assumption that the contribution to the depth dependence of the density that is not due to self-compression is due to compositional change. This contribution is referred to as 'non-adiabatic'. We evaluate the effects of a non-adiabatic density jump (NADJ) for the 670  km discontinuity and the NADG in the upper mantle by adopting a compressible earth model with both a compositional density gradient and a density jump. Numerical calculations based on these models indicate that the magnitude of the Earth's response associated with the NADG is much smaller than that associated with the NADJ at 670  km depth. It is also confirmed that the higher modes associated with the NADJ and the NADG are much more sensitive to the existence of an elastic lithosphere than the fundamental modes associated with the density jumps at the surface and core–mantle boundary.  相似文献   

12.
It is now widely accepted that elastic properties of the continental lithosphere and the underlying sublithospheric mantle are both anisotropic and laterally heterogeneous at a range of scales. To fully exploit modern three-component broad-band array data sets requires the use of comprehensive modelling tools. In this work, we investigate the use of a wide-angle, one-way wave equation to model variations in teleseismic 3-D waveforms due to 2-D elastic heterogeneity and anisotropy. The one-way operators are derived based on a high-frequency approximation of the square-root operator and include the effects of wave propagation as well as multiple scattering. Computational cost is reduced through a number of physically motivated approximations. We present synthetic results from simple 1-D (layer over a half-space) and 2-D (subduction zone) models that are compared with reference solutions. The algorithm is then used to model data from an array of broad-band seismograph stations deployed in northwestern Canada as part of the IRIS-PASSCAL/LITHOPROBE CANOE experiment. In this region radial-component receiver functions show a clear continental Moho and the presence of crustal material dipping into the mantle at the suture of two Palaeo-Proterozoic terranes. The geometry of the suture is better defined on the transverse component where subduction is associated with a ∼10 km thick layer exhibiting strong elastic anisotropy. The modelling reproduces the main features of the receiver functions, including the effects of anisotropy, heterogeneity and finite-frequency scattering.  相似文献   

13.
20 magnetotelluric (MT) soundings were collected on the Isle of Skye, Scotland to provide a high-resolution three-dimensional (3-D) electrical resistivity model of a volcanic province within the framework of a project jointly interpreting gravity, seismic, geological and MT data. The full 3-D inversion of the MT data jointly interpreted with gravity data reveals upper crustal structure. The main features of the model are interpreted in conjunction with previous geological mapping and borehole data. Our model extends to 13 km depth, several kilometres below the top of the Lewisian basement. The top of the Lewisian basement is at approximately 7–8 km depth and the topography of its surface was controlled by Precambrian rifting, during which a 4.5 km thick sequence of Torridonian sediments was deposited. The Mesozoic sediments above, which can reach up to 2.2 km thick, have small-scale depocentres and are covered by up to 600 m of Tertiary lava flows. The interpretation of the resistivity model shows that 3-D MT inversion is an appropriate tool to image sedimentary structures beneath extrusive basalt units, where conventional seismic reflection methods may fail.  相似文献   

14.
Surface wave tomography of the Barents Sea and surrounding regions   总被引：1，自引：0，他引：1  
The goal of this study is to refine knowledge of the structure and tectonic history of the European Arctic using the combination of all available seismological surface wave data, including historical data that were not used before for this purpose. We demonstrate how the improved data coverage leads to better depth and spatial resolution of the seismological model and discovery of intriguing features of upper-mantle structure. To improve the surface wave data set in the European Arctic, we extensively searched for broad-band data from stations in the area from the beginning of the 1970s until 2005. We were able to retrieve surface wave observations from regional data archives in Norway, Finland, Denmark and Russia in addition to data from the data centres of IRIS and GEOFON. Rayleigh and Love wave group velocity measurements between 10 and 150 s period were combined with existing data provided by the University of Colorado at Boulder. This new data set was inverted for maps showing the 2-D group-velocity distribution of Love and Rayleigh waves for specific periods. Using Monte Carlo inversion, we constructed a new 3-D shear velocity model of the crust and upper mantle beneath the European Arctic which provides higher resolution and accuracy than previous models. A new crustal model of the Barents Sea and surrounding areas, published recently by a collaboration between the University of Oslo, NORSAR and the USGS, constrains the 3-D inversion of the surface wave data in the shallow lithosphere. The new 3-D model, BARMOD, reveals substantial variations in shear wave speeds in the upper mantle across the region with a nominal resolution of 1°× 1°. Of particular note are clarified images of the mantle expression of the continent-ocean transition in the Norwegian Sea and a deep, high wave speed lithospheric root beneath the Eastern Barents Sea, which presumably is the remnant of several Palaeozoic collisions.  相似文献   

15.
A new algorithm is presented for the integrated 2-D inversion of seismic traveltime and gravity data. The algorithm adopts the 'maximum likelihood' regularization scheme. We construct a 'probability density function' which includes three kinds of information: information derived from gravity measurements; information derived from the seismic traveltime inversion procedure applied to the model; and information on the physical correlation among the density and the velocity parameters. We assume a linear relation between density and velocity, which can be node-dependent; that is, we can choose different relationships for different parts of the velocity–density grid. In addition, our procedure allows us to consider a covariance matrix related to the error propagation in linking density to velocity. We use seismic data to estimate starting velocity values and the position of boundary nodes. Subsequently, the sequential integrated inversion (SII) optimizes the layer velocities and densities for our models. The procedure is applicable, as an additional step, to any type of seismic tomographic inversion.
We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently.
An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area.  相似文献   

16.
The C -response connects the magnetic vertical component and the horizontal gradient of the horizontal components of electromagnetic variations and forms the basis for deriving the conductivitydepth profile of the Earth. Time-series of daily mean values at 42 observatories typically with 50 years of data are used to estimate C -responses for periods between 1 month and 1  yr. The Z : Y method is applied, which means that the vertical component is taken locally whereas the horizontal components are used globally by expansion in a series of spherical harmonics.
In combination with results from previous analyses, the method yields consistent results for European observatories in the entire period range from a few hours to 1  yr, corresponding to penetration depths between 300 and 1800  km.
1-D conductivity models derived from these results show an increase in conductivity with depth z to about 2  S  m-1 at z =800  km, and almost constant conductivity between z =800 and z =2000  km with values of 310  S  m-1, in good agreement with laboratory measurements of mantle material. Below 2000  km the conductivity is poorly resolved. However, the best-fitting models indicate a further increase in conductivity to values between 50 and 150  S  m-1.  相似文献   

17.
On the density distribution within the Earth   总被引：1，自引：0，他引：1  
The distribution of density as a function of position within the Earth is much less well constrained than the seismic velocities. The primary information comes from the mass and moment of inertia of the Earth and this information alone requires that there be a concentration of mass towards the centre of the globe. Additional information is to be found in the frequencies of the graver normal modes of the Earth which are sensitive to density through self-gravitation effects induced in deformation.
  The present generation of density models has been constructed using linearized inversion techniques from earlier models, which ultimately relate back to models developed by Bullen and based in large part on physical arguments. A number of experiments in non-linear inversion have been conducted using the PREM reference model, with fixed velocity and attenuation, but with the density model constrained to lie within fixed bounds on both density and density gradient. A set of models is constructed from a uniform probability density within the bound and slope constraints. Each of the resultant density models is tested against the mass and moment of inertia of the Earth, and for successful models a comparison is made with observed normal mode frequencies. From the misfit properties of the ensemble of models the robustness of the density profile in different portions of the Earth can be assessed, which can help with the design of parametrization for future reference models. In both the lower mantle and the outer core it would be desirable to allow a more flexible representation than the single cubic polynomial employed in PREM.  相似文献   

18.
Inference of mantle viscosity from GRACE and relative sea level data   总被引：12，自引：0，他引：12  
Gravity Recovery And Climate Experiment (GRACE) satellite observations of secular changes in gravity near Hudson Bay, and geological measurements of relative sea level (RSL) changes over the last 10 000 yr in the same region, are used in a Monte Carlo inversion to infer-mantle viscosity structure. The GRACE secular change in gravity shows a significant positive anomaly over a broad region (>3000 km) near Hudson Bay with a maximum of ∼2.5 μGal yr−1 slightly west of Hudson Bay. The pattern of this anomaly is remarkably consistent with that predicted for postglacial rebound using the ICE-5G deglaciation history, strongly suggesting a postglacial rebound origin for the gravity change. We find that the GRACE and RSL data are insensitive to mantle viscosity below 1800 km depth, a conclusion similar to that from previous studies that used only RSL data. For a mantle with homogeneous viscosity, the GRACE and RSL data require a viscosity between  1.4 × 1021  and  2.3 × 1021  Pa s. An inversion for two mantle viscosity layers separated at a depth of 670 km, shows an ensemble of viscosity structures compatible with the data. While the lowest misfit occurs for upper- and lower-mantle viscosities of  5.3 × 1020  and  2.3 × 1021  Pa s, respectively, a weaker upper mantle may be compensated by a stronger lower mantle, such that there exist other models that also provide a reasonable fit to the data. We find that the GRACE and RSL data used in this study cannot resolve more than two layers in the upper 1800 km of the mantle.  相似文献   

19.
Upper-mantle structure under the Baltic Shield is studied using non-linear high resolution teleseismic P -phase tomography. Observed relative arrival-time residuals from 52 teleseismic earthquakes recorded by the Swedish National Seismological Network (SNSN) are inverted to delineate the structure of the upper mantle. The network consists of 47 (currently working) three-component broad-band stations located in an area about 450 km wide and 1450 km long. In order to reduce complications due to possible significant three-dimensionality of Earth structure, events chosen for this study lay close to in-line with the long-axis of the array  (±30°)  . Results indicate P -wave velocity perturbations of ±3 per cent down to at least 470 km below the network. The size of the array allows inversion for structures even at greater depths, and lateral variations of velocity at depths of up to 680 km appear to be resolved. Below the central part of the array (60°–64° N), where ray coverage is best, the data reveals a large region of relatively low velocity at depths of over about 300 km. At depths less than about 250–300 km, the models include a number of features, including an apparent slab-like structure dipping gently towards the north.  相似文献   

20.
Signature of remnant slabs in the North Pacific from P-wave tomography   总被引：1，自引：0，他引：1  
A 3-D ray-tracing technique was used in a global tomographic inversion in order to obtain tomographic images of the North Pacific. The data reported by the Geophysical Survey of Russia (1955–1997) were used together with the catalogues of the International Seismological Center (1964–1991) and the US Geological Survey National Earthquake Information Center (1991–1998), and the recompiled catalogue was reprocessed. The final data set, used for following the inversion, contained 523 430 summary ray paths. The whole of the Earth's mantle was parametrized by cells of 2° × 2° and 19 layers. The large and sparse system of observation equations was solved using an iterative LSQR algorithm.
A subhorizontal high-velocity anomaly is revealed just above the 660 km discontinuity beneath the Aleutian subduction zone. This high-velocity feature is observed at latitudes of up to ~70°N and is interpreted as a remnant of the subducted Kula plate, which disappeared through ridge subduction at about 48 Ma. A further positive velocity perturbation feature can be identified beneath the Chukotka peninsula and Okhotsk Sea, extending from ~300 to ~660 km depth and then either extending further down to ~800 km (Chukotka) or deflecting along the 660 km discontinuity (Okhotsk Sea). This high-velocity anomaly is interpreted as a remnant slab of the Okhotsk plate accreted to Siberia at ~55 Ma.  相似文献   

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

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