Stress in the lithosphere: Inferences from steady state flow of rocks

Mechanical data and flow processes from steady state deformation experiments may be used to infer the state of stress in the lithosphere and asthenosphere. Extrapolations of flow equations to a representative geologic strain rate of 10^–14/sec. for halite, marble, quartzite, dolomite, dunite and enstatolite are now warranted because the steady state flow processes in the experiments are identical to those in rocks and because the geotherms are reasonably well established. More direct estimates are obtained from free dislocation densities, subgrain sizes and recrystallized grain sizes all of which are functions only of stress. Using the last of these techniques, we have estimated stress profiles as a function of depth from xenoliths in basalts and kimberlites, whose depths of equilibration were determined by pyroxene techniques, from four different areas of subcontinental and suboceanic upper mantle. The results are similar and indicate stress differences of about 200 to 300 bars at 40 to 50 km, decaying to a few tens of bars at depths betow 100 km. These stresses are reasonable and are in accord with extrapolations of the mechanical data provided that allowance is made for a general increase in strain rate and decrease in viscosity with depth.     

Distribution and geochronology of Oregon Plateau (U.S.A.) flood basalt volcanism: The Steens Basalt revisited

The timing and petrogenesis of mid-Miocene flood basalt volcanism in the northwest United States has been extensively addressed, yet the chemical characteristics and temporal details of the Steens Basalt, exposed on the Oregon Plateau, are poorly defined. Steens Basalt volcanism has generally been accepted to have occurred at ∼ 16.6 Ma, coeval and/or just prior to the onset of Columbia River Basalt Group volcanism to the north. New major and trace element analyses and nine ⁴⁰Ar/³⁹Ar ages ranging from 15.51 ± 0.28 to 16.58 ± 0.18 Ma were obtained on Oregon Plateau flood basalt lava flows from stratigraphic sections in close proximity to Steens Mountain. Additionally, new ⁴⁰Ar/³⁹Ar ages were obtained on the uppermost and thirty-first lava flow down from the top of the ∼ 1 km section of Steens Basalt exposed at Steens Mountain and yield eruption ages of 16.59 ± 0.10 and 16.55 ±0.10 Ma, respectively. Field relations between these basalt sections suggest that multiple eruptive centers were present in the vicinity of Steens Mountain.     

Characterization of a mud deposit offshore of the Patos Lagoon,southern Brazil

Rapid deposition of mud on the beach along the shoreface of Rio Grande do Sul, Brazil dramatically influences the normal operations in the littoral zone. In the surf zone, fluid and suspended mud opposes water-wave movement and dissipates water-wave energy; on the beach, mud limits trafficability. As part of a multinational, multidisciplinary program to evaluate the influence of mud strength, density and viscosity on water-wave attenuation, sediments were evaluated in situ or collected for evaluation from an area offshore of Cassino Beach, slightly south of the Patos Lagoon mouth. Shear strength of deposited sediments ranged from 0.6 kPa at the seafloor to 3.4 kPa at ∼1 m below the seafloor. Mud sediments were also collected to simulate the in situ response of fluid mud to shear stresses. For this determination, rheological evaluations were made using a strain-controlled Couette viscometer on numerous remixed samples that ranged in density from 1.05 to 1.30 g/cm³. It was determined that this mud is a non-ideal Bingham material in that it has a true initial yield stress as well as a upper Bingham yield stress. Initial yield stress ranged from 0.59 to 2.62 Pa, upper Bingham yield stress ranged from 1.05 to 7.6 Pa. Apparent viscosity ranged from 0.02 to 4.7 Pa s with the highest viscosities occurring between the two yield stresses. Sediment strength in the remixed samples is 2 to 3 orders of magnitude lower than the horizontal shear strength of the sediment bed as determined by shear vane or predicted from penetrometer measurements. This difference is partially due to the fact that rheological evaluations are made on fully remixed sediments, whereas horizontal shear strength is determined within relatively undisturbed sediments. Similar values of viscosity and shear strength are comparable to those determined for mud in other coastal areas where fluid mud deposits occur.     

