首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 109 毫秒
1.
Glaciation and deglaciation in Fennoscandia during the last glacial cycles has significantly perturbed the Earth's equilibrium figure. Changes in the Earth's solid and geoidal surfaces due to external and internal mass redistributions are recorded in sequences of ancient coastlines, now either submerged or uplifted, and are still visible in observations of present‐day motions of the surface and glacially induced anomalies in the Earth's gravitational field. These observations become increasingly sophisticated with the availability of GPS measurements and new satellite gravity missions.
Observational evidence of the mass changes is widely used to constrain the radial viscosity structure of the Earth's mantle. However, lateral changes in earth model properties are usually not taken into account, as most global models of glacial isostatic adjustment assume radial symmetry for the earth model. This simplifying assumption contrasts with seismological evidence of significant lateral variations in the Earth's crust and upper mantle throughout the Fennoscandian region.
We compare predictions of glacial isostatic adjustment based on an ice model over the Fennoscandian region for the last glacial cycle for both radially symmetric and fully 3‐D earth models. Our results clearly reveal the importance of lateral variations in lithospheric thickness and asthenospheric viscosity for glacially induced model predictions. Relative sea‐level predictions can differ by up to 10–20 m, uplift rate predictions by 1–3 mm yr−1 and free‐air gravity anomaly predictions by 2–4 mGal when a realistic 3‐D earth structure as proposed by seismic modelling is taken into account.  相似文献   

2.
For more than 30 years, Sauramo's (1958) shoreline diagram of the Fennoscandian uplift has been used in geophysical studies for estimates of the glacial-isostatic decay spectrum in order to infer from it the viscosity stratification in the Earth's mantle below Fennoscandia. The intent of the present note is to point out that more recent geological studies suggest that Sauramo's shoreline diagram is an incorrect representation of the Fennoscandian uplift. Geophysical interpretations based on the diagram may therefore require revision.  相似文献   

3.
A wide-angle seismic profile across the western peninsulas of SW Ireland was performed. This region corresponds to the northernmost Variscan thrust and fold deformation. The dense set of 13 shots and 109 stations along the 120  km long profile provides a detailed velocity model of the crust.
  The seismic velocity model, obtained by forward and inverse modelling, defines a five-layer crust. A sedimentary layer, 5–8  km thick, is underlain by an upper-crustal layer of variable thickness, with a base generally at a depth of 10–12  km. Two mid-crustal layers are defined, and a lower-crustal layer below 22  km. The Moho lies at a depth of 30–32  km. A low-velocity zone, which coincides with a well-defined gravity low, is observed in the central part of the region and is modelled as a Caledonian granite which intruded upper-crustal basement. The granite may have acted as a buffer to northward-directed Variscan thrusting. The Dingle–Dungarvan Line (DDL) marks a major change in sedimentary and crustal velocity and structure. It lies immediately to the north of the velocity and gravity low, and shows thickness and velocity differences in many of the underlying crustal layers and even in the Moho. This suggests a deep, pre-Variscan control of the structural development of this area. The model is compatible with thin-skinned tectonics, which terminated at the DDL and which incorporated thrusts involving the sedimentary and upper-crustal layers.  相似文献   

4.
张赤军  陆洋 《极地研究》1998,9(2):71-75
１ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅＡｎｔａｒｃｔｉｃｉｃｅｓｈｅｔｎｅａｒｌｙｏｃｕｐｉｅｓ９０％ｏｆｔｈｅｇｌｏｂａｌｏｎｅｓ,ｔｈｅｆｏｒｍａｔｉｏｎａｎｄａｂｌａｔｉｏｎｏｆｗｈｉｃｈｈａｖｅａｓｔｒｏｎｇｉｍｐａｃｔｏｎｔｈｅｇｌｏｂａｌｇｅｏ...  相似文献   

