Stress and viscosity in the asthenosphere

Stresses and effective viscosities in the asthenosphere to a depth of 400 km are calculated on the basis of Weertmans “temperature method” i.e., on relating viscosity to the ratio of the temperature to the melting point (=homologous temperature). Some oceanic and continental geotherms and two melting point—depth curves, the dry pyrolite solidus and the forsterite₉₀ melting curve are used for the conversion of the homologous temperature to the effective viscosity. Two creep laws are considered, the linear, grain-size-dependent Nabarro—Herring (NH) creep law, and a power creep law, in which the creep rate is proportional to the third power of the stress. A plate tectonic model yields creep rates of 2 · 10⁻¹⁴ s⁻¹ for the oceanic and 3 · 10⁻¹⁵ s⁻¹ for the continental asthenosphere. These values are held constant for the calculations and may be valid for regions inside plates.The dry pyrolite mantle model results in high homologous temperatures in the asthenosphere below oceans (0.9), very low stresses (a few bars and lower) and shows a low viscosity “layer” of about 200-km thickness. Below continental shields the homologous temperature has a maximum value of 0.73, stresses are around 5–20 bar and the low-viscosity region is thicker and less pronounced than in the oceanic case. The Fo₉₀ mantle model generally gives lower homologous temperatures (maximum value below oceans beside active ridges 0.75). The stresses in the asthenosphere beneath oceans vary from a few bars to about 50 bar and below continents to about 100 bar. The low-viscosity region seems to reach great depths without forming a “channel”. The Figs. 1 and 2 show the approximate viscosity—depth distribution for the two mantle models under study.Assuming a completely dry mantle and a mean grain size of 5 mm, power law creep will be the dominating creep process in the asthenosphere. However, grains may grow in a high-temperature—low-stress regime (i.e., below younger oceans), an effect which will further diminish the influence of NH creep. In the upper 100–150 km of the earth some fluid phases may affect considerably creep processes.     

The phase diagram calcite-aragonite as derived from the crystallographic properties

A general method to calculate phase diagrams from the fundamental crystallographic properties has been suggested, using the calcite-aragonite reaction as an example. This procedure requires the measurement of the phonon spectra at different temperatures and pressures, calculation of a reasonable Grüneisen parameter to obtain C _p values from the C _v values and the consideration of the order—disorder phenomena and its contribution to the Gibb's energy. The remarkable agreement seen between the calculated diagram and that determined experimentally speaks for the reliability of this method.     

First principles molecular dynamics simulations of diopside (CaMgSi2O6) liquid to high pressure

We use first principles molecular dynamics simulations based on density functional theory in the local density approximation to investigate CaMgSi₂O₆ liquid over the entire mantle pressure regime. We find that the liquid structure becomes much more densely packed with increasing pressure, with the mean Si-O coordination number increasing nearly linearly with volume from fourfold near ambient pressure to sixfold at the base of the mantle. Fivefold Si-O coordination environments are most abundant at intermediate compression. The properties of Mg and Si coordination environments are nearly identical to those in MgSiO₃ liquid, whereas Ca is more highly coordinated with larger mean Ca-O bond length as compared with Mg. The density increases smoothly with increasing pressure over the entire range studied. The Grüneisen parameter increases by a factor of three on twofold compression. The density contrast between diopside composition liquid and the isochemical crystalline assemblage is less than 2% at the core mantle boundary, less than that in the case of MgSiO₃. Thermodynamic properties are described in terms of a liquid-state fundamental thermodynamic relation.     

Elastic and inelastic triggering of earthquakes in the North Anatolian Fault zone

