首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 390 毫秒
1.
Summary. Models of shallow, global mantle circulation due to the accretion and subduction of lithospheric plates are formulated as potential theory problems on a sphere. Subducting and accreting plate boundaries represent sources and sinks respectively for the sublithospheric flow. Solutions, which are obtained by finite difference approximations, give the instantaneous flow velocities within the asthenosphere compatible with plate boundaries and relative plate motions. Results are presented for present-day plate boundaries and relative plate motions for the case of a uniform viscosity asthenosphere and for that of a low viscosity zone at the base of the lithosphere. These results are discussed in terms of available geophysical data. Some of the implications of a shallow, mantle-wide circulation are also considered.  相似文献   

2.
Summary. A microplate is modelled as an elastic plate with two long strike-slip boundaries, lying over a Maxwell-type viscoelastic asthenosphere. The microplate is subjected to a constant and uniform shear strain rate by the opposite motions of two adjoining larger plates. After the occurrence of an earthquake at one of the microplate boundaries, the time evolution of shear stress at the other boundary is studied. It is found that stress build-up at the second boundary is delayed due to stress diffusion governed by the asthenosphere relaxation. Earthquake occurrence at this latter boundary would be delayed depending upon both the microplate width and the ratio between the Maxwell relaxation time of the asthenosphere and a characteristic time required for tectonic strain to recover rupture conditions. It turns out that the parameters which determine the occurrence of seismic activity along the microplate boundaries are more strictly constrained in the presence of a viscoelastic asthenosphere than in the case of an elastic half-pace model.  相似文献   

3.
Summary. Present-day plate motions imply that about 240 km3 of oceanic lithosphere is created by sea-floor spreading and destroyed by subduction per year. A greater volume of asthenosphere will be dragged along by plate motions. Given the fluxes generated at plate boundaries, the horizontal direction and net rate of counterflow required to maintain mass balance is determined globally by a simple analytical model. Time-dependent calculations indicate that the motions are approximately valid in the hotspot reference frame over the past 5 Myr. Under most plates, the model return flow is opposite to the lithospheric motion in the hotspot frame. The counterflow dominates the resisting stresses to plate motion, so driving force models based on plate drag alone are not valid where the directions of plate motion and counterflow differ. The most marked departure of the two directions is under the North American plate. The model counterflow directions indicate that the sources of mantle hotspots are not located within the asthenosphere. Model flux balances demonstrate exchange of material between asthenospheric reservoirs located beneath different plates. Suggestions of southward asthenospheric motion under the North Atlantic, based on physical features around Iceland and strontium isotope geochemistry, are consistent with the direction of flow predicted by the model.  相似文献   

4.
Summary. The thickening plate theory proposed by Yoshii and Parker & Oldenburg for the oceanic lithosphere is extended to include the continental lithosphere. The theory is based on the assumption that the lithosphere—asthenosphere boundary is a solidus and that as a result solidification of the top of the asthenosphere is occurring. Observational data imply that the relationship between the plate thickness and basement age for the North American continent is y = 1.7 √ t + (50 ± 10), where y (km) is the plate thickness and t (Myr) is the basement age.
The theory is tested against changes with basement age of the observed surface heat-flow and seismic estimate of plate thickness. The following conclusions are inferred:
(1) The changes both of the observed heat flow and plate thickness with basement age are explained by this theory.
(2) The surface erosion and vertical distribution of radiogenic heat sources are important factors in controlling the thickening process of the continental lithosphere.
(3) The equality of the average surface heat-flow over the oceans and over the continents is a consequence of a faster release of latent heat at the lithosphere—asthenosphere boundary under the oceans, instead of a higher heat production in the continental crust.  相似文献   

5.
Summary. Because there is secondary sea-floor spreading in the Tonga and Mariana subduction systems, the island arcs are separate plates. Horizontal forces on the two sides of the arc must balance, and the maximum force on the back-arc side can be calculated from a lithostatic ridge model. This, in combination with gravity data, allows calculation of the average shear stress in the top 100 km of the subduction shear zone. Stress in Tonga is 220±100 bar, and in the Mariana it is 165±75 bar. These low stresses are probably made possible by a fluid pore pressure almost equal to the least compressive stress.
Knowledge of stress allows approximate calculation of temperature in the shear zone by integration of a single differential equation. These temperatures are too low to activate most dehydration reactions in the subducted crust. As it approaches the volcanic line, this crust is at 150–350°C in Tonga and 150–300°C in the Mariana. Shear melting of the crust is ruled out, and conductive melting of the slab by contact with the asthenosphere meets with geochemical objections. Magmas in these systems are probably produced by partial melting of asthenosphere, triggered by a sudden release of water from the slab.  相似文献   

