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.     

Reaction mechanism of uranyl in the presence of zero-valent iron nanoparticles

The current study provides an investigation of abiotic reduction of an oversaturated uranyl solution driven by iron nanoparticle oxidation. The reactivity of nano-scale zero-valent iron (ZVI) under mildly oxic conditions (1.2% O₂ and 0.0017% CO₂) was studied in 1000 ppm uranyl solution in the pH range 3-7, at reaction times from 10 min to 4 h. Reductive precipitation of UO₂ was observed as the main process responsible for the removal of uranium from solution with the kinetics of reaction becoming increasingly favourable at higher pH. Despite working with an oversaturated uranium solution, the precipitation of UO₂ occurred in preference to precipitation of UO₃·2H₂O (metaschoepite) at reaction times between 1 and 4 h and for uranyl solutions initially set up at pH ?5. Characterisation of both solid and solution phases was performed using X-ray photoelectron spectroscopy (XPS), focused ion beam (FIB) imaging, X-ray diffraction (XRD) and inductively coupled plasma atomic emission spectroscopy (ICP-AES).     

Anaerobic oxidation of methane (AOM) in marine sediments from the Skagerrak (Denmark): II. Reaction-transport modeling

A steady-state reaction-transport model is applied to sediments retrieved by gravity core from two stations (S10 and S13) in the Skagerrak to determine the main kinetic and thermodynamic controls on anaerobic oxidation of methane (AOM). The model considers an extended biomass-implicit reaction network for organic carbon degradation, which includes extracellular hydrolysis of macromolecular organic matter, fermentation, sulfate reduction, methanogenesis, AOM, acetogenesis and acetotrophy. Catabolic reaction rates are determined using a modified Monod rate expression that explicitly accounts for limitation by the in situ catabolic energy yields. The fraction of total sulfate reduction due to AOM in the sulfate-methane transition zone (SMTZ) at each site is calculated. The model provides an explanation for the methane tailing phenomenon which is observed here and in other marine sediments, whereby methane diffuses up from the SMTZ to the top of the core without being consumed. The tailing is due to bioenergetic limitation of AOM in the sulfate reduction zone, because the methane concentration is too low to engender favorable thermodynamic drive. AOM is also bioenergetically inhibited below the SMTZ at both sites because of high hydrogen concentrations (∼3-6 nM). The model results imply there is no straightforward relationship between pore water concentrations and the minimum catabolic energy needed to support life because of the highly coupled nature of the reaction network. Best model fits are obtained with a minimum energy for AOM of ∼11 kJ mol⁻¹, which is within the range reported in the literature for anaerobic processes.     

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     

Comparison of arsenic concentrations in simultaneously-collected groundwater and aquifer particles from Bangladesh,India, Vietnam,and Nepal

One of the reasons the processes resulting in As release to groundwater in southern Asia remain poorly understood is the high degree of spatial variability of physical and chemical properties in shallow aquifers. In an attempt to overcome this difficulty, a simple device that collects groundwater and sediment as a slurry from precisely the same interval was developed in Bangladesh. Recently published results from Bangladesh and India relying on the needle-sampler are augmented here with new data from 37 intervals of grey aquifer material of likely Holocene age in Vietnam and Nepal. A total of 145 samples of filtered groundwater ranging in depth from 3 to 36 m that were analyzed for As (1–1000 μg/L), Fe (0.01–40 mg/L), Mn (0.2–4 mg/L) and S (0.04–14 mg/L) are compared. The P-extractable (0.01–36 mg/kg) and HCl-extractable As (0.04–36 mg/kg) content of the particulate phase was determined in the same suite of samples, in addition to Fe(II)/Fe ratios (0.2–1.0) in the acid-leachable fraction of the particulate phase. Needle-sampler data from Bangladesh indicated a relationship between dissolved As in groundwater and P-extractable As in the particulate phase that was interpreted as an indication of adsorptive equilibrium, under sufficiently reducing conditions, across 3 orders of magnitude in concentrations according to a distribution coefficient of 4 mL/g. The more recent observations from India, Vietnam and Nepal show groundwater As concentrations that are often an order of magnitude lower at a given level of P-extractable As compared to Bangladesh, even if only the subset of particularly reducing intervals characterized by leachable Fe(II)/Fe >0.5 and dissolved Fe >0.2 mg/L are considered. Without attempting to explain why As appears to be particularly mobile in reducing aquifers of Bangladesh compared to the other regions, the consequences of increasing the distribution coefficient for As between the particulate and dissolved phase to 40 mL/g for the flushing of shallow aquifers of their initial As content are explored.     

Sedimentation settings of the Vendian Vanavara formation,the Siberian platform

A model of sedimentation settings is elaborated for siliciclastic deposits of the Vendian Vanavara Formation, the Katanga saddle, inner areas of the Siberian platform. Four lithologic complexes are distinguished in the formation. The lower complex is composed of proluvial continental deposits exemplifying a dejection cone of ephemeral streams. Its eroded surface is overlain by second complex largely represented by sandstones of coastal zone, which grade upward into siltstones and shales of deeper sedimentation settings (third complex). Sea transgression advanced in northeastern direction. The fourth complex resting with scouring on the third one was deposited in settings of a spacious shallow-water sea zone: in a tidal flat, sand shoals and islands. Sedimentological data are used to correlate more precisely the Vendian siliciclastic deposits of the Katanga saddle and northeastern Nepa-Botuoba anteclise, and to verify subdivision of the Vanavara Formation into subformations and character of its boundary with the overlying Oskoba Formation.