Deformation models used to explain the triggering mechanism often assume pure elastic behaviour for the crust and upper mantle. In reality however, the mantle and possibly the lower crust behave viscoelastically, particularly over longer time scales. Consequently, the stress field of an earthquake is in general time-dependent. In addition, if the elastic stress increase were enough to trigger a later earthquake, this triggered event should occur instantaneously and not many years after the triggering event. Hence, it is appropriate to include inelastic behaviour when analysing stress transfer and earthquake interaction.In this work, we analyse a sequence of 10 magnitude M_s > 6.5 events along the North Anatolian Fault between 1939 and 1999 to study the evolution of the regional Coulomb stress field. We investigate the triggering of these events by stress transfer, taking viscoelastic relaxation into account. We evaluate the contribution of elastic stress changes, of post-seismic viscoelastic relaxation in the lower crust and mantle, and of steady tectonic loading to the total Coulomb stress field. We analyse the evolution of stress in the region under study, as well as on the rupture surfaces of the considered events and their epicentres. We study the state of the Coulomb stress field before the 1999 İzmit and Düzce earthquakes, as well as in the Marmara Sea region.In general, the Coulomb stress failure criterion offers a plausible explanation for the location of these events. However, we show that using a purely elastic model disregards an important part of the actual stress increase/decrease. In several cases, post-seismic relaxation effects are important and greater in magnitude than the stress changes due to steady tectonic loading. Consequently, viscoelastic relaxation should be considered in any study dealing with Coulomb stress changes.According to our study, and assuming that an important part of the rupture surface must be stressed for an earthquake to occur, the most likely value for the viscosity of the lower crust or mantle in this region is 5 · 10¹⁷–10¹⁸ Pa · s. Our results cannot rule out the possibility of other time-dependent processes involved in the triggering of the 1999 Düzce event. However, the stress increase due to viscoelastic relaxation brought 22% of the 1999 Düzce rupture area over the threshold value of Δσ_c ≥ 0.01 MPa (0.1 bar), and took the whole surface closer to failure by an average of 0.2 MPa. Finally, we argue that the Marmara Sea region is currently being loaded with positive Coulomb stresses at a much faster rate than would arise exclusively from steady tectonic loading on the North Anatolian Fault.     

Analysis of Izmit aftershocks 25 days before the November 12th 1999 Düzce earthquake, Turkey

We investigate spatial clustering of 2414 aftershocks along the Izmit M_w = 7.4 August 17, 1999 earthquake rupture zone. 25 days prior to the Düzce earthquake M_w = 7.2 (November 12, 1999), we analyze two spatial clusters, namely Sakarya (SC) and Karadere–Düzce (KDC). We determine the earthquake frequency–magnitude distribution (b-value) for both clusters. We find two high b-value zones in SC and one high b-value zone in KDC which are in agreement with large coseismic surface displacements along the Izmit rupture. The b-values are significantly lower at the eastern end of the Izmit rupture where the Düzce mainshock occurred. These low b-values at depth are correlated with low postseismic slip rate and positive Coloumb stress change along KDC. Since low b-values are hypothesized with high stress levels, we propose that at the depth of the Düzce hypocenter (12.5 km), earthquakes are triggered at higher stresses compared to shallower crustal earthquake. The decrease in b-value from the Karadere segment towards the Düzce Basin supports this low b-value high stress hypothesis at the eastern end of the Izmit rupture. Consequently, we detect three asperity regions which are correlated with high b-value zones along the Izmit rupture. According to aftershock distribution the half of the Düzce fault segment was active before the 12 November 1999 Düzce mainshock. This part is correlated with low b-values which mean high stress concentration in the Düzce Basin. This high density aftershock activity presumably helped to trigger the Düzce event (M_w = 7.2) after the Izmit M_w 7.4 mainshock.     

Stress field in the Ryukyu–Kyushu Wadati–Benioff zone by inversion of earthquake focal mechanisms

