首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Summary. Quartz eclogite as source rock seems capable of explaining the chemistry of calc-alkaline volcanics. This model requires partial melting of quartz eclogite in the depth range 100–200 km. Shear heating of the subducting crust is one of the most debated models. The problem may be treated with two different boundary conditions: constant strain rate or constant shear stress. In the former case an increase in temperature tends to reduce the heat production thus stabilizing the shear flow; the result is a moderate increase in temperature which will remain constant during flow. The latter case may lead to a thermal feedback instability and thus temperature and strain rate may suddenly grow to very high values. This phenomenon is termed a runaway and will be discussed in this paper using an adiabatic approximation. It is shown that for a runaway to occur the local energy density must amount to a common value independent of the rheology. In contrast to the constant strain rate case, shear heating is negligible until just before the instability occurs. When the melting point is reached shear stress will break down but the stored local energy will be set free and supply the latent heat of melting. The possibility of a runaway occurring is strongly dependent on the ambient temperature. In subduction shear zones neither shear stress nor strain rate are likely to be constant throughout, but if the former is constant or changes little on a 1–10 km scale a runaway is liable to occur at a depth of around 150 km thus possibly being the cause of calcalkaline volcanism.  相似文献   

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

3.
A lower mantle S-wave triplication and the shear velocity structure of D"   总被引：6，自引：0，他引：6  
Summary. A lower mantle S-wave triplication detected with short- and long-period WWSSN and CSN recordings indicates a substantial shear velocity discontinuity near 280 km above the core–mantle boundary. The triplication can be observed in rotated SH seismograms from intermediate and deep focus events throughout the distance range from 70° to 95°. Three distinct source region–receiver array combinations that have been investigated in detail demonstrate consistent travel time and relative amplitude behaviour of the triplication, with slight systematic shifts in the triplication indicating up to 40 km variations in the depth of the discontinuity. Modelling of the observations with synthetic seismograms produced with the Cagniard de Hoop and reflectivity methods constrains the shear velocity increase to be 235 ± 0.25 per cent, comparable to upper mantle discontinuities. Short-period observations indicate that the velocity increase may be a sharp first-order discontinuity, or may extend over a transition zone no more than 50 km thick. The shear velocity gradient below the discontinuity, within the D" layer, is not well-constrained by the SH data, but slightly positive or near zero velocity gradients are consistent with the long-period amplitude ratios of ScSH/SH .  相似文献   

4.
Summary. Most crustal earthquakes of the world are observed to occur within a seismogenic layer which extends from the Earth's surface to a depth of a few tens of kilometres at most. A model is proposed in which the shear zone along a transcurrent plate margin is represented as a viscoelastic medium with depth-dependent power-law rheology. A frictional resistance linearly increasing with depth is assumed on a vertical transcurrent fault within the shear zone. Such a model is able to reproduce a continuous transition from the brittle behaviour of the upper crust to the ductile behaviour at depth. Assuming that the shear zone is subjected to a constant strain rate from the opposite motions of the two adjacent plates, it is found that there exists a maximum depth H below which tectonic stress can never reach the frictional threshold: this may be identified as the maximum depth of earthquake nucleation. The value of H is consistent with observations for plausible values of the model parameters. The stress evolution in the shear zone is calculated in the linear approximation of the constitutive equation. A change in rigidity with depth, which is also introduced in the model, may reproduce the high vertical gradient of shear stress, which has been measured across the San Andreas fault, and the fact that most earthquakes are nucleated at some depth in the seismogenic layer. A crack which drops the ambient stress to the dynamic frictional level is then introduced in the model. To this aim, a crack solution is employed without a stress singularity at its edges, which is compatible with a frictional stress threshold criterion for fracture. A constraint on the vertical friction gradient is obtained if such cracks are assumed to be entirely confined within the seismogenic layer.  相似文献   

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. Fault zones in wet Westerly granite deformed at temperatures of 300° and 400°C require markedly lower shear stresses for sliding than when dry, even when the effective confining pressure is held constant between the wet and dry tests, provided that the strain rate is lower than 10−7s−1. The rate of strength reduction is enhanced by increasing the pore water pressure. The deformation rate is a power function of the applied stress where the stress exponent is approximately 7 for pore water pressure of 100 MPa and 21 for pore water pressure of 20 MPa.
The experimental results are extrapolated to conditions believed to occur at depths of 10 km along the San Andreas Fault Zone. It is suggested that for slow tectonic deformation at strain rates of 10−11 and 10−14s−1 the shear stress for sliding on faults in granite would be approximately 60 and 20 MPa, respectively, at pore water pressures equal to the hydrostatic head. Fluid overpressures of c. 0.8 lithostatic pressure are required to lower the shear stress for sliding into the 10 MPa range at the slower strain rate.  相似文献   

