The November 2002 eruption of Piton de la Fournaise,Réunion: tracking the pre-eruptive thermal evolution of magma using melt inclusions

The November 2002 eruption of Piton de la Fournaise in the Indian Ocean was typical of the activity of the volcano from 1999 to 2006 in terms of duration and volume of magma ejected. The first magma erupted was a basaltic liquid with a small proportion of olivine phenocrysts (Fo₈₁) that contain small numbers of melt inclusions. In subsequent flows, olivine crystals were more abundant and richer in Mg (Fo_83–84). These crystals contain numerous melt and fluid inclusions, healed fractures, and dislocation features such as kink bands. The major element composition of melt inclusions in this later olivine (Fo_83–84) is out of equilibrium with that of its host as a result of extensive post-entrapment crystallization and Fe²⁺ loss by diffusion during cooling. Melt inclusions in Fo₈₁ olivine are also chemically out of equilibrium with their hosts but to a lesser degree. Using olivine–melt geothermometry, we determined that melt inclusions in Fo₈₁ olivine were trapped at lower temperature (1,182 ± 1°C) than inclusions in Fo_83–84 olivine (1,199–1,227°C). This methodology was also used to estimate eruption temperatures. The November 2002 melt inclusion compositions suggest that they were at temperatures between 1,070°C and 1,133°C immediately before eruption and quenching. This relatively wide temperature range may reflect the fact that most of the melt inclusions were from olivine in lava samples and therefore likely underwent minor but variable amounts of post-eruptive crystallization and Fe²⁺ loss by diffusion due to their relatively slow cooling on the surface. In contrast, melt inclusions in tephra samples from past major eruptions yielded a narrower range of higher eruption temperatures (1,163–1,181°C). The melt inclusion data presented here and in earlier publications are consistent with a model of magma recharge from depth during major eruptions, followed by storage, cooling, and crystallization at shallow levels prior to expulsion during events similar in magnitude to the relatively small November 2002 eruption.     

Cooling and crystallization of lava in open channels, and the transition of Pāhoehoe Lava to 'A'ā

 Samples collected from a lava channel active at Kīlauea Volcano during May 1997 are used to constrain rates of lava cooling and crystallization during early stages of flow. Lava erupted at near-liquidus temperatures (∼1150  °C) cooled and crystallized rapidly in upper parts of the channel. Glass geothermometry indicates cooling by 12–14  °C over the first 2 km of transport. At flow velocities of 1–2 m/s, this translates to cooling rates of 22–50  °C/h. Cooling rates this high can be explained by radiative cooling of a well-stirred flow, consistent with observations of non-steady flow in proximal regions of the channel. Crystallization of plagioclase and pyroxene microlites occurred in response to cooling, with crystallization rates of 20–50% per hour. Crystallization proceeded primarily by nucleation of new crystals, and nucleation rates of ∼10⁴/cm³s are similar to those measured in the 1984 open channel flow from Mauna Loa Volcano. There is no evidence for the large nucleation delays commonly assumed for plagioclase crystallization in basaltic melts, possibly a reflection of enhanced nucleation due to stirring of the flow. The transition of the flow surface morphology from pāhoehoe to 'a'ā occurred at a distance of 1.9 km from the vent. At this point, the flow was thermally stratified, with an interior temperature of ∼1137  °C and crystallinity of ∼15%, and a flow surface temperature of ∼1100  °C and crystallinity of ∼45%. 'A'ā formation initiated along channel margins, where crust was continuously disrupted, and involved tearing and clotting of the flow surface. Both observations suggest that the transition involved crossing of a rheological threshold. We suggest this threshold to be the development of a lava yield strength sufficient to prevent viscous flow of lava at the channel margin. We use this concept to propose that 'a'ā formation in open channels requires both sufficiently high strain rates for continued disruption of surface crusts and sufficient groundmass crystallinity to generate a yield strength equivalent to the imposed stress. In Hawai'i, where lava is typically microlite poor on eruption, these combined requirements help to explain two common observations on 'a'ā formation: (a) 'a'ā flow fields are generated when effusion rates are high (thus promoting crustal disruption); and (b) under most eruption conditions, lava issues from the vent as pāhoehoe and changes to 'a'ā only after flowing some distance, thus permitting sufficient crystallization. Received: 3 September 1998 / Accepted: 12 April 1999     

Thermal structure and heat loss at the summit crater of an active lava dome

