Basaltic explosive volcanism: Constraints from deposits and models

Basaltic pyroclastic volcanism takes place over a range of scales and styles, from weak discrete Strombolian explosions (～10²–10³ kg s^?1) to Plinian eruptions of moderate intensity (10⁷–10⁸ kg s^?1). Recent well-documented historical eruptions from Etna, Kīlauea and Stromboli typify this diversity. Etna is Europe's largest and most voluminously productive volcano with an extraordinary level and diversity of Strombolian to subplinian activity since 1990. Kīlauea, the reference volcano for Hawaiian fountaining, has four recent eruptions with high fountaining (>400 m) activity in 1959, 1960, 1969 (–1974) and 1983–1986 (–2008); other summit (1971, 1974, 1982) and flank eruptions have been characterized by low fountaining activity. Stromboli is the type location for mildly explosive Strombolian eruptions, and from 1999 to 2008 these persisted at a rate of ca. 9 per hour, briefly interrupted in 2003 and 2007 by vigorous paroxysmal eruptions. Several properties of basaltic pyroclastic deposits described here, such as bed geometry, grain size, clast morphology and vesicularity, and crystal content are keys to understand the dynamics of the parent eruptions.The lack of clear correlations between eruption rate and style, as well as observed rapid fluctuations in eruptive behavior, point to the likelihood of eruption style being moderated by differences in the fluid dynamics of magma and gas ascent and the mechanism by which the erupting magma fragments. In all cases, the erupting magma consists of a mixture of melt and gaseous bubbles. The depth and rate of degassing, melt rheology, bubble rise and coalescence rates, and extent of syn-eruptive microlite growth define complex feedbacks that permit reversible shifts between fragmentation mechanisms and in eruption style and intensity. However, many basaltic explosive eruptions end after an irreversible shift to open-system outgassing and microlite crystallization in melt within the conduit.Clearer understanding of the factors promoting this diversity of basaltic pyroclastic eruptions is of fundamental importance in order to improve understanding of the range of behaviors of these volcanoes and assess hazards of future explosive events at basaltic volcanoes. The three volcanoes used for this review are the sites of large and growing volcano-tourism operations and there is a public need both for better knowledge of the volcanoes’ behavior and improved forecasting of the likely course of future eruptions.     

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.     

Fast kimberlite ascent rates estimated from hydrogen diffusion profiles in xenolithic mantle olivines from southern Africa

Fourier transform infrared spectrometry (FTIR) analyses of olivines from peridotite xenoliths found in southern African kimberlites indicate 0 to 80 ppm H₂O concentrations. OH absorbance profiles across olivine grains show homogeneous H contents from core to edge for most samples. In one sample the olivines are H-free, while another has olivines characterized by lower H contents at the grain edges compared to the cores, indicating H loss during transport of the xenolith to the surface. Flat or near-flat H profiles place severe constraints on the duration of H loss from olivine grains, with implications for kimberlite magma ascent rates. Diffusion equations were used to estimate times of H loss of about 4 h for the sample with heterogeneous olivine H contents. Resulting kimberlite ascent rates are calculated to be 5-37 m s⁻¹ minimum, although these estimates are highly dependent on volatile contents and degassing behavior of the host kimberlite magma. Xenolithic olivines from alkali basalts generally have lower H contents and more pronounced H diffusion profiles than do those from kimberlites. This difference is likely caused by higher magma temperatures and lower ascent rates of alkali basalts compared to kimberlites.     

Calcite precipitation instability under laminar, open-channel flow

We present a 2D numerical model for the growth of calcite from supersaturated aqueous solutions under laminar, open-channel flow conditions. The model couples solution chemistry, precipitation at solution/calcite interfaces, hydrodynamics, diffusion and degassing. The model output is compared with experimental results obtained using an oversaturated calcite solution produced by mixing CaCl₂ and Na₂CO₃. The precipitation rate is observed to increase when the supersaturated solution flows over an obstruction, leading to a growth instability that causes the formation of terraces. At relatively high flow rates, the most important mechanism for this behaviour seems to be hydrodynamic advection of dissolved species either towards or away from the calcite surface, depending on location relative to the obstruction, which deforms the concentration gradients. At lower flow rates, steepening of diffusion gradients around protrusions becomes important. Enhanced degassing over the obstruction due to shallowing and pressure drop is not important on small scales. Diffusion controlled transport close to the calcite surface can lead to a fingering-type growth instability, which generates porous textures. Our results are consistent with existing diffusive boundary layer theory, but for flow over non-smooth surfaces, simple calcite precipitation models that include empirical correlations between fluid flow rate and calcite precipitation rate are inaccurate.     

Magma Ascent and Crustal Accretion at Ultraslow-Spreading Ridges: Constraints from Plagioclase Ultraphyric Basalts from the Arctic Mid-Ocean Ridge