5.
The southern South African continental margin documents a complex margin system that has undergone both continental rifting and transform processes in a manner that its present‐day architecture and geodynamic evolution can only be better understood through the application of a multidisciplinary and multi‐scale geo‐modelling procedure. In this study, we focus on the proximal section of the larger Bredasdorp sub‐basin (the westernmost of the five southern South African offshore Mesozoic sub‐basins), which is hereto referred as the Western Bredasdorp Basin. Integration of 1200 km of 2D seismic‐reflection profiles, well‐logs and cores yields a consistent 3D structural model of the Upper Jurassic‐Cenozoic sedimentary megasequence comprising six stratigraphic layers that represent the syn‐rift to post‐rift successions with geometric information and lithology‐depth‐dependent properties (porosities and densities). We subsequently applied a combined approach based on Airy's isostatic concept and 3D gravity modelling to predict the depth to the crust‐mantle boundary (Moho) as well as the density structure of the deep crust. The best‐fit 3D model with the measured gravity field is only achievable by considering a heterogeneous deep crustal domain, consisting of an uppermost less dense prerift meta‐sedimentary layer [ρ = 2600 kg m?3] with a series of structural domains. To reproduce the observed density variations for the Upper Cenomanian–Cenozoic sequence, our model predicts a cumulative eroded thickness of ca. 800–1200 m of Tertiary sediments, which may be related to the Late Miocene margin uplift. Analyses of the key features of the first crust‐scale 3D model of the basin, ranging from thickness distribution pattern, Moho shallowing trend, sub‐crustal thinning to shallow and deep crustal extensional regimes, suggest that basin initiation is typical of a mantle involvement deep‐seated pull‐apart setting that is associated with the development of the Agulhas‐Falkland dextral shear zone, and that the system is not in isostatic equilibrium at present day due to a mass excess in the eastern domain of the basin that may be linked to a compensating rise of the asthenospheric mantle during crustal extension. Further corroborating the strike‐slip setting is the variations of sedimentation rates through time. The estimated syn‐rift sedimentation rates are three to four times higher than the post‐rift sedimentation, thereby indicating that a rather fast and short‐lived subsidence during the syn‐rift phase is succeeded by a significantly poor passive margin development in the post‐rift phase. Moreover, the derived lithospheric stretching factors [β = 1.5–1.75] for the main basin axis do not conform to the weak post‐rift subsidence. This therefore suggests that a differential thinning of the crust and the mantle‐lithosphere typical for strike‐slip basins, rather than the classical uniform stretching model, may be applicable to the Western Bredasdorp Basin.  相似文献   

6.
About 50 000 P and S arrival times and 25 000 values of t * recorded at seismic arrays operated in the Central Andes between 20°S and 25°S in the time period from 1994 to 1997 have been used for locating more than 1500 deep and crustal earthquakes and creating 3-D P , S velocity and Qp models. The study volume in the reference model is subdivided into three domains: slab, continental crust and mantle wedge. A starting velocity distribution in each domain is set from a priori information: in the crust it is based on the controlled sources seismic studies; in slab and mantle wedge it is defined using relations between P and S velocities, temperature and composition given by mineral physics. Each iteration of tomographic inversion consists of the following steps: (1) absolute location of sources in 3-D velocity model using P and S arrival times; (2) double-difference relocation of the sources and (3) simultaneous determination of P and S velocity anomalies, P and S station corrections and source parameters by inverting one matrix. Velocity parameters are computed in a mesh with the density of nodes proportional to the ray density with double-sided nodes at the domain boundaries. The next iteration is repeated with the updated velocity model and source parameters obtained at the previous step. Different tests aimed at checking the reliability of the obtained velocity models are presented. In addition, we present the results of inversion for Vp and Vp/Vs parameters, which appear to be practically equivalent to Vp and Vs inversion. A separate inversion for Qp has been performed using the ray paths and source locations in the final velocity model. The resulting Vp , Vs and Qp distributions show complicated, essentially 3-D structure in the lithosphere and asthenosphere. P and S velocities appear to be well correlated, suggesting the important role of variations of composition, temperature, water content and degree of partial melting.  相似文献   

7.
A palaeomagnetic investigation has been carried out of rocks from the eastern part of the Voronezh Massif, which constitutes, together with the Ukrainian Shield, the Sarmatian segment in the southern part of the East European Craton. The samples were collected in a quarry close to the town of Pavlovsk (50.4°N, 40.1°E), where a syenitic-granitic body intrudes Archaean units. U–Pb (zircon) dating has yielded an age of 2080  Ma for the intrusion.
  Two characteristic magnetic components, A and B, were isolated by thermal and alternating-field demagnetization. Component A was obtained from granites and quartz syenites (11 samples) and has a mean direction of D = 229°, I = 28°, and a pole position at 12°N, 172°E. This pole is close to a contemporary mean pole (9°N, 187°E) for the Ukrainian Shield, which implies that the Voronezh Massif and the Shield constituted a single entity at 2.06  Ga. These poles differ from contemporaneous poles of the Fennoscandian Shield, indicating that the relative positions of the two shields were different from their present configuration about 2100  Myr ago.
  A component B, isolated only in quartz monzonites (five samples), has a mean direction D = 144°, I = 49°, and a pole position at 4°N, 251°E, which is close to late Sveconorwegian (approximately 900  Ma) poles for Baltica. This suggests that the East European Craton was consolidated some time between 2080 and 900  Ma. Comparison with other palaeomagnetic data permit us to narrow this time span to 1770–1340  Ma.  相似文献   