Forward-Looking Infrared (FLIR) nighttime thermal images were used to extract the thermal and morphological properties for the surface of a blocky-to-rubbley lava mass active within the summit crater of the Caliente vent at Santiaguito lava dome (Guatemala). Thermally the crater was characterized by three concentric regions: a hot outer annulus of loose fine material at 150–400°C, an inner cold annulus of blocky lava at 40–80°C, and a warm central core at 100–200°C comprising younger, hotter lava. Intermittent explosions resulted in thermal renewal of some surfaces, mostly across the outer annulus where loose, fine, fill material was ejected to expose hotter, underlying, material. Surface heat flux densities (radiative + free convection) were dominated by losses from the outer annulus (0.3–1.5 × 10⁴  s⁻¹m⁻²), followed by the hot central core (0.1–0.4 × 10⁴ J s⁻¹m⁻²) and cold annulus (0.04–0.1 × 10⁴ J s⁻¹m⁻²). Overall surface power output was also dominated by the outer annulus region (31–176 MJ s⁻¹), but the cold annulus contributed equal power (2.41–7.07 MJ s⁻¹) as the hot central core (2.68–6.92 MJ s⁻¹) due to its greater area. Cooled surfaces (i.e. the upper thermal boundary layer separating surface temperatures from underlying material at magmatic temperatures) across the central core and cold annulus had estimated thicknesses, based on simple conductive model, of 0.3–2.2 and 1.5–4.3 m. The stability of the thermal structure through time and between explosions indicates that it is linked to a deeper structural control likely comprising a central massive plug, feeding lava flow from the SW rim of the crater, surrounded by an arcuate, marginal fracture zone through which heat and mass can preferentially flow.     

Estimation of paleo-firing temperatures using luminescence signals for the volcanic lava baked layer in Datong,China

Eight paleo-fired samples from the baked layer in different depths under the lava and one unfired sample were collected from Datong, China. Fine quartz grains (4–11 μm) from samples were used for probing into relationship between luminescence signals and paleo-firing temperatures. Findings from the re-firing experiments indicated that using thermoluminescence (TL) and optically stimulated luminescence (OSL) sensitivity changes could estimate the paleo-firing temperatures of samples: (1) 110 °C TL sensitivity change rate against the re-firing temperature can tell whether the sample has been fired to temperatures above 500 °C or not; (2) 150 °C TL sensitivity against the re-firing temperature can indicate whether the sample has been fired to temperatures above 900 °C or not; (3) the more specific paleo-firing temperatures can be estimated by comparing the ratio of OSL and 150 °C TL sensitivities against re-firing temperatures. Results showed that the paleo-firing temperatures of the eight lava-baked samples decreased exponentially with the distance from the lava. Based on the estimated temperature profile, the temperature of the lava was estimated to be about 1100 °C.     

Morphology and composition of spinel in Pu’u ’O’o lava (1996–1998), Kilauea volcano, Hawaii

The morphology and composition of spinel in rapidly quenched Pu’u ’O’o vent and lava tube samples are described. These samples contain glass, olivine phenocrysts (3–5 vol.%) and microphenocrysts of spinel (0.05 vol.%). The spinel surrounded by glass occurs as idiomorphic octahedra 5–50 μm in diameter and as chains of octahedra that are oriented with respect to each other. Spinel enclosed by olivine phenocrysts is sometimes rounded and does not generally form chains. The temperature before quenching was calculated from the MgO content of the glass and ranges from 1150°C to 1180°C. The oxygen fugacity before quenching was calculated by two independent methods and the log f_O2 ranged from −9.2 to −9.9 (delta QFM=−1). The spinel in the Pu’u ’O’o samples has a narrow range in composition with Cr/(Cr+Al)=0.61 to 0.73 and Fe²⁺/(Fe²⁺+Mg)=0.46 to 0.56. The lower the calculated temperature for the samples, the higher the average Fe²⁺/(Fe²⁺+Mg), Fe³⁺ and Ti in the spinel. Most zoned spinel crystals decrease in Cr/(Cr+Al) from core to rim and, in the chains, the Cr/(Cr+Al) is greater in the core of larger crystals than in the core of smaller crystals. The occurrence of chains and hopper crystals and the presence of Cr/(Cr+Al) zoning from core to rim of the spinel suggest diffusion-controlled growth of the crystals. Some of the spinel crystals may have grown rapidly under the turbulent conditions of the summit reservoir and in the flowing lava, and the crystals may have remained in suspension for a considerable period. The rapid growth may have caused very local (μm) gradients of Cr in the melt ahead of the spinel crystal faces. The crystals seem to have retained the Cr/(Cr+Al) ratio that developed during the original growth of the crystal, but the Fe²⁺/(Fe²⁺+Mg) ratio may have equilibrated fairly rapidly with the changing melt composition due to olivine crystallization. Six of the samples were collected on the same day at various locations along a 10-km lava tube and the calculated pre-collection temperatures of the samples show a 5°C drop with distance from the vent. The average Fe²⁺/(Fe²⁺+Mg) of the spinel in these samples shows a weak positive correlation with decreasing MgO in the glass of these samples. The range in Cr₂O₃ (0.041–0.045 wt.%) of the glass for these six samples is too small to distinguish a consistent change along the lava tube. The spinel in the Pu’u ’O’o samples shows a zoning trend in a Cr–Al–Fe³⁺ diagram almost directly away from the Cr apex. This compares with a zoning trend in rapidly quenched MORB samples away from Cr coupled with decreasing Fe³⁺. The trend away from Cr displayed by spinel in rapidly quenched samples is in marked contrast to the trend of increasing Fe³⁺ shown by spinel in slowly cooled lava.     