6.
Summary We consider a long strike-slip fault in a lithosphere modelled as an elastic slab. To the base of the slab a shear stress distribution is applied which simulates the viscous drag exerted by the asthenosphere. The resulant stress on the fault plane may directly fracture the lithosphere in its brittle upper portion; alternatively it may give rise at first to a stable aseismic sliding in the lower portion. In the latter case, stress concentration due to the deep aseismic slip is the relevant feature of the pre-seismic stress acting on the upper section of the lithosphere. The two cases are examined by use of dislocation theory and their observable effects compared. Different depths of the aseismic slip zone and the presence or absence of a uniform friction on the seismic fault are allowed for. If the model is applied to the San Andreas fault region, where a steady sliding condition actually seems to be present at shallow depth, it turns out that the slip amplitudes commonly associated with large earthquakes are consistent with average basal stress values which can be substantially lower than a few bars, a value often quoted as the steady state basal stress due to a velocity gradient in the upper asthenosphere.  相似文献   

7.
Summary . Frictional heating in upper mantle shear flows may lead to localized thermal runaway and partial melting in the asthenosphere, but only as the result of a finite-amplitude disturbance. A rigorous two-dimensional stability analysis shows that asthenospheric shear flows are stable to small-amplitude perturbations whether such flows are supercritical (shear stress decreases with increasing plate velocity) or subcritical (shear stress increases with increasing plate velocity). Disturbances which maintain a shear stress larger than the critical value for sufficiently long will lead to runaway. The response of the asthenosphere to events which do not satisfy this criterion must be determined by a non-linear analysis. Reasonable models of flow in the asthenosphere could be driven to runaway, at a superexponential growth rate, by sudden increases in shear stress of less than 10 bar. Disturbances resulting from plate collisions may maintain large enough stresses for sufficiently long times to initiate runaways, while stress changes associated with large earthquakes probably occur too rapidly to do so.  相似文献   

8.
Viscous and viscoelastic models for a subduction zone with a faulted lithosphere and internal buoyancy can self-consistently and simultaneously predict long-wavelength geoid highs over slabs, short-wavelength gravity lows over trenches, trench-forebulge morphology, and explain the high apparent strength of oceanic lithosphere in trench environments. The models use two different free-surface formulations of buoyancy-driven flows (see, for example, Part I): Lagrangian viscoelastic and pseudo-free-surface viscous formulations. The lower mantle must be stronger than the upper in order to obtain geoid highs at long wavelengths. Trenches are a simple consequence of the negative buoyancy of slabs and a large thrust fault, decoupling the overriding from underthrusting plates. The lower oceanic lithosphere must have a viscosity of less than to24 Pa s in order to be consistent with the flexural wavelength of forebulges. Forebulges are dynamically maintained by viscous flow in the lower lithosphere and mantle, and give rise to apparently stiffer oceanic lithosphere at trenches. With purely viscous models using a pseudo-free-surface formulation, we find that viscous relaxation of oceanic lithosphere, in the presence of rapid trench rollback, leads to wider and shallower back-arc basins when compared to cases without viscous relaxation. Moreover, in agreement with earlier studies, the stresses necessary to generate forebulges are small (∼ 100 bars) compared to the unrealistically high stresses needed in classic thin elastic plate models.  相似文献   

