共查询到20条相似文献，搜索用时 46 毫秒
1.
Physical models of subduction investigate the impact of regional mantle flow on the structure of the subducted slab and deformation of the downgoing and overriding plates. The initial mantle flow direction beneath the overriding plate can be horizontal or vertical, depending on its location with respect to the asthenospheric flow field. Imposed mantle flow produces either over or underpressure on the lower surface of the slab depending on the initial mantle flow pattern (horizontal or vertical, respectively). Overpressure promotes shallow dip subduction while underpressure tends to steepen the slab. Horizontal mantle flow with rates of 1–10 cm yr−1 provides sufficient overpressure on a dense subducting lithosphere to obtain a subduction angle of ∼60° , while the same lithospheric slab sinks vertically when no flow is imposed. Vertical drag force (due to downward mantle flow) exerted on a slab can result in steep subduction if the slab is neutrally buoyant but fails to produce steep subduction of buoyant oceanic lithosphere. The strain regime in the overriding plate due to the asthenospheric drag force depends largely on slab geometry. When the slab dip is steeper than the interplate zone, the drag force produces negative additional normal stress on the interplate zone and tensile horizontal stress in the overriding plate. When the slab dip is shallower than the interplate zone, an additional positive normal stress is produced on the interplate zone and the overriding plate experiences additional horizontal compressive stress. However, the impact of the mantle drag force on interplate pressure is small compared to the influence of the slab pull force since these stress variations can only be observed when the slab is dense and interplate pressure is low. 相似文献
2.
D. J. Woodward 《Geophysical Journal International》1977,50(2):459-472
Summary. An empirical equation of state, assuming mineralogical equilibrium, is developed for the top 700 km of the mantle. Assuming a uniform viscosity, this equation of state is used to show that the stresses due to the changes in phase induced in a descending lithospheric plate in a subduction zone are an order of magnitude larger than those due to the negative buoyancy of the slab in the asthenosphere. The stresses predicted are well within the power law creep region for likely mantle materials and so the effective viscosity will vary within the slab. Consequently the stresses will be smaller than those of 7.0 × 108 N/m2 obtained here using uniform viscosities. These stresses are relatively compressional near the sides of the slab and tensional in the centre. 相似文献
3.
Summary. Bulletins of the International Seismological Centre (ISC) show very large residuals, up to 15 s early, for arrivals from events in the Tonga–Kermadec subduction zone to the New Zealand network of seismometers. The very early arrivals are confined to events south of about 22°S, and shallower than about 350 km. The waveforms show two distinct phases: an early, emergent, first phase with energy in the high-frequency band 2–10 Hz, and a distinct second phase, containing lower frequency energy, arriving at about the time predicted by JB tables.
The residuals are attributed to propagation through the cold, subducted lithosphere, which has a seismic velocity 5 per cent faster, on average, than normal. Ray tracing shows that the ray paths lie very close to the slab for events south of 22°S, but pass well beneath the slab for events further north, corresponding to the change in residual pattern. This characteristic of the ray paths is due to the curved shape of the seismic zone, and in particular to the bend in the zone where the Louisville ridge intersects the trench at 25°S.
The residuals can only be explained if the high velocity anomaly extends to a depth of 450 km in the region of the gap in deep seismicity from 32 to 36°S. The very high-frequency character of the first phase requires the path from the bottom of the slab to the stations to be of high Q , and to transmit 2–10 Hz energy with little attenuation.
The absence of low-frequency energy in the first phase is due to the narrowness of the high-velocity slab, which transmits only short-wavelength waves. The second phase, which contains low frequencies, is identified as a P -wave travelling beneath the subducted slab in normal mantle. There is no need to invoke any special structures, such as low-velocity waveguides or reflectors, to explain any of the observations. The S -wave arrivals show similar effects. 相似文献
The residuals are attributed to propagation through the cold, subducted lithosphere, which has a seismic velocity 5 per cent faster, on average, than normal. Ray tracing shows that the ray paths lie very close to the slab for events south of 22°S, but pass well beneath the slab for events further north, corresponding to the change in residual pattern. This characteristic of the ray paths is due to the curved shape of the seismic zone, and in particular to the bend in the zone where the Louisville ridge intersects the trench at 25°S.
The residuals can only be explained if the high velocity anomaly extends to a depth of 450 km in the region of the gap in deep seismicity from 32 to 36°S. The very high-frequency character of the first phase requires the path from the bottom of the slab to the stations to be of high Q , and to transmit 2–10 Hz energy with little attenuation.
The absence of low-frequency energy in the first phase is due to the narrowness of the high-velocity slab, which transmits only short-wavelength waves. The second phase, which contains low frequencies, is identified as a P -wave travelling beneath the subducted slab in normal mantle. There is no need to invoke any special structures, such as low-velocity waveguides or reflectors, to explain any of the observations. The S -wave arrivals show similar effects. 相似文献
4.
Megan Anderson Patricia Alvarado † George Zandt Susan Beck 《Geophysical Journal International》2007,171(1):419-434
We use data from the Chile Argentina Geophysical Experiment (CHARGE) broad-band seismic deployment to refine past observations of the geometry and deformation within the subducting slab in the South American subduction zone between 30°S and 36°S. This region contains a zone of flat slab subduction where the subducting Nazca Plate flattens at a depth of ∼100 km and extends ∼300 km eastward before continuing its descent into the mantle. We use a grid-search multiple-event earthquake relocation technique to relocate 1098 events within the subducting slab and generate contours of the Wadati-Benioff zone. These contours reflect slab geometries from previous studies of intermediate-depth seismicity in this region with some small but important deviations. Our hypocentres indicate that the shallowest portion of the flat slab is associated with the inferred location of the subducting Juan Fernández Ridge at 31°S and that the slab deepens both to the south and the north of this region. We have also determined first motion focal mechanisms for ∼180 of the slab earthquakes. The subhorizontal T -axis solutions for these events are almost entirely consistent with a slab pull interpretation, especially when compared to our newly inferred slab geometry. Deviations of T -axes from the direction of slab dip may be explained with a gap within the subducting slab below 150 km in the vicinity of the transition from flat to normal subducting geometry around 33°S. 相似文献
5.
Experiments about working mechanism and mechanical characteristics of the DX model pile foundation under lateral dynamic and static loading were conducted by using a model system of the dynamic frozen ... 相似文献
6.
7.
Simple analytic model for subduction zone thermal structure 总被引：4，自引：0，他引：4
J. Huw Davies 《Geophysical Journal International》1999,139(3):823-828
A new analytic model is presented for the thermal structure of subduction zones. It applies to the deeper regions of a subduction zone, where the overriding mantle is no longer rigid but flows parallel to the slab surface. The model captures the development of one thermal boundary layer out into the mantle wedge, and another into the subducting slab. By combining this model with the analytic model of Royden (1993a , b ), which applies to regions in which the overriding plate is rigid, a nearly complete analytic model for the thermal structure of a steady-state subduction zone can be achieved. A good agreement is demonstrated between the output of the combined analytic model and a numerical finite element calculation. The advantages of this analytic approach include (1) efficiency (only limited computing resources are needed); (2) flexibility (non-linear slab shape, and processes such as erosion, and shear heating are easily incorporated); and (3) transparency (the effect of changes in input variables can be seen directly). 相似文献
8.
基于欧拉双流体模型,运用商用CFD软件FLUENT,对挟雨风条件下高速列车的气动特性和运行稳定性进行了三维模拟研究。分析了不同横摆角和暴雨强度下气动力和倾覆力矩的变化特征,并与无降雨条件下的相应结果进行比较。结果表明：受降雨影响,列车受到的阻力、侧向力和倾覆力矩均增大,列车气动性能损耗增加,运行稳定性降低。在模拟的降雨强度和横摆角条件下,阻力系数将提高1.8%~38.6%,侧向力系数提高0.1%~9.16%,导致列车的能耗增加;倾覆力矩系数提高1.8%~10.0%,意味着与净侧风相比,挟雨风进一步降低了列车的运行稳定性。 相似文献
9.
Toshiro Tanimoto 《Geophysical Journal International》1998,134(1):199-206
The state of stress within a bending spherical shell has some special features that are caused by sphericity. While most lithospheres are more like spherical shells than flat plates, our ideas of the state of stress have been dominated by flat-plate models. As a consequence, we might be missing some important aspects of the state of stress within subducting lithospheres. In order to examine this problem, we analyse spherical-shell bending problems from basic equations. We present two approaches to solve spherical-shell bending problems: one by the variational approach, which is suitable for global-scale problems, and the other by the asymptotic equation, which is valid to first order in h/R , where h is the thickness of the lithosphere and R is its curvature radius (i.e. under the assumption of small curvature). The form of the equation for displacement shows that wavelengths of deformation are determined by the spherical (elastic) effect and the gravitational buoyancy effect, for which only the latter effect is included in the usual flat-plate formulations. In the case of the Earth, the buoyancy force is dominant and, consequently, spherical effects are suppressed to a large extent; this explains why flat-plate models have been successful for Earth's lithospheric problems. On the other hand, the state of stress shows interesting spherical effects: while bending (fibre) stress along the subduction zone is always important, bending stress along the trench-strike direction can also be important, in particular when the subduction zone arc is small. Numerical results also indicate that compressive normal stress along the trench-strike direction is important when a subduction zone arc is large. These two stresses, the bending stress and the compressive normal stress, both along the trench-strike direction, may have important implications for intraplate earthquakes at subduction zones. 相似文献
10.
Mikael Beuthe 《Geophysical Journal International》2008,172(2):817-841
Planetary topography can either be modelled as a load supported by the lithosphere, or as a dynamic effect due to lithospheric flexure caused by mantle convection. In both cases the response of the lithosphere to external forces can be calculated with the theory of thin elastic plates or shells. On one-plate planets the spherical geometry of the lithospheric shell plays an important role in the flexure mechanism. So far the equations governing the deformations and stresses of a spherical shell have only been derived under the assumption of a shell of constant thickness. However, local studies of gravity and topography data suggest large variations in the thickness of the lithosphere. In this paper, we obtain the scalar flexure equations governing the deformations of a thin spherical shell with variable thickness or variable Young's modulus. The resulting equations can be solved in succession, except for a system of two simultaneous equations, the solutions of which are the transverse deflection and an associated stress function. In order to include bottom loading generated by mantle convection, we extend the method of stress functions to include loads with a toroidal tangential component. We further show that toroidal tangential displacement always occurs if the shell thickness varies, even in the absence of toroidal loads. We finally prove that the degree-one harmonic components of the transverse deflection and of the toroidal tangential displacement are independent of the elastic properties of the shell and are associated with translational and rotational freedom. While being constrained by the static assumption, degree-one loads can deform the shell and generate stresses. The flexure equations for a shell of variable thickness are useful not only for the prediction of the gravity signal in local admittance studies, but also for the construction of stress maps in tectonic analysis. 相似文献
11.
Steven N. Ward 《Geophysical Journal International》1985,81(2):409-428
Summary. This paper explores the middle ground between complex thermally-coupled viscous flow models and simple corner flow models of island arc environments. The calculation retains the density-driven nature of convection and relaxes the geometrical constraints of corner flow, yet still provides semianalytical solutions for velocity and stress. A novel aspect of the procedure is its allowance for a coupled elastic lithosphere on top of a Newtonian viscous mantle. Initially, simple box-like density drivers illustrate how vertical and horizontal forces are transmitted through the mantle and how the lithosphere responds by trench formation. The flexural strength of the lithosphere spatially broadens the surface topography and gravity anomalies relative to the functional form of the vertical flow stresses applied to the plate base. I find that drivers in the form of inclined subducting slabs cannot induce self-driven parallel flow; however, the necessary flow can be provided by supplying a basal drag of 1–5 MPa to the mantle from the oceanic lithosphere. These basal drag forces create regional lithospheric stress and they should be quantifiable through seismic observations of the neutral surface. The existence of a shallow elevated phase transition is suggested in two slab models of 300 km length where a maximum excess density of 0.2 g cm−3 was needed to generate an acceptable mantle flow. A North New Hebrides subduction model which satisfies flow requirements and reproduces general features of topography and gravity contains a high shear stress zone (75 MPa) around the upper slab surface to a depth of 150 km and a deviatoric tensional stress in the back arc to a depth of 70 km. The lithospheric stress state of this model suggests that slab detachment is possible through whole plate fracture. 相似文献
12.
Mathieu Dumberry 《Geophysical Journal International》2008,172(3):903-920
A decadal polar motion with an amplitude of approximately 25 milliarcsecs (mas) is observed over the last century, a motion known as the Markowitz wobble. The origin of this motion remains unknown. In this paper, we investigate the possibility that a time-dependent axial misalignment between the density structures of the inner core and mantle can explain this signal. The longitudinal displacement of the inner core density structure leads to a change in the global moment of inertia of the Earth. In addition, as a result of the density misalignment, a gravitational equatorial torque leads to a tilt of the oblate geometric figure of the inner core, causing a further change in the global moment of inertia. To conserve angular momentum, an adjustment of the rotation vector must occur, leading to a polar motion. We develop theoretical expressions for the change in the moment of inertia and the gravitational torque in terms of the angle of longitudinal misalignment and the density structure of the mantle. A model to compute the polar motion in response to time-dependent axial inner core rotations is also presented. We show that the polar motion produced by this mechanism can be polarized about a longitudinal axis and is expected to have decadal periodicities, two general characteristics of the Markowitz wobble. The amplitude of the polar motion depends primarily on the Y 1 2 spherical harmonic component of mantle density, on the longitudinal misalignment between the inner core and mantle, and on the bulk viscosity of the inner core. We establish constraints on the first two of these quantities from considerations of the axial component of this gravitational torque and from observed changes in length of day. These constraints suggest that the maximum polar motion from this mechanism is smaller than 1 mas, and too small to explain the Markowitz wobble. 相似文献
13.
State of stress in the Southern Tyrrhenian subduction zone from fault-plane solutions 总被引：3，自引：0，他引：3
In this paper we present revised locations and original focal mechanisms computed for intermediate and deep earthquakes that occurred within the Southern Tyrrhenian subduction zone between 1988 and 1994, in order to improve our knowledge of the state of stress for this compressional margin. In particular, we define the stress distribution within a large portion of the descending slab, between 40 and about 450 km depth. The seismicity distribution reveals a continuous 40–50 km thick slab that abruptly increases its dip from subhorizontal in the Ionian Sea to a constant 70° dip in the Tyrrhenian. We computed focal mechanisms for events with magnitudes ranging from 2.7 and 5.7, obtaining the distribution of P - and T -axes for many events for which centroid moment tensor (CMT) solutions are not available, thus enabling the sampling of a larger depth range compared to previous studies. We define three portions of the slab characterized by different distributions of P - and T -axes. A general down-dip compression is found between 165 and 370 km depth, whereas in the upper part of the slab (40–165 km depth) the fault-plane solutions are strongly heterogeneous. Below 370 km the P -axes of the few deep events located further to the north have a shallower dip and are not aligned with the 70° dipping slab, possibly suggesting that they belong to a separated piece of subducted lithosphere. There is a good correspondence between the depth range in which the P -axes plunge closer to the slab dip (∼ 70°) and the interval characterized by the highest seismic energy release (190–370 km). 相似文献
14.
Norman H. Sleep 《Geophysical Journal International》1975,42(3):827-857
Summary. The flow pattern, stress distribution, topography, and gravity anomalies were computed from numerical models having density and viscosity distributions resemblant to the Aleutian arc. The results were compatible with the hypothesis that the excess density of the slab drives its descent and that hydrodynamic forces are responsible for topographic and gravity highs over the outer rise seaward of the trench and the frontal arc and lows over the trench. In models with simple distributions of rheological parameters, the force from the slab was transmitted directly upward producing a negative gravity anomaly over the arc. Material with low resistance to flow was needed along the fault plane above the slab or within the crust of the frontal arc and within the wedge of asthenosphere above the slab to reduce that force and to allow the horizontal lithosphere to move with the slab. Models with the resistance to flow thus reduced had outer rises, deep trenches, horizontal tension seaward of the trench, horizontal compression under the trench, and downdip tension in the slab. Free air gravity anomalies, which are the sum of between deflections of the free surface due to hydrodynamic forces and direct attractions from the masses driving the flow, were not fit excellently by any of the models, in part because the coarse grid used precluded accurate representation of the fault zone above the slab and the frontal arc. An alternate to the hypothesis that about 5 kb of stress on the fault plane is needed to produce an outer rise is offered by these models. Shear stress between the slab and the island arc was always below 700 bars in the more successful models if the density distribution was scaled to match the topography of the trench. This is much less than the 2000 bars stresses needed if frictional heating causes island arc volcanism. 相似文献
15.
We investigate the effect of laterally varying earth structure on centroid moment tensor inversions using fundamental mode mantle waves. Theoretical seismograms are calculated using a full formulation of surface wave ray theory. Calculations are made using a variety of global tomographic earth models. Results are compared with those obtained using the so-called great-circle approximation, which assumes that phase corrections are given in terms of mean phase slowness along the great circle, and which neglects amplitude effects of heterogeneity. Synthetic tests suggest that even source parameters which fit the data very well may have large errors due to incomplete knowledge of lateral heterogeneity. The method is applied to 31 shallow, large earthquakes. For a given earthquake, the focal mechanisms calculated using different earth models and different forward modelling techniques can significantly vary. We provide a range of selected solutions based on the fit to the data, rather than one single solution. Difficulties in constraining the dip-slip components of the seismic moment tensor often produce overestimates of seismic moment, leading to near vertical dip-slip mechanisms. This happens more commonly for earth models not fitting the data well, confirming that more accurate modelling of lateral heterogeneity can help to constrain the dip-slip components of the seismic moment tensor. 相似文献
16.
Summary. The Wiener—Hopf technique is used to obtain an exact analytical solution for the problem of H -polarization induction over the edge of a perfectly conducting thin sheet, representing an ocean, electrically connected to a perfectly conducting mantle through a slab of finite conductivity and thickness, which represents the Earths crust. It is shown that the induced currents resulting from this type of induction process are drawn up into the sea from the cust and mantle with the greatest concentration of current near the ocean edge. The surface impedance over the land surface, is calculated for various mantle depths and is shown to increase sharply as the coastline is approached. The magnetic field along the ocean floor is also plotted as a function of distance from the coastline, and the results are found to agree very well with those calculated previously by approximate and numerical methods. 相似文献
17.
Multilayer, spherically stratified, self-gravitating relaxation models with a large number of layers (more than 100) can be dealt with analytically. Relaxation processes are studied for both Heaviside surface loads and tidal forcings. Simulations of the relaxation process of a realistic earth model with an incompressible Maxwell rheology show that models containing about 30 to 40 layers have reached continuum limits on all timescales and for all harmonic degrees up to at least 150 whenever an elastic lithosphere is present, irrespective of the viscosity profile in the mantle. In particular, fine-graded stratification of the shallow layers proves to be important for high harmonic degrees in these models. The models produce correct long-time (fluid) limits. It is shown that differences in the transient behaviour of the various models are due to the applied volume-averaging procedure of the rheological parameters. Our earlier proposed hypothesis that purported shortcomings in the fundamental physics of (discrete) normal-mode theory are artificial consequences of numerical inaccuracies, theoretical misinterpretations and the use of incomplete sets of normal modes is reinforced by the results presented. We show explicitly that the models produce both continuous behaviour resulting from continuous rheological stratifications and discrete behaviour resulting from sharp density contrasts, as at the outer surface and the core-mantle boundary. The differences between volume-averaged models and fixed-boundary contrast models are outlined. Reducing many-layer models with a volume-averaging procedure before employing a normal mode analysis is both economical and highly accurate on all timescales and for all spherical harmonic degrees. The procedure minimizes the chances of missing contributing modes, while using models with more layers will not result in any substantial increase of accuracy. 相似文献
18.
Summary. An inversion of ISC travel-time data from selected earthquakes in the distance range 30°-90° to 53 stations in Central Europe has been used to model velocity down to 600 km depth. The model explains 0.1–0.2s of the residuals, as for other array studies, leaving 0.5 s unexplained as noise. The uppermost 100 km of the mantle and crust contains inhomogeneities that correlate remarkably well with the geology. This may be due to deep-seated thermal anomalies or, in some areas, to delays introduced by passage of the rays through sedimentary cover. The deeper anomalies are smaller and unrelated to those in the lithosphere, which suggests that the asthenosphere is decoupled from the rigid lithosphere. The structure at 600 km depth is again quite inhomogeneous and might be due to undulations of the 650 km discontinuity. The models show some suggestion of a high velocity slab trending from east to west beneath the Alps. 相似文献
19.
Summary A laboratory and theoretical study of the stability of conduits of buoyant fluid in a viscous shear flow has been conducted. The object of the study is to explain the formation of discrete islands in island chains such as the Hawaiian Emperor seamount chain, and to investigate a new method by which the variation of shear with depth in the mantle may be determined. The conduits were made by injecting oil into a more viscous oil of greater density. Initially a growing chamber of lower viscosity oil formed near the injector, but when the chamber got sufficiently large it rose as a buoyant spheroid. Behind this trailed a vertical cylindrical conduit through which fluid could continue to rise to the surface as long as the source continued. If the more viscous fluid was sheared laterally the conduit was gradually rotated to a more horizontal position. The diameter of the conduit increased with time due to a decreasing component of gravitational force along the axis of the conduit. When the conduit was tilted to more than 60° with the vertical, it began to go unstable by developing bumps which ultimately initiated a new chamber which rose to a new spot. In addition, if the Reynolds number of the conduit was greater than approximately ten, an axisymmetric wavy instability appeared in the walls of the conduit and the conduit had to be tilted less before a new chamber was initiated. If shear under the Pacific plate has to tilt buoyant mantle plumes to as much as 60° to form the relatively regular island chains associated with hot spots, most of the shear would be found in a zone with a vertical extent of less than 200 km. 相似文献
20.
Masao Nakada 《Geophysical Journal International》1999,137(3):663-674
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. 相似文献