Synergy of multiple geophysical approaches to unravel explosive eruption conduit and source dynamics - A case study from Stromboli

An effective approach to understanding the dynamics of explosive volcanic eruptions and the conduit systems that drive them is through synergy of multiple data sets. Three data sets that lend themselves to ease of integration are seismic, infrasonic and thermal. Although approaches involving these data have been used to record volcanological phenomena since 1862, 1955 and 1965, respectively, their integrated use has only developed since 1999. When combined, these three data sets allow constraint of shallow system geometry and the dynamics of the explosive events that occur within that system. Using Stromboli volcano (Italy) as a case study, we review the complete range of geochemical and geophysical studies that can be applied. In doing so, we aim to show how integration of these diverse studies allows insights into a plumbing system and the dynamics of the eruptive activity that the system feeds. When combined at Stromboli, these data provide constraint of multiple system parameters including chamber depths, gas and magma fluxes, shallow system magma residence times, explosion source depths, and the rise/ejection velocities of ascending gas slugs and ejecta. In turn, these results allow various conduit and eruption dynamic models to be applied and tested.The persistent and repeated mildly explosive events that characterize Stromboli have been modeled in terms of the coalescence of gas within the magma to form large gas slugs that ascend the remaining portion of the conduit to burst at the free surface. Our integrated seismic, infrasonic and thermal data sets indicate that gas coalescence occurs at a depth of ∼260 m, with a typical event frequency of ∼9/h. Infrasonic and thermal data show the explosion source to be located 20-220 m below the vent. Thermal data give emission velocities for the ejected fragments of 8-20 m/s, which converts to gas jet velocities of 23-39 m/s. Tracking these parameters in space and time shows that, although eruptions at Stromboli can be grouped into two characteristic types (simple and complex-each of which characterizes a particular crater, NE and SW, respectively), events within each type show significant short-term variability. The system does, however, appear robust, maintaining its characteristic strombolian eruption style after significant effusive phases and more energetic explosive events.     

Phosphate defects and apatite inclusions in coral skeletal aragonite revealed by solid-state NMR spectroscopy

Recent development of paleo-nutrient proxies based on the phosphorus/calcium (P/Ca) ratio in tropical- and deep-water corals (also known as cold-water corals) require an understanding of the processes by which P is incorporated into the coral skeletal aragonite. Here, we apply single- and double-resonance solid-state nuclear magnetic resonance (NMR) spectroscopy to determine the speciation of P in coral aragonite. The results show that the majority of P occurs as phosphate defects in the aragonite structure, but in many samples a significant fraction of the P occurs also in crystalline hydroxylapatite inclusions. Quantification of the amount of hydroxylapatite indicates that its presence is not related simply to external environmental factors and that it can occur at varying abundances in different parts of the same corallite. Since there is currently no model available to describe the relationship between dissolved inorganic phosphate and its incorporation as apatite inclusions into carbonates, careful screening of samples which contain only phosphate in the aragonite structure or selective microsampling could improve proxy development.     

Oxygen isotope composition of phenocrysts from Tristan da Cunha and Gough Island lavas: variation with fractional crystallization and evidence for assimilation

We have measured the δ¹⁸O values of the major phenocrysts (olivine, clinopyroxene and plagioclase) present in lavas from Tristan da Cunha and Gough Island. These islands, which result from the same mantle plume, have enriched radiogenic isotope ratios and are, therefore, prime candidates for an oxygen isotope signature that is distinct from that of MORB. Consistent differences between the δ¹⁸O values of olivine, pyroxene and feldspar in the Gough lavas show that the phenocrysts in the mafic Gough Island lavas are in oxygen isotope equilibrium. The olivines in lavas with SiO₂ <50 wt% have a mean δ¹⁸O value of 5.19‰, consistent with crystallization from a magma having the same oxygen isotope composition as MORB. Phenocrysts in all the Gough lavas show a systematic increase in δ¹⁸O value as silica content increases, which is consistent with closed-system fractional crystallization. The lack of enrichment in δ¹⁸O of the Gough magmas suggests that the mantle source contained <2% recycled sediment. In contrast, the Tristan lavas with SiO₂ >48 wt% contain phenocrysts which have δ¹⁸O values that are systematically ∼0.3‰ lower than their counterparts from Gough. We suggest that the parental mafic Tristan magmas were contaminated by material from the volcanic edifice that acquired low δ¹⁸O values by interaction with water at high temperatures. The highly porphyritic SiO₂-poor lavas show a negative correlation between olivine δ¹⁸O value and whole-rock silica content rather than the expected positive correlation. The minimum δ¹⁸O value occurs at an SiO₂ content of about 45 wt%. Below 45 wt% SiO₂, magmas evolved via a combination of assimilation, fractionational crystallization and crystal accumulation; above 45 wt% SiO₂, magmas appeared to have evolved via closed-system fractional crystallization. Received: 23 November 1998 / Accepted: 27 September 1999     