8.
Summary. The analysis of data of seismic crustal studies in the USSR, obtained from waves propagating at different azimuths, reveals considerable horizontal and vertical inhomogeneity of the crust. Against this background it is difficult to predict what kind of velocity anisotropy can be expected in the continental crust. The rare cases of disagreement in velocities on intersecting profiles can be attributed both to anisotropy and to horizontal crustal inhomogeneity. There is a definite disagreement in layer velocities measured by reflected waves: fine layers in the crust and upper mantle have been found to have anomalously high velocities. The role of anisotropy in these events is not clear. The frequently observed splitting of S -wave with different polarization, however, positively implies anisotropy in the Earth's crust.  相似文献   

9.
We study the tectonic setting and lithospheric structure of the greater Barents Sea region by investigating its isostatic state and its gravity field. 3-D forward density modelling utilizing available information from seismic data and boreholes shows an apparent shift between the level of observed and modelled gravity anomalies. This difference cannot be solely explained by changes in crustal density. Furthermore, isostatic calculations show that the present crustal thickness of 35–37 km in the Eastern Barents Sea is greater than required to isostatically balance the deep basins of the area (>19 km). To isostatically compensate the missing masses from the thick crust and deep basins and to adequately explain the gravity field, high-density material (3300–3350 kg m−3) in the lithospheric mantle below the Eastern Barents Sea is needed. The distribution of mantle densities shows a regional division between the Western and Eastern Barents and Kara Seas. In addition, a band of high-densities is observed in the lower crust along the transition zone from the Eastern to Western Barents Sea. The distribution of high-density material in the crust and mantle suggests a connection to the Neoproterozoic Timanide orogen and argues against the presence of a Caledonian suture in the Eastern Barents Sea. Furthermore, the results indicate that the basins of the Western Barents Sea are mainly affected by rifting, while the Eastern Barents Sea basins are located on a stable continental platform.  相似文献   

10.
Anisotropy in multi-offset deep-crustal seismic experiments   总被引：1，自引：0，他引：1  
Modelling of deep-seismic wide-angle data commonly assumes that the Earth is heterogeneous and isotropic. It is important to know the magnitudes of errors that may be introduced by isotropic-based wide-angle models when the Earth is anisotropic. It is equally important to find ways of detecting anisotropy and determining its properties.
  This paper explores the errors introduced by interpreting anisotropic seismic data with isotropic models. Errors in P -wave reflector depths are dependent on the magnitude of the velocity anisotropy and the direction of the fast axis. The interpreted, isotropic, model velocity function is found to correspond closely to the horizontal velocity of the anisotropic medium. An additional observed parameter is the time mismatch , which we define to be the difference between the vertical two-way traveltime to a reflector and the time-converted wide-angle position of the reflector. The magnitude of the time mismatch is typically <1.0  s (when the whole crust is anisotropic) and is found to be closely related to the magnitude and sign of the anisotropic anellipticity. The relationships are extendible to more complicated models, including those with vertical velocity gradients, crustal zonation, and lower symmetry orders.
  A time mismatch may be symptomatic of the presence of anisotropy. We illustrate the observation of a time mismatch for a real multi-offset seismic data set collected north of Scotland and discuss the implications for crustal anisotropy in that region.  相似文献   

11.
b
As a supplement to seismic profiling surveys, crustal thicknesses have been estimated for 11 Fennoscandian seismograph stations equipped with three-component long period instruments, using the so-called spectral ratio technique of Phinney. The largest Moho depths, of the order of 45 km, were found for stations located in the north-east areas of Norway and Sweden and in Finland, with a local maximum in the Bothnian Bay. The coastal area of south-east Norway and Zealand, Denmark exhibit crustal thicknesses in the range 28–33 km. The agreement between our results and those obtained by conventional refraction profiling is good, when this comparison is restricted to profiles of lengths 300 km or more, and when the associated crustal thickness estimate is averaged over the central parts of the profiles in question. Also, a comparison between our results and other available geophysical information gives that the oldest tectonic provinces of the Baltic Shield also are characterized by relatively modest heat flow, and exhibit the greatest crustal thicknesses. Post-glacial uplift data and large wavelength free air gravity data appear to be uncorrelated with crustal thickness. The same partly applies to Bouguer gravity anomalies, thus implying that the isostatic compensation mechanism in Fennoscandia is of both Airy and Pratt type.  相似文献   

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