Precise revision of the garnet-muscovite geothermometer

The garnet-muscovite geothermometer was refined through empirical calibration by using natural rocks metamorphosed under the physical conditions of 238—1306 MPa and 490—700℃. Input temperatures and pressures were determined through simultaneously applying the garnet-biotite geothermometer and the garnet-biotite-plagioclase-quartz barometer, assuming that all FeO in muscovite and garnet be ferrous. Garnet was treated as the asymmetric quaternary solid solution, and muscovite as the symmetric binary solid solution. Input muscovite compositions include Fe atoms between 0.03—0.19 and Mg atoms between 0.04—0.16 on the basis of 11 oxygen atoms, and input garnet compositions include spessartine fractions between 0.01—0.289, grossular fractions between 0.028—0.273, and the Fe/Mg ratio between 3.387-18.986. The resulting garnet-muscovite geothermometer reproduces temperatures within (50℃ compared with the garnet-biotite thermometer. Total random error of ±37℃ of the new thermometer may stem from the pressure uncertainty of ±200 MPa, and uncertainties of ±5% of Fe and Mg components in muscovite, and ±5% of Fe, Mg, Mn and Ca components in garnet, altogether. When there exist 10%, 20%, 30%, 40% and 50% Fe3+ in muscovite, respectively, the computed garnet-muscovite temperatures will be 1—6℃, 2—12℃, 3—16℃, 5—24℃ and 7—29℃, respectively, lower than those obtained when assuming that all FeO be ferrous. The new garnet-muscovite geothermometer can efficiently reflect temperature change of typical prograde sequences and contact aureole rocks, and may be applied to low- to high-grade and low- to high-pressure metamorphic rocks.     

Experimental constraints on the crystallization of natrocarbonatitic lava flows

Natrocarbonatitic magmas are characterized by their extremely low viscosities and fast elemental diffusion, and as a consequence of this, their chemistry and crystallinity can change significantly during residence in shallow reservoirs or even due to cooling during lava flow emplacement. Here, we present the results of a series of crystallization experiments conducted at 1-atm confining pressure and in a temperature range between 630°C and 300°C. The experiments were set up to characterize the chemistry and growth processes of the phenocryst phases present in natrocarbonatites. The results are applicable to (1) processes occurring during residence in shallow magma reservoirs and/or (2) during lava flow emplacement. We show that during crystallization of natrocarbonatites at atmospheric pressure, gregoryite is the first mineral to crystallize at 630°C, followed by nyerereite at 595°C. Crystal size distributions of the gregoryites show that the crystals grow rapidly by textural coarsening (i.e., Ostwald ripening). As the crystallization is a continuous process at this pressure, the composition of the residual melt changes in response to the crystallization. However, the experiments also show that individual crystals completely reequilibrate with the changes in melt composition in as little time as <11 min. We therefore conclude that crystallization and diffusion are extremely fast processes in the natrocarbonatitic system and that the measured chemical variations in phenocrysts from Oldoinyo Lengai can be explained by different cooling histories. Finally, we model the rheological control on the emplacement of highly crystallized natrocarbonatitic lavas at Oldoinyo Lengai.     

Mineralogy and petrology of sample 67075 and the origin of lunar anorthosites

Plagioclase in cataclastic anorthosite 67075 occurs as angular matrix grains and as recrystallized clasts of micro-anorthosite. Olivines are Fe-rich and fall into two compositional groupings. Large grains of pyroxene show exceptionally well-developed exsolution lamellae analogous to those observed in pyroxenes from layered complexes. The low-Ca component in both pigeonites and augites shows varying degrees of inversion to orthopyroxene. The lattices of host and lamellae may deviate slightly (up to 2°) from the ideal orientation. Very slow cooling from magmatic temperatures is required to produce the coarse exsolution textures and inversion features. Augite macrocrystals are distinctly subcalcic indicating crystallization at temperatures around1100 ± 50°C while host-lamellae pairs and small grains in lithic clasts and matrix indicate reequilibration on a micron scale to temperatures less than 800°C. Pyroxene compositions tend to cluster into two groups both of which are among the most Fe-rich reported for highland pyroxenes. Ti and Al contents of pyroxenes are very low and Ti, Cr, and Mn follow well-established magmatic differentiation trends. The high Cr content may reflect low?O₂ conditions and/or early crystallization of olivine and plagioclase.The⁸⁷Sr/⁸⁶Sr ratios in lunar anorthosites are the lowest reported for any lunar rock. It is likely that anorthosites formed as cumulates during the major differentiation episode which occurred prior to～4.3AE. Recrystallization features are common and³⁹Ar/⁴⁰Ar ages cluster around 4.0 AE. This may be the result of the intense bombardment prior to 4.0 AE which caused repeated cycles of in-situ fracturing and granulation followed by recrystallization. The low siderophile element content and the inferred slow cooling indicate a plutonic source region (10km) not frequently plumbed by impact events. The Fe-rich silicates indicate crystallization from a melt at an advanced stage of fractionation. However, the low REE abundances are not consistent with late-stage crystallization. Plagioclase apparently crystallized relatively early and was concentrated by flotation and/or convection currents while the mafic minerals crystallized from a fractionated trapped liquid. The chemical, isotopic, and mineralogical data place stringent constraints on the nature of genetically related rocks and the relationship of anorthosites to other members of the ANT suite does not appear to be one ofsimple fractionation. The data presented in this paper are consistent with the Taylor-Jake?model of lunar evolution.     

Eruptive style of the young high-Mg basaltic-andesite Pelagatos scoria cone,southeast of México City

The eruption of the Pelagatos scoria cone in the Sierra Chichinautzin monogenetic field near the southern suburbs of Mexico City occurred less than 14,000 years ago. The eruption initiated at a fissure with an effusive phase that formed a 7-km-long lava flow, and continued with a phase of alternating and/or simultaneous explosive and effusive activity that built a 50-m-high scoria cone on the western end of the fissure and formed a compound lava flow-field near the vent. The eruption ended with the emplacement of a short lava flow that breached the cone and was accompanied by weak explosions at the crater. Products consist of a microlite-rich high-Mg basaltic andesite. Samples were analyzed to determine the magma’s initial properties as well as the effects of degassing-induced crystallization on eruptive style. Although distal ash fallout deposits from this eruption are not preserved, a recent quarry exposes a large section of the scoria cone. Detailed study of exposed layers allows us to elucidate the mode of cone-building activity. Petrological and textural data, combined with models calibrated by experimental work and melt-inclusion analyses of similar magmas elsewhere, indicate that the magma was initially hot (>1,200°C), gas-rich (up to 5 wt.% H₂O), crystal-poor (~10 vol.% Fo₉₀ olivine phenocrysts) and thus poorly viscous (40–80 Pa s). During the early phase, low magma ascent velocity at the fissure vent allowed low-viscosity magma to degas and crystallize during ascent, producing lava flows with elevated crystal contents at T < 1,100°C, and blocky surfaces. Later, the closure of the fissure by cooling dikes focused the magma flow at a narrow section of the fissure. This led to an increased magma ascent velocity. Rapid and shallow degassing (<3 km deep) triggered ~40 vol.% microlite crystallization. Limited times for gas-escape and higher magma viscosity (6 × 10⁵–4 × 10⁶ Pa s) drove strong explosions of highly (60–80 vol.%) and finely vesicular magma. Coarse clasts broke on landing, which implies brittle behavior due to complete solidification. This requires sufficient time to cool and in turn implies ejection heights of over 1 km, which is much higher than “normal” Strombolian activity. Hence, magma viscosity significantly impacts eruption style at monogenetic volcanoes because it affects the kinetics of shallow degassing. The long-lasting eruptions of Jorullo and Paricutin, which produced similar magmas in western México, were more explosive. This can be related to higher magma fluxes and total erupted volumes. Implications of this study are important because basaltic andesites are commonly erupted to form monogenetic scoria cones of the Trans-Mexican Volcanic Belt.     

Cation distribution in low-calcium pyroxenes: Dependence on temperature and calcium content and the thermal history of lunar and terrestrial pigeonites

Four pyroxenes with compositions En₄₈Fs₄₈Wo₄, En_47·5Fs_47·5Wo₅, En₄₅Fs₄₅Wo₁₀ and En₄₀Fs₄₀Wo₂₀, synthesized at 1200°C at atmospheric pressure, were heat-treated at 500, 600, 700, and 800°C for various lengths of time. These pyroxenes are variously ordered with respect to Fe²⁺ and Mg²⁺ without unmixing. The Fe²⁺-Mg²⁺ distribution over the two nonequivalent sites M1 and M2, determined through Mössbauer spectroscopy, is found to be a function of both temperature and concentration of Ca²⁺ at the M2 site. The preference of Fe²⁺ for the M2 site increases with decreasing temperature and increasing Ca²⁺. These data can be used to determine cation equilibration temperatures of lunar and terrestrial pigeonites. The lunar pigeonites usually indicate equilibration temperatures of 700–860°C, except the pigeonite from rock 14053, which may have been subjected to shock heating due to meteoritic impact.     

Magnetic studies on Kilauea Iki lava lake,Hawaii

Ground surveys made during August, 1961, show large vertical magnetic intensity anomalies associated with the partly lava filled crater of Kilauea Iki. A vertical magnetic variation of 11,600 gammas occurs along a north-south profile across the crater, the maximum being on the north rim of the crater and the minimum on the south edge of the encrusted lava lake below the south rim. An east-west profile shows less vertical magnetic variation, with lake-surface measurements 1500 to 2500 gammas lower than measurements on the east rim of the crater. Computed anomalies using two-dimensional potential field graticules are in good agreement with the observed anomalies and support the following conclusions: 1) Average measured values of remanent magnetization of 10^?2 cgs units and susceptibilities of 10^?3 cgs units give reasonable magnitudes to the computed anomalies. 2) The remanent magnetization is parallel to the earth’s present magnetic field. 3) The maximum vertical magnetic field value in the north-south profile is the result of reinforcement of the positive terrain effect of the north rim of the crater and the positive edge effect of the north side of the lava lake. 4) The minimum value in the same profile is the result of reinforcement of the negative terrain effect at the base of the south rim of the crater and the negative edge effect of the south side of the lava lake. 5) Variation in the east-west magnetic profile is less because the terrain and edge effects of the horizontal components of the earth’s magnetic field and remanent magnetization approach zero. Changes in vertical magnetic field values as the lake solidifies will be maximum at the north edge of the lava lake, but more consistent changes of the opposite sign will occur on the south side of the lava lake. Total change will be approximately + 2300 gammas between the August 1961 measurement at station S6 and the value at that point when the entire lava lake has cooled below 400° C. The maximum rate of change at station S6 will occur when the 500° C isotherm is 35 to 65 meters below the surface and will be about 28 gammas per meter of lowering of the 500°C surface. Because of the steep magnetic anomalies associated with the lava lake and crater rims, the permanent magnetization presently forming in the cooling lake crust will have inclinations as much as 12° less than the average 37.5° inclination in the Kilauea area.     