7.
Summary. This paper extends an earlier study (Sengupta & Julian) of travel times of P waves of deep-focus earthquakes to include shear waves. Primary advantage of deep-focus earthquakes is the reduction of anomalies caused by complex structures near the source. The standard deviations of travel times and station anomalies of this study are about half as large as those determined from the data of shallow-focus earthquakes (e.g. Herrin et al.; Hales & Roberts). Spherically-symmetric velocity models derived from the travel times by a linearized inverse technique have resolving lengths of about 70 km for standard errors in velocity of about 0.02 km/s. No pronounced reversal of either compressional or shear velocity was required at the base of the mantle to satisfy the data, though a small velocity decrease could not be entirely ruled out. Some anomalous rapid changes in compressional velocity gradient were, however, found centred around the depths of 2400 and 2600 km. The models derived in this study agree most closely with that of Herrin et al . for compressional velocity and the model 1066B of Gilbert & Dziewonski for shear velocity.  相似文献   

8.
We develop an approach that allows us to invert for the mantle velocity structure within a finely parametrized region as a perturbation with respect to a low-resolution, global tomographic model. We implement this technique to investigate the upper-mantle structure beneath Eurasia and present a new model of shear wave velocity, parametrized laterally using spherical splines with ∼2.9° spacing in Eurasia and ∼11.5° spacing elsewhere. The model is obtained from a combined data set of surface wave phase velocities, long-period waveforms and body-wave traveltimes. We identify many features as narrow as few hundred kilometres in diameter, such as subducting slabs in eastern Eurasia and slow-velocity anomalies beneath tectonically active regions. In contrast to regional studies in which these features have been identified, our model encompasses the structure of the entire Eurasian continent. Furthermore, including mantle- and body-wave waveforms helped us constrain structures at depths larger than 250 km, which are poorly resolved in earlier models. We find that up to +9 per cent faster-than-average anomalies within the uppermost ∼200 km of the mantle beneath cratons and some orogenic regions are separated by a sharp gradient zone from deeper, +1 to +2 per cent anomalies. We speculate that this gradient zone may represent a boundary separating the lithosphere from the continental root, which might be compositionally distinct from the overlying lithosphere and remain stable either due to its compositional buoyancy or due to higher viscosity compared with the suboceanic mantle. Our regional model of anisotropy is not significantly different from the global one.  相似文献   

9.
Summary. In order to examine the development of the oceanic crust in the neighbourhood of a slowly spreading ridge, a seismic refraction experiment was carried out at 59° 30'N on the Reykjanes Ridge. Three 120 km long overlapped split profiles were shot parallel to the trend of the ridge, on the eastern flank, and recorded on up to five recording sonobuoys. The profiles were at distances of 0, 30 and 90km from the ridge axis, corresponding to approximate crustal ages of 0, 3 and 9 Myr. Data from the main profiles were supplemented by using a large chamber air gun during recovery of the buoys.
The analysis of the data combined standard travel-time interpretation, the 'tau' method of systematic travel-time inversion and detailed amplitude modelling using the Reflectivity Method to calculate synthetic seismograms. Detailed velocity-depth models were constructed for each of the profiles.
There is no indication of a significant magma chamber at the ridge crest, although a slight velocity inversion in layer 3 suggests a zone of elevated temperature. Away from the crest there was a slight positive velocity gradient in layer 3. Layer 2 was most effectively modelled by a region of varying velocity gradients, which thinned with age and the transition to layer 3 is marked by a sharp change in velocity gradient. The transition to mantle velocities is also best modelled by a high velocity gradient rather than an interface.
Although some lateral variation in properties is apparent along the profiles, the lateral velocity gradients were sufficiently weak to allow an effective analysis in terms of laterally uniform models.  相似文献   