The study addresses the space distribution of the stress field in the Kyushu–Ryukyu Wadati–Benioff zone based on homogeneous data of earthquake focal mechanisms and the inverse technique by Gephart and Forsyth [J. Geophys. Res. 89 (1984) 9305]. The used data set consists of 148 Harvard CMT solutions and 22 earthquake focal mechanisms listed in previous studies. The stress field parameters are determined for 0–40, 41–100 and h>100 km depth ranges. The top 100-km layer of the Wadati–Benioff zone (WBZ) is characterized by strike normal maximum compression σ₁ and steeper than the slab minimum compression σ₃, the last indicating for unbalanced slab pull force. The Tokara channel ‘divides’ the subduction into two parts of different stress regime at depth greater than 100 km. To the south of the channel the slab is under slab parallel σ₁ and slab normal σ₃ while its northern part, beneath Kyushu, is under slab parallel extension and slab normal compression. The results of recent studies on the regional velocity structure and geochemistry of the volcanic lava indicate that the most plausible reason for the observed stress field difference below 100 km in the northern and rest part of the arc is the presence of hot low viscosity upper mantle west of Kyushu.The results of this study indicate that the forces involved in the contemporary subduction dynamics in the Ryukyu–Kyushu Wadati–Benioff zone are related to the convergence between the Philippine Sea Plate and the Eurasian plate, the trench suction force, slab pull, the slab anchor force and, in the southern-central part of the arc, mantle resistance.     

Progress report on rhenish massif seismic experiment

This paper presents some data and results from a seismic refraction experiment, completed mainly in 1979 in the Rhenish Massif, Federal Republic of Germany and extending through Luxembourg and Belgium into the Paris Basin in France.Velocity-depth functions have been derived for each record section independently, based on the assumption that velocity varies only with depth: these models are being improved upon by time-term and ray-tracing methods capable of handling laterally varying velocity structures and by calculating synthetic seismograms.The P_g phase which is observed very clearly on all record sections represents a refracted wave, with velocity generally > 6 km/s, from depths below 1.5–5.5 km. Along the 600 km long main profile one intracrustal reflection can usually be recognized, while from the three shorter crossing profiles in the massif two intracrustal reflectors can always be seen. Beneath much of the main profile the crust-mantle boundary is either a first order discontinuity or thin (< 1 km) transition zone at ~30 km depth. However, beneath the Ardennes and West Eifel there is a 6–12 km thick transition zone before a velocity of 8.1 km/s is reached at ~36 km depth. Beneath the crossing profiles, there is generally a transition zone < 3 km thick between crust and mantle. In some cases, there can be recognized at the top of the mantle a thin high velocity layer which is underlain by a low velocity layer which, in turn, is underlain by a reflector 4–11 km below the crust-mantle boundary.     

Imaging of Vp, Vs, and Poisson’s ratio anomalies beneath Kyushu, southwest Japan: Implications for volcanism and forearc mantle wedge serpentinization

We determine detailed 3-D V_p and V_s structures of the crust and uppermost mantle beneath the Kyushu Island, southwest Japan, using a large number of arrival times from local earthquakes. From the obtained V_p and V_s models, we further calculate Poisson’s ratio images beneath the study area. By using this large data set, we successfully image the 3-D seismic velocity and Poisson’s ratio structures beneath Kyushu down to a depth of 150 km with a more reliable spatial resolution than previous studies. Our results show very clear low V_p and low V_s anomalies in the crust and uppermost mantle beneath the northern volcanoes, such as Abu, Kujyu and Unzen. Low-velocity anomalies are seen in the mantle beneath most other volcanoes. In contrast, there are no significant low-velocity anomalies in the crust or in the upper mantle between Aso and Kirishima. The subducting Philippine Sea slab is imaged generally as a high-velocity anomaly down to a depth of 150 km with some patches of normal to low seismic wave velocities. The Poisson’s ratio is almost normal beneath most volcanoes. The crustal seismicity is distributed in both the high- and low-velocity zones, but most distinctly in the low Poisson’s ratio zone. A high Poisson’s ratio region is found in the forearc crustal wedge above the slab in the junction area with Shikoku and Honshu; this high Poisson’s ratio could be caused by fluid-filled cracks induced by dehydration from the Philippine Sea slab. The Poisson’s ratio is normal to low in the forearc mantle in middle-south Kyushu. This is consistent with the absence of low-frequency tremors, and may indicate that dehydration from the subducting crust is not vigorous in this region.     