Tetrahedral occupancy of ferric iron in (Mg,Fe)O: Implications for point defects in the Earth's lower mantle

We investigated the concentration and site occupation of ferric iron (Fe³⁺) in (Mg,Fe)O to understand the influence of point defects on transport properties such as atomic diffusion, electrical conductivity and viscosity. We conducted Mössbauer spectroscopy of (Mg_0.8Fe_0.2)O single crystals synthesized at temperatures from 1673 to 2273 K and pressures from 5 to 15 GPa with Re–ReO₂ and Mo–MoO₂ oxygen fugacity buffers. The isomer shift of the Mössbauer spectra suggests that Fe³⁺ occupies mostly the tetrahedral site at reduced conditions and both the octahedral and tetrahedral sites at oxidized conditions. We formulate a thermodynamic model of point defect dissolution in (Mg,Fe)O which suggests that unassociated tetrahedral Fe³⁺ is more stable than unassociated octahedral Fe³⁺ at high-pressure and low oxygen fugacity due to the effect of configurational entropy. The pressure dependence of Fe³⁺ concentration indicates a change in the dominant site occupancy of Fe³⁺: (1) Fe³⁺ in the tetrahedral site, (2) Fe³⁺ in the octahedral site, and (3) defect clusters of Fe³⁺ and cation vacancy, in the order of increasing oxygen fugacity and decreasing pressure. This is in reasonable agreement with previously reported experiments on Fe³⁺ concentration, Mg–Fe interdiffusivity and electrical conductivity. We consider it plausible that (Mg,Fe)O accommodates Fe³⁺ in the tetrahedral site down to the lower mantle. Based on our results and available experimental data, we discuss the solubility competition between Fe³⁺ and protons (H⁺), and its implications for transport properties in the lower mantle.     