9.
In Paper I (Breuer & Wolf 1995), a preliminary interpretation of the postglacial land emergence observed at a restricted set of six locations in the Svalbard Archipelago was given. The study was based on a simple model of the Barents Sea ice sheet and suggested increases in lithosphere thickness and asthenosphere viscosity with increasing distance from the continental margin.
In the present paper, the newly developed high-resolution load model. BARENTS-2, and land-uplift observations from an extended set of 25 locations are used to study further the possibility of resolving lateral heterogeneity in the upper mantle below the northern Barents Sea. A comparison of the calculated and observed uplift values shows that the lithosphere thickness is not well resolved by the observations, although values above 110 km are most common for this parameter. In contrast to this, there are indications of a lateral variation of asthenosphere viscosity. Whereas values in the range 1018-1020Pas are inferred for locations close to the continental margin, 1020-1021 Pa s are suggested further away from the margin.
A study of the sensitivity of the values found for lithosphere thickness and asthenosphere viscosity to modifications of load model BARENTS-2 shows that such modifications can be largely accommodated by appropriate changes in lithosphere thickness, whereas the suggested lateral variation of asthenosphere viscosity is essentially unaffected. An estimate of the influence of the Fennoscandian. ice sheet leads to the conclusion that its neglect results in an underestimation of the thickness of the Barents Sea ice sheet by about 10 per cent.  相似文献   

10.
Summary. If intermediate and deep earthquakes occur in the coldest portions of the downgoing slabs of lithosphere, then different lengths of seismic zones represent different temperatures in the slabs. As the slab descends through the aesthenosphere, it warms primarily by conduction of heat through its upper and lower surfaces. Isotherms are advected downwards to distances approximately proportional both to the rate of subduction and to the square of the thickness of the lithosphere. Consequently, the lengths of seismic zones should be approximately proportional to the product of the rates and the squares of these thicknesses. As these thicknesses are approximately proportional to the square root of the age of the lithosphere, the lengths ought therefore to be approximately proportional to the product of the convergence rates times the ages. Although there is considerable scatter, observed lengths are approximately proportional to such products, and are not simply related to the rate, the age or the thickness alone. The data crudely fit the relationship: length = rate × age/10. Using this relationship, we infer that the Philippine sea and Pacific plates move slowly, if at all, with respect to one another and that the Farallon plate may have been too young to be subducted to a great distance beneath western North America in the Palaeogene. Calculations of temperatures at the depths of the deepest events suggest that these cut-off temperatures increase from about 600 ± 100°C at 200 km to 830 ± 50°C at 650 km depth, but the cut-off potential temperature is approximately constant. Assuming that the strength is a thermally activated parameter, and using the activation energy for olivine, a crude estimate of activation volume is obtained from the dependence of the cut-off temperature on depth.  相似文献   

11.
Summary. Various factors can make it difficult to explain observations of Love- and Rayleigh-wave dispersion with the same relatively simple isotropic model. These factors include systematic errors which might occur in determinations of observed group and phase velocities, lateral variations in structure along the path of travel, and the attempt to explain observations with a model comprised of only a small number of thick layers. The last of these factors is illustrated by an inversion of dispersion data in the central United States where shear-wave anisotropy had previously been invoked as one way to explain incompatible Love- and Rayleigh-wave velocities. It is shown that the data can be satisfied equally well by an isotropic model consisting of several thin layers.
In cases where the incompatibility of Love- and Rayleigh-wave data might be produced by intrinsic anisotropy, it is necessary to invert those data using an anisotropic theory rather than by separate isotropic inversions of Love and Rayleigh waves. Inversions of fundamental-mode data for a region of the Pacific, assuming anisotropic media in which the layers are transversely isotropic with a vertical axis of symmetry, lead to models which are highly non-unique. Even if the inversions solve only for shear velocities in the litho-sphere and asthenosphere it is not possible, without supplementary information, to ascertain the depth interval over which anisotropy occurs or to determine the thickness of the lithosphere or asthenosphere with much precision.  相似文献   

12.
Summary. We give the analytical formulation for calculating the transient displacement of fields produced by earthquakes in a stratified, selfgravitating, incompressible, viscoelastic earth. We have evaluated the potential of viscous creep in the asthenosphere in exciting the Chandler wobble by a four-layer model consisting of an elastic lithosphere, a two-layer Maxwell viscoelastic mantle, and an inviscid core. The seismic source is modelled as an inhomogeneous boundary condition, which involves a jump condition of the displacement fields across the fault in the lithosphere. The response fields are derived from the solution of a two-point boundary value problem, using analytical propagator matrices in the Laplace-transformed domain. Transient flows produced by post-seismic rebound are found to be confined within the asthenosphere for local viscosity values less than 1020P. The viscosity of the mantle below the low-viscosity channel is kept at 1022P. For low-viscosity zones with widths greater than about 100 km and asthenospheric viscosities less than 1018P, we find that viscoelasticity can amplify the perturbations in the moment of inertia by a factor of 4–5 above the elastic contribution within the time span of the wobble period. We have carried out a comparative study on the changes of the inertia tensor from forcings due to surface loading and to faulting. In general the global responses from faulting are found to be much more sensitive to the viscosity structure of the asthenosphere than those produced from surface loading.  相似文献   

