首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 546 毫秒
1.
The deformation of the lithosphere due to temperature anomalies caused by a heat source located below or within the layer is usually modelled as one of Pratt local isostasy. A more appropriate model is one of rheologically layered lithosphere comprising of a stiff viscoelastic or elastic layer overlying a weaker viscoelastic layer. The surface deformations are a result of not only the perturbations in body forces due to density changes, but thermal bending moments. In geophysically realistic situations the former contribution dominates. Pratt isostasy is attained if the stresses in the entire lithosphere are allowed to relax and this end state is not contingent upon the lithosphere being confined against horizontal deformation. In a rheologically layered lithosphere, even though the non-isostatic thermal stresses persist in the upper layer, the surface deformations are indistinguishable from that of local isostasy if the horizontal dimension of the heat source exceeds about three times the effective elastic thickness of the lithosphere.  相似文献   

2.
In this paper, the lithosphere is considered to be a homogeneous elastic spherical shell for the sake of simplicity and the stress equations for the base of the lithosphere are taken as boundary conditions. Then the stress equations are obtained for use in the computation of the stress field in the lithosphere with the satellite gravitational harmonic coefficients. The 5 × 5° global stress field in the lithosphere is computed from harmonics of 2–30°. The directions of principal stresses of this stress field agree favourably with the directions of principal stresses indicated by mid-plate earthquake mechanisms, in situ stress measurements and sensitive geological features. This result indicates that the drag forces exerted on the base of the lithosphere, due to gravitational mantle convection, may be among the driving forces of plate motion and a major source for the stress field in the lithosphere.  相似文献   

3.
A linear analysis of the stability of the lithosphere considered as a viscoelastic layer with an equilibrium vertical gradient of temperature is carried out. The problem is solved with a complete system of linearized equations of a continuous medium represented in the dimensionless form and containing a set of dimensionless parameters that determine thermomechanical properties of the lithosphere. As a result of the stability analysis, decrements are found that give the time dependence of perturbations and correspond to high-frequency seismic waves and low-frequency tectonic waves. The frequency and velocity of seismic waves are determined by the elasticity and inertial properties of the lithosphere, and their attenuation, by viscous properties of the lithosphere. The temperature gradient existing in the lithosphere influences seismic waves very weakly. On the contrary, the pattern of tectonic waves is controlled by the temperature gradient and viscous properties, while the effect of elastic and inertial properties on these waves is negligibly small. The stability of a viscoelastic lithosphere is examined using such rheological models as the Maxwell, standard linear, and Andrade media (the frequency of tectonic waves is zero in the Maxwell medium).  相似文献   

4.
In this paper the mid-ocean ridge axial valley is modelled as a steady-state lithospheric neck in which lithospheric stretching balances lithospheric accretion. Conversely, the axial high is a steady-state lithospheric bulge. The lithosphere is modelled as a thin plate with a Newtonian rheology. It is shown that an axial valley will occur if the rate of viscosity increase away from the ridge axis is faster than the rate at which accretion decreases. An axial high will occur if the opposite condition holds. This is consistent with the observation that axial valleys occur at low spreading rates and axial highs at high spreading rates. By fitting our model to profiles across the Mid-Atlantic Ridge and the East Pacific Rise and assuming the lithospheric thickness at the ridge axis to be 5 km, we find accretion widths of 6–8 km. We find the width over which there is a significant increase in lithospheric viscosity to be also 6–8 km.  相似文献   

5.
In this paper, the effects of a mass damper on the rocking motion of a non‐symmetric rigid block‐like structure, subject to different seismic excitation, are investigated. The damper is modelled as a single degree of freedom oscillating mass, running at the top of the block and connected to it by a linear visco‐elastic device. The equations of rocking motion, the uplift and the impact conditions are derived. A nondimensionalisation of the governing equations is performed with the aim to obtain an extensive parametric analysis. The results are achieved by numerical integration of these equations. The slenderness and the base of the rigid block, and the eccentricity of the centre of mass are taken as variable parameters in the analyses. The main objective of the study is to check the performance of the damper versus the spectral characteristics of the seismic input. Three earthquake registrations with different frequency contents are used in the analyses. The results show that the presence of the mass damper leads to different levels of improvement of the response of the system, depending on the spectral characteristics of the seismic input. Curves providing the overturning slenderness of blocks of specific sizes versus the characteristics of the TMD are obtained. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