Effects of buoyancy and mechanical layering on collisional deformation of continental lithosphere: Results from physical modeling

We use two suites of lithospheric-scale physical experiments to investigate the manner in which deformation of the continental lithosphere is affected by both (1) variations of lithospheric density (quantified by the net buoyant mass per area in the lithospheric mantle layer, M_B), and (2) the degree of coupling between the crust and lithospheric mantle (characterized by a modified Ampferer ratio, Am). The dynamics of the experiments can be characterized with a Rayleigh–Taylor type ratio, _CLM. Models with a positively buoyant lithospheric mantle layer (M_B > 0 and _CLM > 0) result in distributed root formation and a wide deformation belt. In contrast, models with a negatively buoyant lithospheric mantle layer strongly coupled to the crust (M_B < 0, 0 > _CLM > ≈ − 0.2, and Am > ≈ 10^{− 3}) exhibit localized roots and narrow deformation belts. Syncollisional delamination of the model lithospheric mantle layer and a wide deformation belt is exhibited in models with negatively buoyant lithospheric mantle layers weakly coupled to the crust (M_B < 0, _CLM < 0, and Am < ≈ 10^{− 3}). Syncollisional delamination of the continental lithosphere may initiate due to buoyancy contrasts within the continental plate, instead of resulting from wedging by the opposing plate. Rayleigh–Taylor instabilities dominate the style of deformation in models with a negatively buoyant lithospheric mantle layer strongly coupled to the crust and a slow convergence rate (M_B < 0 and _CLM > ≈ − 0.2). The degree of coupling (Am) between the model crust and lithospheric mantle plays a lesser role in both the style of lower-lithospheric deformation and the width of the crustal deformed zone with increasing density of the lithospheric mantle layer.     

Physical conditions producing slab stagnation: Constraints of the Clapeyron slope, mantle viscosity, trench retreat, and dip angles

Recent seismic tomography has revealed various morphologies in the subducted lithosphere. In particular, significant flattening and stagnation of slabs around the 660-km boundary are seen in some areas beneath the northwestern Pacific subduction zones. We examined the cause of slab stagnation in terms of the Clapeyron slope of the phase transformation from ringwoodite to perovskite + magnesiowüstite, trench retreat velocity, dip angles, and high viscosity of the lower mantle based on two-dimensional (2-D) numerical simulations of thermal convection. In particular, we examined the conditions necessary for slab stagnation assuming a very small absolute value of the Clapeyron slope, which were proposed based on recent high-pressure, high-temperature (high P–T) experiments. Our calculations show that slabs tend to stagnate above the 660-km boundary with an increasing absolute value of the Clapeyron slope, viscosity jump at the boundary, and trench retreat velocity and a decreasing initial dip angle. Stagnant slabs could be obtained numerically for a realistic range of parameters obtained from high P–T experiments and other geophysical observations combining buoyancy, high lower-mantle viscosity, and trench retreat. We found that a low dip angle of a descending slab at the bottom of the upper mantle plays an important role in slab stagnation. Two main regimes underlie slab stagnation: buoyancy-dominated and viscosity-dominated regimes. In the viscosity-dominated regime, it is possible for slabs to stagnate above the 660-km boundary, even when the value of the Clapeyron slope is 0 MPa/K.     

Deep structure of the baikal and other continental rift zones from seismic data

Abyssal variations beneath the Baikal rift zone are revealed in an irregular seismic stratification of the crust, the presence of an intracrust waveguide and by the vast (> 200,000 km²) underlying area of anomalously low velocity (P_n = 7.6−7.8 km/sec) uppermost mantle. In its abyssal structure the Baikal rift is heterogeneous along the strike, with sharp changes in crustal thickness (35–50 km).Comparison of first-arrival seismic-velocity curves and also the respective velocity columns reveals the essential similarity of upper-mantle seismic cross-sections for all continental rift zones. The anomalous upper layer of the mantle (ca. 7.7 km/sec) is relatively thin (15-13 km) and can be linked with the mantle waveguide only locally.     