13.
Alan Trenhaile   《Geomorphology》2006,76(3-4):316-331
Laboratory experiments have been conducted to measure rock expansion and contraction induced by wetting and drying on intertidal shore platforms in eastern Canada. Almost 800 rock cores and cubes were subjected to between 800 and 1000 tidal cycles. The rocks were put into tidal simulators and inundated in de-ionized water for 1, 6, or 11 h of each 12 h cycle, replicating conditions at the high, mid-, and low tidal levels, respectively. Micro-erosion meter (MEM) stations were also installed in slabs of argillite, sandstone, and basalt. The slabs were subjected to conditions at the high, mid-, and low tidal levels, and a traversing MEM was used to measure rock expansion and contraction as the slabs dried in air with variable temperature and humidity during low tide. The slabs were also exposed for long periods to air with variable humidity. Cores and cubes often gained in weight when they absorbed or adsorbed water and retained it for periods ranging from 1 month to 1 year or more, although this did not appear to be a prelude to rock breakdown. Tidal wetting and drying and exposure to humid air caused the argillite and basalt slabs to expand and contract by up to 0.14 and 0.04 mm, respectively. Wetting and drying did not induce expansion or contraction of the sandstone, but the sandstone did respond slightly to changes in air humidity. The greatest amount of expansion and contraction was at the high tidal level in the argillite and at the mid-tidal level in the basalt, which is consistent with previously measured rates of downwearing.  相似文献   

14.
Summary. An interpretation of the geomagnetic inductive response function, C (ω, 0), observed at Kiruna in northern Sweden, is herein undertaken. The bounds of acceptable solutions are initially discovered by a Monte-Carlo random search procedure, and the best-fitting solutions are examined by the application of linear theory to the problem. The data are shown to have a higher degree of internal consistency than that described by the estimated variances of each datum. A further Monte-Carlo inversion of the variance- reduced data set gives solutions with well defined model parameters.
The two major features of the models are: (1) a small, or non-existent, electrical conductivity variation across the seismic Moho boundary, and (2) the unequivocal existence of an electrical asthenosphere, under the Fenno- scandian shield, beginning at a depth of between 155–185 km, and of 60km minimum thickness. Both of these observations have seismic counterparts.
Finally, possible mantle temperature profiles are deduced which depend on the assumptions and laboratory data employed.  相似文献   

15.
Summary. Data from Japanese local seismograph networks suggest that the stresses in double seismic zones are in-plate compression for the upper zone and in-plate tension for the lower zone; the stresses do not necessarily appear to be down-dip. It may therefore be possible to identify other double seismic zones on the basis of data which indicate that events with differing orientations of in-plate stresses occur in a given segment of slab.
A global survey of published focal mechanisms for intermediate depth earthquakes suggests that the stress in the slab is controlled, at least in part, by the age of the slab and the rate of convergence. Old and slow slabs are under in-plate tensile stresses and the amount of in-plate compression in the slab increases with increasing convergence rate or decreasing slab age. Young and fast slabs are an exception to this trend; all such slabs are down-dip tensile. Since these slabs all subduct under continents, they may be bent by continental loading. Double seismic zones are not a feature common to all subduction zones and are only observed in slabs which are not dominated by tensile or compressive stresses.
Unbending of the lithosphere and upper mantle phase changes are unlikely to be the causes of the major features of double zones, although they may contribute to producing some of their characteristics. Sagging or thermal effects, possibly aided by asthenospheric relative motion, may produce the local deviatoric stresses that cause double zones.  相似文献   