10.
Summary. The seismic structure has been measured to a depth of about 3 km along a 30 km seismic profile in east central Ireland. This profile is unusual in that it is the S -wave velocity—depth structure that has been measured to a degree of precision more normally associated with P -wave results. One reason for this is that the sources used were quarry blasts which generated strong S -waves and short-period surface waves but rather weak P -waves.
The results show a layer of Carboniferous limestone with shear velocity 2.65 km−1 s overlying a layer with a velocity of 3.06 km s−1. This second layer was interpreted as Lower Palaeozoic strata (Silurian/Ordovician) since this velocity was evident in an inlier seen at the surface at the northern end of the line. A third refraction horizon, shear velocity 3.45 km s−1 and displaying a basinal structure, was also recognized. This may be Cambrian or Precambrian basement.  相似文献   

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

12.
Summary. The ascent of a hot spherical body through a fluid with a strongly temperature-dependent viscosity has been studied using an axisymmetric finite element method. Numerical solutions range over Peclet numbers of 10−1– 103 from constant viscosity up to viscosity variations of 105. Both rigid and stress-free boundary conditions were applied at the surface of the sphere. The dependence of drag on viscosity variation was shown to have no dependence on the stress boundary condition except for a Stokes flow scaling factor. A Nusselt number parameterization based on the stress-free constant viscosity functional dependence on the Peclet number scaled by a parameter depending on the viscosity structure fits both stress-free and rigid boundary condition data above viscosity variations of 100. The temperature scale height was determined as a function of sphere radius. For the simple physical model studied in this paper pre-heating is required to reduce the ambient viscosity of the country rock to less than 1022 cm2 s−1 in order for a 10 km diapir to penetrate a distance of several radii.  相似文献   

13.
Upper mantle shear structure of North America   总被引：5，自引：0，他引：5  
Summary. The waveforms and travel times of S and SS phases in the range 10°–60° have been used to derive upper mantle shear velocity structures for two distinct tectonic provinces in North America. Data from earthquakes on the East Pacific Rise recorded at stations in western North America were used to derive a tectonic upper mantle model. Events on the north-west coast of North America and earthquakes off the coast of Greenland provided the data to investigate the upper mantle under the Canadian shield. All branches from the triplications due to velocity jumps near 400 and 660 km were observed in both areas. Using synthetic seismograms to model these observations placed tight constraints on heterogeneity in the upper mantle and on the details of its structure. SS–S travel-time differences of 30 s along with consistent differences in waveforms between the two data sets require substantial heterogeneity to at least 350 km depth. Velocities in the upper 170 km of the shield are about 10 per cent higher than in the tectonic area. At 250 km depth the shield velocities are still greater by about 4.5 per cent and they gradually merge near 400 km. Below 400 km no evidence for heterogeneity was found. The two models both have first-order discontinuities of 4.5 per cent at 405 km and 7.5 per cent at 695 km. Both models also have lids with lower velocities beneath. In the western model the lid is very thin and of relatively low velocity. In the shield the lid is 170 km thick with very high elocity (4.78 km s-1); below it the velocity decreases to about 4.65 km s-1. Aside from these features the models are relatively smooth, the major difference between them being a larger gradient in the tectonic region from 200 to 400 km.  相似文献   

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

15.
Summary. Lateral heterogeneity exists in the Earth's mantle, and may result in seismic velocity anomalies up to several per cent. If convection cells and plumes extend down to the core, then these features may be associated with local inhomogeneities observed in the lower mantle.
Published data for direct and core-reflected P -wave residuals are used to delineate velocity anomalies in the middle—lower mantle under the North Atlantic. Differential ( PcP — P ) residuals indicate travel-time anomalies near the core—mantle transition, and may be due to core topography or lateral variations in velocity. It is assumed that the anomalies occur near the midpoints of the ray paths. The main source of error in the data set may arise from phases which have been identified incorrectly. Hence trend-surfaces are fitted to the residual data to show only the large-scale trends in anomaly values, with wavelengths of the order of 1000 km.
The Azores and Colorado hot spots occur in a region covered by the data. A possible interpretation of the trend maps is that an anomalous zone extends from a relatively fast region at the core boundary at 35° N, 50° W up to these hot spots, at about 30 degrees from the vertical. This may agree with the suggestion of Anderson that plumes are chemical rather than thermal in origin. If inclined plumes do exist, the deviation from the ideal vertical plume or convection cell boundary may imply that lateral shear or other distortion effects exist in the mantle.  相似文献   