Upper mantle structure of the Northern Eurasia from peaceful nuclear explosion data

Several long-range seismic profiles were carried out in Russia with Peaceful Nuclear Explosions (PNE). The data from 25 PNEs recorded along these profiles were used to compile a 3-D upper mantle velocity model for the central part of the Northern Eurasia. 2-D crust and upper mantle models were also constructed for all profiles using a common methodology for wavefield interpretation. Five basic boundaries were traced over the study area: N1 boundary (velocity level, V = 8.35 km/s; depth interval, D = 60–130 km), N2 (V = 8.4 km/s; D = 100–140 km), L (V = 8.5 km/s; D = 180–240 km) and H (V = 8.6 km/s; D = 300–330 km) and structural maps were compiled for each boundary. Together these boundaries describe a 3-D upper mantle model for northern Eurasia. A map characterised the velocity distribution in the uppermost mantle down to a depth of 60 km is also presented. Mostly horizontal inhomogeneity is observed in the uppermost mantle, and the velocities range from the average 8.0–8.1 km/s to 8.3–8.4 km/s in some blocks of the Siberian Craton. At a depth of 100–200 km, the local high velocity blocks disappear and only three large anomalies are observed: lower velocities in West Siberia and higher velocities in the East-European platform and in the central part of the Siberian Craton. In contrast, the depths to the H boundary are greater beneath the craton and lower beneath in the West Siberian Platform. A correlation between tectonics, geophysical fields and crustal structure is observed. In general, the old and cold cratons have higher velocities in the mantle than the young platforms with higher heat flows.Structural peculiarities of the upper mantle are difficult to describe in form of classical lithosphere–asthenosphere system. The asthenosphere cannot be traced from the seismic data; in contrary the lithosphere is suggested to be rheologically stratified. All the lithospheric boundaries are not simple discontinuities, they are heterogeneous (thin layering) zones which generate multiphase reflections. Many of them may be a result of fluids concentrated at some critical P–T conditions which produce rheologically weak zones. The most visible rheological variations are observed at depths of around 100 and 250 km.     

Fluid source-based modeling of melt initiation within the subduction zone mantle wedge: Implications for geochemical trends in arc lavas

The GyPSM-S (Geodynamic and Petrological Synthesis Model for Subduction) scheme couples a petrological model with a 2-D thermal and variable viscosity flow model to describe and compare fundamental processes occurring within the subduction mantle wedge, including the development of a low-viscosity channel (LVC) (Hebert et al., 2009, Earth and Planetary Science Letters, v. 278, p. 243–256). Here we supplement the basic coupled model result with more sophisticated treatments of trace element partitioning in the fluid phase and melt transport regimes. We investigate the influences of slab fluid source lithology and fluid transport mechanisms on melt geochemistry, the implications of mantle source depletion related to fluid fluxing, and potential melt migration processes. This study describes two model cases that can be compared to geochemical datasets for the Izu–Bonin intra-oceanic subduction system and the Central Costa Rican part of the Central American arc. We find that there is a progression of geochemical characteristics described in studies of cross-arc and along-arc lavas that can be approximated assuming (i) limited fluid–rock interaction within the mantle wedge and (ii) that melt migration preserves the spatial distinction among melts initiated in different areas of the wedge. Specifically, volcanic front lavas have significant contributions from shallower slab fluid sources, and rear-arc lavas have significant contributions from deeper slab fluid sources. Evidence for limited fluid–rock interaction could imply either a rapid fluid transport mechanism or a fluid-dominated trace element budget within the LVC. Although we do not include a back-arc in these models, interpretations of the results lead to several potential mechanisms to explain hydrous inputs to back-arc source regions.     