6.
The formation of the thermal cross section of the lithosphere and mantle upon the interaction between the mantle convection and the immobile continent surrounded by the oceanic lithosphere is studied by numerical modeling. The convective temperature and velocity fields and then the averaged geotherms for subcontinental and suboceanic regions up to the boundary with the core are calculated from the solution of convection equations with a jump in viscosity in the continental zone. Using the experimental data on the solidus temperature in the rocks of the upper mantle, the average thickness of the continental and oceanic lithosphere is estimated at 190 and 30 km, respectively. The effect of a hot spot formed in the subcontinental upper mantle at a depth of 250–500 km, which has not been previously noted, is revealed. Although the temperature in this zone is typically assumed to be close to adiabatic, the calculations show that it is actually higher than adiabatic by up to 200°C. The physical mechanism responsible for this effect is associated with the accumulation of convective heat beneath the thermally insulating layer of the continental lithosphere. The revealed anomalies can be important in studying the phase and mineral transformations at the base of the lithosphere and in the regional geodynamical reconstructions.  相似文献   

7.
We modelled the thickness and seismic anisotropy of the subcrustal lithosphere from the variations of P-wave delay times and the shear-wave splitting observed at seismological observatories and portable stations in the western part of the Bohemian Massif. The Saxothuringian lithosphere is characterized by a total thickness between 90 and 120 km, the Moldanubian lithosphere is generally thicker –120-140 km, on the average. The subcrustal lithosphere of both units is characterised by divergently dipping anisotropic structures and the suture between them is marked by a lithosphere thinning to about 80km. Within the subcrustal lithosphere a complex structure of the transition of both units extends to about 150 km toward the south. We suggest that the Saxothuringian-Moldanubian suture has created a zone of mechanical predisposition for the Tertiary Ohe (Eger) Graben, as well as for the occurrence of earthquake swarms in the region. Most earthquakes occur within the brittle part of the upper crust above the crossing of the suture between the Saxothuringian in the north and the Moldanubian and the Tepl´-Barrandian in the south, with the tectonically active Mariánské Lázn fault.  相似文献   

8.
In order to study the relationship between mantle flow and global tectogenesis, we present a 3-D spherical shell model with incompressible Newtonian fluid medium to simulate mantle flow which fits the global tectogenesis quite well. The governing equations are derived in spherical coordinates. Both the thermal buoyancy force and the self-gravitation are taken into account. The velocity and pressure coupled with temperature are computed, using the finite-element method with a punitive factor. The results show that the lithosphere, as the boundary layer of the earth's thermodynamic system, moves with the entire mantle. Both its horizontal and vertical movements are the results of the earth's thermal motion. The orogenesis occurs not only in the collision zones at the plates' boundaries, but also occurs within the plates. If the core-mantle boundary is impermeable and the viscosity of the lower mantle is considerable, the vertical movement is mostly confined to the upper mantle. The directions of the asthenospheric movements are not fully consistent with those of the lithospheric movements. The depths of spreading movements beneath all ridges are less than 220 km. In some regions, the shear stresses, acting on the base of the lithosphere by the asthenosphere, are the main driving force; but in other regions, the shear stresses are the resisting force.  相似文献   

9.
We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.  相似文献   

10.
11.
Vertical dynamic response of a disk on a saturated poroelastic half-space   总被引:2,自引:0,他引:2  
This paper considers the vertical dynamic response of a disk on a saturated poroelastic half-space. Firstly the pressure-solid displacement form of the harmonic equations of motion for a poroelastic solid are developed from the form of the equations originally presented by Biot. These equations are solved by a new method. Then the mixed boundary value problem for the vertical harmonic vibration of a disk on a poroelastic half-space is studied. The two types of drainage conditions at the surface of the poroelastic half-space are considered: (a) the surface of the poroelastic half-space is assumed to be completely pervious both within and exterior to the plate; (b) The interface between the plate and the poroelastic half-space is assumed to be impervious and the exterior region is assumed to be pervious. By using the Hankel transform techniques, the paper develops the governing dual integral equations. These governing integral equations are further reduced to systems of standard Fredholm integral equations of the second kind by Abel transform.  相似文献   