Iron oxidation state in lower mantle mineral assemblages: II. Inclusions in diamonds from Kankan, Guinea

Inclusions of ferropericlase and former (Mg,Fe)(Si,Al)O₃ perovskite in diamonds from Kankan, Guinea believed to originate in the lower mantle were studied using Mössbauer spectroscopy to determine Fe³⁺/ΣFe. Fe³⁺ concentration in the (Mg,Fe)(Si,Al)O₃ inclusion is consistent with empirical relations relating Fe³⁺/ΣFe to Al concentration, supporting the inference that it crystallised in the perovskite structure at lower mantle conditions. In ferropericlase there is a nearly linear variation of trivalent cation abundance with monovalent cation abundance, suggesting a substitution of the form Na_0.5M_0.5³⁺O (M=Fe³⁺, Cr³⁺, Al³⁺). Excess positive charge is likely balanced by cation vacancies, where their abundance is observed to increase with increasing iron concentration, consistent with high-pressure experiments. The abundance of cation vacancies is related to oxygen fugacity, where ferropericlase inclusions from Kankan and São Luiz (Brazil) are inferred to have formed at conditions more oxidising than Fe-(Mg,Fe)O equilibrium, but more reducing than Re-ReO₂ equilibrium. Fe²⁺/Mg partition coefficients between (Mg,Fe)(Si,Al)O₃ and ferropericlase were calculated for inclusions co-existing in the same diamond using Mössbauer data and empirical relations based on high-pressure experimental work. Most values are consistent with high-pressure experiments, suggesting that these inclusions equilibrated at lower mantle conditions. The measured ferropericlase Fe³⁺ concentrations are consistent with diamond formation in a region of redox gradients, possibly arising from the subduction of oxidised material into reduced lower mantle. Reduction of carbonate to form ferropericlase and diamond is consistent with a slight shift of Kankan δ¹³C values to isotopically heavy compositions compared to the worldwide dataset, and could supply the oxygen necessary to satisfy the high Fe³⁺ concentration in (Mg,Fe)(Si,Al)O₃ perovskite, as well as account for the high proportion of ferropericlase in the lower mantle paragenesis. The heterogeneity of lower mantle diamond sources indicates that the composition of lower mantle diamonds do not necessarily reflect those of the bulk mantle.     