Eruption and emplacement of a basaltic welded ignimbrite during caldera formation on Gran Canaria

The 14.1 Ma old composite ignimbrite cooling unit P1 (45 km³) on Gran Canaria comprises a lower mixed rhyolite-trachyte tuff, a central rhyolite-basalt mixed tuff, and a slightly rhyolite-contaminated basaltic tuff at the top. The basaltic tuff is compositionally zoned with (a) an upward change in basalt composition to higher MgO content (4.3–5.2 wt.%), (b) variably admixed rhyolite or trachyte (commonly <5 wt.%), and (c) an upward increasing abundance of basaltic and plutonic lithic fragments and cognate cumulate fragments. The basaltic tuff is divided into three structural units: (I) the welded basaltic ignimbrite, which forms the thickest part (c. 95 vol.%) and is the main subject of the present paper; (II) poorly consolidated massive, bomb- and block-rich beds interpreted as phreatomagmatic pyroclastic flow deposits; and (III) various facies of reworked basaltic tuff. Tuff unit I is a basaltic ignimbrite rather than a lava flow because of the absence of top and bottom breccias, radial sheet-like distribution around the central Tejeda caldera, thickening in valleys but also covering higher ground, and local erosion of the underlying P1 ash. A gradual transition from dense rock in the interior to ash at the top of the basaltic ignimbrite reflects a decrease in welding; the shape of the welding profile is typical for emplacement temperatures well above the minimum welding temperature. A similar transition occurs at the base where the ignimbrite was emplaced on cold ground in distal sections. In proximal sections the base is dense where it was emplaced on hot felsic P1 tuff. The intensity of welding, especially at the base, and the presence of spherical particles and of mantled and composite particles formed by accretion and coalescence in a viscous state imply that the flow was a suspension of hot magma droplets. The flow most likely had to be density stratified and highly turbulent to prevent massive coalescence and collapse. Model calculations suggest eruption through low pyroclastic fountains (<1000 m high) with limited cooling during eruption and turbulent flow from an initial temperature of 1160°C. The large volume of 26 km³ of erupted basalt compared with only 16 km³ of the evolved P1 magmas, and the extremely high discharge rates inferred from model calculations are unusual for a basaltic eruption. It is suggested that the basaltic magma was erupted and emplaced in a fashion commonly only attributed to felsic magmas because it utilized the felsic P1 magma chamber and its ring-fissure conduits. Evolution of the entire P1 eruption was controlled by withdrawal dynamics involving magmas differing in viscosity by more than four orders of magnitude. The basaltic eruption phase was initially driven by buoyancy of the basaltic magma at chamber depth and continued degassing of felsic magma, but most of the large volume of basalt magma was driven out of the reservoir by subsidence of a c. 10 km diameter roof block, which followed a decrease in magma chamber pressure during low viscosity basaltic outflow.     

Quantifying the effect of rheology on lava-flow margins using fractal geometry

This study aims at quantifying the effect of rheology on plan-view shapes of lava flows using fractal geometry. Plan-view shapes of lava flows are important because they reflect the processes governing flow emplacement and may provide insight into lava-flow rheology and dynamics. In our earlier investigation (Bruno et al. 1992), we reported that flow margins of basalts are fractal, having a characteristic shape regardless of scale. We also found we could use fractal dimension (D, a parameter which quantifies flow-margin convolution) to distinguish between the two endmember types of basalts: a a (D: 1.05–1.09) and pahoehoe (D: 1.13–1.23). In this work, we confirm those earlier results for basalts based on a larger database and over a wider range of scale (0.125 m–2.4 km). Additionally, we analyze ten silicic flows (SiO₂: 52–74%) over a similar scale range (10 m–4.5 km). We note that silicic flows tend to exhibit scale-dependent, or non-fractal, behavior. We attribute this breakdown of fractal behavior at increased silica contents to the suppression of small-scale features in the flow margin, due to the higher viscosities and yield strengths of silicic flows. These results suggest we can use the fractal properties of flow margins as a remote-sensing tool to distinguish flow types. Our evaluation of the nonlinear aspects of flow dynamics indicates a tendency toward fractal behavior for basaltic lavas whose flow is controlled by internal fluid dynamic processes. For silicic flows, or basaltic flows whose flow is controlled by steep slopes, our evaluation indicates non-fractal behavior, consistent with our observations.     

