Phase relations of iron and iron–nickel alloys up to 300 GPa: Implications for composition and structure of the Earth's inner core

We have investigated the phase relations of iron and iron–nickel alloys with 18 to 50 wt.% Ni up to over 300 GPa using a laser-heated diamond-anvil cell. The synchrotron X-ray diffraction measurements show the wide stability of hcp-iron up to 301 GPa and 2000 K and 319 GPa and 300 K without phase transition to dhcp, orthorhombic, or bcc phases. On the other hand, the incorporation of nickel has a remarkable effect on expanding the stability field of fcc phase. The geometry of the temperature–composition phase diagram of iron–nickel alloys suggests that the hcp–fcc–liquid triple point is located at 10 to 20 wt.% Ni at the pressure of the inner core boundary. The fcc phase could crystallize depending on the nickel and silicon contents in the Earth's core, both of which are fcc stabilizer.     

Pre-eruptive and syn-eruptive conditions in the Black Butte, California dacite: Insight into crystallization kinetics in a silicic magma system

A series of experiments and petrographic analyses have been run to determine the pre-eruption phase equilibria and ascent dynamics of dacitic lavas composing Black Butte, a dome complex on the flank of Mount Shasta, California. Major and trace element analyses indicate that the Black Butte magma shared a common parent with contemporaneously erupted magmas at the Shasta summit. The Black Butte lava phenocryst phase assemblage (20 v.%) consists of amphibole, plagioclase (core An_77.5), and Fe–Ti oxides in a fine-grained (< 0.5 mm) groundmass of plagioclase, pyroxene, Fe–Ti oxides, amphibole, and cristobalite. The phenocryst assemblage and crystal compositions are reproduced experimentally between 890 °C and 910 °C, ≥ 300 MPa, X_H2O = 1, and oxygen fugacity = NNO + 1. This study has quantified the extent of three crystallization processes occurring in the Black Butte dacite that can be used to discern ascent processes. Magma ascent rate was quantified using the widths of amphibole breakdown rims in natural samples, using an experimental calibration of rim development in a similar magma at relevant conditions. The majority of rims are 34 ± 10 μm thick, suggesting a time-integrated magma ascent rate of 0.004–0.006 m/s among all four dome lobes. This is comparable to values for effusive samples from the 1980 Mount St. Helens eruption and slightly faster than those estimated at Montserrat. A gap between the compositions of plagioclase phenocryst cores and microlites suggests that while phenocryst growth was continuous throughout ascent, microlite formation did not occur until significantly into ascent. The duration of crystallization is estimated using the magma reservoir depth and ascent rate, as determined from phase equilibria and amphibole rim widths, respectively. Textural analysis of the natural plagioclase crystals yields maximum growth rates of plagioclase phenocryst rims and groundmass microlites of 8.7 × 10^− 8 and 2.5 × 10^− 8 mm/s, respectively. These rates are comparable to values determined from time-sequenced samples of dacite erupted effusively from Mount St. Helens during 1980 and indicate that syneruptive crystallization processes were important during the Black Butte eruptive cycle.     

Sedimentological constraints on hydrometeor-enhanced particle deposition: 1992 Eruptions of Crater Peak, Alaska

Water is a dominant component of volcanic clouds and has fundamental control on very fine particle deposition. Particle size characteristics of distal tephra-fall (100s km from source volcano) have a higher proportion of very fine particles compared to predictions based on single particle settling rates. In this study, sedimentological analyses of fallout from for the 18 August and 16–17 September 1992 eruptions of Crater Peak, Alaska, are combined with satellite observations, and cloud trajectory and microphysics modeling to investigate meteorological influences on particle sedimentation. Total grain size distributions of tephra fallout were reconstructed for both Crater Peak eruptions and indicate a predominance of fine particles < 125 μm. Polymodal analysis of the deposits has identified a particle subpopulation with mode ~ 15–18 μm involved in particle aggregation. Accounting for the magmatic water source only, calculated ice water content of the 3.7 hour old September 1992 Spurr cloud was ~ 4.5 × 10^− 2 g m^− 3 (based on an estimated cloud thickness of ~ 1000 m from trajectory modeling). Hydrometeor formation on particles in the volcanic cloud and subsequent sublimation may induce a cloud base instability that leads to rapid bulk (en masse) sedimentation of very fine particles through a mammatus-like mechanism.     

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.     