16.
The McKenzie model proposed in 1978, which is widely used in calculating the thermal history of rift basins and other extensional basins, incorrectly assumes that all heat passing through the lithosphere originates below the lithosphere. In reality, heat from radiogenic sources within the lithosphere, especially in the upper crust, may represent more than half the heat flow at the top of basement. Thinning of the lithosphere during extension does indeed result in an increase of heat flowing from the asthenosphere, but this thinning also reduces the radiogenic heat from within the lithosphere. Because these two effects cancel to a large degree, the direct effects of lithospheric extension on heat flow at the top of basement are smaller than those predicted by the McKenzie model. Because of permanent loss of radiogenic material by lithospheric thinning, the heat flow at the top of basement long after rifting will be lower than the pre-rift heat flow.The McKenzie model predicts an instantaneous increase in heat flow during rifting. The Morgan model proposed in 1983, however, predicts a substantial time delay in the arrival of the higher heat flow from the asthenosphere at the top of basement or within sediments. Using the Morgan model, heat flow during the early stages of rifting will actually be lower than prior to rifting, because the time delay in the loss of radiogenic heat is less than the time delay in arrival of new heat from the asthenosphere.  相似文献   

17.
Summary. Motion of the lithosphere over a low viscosity asthenosphere concentrates shear and thus energy dissipation in the asthenosphere. This heat source warms the asthenosphere and, in extreme circumstances, may lead to thermal instabilities. The conditions for thermal stability have been investigated by Melosh who supposed that constant stress acted on the plate, and by Yuen & Schubert who assumed constant velocity boundary conditions. In this paper we investigate a simple analytical model which behaves qualitatively like the more complex systems. This model reproduces the results of Melosh for constant stress and of Yuen & Schubert for constant velocity. The velocity—shear stress characteristic curve for this model shows three branches. The stability of solutions on each branch is a function of the boundary conditions, whether constant stress or constant velocity. The simplicity of the model allows us to investigate stability when neither constant stress nor constant velocity apply and to study the structure of the solutions as these limits are approached. A relation between the velocity of a plate and the driving force is constructed. A loading-line analysis specifies the actual stress and velocity of the plate. Although the solutions are unique for many combinations of the loading-line parameters, there is a region of multiple solutions. These solutions exhibit the characteristics of a 'cusp catastrophe' both a low velocity and a high velocity state are stable, while an intermediate state is unstable. Continental lithosphere may lie in this region, leading to epirogenic movements when the plate changes its velocity with respect to the mantle. Oceanic lithosphere almost certainly moves in the low velocity state.  相似文献   

18.
Summary. Observations of surface waves crossing ocean basins indicate that, at least in some cases, a better interpretation of what was previously considered to be the simultaneous arrival of fundamental Love and second Rayleigh modes, with similar group velocities, is the coupled motion of the second-generalized mode surface wave propagating in an anisotropic upper-mantle. This mode displays the polarization anomalies expected from anisotropic alignments caused by syntectonic recrystallization of olivine in a zone of horizontal shear. The polarity of the pattern of anomalies, found in waves crossing the Pacific Basin, suggests that the lithosphere is dragging the asthenosphere.  相似文献   

19.
We present the first results of a high-resolution teleseismic traveltime tomography and seismic anisotropy study of the lithosphere–asthenosphere system beneath the western Bohemian Massif. The initial high-resolution tomography down to a depth of 250 km did not image any columnar low-velocity anomaly which could be interpreted as a mantle plume anticipated beneath the Eger Rift, similar to recent findings of small plumes beneath the French Massif Central and the Eifel in Germany. Alternatively, we interpret the broad low-velocity anomaly beneath the Eger Rift by an upwelling of the lithosphere–asthenosphere transition. We also map lateral variations of seismic anisotropy of the mantle lithosphere from spatial variations of P -wave delay times and the shear wave splitting. Three major domains characterised by different orientations of seismic anisotropy correspond to the major tectonic units—Saxothuringian, Moldanubian and the Teplá-Barrandian—and their fabrics fit to those found in our previous studies of mantle anisotropy on large European scales.  相似文献   

20.
The deformation at the core–mantle boundary produced by the 2004 Sumatra earthquake is investigated by means of a semi-analytic theoretical model of global coseismic and postseismic deformation, predicting a millimetric coseismic perturbation over a large portion of the core–mantle boundary. Spectral features of such deformations are analysed and discussed. The time-dependent postseismic evolution of the elliptical part of the gravity field ( J 2) is also computed for different asthenosphere viscosity models. Our results show that, for asthenospheric viscosities smaller than 1018 Pa s, the postseismic J 2 variation in the next years is expected to leave a detectable signal in geodetic observations.  相似文献   

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

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