Thick lava flows of Karisimbi Volcano, Rwanda: insights from SIR-C interferometric topography

 We use a digital elevation model (DEM) derived from interferometrically processed SIR-C radar data to estimate the thickness of massive trachyte lava flows on the east flank of Karisimbi Volcano, Rwanda. The flows are as long as 12 km and average 40–60 m (up to >140 m) in thickness. By calculating and subtracting a reference surface from the DEM, we derived a map of flow thickness, which we used to calculate the volume (up to 1 km³ for an individual flow, and 1.8 km³ for all the identified flows) and yield strength of several flows (23–124 kPa). Using the DEM we estimated apparent viscosity based on the spacing of large folds (1.2×10¹² to 5.5×10¹² Pa s for surface viscosity, and 7.5×10¹⁰ to 5.2×10¹¹ Pa s for interior viscosity, for a strain interval of 24 h). We use shaded-relief images of the DEM to map basic flow structures such as channels, shear zones, and surface folds, as well as flow boundaries. The flow thickness map also proves invaluable in mapping flows where flow boundaries are indistinct and poorly expressed in the radar backscatter and shaded-relief images. Received: 6 September 1997 / Accepted: 15 May 1998     

Two-dimensional kinematic and rheological modeling of the 1912 pyroclastic flow,Katmai, Alaska

Pyroclastic flow emplacement is strongly influenced by eruption column height. A surface along which kinetic energy is zero theoretically connects the loci of eruption column collapse with all coeval ignimbrite termini. This surface is reconstructed as a two-dimensional energy line for the 1912 Katmai pyroclastic flow in the Valley of Ten Thousand Smokes from mapped flow termini and the runup of the ignimbrite onto obstructions and through passes. Extrapolation of the energy line to the vicinity of the source vent at Novarupta suggests the eruption column which generated the ignimbrite eruption was approximately 425 m high. The 1912 pyroclastic flow travelled about 25 km downvalley. Empirical velocity data calculated from runup elevations and surveyed centrifugal superelevations indicate initial velocities near Novarupta were greater than 79–88 m s^–1. The flow progressively decelerated and was travelling only 2–8 m s^–1 when it crossed a moraine 16 km downvalley. The constant slope of the energy line away from Novarupta suggests the flow was systematically slowed by internal and basal friction. Using a simple physical model to calculate flow velocities and a constant kinetic friction coefficient (Heim coefficient) of 0.04 derived from the reconstructed energy line, the flow is estimated to have decelerated at an average rate of –0.16 m s^–2 and to have taken approximately 9.5 minutes to travel 25 km down the Valley of Ten Thousand Smokes. The shear strength of the flowing ignimbrite at the moraine was approximately 0.5 kPa, and its Bingham viscosity when it crossed the moraine was 3.5 × 10³ P. If the flow was Newtonian, its viscosity was 4.2 × 10³ P. Reynolds and Froude numbers at the moraine were only 41–62 and 0.84–1.04, respectively, indicating laminar, subcritical flow.     

Pahoehoe and aa in Hawaii: volumetric flow rate controls the lava structure