Plagioclase ultraphyric basalts (PUBs) with up to 54% plagioclasephenocrysts were dredged in the rift valley and adjacent flanksof the ultraslow-spreading Mohns and Knipovich ridges. The PUBsshow large variations in crystal morphologies and zoning. Thelarge variations suggest that single basalt samples containa mixture of plagioclase crystals that aggregated at differentlevels in the magma conduits. Resorbed crystals and repeatedreverse zones suggest that the magma reservoirs were replenishedand heated several times. Thin concentric zones with melt inclusions,and sharp reductions in the anorthite content of 3–7%,are common between the reverse zones. These zones, and skeletalcrystals with distinctly lower anorthite contents than massivecrystals, are interpreted to be the result of rapid crystalliztionduring strong undercooling. The changes between short periodsof cooling and longer periods with reheating are explained bymultiple advances of crystal-rich magma into cool regions followedby longer periods of gradual magma inflow and temperature increase.The porphyritic basalts are characterizd by more depleted andmore fractionated compositions than the aphyric basalts, withlower (La/Sm)_N, K₂O and Mg-numbers. This relationship, and theobservation that PUBs are sampled only close to segment centresalong these ridges, suggests that the PUBs formed by higherdegrees of melting and evolved in more long-lived magma reservoirs.We propose that the zoning patterns of plagioclase crystalsand crystal morphologies of these PUBs reflect the developmentand flow of magma through a stacked sill complex-like conduitsystem, whereas the aphyric equivalents represent later flowof magma through the conduit. The formation of voluminous higher-degreemelts may trigger the development of the magma conduits andexplain the generally depleted compositions of PUB magmas. KEY WORDS: basalt; mineral chemistry; MORB; magma mixing; magma chamber; major element     

National-scale estimation of potentially harmful element ambient background concentrations in topsoil using parent material classified soil:stream–sediment relationships

Regulatory authorities require estimates of ambient background concentrations (ABCs) of potentially harmful elements (PHEs) in topsoil; such data are currently not available in many countries. High resolution soil geochemical data exist for only part of England and Wales, whilst stream sediment data cover the entire landscape. A novel methodology is presented for estimating soil equivalent ABCs for PHEs from high-resolution (HR) stream sediment geochemical data grouped by common parent materials (PM), using arsenic (As) as an example. Geometric mean (GM) values for local PM groups are used to investigate different approaches for transforming sediment to soil equivalent concentrations. Holdout validation is used to assess: (i) the optimum number of samples for calculating local GM values, and (ii) the optimum scale at which to group data when using linear regression analysis to estimate GM soil ABCs from local sediment geochemical values. Holdout validation showed that the smallest differences were generally observed when five observations were used to calculate the GM and that these should be grouped over the smallest possible area in order to encompass soils over PMs with elevated GM As concentrations. Geometric mean ABCs are estimated and mapped for As in mineral soil across all of England and Wales within delineations of PM polygons. Errors for the estimation of soil equivalent GM As ABCs based on sediment data for an independent validation set were of a similar magnitude to those from holdout validation applied to the original data suggesting the approach is robust. The estimates of soil equivalent ABCs suggest that As exceeds the regulatory threshold used in risk assessments for residential land use (20 mg kg⁻¹) across 16% of the landscape of England and Wales. The applicability of the method for cognate landscapes, and potential refinements is discussed.     

Methodology of hydrogeological characterization of deep carbonate aquifers as potential reservoirs of groundwater. Case of study: the Jurassic aquifer of <Emphasis Type="Italic">El Maestrazgo</Emphasis> (Castellón,Spain)

A methodology for the characterization of deep carbonate aquifers has been developed and applied to El Maestrazgo Jurassic aquifer in Castellón, Spain. Characterization of these aquifer formations, located at more than 300 m deep, consisted of a previous phase of compilation, analysis and synthesis of the existing information about the area, followed by a coordinated combination of different speciality studies: geology, stratigraphy, structural analysis, hydrogeology, hydrochemistry, geophysics and remote sensing. Geological studies included geological mapping, definition of stratigraphical units and facies and structural analysis. The aim of the hydrogeology study was to define aquifer formations, recharge area, aquifer points inventory and groundwater flow directions for the establishment of piezometric and water quality observation nets. Special techniques were applied, like thermal infrared aerial images and the evaluation of submarine groundwater discharge by means of natural radium isotopes. Hydrochemical techniques, including majority elements characterization and stable isotopes (¹⁸O, ²H and ³H) determination, allowed classifying hydrochemical facies and establishing a renewal pattern for water within the system. Geophysics was useful in determining the aquifer geometry, the features of the basement and the petrophysical characteristics of the geological formations. Preliminary results show an important tectonic complexity and the possibilities for groundwater uses in the area of study.     

The role of the volatile phase for REE and Y fractionation in low-silica carbonate magmas: implications from natural carbonatites, Namibia

Summary Two bimodal carbonatite complexes in Namibia of Cretaceous age are explored as to the presence and composition of a coexisting carbonatitic fluid. The Kalkfeld and Ondurakorume complexes contain both Ca- and Mg/Fe-carbonatites, composed of calcite alone or calcite with ferroan dolomite, fluorapatite and strontianite. The major element evolution in the bulk rocks from s?vites to beforsites is due to crystallization of calcite and fluorapatite. All carbonatites show a negative Y anomaly in normalised REE plots. Fractionation is accompanied by successively lower HREE contents between Tb and Yb, expressed by the ratios Nd/Ho and Ho/Lu. The evolution of this downward-facing hump goes along with decreasing Y contents in bulk rocks and minerals. All this requires an additional phase coexisting with the carbonate liquid during fractionation. Comparison between the bulk rocks and the expelled fluid shows that the latter had preferentially accumulated the HREE and Y. Further evidence for this process are hydrothermal, HREE, Y-rich fluorites in other carbonatite complexes which reflect the composition of the expelled fluid. The high strength of fluoride complexes suggests that fluoride complexing in the carbonatitic fluid is the process responsible for extracting HREE and Y from the carbonatite magma, leaving fractionated carbonatite rocks depleted in these elements. The geochemical evolution of carbonatite magmas along fractionation has therefore to be considered in a melt-mineral-fluid system. Correspondence: B. Bühn, Instituto de Geociências, Universidade de Brasilia, Campus Universitário Darcy Ribeiro, Brasilia 70910-900, Brazil