Ueber Herde elastischer Wellen in isotropen,homogenen Medien

Zusammenfassung 1) Es werden Multipollösungen der skalaren Wellengleichung ² f/t ² – c² div gradf=0 betrachtet. Einerseits kann man solche Lösungen direkt durch Kugelfunktionenn-ter Ordnung ausdrücken, anderseits aus der Einpollösungf=1/p F(t–p/c) durch Differentiation nachn Richtungen erhalten. Es wird der Zusammenhang zwischen den Ergebnissen der beiden Verfahren gezeigt. — 2) Für die Energiedichte und den Energiefluss durch Kugelflächen bei kleinen elastischen Verschiebungen werden Ausdrücke in Kugelkoordinaten angegeben. — 3) Für die Wellengleichung grad div –b ² rot rot werden rotationsfreie Multipollösungen angegeben und Ausdrücke für Energiedichte und Energiefluss hergeleitet. — 4) Das gleiche wird für divergenzfreie Multipollösungen durchgeführt. — 5) Es werden Multipole betrachtet, die weder rotationsfrei noch divergenzfrei sind. Als Spezialfälle werden Multipole mit zeitlich begrenzter und solche mit periodischer Erregung gezeigt, ferner Lösungen der Wellengleichung, die sowohl rotationsfrei wie divergenzfrei sind. — 6) Es wird gezeigt, wie man die elastischen Wellen, die im Sinne vonStokes von einem Herdgebiet endlicher Ausdehnung ausgehen, näherungsweise durch elastische Multipole darstellen kann. — 7) Es wird angedeutet, wie man durch Messung von Komponenten von oder u.s.w. in Punkten im Innern des Mediums die Erregung und Energie von elastischen Multipolen bestimmen kann. Ferner wird auf den Fall hingewiesen, wo ein rotationsfreier Einpol sich im Innern eines Halbraumes befindet und die Messungen an seiner Oberfläche ausgeführt werden.

---

Summary (On foci of elastic waves in isotropic homogeneous media) — 1) Multiplets as solutions of the scalar wave equation ² f/t ² – c² div gradf=0 are considered. Such solutions can be obtained either directly by aid of spherical harmonics of ordern, or by differentiating the single polef=1/p F(t–p/c) with respect ton directions. The relations between the results of those two procedures are shown. — 2) In the case of small elastic displacements , the density of energy and the flow of energy through spherical surfaces are expressed by spherical coordinates. — 3) Multiplets which satisfy the equation of motion =a ² grad div b ² curl curl and the equation curl = 0 are given, and expressions for the density and flow of energy are found. — 4) The same is done with multiplets satisfying the equation of motion and the equation div = 0. — 5) General multiplets which satisfy the equation of motion are treated. As special cases, multiplets with excitation of finite length and multiplets with periodic excitation are considered, furthermore solutions of the equation of motion and of the equations curl = 0 and div = 0 are given. — 6) It is shown how elastic waves whose origin is a region of finite extension in the sense given byStokes, can be approximated by elastic multiplets. — 7) Some indications are given on the problem of how to find the functions of excitation and the energy of an elastic multiplet by measuring components of or etc., at points in the interior of the medium. The same problem is considered in the case of the single elastic pole. = grad 1/p F (t–p/a), if the measurements are made at the surface of an elastic half space.

     

Experimental and theoretical study of stability of dense hydrous magnesium silicates in the deep upper mantle

Fluid-undersaturated experiments were conducted to determine the phase relations in the simplified peridotite system MgO-SiO₂-H₂O (MSH) at 11.0-14.5 GPa and 800-1400 °C. Stability relations of dense hydrous magnesium silicates (DHMSs) under fluid-undersaturated conditions were experimentally examined. From the fluid-absent experimental results, we retrieved thermodynamic data for clinohumite, phase A, phase E, and hydrous wadsleyite, consistent with the published data set for dry mantle minerals. With this new data set, we have calculated phase equilibria in the MSH system including dehydration reactions. The dehydration reactions calculated with lower water activities of 0.68-0.60 match the fluid-present experiments of this study above 11.0 GPa and 1000 °C, indicating that considerable amounts of silicate component were dissolved into the fluid phase. The calculated phase equilibria illustrate the stability of the post-antigorite phase A-bearing assemblages. In the cold subducting slab peridotite, phase A + enstatite assemblage survives into the transition zone, whereas phase A + forsterite + enstatite assemblage forms hydrous wadsleyite at a much shallower depth of about 360-km. The slab is subducted with no dehydration reactions occurring when entering the transition zone. The phase equilibria also show the high temperature stability of phase E. Phase E is stable up to 1200 °C at 13.5 GPa, a plausible condition in the mantle of relatively low temperature, i.e., beneath subduction zones. Phase E is a possible water reservoir in the mantle as well as wadsleyite and ringwoodite.     