13.
On crustal corrections in surface wave tomography   总被引：1，自引：0，他引：1  
Mantle models from surface waves rely on good crustal corrections. We investigated how far ray theoretical and finite frequency approximations can predict crustal corrections for fundamental mode surface waves. Using a spectral element method, we calculated synthetic seismograms in transversely isotropic PREM and in the 3-D crustal model Crust2.0 on top of PREM, and measured the corresponding time-shifts as a function of period. We then applied phase corrections to the PREM seismograms using ray theory and finite frequency theory with exact local phase velocity perturbations from Crust2.0 and looked at the residual time-shifts. After crustal corrections, residuals fall within the uncertainty of measured phase velocities for periods longer than 60 and 80 s for Rayleigh and Love waves, respectively. Rayleigh and Love waves are affected in a highly non-linear way by the crustal type. Oceanic crust affects Love waves stronger, while Rayleigh waves change most in continental crust. As a consequence, we find that the imperfect crustal corrections could have a large impact on our inferences of radial anisotropy. If we want to map anisotropy correctly, we should invert simultaneously for mantle and crust. The latter can only be achieved by using perturbation theory from a good 3-D starting model, or implementing full non-linearity from a 1-D starting model.  相似文献   

14.
The determination of the lower edges of magnetized bodies in the Earth's crust is a complex geophysical problem, although these values can be estimated by using geothermal data. An analysis of the temperature regime and location of the lower edges of magnetized bodies has been carried out for the geosynclinal region of the southern Caucasus and the area joining the ancient platform with the Arabian Shield in Israel. Geothermal calculations for Israel have been performed for three models of the thermal regime for the Earth's crust and upper mantle. The process of ultrabasic rock serpentinization is accompanied by the transformation of iron suboxide to iron oxide. Both these processes run under identical thermodynamic conditions within an average temperature interval of 200°-400°C. The Curie surface controls the position of lower edges only in fault zones where oxidation conditions hold up to great depths.  相似文献   

15.
The inverse tomography method has been used to study the P - and S -waves velocity structure of the crust and upper mantle underneath Iran. The method, based on the principle of source–receiver reciprocity, allows for tomographic studies of regions with sparse distribution of seismic stations if the region has sufficient seismicity. The arrival times of body waves from earthquakes in the study area as reported in the ISC catalogue (1964–1996) at all available epicentral distances are used for calculation of residual arrival times. Prior to inversion we have relocated hypocentres based on a 1-D spherical earth's model taking into account variable crustal thickness and surface topography. During the inversion seismic sources are further relocated simultaneously with the calculation of velocity perturbations. With a series of synthetic tests we demonstrate the power of the algorithm and the data to reconstruct introduced anomalies using the ray paths of the real data set and taking into account the measurement errors and outliers. The velocity anomalies show that the crust and upper mantle beneath the Iranian Plateau comprises a low velocity domain between the Arabian Plate and the Caspian Block. This is in agreement with global tomographic models, and also tectonic models, in which active Iranian plateau is trapped between the stable Turan plate in the north and the Arabian shield in the south. Our results show clear evidence of the mainly aseismic subduction of the oceanic crust of the Oman Sea underneath the Iranian Plateau. However, along the Zagros suture zone, the subduction pattern is more complex than at Makran where the collision of the two plates is highly seismic.  相似文献   

16.
The dispersive properties of surface waves are used to infer earth structure in the Eastern Mediterranean region. Using group velocity maps for Rayleigh and Love waves from 7 to 100 s, we invert for the best 1-D crust and upper-mantle structure at a regular series of points. Assembling the results produces a 3-D lithospheric model, along with corresponding maps of sediment and crustal thickness. A comparison of our results to other studies finds the uncertainties of the Moho estimates to be about 5 km. We find thick sediments beneath most of the Eastern Mediterranean basin, in the Hellenic subduction zone and the Cyprus arc. The Ionian Sea is more characteristic of oceanic crust than the rest of the Eastern Mediterranean region as demonstrated, in particular, by the crustal thickness. We also find significant crustal thinning in the Aegean Sea portion of the backarc, particularly towards the south. Notably slower S -wave velocities are found in the upper mantle, especially in the northern Red Sea and Dead Sea Rift, central Turkey, and along the subduction zone. The low velocities in the upper mantle that span from North Africa to Crete, in the Libyan Sea, might be an indication of serpentinized mantle from the subducting African lithosphere. We also find evidence of a strong reverse correlation between sediment and crustal thickness which, while previously demonstrated for extensional regions, also seems applicable for this convergence zone.  相似文献   