The historical records of Kilauea and Mauna Loa volcanoes reveal that the rough-surfaced variety of basalt lava called aa forms when lava flows at a high volumetric rate (>5–10 m³/s), and the smooth-surfaced variety called pahoehoe forms at a low volumetric rate (<5–10 m³/s). This relationship is well illustrated by the 1983–1990 and 1969–1974 eruptions of Kilauea and the recent eruptions of Mauna Loa. It is also illustrated by the eruptions that produced the remarkable paired flows of Mauna Loa, in which aa formed during an initial short period of high discharge rate (associated with high fountaining) and was followed by the eruption of pahoehoe over a sustained period at a low discharge rate (with little or no fountaining). The finest examples of paired lava flows are those of 1859 and 1880–1881. We attribute aa formation to rapid and concentrated flow in open channels. There, rapid heat loss causes an increase in viscosity to a threshold value (that varies depending on the actual flow velocity) at which, when surface crust is torn by differential flow, the underlying lava is unable to move sufficiently fast to heal the tear. We attribute pahoehoe formation to the flowage of lava at a low volumetric rate, commonly in tubes that minimize heat loss. Flow units of pahoehoe are small (usually <1 m thick), move slowly, develop a chilled skin, and become virtually static before the viscosity has risen, to the threshold value. We infer that the high-discharge-rate eruptions that generate aa flows result from the rapid emptying of major or subsidiary magma chambers. Rapid near-surface vesiculation of gas-rich magma leads to eruptions with high discharge rates, high lava fountains, and fast-moving channelized flows. We also infer that long periods of sustained flow at a low discharge rate, which favor pahoehoe, result from the development of a free and unimpeded pathway from the deep plumbing system of the volcano and the separation of gases from the magma before eruption. Achievement of this condition requires one or more episodes of rapid magma excursion through the rift zone to establish a stable magma pathway.     

Zenith angular distributions of atmospheric high-energy neutrinos

In this paper, we present a new calculation of the atmospheric neutrino flux in the energy range 10–10⁷ GeV, which reveals sizable differences in muon neutrino flux predictions obtained with known hadronic models. The calculation is based on the method of solving nuclear cascade equations in the atmosphere, which takes into account the nonscaling behavior of inclusive cross sections for particle production, the increase in the total inelastic hadron-nucleus cross sections, and the non-power-law character of the primary cosmic ray spectrum. The efficiency of the method was recently tested in atmospheric muon flux calculations. The results of neutrino spectrum calculations have been compared with Frejus, AMANDA-II, and IceCube measurement data.     

The reaction Mg(n,α)Ne at 14.1 and 14.7 MeV: cross sections and implications for meteorites

Mass spectrometric analyses of neutron-irradiated targets of natural magnesium yield cross sections of59 ± 14,160 ± 8, and11.0 ± 3.3mb for²⁰Ne,²¹Ne, and²²Ne, respectively, at 14.1 MeV and of94 ± 8,152 ± 12, and13.0 ± 2.0mb at 14.7 MeV. With the incorporation of these cross sections, calculations modeling cosmic-ray interactions in stony meteoroids of radii 20 and 26 cm predict that between the surface and center the²²Ne/²¹Ne ratio falls more than 10% while the²¹Ne production rate rises about 30%. The reaction²⁴Mg(n,α)²¹Ne predominantly controls these trends: the²²Ne/²¹Ne ratio due to magnesium decreases over 15% while that due to silicon remains constant with increasing depth. The calculations are compared with published neon measurements for the Keyes and St. Séverin meteorites.     

Toothpaste lava: Characteristics and origin of a lava structural type transitional between pahoehoe and aa

Toothpaste lava, an important basalt structural type which illustrates the transition from pahoehoe to aa, is particularly well displayed on the 1960 Kapoho lava of Kilauea Volcano. Its transitional features stem from a viscosity higher than that of pahoehoe and a rate of flow slower than that of aa. Viscosity can be quantified by the limited settling of olivine phenocrysts and rate of flow by field observations related to the low-angle slope on which the lava flowed. Much can be learned about the viscosity, rheologic condition, and flow velocity of lavas long after solidification by analyses of their structural characteristics, and it is possible to make at least a semiquantitative assessment of the numerical values of these parameters.     

Field observations,rheological testing and numerical modelling of a debris‐flow event