Luminescence dating of scoria fall and lahar deposits from Somma–Vesuvius,Italy

Luminescence dating has been applied to scoria and lahar deposits from Somma–Vesuvius, Italy. Samples include scoria from the AD472 and 512 (or 536) eruptions and lahar deposits. In order to find a stable luminescence signal which is less affected by anomalous fading, infrared stimulated luminescence (IRSL) signals at elevated temperatures after bleaching with IR at 50 °C (termed post-IR IRSL; pIRIR) were tested at different preheat and elevated stimulation temperatures. The fading rates of both IRSL and pIRIR signals reduced dramatically with increasing preheat and pIRIR stimulation temperatures. A pIRIR signal measured at 290 °C after a preheat at 320 °C (60 s) and an IR stimulation at 50 °C (100 s) was selected to calculate the equivalent dose (D_e). The gamma spectrometry results indicate that the U-series nuclides are not in equilibrium and there is a large ²²⁶Ra excess. The dose rates and ages were calculated by assuming a ²²⁶Ra excess (over its parent ²³⁰Th) at deposition, and that this unsupported excess then decayed to the present level. The resulting luminescence ages of the two scoria samples agreed with the expected ages, and the ages of the lahar deposits indicate that they are associated with the AD1631 eruption.     

Mafic and ultramafic xenoliths in San Bartolo lava field: New insights on the ascent and storage of Stromboli magmas

Mafic and ultramafic xenoliths are well represented within a large basaltic lava field of Stromboli. These basalts, known as San Bartolo lavas, show a high-K calc-alkaline (HKCA) affinity and were erupted <5 ka BP. Xenoliths consist of olivin-gabbro, gabbronorite, anorthosite, dunite, wehrlite and clinopyroxenite. Thermobarometric estimates for the crystallization of gabbroic materials show minima equilibration pressures of 0.17–0.24 GPa, at temperatures ranging from 940 to 1,030°C. These materials interacted with hydrous ascending HKCA basaltic magmas (with temperatures of 1,050–1,100°C) at pressures of about 0.2–0.4 GPa. These pressure regimes are nearly identical to those found for the crystallization of phenocrystic phases within HKCA basaltic lavas. Gabbroic inclusions are regarded as cumulates and represent crystallized portions of earlier HKCA Strombolian basalts.Dunite and wehrlite show porphyroclastic-heterogranular textures, whereas the clinopyroxenite exhibit a mosaic-equigranular texture typical of mantle peridotites. These ultramafic materials are in equilibrium with more primitive basaltic magmas (under moderately hydrous and anhydrous conditions) at pressures of 0.8–1.2 GPa, which is below the crust-mantle transition, located at about 20 km depth under Stromboli.Major and trace element distributions indicate comagmatism between the host basaltic lava and the mafic and ultramafic inclusions. REE patterns for mafic nodules are relatively regular and overlap the field of basaltic lavas (HKCA). They show moderate to high LREE enrichments and moderate enrichments in HREE relative to chonrites. Spider diagrams also show significant similarities between the lavas and the mafic-ultramafic xenoliths as well.During their ascent, primitive Strombolian magmas may be stored in upper-mantle regions where they interact with peridotitic materials and partly differentiate (to give dunite and wehrlite) before migrating to upper crustal levels. In this region, hydrous basaltic magmas (with estimated water contents of 2–3.5 wt%) are stored in the subvolcanic environment, and are allowed to crystallize the gabbroic materials before reaching the surface under nearly anhydrous conditions.An erratum to this article can be found at     