Estimation of uranium endowment by subjective geological analysis—a comparison of methods and estimates for the San Juan Basin,New Mexico

This study investigates how estimates of uranium endowment made by a geologist using an appraisal system that is based upon a formalization of geoscience and decision rules compare with estimates made by informal and unconstrained intuitive processes. The motivation for this study derives from the premise that formalization of decisions would mitigate the heuristic biases and hedging that may result from the use of unconstrained intuitive processes. Estimates of the uranium endowment of the San Juan Basin of New Mexico by four methodologies are compared in this study. These methods, ranked from top to bottom by degree of decomposition (mitigating of heuristic bias)and control on hedging, are as follows Implicit 2 1.5 × 10⁶ s.t. of U₃O₈ Implicit 1 1.6 × 10⁶ s.t. of U₃O₈ NURE (1980) 2.4 × 10⁶ s.t. of U₃O₈ Appraisal system 3.9 × 10⁶ s.t. of U₃O₈ The magnitude of expected uranium endowment estimated by these methods, ranked from smallest to largest, is in this same order. With the exception of the NURE estimates, the magnitude of the variance (uncertainty)of uranium endowment, ranked from smallest to largest, also is in this same order. These results prompt the suggestion that the more decomposed and formalized the estimation procedure, the greater the expected value and the variance of uranium endowment. Equivalently, predicating U ₃ O ₈ endowment estimation strictly upon that part of the geologist's geoscience that is useful in making U ₃ O ₈ endowment estimates and upon his understanding of the region's history produced larger estimates than have previously been reported. However, this method of estimation also shows that uncertainty about the actual state of U ₃ O ₈ endowment is much greater than previously described.     

Application of precise 142Nd/144Nd analysis of small samples to inclusions in diamonds (Finsch,South Africa) and Hadean Zircons (Jack Hills,Western Australia)

¹⁴⁶Sm–¹⁴²Nd and ¹⁴⁷Sm–¹⁴³Nd systematics were investigated in garnet inclusions in diamonds from Finsch (S. Africa) and Hadean zircons from Jack Hills (W. Australia) to assess the potential of these systems as recorders of early Earth evolution. The study of Finsch inclusions was conducted on a composite sample of 50 peridotitic pyropes with a Nd model age of 3.3 Ga. Analysis of the Jack Hills zircons was performed on 790 grains with ion microprobe ²⁰⁷Pb/²⁰⁶Pb spot ages from 3.95 to 4.19 Ga. Finsch pyropes yield 100 × ?¹⁴²Nd = ? 6 ± 12 ppm, ?¹⁴³Nd = ? 32.5, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.1150. These results do not confirm previous claims for a 30 ppm ¹⁴²Nd excess in South African cratonic mantle. The lack of a ¹⁴²Nd anomaly in these inclusions suggests that isotopic heterogeneities created by early mantle differentiation were remixed at a very fine scale prior to isolation of the South African lithosphere. Alternatively, this result may indicate that only a fraction of the mantle experienced depletion during the first 400 Myr of its history. Analysis of the Jack Hills zircon composite yielded 100 × ?¹⁴²Nd = 8 ± 10 ppm, ?¹⁴³Nd = 45 ± 1, and ¹⁴⁷Sm/¹⁴⁴Nd = 0.5891. Back-calculation of this present-day ?¹⁴³Nd yields an unrealistic estimate for the initial ?¹⁴³Nd of ? 160 ?-units, clearly indicating post-crystallization disturbance of the ¹⁴⁷Sm–¹⁴³Nd system. Examination of ^146,147Sm–^142,143Nd data reveals that the Nd budget of the Jack Hills sample is dominated by non-radiogenic Nd, possibly contained in recrystallized zircon rims or secondary subsurface minerals. This secondary material is characterized by highly discordant U–Pb ages. Although the mass fraction of altered zircon is unlikely to exceed 5–10% of total sample, its high LREE content precludes a reliable evaluation of ¹⁴⁶Sm–¹⁴²Nd systematics in Jack Hills zircons.