Continental rifting as a function of lithosphere mantle strength

The role of the uppermost mantle strength in the pattern of lithosphere rifting is investigated using a thermo-mechanical finite-element code. In the lithosphere, the mantle/crust strength ratio (S_M/S_C) that decreases with increasing Moho temperature T_M allows two strength regimes to be defined: mantle dominated (S_M > S_C) and crust dominated (S_M < S_C). The transition between the two regimes corresponds to the disappearance of a high strength uppermost mantle for T_M > 700 °C. 2D numerical simulations for different values of S_M/S_C show how the uppermost mantle strength controls the style of continental rifting. A high strength mantle leads to strain localisation at lithosphere scale, with two main patterns of narrow rifting: “coupled crust–mantle” at the lowest T_M values and “deep crustal décollement” for increasing T_M values, typical of some continental rifts and non-volcanic passive margins. The absence of a high strength mantle leads to distributed deformations and wide rifting in the upper crust. These numerical results are compared and discussed in relation with series of classical rift examples.     

Influence of a low-viscosity layer between rigid inclusion and viscous matrix on inclusion rotation and matrix flow: A numerical study

We have used 2-D finite element modelling to investigate the influence of a permanent low-viscosity layer between matrix and inclusion on matrix flow and inclusion rotation under viscous simple shear flow. Rigid inclusions of different shape (circle, square, ellipse, lozenge, rectangle and skewed rectangles) and aspect ratio (R) were used. The calculated matrix flow pattern is neither bow tie nor eye-shaped. It is a new flow pattern that we call cat eyes-shaped, which is characterized by: (i) straight streamlines that slightly bend inwards at the inclusion's crests; (ii) elongate eye-shaped streamlines on each side of the inclusion; (iii) stagnation points in the centre of the eyes; (iv) absence of closed streamlines surrounding the inclusion; (v) changes in flow configuration with inclusion orientation; the lines of flow reversal bend and tilt, closed streamline circuits may disappear, and streamlines may bend outwards at the inclusion's crests.Concerning inclusion rotation, the numerical results show that: (i) a low-viscosity layer (LVL) makes inclusions with R = 1 rotate synthetically, but the rotation rate depends upon shape (circle or square) and orientation. Therefore, shape matters in the slipping mode. (ii) All studied shapes with R > 1 rotate antithetically when starting with the greatest principal axis (e₁) parallel to the shear direction ( = 0°); (iii) rotation is limited because there is a stable equilibrium orientation (_se) for all studied shapes with R > 1. (iii) There is also an unstable equilibrium orientation (_ue), and both _se and _ue depend upon inclusion's R and shape.The present numerical results closely agree with previous results of analogue experiments with a permanent low viscosity interface. Only minor deviations related with small shape differences were detected.     

Structure, thermodynamics, and diffusion in CaAl2Si2O8 liquid from first-principles molecular dynamics

We have performed first-principles molecular dynamics simulations of CaAl₂Si₂O₈ (anorthite) liquid at pressures up to 120 GPa and temperatures of 3000, 4000 and 6000 K. At the lowest degrees of compression the liquid is seen to accommodate changes in density through decreasing the abundance of 3- and 4-membered rings, while increases in coordination of network forming cations take effect at somewhat higher degrees of compression. Results are fit to a fundamental thermodynamic relation with 4th order finite strain and 1st order thermal variable expansions. Upon compression by a factor of two, the Grüneisen parameter (γ) is found to increase continuously from 0.35 to 1.10. Weak temperature dependence in γ is thermodynamically consistent with a slight decrease in isochoric heat capacity (C_V), for which values of between 4.4 and 5.2 Nk_B are obtained, depending on the temperature. Pressure and temperature dependence of self-diffusivities is found to be well represented by an Arrhenius relation, except at 3000 K and pressures lower than 5 GPa, where self-diffusivities of Si, Al, and O increase with pressure. Analysis of the lifetimes of individual coordination species reveals that this phenomenon arises due to the disproportionately high stability of 4-fold coordinated Si, and to a lesser extent 4-fold coordinated Al. Our results represent a marked improvement in accuracy and reliability in describing the physics of CaAl₂Si₂O₈ liquid at deep mantle pressures, pointing the way to a general thermodynamic model of melts at extreme pressures and temperatures relevant to planetary-scale magma oceans and deep mantle partial melting.     