Thermal equation of state of magnesite to 32 GPa and 2073 K

Pressure–volume–temperature relations have been measured to 32 GPa and 2073 K for natural magnesite (Mg_0.975Fe_0.015Mn_0.006Ca_0.004CO₃) using synchrotron X-ray diffraction with a multianvil apparatus at the SPring-8 facility. A least-squares fit of the room-temperature compression data to a third-order Birch–Murnaghan equation of state (EOS) yielded K₀ = 97.1 ± 0.5 GPa and K′ = 5.44 ± 0.07, with fixed V₀ = 279.55 ± 0.02 Å³. Further analysis of the high-temperature compression data yielded the temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.013 ± 0.001 GPa/K and zero-pressure thermal expansion α = a₀ + a₁T with a₀ = 4.03 (7) × 10⁻⁵ K⁻¹ and a₁ = 0.49 (10) × 10⁻⁸ K⁻². The Anderson–Grüneisen parameter is estimated to be δ_T = 3.3. The analysis of axial compressibility and thermal expansivity indicates that the c-axis is over three times more compressible (K_Tc = 47 ± 1 GPa) than the a-axis (K_Tc = 157 ± 1 GPa), whereas the thermal expansion of the c-axis (a₀ = 6.8 (2) × 10⁻⁵ K⁻¹ and a₁ = 2.2 (4) × 10⁻⁸ K⁻²) is greater than that of the a-axis (a₀ = 2.7 (4) × 10⁻⁵ K⁻¹ and a₁ = −0.2 (2) × 10⁻⁸ K⁻²). The present thermal EOS enables us to accurately calculate the density of magnesite to the deep mantle conditions. Decarbonation of a subducting oceanic crust containing 2 wt.% magnesite would result in a 0.6% density reduction at 30 GPa and 1273 K. Using the new EOS parameters we performed thermodynamic calculations for magnesite decarbonation reactions at pressures to 20 GPa. We also estimated stability of magnesite-bearing assemblages in the lower mantle.     

Geoelectric structure of the Ein Gedi sinkhole occurrence site at the Dead Sea shore in Israel

Hazardous sinkholes started to appear in alluvial fans and unconsolidated sediments along the western Dead Sea coast in 1990. Since then hundreds of sinkholes have appeared from north to south along the shoreline. The Electrical Resistivity Tomography (ERT) method was used to achieve a better understanding of the subsurface geoelectric structure at the sinkhole development sites, taking into account that electric parameters (such as resistivity or conductivity) are very sensitive to formation properties and their variations in time. Fifteen image lines were surveyed at the Ein Gedi area during a period of active sinkhole development (in 2001–2002) over an area of 300 × 550 m². Resistivity cross-sections and maps were constructed from 2-D linear surveys. The process of sinkhole formation in the surveyed area is located in a strip 50–70 m wide and 300–500 m long, extending approximately in a north–south direction. The sinkholes are arranged along a tortuous line within this strip. On resistivity maps and sections this U-shaped zone appears as an alternation of high resistivity anomalies of 350–1000 Ωm (at sinkhole group locations) with narrow background resistivity zones of 50–100 Ωm. The large size of resistivity anomalies (250 × 300 m²), which are considerably greater than those of the sinkholes, form one of the features of the sinkhole sites in the Ein Gedi area. The anomalies continue down to the water table or even deeper (maximum of 25–35 m depth). A low resistivity layer of 1–8 Ωm underlies them. The combined analysis of the image results and other geophysical data shows that high resistivity anomalies are associated with the decompaction of the soil mass at the sinkhole development sites and surrounding areas. Recent studies have shown that sinkholes in the Ein Gedi area are developing along the salt western edge located at a depth of 50 m. The subsurface high resistivity anomaly conforms to the sinkhole line (and salt boundary). They are presumably located above the great dissolution caverns at the salt edge. The heterogeneity of the resistivity structure within the high resistivity anomaly (seen in both lateral and vertical planes) confirms that a disintegration of internal formation structure takes place. Away from the sinkhole sites the subsurface resistivity distribution is homogeneous.     