Microlites and "nanolites" in rhyolitic glass: microstructural and chemical characterization

 We used transmission electron microscopy (TEM) and scanning electron microscopy (SEM) to study magmatic crystals in the Ben Lomond rhyolite lava dome, Taupo Volcanic Center, New Zealand. Using TEM and SEM to investigate the size distributions of these crystals, we identified three size populations: microphenocrysts (>1.2 μm wide), microlites (>0.6 μm wide), and smaller crystals (<0.6 μm wide) which we term "nanolites". The predominant mineral phases of the microlites and nanolites are augites, pigeonites, and hypersthenes. The compositions and microstructures within these pyroxenes indicate disequilibrium crystallization at approximately 850–900  °C and undercoolings as high as 300  °C from equilibrium crystallization temperatures. Complex microstructures resulting from subsolidus reactions in augite and pigeonite are consistent with moderate cooling rates within the upper obsidian layer of the Ben Lomond rhyolite dome. This study demonstrates the existence of sub-micron magmatic crystals in a rhyolite and illustrates the potential of TEM to provide unique information about the crystallization and cooling histories of glassy volcanic rocks. Received: May 8, 1995 / Accepted: November 27, 1995     

Systematic change of foraminiferal Mg/Ca ratios across a strong salinity gradient

The Mg/Ca ratio of foraminiferal calcite is an important proxy for estimating past ocean temperatures. Used in conjunction with δ¹⁸O of foraminiferal calcite it allows deconvolution of temperature and ice-volume signals to infer past ocean temperatures and salinities (assuming the relationship between seawater δ¹⁸O and salinity is known). Such work assumes that temperature is the only, or at least the dominant, environmental controller of foraminiferal Mg/Ca. The semi-enclosed Mediterranean Sea, where salinity varies from 36 to 40 psu over a seasonal temperature range of between only 5 °C to 8 °C, provides a relevant setting to test this assumption outside the laboratory. In this study, planktonic foraminifera (O. universa, G. siphonifera, G. bulloides and G. ruber (white) and (pink)) were picked from 11 box core tops spanning the Mediterranean salinity gradient and analysed for their trace-element concentrations. Mg/Ca ratios are higher, for the associated calcification temperatures, than in other regions where calibrations have been conducted and correlate poorly with calcification temperature. Mg/Ca ratios are particularly high for samples from the Eastern Mediterranean where salinity is unusually high. Correlations of Mg/Ca with the calcification salinity are statistically significant with Mg/Ca changing by 15–59% per psu, suggesting that salinity may act as a control on Mg/Ca ratios in addition to the dominant temperature control. We show that contamination by non-carbonate material and diagenetic high-Mg carbonate overgrowths cannot account for the observed trend of increasing Mg/Ca with salinity. A relationship between Mg/Ca and salinity is also suggested by re-analysis of calibrations from open-ocean settings. These new Mediterranean results are from a region with unusually high salinity but suggest that the effects of salinity on the Mg/Ca palaeothermometer should be considered even in open-ocean settings, particularly where large salinity changes occurred in the past.     

Columnar jointing in vapor-phase-altered,non-welded Cerro Galán Ignimbrite,Paycuqui, Argentina

Columnar jointing is thought to occur primarily in lavas and welded pyroclastic flow deposits. However, the non-welded Cerro Galán Ignimbrite at Paycuqui, Argentina, contains well-developed columnar joints that are instead due to high-temperature vapor-phase alteration of the deposit, where devitrification and vapor-phase crystallization have increased the density and cohesion of the upper half of the section. Thermal remanent magnetization analyses of entrained lithic clasts indicate high emplacement temperatures, above 630°C, but the lack of welding textures indicates temperatures below the glass transition temperature. In order to remain below the glass transition at 630°C, the minimum cooling rate prior to deposition was 3.0 × 10⁻³–8.5 × 10⁻²°C/min (depending on the experimental data used for comparison). Alternatively, if the deposit was emplaced above the glass transition temperature, conductive cooling alone was insufficient to prevent welding. Crack patterns (average, 4.5 sides to each polygon) and column diameters (average, 75 cm) are consistent with relatively rapid cooling, where advective heat loss due to vapor fluxing increases cooling over simple conductive heat transfer. The presence of regularly spaced, complex radiating joint patterns is consistent with fumarolic gas rise, where volatiles originated in the valley-confined drainage system below. Joint spacing is a proxy for cooling rates and is controlled by depositional thickness/valley width. We suggest that the formation of joints in high-temperature, non-welded deposits is aided by the presence of underlying external water, where vapor transfer causes crystallization in pore spaces, densifies the deposit, and helps prevent welding.     