Phase transformations in descending plates and implications for mantle dynamics

The role of phase transformations in a mantle of pyrolite composition is reviewed in the light of recent experimental data. The pyroxene component of pyrolite transforms to the garnet structure at 300–350 km whilst olivine transforms to beta-Mg₂SiO₄ near 400 km. Between about 500 and 550 km, beta-Mg₂SiO₄ probably transforms to a partially inverse spinel structure whilst the CaSiO₃ component of the complex garnet solid solution exsolves and transforms to the perovskite structure. The major seismic discontinuity near 650–700 km is probably caused by disproportionation of Mg₂SiO₄ spinel into periclase plus stishovite. At a slightly greater depth, the remaining magnesian garnet transforms to the corundum or ilmenite structure. Finally, at a depth probably in the vicinity of 800–1000 km, the (Mg,Fe)SiO₃ component of the ilmenite phase transforms to a perovskite structure whilst stishovite and some of the periclase recombine to form perovskite also. The mineral assemblage so formed is about 4% denser than mixed oxides (MgO + FeO + A1₂O₃ + CaO + stishovite) isochemical with pyrolite. The above sequence of phase transformations in pyrolite provides a satisfactory general explanation of the elastic properties and density distribution in the mantle. In particular, there is no evidence requiring an increase of FeO/(FeO + MgO) ratio with depth.The depths at which major phase transformations occur in subducted lithosphere differ from those in ‘normal’ mantle. These differences are caused by two factors: (1) Temperatures within sinking plates are much lower than in surrounding mantle to depths of 700 km or more. (2) Irreversible chemical differentiation of pyrolite occurs at oceanic ridges. Lithosphere plates so formed consist of a layer of basaltic rocks underlain successively by layers of harzburgite, lherzolite, and pyrolite slightly depleted in highly incompatible elements (e.g. La, Ba, Rb, U). The phase-transformation behaviour of the first three of these layers differs from that of pyrolite.The effects of these and other factors connected with phase transformations on the dynamics of plate subsidence are discussed. It appears quite likely that plates penetrate the 650–700 km discontinuity, largely because the slope of the spinel disproportionation is probably positive, not negative as generally supposed. The former basaltic oceanic crust probably sinks deeply into the lower mantle, whilst the former harzburgite component of the plate may collect above the perovskite transition boundary. Phase transformations may thus serve as a kind of filter, leading to increased and irreversible mantle heterogeneity with time.The possible roles of phase transformations in causing deep-focus earthquakes and introducing water into the mantle in subduction zones are also briefly discussed.     

A two-phase model for the description of the influence of cracks on the P- and S-wave velocities in dry and saturated rock samples

The premonitory variations of seismic-wave velocities before earthquakes originate from various cracking processes before the fracture. It can be shown that these variations are comparable with wave-velocity variations in porous model samples with defined pore sizes. Furthermore, it is possible to describe analytically the variation of wave velocities as a function of the parameter K₀, which describes the fracturing process, and a material/depth parameter AP. On the basis of the wave velocity vs. pressure curves of rocks, it is possible to determine K₀ and A. Using the material/depth parameter AP sediments in covering strata and eruptive rocks in regions of earthquakes of shallow to medium depth can be descirbed. A relationship between Δυ_D and ΔK₀ can be established. For acid to ultrabasic rocks, a variation of up to 2 km/s indicates a variation of K₀ of 0.1–1.0.Moreover, it is possible to establish a relationship between K₀ and the number as well as the mean length of cracks in the rock. The solutions differ depending on the cracks being closed or open. For closed cracks a wave-velocity minimum of 6% results. For oper cracks the variations of the number and mean length of cracks are taken into account by means of a stochastic process; the resulting variations of K₀ and the wave velocity car adequately explain the variations in seismic-wave velocity. Variations of the pore pressure have an influence on K₀and the wave velocity only under most favourable geologica conditions; generally they are insignificant. Also for S-wavc velocities and for the ratio υ_p/υ_s the wave velocity vs. pressure equations are valid; it is possible to state K₀ and A-values.     