Modern and Pleistocene boron isotope composition of the benthic foraminifer Cibicidoides wuellerstorfi

Here we present the first species-specific study of boron isotopes in the epibenthic foraminifer species Cibicidoides wuellerstorfi. Coretop samples from a water depth profile from 1000 to 4500 m on the northern flank of the Walvis Ridge are 4.4‰ lower than the values expected, based on calculations of the δ¹¹B_borate of ambient seawater. Similar values for this foraminifer species are presented from ODP site 668B at the Sierra Leone Rise, in the equatorial Atlantic. The consistency between data of the same species suggests the offsets are primary, rather than diagenetic. Glacial C. wuellerstorfi from ODP 668B and Walvis Ridge have boron isotope compositions only slightly different to interglacial samples, that is no larger than + 0.10 pH units, or + 23 µmol kg^− 1 in [CO₃²⁻] above the reconstructed glacial lysocline, and − 0.07 pH units, or − 14 µmol kg^− 1 in [CO₃²⁻] below. We use these results to suggest that glacial deep water pH in the Atlantic was similar to interglacial pH. The new data resolve the inconsistency between the previously reported high bottom water pH and the lack of significant carbonate preservation of the glacial deep ocean.     

Determination of phase transition pressures of ZnTe under quasihydrostatic conditions

Pressure behavior of ZnTe at room temperature was studied using an X-ray energy dispersive method on a DIA type cubic anvil apparatus (SAM-85) at NSLS-X17B1. By using powdered polyethylene, the sample and NaCl for a pressure scale were held under quasihydrostatic conditions, which were confirmed by X-ray diffraction method. Two high-pressure phase transitions were confirmed using X-ray powder diffraction simultaneously with electrical resistance measurements. The phase transition pressures under quasihydrostatic conditions were determined to be 9.6 GPa, at which the resistance increased, and 12.0 GPa, which was the midpoint of a large resistance decrease. Errors in the pressure determinations were estimated to be less than 0.2 GPa. These pressure values may depend on grain size and anisotropic stress effects on the calibrant. From X-ray observation of ZnTe, the bulk modulus of the zinc blende structure was calculated to beK ₀=51(3) GPa andK ₀ =3.6(0.8), and the first transition at 9.6 GPa was found to have about 9% volume change. It was consistent with an anomaly in the pressure generating curves.     

High-pressure phases of CaCO3: Crystal structure prediction and experiment

Post-aragonite phase of CaCO₃, experimentally known to be stable above 40 GPa [S. Ono, T. Kikegawa, Y. Ohishi, J. Tsuchiya, Post-aragonite phase transformation in CaCO₃ at 40 GPa, Am. Mineral. 90 (2005) 667–671], is believed to be a major carbon-containing mineral in the Earth's mantle. Crystal structure of this mineral phase could not be solved using experimental data or traditional theoretical simulation methods and remained a controversial issue. Using a combination of advanced ab initio simulation techniques and high-pressure experiment, we have been able to determine the crystal structure of CaCO₃ post-aragonite. Here, we performed simulations with the USPEX code [C.W. Glass, A.R. Oganov, and N. Hansen, (in preparation). USPEX: a universal structure prediction program], which is based on an evolutionary algorithm using ab initio free energy as the fitness function. This novel methodology for crystal structure prediction, which uses only the chemical composition as input, is described in detail. For CaCO₃, we identify a number of energetically competitive structures, the most stable of which closely matches the experimental powder diffraction pattern and, in agreement with experiment, becomes more stable than aragonite above 42 GPa. This structure belongs to a new structure type, which is also adopted by the high-pressure post-aragonite phases of SrCO₃ and BaCO₃. It has 2 formula units in the orthorhombic unit cell (space group Pmmn) and contains triangular CO₃²⁻ ions and Ca²⁺ ions in the 12-fold coordination. Above 137 GPa, a pyroxene-type structure (space group C222₁) with chains of CO₄⁴⁻ tetrahedra becomes more stable than post-aragonite. For MgCO₃, this structure becomes more stable than magnesite above 106 GPa and is a good candidate structure for MgCO₃ post-magnesite.     