Debris flows generated from landslides are common processes and represent a severe hazard in mountain regions due to their high mobility and impact energy. We investigate the dynamics and the rheological properties of a 90 000 m³ debris‐flow event triggered by a rapid regressive landslide with high water content. Field evidence revealed a maximum flow depth of 7–8 m, with an estimated peak discharge of 350–400 m³ s^?1. Depositional evidence and grain‐size distribution of the debris pose the debris flow in an intermediate condition between the fluid‐mud and grain‐flow behaviour. The debris‐flow material has silt–clay content up to 15 per cent. The rheological behaviour of the finer matrix was directly assessed with the ball measuring system. The measurements, performed on two samples at 45–63 per cent in sediment concentration by volume, gave values of 3·5–577 Pa for the yield strength, and 0·6–27·9 Pa s for the viscosity. Based on field evidence, we have empirically estimated the yield strength and viscosity ranging between 4000 ± 200 Pa, and 108–134 Pa s, respectively. We used the Flo‐2D code to replicate the debris‐flow event. We applied the model with rheological properties estimated by means of direct measurements and back‐analyses. The results of these models show that the rheological behaviour of a debris‐flow mass containing coarse clasts can not be assessed solely on the contribution of the finer matrix and thus neglecting the effects of direct grain contacts. For debris flows composed by a significant number of coarse clasts a back‐analysis estimation of the rheological parameters is necessary to replicate satisfactorily the depositional extent of debris flows. In these cases, the back‐estimated coefficients do not adequately describe the rheological properties of the debris flow. Copyright © 2006 John Wiley & Sons, Ltd.     

Crustal thermal models for South Siberia

Summary The results of geothermal research carried out in South Siberia (West-Siberian and Siberian Platforms, Altai-Sayan folded area and the Baikalian arched uplift zone) were employed to calculate the Earth's crustal temperatures. It is shown that temperature is a function of the heat flow value. The maps of the surface heat flow and the temperatures at the Moho discontinuity are presented and characteristic temperature-depth profiles were compiled for specific tectonic areas. The West-Siberian and Siberian Platforms, as well as the Altai-Sayan folded region are characterized by rather low heat flow, 42–50mW m ^–2 (1.0–1.2 µcal/cm ² s), and low temperatures at the Moho boundary.400–500 °C. The thermal conditions in the area of Lake Baikal are extreme, the heat flow reaches here100–142 mW m ^–2 (2.4–3.4 µcal/cm ² s), and the temperature at the crust's bottom may exceed1000 °C. The Moho-surface according to our calculations is nonisothermal.     

Two-layer Ekman model for wind-induced shear flow near the sea surface

A simple two-layered model for steady wind-induced shear flow near the sea surface has been formulated. Basic assumptions of Ekman's theory are retained, including horizontal uniformity, infinite depth and constant (but differing) vertical eddy viscosities in the respective layers. Employing information coming from observational data, the parameters of the model (depth of the surface layer and the two eddy viscosities) are evaluated and optimized. The results thereby obtained favour the presence of a high shear layer, about 1 m thick, immediately below the sea surface. The eddy viscosity in that layer increases approximately linearly with wind speed but remains comparatively small, being generally less than 20 cm² s⁻¹ for wind speeds less than 20 m s⁻¹. In contrast, the eddy viscosity directly below the layer is two orders of magnitude larger and increases more steeply with wind speed.     

Rheology of the lower crust beneath the northern part of North China: Inferences from lower crustal xenoliths from Hannuoba basalts, Hebei Province, China

Lower crustal xenoliths brought up rapidly by basaltic magma onto the earth surface may provide di-rect information on the lower crust. The main purpose of this research is to gain an insight into the rheology of the lower crust through the detailed study of lower crustal xenoliths collected from the Hannuoba basalt, North China. The lower crustal xenoliths in this area consist mainly of two pyroxene granulite, garnet granulite, and light-colored granulite, with a few exception of felsic granulite. The equilibration temperature and pressure of these xenoliths are estimated by using geothermometers and geobarometers suitable for lower crustal xenoliths. The obtained results show that the equilibration temperature of these xenoliths is within the range of 785―900℃, and the equilibrium pressure is within the range of 0.8―1.2 GPa, corresponding to a depth range of 28―42 km. These results have been used to modify the previously constructed lower crust-upper mantle geotherm for the studied area. The dif-ferential stress during the deformation process of the lower crustal xenoliths is estimated by using recrystallized grain-size paleo-piezometer to be in the range of 14―20 MPa. Comparing the available steady state flow laws for lower crustal rocks, it is confirmed that the flow law proposed by Wilks et al. in 1990 is applicable to the lower crustal xenoliths studied in this paper. The strain rate of the lower crust estimated by using this flow law is within the range of 10-13―10-11 s-1, higher than the strain rate of the upper mantle estimated previously for the studied area (10-17―10-13 s-1); the equivalent viscosity is estimated to be within the range of 1017―1019Pa·s, lower than that of the upper mantle (1019―1021 Pa·s). The constructed rheological profiles of the lower crust indicate that the differential stress shows no significant linear relation with depth, while the strain rate increases with depth and equivalent vis-cosity decrease with depth. The results support the viewpoint of weak lower continental crust.     