Geochemistry of Cretaceous mafic rocks from the Lower Yangtze region, eastern China: Characteristics and evolution of the lithospheric mantle

Geochemical and isotopic data for Cretaceous mafic rocks (basalt, gabbro, and diorite) from the Lower Yangtze region, northern Yangtze block, constrain the evolution of the lithospheric mantle. The mafic rocks, separated into the northeast and southwest groups, are alkaline and evolved, with low Mg# values (44–58) and variable SiO₂ contents (47.6–57.4 wt%). Enriched LREEs, LILEs, and Pb, together with depleted Nb, Zr, and Ti, suggest that the mantle sources were metasomatized by slab-derived fluid/melt. All samples show high radiogenic ²⁰⁷Pb/²⁰⁴Pb(t) (15.41–15.65) and ²⁰⁸Pb/²⁰⁴Pb(t) (37.66–38.51) ratios at given ²⁰⁶Pb/²⁰⁴Pb(t) (17.65–19.00) ratios, consistent with the mantle sources having been metasomatized by ancient slab-derived material. Mafic rocks of the southwest group show enriched Sr–Nd isotopic characteristics, with ⁸⁷Sr/⁸⁶Sr(t) ranging from 0.7056 to 0.7071 and ε_Nd(t) ranging from −5.3 to −8.3, indicating an origin from enriched lithospheric mantle. Mafic rocks of the northeast group, which record ⁸⁷Sr/⁸⁶Sr(t) ratios of between 0.7044 and 0.7050 and ε_Nd(t) of −2.8 to −0.7, possibly formed by the mixing of melts from isotopically enriched lithospheric mantle and isotopically depleted asthenospheric mantle. Taking into consideration the geochemical and isotopic characteristics of Cretaceous mafic rocks, Cenozoic basalts, and basalt-hosted peridotite xenoliths from the Lower Yangtze region, we propose that an isotopically enriched, subduction-modified lithospheric mantle was replaced by or transformed into an isotopically depleted “oceanic-type” mantle. Such a process appears to have occurred in the eastern North China Craton as well as the eastern Yangtze block, probably in response to subduction of the paleo-Pacific plate beneath East Asia.     

On the possibility of Newtonian flow in the upper mantle

The possibility of a linear-creep (Newtonian-viscosity) upper mantle is reexamined on the basis of present knowledge on flow mechanisms in olivine, including differences between activation parameters for creep and for diffusion, and revised estimates of grain boundary width. Results of the comparison between linear superplastic creep and power-law creep are presented as crossover temperature between lattice and grain-boundary diffusion, and crossover stress between non-linear and linear creep, as a function of temperature (depth), grain size, grain-boundary width, grain-boundary diffusion activation energy, and rate-controlling species (silicon or oxygen). For the most realistic values of the parameters, linear creep is well within the range of possibilities. There is no major objection from rheology to the idea of a Newtonian-viscosity upper mantle.Viscosities for the two kinds of creep are compared for laboratory, upper mantle, and—tentatively—lower mantle conditions. This results in a prediction of the grain size at which Coble creep could be observed in olivine polycrystals in the laboratory (~ 10 μm or less). The upper mantle viscosity is in the 10²⁰–10²¹ Pa s range. Minimum estimates for the viscosity of the lower mantle are of the same order.The pictures of the rheology of the mantle derived from microphysical models of flow and from geophysical observation can therefore be made compatible.