Flux and dispersion of gases from the “Drachenschlund” hydrothermal vent at 8°18' S, 13°30' W on the Mid-Atlantic Ridge

Hydrothermal emission of mantle helium appears to be directly related to magma production rate, but other processes can generate methane and hydrogen on mid-ocean ridges. In an on-going effort to characterize these processes in the South Atlantic, the flux and distribution of these gases were investigated in the vicinity of a powerful black smoker recently discovered at 8°17.9' S, 13°30.4' W. The vent lies on the shoulder of an oblique offset in the Mid-Atlantic Ridge and discharges high concentrations of methane and hydrogen. Measurements during expeditions in 2004 and 2006 show that the ratio of CH₄ to ³He in the neutrally buoyant plume is quite high, 4 × 10⁸. The CTD stations were accompanied by velocity measurements with lowered acoustic Doppler current profilers (LADCP), and from these data we estimate the methane transport to have been 0.5 mol s^− 1 in a WSW-trending plume that seems to develop during the ebb tidal phase. This transport is an order of magnitude greater than the source of CH₄ calculated from its concentration in the vent fluid and the rise height of the plume. From this range of methane fluxes, the source of ³He is estimated to be between 0.14 and 1.2 nmol s^− 1. In either case, the ³He source is significantly lower than expected from the spreading rate of the Mid-Atlantic Ridge. From the inventory of methane in the rift valley adjacent to the vent, it appears that the average specific rate of oxidation is 2.6 to 23 yr^− 1, corresponding to a turnover time between 140 and 16 days. Vertical profiles of methane in the surrounding region often exhibited Gaussian-like distributions, and the variances appear to increase with distance from the vent. Using a Gaussian plume model, we obtained a range of vertical eddy diffusivities between 0.009 and 0.08 m²m² s^− 1. These high values may be due to tidally driven internal waves across the promontory on which the vent is located.     

Melting of enstatite from 13 to 18 GPa under hydrous conditions

Experiments on MgSiO₃ enstatite were conducted in the pressure range from 13 to 18 GPa under hydrous conditions in order to clarify the effect of water on the melting phase relations of enstatite at pressures corresponding to the Earth’s mantle transition zone. In some previous experiments [Geol. Soc. Am. Bull. 79 (1968) 1685; Phys. Earth Planet. Inter. 85 (1994) 237], incongruent melting behavior to form Mg₂SiO₄ forsterite and SiO₂ enriched liquid up to 5 GPa was observed, and congruent melting behavior at pressures up to 12 GPa was observed. Under hydrous conditions, we found that the melting reaction changes from congruent to incongruent at around 13.5 GPa. Liquid formed above 13.5 GPa is enriched in MgO component relative to MgSiO₃ because it coexists with stishovite (SiO₂). Moreover, the solidus temperature decreases drastically at around 13.5 GPa, in unison with the change in the melting reaction. The solidus temperature is about 1400 °C at 13 GPa, but approximately 900 °C at 15 GPa. Our results show that the liquidus phase changes from clinoenstatite to stishovite with increasing pressure and water content above 13.5 GPa. MgSiO₃ enstatite is one of the major constituent minerals in the Earth’s mantle, and it is expected that MgO-enriched liquid will be generated in the transition zone if water is present.     

Hydrothermal activity on the southern Mid-Atlantic Ridge: Tectonically- and volcanically-controlled venting at 4–5°S

We report results from an investigation of the geologic processes controlling hydrothermal activity along the previously-unstudied southern Mid-Atlantic Ridge (3–7°S). Our study employed the NOC (UK) deep-tow sidescan sonar instrument, TOBI, in concert with the WHOI (USA) autonomous underwater vehicle, ABE, to collect information concerning hydrothermal plume distributions in the water column co-registered with geologic investigations of the underlying seafloor. Two areas of high-temperature hydrothermal venting were identified. The first was situated in a non-transform discontinuity (NTD) between two adjacent second-order ridge-segments near 4°02′S, distant from any neovolcanic activity. This geologic setting is very similar to that of the ultramafic-hosted and tectonically-controlled Rainbow vent-site on the northern Mid-Atlantic Ridge. The second site was located at 4°48′S at the axial-summit centre of a second-order ridge-segment. There, high-temperature venting is hosted in an  18 km² area of young lava flows which in some cases are observed to have flowed over and engulfed pre-existing chemosynthetic vent-fauna. In both appearance and extent, these lava flows are directly reminiscent of those emplaced in Winter 2005−06 at the East Pacific Rise, 9°50′N and reference to global seismic catalogues reveals that a swarm of large (M 4.6−5.6) seismic events was centred on the 5°S segment over a  24 h period in late June 2002, perhaps indicating the precise timing of this volcanic eruptive episode. Temperature measurements at one of the vents found directly adjacent to the fresh lava flows at 5°S MAR (Turtle Pits) have subsequently revealed vent-fluids that are actively phase separating under conditions very close to the Critical Point for seawater, at  3000 m depth and 407 °C: the hottest vent-fluids yet reported from anywhere along the global ridge crest.     