Nature and origin of cone-forming volcanic breccias in the Te Herenga Formation, Ruapehu, New Zealand

 Volcanic breccias form large parts of composite volcanoes and are commonly viewed as containing pyroclastic fragments emplaced by pyroclastic processes or redistributed as laharic deposits. Field study of cone-forming breccias of the andesitic middle Pleistocene Te Herenga Formation on Ruapehu volcano, New Zealand, was complemented by paleomagnetic laboratory investigation permitting estimation of emplacement temperatures of constituent breccia clasts. The observations and data collected suggest that most breccias are autoclastic deposits. Five breccia types and subordinate, coherent lava-flow cores constitute nine, unconformity-bounded constructional units. Two types of breccia are gradational with lava-flow cores. Red breccias gradational with irregularly shaped lava-flow cores were emplaced at temperatures in excess of 580  °C and are interpreted as aa flow breccias. Clasts in gray breccia gradational with tabular lava-flow cores, and in some places forming down-slope-dipping avalanche bedding beneath flows, were emplaced at varying temperatures between 200 and 550  °C and are interpreted as forming part of block lava flows. Three textural types of breccia are found in less intimate association with lava-flow cores. Matrix-poor, well-sorted breccia can be traced upslope to lava-flow cores encased in autoclastic breccia. Unsorted boulder breccia comprises constructional units lacking significant exposed lava-flow cores. Clasts in both of these breccia types have paleomagnetic properties generally similar to those of the gray breccias gradational with lava-flow cores; they indicate reorientation after acquisition of some, or all, magnetization and ultimate emplacement over a range of temperatures between 100 and 550  °C. These breccias are interpreted as autoclastic breccias associated with block lava flows. Matrix-poor, well-sorted breccia formed by disintegration of lava flows on steep slopes and unsorted boulder breccia is interpreted to represent channel-floor and levee breccias for block lava flows that continued down slope. Less common, matrix-rich, stratified tuff breccias consisting of angular blocks, minor scoria, and a conspicuously well-sorted ash matrix were generally emplaced at ambient temperature, although some deposits contain clasts possibly emplaced at temperatures as high as 525  °C. These breccias are interpreted as debris-flow and sheetwash deposits with a dominant pyroclastic matrix and containing clasts likely of mixed autoclastic and pyroclastic origin. Pyroclastic deposits have limited preservation potential on the steep, proximal slopes of composite volcanoes. Likewise, these steep slopes are more likely sites of erosion and transport by channeled or unconfined runoff rather than depositional sites for reworked volcaniclastic debris. Autoclastic breccias need not be intimately associated with coherent lava flows in single outcrops, and fine matrix can be of autoclastic rather than pyroclastic origin. In these cases, and likely many other cases, the alternation of coherent lava flows and fragmental deposits defining composite volcanoes is better described as interlayered lava-flow cores and cogenetic autoclastic breccias, rather than as interlayered lava flows and pyroclastic beds. Reworked deposits are probably insignificant components of most proximal cone-forming sequences. Received: 1 October 1998 / Accepted: 28 December 1998     

Cone-building block-and-ash flows: the Senyama volcanic products of O’e Takayama volcano,SW Japan

The Senyama volcanic products of the late Pliocene to early Pleistocene O’e Takayama volcano overlie a 100-m-thick, late Pliocene coastal quartz-sandstone and are intruded by an early Pleistocene dacite dome. The Senyama volcanic products are the remains of a cone that retains a basal part 1.5 km across and 150–250 m high from the substrate. The cone comprises dacite block-and-ash flow deposits and minor base-surge deposits occur at the base. Single beds of the block-and-ash flow deposits are 1–16 m thick and dip inward 20–40° at the base of the cone and inward or outward 10–20° at the summit. Juvenile fragments in the block-and-ash flow deposits are non- to poorly vesicular and commonly have curviplanar surfaces and prismatic joints extending inward from the surfaces, which imply quenching and brittle fracturing of dacite lava. They are variably hydrothermally altered. Nevertheless, juvenile blocks appear to retain a uniform direction of the magnetization vector residual during thermal demagnetization between 280°C and 625°C. At the time of the eruption, the well-sorted sand of the substrate was at the coast and a good aquifer that facilitated explosive interaction of water and the ascending dacite lava. The mechanism of the explosion perhaps involved thermal contraction cracking of the dacite lava, water-inflow into the interior of the lava, and explosive expansion of the water. Initial phreatomagmatic explosions opened the vent. Succeeding phreatomagmatic or phreatomagmatic–vulcanian explosions produced block-and-ash flow deposits around the vent. Hydrothermal silver-ore deposits and manganese-oxide deposits occur in the Senyama volcanic products and the underlying sandstone, respectively. They could represent post-eruptive activity of the hydrothermal system developed in and around the cone.