Viscosity measurements of subliquidus magmas: Alkali olivine basalt from the Higashi-Matsuura district, Southwest Japan

We carried out viscosity measurements and sampling of a crystal suspension derived from alkali olivine basalt from the Matsuura district, SW Japan, at subliquidus temperatures from 1230 °C to 1140 °C under 1 atm with NNO oxygen buffered conditions. Viscosity increased from 31 to 1235 Pa s with a decrease in temperature from 1230 to 1140 °C. On cooling, olivine first appeared at 1210 °C, followed by plagioclase at 1170 °C. The crystal content of the sample attained 31 vol.% at 1140 °C (plagioclase 22%, olivine 9%). Non-Newtonian behaviors, including thixotropy and shear thinning, were pronounced in the presence of tabular plagioclase crystals. The cause of such behavior is discussed in relation to shear-induced changes in melt–crystal textures. Relative viscosities, η_r (= η_s / η_m, where η_s and η_m are the viscosities of the suspension and the melt, respectively), were obtained by calculating melt viscosities from the melt composition and temperature at 1 atm using the equation proposed by Giordano and Dingwell [Giordano, D., Dingwell, D.B., 2003. Non-Arrhenian multicomponent melt viscosity: a model. Earth and Planetary Science Letters, 208, 337–349.]. The obtained relative viscosities are generally consistent with the Einstein–Roscoe relation, which represents η_r for suspensions that contain equant and equigranular crystals, even though the crystal suspension analyzed in the present experiments contained tabular plagioclase and granular olivine of various grain sizes. This consistency is attributed to the fact that the effect of crystal shape was counterbalanced by the effect of the dispersion of crystal size. The applicability of the Einstein–Roscoe equation with respect to crystal shape is discussed on the basis of the present experimental results. Our experiments and those of Sato [Sato, H., 2005. Viscosity measurement of subliquidus magmas: 1707 basalt of Fuji volcano. Journal of Mineralogical and Petrological Sciences, 100, 133–142.] show that the relationship between relative viscosity and crystal fraction is consistent with the Einstein–Roscoe relationship for axial ratios that are smaller than the critical value of 4–6.5, but discrepancies occur for higher ratios.     

Ultrasonic measurements of Vp and Qp: relaxation spectrum of complex modulus on basalt melts

Ultrasonic compressional wave velocity V_p and quality factor Q_p have been measured in alkali basalt, olivine basalt and basic andesite melts in the frequency range of 3.4–22 MHz and in the temperature range of 1100–1400°C. Velocity and attenuation of the melts depend on frequency and temperature, showing that there are relaxation mechanisms in the melts. Complex moduli are calculated from the ultrasonic data. The results fit well a complex modulus of Arrhenius temperature dependence with log-normal Gaussian distribution in relaxation times of attenuation. The analysis yields average relaxation time, its activation energy, relaxed modulus, unrelaxed modulus and width of Gaussian distribution in relaxation times. Relaxed modulus is smaller (17.5 GPa) for basic andesite melt of high silica and high alumina contents than for the other two basalt melts (18.1–18.4 GPa). The most probable relaxation times decrease from ～ 3 × 10^?10 s for basic andesite to ～ 10^?11 s for alkali basalt at 1400°C. Activation energies of attenuation, ranging from 270 to 340 kJ mol^?1 in the three melts, are highest in basic andesite. Longitudinal viscosity values and their temperature dependences are also calculated from V_p and Q_p data. The volume viscosity values are estimated from the data using the shear viscosity values. Longitudinal, volume and shear viscosities and their activation energies are highest in the basic andesite melt of the most polymerized structure.     