17.
When interpreting electromagnetic fields observed at the Earth's surface in a realistic geophysical environment it is often necessary to pay special attention to the effects caused by inhomogeneities of the subsurface sedimentary and/or water layer and by inhomogeneities of the Earth's crust. The inhomogeneities of the Earth's crust are expected to be especially important when the electromagnetic field is generated by a source located in a magma chamber of a volcano. The simulation of such effects can be carried out using generalized thin-sheet models, which were independently introduced by Dmitriev (1969 ) and Ranganayaki & Madden (1980 ). In the first part of the paper, a system of integral equations is derived for the horizontal current that flows in the subsurface inhomogeneous conductive layer and for the vertical current crossing the inhomogeneous resistive layer representing the Earth's mantle. The terms relating to the finite thickness of the laterally inhomogeneous part of the model are retained in the equations. This only marginally complicates the equations, whilst allowing for a significant expansion of the approximation limits.
  The system of integral equations is solved using the iterative dissipative method developed by the authors in the period from 1978 to 1988. The method can be applied to the simulation of the electromagnetic field in an arbitrary inhomogeneous medium that dissipates the electromagnetic energy. When considered on a finite numerical grid, the integral equations are reduced to a system of linear equations that possess the same contraction properties as the original equations. As a result, the rate at which the iterative-perturbation sequence converges to the solution remains independent of the numerical grid used for the calculations. In contrast to previous publications on the method, aspects of the algorithm implementation that guarantee its effectiveness and robustness are discussed here.  相似文献   

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

19.
Summary. Numerical convection models are presented in which plates are simulated by imposing piecewise constant horizontal velocities on the upper boundary. A 4 × 1 box of constant viscosity fluid and two-dimensional (2-D) flow is assumed. Four heating modes are compared: the four combinations of internal or bottom heating and prescribed bottom temperature or heat flux. The case with internal heating and an isothermal base is relevant to lower mantle or whole mantle convection, and it yields a lower thermal boundary layer which is laterally variable and can be locally reversed, corresponding to heat flowing back into the core locally. When scaled to the whole mantle, the surface deflections and gravity and geoid perturbations calculated from the models are comparable to those observed at the Earth's surface. For models with migrating ridges and trenches, the flow structure lags well behind the changing surface 'plate'configurations. This may help to explain the poor correlation between the main geoid features and plate boundaries. Trench migration substantially affects the dip of the cool descending fluid because of induced horizontal shear in the vicinity of the trench. Such shear is small for whole mantle convection, but is large for upper mantle convection, and would probably result in the Tonga Benioff zone dipping to the SE, opposite to the observed dip, for the case of upper mantle convection.  相似文献   

20.
Joint inversion of receiver function and surface wave dispersion observations   总被引：16，自引：0，他引：16  
We implement a method to invert jointly teleseismic P wave receiver functions and surface wave group and phase velocities for a mutually consistent estimate of earth structure. Receiver functions are primarily sensitive to shear wave velocity contrasts and vertical traveltimes, and surface wave dispersion measurements are sensitive to vertical shear wave velocity averages. Their combination may bridge resolution gaps associated with each individual data set. We formulate a linearized shear velocity inversion that is solved using a damped leastsquares scheme that incorporates a priori smoothness constraints for velocities in adjacent layers. The data sets are equalized for the number of data points and physical units in the inversion process. The combination of information produces a relatively simple model with a minimal number of sharp velocity contrasts. We illustrate the approach using noisefree and realistic noise simulations and conclude with an inversion of observations from the Saudi Arabian Shield. Inversion results for station SODA, located in the Arabian Shield, include a crust with a sharp gradient near the surface (shear velocity changing from 1.8 to 3.5 km s1 in 3 km) underlain by a 5kmthick layer with a shear velocity of 3.5 km s1 and a 27kmthick layer with a shear velocity of 3.8 km s1, and an upper mantle with an average shear velocity of 4.7 km s1. The crustmantle transition has a significant gradient, with velocity values varying from 3.8 to 4.7 km s1 between 35 and 40 km depth. Our results are compatible with independent inversions for crustal structure using refraction data.  相似文献   

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

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