12.
The scaled boundary finite‐element method is extended to simulate time‐harmonic responses of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. The unbounded domains and the elasticity matrices are transformed to the scaled boundary coordinates. The scaled boundary finite‐element equation in displacement amplitudes are derived directly from the governing equations of elastodynamics. To enforce the radiation condition at infinity, an asymptotic expansion of the dynamic‐stiffness matrix for high frequency is developed. The dynamic‐stiffness matrix at lower frequency is obtained by numerical integration of ordinary differential equations. Only the boundary is discretized yielding a reduction of the spatial dimension by one. No fundamental solution is required. Material anisotropy is modelled without additional efforts. Examples of two‐ and three‐dimensional non‐homogeneous isotropic and transversely isotropic unbounded domains are presented. The results demonstrate the accuracy and simplicity of the scaled boundary finite‐element method. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

13.
John Z. Shi  Li‐Feng Lu 《水文研究》2007,21(13):1780-1786
A model of the wave and current boundary‐layer structure was developed using the k–ε turbulent closure model. The finite‐difference method was used to solve the governing equations. Vertical logarithmic grids and equal time steps were adopted. The following modelled simulations were obtained: (1) vertical profiles of wave velocity amplitude, eddy viscosity coefficient and turbulent kinetic energy with waves only; (2) vertical profiles of wave velocity amplitude, mean current velocity, eddy viscosity coefficient and turbulent kinetic energy with waves having a following current. To test the validity and the rationality of the present model, vertical profiles of modelled wave velocity amplitude and mean velocity were compared with corresponding experimental results available in the literature. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

14.
Locations of the Eger Rift, Cheb Basin, Quaternary volcanoes, crustal earthquake swarms and exhalation centers of CO2 and 3He of mantle origin correlate with the tectonic fabric of the mantle lithosphere modelled from seismic anisotropy. We suggest that positions of the seismic and volcanic phenomena, as well as of the Cenozoic sedimentary basins, correlate with a “triple junction” of three mantle lithospheres distinguished by different orientations of their tectonic fabric consistent within each unit. The three mantle domains most probably belong to the originally separated microcontinents – the Saxothuringian, Teplá-Barrandian and Moldanubian – assembled during the Variscan orogeny. Cenozoic extension reactivated the junction and locally thinned the crust and mantle lithosphere. The rigid part of the crust, characterized by the presence of earthquake foci, decoupled near the junction from the mantle probably during the Variscan. The boundaries (transitions) of three mantle domains provided open pathways for Quaternary volcanism and the ascent of 3He- and CO2-rich fluids released from the asthenosphere. The deepest earthquakes, interpreted as an upper limit of the brittle–ductile transition in the crust, are shallower above the junction of the mantle blocks (at about 12 km) than above the more stable Saxothuringian mantle lithosphere (at about 20 km), probably due to a higher heat flow and presence of fluids.  相似文献   

15.
《Journal of Geodynamics》2008,45(3-5):149-159
Locations of the Eger Rift, Cheb Basin, Quaternary volcanoes, crustal earthquake swarms and exhalation centers of CO2 and 3He of mantle origin correlate with the tectonic fabric of the mantle lithosphere modelled from seismic anisotropy. We suggest that positions of the seismic and volcanic phenomena, as well as of the Cenozoic sedimentary basins, correlate with a “triple junction” of three mantle lithospheres distinguished by different orientations of their tectonic fabric consistent within each unit. The three mantle domains most probably belong to the originally separated microcontinents – the Saxothuringian, Teplá-Barrandian and Moldanubian – assembled during the Variscan orogeny. Cenozoic extension reactivated the junction and locally thinned the crust and mantle lithosphere. The rigid part of the crust, characterized by the presence of earthquake foci, decoupled near the junction from the mantle probably during the Variscan. The boundaries (transitions) of three mantle domains provided open pathways for Quaternary volcanism and the ascent of 3He- and CO2-rich fluids released from the asthenosphere. The deepest earthquakes, interpreted as an upper limit of the brittle–ductile transition in the crust, are shallower above the junction of the mantle blocks (at about 12 km) than above the more stable Saxothuringian mantle lithosphere (at about 20 km), probably due to a higher heat flow and presence of fluids.  相似文献   