Stability of the oceanic lithosphere with variable viscosity: an initial-value approach

We have studied the problem concerning the onset of convective instabilities below the oceanic lithosphere. A system of linear partial differential equations, in which the background temperature field is time-dependent, is integrated in time to monitor the evolution of incipient disturbances. Two types of rheologies have been examined. One depends strongly on temperature. The other involves a viscosity which is both temperature- and pressure-dependent. The results from this initial-value approach, in which the viscosity profiles migrate downward with time, reveal the importance of considering temperature- and pressure-dependent rheology in issues regarding the development of local instabilities in upper mantle convection. For temperature-dependent viscosity, viscosities of 0(10²⁰P) are required to produce instabilities with growth-rates of 0(.1/Ma). In contrast, these same growth rates can be attained for a temperature- and pressure-dependent viscosity profile with a mean value close to 0(10²⁰P) in the upper mantle, owing to the presence of a low viscosity zone, 0(10²⁰P), existing right below the lithosphere. Unlike the results of temperature-dependent viscosity, whose growth-rates increase with time, the amplification of disturbances in a fluid medium with temperature- and pressure-dependent rheology reaches a maximum at an early age, < 50 Ma, and decreases thereafter with time. This suggests the potential importance played by initial disturbances in the evolution of the oceanic lithosphere.     

Lateral variation in upper mantle viscosity: role of water

Differences in the viscosity of the earth's upper mantle beneath the western US (∼10¹⁸-10¹⁹ Pa s) and global average values based on glacial isostatic adjustment and other data (∼10²⁰-10²¹ Pa s) are generally ascribed to differences in temperature. We compile geochemical data on the water contents of western US lavas and mantle xenoliths, compare these data to water solubility in olivine, and calculate the corresponding effective viscosity of olivine, the major constituent of the upper mantle, using a power law creep rheological model. These data and calculations suggest that the low viscosities of the western US upper mantle reflect the combined effect of high water concentration and elevated temperature. The high water content of the western US upper mantle may reflect the long history of Farallon plate subduction, including flat slab subduction, which effectively advected water as far inland as the Colorado Plateau, hydrating and weakening the upper mantle.     

Mantle rheology: radial and lateral viscosity variations inferred from microphysical creep laws

Inferences on the rheology of the mantle based on theoretical and experimental rate equations for steady state creep are discussed and compared with results from geophysical models. The radial increase of viscosity by one to three orders of magnitude across the mantle, required by inversion of postglacial rebound and geodynamic data, is confirmed by microphysical models based on the estimation of continuous and discontinuous changes of creep parameters with depth. The upper mantle (viscosity 10²⁰–10²¹ Pa s) is likely to show non-Newtonian rheology (power-law creep) for average grain sizes larger than 0.1 mm as an order of magnitude. Given the variability of both grain size and stress conditions, local regions of linear rheology can be present. The rheology of transition zone and lower mantle (viscosity 10²²–10²⁴ Pa s) cannot be definitely resolved at present. Estimation of creep parameters leads to possible nonlinear or mixed rheology, if grain sizes are not lower than 0.1 mm and flow conditions can be approximated by a constant strain rate of about 10⁻¹⁵ s⁻¹. This conclusion can be modified by different flow conditions (e.g. a decrease in strain rate or constant viscous dissipation). Furthermore, experiments on fine-grained garnetites and perovskite analogues have shown that diffusion creep is predominant at laboratory conditions. However, the pressure dependence of creep in these phases is unknown, and therefore direct extrapolation to lower mantle conditions is necessarily speculative. Lateral variations of viscosity, largest in the upper and lowermost mantle (up to 2–4 orders of magnitude) are predicted by models based on lateral temperature anomalies derived from seismic tomographic models.