16.
The earthquakes in the seismicity belt extending through Indonesia, New Guinea, Vanuatu and Fiji to the Tonga–Kermadec subduction zone recorded at the 65 portable broad-band stations deployed during the Skippy experiment from 1993–1996 provide good coverage of the lithosphere and mantle under the Australian continent, Coral Sea and Tasman Sea.
The variation in structure in the upper part of the mantle is characterized by deter-mining a suite of 1-D structures from stacked record sections utilizing clear P and S arrivals, prepared for all propagation paths lying within a 10° azimuth band. The azimuth of these bands is rotated by 20° steps with four parallel corridors for each azimuth. This gives 26 separate azimuthal corridors for which 15 independent 1-D seismic velocity structures have been derived, which show significant variation in P and S structure.
The set of 1-D structures is combined to produce a 3-D representation by projecting the velocity values along the ray path using a turning point approximation and stacking into 3-D cells (5° by 50 km in depth). Even though this procedure will tend to underestimate wave-speed perturbations, S -velocity deviations from the ak135 reference model exceed 6 per cent in the lithosphere.
In the uppermost mantle the results display complex features and very high S -wave speeds beneath the Precambrian shields with a significant low-velocity zone beneath. High velocities are also found towards the base of the transition zone, with high S -wave speeds beneath the continent and high P -wave speeds beneath the ocean. The wave-speed patterns agree well with independent surface wave studies and delay time tomography studies in the zones of common coverage.  相似文献   

17.
We describe results of an active-source seismology experiment across the Chilean subduction zone at 38.2°S. The seismic sections clearly show the subducted Nazca plate with varying reflectivity. Below the coast the plate interface occurs at 25 km depth as the sharp lower boundary of a 2–5 km thick, highly reflective region, which we interpret as the subduction channel, that is, a zone of subducted material with a velocity gradient with respect to the upper and lower plate. Further downdip along the seismogenic coupling zone the reflectivity decreases in the area of the presumed 1960 Valdivia hypocentre. The plate interface itself can be traced further down to depths of 50–60 km below the Central Valley. We observe strong reflectivity at the plate interface as well as in the continental mantle wedge. The sections also show a segmented forearc crust in the overriding South American plate. Major features in the accretionary wedge, such as the Lanalhue fault zone, can be identified. At the eastern end of the profile a bright west-dipping reflector lies perpendicular to the plate interface and may be linked to the volcanic arc.  相似文献   

18.
Summary. A tripartite ocean-bottom seismograph array at the junction of the East Pacific Rise and Rivera Fracture Zone recorded an eathquake sequence, consisting of three main shocks ( m B= 4.3, 4.3 and 4.8) and numerous aftershocks from the fracture zone, in the distance range 35–50 km. Delineation of the rupture zones by aftershocks indicates that the first two main shocks took place on overlapping fault areas, while the third occurred over a fault area separated from the first by several kilometres. Both rupture zones were about 4 km long. Surface wave spectra indicate a shallow (about 3 km below the sea floor) source, as does OBS array phase velocity data. The seismic moments, obtained from teleseismic surface wave data, of 1.3, 2.1 and 2.8 × 1023 dyn cm, with the fault areas as delineated by aftershocks, imply a stress drop of about 8 bars for the main shocks. Aftershock sequences of each of the main shocks are similar, with a b -value of about 0.65. Teleseismic P travel times are similar to those from near-surface sources in Nevada.  相似文献   

19.
Summary. Monthly sea-levels from an extensive array of North Atlantic tide gauges (26-50N) are examined. The spatial scale of the sea-level variations, and the reasons for them, are discussed; one application of such a study is clearly in the design of a tide gauge network for monitoring eustatic changes of sea-level.
The spatial scale of the sea-level changes is large. There is a coherent sea-level signal which can be traced along the eastern boundary of the North Atlantic from Newlyn (50N) to Tenerife (28N). There are also two distinct groupings of tide gauges along the western boundary, separated by Cape Hatteras.
The contribution of local air pressure and wind stress is quantified at each gauge through multiple regression techniques and the gains are then interpreted in terms of recent theoretical and numerical modelling studies. For example, the gains suggest that the wind-forced boundary current along the Nova Scotian shelf is trapped to within about 16 km of the coast.
The influence of local meteorology cannot account for the large-scale modes of variability. The coherent signal along the eastern boundary is correlated with changes in the Sverdrup transport of the North Atlantic and hence the large-scale wind field. The two modes on the western boundary appear to be related to baroclinic boundary current variations.
The Newlyn sea-level record is finally 'corrected' for some of the above effects to illustrate the utility of such a residual series in the identification of eustatic changes and vertical crustal movement.  相似文献   

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

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

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