16.
This study considers two-dimensional mantle flow beneath a rigid lithosphere. The lithosphere which forms the upper boundary of a convecting region moves with a prescribed uniform horizontal velocity, and thickens with distance from the accreting plate boundary as it cools. Beneath the lithosphere, the mantle deforms viscously by diffusion creep and is heated radiogenically from within. Solutions for thermal convection beneath the lithosphere are obtained by finite-difference methods. Two important conclusions have resulted from this study: (1) convective patterns of large aspect ratio are stable beneath a rigid moving lithosphere; (2) even for a lithosphere velocity as small as 3 cm/yr. and a Rayleigh number as large as 106, mantle circulation with large aspect ratio is driven dominantly by the motion of the lithosphere rather than by temperature gradients within the flow. Gravity, topography and heat flow are determined and implications for convection in the upper mantle are discussed.  相似文献   

17.
Using thin elastic plate theory and neglecting horizontal applied forces, a universal deflection profile applicable to many oceanic trenches is derived. This theoretical profile is compared with bathymetric profiles from the central Aleutian, Kuril, Bonin, and Mariana trench-outer rise regions. The profiles were corrected for sediment thickness and age variation of the lithosphere. Good agreement between theory and observation is found. The distance from the first point of zero deflection seaward of the trench to the point of maximum height of the outer rise is directly related to the flexural rigidity of the lithosphere. The thickness of the elastic lithosphere is found to vary between 20 and 29 km for the trench profiles considered. The good agreement obtained shows that horizontal forces may be neglected and that the bending lithosphere behaves elastically in the cases considered. The analysis shows that only unreasonably large horizontal forces would affect the universal deflection curve. It is concluded that although the near-surface lithosphere may be subject to brittle fracture, the deeper lithosphere is capable of transmitting elastic stresses as high as 9 kbar.  相似文献   

18.
The modern concepts of the rheology of viscous mantle and brittle lithosphere, as well as the results of the numerical experiments on the processes in a heated layer with a viscosity dependent on pressure, temperature, and shear stress, are reviewed. These dependences are inferred from the laboratory studies of olivine and measurements of postglacial rebound (glacial isostatic adjustment) and geoid anomalies. The numerical solution of classical conservation equations for mass, heat, and momentum shows that thermal convection with a highly viscous rigid lithosphere develops in the layer with the parameters of the mantle with the considered rheology under a temperature difference of 3500 K, without any special additional conditions due to the self-organization of the material. If the viscosity parameters of the lithosphere correspond to dry olivine, the lithosphere remains monolithic (unbroken). At a lower strength (probably due to the effects of water), the lithosphere splits into a set of separate rigid plates divided by the ridges and subduction zones. The plates submerge into the mantle, and their material is involved in the convective circulation. The results of the numerical experiment may serve as direct empirical evidence to validate the basic concepts of the theory of plate tectonics; these experiments also reveal some new features of the mantle convection. The probable structure of the flows in the upper and lower mantle (including the asthenosphere), which shows the primary role of the lithospheric plates, is demonstrated for the first time.  相似文献   

19.
An analytical approach is developed to study the dynamic response of a flexible plate on single-layered saturated soil. The analysis is based on Biot's two-phased theory of poroelasticity and also on the classical thin-plate theory. First, the governing differential equations for saturated soil are solved by the use of Hankel transform. The general solutions of the skeleton displacements, stresses, and pore pressures, derived in the transformed domain, are subsequently incorporated into the imposed boundary conditions, which leads to a set of dual integral equations describing the corresponding mixed boundary value problem. These governing integral equations are finally reduced to the Fredholm integral equations of the second kind and solved by standard numerical procedures. The accuracy of the present solution is validated via comparisons with existing solutions for an ideal elastic half-space. Furthermore, some numerical results are presented to show the influences of the layer depth, the plate flexibility, and the soil porosity on the dynamic compliances.  相似文献   

20.
Research on high-resolution tomography, three-dimensional anisotropy modelled from shear-wave splitting and P-residual analysis, as well as receiver function interpretation, are in progress with the aim to image the crust and uppermost mantle velocities in the Bohemian Massif. Structure of the deep lithosphere and location of boundaries of mantle domains will be compared with tectonics of the crust and limits of terranes derived by geologic and palaeomagnetic methods.  相似文献   

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

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