Eocene melting of subducting continental crust and early uplifting of central Tibet: Evidence from central-western Qiangtang high-K calc-alkaline andesites, dacites and rhyolites

Changes in oceanic O–Sr isotopic compositions and global cooling beginning in the Eocene are considered to have been caused by the uplift of the Tibetan Plateau. The specific timing and uplift mechanism, however, have long been subjects of debate. We investigated the Duogecuoren lavas of the central-western Qiangtang Block, which form the largest outcrops among Cenozoic lavas in northern-central Tibet and have widely been considered as shoshonitic. Our study demonstrates, however, that most of these lavas are high-K calc-alkaline andesites, dacites and rhyolites. Moreover, they are characterized by high Sr (367–2472 ppm) and Al₂O₃ (14.55–16.86 wt.%) and low Y (3.05–16.9 ppm) and Yb (0.31–1.48 ppm) contents and high La/Yb (27–100) and Sr/Y (48–240) ratios, similar to adakitic rocks derived by partial melting of an eclogitic source. They can be further classified as either peraluminous and metaluminous subtypes. The peraluminous rocks have relatively high SiO₂ (> 66 wt.%) contents, and low MgO (< 1.0 wt.%), Cr (4.94–23.3 ppm) and Ni (2.33–17.0 ppm) contents and Mg^# (20–50) values, while the metaluminous rocks exhibit relatively low SiO₂ (55–69 wt.%) contents, and high MgO (1.41–6.34), Cr (25.7–383 ppm), Ni (14.13–183 ppm) and Mg^# (46–69) values, similar to magnesian andesites. ⁴⁰Ar/³⁹Ar and SHRIMP zircon U–Pb dating reveal that both peraluminous and metaluminous adakitic rocks erupted in the Eocene (46–38 Ma). Paleocene–Early Miocene thrust faults and associated syn-contractional basin deposits in the Qiangtang Block suggest that this region was undergoing crustal shortening within a continent during the Eocene. The low ε_Nd (− 2.81 to − 6.91) and high ⁸⁷Sr/⁸⁶Sr (0.7057–0.7097), Th (11.2–32.3 ppm) and Th/La (0.23–0.88) values in the Duogecuoren adakitic rocks further indicate that they were not derived by partial melting of subducted oceanic crust. Taking into account tectonic and geophysical data and the compositions of xenoliths in Cenozoic lava in northern-central Tibet, we suggest that the peraluminous adakitic rocks were most probably derived by partial melting of subducted sediment-dominated continent of the Songpan-Ganzi Block along the Jinsha suture to the north at a relatively shallow position (the hornblende + garnet stability field), but the metaluminous adakitic rocks likely originated from the interaction between peraluminous adakitic melts generated at greater depths (the garnet + rutile stability field) and mantle. Because the Duogecuoren adakitic rocks must have originated from a garnet-bearing (namely, eclogite facies) source, Eocene continental subduction along the Jinsha suture caused the thickening of the Qiangtang crust. Given that crustal thickening generally equates with elevation, the uplift of the Central Tibetan Plateau probably began as early as 45–38 Ma, which provides important evidence for tectonically driven models of oceanic O–Sr isotope evolution during global cooling and Asian continental aridification beginning in the Eocene.     

Simulating the dispersal of tephra from the 1991 Pinatubo eruption: Implications for the formation of widespread ash layers

PUFF and HAZMAP, two tephra dispersal models developed for volcanic hazard mitigation, are used to simulate the climatic 1991 eruption of Mt. Pinatubo. PUFF simulations indicate that the majority of ash was advected away from the source at the level of the tropopause (~ 17 km). Several eruptive pulses injected ash and SO₂ gas to higher altitudes (~ 25 km), but these pulses represent only a small fraction (~ 1%) of the total erupted material released during the simulation. Comparison with TOMS images of the SO₂ cloud after 71 and 93 h indicate that the SO₂ gas originated at an altitude of ~ 25 km near the source and descended to an altitude of ~ 22 km as the cloud moved across the Indian Ocean. HAZMAP simulations indicate that the Pinatubo tephra fall deposit in the South China Sea was formed by an eruption cloud with the majority of the ash concentrated at a height of 16–18 km. Results of this study demonstrate that the largest concentration of distal ash was transported at a level significantly below the maximum eruption column height (~ 40 km) and at a level below the calculated height of neutral buoyancy (~ 25 km). Simulations showed that distal ash transport was dominated by atmospheric circulation patterns near the regional tropopause. In contrast, the movement of the SO₂ cloud occurred at higher levels, along slightly different trajectories, and may have resulted from gas/particle segregations that took place during intrusion of the Pinatubo umbrella cloud as it moved away from source.     

Stability relation of (Mg,Fe)SiO3 garnets, major constituents in the Earth's interior

High-pressure and high temperature experiments at 20 GPa on (Mg,Fe)SiO₃ have revealed stability fields of two types of aluminium-free ferromagnesian garnets; non-cubic garnet and cubic garnet (majorite). Majorite garnet is stable only within a limited compositional variation, 0.2 < Fe/(Mg + Fe)< 0.4, and in the narrow temperature interval of 200°C around 2000°C, while the stability of non-cubic garnet with more iron-deficient compositions persists up to higher temperatures. These two garnets show fractional melting into iron-deficient garnet and iron-rich liquid, and the crystallization field of cubic garnet extends over Fe/(Mg + Fe)= 0.5. The assemblage silicate spinel and stishovite is a low-temperature phase, which also occurs in the iron-rich portion of the MgSiO₃—FeSiO₃ system. The sequence as given by the Fe/(Mg + Fe) value for the coexisting phases with the two garnets at 2000°C and 20 GPa is: silicate modified spinel aluminium-free garnets silicate spinel.Natural majorite in shock-metamorphosed chondrites is clarified to be produced at pressures above 20 GPa and temperatures around 2000°C. Similar shock events may cause the occurrence of non-cubic garnet in iron-deficient meteorites. Non-cubic garnet could be a stable phase in the Earth's mantle if a sufficiently low concentration of aluminium is present in the layer corresponding to the stable pressure range of non-cubic garnet. The chemical differentiation by melting in the deep mantle is also discussed on the basis of the present experimental results and the observed coexistence of majorite garnet with magnesiowüstite in chondrites.     

Optically stimulated luminescence dating of young (<500 years old) sediments: Testing estimates of burial dose

The reliability of single-grain optically stimulated luminescence (OSL) dating of young (<500 years old) samples is assessed by examining the results of 166 single-grain OSL analyses of Australian alluvial, lacustrine, estuarine, coastal and marine sediment samples. All samples have been analysed using the modified SAR protocol of Olley et al. [2004b. Optical dating of Holocene sediments from a variety of geomorphic settings using single grains of quartz. Geomorphology 60, 337–358], with burial doses (D_b) calculated by fitting a single Gaussian curve to the peak of a multi-Gaussian summed probability distribution (PDF_Gaussian). The near-equivalence of this approach to currently accepted methods of palaeodose estimation – namely the minimum age and central age models of Galbraith et al. [1999. Optical dating of single and multiple grains of quartz from Jinmium rock shelter, northern Australia: part I, experimental design and statistical models. Archaeometry 41, 339–364] – is demonstrated. Considered as a whole, OSL ages calculated using the PDF_Gaussian D_b determination method are consistent with other chronometric, geomorphological, botanical and historical indicators of depositional age, while four alternative (previously published) methods are generally less reliable when applied to the samples analysed here. The single-grain data from 12 samples of known age (reported by Olley et al. [2004b. Optical dating of Holocene sediments from a variety of geomorphic settings using single grains of quartz. Geomorphology 60, 337–358]) are re-analysed using the PDF_Gaussian approach, with all but one yielding a D_b consistent with the independently known age. Examples of luminescence data (dose recovery tests, preheat plateaux, continuous-wave and linear-modulated OSL, growth curves, single-grain dose distributions) are provided that together indicate the suitability of Australian quartz for dating young sediments.